OR Oncology Reports 1021-335X 1791-2431 D.A. Spandidos 10.3892/or.2018.6376 or-39-06-2978 Articles HSPB8 promotes cancer cell growth by activating the ERK-CREB pathway and is indicative of a poor prognosis in gastric cancer patients ShenJian 1 LiMinzhe 1 MinLi 2 Department of Surgery, Beijing Chaoyang Hospital, Capital Medical University, Beijing 100020, P.R. China Department of Gastroenterology, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory for Precancerous Lesion of Digestive Diseases, National Clinical Research Center for Digestive Diseases, Beijing 100050, P.R. China Correspondence to: Dr Li Min, Department of Gastroenterology, Beijing Friendship Hospital, Capital Medical University, 95 Yong'an Road, Xicheng, Beijing 100050, P.R. China, E-mail: minli@ccmu.edu.cn Dr Jian Shen, Department of Surgery, Beijing Chaoyang Hospital, Capital Medical University, 5 Gongtinan Road, Chaoyang, Beijing 100050, P.R. China, E-mail: 410078@163.com 062018 18042018 39 6 2978 2986 26102017 11042018 Copyright © 2018, Spandidos Publications 2018

Gastric cancer (GC) is one of the most commonly diagnosed malignancies worldwide, especially in East Asia. Discovery of new biomarker and the elucidation of the molecular mechanisms involved in GC development and progression continue to be important issues for both researchers and clinicians. In the present study, we report that siRNA knockdown of heat shock protein family B (small) member 8 (HSPB8) inhibited the proliferation of GC cells and promoted their apoptosis. Analysis of TCGA dataset indicated that the HSPB8 expression level was strongly positively correlated with the KEGG MAPK signaling pathway (P<0.001, FDR=0.006) and BIOCARTA CREB pathway (P=0.006, FDR=0.043). The association between HSPB8 and the ERK-CREB pathway was confirmed by western blot analysis and we found that pERK and pCREB were significantly decreased following HSPB8 knockdown. Downstream genes of the ERK-CREB pathway were all significantly decreased following HSPB8 knockdown. By evaluating the survival of TCGA GC patients, we found that patients with a high HSPB8 level exhibited significantly worse prognosis than those with low HSPB8 in both overall survival (OS) (log-rank χ2=10.60, P=0.001) and disease-free survival (DFS) (log-rank χ2=11.31, P<0.001). The methylation level of HSPB8 DNA was significantly negatively associated with its expression (R=−0.1368, P=0.041), and positively associated with OS (log-rank χ2=10.60, P=0.001). In conclusion, we provide evidence that HSPB8 promotes cancer cell growth by activating the ERK-CREB pathway and may serve as a potential prognostic factor in GC patients.

 HSPB8 gastric cancer cell growth apoptosis prognosis Introduction

Gastric cancer (GC) is one of the most aggressive cancer worldwide, with approximately 951,600 new cases and 723,100 new deaths in 2012 (1). Even though the global GC incidence is decreasing, the overall disease burden still ranks 5th in incidence and 3rd in mortality (1). In particular, 40% of all GC cases occur in China and other East Asian countries, which may be associated with various unhealthy living habits (2). Recently, endoscopy examination is the main approach for GC screening, but the high medical cost of endoscopy and the shortage of endoscopic professionals in primary medical centers largely hampers the early diagnosis of GC (3,4). Advances in the identification of biomarkers for detecting precancerous and cancerous gastric lesions have offered alternative strategies for GC screening (5). For example, serum levels of pepsinogen I and II have been used in the identification of high GC risk individuals before endoscopic examination in numerous studies (5,6).

GC is considered as a disease with high heterogeneity, yet the majority of GC patients are treated with similar chemotherapeutic drugs and surgical techniques, resulting in unfavorable combined sequelae and side effects (7,8). Biomarkers also hold high expectation in the preoperative classification of GC. Actually, protein expression levels, non-coding RNAs, gene copy variations, and single nucleotide polymorphisms (SNPs) are all potential biomarkers of GC (9–11). For example, high expression of serpin A1 was found to be an indicator of a poor prognosis (12), and rs629367 was identified as correlated with poor survival in GC patients (13). However, apart from the well-known HER2 amplification, there are few biomarkers that can be utilized in GC patient classification (14). It is still urgent to discover new biomarkers to facilitate clinical decision-making and avoid unnecessary over-treatment.

