Abstract
Yanhua Jiang1,*, Xiandi Wang2,*, Yi Cheng1, Jun Peng1, Jiwei Xiao1, Dingbo Tang1 and Ying Yi1
1Department of Orthopedics, Sichuan Cancer Hospital, Chengdu, China
2Department of Orthopedics, Huashan Hospital, Fudan University, Shanghai, China
*These authors contributed equally to this work
Correspondence to:
Yanhua Jiang, email: [email protected]
Keywords: osteosarcoma, inflammatory gene, polymorphism, meta-analysis
Received: April 27, 2017 Accepted: June 12, 2017 Published: June 29, 2017
ABSTRACT
Associations between inflammatory gene polymorphisms (TNF-α 308G/A, TNF-α 238G/A, TNF-β 252A/G, TGF-β1 29T/C, IL-6 174G/C and IL-10 1082A/G) and osteosarcoma (OS) risk remain unclear. We conducted a systematic search to retrieve studies that investigated associations between inflammatory gene polymorphisms and OS risk. Nine studies that met the inclusion criteria were finally recruited in this meta-analysis. Overall, there was a significant association between TNF-α 308G/A, IL-10 1082A/G and OS risk, while there was no significant association between TNF-α 238G/A, TNF-β 252A/G and IL-6 174G/C and OS risk. Our subgroup analysis showed a significant association between IL-6 174G/C and IL-10 1082A/G and OS risk in Asians, while no such significant correlation was observed with TNF-α 308G/A, TNF-α 238G/A, TNF-β 252A/G and TGF-β1 29T/C polymorphisms. In Caucasians, there was a significant association between TNF-α 238G/A and the decreased incidence of OS. In conclusion, inflammatory gene polymorphisms play a key role in the occurrence and progression of OS. IL-6 174G/C polymorphism was obviously associated with OS risk in Asians, while TNF-α 238G/A polymorphism seemed to be associated with the decreased susceptibility to OS in Caucasians as Altman and Bland test indicated. Although controversial results were observed between IL-10 1082A/G and OS risk in Asians and Caucasians, it is difficult to make a definite conclusion about the role of IL-10 1082A/G polymorphism in the etiology of OS because our Altman and Bland test showed no good evidence to support a different effect in Asians and Caucasians.
INTRODUCTION
Osteosarcoma (OS) is one of the most common primary malignant bone tumors probably due to repeated production of osteoid by malignant cells [1, 2]. OS is a high-grade malignant mesenchymal tumor in the bone tissue that often leads to high disability and mortality in both children and adults, causing great burdens to families and society [1, 3, 4]. According to a survey conducted by the Surveillance, Epidemiology, and End Results (SEER), about 3010 people were diagnosed with primary malignant tumors of bone and joints in the United States in 2013 [5]. In China, OS is also the most common primary malignant bone and soft tissue tumor. According to the database of Beijing Jishuitan Hospital (one of the most famous orthopedics hospitals in China) [6], OS is the most common histologic type of bone tumors, accounting for 22.8% of all bone tumors, followed by giant cell tumor (16.7%), osteochondroma (9.1%), and primary chondrosarcoma (5.6%).
The etiology of OS is complex involving multiple factors, but its actual pathogenesis remains unclear [1, 2, 7, 8]. Abnormal bone growth is reported to be significantly associated with the development and progression of OS [9], and diverse strategies have been attempted to check the progression of OS during the rapid stage of bone growth [10]. Environmental factors have also been found to be associated with OS risk [11, 12]. In addition, age [13], gender [14], ethnicity [15], fluorinated drinking water [16], height [17, 18], and other factors [1] have also been implied as risk factors of OS.
In addition to the above-mentioned risk factors, more studies have verified the importance of genes and gene polymorphisms in the etiology, development and complexity of OS [1, 4, 8]. Copy number gains at chromosomes 1p, 1q, 6p, 8q and 17p, and copy number losses at chromosomes 3q, 6q, 9, 10, 13, 17p and 18q have been reported by many research groups, including the systematic review regarding the role of genes in the etiology of OS conducted by Martin et al. [8]. In addition, more gene polymorphisms have been found to play a key role in the etiology and development of OS, highlighting the importance of gene polymorphisms in this disease [19–22].
Inflammation has been increasingly recognized as an important factor in the etiology and development of cancers, for many studies have demonstrated that inflammatory cells play an important role in the tumor microenvironment, where they foster tumor cell proliferation, survival and migration [23–25] and that gene polymorphisms also play important roles in various cancers and diseases including OS [26–28]. Zhao et al. [29] showed that the AA genotype of TNF-α-308G/A locus of TNF-α gene was a risk factor for OS. However, Oliveira et al. [26] reported that there was no significant association between them. The result of a previous meta-analysis performed to detect the correlation of TNF-α 308G/A, TNF-β (Tumor Necrosis Factor-β) 252A/G and TGF-β1 (Transforming Growth Factor-β1) 29T/C polymorphisms with OS risk [30] showed that TGF-β1 29T/C variants were significantly associated with OS while no association was found between TNF-α -308G/A or TNF-β 252A/G polymorphism and OS risk. But as their meta-analysis only focused on the association of TNF-α 308G/A, TNF-β 252A/G and TGF-β1 29T/C polymorphisms with OS susceptibility, whether other gene polymorphisms such as IL-6 (Interleukin 6), IL-10 (Interleukin 10) and TNF-α 238G/A play a key role in the etiology and development of OS remains unclear. Although increased numbers of inflammatory gene polymorphisms have been found to be associated with susceptibility to OS [26, 31–34], no meta-analysis has been conducted to explore the associations between these novel inflammatory gene polymorphisms and OS risk. The aim of our meta-analysis is to explore the associations of inflammatory gene polymorphisms (TNF-α 308G/A, TNF-α 238G/A, TNF-β 252A/G, TGF-β1 29T/C, IL-6 174G/C and IL-10 1082A/G) with OS risk.
MATERIALS AND METHODS
Literature search
We conducted a systematic online search using ‘PubMed’, ‘EMBASE’, ‘Web of Science’, ‘the Cochrane Library’, ‘Science Direct’, ‘Karger’, ‘Wiley Online Library’ and ‘Springer’ and ‘China WeiPu Library’ to identify case-control studies that investigated the relationship between inflammatory gene polymorphisms and susceptibility to OS. First, the following search terms were used to find out how many inflammatory gene polymorphisms were studied: (‘osteosarcoma’ OR ‘OS’ OR ‘osteosarcoma tumor’) AND (‘polymorphism’ OR ‘single nucleotide polymorphism’ OR ‘SNP’ OR ‘variation’), and then, one-by-one screening was performed by two authors. We found that inflammatory gene polymorphisms TNF-α, TNF-β, TGF-β1, IL-1, IL-6, IL-8, IL-10, IL-12 and IL-16 were studied in case-control studies. Next, we used search terms one-by-one to find out whether there were sufficient articles that reported at least one of the gene polymorphisms identified in our previous search stage. The following search terms were used to find out the eligible studies regarding the association between TNF polymorphisms and OS risk: (‘TNF’ OR ‘Tumor Necrosis Factor) AND (‘osteosarcoma’ OR ‘OS’ OR ‘osteosarcoma tumor’) AND (‘polymorphism’ OR ‘single nucleotide polymorphism’ OR ‘SNP’ OR ‘variation’). The search terms, (‘TGF’ OR ‘Transforming Growth Factor’) AND (‘osteosarcoma’ OR ‘OS’ OR ‘osteosarcoma tumor’) AND (‘polymorphism’ OR ‘single nucleotide polymorphism’ OR ‘SNP’ OR ‘variation’) were used to find the associations between TGF and OS risk. As to IL-6, the following search terms were used: (‘IL 6’ OR ‘Interleukin 6’) AND (‘osteosarcoma’ OR ‘OS’ OR ‘osteosarcoma tumor’) AND (‘polymorphism’ OR ‘single nucleotide polymorphism’ OR ‘SNP’ OR ‘variation’). As to IL-10, the following search terms were used: (‘IL 10’ OR ‘Interleukin 10’) AND (‘osteosarcoma’ OR ‘OS’ OR ‘osteosarcoma tumor’) AND (‘polymorphism’ OR ‘single nucleotide polymorphism’ OR ‘SNP’ OR ‘variation’). And then, these studies were further selected. Furthermore, no language restrictions were applied. Secondary searches of the unpublished literature were conducted by searching the reference lists of the selected studies, reviews and conference reports. Reviews and comments were also examined to further identify eligible studies.
Inclusion and exclusion criteria
The inclusion criteria of our meta-analysis were as follows: (1) case-control studies; (2) evaluation of OS risk and at least one of the identified inflammatory gene polymorphisms were reported; (3) having detailed genotype frequency or numbers of alleles and genotypes between cases and controls or effect size reported or estimated from the articles. The exclusion criteria were: (1) reviews or case reports that were not case-control studies; (2) no available data reported; (3) duplicated reports.
Data extraction
Data from the eligible studies were extracted according to the inclusion and exclusion criteria by two authors, and a consensus was reached. For each study, the following data were collected: author list, year of publication, ethnicity, sample size, alleles, genotype of each gene polymorphism, and HWE (Hardy-Weinberg Equilibrium). Incomplete data were also obtained by contacting the respective corresponding author.
Data synthesis and statistical analysis
Odds ratios (OR) and 95% confidence interval (CI) were calculated to evaluate the association between inflammatory gene polymorphisms and OS. An allele contrast model (A vs. G), heterozygote model (AG vs. GG), homozygote model (AA vs. GG), dominant (AA/AG vs. GG), and recessive (AA vs. AG/GG) model were used for TNF-α 308G/A and TNF-α 238G/A polymorphism. An allele contrast model (G vs. A), heterozygote model (AG vs. AA), homozygote model (GG vs. AA), dominant (GG/AG vs. AA), and recessive (GG vs. AG/AA) model were used for TNF-β 252A/G and IL-10 1082A/G polymorphism. An allele contrast model (C vs. T), heterozygote model (TC vs. TT), homozygote model (CC vs. TT), dominant (CC/TC vs. TT), and recessive (CC vs. TC/TT) model were used for TGF-β1 29T/C polymorphism. And the strength of association between IL-6 174G/C polymorphism and OS susceptibility was evaluated by OR and 95% CI according to allele contrast model (C vs. G), heterozygote model (GC vs. GG), homozygote model (CC vs. GG), dominant (GC/CC vs. GG), and recessive (CC vs. GC/GG) model.
The assumption of heterogeneity was verified by a chi-squared-based Q statistical test and quantified by I2 metric value. ORs were pooled by random effect model across the study, which considers both within-study and between-study variation.
Sensitivity analysis was performed to assess the impact of each study on the combined effect of the present meta-analysis, and subgroup analysis was performed according to the ethnicity of the study populations. Altman and Bland test [35] was conducted to evaluate the difference between Asians and Caucasians. The ratio of relative risk (RRR) with 95% CI and z score were used to assess the interaction. HWE was calculated in control population to evaluate the quality of the data by using chi-square test. All meta-analyses were conducted using Stata 12.0 software (StataCorp, College Station, TX, USA) and a P value below 0.05 was considered statistically significant.
RESULTS
A total of nine studies [26, 29, 32, 33, 36–40] with 1119 cases and 1358 controls were recruited in this meta-analysis. The process of study selection and the inclusion process are shown in Figure 1. Four studies [26, 29, 36, 37] focused on the association between TNF-α 308G/A polymorphism and OS risk, among which two studies [29, 36] also reported both alleles and genotypes of TNF-α 238G/A polymorphism. Oliveira et al. [26] and Xie et al. [37] studied the relationship between TNF-β 252A/G and susceptibility to OS. Ma et al. [38], Xu et al. [39] and Wu et al. [40] reported the alleles and genotypes of TGF-β1 29T/C polymorphism. Oliveira et al. [26] reported the alleles and genotypes of IL-6 174G/C and IL-10 1082A/G polymorphism, and these two polymorphisms were also studied by Qi et al. [32] and Cui et al. [33]. Furthermore, seven studies [29, 32, 33, 37–40] were performed in Asian populations, and the other two studies [26, 36] were conducted in Caucasian populations. In addition, all these studies had complied with HWE. The general characteristics of the recruited studies are summarized in Table 1.
Figure 1: The study selection and inclusion process.
Table 1: General characteristics of studies included in the meta-analysis
Meta-analysis results
In the overall populations, TNF-α 308G/A and IL-10 1082A/G polymorphisms were significantly associated with the increased susceptibility to OS (TNF-α 308G/A: AA vs. GG: OR = 2.99, 95% CI = 1.34–6.68, P = 0.007; AA vs. AG/GG: OR = 2.90, 95% CI = 1.37–6.13, P = 0.005; IL-10 1082A/G (G vs. A: OR = 1.32, 95% CI = 1.06–1.65, P = 0.014; GG vs. AA: OR = 1.72, 95% CI = 1.11–2.67, P = 0.016; GG vs. AA/AG: OR = 1.52, 95% CI = 1.05–2.20, P = 0.028), while no significant association was found between TNF-α 238G/A, TNF-β 252A/G and IL-6 174G/C and OS risk (all P > 0.05) (Tables 2 and 3).
Table 2: Results of genetic models for TNF-α 308G/A, TNF-α 238G/A, TGF-β 252A/G and TGF-β1 29T/C polymorphisms and osteosarcoma
Table 3: Results of genetic models for IL-6 174G/C and IL-10 1082A/G polymorphisms and osteosarcoma
Our subgroup analysis showed that there was no significant association between TNF-α 308G/A and TNF-β 252A/G polymorphisms and susceptibility to OS in either Asian or Caucasian populations (both P > 0.05) (Table 2). Altman and Bland test showed consistency of our subgroup results regarding the role of TNF-α 308G/A and TNF-β 252A/G polymorphisms in OS etiology (Table 2).
It was found in our meta-analysis that TNF-α 238G/A polymorphism was a protective contributor to OS susceptibility in Caucasian populations (A vs. G: OR = 0.18, 95% CI = 0.04–0.80, P = 0.024; AG vs. GG: OR = 0.17, 95% CI = 0.04–0.76, P = 0.020; AA/AG vs. GG: OR = 0.17, 95% CI = 0.04–0.76, P = 0.020, Table 2). However, we could not find any association between TNF-α 238G/A polymorphism and susceptibility to OS in the Asian populations (Table 2), suggesting that TNF-α 238G/A polymorphism might play a different role in the etiology of OS in different ethnicities, as our Altman and Bland test showed that there was good evidence to support a different effect in Asians and Caucasians [A vs. G: ratio of relative risk (95% CI) = 0.19 (0.04–0.96), z score = -2.01; AA/AG vs. GG: ratio of relative risk (95% CI) = 0.20 (0.04–1.00), z score = −2.00]. All these data are shown in Table 2.
As shown in Table 3, although no significant association was observed between IL-6 174G/C polymorphism and OS risk in the overall and Caucasian populations (both P > 0.05), IL-6 174G/C mutation might increase the risk of OS in Asian populations (C vs. G: OR = 1.74, 95% CI = 1.30–2.33, P < 0.001; CC vs. GG: OR = 2.25, 95% CI = 1.27–3.99, P = 0.006; GC vs. GG: OR = 1.85, 95% CI = 1.21–2.83, P = 0.004; CC/GC vs. GG: OR = 1.96, 95% CI = 1.34–2.88, P = 0.001; CC vs. GC/GG: OR = 1.76, 95% CI = 1.02–3.05, P = 0.043). This controversial result might be attributed to the different genetic backgrounds in Asians and Caucasians, as shown by the results of Altman and Bland test [CC vs. GG: ratio of relative risk (95% CI) = 0.29 (0.09–0.91), z score = −2.13; CC vs. GC/GG: ratio of relative risk (95% CI) = 0.45 (0.21–0.97), z score = −2.04], while other factors such as the limited number of included studies or the small number of included cases in individual studies could not be ignored.
A significant association was also observed between IL-10 1082A/G polymorphism and OS risk in Asian populations (G vs. A: OR = 1.39, 95% CI = 1.06–1.81, P = 0.017; GG vs. AA: OR = 1.85, 95% CI = 1.10–3.09, P = 0.020; GG vs. AA/AG: OR = 1.80, 95% CI = 1.10–2.97, P = 0.020, Table 3); however, we did not find any significant association in Caucasian populations (all P > 0.05) (Table 3). We could hardly make a certain conclusion regarding the role of IL-10 1082A/G polymorphism in the etiology of OS, for our Altman and Bland test showed that there was no good evidence to support a different effect in Asians and Caucasians [G vs. A: ratio of relative risk (95% CI) = 0.86 (0.53–1.38), z score = −0.64; CC vs. GG vs. AA: ratio of relative risk (95% CI) = 0.77 (0.29–2.06), z score = −0.52; AG vs. AA: ratio of relative risk (95% CI) = 1.12 (0.45–2.82), z score = 0.24; AG/GG vs. AA: ratio of relative risk (95% CI) = 1.01 (0.43–2.38), z score = 0.02; GG vs. AA/AG: ratio of relative risk (95% CI) = 0.68 (0.32–1.43), z score = −1.01].
As to TGF-β1 29T/C polymorphism, no significant association with OS was observed in Asian populations (P > 0.05) (Table 2). As no study reported the association between TGF-β1 29T/C polymorphism and OS risk, whether it plays a key role in Caucasian populations remains unclear.
Sensitivity analysis and publication bias
A leave-one-out analysis was performed to estimate the sensitivity of our meta-analysis. Any single study could be omitted, without any effect on the overall statistical significance, indicating that the results are stable. Furthermore, we did not assess publication bias because the number of studies analyzed in each polymorphism is smaller than 10.
DISCUSSION
Inflammation is the body’s immune response to infection, injury, insult or any other acute and chronic factors [41]. Immune cells and inflammatory cytokines such as tumor necrosis factor (TNF), tumor growth factor and interleukin (IL) are commonly involved in the process of inflammation [41]. The continuous presence of these inflammatory cytokines may alter the tissue microenvironment and cause mutation of normal cells. In addition, increasing attention has been paid to the role of immune cells and inflammatory cytokines in cancer etiology and progression [42–44]. Several immune-based therapies have been proposed to treat different cancers in clinical practice [45, 46]. More researchers have recognized the importance of inflammatory cytokines and immune cells in epithelial cancers [41], especially in mesenchymal stem cell-derived tumors such as OS [27, 47].
Multiple factors are involved in the etiology of OS, among which genetics and gene polymorphisms have aroused increased attention [1, 4, 8]. Given the important roles of inflammation and gene polymorphisms in cancer etiology and progression, more studies have been performed to explore the associations between inflammatory gene polymorphisms and susceptibility to OS [26, 29, 32, 33, 36–40, 48, 49]. Although several studies including Patio-Garcia et al. [36], Oliveira et al. [26], Zhao et al. [29] and Xie et al. [37] have explored the possible association between TNF-α 308G/A polymorphism and OS risk, there results are inconsistent. Other TNF polymorphisms have also been genotyped in OS patients vs. healthy subjects in an attempt to define the association of these polymorphisms with OS risk [26, 29, 36, 37]. Although a previous meta-analysis [30] had investigated the role of TNF-α 308G/A and TNF-β 252A/G polymorphisms in OS risk, it had the following demerits: 1) it only included TNF-α 308G/A, TNF-β 252A/G and TGF-β1 29T/C polymorphisms, devoid of other inflammatory gene polymorphisms such as IL-6 and IL-10; 2) it did not include the results reported by Zhao et al. [29], thus weakening its statistical power; 3) it did not discuss the impact of different ethnicities on the obtained result, knowing that ethnicity is an important factor affecting the association between gene polymorphisms and the disease. In addition, the results about IL-6 174G/C polymorphism are inconsistent between the study of Oliveira et al. [26] and that of Qi et al. [32], and the results about IL-10 1082A/G polymorphism are inconsistent between the study of Oliveira [26] and that of Cui et al. [33]. To the best of our knowledge, this is the first meta-analysis reporting the role of IL-6 174G/C and IL-10 1082A/G polymorphism in the etiology and progression of OS. This is also the first study reporting the association of TNF-β 252A/G, IL-6 174G/C and IL-10 1082A/G polymorphisms with OS, and the impact of ethnicity on OS susceptibility. Compared with previous meta-analyses, we retrieved more articles [29, 32, 33, 40] to explore the association of inflammatory gene polymorphisms with OS susceptibility.
TNF-α 308G/A polymorphism and OS risk
TNF-α 308G/A polymorphism is a change G-to-A single nucleotide polymorphism (SNP) at position 308, affecting TNF-α gene regulation and associated with altered transcriptional activity in many diseases [50]. It has been found to play an important role in the etiology of many diseases such as systemic lupus erythematosus (SLE) [50], coronary artery disease (CAD) [51], and inflammatory bowel disease (IBD) [52]. The association between TNF-α 308G/A polymorphism and OS susceptibility has also been explored by many researchers. The results of studies performed by Oliveira et al. [26] and Patio-Garcia et al. [36] in Caucasian populations showed no significant association between them, while Zhao et al. [29] reported controversial results in Asian populations. Given these conflicting results, we performed this meta-analysis to determine whether TNF-α 308G/A polymorphism was truly associated with OS susceptibility in both Caucasian and Asian populations. Four relevant studies [26, 29, 36, 37] were included in this meta-analysis and the results showed that TNF-α 308G/A polymorphism was associated with an increased OS risk in the overall populations. Our study provided relatively objective results regarding the role of TNF-α 308G/A polymorphism and OS risk. Although these results are still inconsistent with the result reported by Bian et al. [30], they seem more creditable because they were retrieved from more studies. We believe that the most important factor contributing to this difference is the number of studies involved in our research. Our study was based on four eligible studies including Zhao’s study [29], which was not included in Bian’s study [30]. Another advantage of our study over Bian’s [30] is the sub-group analysis on ethnicity. As no sub-group analysis was performed in Bian’s study [30], their result only represents the overall population and cannot be used in either Caucasian or Asian population. In addition, no association was observed in either Asian or Caucasion population in sub-group analysis, although contradictory results were found in the overall populations. The combination of different original data in each study might have great impact on the pooled distribution of each genotype, which might be an important contributor to the different results of overall populations and Asians and Caucasians.
TNF-α 238G/A polymorphism and OS risk
Another common SNP of TNF-α is TNF-α 238G/A polymorphism. It is recently believed that the A-238 allele at position 238 can significantly reduce the TNF-α production; on the other side, it may be associated with down-regulation of tissue inflammation [53]. Interestingly, we found that TNF-α 238G/A polymorphism played a protective role in the etiology of OS in Caucasian populations. Multiple factors are believed to be involved in the etiology of cancers or tumors, including gene, gene polymorphism, immune system, tissue microenvironment and eating habit. Inflammatory cytokines such as TNF-α occupy much of the tumor microenvironment, fostering tumor cell proliferation, survival and migration [23–25]. Therefore, inflammatory cytokines are considered important stimulators or factors contributing to many types of cancer and sarcoma. According to this consensus, it is easy to understand why TNF-α 238G/A polymorphism could significantly decrease the susceptibility to OS in Caucasian populations: A-238 allele of TNF-α 238G/A could significantly reduce the TNF-α production, and then it would indirectly reduce the immune response and change the immune microenvironment in tissues, thus reducing OS risk in Caucasian populations. However, the exact mechanism about how TNF-α 238G/A polymorphism affects TNF-α function and expression and plays a protective role in OS risk is not fully understood and needs further exploration [36]. In addition, there is no significant association between TNF-α 238G/A polymorphism and OS risk in either the overall population and Asian population, which further demonstrates the different roles of gene polymorphisms in various ethnicities. We could think that TNF-α 238G/A polymorphism might play different roles in the etiology of OS in different ethnicities as our Altman and Bland test showed that there was good evidence to support a different effect in Asians and Caucasians. However, what we should not ignore is that the limited number of included studies or the small number of included cases in individual studies might also contribute to this unremarkable relationship in Asian populations, which might not have enough statistical power to find out this difference.
TNF-β 252A/G polymorphism and OS risk
TNF-β is a sub-type of the TNF family and believed to be involved in a variety of inflammatory, immune-stimulatory and antiviral responses [54]. Meanwhile, it is generally accepted that a SNP at position 252 in the intron of TNF-β called TNF-β 252A/G polymorphism leads to the overexpression of TNF-β, and it is considered a candidate biomarker of gastric cancer susceptibility [54]. Due to the function of TNF-β 252A/G polymorphism in many diseases [54–56], the relationship between TNF-β 252A/G polymorphism and OS susceptibility has also been determined [26, 37]. Oliveira ID et al. [26] studied TNF-β 252A/G polymorphism in 80 osteosarcoma patients and 160 control individuals using polymerase chain reaction-restriction fragment length polymorphism method. Although no remarkable relationship was observed in their study, they found that the Caucasian patients with a variant genotype (GG) of the 252A/G TNF-β gene showed an event-free survival rate of 20% at 100 months. In their opinion, the high-dose local TNF management selectively destroyed tumor blood vessels and had powerful anticancer action, but when chronically produced these cytokines might act as an endogenous tumor promoter, contributing to the tissue remodeling and stromal development necessary for tumor growth and spread [26]. This may be the reason why TNF-β 252A/G polymorphism had impact on the event-free survival rate. However, no significant association of TNF-β 252A/G polymorphism with the survival rate was found in Asian populations [37]. The inconsistent backgrounds, the relatively small sample size and different genotyping technologies might be important contributors to this difference. Our study showed that TNF-β 252A/G polymorphism was not significantly associated with OS susceptibility in either overall populations, Asian populations or Caucasian populations, which is consistent with Oliveira’s study [26] and Xie’s study [37]. However, whether TNF-β 252A/G polymorphism affects the survival rate in Asian populations remains unknown and needs further investigation.
TGF-β1 29T/C polymorphism and OS risk
Three studies [38–40] were performed to detect the relationship between TGF-β1 29T/C polymorphism and OS susceptibility, hoping to find the important roles of TGF-β1 29T/C polymorphism in the etiology of OS, knowing that this cytokine plays a key role in cancer etiology [57]. TGF-β1 is known to inhibit cellular proliferation or promote cellular differentiation and apoptosis in normal cells [57]. However, like TNF, the expression of TGF-β1 appears to be increased in cancer cells [57]. Wu et al. [40] reported that patients carrying TT genotype had a higher risk to get OS than patients carrying CC genotype, and that the morbidity and metastasis of OS were relevant to TGF-β1 29T/C polymorphism. These results are consistent with Xu et al. [39], who reported that the prevalence of TT genotype and T allele was significantly elevated in patients as compared with healthy controls. According to Ma et al. [38], the TT genotype among 29TC polymorphism of TGF-β1 gene might also be a genetic factor that aroused OS. However, no significant association was observed between TGF-β1 29T/C polymorphism and OS susceptibility. Our meta-analysis recruited the original data of Wu et al. [40], Xu et al. [39] and Ma et al. [38], based on which we performed a pooled analysis. When all the data were merged, the distribution of TT, CT and CC genotypes might be changed among OS patients and healthy subjects, which might be the most important contributor to the different results of our meta-analysis and the above three studies. Although the meta-analysis of Bian et al. [30] found the same result as that of the above three studies, our study might be more creditable than theirs, as we included more studies with larger sample sizes. In addition, all the studies included in our meta-analysis were performed in Asian populations, and whether TGF-β1 29T/C polymorphism increases or decreases OS susceptibility in Caucasian populations remains unclear and more studies are required.
IL-6 174 G/C and IL-10 1082A/G polymorphisms and OS risk
Interleukins are important members of inflammatory cytokines that have been regarded as key contributors to tumor development with respect to proliferation, apoptosis, angiogenesis and differentiation [32, 33, 58, 59]. IL-6 and IL-10 are the most commonly studied interleukins and always studied together in many studies, knowing that both IL-6 and IL-10 are pro-inflammatory cytokines that may enhance the pro-inflammatory status in different physiologic and pathophysiologic processes. It has been well documented that IL-6 and IL-10 polymorphisms play a key role in many diseases [58, 60–62]. Many studies have also been conducted to explore the association between IL-6 and IL-10 polymorphisms and OS [26, 32, 33]. Oliveira [26] performed the first study to determinate whether there existed difference in the distribution of IL-6 and IL-10 polymorphisms between OS patients and health controls but found no remarkable correlation in Caucasian populations. In contrast, Qi et al. [32] reported that the genotype CC of IL-6 174 G/C carried a higher risk of OS metastasis and later Enneking stages, compared with the GG genotype. They thought that the IL-6 −174 G/C genotype was associated with the risk for OS development and metastasis in Chinese Han population. The ethnicity, sample size, selection bias of cases and control subjects, and genotyping technologies might contribute to this difference. Our pooled results showed that there was no significant association between IL-6 174 G/C polymorphism and OS susceptibility in Caucasian populations, and the controversial results were observed in Asian population, which is consistent with Oliveira et al.’s [26] and Qi et al.’s study [32]. These controversial results might be attributed to the different genetic backgrounds in Asians and Caucasians because the results of Altman and Bland test in our study showed that there was good evidence to support a different effect in Asians and Caucasians, while we could not ignore the influences of other factors such as the limited number of included studies or the small number of included cases in individual studies. In addition, we also found an insignificant association in the overall populations, which might result from the changed gene distributions in the pooled analysis. As to IL-10 1082A/G polymorphism, two studies were included in our meta-analysis [26, 33]. In the overall populations, our meta-analysis showed a significant association between 1082A/G polymorphism and OS risk, indicating that the IL-10 1082G allele might be the genetic risk and susceptibility factor for developing OS. It has been verified that the 1082A allele could down-regulate the IL-10 transcription, whereas the IL-10 1082G allele has an increasing effect on IL-10 expression [63]. Moreover, IL-10 has impact on pro-inflammatory cytokine gene transcription and stability by reducing protein translation, which might be an important contributor to the etiology and development of many diseases [63]. Therefore, it is easy to understand why the IL-10 1082G allele in OS patients significantly differed from that in the control group. Besides, Cui et al. [33] found that there was also a statistically significant association between the IL-10 1082A/G genotype and poor survival, further verifying the role of 1082G allele in the etiology and development of OS and survival of OS patients. Our sub-group analysis revealed a remarkable relationship between IL-10 1082A/G polymorphism and OS risk in Asian populations, which is consistent with Cui’s study [33], while the IL-10 1082G allele might not have significant impact on the etiology of OS in Caucasian populations. However, we could not only consider the different role of IL-10 1082A/G in OS etiology as the contributor to the contrast results in Asian and Caucasian populations, as our Altman and Bland test showed that there was no good evidence to support a different effect in Asians and Caucasians. Other factors, such as the limited number of included studies, the small number of included cases in individual studies, and different genotyping techniques might also have great impact. Therefore, the effect of 1082A/G polymorphism on OS in different ethnicities needs to be further explored.
Despite the comprehensive analysis on the association between inflammatory gene polymorphisms and OS, our meta-analysis still has some limitations. First, of the nine studies that eventually satisfied the eligibility criteria of the present study, only two studies were included in the pooled analysis of TNF-a 238G/A, TNF-β 252A/G, IL-6 174G/C, and IL-10 1082A/G polymorphisms. Second, the sample size of some studies used in this meta-analysis is still relatively small, which might contribute to the inconsistent results and affect the conclusions drawn. Third, no study was conducted to explore the association between TGF-β1 29T/C polymorphism in Caucasian populations, which may not be sufficient to make a convincing conclusion on all ethnicities and indicate a race-specific effect. Therefore, larger-scale and better-designed studies are necessary to determine the association between inflammatory gene polymorphisms (TNF-α 308G/A, TNF-α 238G/A, TNF-β 252A/G, TGF-β1 29T/C, IL-6 174G/C and IL-10 1082A/G) and OS susceptibility.
CONCLUSIONS
In summary, this is the first meta-analysis to determinate the association of TNF-α 308G/A, TNF-α 238G/A, TNF-β 252A/G, TGF-β1 29T/C, IL-6 174G/C and IL-10 1082A/G polymorphisms with OS susceptibility. Inflammatory gene polymorphisms play a key role in the occurrence and progression of OS. IL-6 174G/C polymorphism was obviously associated with OS risk in Asian populations, while TNF-α 238G/A polymorphism might be associated with the decreased susceptibility to OS in Caucasians. Although controversial associations were observed between IL-10 1082A/G and OS risk in Asian and Caucasian populations, it is difficult to make a definite conclusion about the role of IL-10 1082A/G polymorphism in the etiology of OS because our Altman and Bland test showed that there was no good evidence to support a different effect in Asians and Caucasians.
CONFLICTS OF INTEREST
The authors have no conflicts of interest to declare.
REFERENCES
1. Lindsey BA, Markel JE, Kleinerman ES. Osteosarcoma Overview. Rheumatol Ther. 2017; 4:25–43.
2. Geller DS, Gorlick R. Osteosarcoma: a review of diagnosis, management, and treatment strategies. Clin Adv Hematol Oncol. 2010; 8:705–18.
3. PDQ Pediatric Treatment Editorial Board. Osteosarcoma and Malignant Fibrous Histiocytoma of Bone Treatment (PDQ®): Health Professional Version. 2017. PDQ Cancer Information Summaries. Bethesda (MD): National Cancer Institute (US).
4. Durfee RA, Mohammed M, Luu HH. Review of Osteosarcoma and Current Management. Rheumatol Ther. 2016; 3:221–43.
5. Siegel R, Naishadham D, Jemal A. Cancer statistics, 2013. CA Cancer J Clin. 2013; 63:11–30.
6. Niu X, Xu H, Inwards CY, Li Y, Ding Y, Letson GD, Bui MM. Primary Bone Tumors: Epidemiologic Comparison of 9200 Patients Treated at Beijing Ji Shui Tan Hospital, Beijing, China, With 10165 Patients at Mayo Clinic, Rochester, Minnesota. Arch Pathol Lab Med. 2015; 139:1149–55.
7. Zhou W, Hao M, Du X, Chen K, Wang G, Yang J. Advances in targeted therapy for osteosarcoma. Discov Med. 2014; 17:301–07.
8. Martin JW, Squire JA, Zielenska M. The genetics of osteosarcoma. Sarcoma. 2012; 2012:627254.
9. Byrum SD, Washam CL, Montgomery CO, Tackett AJ, Suva LJ. Proteomic technologies for the study of osteosarcoma. Sarcoma. 2012; 2012:169416.
10. Battaglia S, Dumoucel S, Chesneau J, Heymann MF, Picarda G, Gouin F, Corradini N, Heymann D, Redini F. Impact of oncopediatric dosing regimen of zoledronic acid on bone growth: preclinical studies and case report of an osteosarcoma pediatric patient. J Bone Miner Res. 2011; 26:2439–51.
11. Watanabe A, Yoneyama S, Nakajima M, Sato N, Takao-Kawabata R, Isogai Y, Sakurai-Tanikawa A, Higuchi K, Shimoi A, Yamatoya H, Yoshida K, Kohira T. Osteosarcoma in Sprague-Dawley rats after long-term treatment with teriparatide (human parathyroid hormone (1-34)). J Toxicol Sci. 2012; 37:617–29.
12. Paioli A, Rocca M, Cevolani L, Rimondi E, Vanel D, Palmerini E, Cesari M, Longhi A, Abate ME, Marchesi E, Picci P, Ferrari S. Erratum to: Osteosarcoma follow-up: chest X-ray or computed tomography? Clin Sarcoma Res. 2017; 7:7.
13. Pakos EE, Nearchou AD, Grimer RJ, Koumoullis HD, Abudu A, Bramer JA, Jeys LM, Franchi A, Scoccianti G, Campanacci D, Capanna R, Aparicio J, Tabone MD, et al. Prognostic factors and outcomes for osteosarcoma: an international collaboration. Eur J Cancer. 2009; 45:2367–75.
14. Smeland S, Müller C, Alvegard TA, Wiklund T, Wiebe T, Björk O, Stenwig AE, Willén H, Holmström T, Follerås G, Brosjö O, Kivioja A, Jonsson K, et al. Scandinavian Sarcoma Group Osteosarcoma Study SSG VIII: prognostic factors for outcome and the role of replacement salvage chemotherapy for poor histological responders. Eur J Cancer. 2003; 39:488–94.
15. Pui CH, Boyett JM, Hancock ML, Pratt CB, Meyer WH, Crist WM. Outcome of treatment for childhood cancer in black as compared with white children. The St Jude Children’s Research Hospital experience, 1962 through 1992. JAMA. 1995; 273:633–37.
16. Gandhi D, Naoghare PK, Bafana A, Kannan K, Sivanesan S. Fluoride-Induced Oxidative and Inflammatory Stress in Osteosarcoma Cells: Does It Affect Bone Development Pathway? Biol Trace Elem Res. 2017; 175:103–11.
17. Longhi A, Pasini A, Cicognani A, Baronio F, Pellacani A, Baldini N, Bacci G. Height as a risk factor for osteosarcoma. J Pediatr Hematol Oncol. 2005; 27:314–18.
18. Mirabello L, Pfeiffer R, Murphy G, Daw NC, Patiño-Garcia A, Troisi RJ, Hoover RN, Douglass C, Schüz J, Craft AW, Savage SA. Height at diagnosis and birth-weight as risk factors for osteosarcoma. Cancer Causes Control. 2011; 22:899–908.
19. Qu WR, Wu J, Li R. Contribution of the GSTP1 gene polymorphism to the development of osteosarcoma in a Chinese population. Genet Mol Res. 2016; 15:15.
20. Zhang S, Hou C, Li G, Zhong Y, Zhang J, Guo X, Li B, Bi Z, Shao M. A single nucleotide polymorphism in the 3′-untranslated region of the KRAS gene disrupts the interaction with let-7a and enhances the metastatic potential of osteosarcoma cells. Int J Mol Med. 2016; 38:919–26.
21. Adiguzel M, Horozoglu C, Kilicoglu O, Ozger H, Acar L, Ergen A. MMP-3 gene polymorphisms and Osteosarcoma. Indian J Exp Biol. 2016; 54:175–79.
22. Zhao J, Chen ZQ, Li GW, Yang M, Shao J, Li M. The associations of vascular endothelial growth factor gene polymorphisms with susceptibility to osteosarcoma: evidence from a meta-analysis. Eur J Cancer Care (Engl). 2017; 26: e12513.
23. He Q, Tian L, Jiang H, Zhang J, Li Q, Sun Y, Zhao J, Li H, Liu M. Identification of laryngeal cancer prognostic biomarkers using an inflammatory gene-related, competitive endogenous RNA network. Oncotarget. 2017; 8:9525–34. https://doi.org/10.18632/oncotarget.13627.
24. Zhou ZH, Ji CD, Zhu J, Xiao HL, Zhao HB, Cui YH, Bian XW. The prognostic value and pathobiological significance of Glasgow microenvironment score in gastric cancer. J Cancer Res Clin Oncol. 2017; 143:883–94.
25. Kwak T, Drews-Elger K, Ergonul A, Miller PC, Braley A, Hwang GH, Zhao D, Besser A, Yamamoto Y, Yamamoto H, El-Ashry D, Slingerland JM, Lippman ME, Hudson BI. Targeting of RAGE-ligand signaling impairs breast cancer cell invasion and metastasis. Oncogene. 2017; 36:1559–72.
26. Oliveira ID, Petrilli AS, Tavela MH, Zago MA, de Toledo SR. TNF-alpha, TNF-beta, IL-6, IL-10, PECAM-1 and the MPO inflammatory gene polymorphisms in osteosarcoma. J Pediatr Hematol Oncol. 2007; 29:293–97.
27. Hong JT, Son DJ, Lee CK, Yoon DY, Lee DH, Park MH. Interleukin 32, inflammation and cancer. Pharmacol Ther. 2017; 174:127–37.
28. Liu B, Yan S, Qu L, Zhu J. Celecoxib enhances anticancer effect of cisplatin and induces anoikis in osteosarcoma via PI3K/Akt pathway. Cancer Cell Int. 2017; 17:1.
29. Zhao Z, Tang X, Song K, Li X, Zhang Y. Association of -308G/A and -238G/A polymorphisms of TNF-α and osteosarcoma risk. Int J Clin Exp Pathol. 2015; 8:4177–81.
30. Bian Z, He Q, Wang X, Li M, Zhu L. Association of genetic polymorphisms with osteosarcoma risk: a meta-analysis. Int J Clin Exp Med. 2015; 8:8317–28.
31. Chen Y, Yang Y, Liu S, Zhu S, Jiang H, Ding J. Association between interleukin 8 -251 A/T and +781 C/T polymorphisms and osteosarcoma risk in Chinese population: a case-control study. Tumour Biol. 2016; 37:6191–96.
32. Qi Y, Zhao C, Li H, Zhang B, Tada K, Abe H, Tada M. Genetic variations in interleukin-6 polymorphism and the association with susceptibility and overall survival of osteosarcoma. Tumour Biol. 2016; 37:9807–11.
33. Cui Y, Zhu JJ, Ma CB, Cui K, Wang F, Ni SH, Zhang ZY. Interleukin 10 gene -1082A/G polymorphism is associated with osteosarcoma risk and poor outcomes in the Chinese population. Tumour Biol. 2016; 37:4517–22.
34. Tang YJ, Wang JL, Xie KG, Lan CG. Association of interleukin 16 gene polymorphisms and plasma IL16 level with osteosarcoma risk. Sci Rep. 2016; 6:34607.
35. Altman DG, Bland JM. Interaction revisited: the difference between two estimates. BMJ. 2003; 326:219.
36. Patio-Garcia A, Sotillo-Pieiro E, Modesto C, Sierrases-Maga L. Analysis of the human tumour necrosis factor-alpha (TNFalpha) gene promoter polymorphisms in children with bone cancer. J Med Genet. 2000; 37:789–92.
37. Xie J, Liu S, Lu B. [Relationship between polymorphism of tumor necrosis factor and osteosarcoma]. [Article in Chinese]. Zhonghua Shi Yan Wai Ke Za Zhi. 2008; 25:723–24.
38. Ma J, Zhou D, Ma X. Association of the 29T > C polymorphism of TGF-β1gene with osteosarcoma. Orthopedic Journal of China. 2010; 18:1510–13.
39. Xu S, Yang S, Sun G, Huang W, Zhang Y. Transforming growth factor-beta polymorphisms and serum level in the development of osteosarcoma. DNA Cell Biol. 2014; 33:802–06.
40. Wu Y, Zhao J, He M. Correlation between TGF-β1 gene 29 T > C single nucleotide polymorphism and clinicopathological characteristics of osteosarcoma. Tumour Biol. 2015; 36:5149–56.
41. Shivappa N, Hébert JR, Rosato V, Garavello W, Serraino D, La Vecchia C. Inflammatory potential of diet and risk of oral and pharyngeal cancer in a large case-control study from Italy. Int J Cancer. 2017; 141:471–79.
42. Espinoza-Sánchez NA, Chimal-Ramírez GK, Mantilla A, Fuentes-Pananá EM. IL-1β, IL-8, and Matrix Metalloproteinases-1, -2, and -10 Are Enriched upon Monocyte-Breast Cancer Cell Cocultivation in a Matrigel-Based Three-Dimensional System. Front Immunol. 2017; 8:205.
43. Pastorino U, Morelli D, Marchianò A, Sestini S, Suatoni P, Taverna F, Boeri M, Sozzi G, Cantarutti A, Corrao G. Inflammatory status and lung function predict mortality in lung cancer screening participants. Eur J Cancer Prev. 2017.
44. Culig Z, Puhr M. Interleukin-6 and prostate cancer: current developments and unsolved questions. Mol Cell Endocrinol. 2017; S0303–7207:30181–8.
45. Abdel-Wahab N, Alshawa A, Suarez-Almazor ME. Adverse Events in Cancer Immunotherapy. Adv Exp Med Biol. 2017; 995:155–74.
46. Hendriks D, Choi G, de Bruyn M, Wiersma VR, Bremer E. Antibody-Based Cancer Therapy: Successful Agents and Novel Approaches. Int Rev Cell Mol Biol. 2017; 331:289–383.
47. Kelleher FC, O’Sullivan H. Monocytes, Macrophages, and Osteoclasts in Osteosarcoma. J Adolesc Young Adult Oncol. 2017; 6:396–405.
48. Bilbao-Aldaiturriaga N, Patino-Garcia A, Martin-Guerrero I, Garcia-Orad A. Cytotoxic T lymphocyte-associated antigen 4 rs231775 polymorphism and osteosarcoma. Neoplasma. 2017; 64:299–304.
49. Tang YJ, Wang JL, Nong LG, Lan CG, Zha ZG, Liao PH. Associations of IL-27 polymorphisms and serum IL-27p28 levels with osteosarcoma risk. Medicine (Baltimore). 2014; 93:e56.
50. Yang ZC, Xu F, Tang M, Xiong X. Association Between TNF-α Promoter -308 A/G Polymorphism and Systemic Lupus Erythematosus Susceptibility: A Case-Control Study and Meta-Analysis. Scand J Immunol. 2017; 85:197–210.
51. Kazemi E, Jamialahmadi K, Avan A, Mirhafez SR, Mohiti J, Pirhoushiaran M, Hosseini N, Mohammadi A, Ferns GA, Pasdar A, Ghayour-Mobarhan M. Association of tumor necrosis factor-α -308 G/A gene polymorphism with coronary artery diseases: an evidence-based study. J Clin Lab Anal. 2017; 96:e6569.
52. Al-Meghaiseeb ES, Al-Robayan AA, Al-Otaibi MM, Arfin M, Al-Asmari AK. Association of tumor necrosis factor-α and -β gene polymorphisms in inflammatory bowel disease. J Inflamm Res. 2016; 9:133–40.
53. Sghaier I, Zidi S, Mouelhi L, Dabbech R, Ghazouani E, Brochot E, Stayoussef M, Yacoubi-Loueslati B. The relationship between TNF alpha gene polymorphisms (-238/-308), TNF RII VNTR (p75) and outcomes of hepatitis B virus infection in Tunisian population. Gene. 2015; 568:140–45.
54. Xu Z, Shi R, Zhang R, Zhang D, Wang L. Association between tumor necrosis factor β 252 A/G polymorphism and risk of gastric cancer: a meta-analysis. Tumour Biol. 2013; 34:4001–05.
55. Liu R, Song Q, Zhang Z, Yang Z, Liu Y, Zhang X. [TNF-β genetic variant contributes to the risk of non-small cell lung cancer]. [Article in Chinese]. Zhonghua Yu Fang Yi Xue Za Zhi. 2015; 49:41–44.
56. Jeng JE, Wu HF, Tsai MF, Tsai HR, Chuang LY, Lin ZY, Hsieh MY, Chen SC, Chuang WL, Wang LY, Yu ML, Dai CY, Tsai JF. Independent and additive interaction between tumor necrosis factor β +252 polymorphisms and chronic hepatitis B and C virus infection on risk and prognosis of hepatocellular carcinoma: a case-control study. Asian Pac J Cancer Prev. 2014; 15:10209–15.
57. Migita K, Miyazoe S, Maeda Y, Daikoku M, Abiru S, Ueki T, Yano K, Nagaoka S, Matsumoto T, Nakao K, Hamasaki K, Yatsuhashi H, Ishibashi H, Eguchi K. Cytokine gene polymorphisms in Japanese patients with hepatitis B virus infection—association between TGF-beta1 polymorphisms and hepatocellular carcinoma. J Hepatol. 2005; 42:505–10.
58. Giannitrapani L, Soresi M, Balasus D, Licata A, Montalto G. Genetic association of interleukin-6 polymorphism (-174 G/C) with chronic liver diseases and hepatocellular carcinoma. World J Gastroenterol. 2013; 19:2449–55.
59. Mocellin S, Marincola F, Rossi CR, Nitti D, Lise M. The multifaceted relationship between IL-10 and adaptive immunity: putting together the pieces of a puzzle. Cytokine Growth Factor Rev. 2004; 15:61–76.
60. Zhao J, Yang M, Li M. Association of IL-6 and MMP-3 gene polymorphisms with susceptibility to adolescent idiopathic scoliosis: a meta-analysis. J Genet. 2016; 95:573–79.
61. Sghaier I, Mouelhi L, Rabia NA, Alsaleh BR, Ghazoueni E, Almawi WY, Loueslati BY. Genetic variants in IL-6 and IL-10 genes and susceptibility to hepatocellular carcinoma in HCV infected patients. Cytokine. 2017; 89:62–67.
62. Crosbie EJ, Einstein MH, Franceschi S, Kitchener HC. Human papillomavirus and cervical cancer. Lancet. 2013; 382:889–99.
63. Rad R, Dossumbekova A, Neu B, Lang R, Bauer S, Saur D, Gerhard M, Prinz C. Cytokine gene polymorphisms influence mucosal cytokine expression, gastric inflammation, and host specific colonisation during Helicobacter pylori infection. Gut. 2004; 53:1082–89.
