Article Text

Download PDFPDF

Association between COX-2 gene polymorphisms and risk of hepatocellular carcinoma development: a meta-analysis
  1. Si-Cong Lu1,
  2. Jian-Hong Zhong2,
  3. Jun-Tao Tan2,
  4. Hua-Lin Tang2,
  5. Xiao-Guang Liu1,
  6. Bang-De Xiang2,
  7. Le-Qun Li2,
  8. Tao Peng1
  1. 1Hepatobiliary Surgery Department, First Affiliated Hospital of Guangxi Medical University, Nanning, PR China
  2. 2Hepatobiliary Surgery Department, Affiliated Tumor Hospital of Guangxi Medical University, Nanning, PR China
  1. Correspondence to Professor Tao Peng; 3041375003{at}


Objective To investigate the association between cyclo-oxygenase-2 (COX-2) polymorphism and the risk of hepatocellular carcinoma (HCC) development.

Design Systematic review and meta-analysis of COX-2 polymorphism and risk of HCC development among people with or without HCC.

Data sources EMBASE, PubMed, Public Library of Science, SCOPUS, Web of Knowledge and Chinese National Knowledge Infrastructure were searched for all clinical and experimental case–control studies of COX-2 polymorphism and HCC risk. Studies published up to March 2015 were included.

Review method Ten studies were included for data extraction, which were mainly from Asian countries.

Results 2538 people with HCC and 3714 without HCC were found to satisfy the inclusion criteria and included in the review. The associations of specific genotypes in the eight polymorphic variants of COX-2 and the risk of HCC development were analysed. GG genotype at the A-1195G polymorphism may be associated with a reduced risk of HCC development: the OR across all studies was 0.87 (95% CI 0.75 to 1.02) for the G allele versus the A allele, 0.72 (0.53 to 0.97) for GG versus AA, 0.72 (0.57 to 0.92) for GG versus GA+AA and 1.05 (0.77 to 1.44) for AA versus GA+GG. Similar results were found when the meta-analysis was repeated separately for the Chinese subgroup. However, more reliable data are needed to demonstrate associations between variants in G-765C, T+8473C, A-1290G, G-899C and introns 1, 5 and 6 polymorphisms and the risk of HCC development.

Conclusions Only the COX-2 A-1195G gene polymorphism may be associated with a decreased risk of HCC development. These conclusions should be verified in further studies.

  • cyclooxygenase-2
  • hepatocellular carcinoma
  • polymorphism
  • susceptibility

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See:

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • Eight polymorphic variants of cyclo-oxygenase-2 gene were studied.

  • The study was limited by a lack of gene–gene and gene–environment interaction data.

  • This study was also limited by the small sample size in some polymorphic variants.


Hepatocellular carcinoma (HCC) is a significant cause of cancer morbidity and mortality worldwide. The estimated incidence of new HCC cases each year is more than 0.5 million.1 China is one of the regions with highest incidence of HCC (>20 per 100 000 people), which accounts for more than 50% of the total number of cases.2 ,3 Epidemiologically, HCC is strongly associated with hepatitis B or C virus infection, alcohol consumption and metabolic disease. However, not all people with these factors appear to have the same risk of developing HCC. HCC is a multifactorial disease and many studies have shown that gene polymorphisms may also contribute to the risk of hepatocarcinogenesis.4 ,5 Patients with HCC exhibit a high degree of genetic heterogeneity.

Cyclo-oxygenase-2 (COX-2, also known as prostaglandin endoperoxide synthases or prostaglandin H synthases), is an inducible enzyme that converts arachidonic acid to prostaglandins, which are potent mediators of inflammation. COX-2 is not normally present in most tissue cells. It is induced in response to inflammatory cytokines, mitogens, angiogenic growth factors and tumour promoters.6 ,7 Increased COX-2 expression has been associated with the early stages of hepatocarcinogenesis.8 ,9 However, the association of COX-2 genotypes polymorphism with the risk of HCC development has not been established.