Heat shock protein family B (small) member 8 (HSPB8) is a member of the small heat shock protein superfamily, which contains a conservative α-crystallin domain at the C-terminal (15–17). The most well-known function of HSPB8 is acting as a chaperone in association with Bag3 in the regulation of macroautophagy (15,18–20). HSPB8 was found to be associated with many diseases, such as cardiomyopathy (21), amyotrophic lateral sclerosis (22) and Alzheimer's disease (23). HSPB8 was also found to be associated with estrogen-related cancers. Piccolella et al found that HSPB8 modulates the proliferation and migration of breast cancer cells (24). Suzuki et al found that HSPB8 regulates TGF-α-induced ovarian cancer cell migration (25). However, there are few reports concerning the role of HSPB8 in gastrointestinal cancers. Based on our previous study of the molecular signature of GC subtypes, HSPB8 was found to be involved in both diffuse and intestinal GCs (26). In the present study, we conducted a series of further assays in vitro and in silico to reveal the biological role and potential prognostic value of HSPB8 in GC.

 Materials and methods <sec> <title>Cell line culture and siRNA transfection

Human gastric cancer cell lines (AGS, BGC-803, BGC-823, SCG-7901 and N87) and normal gastric cell line GES-1 were provided by the Cancer Institute and Hospital, Chinese Academy of Medical Sciences (Beijing, China) and maintained by our laboratory. All cells were maintained in Dulbecco's modified Eagle's medium (DMEM) with 10% fetal bovine serum (FBS) (purchased from Gibco (Grand Island, NJ, USA), in an incubator containing 5% CO2 at 37°C. A total of 4×105 cells were seeded in each well of 6-well plates. After culturing for 24 h, cells were transfected with the siRNAs using Lipofectamine 2000 (Life Technologies; Thermo Fisher Scientific, Inc., Waltham, MA, USA) according to the manufacturer's protocol. All siRNAs were de novo synthesized by Gene Pharma Co. (Shanghai, China) and all the sequences are listed in Table I.

 Western blot analysis and antibodies

Rabbit antibodies against HSPB8 (1:500 diluted; cat. no. ab96837), CREB (1:1,000 diluted; cat. no. ab31387) and pCREB (Ser113, 1:1,000 diluted; cat. no. ab32096) were purchased from Abcam (Cambridge, MA, USA). Rabbit antibodies against GAPDH (1:10,000 diluted; cat. no. 14C10), Erk (1:2,000 diluted; cat. no. 137F5) and p-Erk1/2 (Thr202/Tyr204; 1:1,000 diluted; cat. no. D13.14.4E) were purchased from Cell Signaling Technology, Inc. (Berkley, MA, USA). Proteins were extracted using lysis buffer (50 mM Tris-HCl, pH 7.4; 10 mM EDTA; 0.5% NP-40; 1% Triton X-100) with a protease inhibitor cocktail (Sigma-Aldrich; Merck KGaA, Darmstadt, Germany) added. SDS-polyacrylamide gel electrophoresis (PAGE) was performed and then proteins were transferred from gels to PVDF membranes. The membranes were incubated with primary antibodies overnight at 4°C, washed with TBST for 5 times, 5 min each time, and then incubated for 45 min with HRP-labeled goat anti-rabbit IgG antibodies (1:2,000 diluted; cat. no. ZDR-5306; ZSGB-Bio Co., Beijing, China) at room temperature, washed with TBST for 3 times, 5 min each time. All proteins were detected with an ECL Plus system (Beyotime Institute of Biotechnology, Jiangsu, China).

 Cell proliferation and apoptosis assays

Cells were seeded in 96-well plates at a density of 2×103 cells per well. The viability of cells was determined using Cell Counting Kit-8 (CCK-8; Dojindo Molecular Technologies, Inc., Kumamoto, Japan). The optical density at 490 nm was measured by Smartspec Model 450 (Bio-Rad Laboratories, Inc., Hercules, CA, USA) after 0, 24, 48 and 72 h, in triplicate. The apoptosis rate of the cells was measured by flow cytometry. Cells were obtained and washed with pre-cooling PBS, and re-suspended in binding buffer. Annexin V-FITC (5 µl) and 7AAD (5 µl) reagents were added to each sample and incubated at 25°C for 15 min away from light. An additional 400 µl binding buffer was added and then the cells were analyzed by a flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA). All assays were repeated three times independently.

 Reverse transcription and real-time quantitative PCR

