Abstract

The polymorphisms in the three main heat shock protein 70 (HSP70-1, HSP70-2, and HSP70-hom) genes were identified to be associated with cancer risk. However, the results are inconsistent. We perform a meta-analysis to evaluate the association between the three HSP70 polymorphisms and cancer risk. Relevant studies were identified using PubMed, Web of Science, Chinese National Knowledge Infrastructure (CNKI), and Wanfang databases up to March 29, 2014. The cancer risk associated with the HSP70 polymorphisms was estimated for each study by odds ratios (OR) together with its 95% confidence interval (CI), respectively. Twenty case-control studies from eighteen publications were included; a significant association was observed for HSP70-2 polymorphism (dominant model: OR = 1.53, 95% CI: 1.11–2.09; recessive model: OR = 1.91, 95% CI: 1.06–3.45; AG versus AA: OR = 1.38, 95% CI: 1.03–1.84; GG versus AA: OR = 2.34, 95% CI: 1.21–4.54), while there was no significant association for HSP70-1 and HSP70-hom polymorphisms. Besides, in stratification analyses by ethnicity, cancer type, and source of control, significant association was detected for HSP70-2 polymorphism, while for HSP70-hom polymorphism, we found a significant association in hospital-based population under homozygote comparison model. This meta-analysis suggests that the HSP70-2 polymorphism rather than HSP70-hom and HSP70-1 polymorphisms was associated with the risk of cancer.

1. Introduction

Cancer is recognized as one of the leading causes of death in economically developed countries as well as in developing countries. According to the estimation of GLOBOCAN, approximately 12.7 million new cases and 7.6 million deaths of cancer had occurred in 2008, it has become a major public health challenge [1]. Although the mechanism of carcinogenesis is still not fully understood, it has been suggested that environmental factors, interplaying with low-penetrance susceptibility genes, may be important in the development of cancer [2, 3].

Heat shock proteins (HSPs) are evolutionarily highly conserved stress proteins expressed and induced by heat shock, infection, inflammation, ischemia hypoxia, oxidative stress, carcinogens, and so on [4, 5]. Mammalian HSPs have been classified into six families according to their molecular weight. Among them, the HSP70 family is one of the most conservative and well-known HSPs. There are three main genes (HSPA1A, HSPA1B, and HSPA1L) in the human HSP70 family and the coding proteins were defined as HSP70-1, HSP70-2, and HSP70-hom, respectively. All three genes were located in class III region of the human major histocompatibility complex (MHC) on chromosome 6. HSPA1A and HSPA1B genes encode an identical protein, the heat inducible HSP70 protein, whereas the HSPA1L gene encodes a non-heat inducible protein that shares 90% sequence identity with HSP70 protein [6]. It has been proposed that the HSP70 plays an important role in tumor development, treatment, and prognosis and has distinct immunologic mechanisms affecting tumor cells and promoting cell growth [7, 8]. In cancer cells, the expression of HSP70 is abnormally high and the protein may participate in oncogenesis and in resistance to chemotherapy [9]. Several single nucleotide polymorphisms (SNPs) have been described in these genes. The most studied regions are located at positions +1267 of HSP70-2 (rs1061581), +2437 of HSP70-hom (rs2227956), and +190 of HSP70-1 (rs1043618). These SNPs could affect HSP70 expression or function and further contribute to disease susceptibility and stress tolerance [10].

To date, several studies have investigated the association between the three HSP70 polymorphisms and risk of cancer [1131]. However, the results remain controversial. Therefore, we conduct a meta-analysis to evaluate the association between the HSP70 polymorphisms and cancer susceptibility.

2. Materials and Methods

2.1. Search Strategy

We searched for relevant studies before March 29, 2014, by using electronic PubMed, Web of Science, Chinese National Knowledge Infrastructure (CNKI), and Wanfang databases with the following terms: “heat shock protein 70 or HSP70,” “genetic polymorphism or polymorphisms or variant or SNP,” and “cancer or carcinoma or tumor.” The search was restricted to humans and without language restrictions. Additional studies were identified by a hand search of references of original or review articles on this topic. If more than one geographic or ethnic heterogeneous group was reported in one report, each was extracted separately.

2.2. Inclusion Criteria and Exclusion Criteria

The inclusion criteria were as follows: (1) studies that evaluated the association between the HSP70 polymorphisms and cancer, (2) a case-control study design, and (3) studies that had detailed genotype frequency of cases and controls or could be calculated from the paper text, while major exclusion criteria were (1) case-only study, case reports, and review articles, (2) studies without the raw data of the HSP70 genotype, and (3) studies that compared the HSP70 variants in precancerous lesions.

