PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of oncotargetLink to Publisher's site
 
Oncotarget. 2016 November 29; 7(48): 79494–79502.
Published online 2016 November 4. doi:  10.18632/oncotarget.13104
PMCID: PMC5346730

Identification of KLK10 as a therapeutic target to reverse trastuzumab resistance in breast cancer

Abstract

Trastuzumab, the first antibody widely used in anti-HER2 targeted therapy, dramatically improved the overall outcome of HER2 positive breast cancer patients. However, trastuzumab resistance emerged as a major problem in its clinical application. In order to explore mechanisms underlying trastuzumab resistance, we performed RNA-Seq to analyze the gene expression variation in trastuzumab resistant breast cancer cell line. The sequencing result was then combined with the relevant data in TCGA database to conduct a co-expression analysis. We found a series of differentially expressed genes with potential contributions to trastuzumab resistance. Among them, KLK10 was verified to be a potential therapeutic target for reversing trastuzumab resistance. In summary, this study provides a new clue to screen molecular targets and predictive biomarkers for trastuzumab resistance.

Keywords: trastuzumab resistance, KLK10, breast cancer, HER2

INTRODUCTION

Breast cancer accounts for the highest morbidity and mortality among all cancers in female globally [1]. Nearly 30% breast cancer patients have HER2 gene amplification [2]. Trastuzumab, the first antibody widely used in anti-HER2 targeted therapy, dramatically improved the overall survival of HER2 positive breast cancer patients. Recently, it has been approved to treat HER2 positive gastric cancer [3]. With its extensive application, trastuzumab resistance emerged as a major problem. Novel targets are expected to reverse trastuzumab resistance. Unfortunately, no effective targets or biomarkers have been approved for trastuzumab resistance. Most of such efforts to identify biomarkers or targets for trastuzumab resistance were based on the molecular mechanism of trastuzumab [4, 5]. More practical alternative approaches would be necessary to identify biomarkers to predict and targets to reverse trastuzumab resistance.

RNA-Seq(RNA sequencing) technology has been commonly used in high throughput analysis of genome-wide gene expression [6]. In addition, the Cancer Genome Atlas (TCGA) project collects high throughput analyses such as gene expression profiling, exon sequencing, SNP genotyping, genomic DNA methylation profiling and microRNA profiling along with clinical data of each patient [7]. In this study, we are trying to combine our RNA-Seq analysis of trastuzumab resistant breast cancer cells with TCGA database to discover potential biomarkers and therapeutic targets for trastuzumab resistance in breast cancer.

RESULTS

Establishment of trastuzumab resistant breast cancer cell line