Total RNA was extracted using TRIzol reagents and reverse transcripted using SuperScript II reverse transcriptase (all from Life Technologies; Thermo Fisher Scientific, Inc., Waltham, MA, USA) according to the manufacturer's protocol. Real-time quantitative PCR (qPCR) was conducted in triplicates using Applied Biosystems 7500 using the SYBR-Green PCR Master Mix (Roche). Relative mRNA levels of each gene were normalized to GAPDH. All primer sequences used in this study are listed in Table II.

 Online patient data acquisition and statistical analysis

Phenotype data with genetic information of the patients were all downloaded from TCGA (http://cancergenome.nih.gov/) (27,28). For Gene Set Enrichment Analysis (GSEA) (29), a Pearson correlation based method (1000 permutations of phenotype were run as recommended) was applied to evaluate the correlation between HSPB8 expression and gene sets of Kyoto Encyclopedia of Genes and Genomes (KEGG) (30) and BIOCARTA (31). A FDR q-value <0.25 was considered significant. For survival analysis, the median level of mRNA/methylation of HSPB8 was chosen as the cut-off to separate two subgroups and the Log-rank test was applied and corresponding Kaplan-Meier plots were drawn to show the results intuitionally. Correlation between mRNA and methylation of HSPB8 was calculated by Pearson Chi-square analysis. For cell proliferation, apoptosis and gene expression analysis, unpaired Student's t-tests were used to compare two groups when the means follow the normal distribution, and P<0.05 was considered statistically significant.

 Results <sec> <title>Knockdown of HSPB8 inhibits gastric cancer cell proliferation

To evaluate the expression of HSPB8 in GC cells and normal gastric epithelial cells, we detected HSPB8 in five GC cell lines (AGS, BGC-803, BGC-823, SGC-7901 and N87) and one gastric epithelial cell line (GES-1). The results demonstrated that HSPB8 was relatively higher in the GC cells than that noted in the normal epithelial cells (Fig. 1A). To examine the proliferation-promoting potential of HSPB8, two siRNAs targeting HSPB8 were designed (sequences are shown in Table I) and transfected into AGS and BGC-803 cells. The HSPB8 protein levels were decreased by both siRNAs successfully (Fig. 1B and C), and the proliferation rates of both AGS (Fig. 1D and F) and BGC-803 (Fig. 1E and G) cell lines were significantly hampered by HSPB8 knockdown. Considering that the efficiency of siRNA1 was higher than that of siRNA2, the subsequent analyses were all performed using siRNA1.

 Knockdown of HSPB8 promotes gastric cancer cell apoptosis

Suppression of cancer cell growth is usually associated with activation of cellular apoptosis. To examine the influence on apoptosis of HSPB8 knockdown, AGS and BGC803 cells transfected with HSPB8-siRNA1 and NC were subjected to Annexin V-FITC/7AAD staining and flow cytometry. The results indicated that in both AGS (Fig. 2A-C) and BGC-803 (Fig. 2D-F) cell lines, HSPB8 knockdown had no significant influence on early apoptosis but strongly increased the percentage of late apoptotic cells.

 HSPB8 is positively correlated with the ERK-CREB pathway

To further reveal the potential mechanism attributed to the oncogenic role of HSPB8 in gastric cancer, the TCGA dataset of gastric cancer expression data was downloaded and analyzed. The GSEA analysis revealed that the HSPB8 expression level was strongly positively correlated with the KEGG MAPK signaling pathway (NES=2.042, P<0.001, FDR=0.006, Fig. 3A) and the BIOCARTA CREB pathway (NES=1.776, P=0.006, FDR=0.043, Fig. 3B). Activation of the ERK-CREB pathway would be a possible key mechanism of the oncogenic activity of HSPB8.

 Knockdown of HSPB8 decreases the activity of the ERK-CREB pathway

To verify our hypothesis of HSPB8's association with the ERK-CREB pathway, we detected levels of phosphorylated and total ERK and CREB proteins under siRNA interference of HSPB8 by western blot analysis. Our results indicated that the total CREB level was stable, while pErk and pCREB were significantly decreased under HSPB8 knockdown (Fig. 4A and B). To reveal whether the changes in pErk and pCREB were time-dependent, we also detected the protein levels of these proteins at 48, 72 and 96 h. The level of HSPB8 returned to a normal level at 72 h after siRNA knockdown, suggesting that HSPB8 is crucial for cell survival and maintains a high synthesis rate in these cells. pErk and pCREB remained suppressed at 96 h in AGS cells, while in BGC-803 cells the levels of pErk and pCREB were restored at 72 h (Fig. 4C and D). Downstream genes of the ERK-CREB pathway were also detected by RT-qPCR (primers are shown in Table II), and we found that all the targeted genes we detected were significantly decreased under HSPB8 knockdown (Fig. 4E).

 Expression and methylation levels of HSPB8 are potential prognostic factors

To explore whether HSPB8 could be used as a biomarker in gastric cancer, we further evaluated the survival of TCGA gastric cancer patients. When divided by the median expression level, the patients with a high HSPB8 level exhibited a significantly worse prognosis than those with low HSPB8 in both overall survival (OS) (log-rank χ2=10.60, P=0.001, Fig. 5A) and disease-free survival (DFS) (log-rank χ2=11.31, P<0.001, Fig. 5B).

DNA methylation data were also co-analyzed with expression data, and indicated that the methylation level of HSPB8 DNA was significantly negatively associated with its expression (R= −0.1368, P=0.041, Fig. 5C). Patients with high HSPB8 methylation level exhibited a significantly better prognosis than those with low methylation in regards to OS (log-rank χ2=10.60, P=0.001, Fig. 5D).

 Discussion

The small heat shock protein (sHsp) superfamily consists of a series of 15–30 kDa proteins with a common α-crystallin domain at the C-terminal. The most well-studied role of sHsp is acting as molecular chaperones (32,33), i.e. preventing the aggregation of enzymes under heat shock conditions and stabilizing the proteins. There are also numerous studies revealing that sHsp promote the functional refolding of proteins after urea denaturation in an ATP-independent manner (34). There are many members of the sHsp superfamily involved in GC, such as HSP70 (35), HSP110 (36) and HSP27 (37). However, HSPB8 has drawn little attention in this field.

In our previous study, we conducted a pilot investigation of expression-based GC biomarker screening and found that HSPB8 is a potential biomarker in both diffuse and intestinal GC (26). In the present study, our results indicated that HSPB8 is a proliferation-promoting protein in GC cell lines, which could also aid in the evasion of apoptosis by GC cells. The oncogenic role of HSPB8 is in accordance with most other sHsp members, such as HSP70 (35), HSP110 (36) and HSP60 (38).

To reveal the possible mechanism of HSPB8 in GC, we identified molecular pathways associated with HSPB8 expression by analyzing the TCGA GC dataset. From all KEGG and BIOCARTA gene sets, we found that the pathways of MAPK and CREB were positively correlated with HSPB8, which was also validated by western blot analysis in two GC cell lines. Furthermore, we found that the targeted genes of the ERK-CREB pathway (MMP9, EGR1, CCNA1, CCND1, CCN1 and BCL2) were significantly decreased under HSPB8 knockdown. The positive correlation between CCN1 and HSPB8 (39) and the association between BCL2 and HSPB8 (40), were both reported in other cell lines, which may partially verify our finding of HSPB8's role in the regulation of the ERK-CREB pathway.

It is worth noting that HSPB8 was restored to a normal level at 72 h after siRNA knockdown, suggested that HSPB8 is crucial for cell survival and maintains a high synthesis rate in these cells. However, pErk and pCREB levels were restored along with HSPB8 at 72 h in BGC-803 cells, but maintained a suppressed level in AGS cells until 96 h. Thus, we suggest that HSPB8 knockdown may have a phenotypic effect on pErk and pCREB, but this appears to be cell line-specific (i.e. this consistent change of HSPB8 and pERK/pCREB was observed for the BGC-803 cell line but not for the AGS cell line at different time-points). It is also possible that HSPB8 knockdown resulted in a more far-reaching effect on AGS cells (i.e. even when HSPB8 was restored quickly after siRNA expired at 72 and 96 h, pErk and pCREB remained suppressed). In all, any phenotypic effects observed after transfection at 72 h should be interpreted with caution, and further investigation of how HSPB8 promotes pErk and pCREB in a cell line-specific manner is warranted.

Additionally, we also explored the potential of HSPB8 as a prognostic biomarker in GC. Based on our results, both the expression level and the methylation level of HSPB8 could be used as unfavorable prognostic factors. Although the expression level and the methylation level of HSPB8 were significantly correlated, the expression factor still achieved a higher discrimination than the methylation factor in GC patients stratified based on survival data. However, further studies with a larger population are needed before developing reagents for clinical use. In our opinion, the effect of HSPB8 in GC is only partially dependent on its methylation level. Transcriptional factor level and activity, histone modification (methylation and acetylation) and chromatin folding are all possible regulators of HSPB8 expression. For TCGA data analysis, most GC patients had a very low HSPB8 level (Fig. 5C), indicating that the differential expression of HSPB8 was only partially determined by its methylation level.

In conclusion, we demonstrated that HSPB8 could promote the proliferation and inhibit the apoptosis of GC cells by activating ERK-CREB signaling. We also demonstrated that high expression of HSPB8 was indicative of a poor prognosis in GC patients. These findings suggest new evidence to improve our understanding of the mechanisms attributed to the carcinogenesis of GC, and provide insight into the development of patient stratification strategies for GC treatment.

 Acknowledgements

Not applicable.

 Funding

This study was supported by the Funding Program for Excellent Talents of Beijing (2017000021469G212).

 Availability of data and materials

The datasets used during the present study are available from the corresponding author upon reasonable request.

 Authors' contributions

LM designed the study. JS and ML performed the experiments. JS and LM analyzed the data. JS and LM draft the manuscript. All authors read and approved the manuscript and agree to be accountable for all aspects of the research in ensuring that the accuracy or integrity of any part of the work are appropriately investigated and resolved.

 Ethics approval and consent to participate

This study does not contain any research using human participants or animals performed by any of the authors.

 Consent for publication

Not applicable.

 Competing interests

The authors state that they have no competing interests.

 References TorreLABrayFSiegelRLFerlayJLortet-TieulentJJemalAGlobal cancer statistics, 2012CA Cancer J Clin6587108201510.3322/caac.2126225651787 ChenWZhengRBaadePDZhangSZengHBrayFJemalAYuXQHeJCancer statistics in china, 2015CA Cancer J Clin66115132201610.3322/caac.2133826808342 ChoiKSJunJKSuhMParkBNohDKSongSHJungKWLeeHYChoiIJParkECEffect of endoscopy screening on stage at gastric cancer diagnosis: Results of the national cancer screening programme in koreaBr J Cancer112608612201510.1038/bjc.2014.608254905284453643 WeiWQYangCXLuSHYangJLiBYLianSYQiaoYLCost-benefit analysis of screening for esophageal and gastric cardiac cancerChin J Cancer30213218201110.5732/cjc.010.10425213526994013318 YehJMHurCWardZSchragDGoldieSJGastric adenocarcinoma screening and prevention in the era of new biomarker and endoscopic technologies: A cost-effectiveness analysisGut65563574201610.1136/gutjnl-2014-308588257795974573370 den HoedCMvan EijckBCCapelleLGvan DekkenHBiermannKSiersemaPDKuipersEJThe prevalence of premalignant gastric lesions in asymptomatic patients: Predicting the future incidence of gastric cancerEur J Cancer4712111218201110.1016/j.ejca.2010.12.01221239166 AlsinaMGulloICarneiroFIntratumoral heterogeneity in gastric cancer: A new challenge to faceAnn Oncol28912913201710.1093/annonc/mdx13428368465 YePZhangMFanSZhangTFuHSuXGavinePRLiuQYinXIntra-tumoral heterogeneity of HER2, FGFR2, cMET and ATM in gastric cancer: Optimizing personalized healthcare through innovative pathological and statistical analysisPLoS One10e0143207201510.1371/journal.pone.0143207265879924654477 ZhangYYangRLianJXuHLncRNA Sox2ot overexpression serves as a poor prognostic biomarker in gastric cancerAm J Transl Res8503550432016279047045126346 HigakiEKuwataTNagatsumaAKNishidaYKinoshitaTAizawaMNittaHNaginoMOchiaiAGene copy number gain of EGFR is a poor prognostic biomarker in gastric cancer: Evaluation of 855 patients with bright-field dual in situ hybridization (DISH) methodGastric Cancer196373201610.1007/s10120-014-0449-925487305 YaoYSYaoRYZhuangLKQiWWLvJZhouFQiuWSYueLMOR1 expression in gastric cancer: A biomarker associated with poor outcomeClin Transl Sci8137142201510.1111/cts.12246254417635350983 KwonCHParkHJLeeJRKimHKJeonTYJoHJKimDHKimGHParkDYSerpin peptidase inhibitor clade A member 1 is a biomarker of poor prognosis in gastric cancerBr J Cancer11119932002201410.1038/bjc.2014.490252116654229634 XuQDongQHeCLiuWSunLLiuJXingCLiXWangBYuanYA new polymorphism biomarker rs629367 associated with increased risk and poor survival of gastric cancer in Chinese by up-regulated miRNA-let-7a expressionPLoS One9e95249201410.1371/journal.pone.0095249247600093997364 QiuMZhouYZhangXWangZWangFShaoJLuJJinYWeiXZhangDLauren classification combined with HER2 status is a better prognostic factor in Chinese gastric cancer patientsBMC Cancer14823201410.1186/1471-2407-14-823253806544289219 CarraSSeguinSJLambertHLandryJHspB8 chaperone activity toward poly(Q)-containing proteins depends on its association with Bag3, a stimulator of macroautophagyJ Biol Chem28314371444200810.1074/jbc.M70630420018006506 HamoudaMABelhaceneNPuissantAColosettiPRobertGJacquelAMariBAubergerPLucianoFThe small heat shock protein B8 (HSPB8) confers resistance to bortezomib by promoting autophagic removal of misfolded proteins in multiple myeloma cellsOncotarget562526266201410.18632/oncotarget.2193250513694171627 KeLMeijeringRAHoogstra-BerendsFMackovicovaKVosMJVan GelderICHenningRHKampingaHHBrundelBJHSPB1, HSPB6, HSPB7 and HSPB8 protect against RhoA GTPase-induced remodeling in tachypaced atrial myocytesPLoS One6e20395201110.1371/journal.pone.0020395217316113123278 CarraSSeguinSJLandryJHspB8 and Bag3: A new chaperone complex targeting misfolded proteins to macroautophagyAutophagy4237239200810.4161/auto.540718094623 HishiyaASalmanMNCarraSKampingaHHTakayamaSBAG3 directly interacts with mutated alphaB-crystallin to suppress its aggregation and toxicityPLoS One6e16828201110.1371/journal.pone.0016828214236623057972 MarshNMWarehamAWhiteBGMiskiewiczEILandryJMacPheeDJHSPB8 and the cochaperone BAG3 are highly expressed during the synthetic phase of rat myometrium programming during pregnancyBiol Reprod92131201510.1095/biolreprod.114.12540125904010 SanbeADaichoTMizutaniREndoTMiyauchiNYamauchiJTanonakaKGlabeCTanoueAProtective effect of geranylgeranylacetone via enhancement of HSPB8 induction in desmin-related cardiomyopathyPLoS One4e5351200910.1371/journal.pone.0005351193991792670514 CrippaVSauDRusminiPBoncoraglioAOnestoEBolzoniEGalbiatiMFontanaEMarinoMCarraSThe small heat shock protein B8 (HspB8) promotes autophagic removal of misfolded proteins involved in amyotrophic lateral sclerosis (ALS)Hum Mol Genet1934403456201010.1093/hmg/ddq25720570967 WilhelmusMMBoelensWCOtte-HöllerIKampsBKustersBMaat-SchiemanMLde WaalRMVerbeekMMSmall heat shock protein HspB8: Its distribution in alzheimer's disease brains and its inhibition of amyloid-beta protein aggregation and cerebrovascular amyloid-beta toxicityActa Neuropathol111139149200610.1007/s00401-005-0030-z16485107 PiccolellaMCrippaVCristofaniRRusminiPGalbiatiMCicardiMEMeroniMFerriNMorelliFFCarraSThe small heat shock protein B8 (HSPB8) modulates proliferation and migration of breast cancer cellsOncotarget81040010415201710.18632/oncotarget.14422280607515354667 SuzukiMMatsushima-NishiwakiRKuroyanagiGSuzukiNTakamatsuRFuruiTYoshimiNKozawaOMorishigeKRegulation by heat shock protein 22 (HSPB8) of transforming growth factor-α-induced ovary cancer cell migrationArch Biochem Biophys5714049201510.1016/j.abb.2015.02.03025731856 MinLZhaoYZhuSQiuXChengRXingJShaoLGuoSZhangSIntegrated analysis identifies molecular signatures and specific prognostic factors for different gastric cancer subtypesTransl Oncol1099107201710.1016/j.tranon.2016.11.003280131685198736 Cancer Genome Atlas Network: Comprehensive molecular portraits of human breast tumoursNature4906170201210.1038/nature11412230008973465532 Cancer Genome Atlas Research NetworkWeinsteinJNCollissonEAMillsGBShawKROzenbergerBAEllrottKShmulevichISanderCStuartJMThe cancer genome atlas pan-cancer analysis projectNature Genet4511131120201310.1038/ng.2764240718493919969 SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAPaulovichAPomeroySLGolubTRLanderESMesirovJPGene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profilesProc Natl Acad Sci USA1021554515550200510.1073/pnas.0506580102161995171239896 KanehisaMGotoSKEGG: Kyoto encyclopedia of genes and genomesNucleic Acids Res282730200010.1093/nar/28.1.2710592173102409 NishimuraDBioCartaBiotech Software & Internet Report2117120200410.1089/152791601750294344 ShimizuMTanakaMAtomiYSmall heat shock protein αB-crystallin controls shape and adhesion of glioma and myoblast cells in the absence of stressPLoS One11e0168136201610.1371/journal.pone.0168136279777385158045 LiuJLuoZZhangLWangLNieQWangZFHuangZHuXGongLArrigoAPThe small heat shock protein αA-crystallin negatively regulates pancreatic tumorigenesisOncotarget765808658242016275884675323194 NakamotoHFujitaKOhtakiAWatanabeSNarumiSMaruyamaTSuenagaEMisonoTSKumarPKGoloubinoffPYoshikawaHPhysical interaction between bacterial heat shock protein (Hsp) 90 and Hsp70 chaperones mediates their cooperative action to refold denatured proteinsJ Biol Chem28961106119201410.1074/jbc.M113.524801244157653937677 BodoorKJalboushSAMatalkaIAbu-SheikhaAWaqRAEbwainiHAbu-AwadAFayyadLAl-ArjatJHaddadYHeat shock protein association with clinico-pathological characteristics of gastric cancer in jordan: HSP70 is predictive of poor prognosisAsian Pac J Cancer Prev1739293937201627644641 KimuraAOgataKAltanBYokoboriTIdeMMochikiEToyomasuYKogureNYanomaTSuzukiMNuclear heat shock protein 110 expression is associated with poor prognosis and chemotherapy resistance in gastric cancerOncotarget71841518423201610.18632/oncotarget.7821269437744951298 GiaginisCDaskalopoulouSSVgenopoulouSSfiniadakisIKouraklisGTheocharisSEHeat shock protein-27, −60 and −90 expression in gastric cancer: Association with clinicopathological variables and patient survivalBMC Gastroenterol914200910.1186/1471-230X-9-14192033812644705 LiXSXuQFuXYLuoWSHeat shock protein 60 overexpression is associated with the progression and prognosis in gastric cancerPLoS One9e107507201410.1371/journal.pone.0107507252076544160299 LiuJXuWSunTWangFPuscheckEBrigstockDWangQTDavisRRappoleeDAHyperosmolar stress induces global mRNA responses in placental trophoblast stem cells that emulate early post-implantation differentiationPlacenta306673200910.1016/j.placenta.2008.10.009190364363343317 KapilaNSharmaAKishoreASodhiMTripathiPKMohantyAKMukeshMImpact of heat stress on cellular and transcriptional adaptation of mammary epithelial cells in riverine buffalo (Bubalus bubalis)PLoS One11e0157237201610.1371/journal.pone.0157237276822565040452

Efficiency of HSPB8 siRNA knockdown and its effect on gastric cancer cell proliferation. (A) HSPB8 expression was higher in five GC cell lines (AGS, BGC-803, BGC-823, SGC-7901 and N87) than the level in a gastric epithelial cell line (GES-1). (B) Efficiency of siRNA knockdown at the HSPB8 protein level in AGS cells. (C) Efficiency of siRNA knockdown at the HSPB8 protein level in BGC-803 cells. (D) Proliferation curves of AGS cells with siRNA knockdown of HSPB8 as determined by CCK-8 assay. (E) Proliferation curves of BGC-803 cells with siRNA knockdown of HSPB8 as determined by CCK-8 assays. (F) OD490 levels of AGS cells transfected with HSPB8-siRNA1 and -siRNA2 and control AGS cells at 48 h. (G) OD490 levels of BGC803 cells transfected with HSPB8-siRNA1 and -siRNA2 and control BGC803 cells at 48 h. Mean ± SD of three independent experiments were displayed with the Student's t-test to calculate the statistical significance. *P<0.05, **P<0.01, ***P<0.001. HSPB8, heat shock protein family B (small) member 8.

Effects of HSPB8 knockdown on gastric cancer cell apoptosis. Apoptosis levels were determined by flow cytometric analysis with Annexin V and 7AAD staining, and all experiments were triplicated. Representative flow cytometric charts of (A) control AGS cells and (B) HSPB8-knockdown AGS cells. (C) Percentages of early and late apoptotic cells in the control and HSPB8-knockdown AGS cells. Representative flow cytometric charts of (D) control BGC-803 cells and (E) HSPB8-knockdown BGC-803 cells. (F) Percentages of early and late apoptotic cells in control and HSPB8-knockdown BGC-803 cells. **P<0.01, ***P<0.001; NS, not significant at P>0.05. HSPB8, heat shock protein family B (small) member 8.

HSPB8 is positively correlated with the MAPK and CREB pathway. TCGA data of 443 gastric cancer patients were used for GSEA analysis. A Pearson correlation based method was adopted and a False Discovery Rate (FDR) <0.25 was considered to be statistically significant. (A) GSEA analysis revealed that the KEGG MAPK signaling pathway is positively correlated with HSPB8. (B) GSEA analysis revealed that the BIOCARTA CREB pathway is positively correlated with HSPB8. Upper panels, enrichment plots of GSEA analysis; lower panels, statistical parameters of the GSEA analysis. HSPB8, heat shock protein family B (small) member 8.

siRNA knockdown of HSPB8 decreases the activity of the ERK-CREB pathway. Protein levels of total and phosphorylated Erk (Erk and pErk) and CREB (CREB and pCREB) in (A) control and HSPB8-knockdown AGS cells and in (B) control and HSPB8-knockdown BGC-803 cells. Protein levels in (C) control and HSPB8-knockdown AGS cells and in (D) control and HSPB8-knockdown BGC-803 cells at different time points. (E) mRNA levels of downstream genes of the ERK-CREB pathway under the siRNA knockdown of HSPB8. Mean ± SD of three independent experiments were displayed with Student's t-test to calculate the statistical significance. *P<0.05, **P<0.01, ***P<0.001. HSPB8, heat shock protein family B (small) member 8.

High expression and low methylation levels of HSPB8 are indicative of a poor prognosis in gastric cancer patients. (A) The Kaplan-Meier curve indicated that gastric cancer patients with high expression levels of HSPB8 had shorter overall survival (OS) than those with low HSPB8 expression. (B) The Kaplan-Meier curve indicated that gastric cancer patients with high expression levels of HSPB8 had shorter disease-free survival (DFS) than those with low HSPB8 expression. (C) The Spearman correlation analysis indicated that mRNA expression and promoter methylation levels of HSPB8 were negatively correlated. (D) The Kaplan-Meier curve indicated that gastric cancer patients with low methylation levels of HSPB8 had shorter OS than those with high methylation levels. HSPB8, heat shock protein family B (small) member 8.

Sequences of the siRNAs used in the present study.

ID Sense (5′-3′) Antisense (5′-3′)
HSPB8 siRNA1 GCATTGTTTCTAAGAACTTCATT TGAAGTTCTTAGAAACAATGCTT
HSPB8 siRNA2 GGTCCCTCCTTACTCAACATT TGTTGAGTAAGGAGGGACCTT
NC siRNA UUCUCCGAACGUGUCACGUTT ACGUGACACGUUCGGAGAATT

Primers for real-time PCR used in the present study.

Gene name Primer sequences
HSPB8 F: CTCCTGCCACTACCCAAGC
R: GGCCAAGAGGCTGTCAAGT
GAPDH F: GGACCTGACCTGCCGTCTAGAA
R: GGTGTCGCTGTTGAAGTCAGAG
MMP9 F: AGACCTGGGCAGATTCCAAAC
R: CGGCAAGTCTTCCGAGTAGT
EGR1 F: GGTCAGTGGCCTAGTGAGC
R: GTGCCGCTGAGTAAATGGGA
CCNA1 F: GAGGTCCCGATGCTTGTCAG
R: GTTAGCAGCCCTAGCACTGTC
CCND1 F: CGGCAAGTCTTCCGAGTAGT
R: CCTCCTTCTGCACACATTTGAA
CCN1 F: GGATGGTAGTTTTGAGTCACCAC
R: CACGAGGATAGCTCTCATACTGT
BCL2 F: CCTCCTTCTGCACACATTTGAA
R: CGGTTCAGGTACTCAGTCATCC

F, forward; R, reverse.