2.3. Data Extraction

The following information was extracted from each eligible publication: the first author’s name, year of publication, country of origin, ethnicity, cancer type, source of control, genotyping methods, number of cases and controls, and Hardy-Weinberg equilibrium (HWE) in controls ( value). All data were extracted by two investigators independently, using the same standard. The results were compared and disagreements were resolved by consensus.

2.4. Statistical Analysis

The risk of cancer associated with the HSP70 polymorphisms was estimated for each study by odds ratio (OR) and 95% confidence interval (95% CI). Four different ORs were calculated: dominant model (the combined variant homozygote and heterozygote versus the wild-type homozygote), recessive model (the variant homozygote versus the combined heterozygote and wild-type homozygote), heterozygote comparison (heterozygote versus the wild-type homozygote), and homozygote comparison (variant homozygote versus the wild-type homozygote). A -test-based statistic test was performed to assess the between-study heterogeneity [32]. We also quantified the effect of heterogeneity by test. When a significant test () or indicated heterogeneity across studies, the random effects model was used [33] or else the fixed effects model was chosen [34]. In addition, we tested whether genotype frequencies of controls were in HWE using test. We performed stratification analyses on ethnicity, cancer type, and source of controls. Analysis of sensitivity was performed to evaluate the stability of the results. Finally, potential publication bias was investigated using Begg’ funnel plot and Egger’s regression test [35, 36]. was considered statistically significant.

All analyses were performed by the Cochrane Collaboration RevMan 5.2 and STATA package version 12.0 (Stata Corporation, College Station, Texas).

3. Results

3.1. Study Characteristics

After an initial search, a total of 132 published articles relevant to the topic were identified. According to the inclusion criteria, 21 studies [1131] with full-text were included in this meta-analysis and 111 studies were excluded. The flow chart of study selection is summarized in Figure 1. As shown in Table 1, because the study by Chouchane et al. [11] included two cancer types and the study by Ucisik-Akkaya et al. [26] included two populations, we treated them separately in this meta-analysis. In addition, we excluded three studies [2931] because they included the overlapped data with those included in the analysis [27]. Therefore, there were 17 case-control studies with 2134 cases and 2818 controls concerning HSP70-2 polymorphism, 10 studies with 2042 cases and 2661 controls concerning HSP70-hom polymorphism, and 5 studies with 1558 cases and 2075 controls concerning HSP70-1 polymorphism. Of the 19 eligible studies, four ethnicities were addressed: 8 studies on Asians, 4 studies on Europeans, 5 studies on Africans, and 3 studies on mixed populations. 3 studies focused on hepatocellular carcinoma, 3 studies on gastric cancer, 3 studies on breast cancer, and 11 studies on others (2 on lung, childhood acute lymphoblastic leukemia (ALL) and one on colorectal, pancreatic, prostate, and Kangri cancer, multiple myeloma, and non-Hodgkin’s lymphoma, resp.). The distribution of genotypes in the controls was consistent with the Hardy-Weinberg equilibrium for all selected studies except for six studies [11, 19, 20, 23, 28] for HSP70-2 polymorphism and four studies [18, 20, 26] for HSP70-hom polymorphism.

3.2. Quantitative Data Synthesis

For HSP70-2 polymorphism, 17 studies with 2134 cases and 2818 controls were identified. Overall, a significant association was found (dominant model: OR = 1.53, 95% CI: 1.11–2.09; recessive model: OR = 1.91, 95% CI: 1.06–3.45; AG versus AA: OR = 1.38, 95% CI: 1.03–1.84; GG versus AA: OR = 2.34, 95% CI: 1.21–4.54) (Figure 2). In stratified analysis by ethnicity, we found that the polymorphism was associated with an increased risk of cancer in Asians (dominant model: OR = 1.96, 95% CI: 1.10–3.51; AG versus AA: OR = 1.67, 95% CI: 1.03–2.71), and Africans (recessive model: OR = 7.06, 95% CI: 2.33–21.41; GG versus AA: OR = 7.56, 95% CI: 2.44–23.39) but not for other populations. In the stratified analysis based on cancer type, a significant association was detected in hepatocellular carcinoma (dominant model: OR = 2.41, 95% CI: 1.50–3.87; recessive model: OR = 4.98, 95% CI: 3.18–7.79; AG versus AA: OR = 1.80, 95% CI: 1.34–2.42; GG versus AA: OR = 6.07, 95% CI: 2.79–13.19), and we failed to detect any association between them among breast and other cancers. Stratification based on the source of controls showed significant associations between the polymorphism and risk of cancer in the population-based subgroup (dominant model: OR = 1.57, 95% CI: 1.06–2.34; recessive model: OR = 2.69, 95% CI: 1.21–5.97; GG versus AA: OR = 3.27, 95% CI: 1.36–7.83); however, no significant association was found in the hospital-based subgroup (Table 2).