Recently, a number of studies10–19 have examined whether an association exists between COX-2 polymorphism and the risk of developing HCC. These studies have reached different conclusions, with some suggesting a significant association and others no association. Since individual case–control studies may fail to detect a complicated genetic relationship because of small sample size, this review aims to comprehensively assess the literature examining a possible link between the COX-2 polymorphism and the risk of HCC development.


Literature search strategy

All clinical and experimental case–control studies of COX-2 polymorphism and the risk of developing HCC published up to 31March 2015 were identified through systematic searches in EMBASE, PubMed, Public Library of Science (, SCOPUS, Web of Knowledge and Chinese National Knowledge Infrastructure. A large number of papers have been published recently by the Public Library of Science and so we also searched this database. No language restriction was imposed. The following search terms were used to identify studies: cyclo-oxygenase-2 or COX-2, gene or polymorphism or variation or genotype or genetic or mutation, “hepatocellular carcinoma” or “liver cancer” or HCC. Detailed database search strategies of EMBASE are shown in table 1. We also searched the Catalog of Published Genome-Wide Association Studies ( of the US National Human Genome Research Institute. Reference lists of these articles and relevant literature from review articles were also examined to identify additional relevant publications.

Table 1

EMBASE search strategies

Inclusion criteria

Only full-length research studies satisfying the following criteria were included in this review: (a) they assessed the association between COX-2 polymorphism and the risk of HCC development; (b) they used a case–control or cohort design in which cases were patients with HCC and controls were healthy people, or people with chronic hepatitis B or C, or with cirrhosis; (c) they focused on human beings; (d) they provided sufficient published data to allow estimation of an OR with a 95% CI. In the case of multiple studies apparently based on the same case or control population, we included only the study with the largest number of participants. Conference abstracts or other summaries were not included. If the data on genotype frequency in the study were incomplete, we tried to contact the authors to collect these data.20

Data extraction

Two authors (S-CL, J-TT) independently searched the literature and identified eligible articles based on our inclusion criteria. These two authors also independently extracted the following data from included studies: first author's family name, year of publication, genotyping methods, source of controls (population-based and hospital-based), numbers and genotypes of cases and controls and Hardy–Weinberg equilibrium (HWE) of controls. Extracted data were compared and discrepancies were resolved by discussion with a third author (J-HZ).

Statistical methods and bias testing

As described previously,20 ,21 the unadjusted OR with 95% CI was used to assess the strength of the association between the COX-2 polymorphism and HCC susceptibility based on the genotype frequencies in cases and controls. The meta-analysis examined the association between different genotypes at different loci of COX-2 and HCC risk by comparing the alleles, comparing homozygous genotypes and applying recessive and dominant genetic models.

A Mantel–Haenszel estimate was used to give a pooled OR using a fixed-effect model, while a DerSimonian–Laird estimate used a random-effect model. The significance of OR was assessed using the Z test, and p<0.05 was considered statistically significant. I2 was used to estimate total variation across studies due to heterogeneity in percentage.22 ,23 A percentage of <25% was considered as a low level of heterogeneity, 25–50% as a moderate level of heterogeneity and >50% as a high level of heterogeneity. I2>50% could suggest heterogeneity and suggest using a random-effect estimate.22 ,23 Otherwise, the fixed-effect model was used to calculate pooled ORs. HWE in the control group was assessed using the χ2 goodness-of-fit test, with p<0.05 considered significant. As far as possible, the meta-analysis was performed according to the PRISMA guidelines.24

As described previously,20 ,21 to detect associations that might be masked in the overall sample, we performed subgroup analyses based on ethnicity. Meta-regression was performed to examine the effect of ethnicity to compare results from the meta-analyses. To assess the reliability of the outcomes in the meta-analysis, a sensitivity analysis was performed by excluding one study at a time.

Publication bias was assessed by visual inspection of Begg's funnel plots. An asymmetric plot suggested possible publication bias, in which case Egger's test was used.25 All statistical tests for this meta-analysis were performed using Stata V.11.0 (Stata-Corp, College Station, USA) and RevMan V.5.3 (Cochrane Collaboration).


Description of studies

Several research databases were searched to identify studies assessing the possible association between the polymorphism in the COX-2 gene and the risk of developing HCC. A total of 562 studies were identified, none of which was a genome-wide association study. This list was reduced to 22 after removing duplicates and screening based on the title and abstract review. These articles were read in full and eight studies were removed because they did not include a control group, while another four studies were removed because overlapping patients were analysed or the data were incomplete. No study which was published in a language other than in Chinese or in English was excluded. Finally, 10 studies were included in the analysis (figure 1).10–19 Four were published in Chinese12 ,14 ,15 ,17 and the other six studies were published in English.10 ,11 ,13 ,16 ,18 ,19 The main characteristics of the included studies are shown in tables 24. All the studies reported that the cases and controls were matched for age and gender.

Table 2

Main characteristics of studies of the cyclo-oxygenase-2 A-1195G polymorphism and the risk of hepatocellular carcinoma

Table 3

Main characteristics of studies of cyclo-oxygenase-2 G-765C polymorphism and the risk of hepatocellular carcinoma

Table 4

Main characteristics of studies of cyclo-oxygenase-2 T+8473C polymorphism and the risk of hepatocellular carcinoma

Figure 1

Flow chart of study selection.

The studies included 2538 people with HCC and 3714 people without HCC. The A-1195G polymorphism in the COX-2 gene and risk of HCC development was reported by eight studies10–17 (table 2), G-765C in six studies10 ,11 ,13 ,14 ,17 ,18 (table 3) and T+8473C in three studies (table 4).10–12

Quantitative data synthesis


Although the polymorphism in the allelic contrast model only slightly affects HCC development risk (OR=0.87, 95% CI 0.75 to 1.02, p=0.09), the GG genotype was significantly associated with a reduced risk across the genetic models tested: the OR across all studies was 0.72 (95% CI 0.53 to 0.97) for the GG versus AA and 0.72 (95% CI 0.57 to 0.92) for GG versus GA+AA (figure 2). However, the AA genotype was not associated with higher or lower HCC development risk: the OR across all studies was 1.05 (95% CI 0.77 to 1.44) for AA versus GA+GG (table 5). The results after deleting each study were similar to those obtained across all studies. We loosely classified the study population as Chinese and non-Chinese based on the ethnicity of the participants. Meta-analyses of subgroups found that Chinese population has the same phenomena as the total population. However, the A-1195G polymorphism in the COX-2 gene was not associated with either increased or reduced risk of HCC development in the non-Chinese population (table 5). Meta-regression also supported our results (table 6).

Table 5

Overall and stratified meta-analyses of the association between COX-2 polymorphisms and risk of hepatocellular carcinoma

Table 6

Ethnicity meta-regression and publication bias of COX-2 A-1195G polymorphisms and risk of hepatocellular carcinoma

Figure 2

Forest plots describing the association of A-1195G COX-2 polymorphism with hepatocellular carcinoma (GG vs GA+AA).


For the COX-2 G-765C polymorphism, no significant association was seen in any of the six studies (C allele vs G allele: OR=1.32, 95% CI 0.76 to 2.30; CC vs GC+GG: OR=0.88, 95% CI 0.16 to 4.75; CC vs GG: OR=0.93, 95% CI 0.16 to 5.35; GG vs CC+GC: OR=0.48, 95% CI 0.14 to 1.59). Since the two non-Chinese studies10 ,13 had a small sample size and the GG genotype was absent in three studies,11 ,14 ,17 subgroup analyses were not performed (table 5).


For the COX-2 T+8473C polymorphism, no significant association was seen in all the three studies (C allele vs T allele: OR=0.99, 95% CI 0.86 to 1.14; CC vs CT+TT: OR=1.31, 95% CI 0.83 to 2.07; CC vs TT: OR=1.25, 95% CI 0.78 to 1.98; TT vs CT+CC: OR=1.05, 95% CI 0.89 to 1.24) (table 5).

Other loci

The study by Chang et al11 also reported three loci polymorphisms in the COX-2 gene: intron 1, intron 5 and intron 6. This study showed that, for each of the six genotypes, no differences in distribution between the HCC and control groups were found. The locus polymorphism of A-1290G was reported by one study with 270 cases and 540 healthy controls.17 This study found no significant association between the COX-2 A-1290G polymorphism and risk of HCC. The locus polymorphism of C-899G in the COX-2 gene was also reported by one study with 300 patients with chronic hepatitis B, 300 patients with liver cirrhosis, 300 patients with HCC and 300 healthy controls.19 This study found that the COX-2-899C genotype may increase the susceptibility of individuals to HCC.

Publication bias and small-study bias

Begg's funnel plots were prepared for the eight studies to assess publication bias for studies of the A-1195G polymorphism of COX-2 and HCC risk. The shape of the funnel plots appeared to be symmetrical for allele contrast, homozygous comparison and recessive and dominant genetic models, suggesting the absence of publication bias. Moreover, Egger's test also suggested no publication bias (table 6).


Some studies reported an association between the COX-2 gene polymorphism and the risk of HCC development, while others found no such association. The most likely reason for the inconsistencies among these studies is the small sample size. To help resolve these conflicting results using a larger sample size, we conducted a systematic review of published studies. In this review, we included 10 studies investigating the association of eight polymorphic variants of COX-2 and the HCC development risk. We found that the GG genotype of A-1195G in the COX-2 gene was associated with a decreased risk of HCC development, especially in a Chinese population. However, we did not find compelling evidence of an association between other COX-2 gene polymorphisms and the risk of HCC development.

As is known, the polymorphisms in the COX-2 promoter may have an important effect on gene transcriptional activity by changing the binding capacity of certain nuclear proteins, thereby affecting COX-2 expression. Even though the exact molecular mechanism still remains unclear, several polymorphisms of COX-2 have been published previously and the results are still conflicting. A previous meta-analysis of eight studies showed that the COX-2 C+202T polymorphism is associated with a lower risk of prostate cancer in Caucasians.26 Another meta-analysis of 25 studies found that the COX-2 A-1195G polymorphism is a low penetrance risk factor for cancer.27 However, COX-2 C-765G and T+8473C polymorphisms are significantly associated with an increased risk of digestive system cancers.28 ,29 The meta-analysis by Bu and Zhao30 included five10–12 ,15 ,17 of the 10 studies included in our review. They found an association between COX-2 A-1195G polymorphism and the risk of HCC development, especially in Asians.

In our updated review, with a larger sample size, another five studies13 ,14 ,16 ,18 ,19 were included. We found that the GG genotype at the A-1195G polymorphism was associated with a reduced risk of HCC development across all studies. We also investigated another seven polymorphic variants (G-765C, T+8473C, intron 1, intron 5, intron 6, A-1290G, C-899G) of COX-2. Although COX-2 C-899G polymorphism may increase the risk of HCC, this result was based on only one study. In order to demonstrate the association between the COX-2 C-899G polymorphism and risk of HCC development, more reliable data with a large sample size are needed.

Development of HCC is a complex, multistep and heterogeneous malignant tumorigenesis process. The aetiology of HCC involves various host and environmental factors. Furthermore, host and environmental factors may interact synergistically in HCC pathogenesis and progression.4 Several studies in this review indicate that COX-2 polymorphisms can interact with environmental factors to modulate the HCC risk. Among people with a history of drinking, COX-2–765 C allele carriers lead to a significantly greater risk of HCC development than the G allele.18 ,31 Though single gene polymorphism and risk of HCC development was not found in the study by Fan et al,12 demographic interactions were seen. Among people aged <55 years, the A allele of the COX-2 A-1195G polymorphism is a high penetrance risk factor for HCC development, whereas among women, the C allele of COX-2 T+8473C is a low penetrance risk factor for HCC. For gene–gene interactions, no significant difference in the frequencies of any combined genotypes was seen between patients with HCC and healthy controls.11 Additionally, no combined effects of COX-2 genotypes and smoking or alcohol drinking were found.11 Moreover, no significant difference in COX-2 C-899G genotype distribution interactions with age, sex or smoking history was found.19 Therefore, whether the gene–gene and gene–environment interactions of the COX-2 polymorphism may contribute to the risk of HCC is unknown.

Our data showed that COX-2 A-1195G gene polymorphism may be a protective factor for hepatocarcinogenesis, but the complete picture is more complex. Seven11 ,12 ,14 ,15 ,17–19 of the 10 studies included Chinese people. China has among the highest incidence of HCC in the world, as well as a high prevalence of hepatitis B virus infection and dietary exposure to aflatoxin B1, which are the two main risk factors for HCC.32–34 Some of the included controls had hepatitis B or C virus infection, or cirrhosis. Because the sample size of these controls is small, a subgroup analysis based on a liver disease background was not performed. In addition, polymorphisms in numerous other genes, such as those encoding microsomal epoxide hydrolase4 and epidermal growth factor5 are associated with the risk of HCC. It may be that any single nucleotide polymorphism such as COX-2 A-1195G or epidermal growth factor 61*A/G is insufficient on its own to cause HCC.

As mentioned, some of the controls had one or more of the following: alcoholic liver disease, hepatitis B or C virus infection and cirrhosis. Since the studies included in this review often did not report detailed statistics on the proportion of patients with HCC or control subjects with these background conditions, we could not perform a subgroup analysis to separate the contribution of COX-2 polymorphism from that of possible confounders like hepatitis B or C infection. In addition, it is hard to assess the quality of the included studies, which may also lead to bias.

Low occurrence of genotypes within the COX-2 G-765C and COX-2 T+8473C polymorphisms may lead to null results in table 5. Therefore, more reliable data with larger sample sizes are needed to determine the relationships involving COX-2 G-765C and COX-2 T+8473C polymorphisms, whose analyses have been underpowered and whose null results must be treated with caution.

Some other limitations of this review should also be considered. Although we searched all the eligible records, the number of studies included was still relatively small. Subgroup stratification analysis of other COX-2 gene polymorphisms was not performed. Moreover, meta-analysis was not carried out for five polymorphic variants of COX-2. Second, the results may be affected by additional confounding factors, such as tumour status, age or gender, but most studies either did not report these baseline data or aggregated them in different ways, making it impossible to include them into the pooled analysis. Moreover, the distribution of genotypes among controls did not show HWE in several studies. Finally, because of the lack of the individual original data, our meta-analysis was based on unadjusted data, and a more precise analysis stratified by clinical manifestation and environmental factors was not performed.

In conclusion, this review suggests that the COX-2 A-1195G gene polymorphism, rather than the other seven polymorphic variants of COX-2, might be a protective factor against HCC development. However, since this review included few studies, large, well-designed studies are warranted to re-evaluate these associations.



  • Contributors S-CL conceived and designed the experiments; S-CL, J-TT, H-LT and J-HZ performed the research; B-DX, L-QL and X-GL performed the statistical analysis; S-CL and J-HZ wrote the manuscript; all authors read and approved the final manuscript.

  • Funding This work was supported by grants in part by the National Nature Science Foundation of China (NSFC 81560535, 81072321, 30760243, 30460143 and 30560133), 2009 Program for New Century Excellent Talents in University, Guangxi Nature Sciences Grant (GuiKeGong 1104003A-7), Guangxi Health Ministry Medicine Grant (Key-Scientific-Research-Grant Z201018), the Innovation Project of Guangxi Graduate Education (YCBZ2015030), grants from Programs for Changjiang Scholars and Innovative Research Team in University (No.IRT1119), and Innovative Research Team in Guangxi Natural Science Foundation(No.2011GXNSFF018005).

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.