BT474 HR (Herceptin Resistant) cells were established by culturing BT474 cells with 1μg/ml Trastuzumab for 6 months and 4 μg/ml Trastuzumab for 3 months. No obvious cellular morphology changes were observed in BT474 and BT474 HR. As expected, trastuzumab could remarkably inhibit the growth of BT474 but not BT474 HR cells (Figure (Figure1A).1A). To determine why trastuzumab can inhibit cell growth, cell apoptosis and cell cycle distribution were determined after trastuzumab treatment. Significant changes were observed in the distribution of cell cycle phases in BT474 after trastuzumab treatment (Figure (Figure1B).1B). Trastuzumab could induce G1 phase arrest strikingly in BT474 cells in a dose-dependent manner, but not in BT474 HR. However, Trastuzumab failed to induce apoptosis in neither BT474 nor BT474 HR (Figure (Figure1C1C and and1D1D).

Figure 1
Establishment of Trastuzumab resistant breast cancer cell line

RNA expression profiling of BT474 and BT474HR cells

We used RNA-Seq to reveal changes of transcriptome in BT474 and BT474 HR cells. 65,677 differentially expressed transcripts from 16,170 genes were received through RNA-Seq analysis. Genes had a mean transcript variant number of about 3 (1–41) (Figure (Figure2A).2A). The volcano plot was used to observe for abnormal signals (Figure (Figure2B).2B). After filtering out these outliers, 246 genes were found to be statistical significance (p < 0.05). Next, quantitative real-time PCR was used to validate differentially expressed genes including MAP9, MET, SPNS2, TCEA3 and UGCG using highly and equally expressed GAPDH, ERBB2 and SQSTM1 as the control. For all these tested genes, the expression determined by quantitative real-time PCR was consistent with RNA-Seq results (Figure (Figure2C).2C). The representative transcripts for each protein coding gene, which had a higher level of expression, were selected for further analysis. The data was plotted with expression ratio vs. average expression (Figure (Figure2D),2D), and there was neither obvious skewed distribution nor abnormal signal after filtering. Finally, differential expression data of 12,228 transcripts was extracted as representatives of effective protein coding genes.

Figure 2
RNA expression profiling of BT474 HR cells

Co-expression analysis

To explore functions of differential genes systematically, gene co-expression network was utilized. In this method, we selected genes both meaningful in our RNA-Seq data and in expression profile from TCGA. A total of 9,913 genes were obtained in two data sets. In TCGA, 444 cases were in accordance with the co-expression analysis criteria. This data set was analyzed by WGCNA clustering and 36 gene sets were finally clustered. The clusters were then correlated with expression features in tumor tissues, ER, PR and HER2 states (Figure (Figure3A).3A). For summarizing such clusters, the principal component of each cluster or module eigengene (ME) was used. For instance, ME0 had no significant correlations with all features, while HER2 status had no significant correlation to any clusters but ME32. Different cluster had various degrees of relevance to tissue types, ER and PR. Highly similar correlation patterns of ER and PR implied the clustering of co-expression was a good indicator for biological functions.

Figure 3
Co-expression analysis of RNA-Seq and TCGA database

If drug resistance-related genes were irrelevant to co-expression cluster genes, selected genes that changed most remarkably in the expression should be uniformly distributed in the co-expression cluster gene sets. In contrast, the relationship between this gene set and drug resistance was significant when a particularly large number of differentially expressed genes were presented in some co-expression gene sets. Therefore, the top 10% differentially expressed genes were selected, and the distributions of their frequency of occurrence in the co-expression gene cluster sets were compared and statistically tested to show whether they consisted more than 10% of a gene set. As shown in Figure Figure3B,3B, ME3 and ME6 gene sets had more top 10% differentially expressed genes. It implied that these gene sets were more significantly related with drug resistance. Also, they were more related to tumor, representing good sources for targets and biomarkers identification.

Target validation

Therefore, KLK10 from ME3 and KLK11 from ME6 were selected as potential targets for further validations. Receptor tyrosine kinase encoding EPHA3 from ME4, which had a low score, was chosen as a control.

Quantitative real-time PCR validated the differential expression of these genes (Figure (Figure4A).4A). To further explore the biological relevance of these genes to drug resistance, trastuzumab induced growth inhibition before and after knock-down of these potential targets were determined by MTS assay. While depletion of EPHA3 or KLK11 had no significant effects on trastuzumab sensitivity (p > 0.05), KLK10 knock-down significantly reversed resistance to trastuzumab (p < 0.05) (Figure (Figure4B4B and and4C4C).

Figure 4
Functional validation of potential targets

To further confirm the relevance of KLK10 to trastuzumab resistance, cell cycle distribution was measured in the presence of trastuzumab and KLK10 depletion. After KLK10 depletion, trastuzumab successfully induced G1 arrest (Figure (Figure5A).5A). As dephosphorylated Retinoblastoma (RB) protein is a well-known marker for G1 phase [8], we detected phosphorylation level of RB in BT474 HR cells treated with trastuzumab and KLK10 siRNA. pRB in BT474 but not BT474 HR cells were indeed decreased after trastuzumab treatment (Figure (Figure5B).5B). Once KLK10 was depleted, trastuzumab successfully decreased pRB in BT474 HR cells (Figure (Figure5B5B and and5C).5C). Trastuzumab treatment led to almost 60% inhibition of pRB in BT474 cells, but only 17% in BT474 HR cells (Figure (Figure5D).5D). However, KLK10 siRNA significantly increased trastuzumab-induced inhibition of pRB. Furthermore, SFIT method [9] was used to distinguish the cell cycle phases. After knock-down of KLK10, trastuzumab-treated cells were accumulated in G1 phase with low pRB levels (Figure (Figure5E).5E). Together, KLK10 siRNA succeeded to reverse trastuzumab resistance.

Figure 5
KLK10 is a potential target to reverse trastuzumab resistance

Given the significant effect of drug resistance on clinical outcome, we explored the relevance of KLK10 expression to the overall survival of breast cancer patients. According to the median expression of KLK10, 434 patients from TCGA database were divided into high expression group and low expression group (Figure (Figure5F).5F). High expression of KLK10 was related to poor prognosis in HER2 positive breast cancer patients (n = 434, log-rank p = 0.0165418). Cox proportional hazards regression analysis showed that high expression of KLK10 led to an elevation of 18.1% hazard risk (95% CI:1.015–1.375, p = 0.031).

DISCUSSION

Despite the successful application of trastuzumab for the treatment of HER2 positive breast cancer, its acquired or intrinsic resistance hurdled improvement of breast cancer patients. In this study, we profiled gene expression in a newly established trastuzumab resistant breast cancer model and combined it with publically available database to identify potential targets to reverse trastuzumab resistance. Indeed, KLK10 was increased in trastuzumab-resistant cells and its depletion succeeded to reverse trastuzumab resistance. Importantly, high KLK10 expression was associated with poor prognosis of breast cancer, indicating that KLK10 is a potential target and prognosis predictor for trastuzumab resistance.

Although there are 246 genes detected to be statistically significant differentially expressed between BT474 and BT474 HR cells, some may have no relevance to trastuzumab resistance per se. Therefore, other approaches are needed to screen potential targets for further functional analyses. Co-expression is models of multiple gene expression pattern associated with certain potential biological functions. Co-expression can determine the functional expression pattern of breast tissue while differential expression analysis can reflect the single gene expression patterns in two cell straits, with both random patterns and meaningful patterns related to trastuzumab resistance. By constructing a reasonable score, a combination of both methods can be achieved to get a list of genes with important functions. By enrichment analysis of clustering gene sets and differential expression data, we obtained ME3, ME6 and other gene sets potentially related with trastuzumab resistance. The persistent differentially expressed genes in such gene sets would have much higher chance to be prediction biomarkers and intervention targets for trastuzumab resistance. Indeed, we have successfully identified KLK10 from ME3 as a relevant target to reverse trastuzumab resistance.

The members of KLK family are exocrine protein readily to be detected in the serum. Interestingly, KLK3 is a tumor-associated biomarker well known as prostate specific antigen (PSA). KLK4, KLK5, KLK6 and KLK7 have been found related to the prognosis of ovarian cancer [1013]. Therefore, it would be interesting to explore whether plasma level of KLK10 are associated with drug resistance and overall survival in breast cancer. In fact, KLK10 also has been found related to the prognosis of breast cancer due to its association with tamoxifen resistance [14]. In summary, by co-expression analysis of TCGA data with gene expression profile of trastuzumab resistant breast cancer, we identified KLK10 as a potential biomarker and intervention target for trastuzumab resistance. We found a series of differentially expressed genes with potential contributions to trastuzumab resistance. Among them, KLK10 was verified to be a potential therapeutic target for reversing trastuzumab resistance.

MATERIALS AND METHODS

Cell lines

Human breast cancer cell line BT474 was obtained from American Type Culture Collection (ATCC, Manassas, VA, USA) and maintained in RPIM 1640 supplemented with 10% FBS at 37°C in a humidified, 5% CO2 incubator. BT474 HR (Herceptin Resistant) cells were maintained in the presence of 10 μg/ml trastuzumab.

Cell viability assay

MTS assay was performed with the CellTiter 96® Aqueous No-Radioactive Cell Proliferation Assay Kit (Promega). The cells were transferred into a 96-well plate and cultured overnight before adding trastuzumab or phosphate buffer saline (PBS). 48 h or 72 h later, the cell viability was measured following the manufacturer's instruction. Samples were prepared in triplicates at least.

RNA extraction and quantitative real-time RT-PCR

Total RNA was extracted by TRIzol reagent (Invitrogen) according to the manufacturer's. RNA concentrations were quantified by NanoDrop™ 2000 (Thermo Scientific). Reverse transcription reaction was performed using 2 μg of total RNA with Quantscript RT Kit (Tiangen biotechnology, Beijing, China). The mRNA expression level was determined by quantitative real-time PCR using Bestar® SybrGreen qPCR mastermix (DBI) and LightCycler 480® II Real-Time PCR System (Roche). Primers used are listed in Table Table11.

Table 1
Primers used in the study

RNA-seq

RNA-Seq was performed with Ion Total RNA-Seq kit v2 of Ion Proton™ (Ion Torrent, Life Technologies). GRCh37 reference genome from the phase 3 of the 1000 Genomes Project was used for RNA-Seq alignment. Gene annotation of GRCh37.p13 GENCODE Release 19 was utilized to determine the splicing site annotation files. STAR v2.4.1d was used [15] and 2-pass strategy preliminary comparison was adopted, annotations splicing site information was utilized, and the parameters were same in ENCODE Projec [16]. To align the unsuccessful sequence of STAR, Bowtie2 v2.2.4 was used with a more sensitive parameter “—local —very-sensitive-local” [17]. Finally, Samtools v1.1 was utilized to combine the results of the two methods.

Obtain the differential expression by annotation database

The numbers of transcript fragments per kilobase of per million mapped reads (FPKM) were calculated to determine the relative amount of mRNA in the cells. Data was normalized with the default FPKM method of cuffdiff v2.2.1 from Cufflink [18]. GRCh37.p13 GENCODE Release 19 was selected for annotation and cummeRbund on the platform R v3.2.4 was used to determine the differential expression according to the comparison of BT474 and BT474 HR data.

TCGA database extraction

All data which included both Level 3 microarray data and clinical data in XML format data set of breast cancer patients by TCGA DATA Portal were selected. XML package of R was used to parse clinical data in XML format. The required information was extracted and merged into the corresponding microarray data. Microarray data downloaded from UNC AgilentG4502 platform including 593 cases of specimens. These specimens were from 529 breast cancer patients' tumor and normal tissues. The patients' clinical profiles including the age, gender, race, follow-up times (days), end event, the method of first time confirmed diagnosis, histological type, ER and PR status, immunohistochemistry and FISH of HER2, pathological stage and grade were extracted.

Co-expression analysis

WGCNA v1.46 of R was used to cluster the dataset [19]. Correlation analyses of other clinical indicators were adopted to identify potential tumor-associated co-expression patterns. For each gene co-expression cluster, the differential expression result was evaluated with enrichment analysis by binomal test and calculated a score for the association of drug resistance.

siRNA transfection

KLK10, KLK11 and EPHA3 depletion were achieved by transfection with siRNA (Gene pharma, Shanghai, China). Cells were seeded overnight in 6-well plates (4 × 105/well) and transfected with siRNA duplexes (20 nM) using LipofectamineTM RNAiMAX transfection reagent (Invitrogen) following the manufacturer's instruction. The sequences of various siRNAs are listed in Table Table22.

Table 2
siRNAs used in the study

Flow cytometry analysis

4 × 105 cells cultured overnight in 6-well plates were treated with or without trastuzumab and harvested after 48 or 72 hours. Cell apoptosis were detected with apoptosis kit (FITC Annexin V Apoptosis Detection Kit I, BD PharmingenTM). Briefly, cells were washed twice with cold PBS and incubated in 100 μl binding buffer with 5 μl of FITC Annexin V and 5 μl PI for 15 minutes in the dark. For cell cycle analysis, cells were resuspended in 200 μl cold PBS, and then added into 1 ml 70% ethanol. After an hour, cells were transferred to 450 μl PBS with 40 μl RNAse (Sigma) and 10 μl PI. For pRB expression analysis [20], room-temperature 1.5% (vol/vol) paraformaldehyde was added for 10 minutes to fix the cells washed twice with cold PBS. Fixed cells were permeabilized by slowly adding cold 100% methanol. 100 μl PBS with 1 μg pRB antibody (Phospho-RB Ser807/811, Cell Signaling Technology) was then added to mark the cells. After the incubation of fluorochrome-conjugated secondary antibody, cells were then stained with PI as previously described.

Survival analysis

434 breast cancer (invasive ductal carcinoma) patients with follow-up information from TCGA DATA Portal were chosen. According to the median of KLK10 expression, patients were divided into high and low expression group, Kaplan-Meier Plot was used to reveal differences in survival between the two groups. The log-rank test was used for the statistical analysis of overall survival. A cox proportional hazards regression analysis was conducted to quantify the risk.

Footnotes

CONFLICTS OF INTEREST

Authors declare no conflicts of interest.

GRANT SUPPORT

This work was supported by High Level Talents Program from the Department of Health, 151 talents program in Zhejiang, Natural Science foundation of China (81572715), and the Department of Health in Zhejiang Province (WKJ-ZJ-1520).

REFERENCES

1. Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108. doi: 10.3322/caac.21262. [PubMed] [Cross Ref]
2. Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL. Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987;235:177–82. doi. [PubMed]
3. Yarden Y, Pines G. The ERBB network: at last, cancer therapy meets systems biology. Nat Rev Cancer. 2012;12:553–63. doi: 10.1038/nrc3309. [PubMed] [Cross Ref]
4. Minuti G, Cappuzzo F, Duchnowska R, Jassem J, Fabi A, O'Brien T, Mendoza AD, Landi L, Biernat W, Czartoryska-Arlukowicz B, Jankowski T, Zuziak D, Zok J, et al. Increased MET, HGF gene copy numbers are associated with trastuzumab failure in HER2-positive metastatic breast cancer. Br J Cancer. 2012;107:793–9. doi: 10.1038/bjc.2012.335. [PMC free article] [PubMed] [Cross Ref]
5. Lipton A, Goodman L, Leitzel K, Cook J, Sperinde J, Haddad M, Kostler WJ, Huang W, Weidler JM, Ali S, Newton A, Fuchs EM, Paquet A, et al. HER3, p95HER2, and HER2 protein expression levels define multiple subtypes of HER2-positive metastatic breast cancer. Breast Cancer Res Treat. 2013;141:43–53. doi: 10.1007/s10549-013-2665-0. [PMC free article] [PubMed] [Cross Ref]
6. Chu Y, Corey DR. RNA sequencing: platform selection, experimental design, and data interpretation. Nucleic Acid Ther. 2012;22:271–4. doi: 10.1089/nat.2012.0367. [PMC free article] [PubMed] [Cross Ref]
7. Cancer Genome Atlas N. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70. doi: 10.1038/nature11412. [PMC free article] [PubMed] [Cross Ref]
8. Goodrich DW, Wang NP, Qian Y-W, Eva Y-HL, Lee W-H. The retinoblastoma gene product regulates progression through the G1 phase of the cell cycle. Cell. 1991;67:293–302. doi. [PubMed]
9. Dean PN. A simplified method of DNA distribution analysis. Cell Tissue Kinet. 1980;13:299–308. doi. [PubMed]
10. Obiezu CV, Scorilas A, Katsaros D, Massobrio M, Yousef GM, Fracchioli S, Rigault de la Longrais IA, Arisio R, Diamandis EP. Higher human kallikrein gene 4 (KLK4) expression indicates poor prognosis of ovarian cancer patients. Clin Cancer Res. 2001;7:2380–6. doi. [PubMed]
11. Kim H, Scorilas A, Katsaros D, Yousef GM, Massobrio M, Fracchioli S, Piccinno R, Gordini G, Diamandis EP. Human kallikrein gene 5 (KLK5) expression is an indicator of poor prognosis in ovarian cancer. Br J Cancer. 2001;84:643–50. doi: 10.1054/bjoc.2000.1649. [PMC free article] [PubMed] [Cross Ref]
12. Diamandis EP, Yousef GM, Soosaipillai AR, Bunting P. Human kallikrein 6 (zyme/protease M/neurosin): a new serum biomarker of ovarian carcinoma. Clin Biochem. 2000;33:579–83. doi. [PubMed]
13. Psyrri A, Kountourakis P, Scorilas A, Markakis S, Camp R, Kowalski D, Diamandis EP, Dimopoulos MA. Human tissue kallikrein 7, a novel biomarker for advanced ovarian carcinoma using a novel in situ quantitative method of protein expression. Ann Oncol. 2008;19:1271–7. doi: 10.1093/annonc/mdn035. [PubMed] [Cross Ref]
14. Luo LY, Diamandis EP, Look MP, Soosaipillai AP, Foekens JA. Higher expression of human kallikrein 10 in breast cancer tissue predicts tamoxifen resistance. Br J Cancer. 2002;86:1790–6. doi: 10.1038/sj.bjc.6600323. [PMC free article] [PubMed] [Cross Ref]
15. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, Batut P, Chaisson M, Gingeras TR. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21. doi: 10.1093/bioinformatics/bts635. [PMC free article] [PubMed] [Cross Ref]
17. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. doi: 10.1038/nmeth.1923. [PMC free article] [PubMed] [Cross Ref]
18. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5. doi: 10.1038/nbt.1621. [PMC free article] [PubMed] [Cross Ref]
19. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC bioinformatics. 2008;9:1. doi. [PMC free article] [PubMed]
20. Rosner M, Schipany K, Hengstschlager M. Merging high-quality biochemical fractionation with a refined flow cytometry approach to monitor nucleocytoplasmic protein expression throughout the unperturbed mammalian cell cycle. Nat Protoc. 2013;8:602–26. doi: 10.1038/nprot.2013.011. [PubMed] [Cross Ref]

Articles from Oncotarget are provided here courtesy of Impact Journals, LLC