For HSP70-hom polymorphism, 10 studies with 2042 cases and 2661 controls were identified. Overall, no significant association was found under all genetic models (dominant model: OR = 0.85, 95% CI: 0.53–1.37; recessive model: OR = 1.05, 95% CI: 0.50–2.18; TC versus TT: OR = 0.84, 95% CI: 0.55–1.29; CC versus TT: OR = 0.98, 95% CI: 0.43–2.24). When the analysis was stratified by ethnicity, similar results were observed among Asian, European, African, and mixed populations. While, in stratified analysis by source of controls, a significant association was found in the hospital-based subgroup (CC versus TT: OR = 3.66, 95% CI: 1.03–13.02) but not in the population-based subgroup (Table 2).

For HSP70-1 polymorphism, 5 studies with 1558 cases and 2075 controls were identified. The pooled results suggested that no significant association was found in overall analysis (dominant model: OR = 1.30, 95% CI: 0.75–2.24; recessive model: OR = 1.02, 95% CI: 0.52–2.00; GC versus GG: OR = 1.26, 95% CI: 0.75–2.12; CC versus GG: OR = 1.13, 95% CI: 0.52–2.45) (Table 2).

3.3. Heterogeneity and Sensitivity Analysis

For all three polymorphisms, substantial heterogeneities were observed among overall studies in all four genetic models (HSP70-2: dominant model: , , recessive model: , , AG versus AA: , , GG versus AA: , ; HSP70-hom: dominant model: , ; recessive model: , ; TC versus TT: , ; CC versus TT: , ; HSP70-1: dominant model: , ; recessive model: , ; GC versus GG: , ; CC versus GG: , ) (Table 2). Therefore, we conducted stratified analysis by ethnicity, cancer type, and source of the controls to find the potential sources of heterogeneity. For HSP70-2 polymorphism, we found that heterogeneity significantly reduced or removed among Africans, mixed populations, hepatocellular, and breast cancers. However, heterogeneity still exists in Asians, Europeans, other cancers, population-based, and hospital-based populations. For HSP70-hom polymorphism, heterogeneity significantly reduced or removed in mixed populations and hospital-based populations, but it still exists among Asians, Africans, and population-based populations. Then, sensitivity analysis, after removing one study at a time, was performed to evaluate the stability of the results. When excluded the study by Partida-Rodríguez et al. or Wang et al. in HSP70-1 polymorphism, the study by Medhi et al. in HSP70-hom polymorphism, the heterogeneity was effectively decreased, which suggests that the particular study may be the source of heterogeneity. In addition, no other single study influenced the pooled OR qualitatively as indicated, suggesting that the results of this meta-analysis are credible.

3.4. Publication Bias

Begg’s funnel plot and Egger’s test were performed to assess the potential publication bias in the available literature. The shape of funnel plots did not reveal any evidence of funnel plot asymmetry (data not shown). Egger’s test also showed that there was no statistical significance for the evaluation of publication bias under dominant model (HSP70-2 polymorphism: , HSP70-hom polymorphism: , and HSP70-1 polymorphism: ).

4. Discussion

To our knowledge, this is the first meta-analysis which comprehensively assessed the associations between HSP70 polymorphisms and cancer risk. In this study, we found significant associations in the overall comparison for HSP70-2 polymorphism. Individuals with the AG/GG genotype could have an increased risk of cancer. However, we failed to detect any association for HSP70-1 and HSP70-hom polymorphisms. Moreover, in the stratified analyses by several variables, including ethnicity, cancer type, and source of the controls, significant association was detected among Asians, Africans, hepatocellular carcinoma, and population-based population for HSP70-2 polymorphism, while for HSP70-hom polymorphism, we observed a significant association in hospital-based population under homozygote comparison model.

