Amplification of LAPTM4B and YWHAZ contributes to chemotherapy resistance and recurrence of breast cancer
Abstract
Post-surgery adjuvant chemotherapy for breast cancer has effectively reduced metastatic recurrence rates1. However, a significant proportion of women suffer recurrent cancer at distant metastatic sites despite adjuvant treatment. Identification of the genes critical for tumor response to specific chemotherapy drugs is a challenge, but necessary to improve outcomes2. Using integrated genomics, we identified a small number of over-expressed and amplified genes from chromosome 8q22 significantly associated with early disease recurrence despite anthracycline-based adjuvant chemotherapy. The association was confirmed in an analysis of multiple independent cohorts. Two of these genes, the anti-apoptotic gene YWHAZ, and LAPTM4B, a novel lysosomal gene, sensitized tumor cells to anthracyclines when either was depleted by siRNA knockdown and induced drug resistance when either was over-expressed. Over-expression of LAPTM4B resulted in sequestration of drug, delaying its appearance in the nucleus. Over-expression of these two genes was associated with poor tumor response to anthracycline treatment in a neo-adjuvant chemotherapy trial in women with primary breast cancer. Our results suggest that 8q22 amplification and over-expression of LAPTM4B and YWHAZ contribute to de novo chemoresistance to anthracyclines, and are permissive for metastatic recurrence. These two genes may predict anthracycline resistance and influence selection of chemotherapy.
Breast cancer recurs at distant sites in a significant number of women who receive adjuvant chemotherapy after surgical removal of the primary breast tumor1. De novo resistance mechanisms present within tumor cells prior to treatment are key factors leading to failure of chemotherapeutic drugs to prevent metastatic recurrence. Therefore, discovery of the genomic alterations and genes contributing to de novo chemo-resistance to specific drugs is an important goal2. Although a number of multidrug resistance genes have been discovered, their over-expression is often induced during drug treatment3,4 and not useful for initial guidance of drug selection. Gene signatures generated from drug responses of tumor cell lines are reported to predict drug response in patients5–7; however, others found cell-line derived signatures are not predictive of response in clinical cases8. Repeatedly observed genomic gains or losses have identified genomic regions that may harbor genes contributing to malignant behavior and poor outcome9–12. Which genomic region(s) harbors genes that may contribute to de novo resistance to therapy is currently unknown.
We analyzed gene expression profiles of 115 breast carcinomas from women diagnosed between 2000 and 2003 and treated according to current guidelines including adjuvant chemotherapy if indicated. We performed predictive analysis of microarrays (PAM; 13) and identified 114 probes, encoding 75 known genes, differentially expressed between cases with early distant metastatic recurrence and cases without distant recurrence (Supplementary Table 1). Fifteen percent of these probes, corresponding to 12 different genes, mapped to chromosome 8q22, the only chromosomal region with statistically significant enrichment (P < 2.1e-09) of probes associated with metastatic recurrence (Fig. 1a). We applied Cox proportional hazard regression 13,14 which also demonstrated differential over-expression of these 8q22 genes in tumors with distant recurrence. These genes included CCNE2 and MTDH which are reported associated with metastatic recurrence and poor prognosis of breast cancer 15,16. The coordinate over-expression of neighboring genes often reflects chromosomal amplification. Indeed, 8q22 amplification was observed by SNP array analysis in 50 breast cancers (Supplementary Fig. 1) and expression of the 8q22 genes correlated with DNA copy number (Supplementary Table 1).
We confirmed 8q22 amplification by DNA interphase fluorescence in situ hybridization (FISH) (Fig. 1b,c) and found it in 21% of 85 breast cancers. Degree of copy gain was correlated with average expression of the 12 recurrence-associated 8q22 genes (8q gene expression index, 8qEI) (Fig. 1d and Supplementary Table 1). Kaplan-Meier analysis showed 8q22 amplification was associated with reduced metastasis-free survival in the entire cohort evaluated by FISH (Fig. 1e), in the ER− cases (Supplementary Fig. 2a), and in the women who had received anthracycline-based adjuvant chemotherapy (Fig. 1f). In multivariate analysis, amplification of 8q22 was a strong independent prognostic factor for breast cancer recurrence (Supplementary Table 2).
We sought validation in a meta analysis of six independent cohorts annotated with treatment and outcome 12,14,17–20. Kaplan-Meier analysis demonstrated a significant difference in disease-free survival between 8qEI low-expression and high-expression groups in either chemo-treated (Fig. 1g) or untreated cases (Supplemental Fig.2b). These results indicate that 8q22 amplification promotes over-expression of 8q22 genes in tumor tissue, which are associated with poor prognosis in untreated cases and inferior disease-free survival despite adjuvant chemotherapy.
To determine if 8q22 genes influence sensitivity to chemotherapy, we treated the breast cancer cell line BT549 harboring 8q gain with siRNA against the 12 candidate genes (Supplementary Fig. 3) and screened for alteration of sensitivity to chemotherapeutic drugs (Fig. 2a). Depletion of two genes significantly increased the sensitivity to anthracyclines (Fig. 2a). One of these, YWHAZ, codes for a known anti-apoptotic protein 14-3-3ζ 21,22. The second gene codes a novel lysosomal protein LAPTM4B (Lysosomal Associated Protein Transmembrane 4B), about which little is known in breast cancer. We examined 16 breast cancer cell lines and found a strong positive correlation between higher endogenous LAPTM4B mRNA level and higher IC50 (relative resistance) to anthracyclines (P < 0.00034, Supplementary Fig. 4a), a weaker or no correlation with IC50 to cisplatin and paclitaxel (P = 0.008 and 0.4, respectively; data not shown). The expression of YWHAZ in cell lines also correlated with the IC50 to doxorubicin (Supplementary Fig 4). Specific knockdown of LAPTM4B and YWHAZ in three cell lines (Fig. 2b) increased sensitivity to the anthracyclines doxorubicin (Fig. 2c upper panels) and daunorubicin (data not shown), but had weaker or no effect on sensitivity to cisplatin and paclitaxel (Fig. 2c lower panels). Knockdown of either LAPTM4B or YWHAZ significantly increased doxorubicin-induced apoptosis (Fig. 2d). Induction of apoptosis was less apparent in response to cisplatin treatment (Fig. 2d).
To investigate mechanism, intracellular localization of anthracyclines was tracked by following the autofluorescence of doxorubicin. LAPTM4B expression in cell lines correlated with both IC50 of doxorubicin (Fig. 3a) and with appearance of doxorubicin in the nucleus within 24 hours (Fig. 3b). Knockdown of LAPTM4B by siRNA in MDA-MB-231 cells resulted in a significant increase in nuclear localization of doxorubicin, detectable within 12 h to 24 h of treatment, maximal at 24 h to 36 h (12 h after withdrawal of drug), and sustained at 48 h (Fig. 3c). Knockdown of LAPTM4B in BT549 cells resulted in a similar increase in nuclear localization at 24 h (data not shown). Decreased distribution of doxorubicin into the nucleus is associated with lower phospho-H2AX (Fig. 3c,d), a marker of DNA damage response, and consistent with reduced doxorubicin-induced apoptosis and increased IC50 (Fig. 2). YWHAZ levels correlated with IC50 of doxorubicin (Fig. 3a) and knockdown of YWHAZ increased slightly the phospho-H2AX levels in drug treated cells (Fig. 3d). These results suggest one mechanism of poorer outcome for women harboring breast cancers with 8q22 amplification is increased expression of LAPTM4B and interference with nuclear accumulation of anthracyclines. YWHAZ may effect drug sensitivity through inhibition of apoptosis, consistent with reports of others21,22.
We introduced HA-tagged full length LAPTM4B and YWHAZ into partially transformed human mammary epithelial cells (HMECs) 23. Expression of either exogenous LAPTM4B or YWHAZ increased the IC50 of doxorubicin (270% increase, P < 0.002 and 394% increase, P = 0.0001, respectively) but had no significant effect on sensitivity to paclitaxel or cisplatin (Fig. 3e, f). The LAPTM4B-induced decrease in drug sensitivity parallels and is consistent with delayed appearance of anthracycline in the nucleus of LAPTM4B-overexpressing HMEC (Fig. 3g).
Kaplan Meier analysis of women treated with adjuvant chemotherapy showed that the expression of YWHAZ and LAPTM4B above median level was associated with shorter disease-free survival (Supplementary Fig. 2c). The association with poor outcome after adjuvant chemotherapy is consistent with either a prognostic effect or a role of these two genes in chemotherapy resistance.
Finally, LAPTM4B and YWHAZ were tested for their association with response to anthracyclines in a neoadjuvant (pre-operative) treatment trial of epirubicin monotherapy. The average of LAPTM4B and YWHAZ expression levels from expression array data of pretreatment tumor biopsies was evaluated for association with pathologic complete response (pCR) to epirubicin. The two-gene expression levels displayed a coherent pattern with higher levels of expression associated with absence of pCR (presence of residual disease) after epirubicin treatment (Fig. 4a, b). We evaluated the capability of the two genes to predict pCR by measuring the area under receiver operator characteristic (ROC) curves24,25(AUC); which demonstrated higher expression of the two genes is associated with absence of pCR after anthracycline chemotherapy in the cohort of 118 breast tumors (AUC 0.315, p < 0.00058, Fig. 4d). The association is more significant in 87 ER− HER2− tumors (AUC 0.241, p < 0.000062, Fig. 4e). When the two genes were analyzed separately, both genes were significantly predictive of poor response in the whole cohort and in the ER− HER2− subset, but only LAPTM4B level was predictive in the ER− HER2+ subset (Supplementary Table 3). In contrast, the expression of these two genes was not associated with treatment response to cisplatin monotherapy in a separate neo-adjuvant clinical trial in triple negative (ER− PR− HER2−) cancer cases26 (AUC 0.675, P > 0.3, Fig. 4c,f; Supplementary Table 3). Although MTDH, one of the twelve 8q22 genes, was reported to induce chemo-resistance to a broad spectrum of drugs in experimental models16, its expression was not predictive for pCR in either the epirubicin or cisplatin trials (data not shown). We analyzed a third single agent neoadjuvant therapy trial of predominantly ER+ tumors treated with docetaxel27 and found that higher levels of expression of LAPTM4B and YWHAZ were not associated with an inferior clinical response to therapy (data not shown). The results support the notion that LAPTM4B and YWHAZ over-expression is associated preferentially with poor response to anthracyclines.
LAPTM4B is similar to its family member LAPTM4A that promotes selective resistance to anthracyclines and not cisplatin in Saccharomyces cerevisiae 28. Our results show that LAPTM4B acts on anthracycline trafficking by reducing drug entry into the nucleus and decreasing drug-induced DNA damage. Higher YWHAZ expression protects cells from drug-induced apoptosis. Because they reside in proximity, amplification produces coordinated up-regulation of their various functions, together resulting in preferential resistance to anthracyclines. The results from three clinical trials support this contention, and suggest clinical options for the treatment of primary breast cancers might depend upon the status of 8q22 amplification and over-expression of these two genes in tumors. Anthracyclines appear to be reasonable treatment in tumors without 8q22 amplification, and alternatives might be selected for those whose cancers harbor amplification.
Microarray data
Microarray data sets are deposited in the NCBI GEO database under the following accession numbers: Gene expression data from the DF/HCC cases, GSE19615; SNP array data from 50 of the DF/HCC cases, GSE19594; gene expression data from neo-adjuvant epirubicin “Trial of Principle”, GSE16446; gene expression data from neo-adjuvant cisplatin trial, GSE18864.
Methods
Cohort
Primary breast tumors of 115 subjects were obtained from the NCI-Harvard Breast SPORE blood and tissue repository under Dana-Farber/Harvard Cancer Center Institutional Review Board approved protocols, with informed consent from subjects. Affymetrix U133 plus 2 gene expression array analysis was performed as described 29,30 A subset of 85 tumors were represented in tissue microarrays and used for FISH analysis (see Supplementary Methods). Fifty of the cases were analyzed by Affymetrix 10K SNP array as described31,32. A portion of the SNP and gene expression data were reported previously 29,30,31,32. Clinical and pathologic characteristics for each sample in the cohort are provided in Supplementary Table 4.
Neoadjuvant Clinical Trials
The neoadjuvant “Trial of Principle” for breast cancer is conducted in European hospitals and coordinated at the Institut Jules Bordet. This trial is registered on the clinical trials site of the US National Cancer Institute website http://clinicaltrials.gov/ct2/show/NCT00162812?term=NCT00162812&rank=1. Single agent epirubicin was given as neo-adjuvant (pre-operative) chemotherapy to 118 reportedly ER− cases. After central review, 4 of the cases were found to be ER+. Of the remaining 114 ER−, 87 cases were classified as HER2− based on low ERBB2 module score33. Pretreatment core biopsy of the primary breast tumor was performed for diagnosis and RNA isolation. At completion of chemotherapy, pathologic response was determined by microscopic examination of the excised tumor and nodes. Pathological complete response (pCR) was defined by the absence of residual invasive breast carcinoma in the breast and axillary nodes. This study has been approved by the medical ethics committee of Institute Jules Bordet and all women given written informed consent prior to study entry. Gene expression data of U133plus 2 were generated from RNA of pretreatment core biopsies.
The trial of single agent cisplatin given as neo-adjuvant chemotherapy to women with triple negative breast cancer26. Gene expression array data from pre-treatment biopsies was available for 24 cases26. Pathological response was determined by microscopy examination after chemotherapy as described above.
ROC (Receiver Operating Characteristic) curve analysis was performed to evaluate the mean level of combined LAPTM4B and YWHAZ expression, or the levels of each individual gene, for their capacity to predict pathological complete response (pCR). The association of ranked gene expression levels with pCR was evaluated by determining the area under the curve (AUC) estimated through the concordance index34; the corresponding p-value is from one-sided Wilcoxon's rank test35. The ROC curves are plotted for prediction of pCR, so a curve below the midline in low-right area of the graph indicates the ranked gene expression is associated with absence of pCR.
Statistical analysis of gene expression arrays
Affymetrix U133plus2.0 array data from tumor samples of 115 cases were classified into those with distant recurrence within 36 months of diagnosis or those without distant recurrence and at least 36 months of follow up. The 115 arrays were log-transformed, normalized and invariable genes were removed by filtering using dChip software (p-value < 0.05). Differentially expressed genes were determined by Prediction Analysis of Microarrays (PAM) 13 using the pamr package implemented in R language (http://cran.r-project.org/web/packages/pamr/index.html). PAMR implements the shrunken nearest centroid method 13. Genes were selected at a false discovery threshold that minimized a 10-fold cross-validation and test errors near the shrinkage parameter Δ = 2. The PAM score indicates the degree of statistical association for expression of each gene and metastatic recurrence. Cox proportional hazard regression analysis 13,14 was used to discover differentially expressed genes associated with time-to-recurrence.
Analysis in independent cohorts
Pooling six independent gene expression array datasets 12,14,17–20 for analysis was performed as described 36. 8q22 gene expression-based Index (8qEI) is calculated as the mean expression value of the twelve 8q22 genes identified in PAM analysis, and was scored in each of the six independent cohorts. We defined a simple median-based classifier: 8qEI level higher than the median as 8qEI high and lower than the median as 8qEI low. The classification was performed in 6 cohorts separately, and then combined into one data sheet with 1348 samples. Of these, 1130 had annotation for treatment and followup, 361 cases received adjuvant chemotherapy with or without hormonal therapy and 725 cases received no adjuvant hormonal or chemotherapy. Follow-up time was constrained to a maximum of 10 years with annotation for disease free survival. Kaplan-Meier analyses were carried out using the survival package within the R statistical package. P-values are derived using Mantel-Cox logrank test.
Additional methods are provided in supplementary material online.
Supplementary Material
Footnotes
AUTHOR CONTRIBUTIONS
A.L.R. and Z.C.W. designed the experiment and supervised the project. Y.L. (Yang Li) performed the in vitro laboratory experiments including generating cDNA vector constructs, gene transfer and knockdown, RT-PCR and Western blot analysis, and drug sensitivity and localization studies. L.Z. performed the PAM data analysis and statistics. Q.L., under supervision of Z.S., performed bioinformatic analysis on validation data sets and the cisplatin trial data set. R.T. performed apoptosis assays. Y.L. (Yan Li) contributed to the preparation of cDNA vector constructs. C.D., C.S., and B.H.-K. performed the epirubicin trial and provided clinical data and analysis. A.L.R. provided the DFHCC and cisplatin trial clinical samples and performed gene expression array analysis. Z.C.W. performed the SNP array analysis and scored the FISH assays. A.L.R., J.D.I. and Z.C.W. wrote the manuscript with comments from Y.L.(Yang Li), C.D., C.S., and Z.S.
ACKNOWLEDGMENTS
We thank Charles Lee, Abha Aggarwal (DF/HCC cytogenetics core facility); Edward Fox, Pamela Hollasch, Maura Berkeley (DF/HCC microarray core facility); Rebecca Gelman, Weixiu Luo, Xin Lu (Biostatisticians); and members of the Richardson-Wang lab for their advice and assistance. We also thank Jean Zhao for kindly providing vectors and cell lines, Aron C. Eklund, Nicolai Juul, and Rou-Li Zhou for helpful discussions and advice, and Daniel Silver and Ji-Young Kim for their critical review of this manuscript. Human mammary epithelial cells (HMECs) carrying a dominant negative allele of p53 were generously provided by Dr. Jean Zhao, Dana-Farber Cancer Institute, Boston, MA. This work was supported by the Breast Cancer Research Foundation (BCRF) in New York. Support also came from the NCI SPORE in Breast Cancer at Harvard (CA89393) and a DOD Concept Award (BC053041). The Trial of Principle (TOP) trial was supported by the Fondation Luxembourgeoise contre le Cancer, by the Fonds National de la Recherche Scientifique (CS, BH, CD), by the Brussels Region, as well as by the European Commission through the “Advancing Clinico-Genomic Trials” project (FP6-2005-IST-026996).
References
- 1. E.B.C.T.C.G.Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trialsLancet200536516871717[PubMed][Google Scholar]
- 2. International Web-based consultation on priorities for translational breast cancer researchBreast Cancer Res20079R81[PubMed][Google Scholar]
- 3. Multidrug resistance in cancer: role of ATP-dependent transportersNat Rev Cancer200224858[PubMed][Google Scholar]
- 4. Gene expression and amplification in breast carcinoma cells with intrinsic and acquired doxorubicin resistanceOncogene20012013001306[PubMed][Google Scholar]
- 5. Genomic signatures to guide the use of chemotherapeuticsNat Med20061212941300[PubMed][Google Scholar]
- 6. Prediction of doxorubicin sensitivity in breast tumors based on gene expression profiles of drug-resistant cell lines correlates with patient survivalOncogene20052475427551[PubMed][Google Scholar]
- 7. A strategy for predicting the chemosensitivity of human cancers and its application to drug discoveryProc Natl Acad Sci U S A20071041308613091[PubMed][Google Scholar]
- 8. Clinical evaluation of chemotherapy response predictors developed from breast cancer cell linesBreast Cancer Res Treat2009[Google Scholar]
- 9. Molecular genetics of breast cancer progressionSemin Cancer Biol19999277288[PubMed][Google Scholar]
- 10. Profiling breast cancer by array CGHBreast Cancer Res Treat200378289298[PubMed][Google Scholar]
- 11. Chromosome aberrations in solid tumorsNat Genet200334369376[PubMed][Google Scholar]
- 12. Genomic and transcriptional aberrations linked to breast cancer pathophysiologiesCancer Cell200610529541[PubMed][Google Scholar]
- 13. Diagnosis of multiple cancer types by shrunken centroids of gene expressionProc Natl Acad Sci U S A20029965676572[PubMed][Google Scholar]
- 14. Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancerLancet2005365671679[PubMed][Google Scholar]
- 15. Gene expression profiling predicts clinical outcome of breast cancerNature2002415484485[PubMed][Google Scholar]
- 16. MTDH activation by 8q22 genomic gain promotes chemoresistance and metastasis of poor-prognosis breast cancerCancer Cell200915920[PubMed][Google Scholar]
- 17. A gene-expression signature as a predictor of survival in breast cancerN Engl J Med200234719992009[PubMed][Google Scholar]
- 18. Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancerCancer Res2006661029210301[PubMed][Google Scholar]
- 19. Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosisJ Natl Cancer Inst200698262272[PubMed][Google Scholar]
- 20. Gene expression profiling spares early breast cancer patients from adjuvant therapy: derived and validated in two population-based cohortsBreast Cancer Res20057R953R964[PubMed][Google Scholar]
- 21. Cellular functions of 14-3-3zeta in apoptosis and cell adhesion emphasize its oncogenic characterOncogene2007[Google Scholar]
- 22. 14-3-3zeta overexpression defines high risk for breast cancer recurrence and promotes cancer cell survivalCancer Res20096934253432[PubMed][Google Scholar]
- 23. The oncogenic properties of mutant p110alpha and p110beta phosphatidylinositol 3-kinases in human mammary epithelial cellsProc Natl Acad Sci U S A20051021844318448[PubMed][Google Scholar]
- 24. Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicineClin Chem199339561577[PubMed][Google Scholar]
- 25. A stroma-related gene signature predicts resistance to neoadjuvant chemotherapy in breast cancerNat Med2009156874[PubMed][Google Scholar]
- 26. Efficacy of Neoadjuvant Cisplatin in Triple-Negative Breast CancerJ Clin Oncol(in press).[Google Scholar]
- 27. Gene expression profiling for the prediction of therapeutic response to docetaxel in patients with breast cancerLancet2003362362369[PubMed][Google Scholar]
- 28. A mammalian lysosomal membrane protein confers multidrug resistance upon expression in Saccharomyces cerevisiaeJ Biol Chem19992741287712882[PubMed][Google Scholar]
- 29. BRCA1 promoter methylation in sporadic breast tumors: relationship to gene expression profilesBreast Cancer Res Treat200591179186[PubMed][Google Scholar]
- 30. Predicting features of breast cancer with gene expression patternsBreast Cancer Res Treat2008108191201[PubMed][Google Scholar]
- 31. X chromosomal abnormalities in basal-like human breast cancerCancer Cell20069121132[PubMed][Google Scholar]
- 32. Loss of heterozygosity and its correlation with expression profiles in subclasses of invasive breast cancersCancer Research2004646471[PubMed][Google Scholar]
- 33. Biological processes associated with breast cancer clinical outcome depend on the molecular subtypesClin Cancer Res20081451585165[PubMed][Google Scholar]
- 34. Evaluating the yield of medical testsJama198224725432546[PubMed][Google Scholar]
- 35. Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimationStat Med20042321092123[PubMed][Google Scholar]
- 36. Pooling breast cancer datasets has a synergetic effect on classification performance and improves signature stabilityBMC Genomics20089375[PubMed][Google Scholar]