The HSP70 family is the most important and best characterized family of stress proteins. It acts as a chaperone molecule for antigenic peptides derived from tumor cells, leading to an antitumor immune recognition by cytotoxic T lymphocytes [3739]. In addition, HSP70 is induced in tumor cells to overcome the stressful conditions faced by the tumor, such as lack of nutrients, oxygen, or antitumor immune response contributing to their survival [40]. Recently, the three main polymorphisms in these genes have investigated the association with many cancers, such as gastric, colorectal, hepatocellular, and breast cancer. As for the HSP70-2 and HSP70-hom polymorphisms, the HSP70-2 A+1267G polymorphism is a synonymous mutation located at the coding region and is likely to affect the secondary structure of mRNA, thus affecting the stability of mRNA and protein expression, and the HSP70-hom T > C polymorphism has nonsynonymous mutations, which leads to a Met to Thr substitution at position 493 in the peptide binding domain and may affect substrate binding specificity and chaperone activity of HSP70 [6]. Some investigations demonstrated that genetic alteration of the HSP70-2 and/or HSP70-hom can modulate cancer susceptibility and that the frequency of the variant genotype was significantly higher in patients when compared with controls [11, 12, 18, 22]. However, the association of Ala variants and cancer risk was not validated by others [16, 21, 25, 27, 28]. Besides, Ucisik-Akkaya et al. [26] reported that the HSP70-2 polymorphism may be a protective factor for acute lymphoblastic leukemia in both Welsh and Mexican populations. With regard to HSP70-1, the HSP70-1 G+190C polymorphism base pair lies upstream of the translation initial site [6]. In a study from Mexico, Partida-Rodríguez et al. [21] reported that HSP70-1 C/G showed significant association with gastric cancer; similarly, Wang et al. [27] suggest that HSP70-1 G+190C may contribute to individual susceptibility to lung cancer in a Chinese Han population; however, Guo et al. [13] found that the HSP70-1 polymorphism was not associated with lung cancer risk.

In this meta-analysis, we found that individuals with AG/GG genotype had a higher risk of developing cancer under all four models in HSP70-2 polymorphism; besides, in the stratified analyses by ethnicity, cancer type, and source of control, we found that G allele carriers had a higher risk of cancer than AA genotype carriers in Asians, Africans, hepatocellular carcinoma, and population-based population. With regard to HSP70-hom and HSP70-1 polymorphisms, the genotype distribution between cancer and control was not of significant difference. The inconsistent results may be attributed to differences in genetic backgrounds, environmental factors, and other factors, such as small sample size or inadequate adjustment for confounding factors. For example, the distribution of the AA genotype is about twice as frequent in the Chinese Han population, and the frequency of the GG genotype is similar, slightly above 25%, in the Costa Rican, Mexican, and Chinese population, whereas it does not reach 10% in the Tunisian, Indian, and Japanese groups [12]. Furthermore, the interaction among some other SNPs might affect the relationship of each polymorphism included with the development of cancer. Some reports suggest possible linkage disequilibrium between HSP70 and TNF SNPs, since TNF and HSP70 gene families are located just 600 kb apart from each other [11, 19], which indicate that TNF and HSP70 may act as endogenous tumor promoters in vivo. HSP70 effects, thereby, may be modulated by the TNF genotype. Additionally, because only few studies on subgroup (such as hepatocellular carcinoma, breast cancer, Africans, and hospital-based) were included, the results should be interpreted with caution, and more studies are needed.

Two significant issues should be addressed in this study, that is, heterogeneity and publication bias, which may influence the results of meta-analysis. We do not detect a significant publication bias in this meta-analysis, suggesting the reliability of our results. With regard to heterogeneity, substantial heterogeneities were observed among overall studies in all four genetic models for all three polymorphisms, when stratified analysis by ethnicity, cancer type and source of the controls were conducted. For HSP70-2 polymorphism, we found that heterogeneity significantly reduced or removed among Africans, mixed populations, and hepatocellular and breast cancers but not among Asians, Europeans, other cancers, population-based, and hospital-based populations. For HSP70-hom polymorphism, heterogeneity significantly reduced or removed in mixed populations and hospital-based populations, but it still exists among Asians, Africans, and population-based populations. When excluded the study by Partida-Rodríguez et al. or Wang et al. in HSP70-1 polymorphism, the study by Medhi et al. in HSP70-hom polymorphism, the heterogeneity was effectively decreased. The results above suggest that different ethnicity, tumor types, control selection, and particular study may be the source of heterogeneity.

This meta-analysis has limitations that must be acknowledged. First, because of incomplete raw data or publication limitations, some relevant studies could not be included in our analysis. Second, the controls included in our analysis were selected variously either from populations or hospitals. Therefore, misclassification bias was possible because these studies may have included control groups who have different risks of developing cancer. Third, our results were based on unadjusted estimates, while lacking of the information (such as age, gender, family history and other risk factors) for the date analysis may cause serious confounding bias.

In summary, this meta-analysis suggested that the HSP70-2 polymorphism rather than HSP70-hom and HSP70-1 polymorphisms was associated with the risk of cancer. However, large and well-designed studies taking into consideration gene-gene and gene-environment interactions are warranted to validate our findings.

Conflict of Interests

The authors have declared that no conflict of interests exists.