Abstract
Several potential functional polymorphisms (Arg194Trp, Arg280His, Arg399Gln) in the DNA base excision repair gene X-ray repair cross-complementing group 1 (XRCC1) have been implicated in cancer risk. Our meta-analysis on total of 11,957 cancer cases and 14,174 control subjects from 38 published case-control studies showed that the odds ratio (OR) for the variant genotypes (Trp/Trp + Arg/Trp) of the Arg194Trp polymorphism, compared with the wild-type homozygote (Arg/Arg), was 0.89 [95% confidence interval (95% CI), 0.81-0.98] for all tumor types without between-study heterogeneity. Similarly, the overall risk for the combined variant genotypes (His/His + Arg/His) of the Arg280His, compared with the wild homozygote (Arg/Arg), was 1.19 (95% CI, 1.00-1.42). However, there was no main effect in either recessive or dominant modeling for the Arg399Gln, and the variant Gln/Gln homozygote was not associated with overall cancer risk (OR, 1.01; 95% CI, 0.90-1.14). The analyses suggest that XRCC1 Arg194Trp, Arg280His polymorphisms may be biomarkers of cancer susceptibility and a single larger study with thousands of subjects and tissue-specific biochemical and biological characterization is warranted to further evaluate potential gene-to-gene and gene-to-environment interactions on XRCC1 polymorphisms and cancer risk.
- X-ray repair cross-complementing group 1 (XRCC1)
- polymorphism
- cancer
- meta-analysis
- susceptibility
Introduction
Human cancer can be initiated by DNA damage caused by UV, ionizing radiation, and environmental chemical agents. To safeguard the integrity of genome, humans have developed a set of complex DNA repair systems. Among the five main DNA maintenance mechanisms operating in mammals, base excision repair is the primary guardian against damage that results from cellular metabolism, including reactive oxygen species, methylation, deamination, and hydroxylation. Therefore, base excision repair is a universal event in the cells and is relevant for preventing mutagenesis.
X-ray repair cross complementing group 1 (XRCC1), one of the >20 genes that participate in base excision repair pathway, encodes a scaffolding protein that functions in the repair of single-strand breaks, the most common lesions in cellular DNA (1). Both biological and biochemical evidence indicates a direct role for XRCC1 in base excision repair because it interacts with a complex of DNA repair proteins, including poly(ADP-ribose) polymerase, DNA ligase 3, and DNA polymerase-β (1-3). There are a total of eight nonsynonymous coding single nucleotide polymorphisms in XRCC1, three of which are common (variant allele frequency > 0.05) and lead to amino acid substitutions in XRCC1 at codon 194 (exon 6, base C to T, amino acid Arg to Trp), codon 280 (exon 9, base G to A, amino acid Arg to His), and codon 399 (exon 10, base G to A, amino acid Arg to Gln) (http://egp.gs.washington.edu). Because the Arg399Gln polymorphism is located in the region of the BRCT-I interaction domain of XRCC1 within a poly(ADP-ribose) polymerase binding region, this polymorphism has been extensively investigated both in its function and in its association with cancer risk. The presence of the variant Gln399 allele has been shown to be associated with measurable reduced DNA repair capacity as assessed by the persistence of DNA adducts (4-6), elevated levels of sister chromatid exchanges (5, 7), increased RBC glycophorin A (4), p53 mutations (8), and prolonged cell cycle delay (9). However, Taylor et al. (10) reported that whereas BRCT-I is critical for XRCC1-dependent SSBR for maintenance of genetic integrity, the Arg399Gln polymorphism in BRCT-I does not have a significant impact on this function and negative findings were also obtained from other individual studies (11-13). A large number of molecular epidemiologic studies have been conducted to evaluate the role of the Arg399Gln polymorphism on cancer risk; however, the results remain conflicting rather than conclusive (6, 14-49).
Both the XRCC1 Arg194Trp and Arg280His variants occur in the newly identified proliferating cell nuclear antigen binding region (50), which consists of polar Pro-, Ser-, and Arg/Lys-rich regions. The transition from a positively charged Arg to a hydrophobic Trp within the conserved region may alter XRCC1 function. Recently, Wang et al. (51) reported that individuals with the variant Trp194 allele had fewer bleomycinor benzo(a)pyrene diol epoxide–induced chromosomal breaks than those with wild-type genotype; however, others did not find a significant association of Arg194Trp with altered levels of DNA adducts (4) and G2 cell cycle delay (9). Molecular epidemiologic studies on the association between this polymorphism and cancer risk also presented contradicting results (14-17, 19-23, 26-29, 31-34, 39-43, 49). To date, there are relatively few studies conducted to examine the association between Arg280His variant and cancer risk (19, 20, 22, 31, 33, 35, 40, 52) and only one study evaluated the association of Arg280His and altered DNA adducts (4). Because a single study may have been underpowered to detect the effect of low-penetrance genes and particularly their dose-response relationships, a quantitative synthesis to accumulate data from different studies may provide evidence on the association of genetic polymorphisms with cancer risk. In this meta-analysis, we aimed to obtain summary risk estimates for the above-mentioned three nonsynonymous coding single nucleotide polymorphisms of XRCC1 associated with cancer risk, as well as to quantify the potential between-study heterogeneity.
Materials and Methods
Identification and Eligibility of Relevant Studies
We have attempted to include all the case-control association studies of cancer with genotyping data for at least one of the three polymorphisms, Arg399Gln, Arg194Trp, and Arg280His. Eligible studies were identified by searching the electronic literature MEDLINE for relevant reports (last search update July 15, 2004, using the search terms “XRCC1 and cancer”). Additional studies were identified by a hand search of references of original studies or review articles on this topic and by personal contact with the authors if necessary.
A total of 43 published studies examined the relationship between XRCC1 polymorphisms and cancer risk, three of which were excluded because they did not have an appropriate case-control design (53) or they did not present detailed genotyping information for the three polymorphisms (54, 55). Another three studies were excluded because they investigated the same or a subset population of reported articles (56-58). Hence, the data for this analysis were available from 38 case-control studies, including 11,542 cancer cases and 13,694 controls for Arg399Gln (from 37 studies), 4,933 cancer cases and 6,775 controls for Arg194Trp (from 22 studies), and 1,688 cancer cases and 2,129 controls for Arg280His (from 8 studies).
Data Extraction
Two investigators independently extracted data and reached a consensus on all of the items. The following information was sought from each article: author, journal and year of publication, country of origin, selection and characteristics of cancer cases and controls, demographics, ethnicity, and genotyping information. For studies including subjects of different ethnicities, data were extracted separately and categorized as European, African, and Asian. However, if the authors did not clearly state the ethnic information or we could not separate them according to genotypes, the term “mixed ethnicity” was used (Table 1).
Number of cases and controls, allele frequencies, and quality score of the studies included in the meta-analysis
Meta-analysis
The risks (odds ratios, ORs) of cancer associated with the XRCC1 polymorphisms were estimated for each study. For the Arg399Gln, we first estimated the risk of the variant genotype Gln/Gln, compared with the wild-type Arg/Arg homozygotes, and then evaluated the risks of Gln/Gln versus (Arg/Gln + Arg/Arg) and (Gln/Gln + Arg/Gln) versus Arg/Arg, which assumed recessive and dominant effects, respectively, of the variant Gln399 allele. For the Arg194Trp and Arg280His polymorphisms, we evaluated only the risk of combined variant genotypes (Trp/Trp + Arg/Trp for Arg194Trp and His/His + Arg/His for Arg280His) versus their wild-type homozygote Arg/Arg, because of the rare variant allele frequencies of these two polymorphisms.
In addition to comparisons for total subjects, studies were categorized into different subgroup analyses according to the ethnicity and tumor type (if the tumor type contains less than three individual studies, it was categorized into the “other cancer” group). For each subgroup, we estimated the between-study heterogeneity across the eligible comparisons using the χ2-based Q test (59) and the heterogeneity was considered significant for P < 0.05. Values from single studies were combined using models of both fixed effects (Mantel-Haenszel) and random effects (DerSimonian and Laird; ref. 60). Random effects incorporate an estimate of the between-study variance and tend to provide wider confidence intervals when the results of the constituent studies differ among themselves. In the absence of between-study heterogeneity, the two methods provide identical results. We also did cumulative meta-analysis to evaluate whether the summary OR for the allele contrasts changed over time as more data accumulated (61). Inverted funnel plots and the Egger's test were used to provide diagnosis of publication bias (linear regression analysis; ref. 62).
All analyses were done in Statistical Analysis System software (v.8.0; SAS Institute, Cary, NC) and Review Manage (v.4.2; Oxford, England). All the P values were two-sided.
Results
Meta-analyses Databases
We established a database according to the extracted information from each article. Table 1 lists the cancer type of the study, ethnicity of the population, and the number of cases and controls for each XRCC1 polymorphisms. There are total 37 case-control studies concerning Arg399Gln polymorphism, 22 for Arg194Trp, and 8 for Arg280His. All studies indicated that the distribution of genotypes in the controls was consistent with Hardy-Weinberg equilibrium [goodness of fit χ2 test, degree of freedom (df) = 1], except for three studies for Arg399Gln (21, 27, 44), one study for Arg194Trp (32), and one study for Arg280His (20). However, only 27% of the studies for Arg399Gln, 14% for Arg194Trp, and none for Arg280His had a statistical power >80% based on the assumption that the cancer risk (OR) associated with the variant genotypes was >1.5. According to quality control on genotyping, three studies obtained DNA from surgically resected “normal” tissues adjacent to the tumor instead of peripheral blood (18, 28, 52). A classic PCR-RFLP assay was done in 81% of the studies, 41% randomly repeated a portion of samples while genotyping, and 22% used other genotyping assay to validate the data. Only 19% of the studies described use of blindness of the case-control status of DNA samples while genotyping.
Quantitative Synthesis
XRCC1 Arg399Gln. The eligible studies included 11,542 cancer patients and 13,694 control subjects. There were significant differences in terms of the variant Gln399 allele frequency between the three major ethnicities [European, 34.7%; 95% confidence interval (95% CI), 33.8-35.6; Asian, 26.5%; 95% CI, 25.6-27.4; African, 15.5%; 95% CI, 13.5-17.7; P < 0.0001; Fig. 1]. Figure 2 shows the cancer risks (ORs) associated with the XRCC1 Gln/Gln genotype compared with the Arg/Arg genotype. Overall, individuals carrying the XRCC1 Gln/Gln genotype did not have elevated cancer risk compared with the individuals with the Arg/Arg genotype (OR, 1.01; 95% CI, 0.90-1.14; P = 0.02 for heterogeneity), and this negative associations were also observed in subgroups stratified by cancer type and ethnicity (Fig. 2). Similarly, no association with cancer risk was found, neither in the recessive (Gln/Gln versus Arg/Gln + Arg/Arg: OR, 1.03; 95% CI, 0.92-1.15, P = 0.02 for heterogeneity) nor in the dominant model of the Gln399 allele (Gln/Gln + Arg/Gln versus Arg/Arg: OR, 1.00; 95% CI, 0.95-1.05; Table 2). In the stratified analyses, however, significantly increased risks were found in the Asian subjects (OR, 1.18; 95% CI, 1.00-1.39) in the recessive model and among the African subjects (OR, 1.45; 95% CI, 1.13-1.88) in the dominant model (Table 2).
Allele frequencies and their 95% CIs of the XRCC1 codon 399 and 194 polymorphisms among control subjects by different ethnicity. Each data point represents a separate study for the indicated association.
ORs (log scale) of cancer associated with XRCC1 codon 399 for the Gln/Gln genotype compared with the Arg/Arg genotype. For each study, the estimate of OR and its 95% CI is plotted with a box and a horizontal line. ⧫, pooled OR and its 95% CI.
Summary ORs for various contrasts of the XRCC1 codon 399 and codon 194 polymorphisms and cancer risk
XRCC1 Arg194Trp. The eligible studies included 4,933 cancer patients and 6,775 controls for this locus. Because of the low frequency of the variant allele in some ethnic subgroups, we only did the analysis with the dominant model of the Trp194 allele. Overall, the Trp194 allele was 31.2% (95% CI, 29.6-32.8) among Asian controls, which was significantly higher than that in Caucasians (6.6%; 95% CI, 5.9-7.4) and in African population (7.3%; 95% CI, 5.7-9.2; P < 0.0001; Fig. 1). Overall, a significantly decreased risk was associated with the variant genotypes (Trp/Trp + Arg/Trp), compared with the wild homozygote Arg/Arg genotype (OR, 0.89; 95% CI, 0.81-0.98; Fig. 3). In stratified analyses, however, no significant associations were found in populations with either different ethnicity or different tumor types (Table 2).
ORs (log scale) of cancer associated with XRCC1 codon 194 for the Trp/Trp and Arg/Trp genotypes compared with the Arg/Arg genotype. For each study, the estimate of OR and its 95% CI is plotted with a box and a horizontal line. ⧫, pooled OR and its 95% CI.
XRCC1 Arg280His. Because only eight studies have investigated the Arg280His polymorphism and cancer risk to date, we did not perform stratification analysis. The eight studies included 1,688 cancer patients and 2,129 control subjects. Similar to the Arg194Trp, we did the analysis only with the dominant model of the His280 allele because of its low frequency. Overall, individuals with the variant genotypes (His/His + Arg/His) had a borderline significantly increased cancer risk, compared with individuals with the Arg/Arg genotype (OR, 1.19; 95% CI, 1.00-1.42) without between-study heterogeneity (Fig. 4).
ORs (log scale) of cancer associated with XRCC1 codon 280 for the His/His and Arg/His genotypes compared with the Arg/Arg genotype. For each study, the estimate of OR and its 95% CI is plotted with a box and a horizontal line. ⧫, pooled OR and its 95% CI.
Test of Heterogeneity
There was substantial heterogeneity among the 37 studies that included the Arg399Gln polymorphism (χ2 = 54.46, df = 36, P = 0.02) but not among the 22 studies that included the Arg194Trp polymorphism (χ2 = 25.50, df = 21, P = 0.23) and other 8 studies that included the Arg280His polymorphism (χ2 = 10.12, df = 7, P = 0.18). Therefore, we evaluated the source of heterogeneity for the Gln/Gln genotype (Gln/Gln versus Arg/Arg) by study sample size, ethnicity, and tumor type. When we dichotomized the 37 studies by the sample size cutoff value of 400 lung cancer cases, the difference of study sample size did not significantly contribute to the observed heterogeneity (χ2 = 0.19, df = 1, P = 0.66), nor did ethnicity (χ2 = 3.41, df = 3, P = 0.33) and tumor type (χ2 = 5.79, df = 5, P = 0.33).
Bias Diagnostics
XRCC1 Arg399Gln. The magnitude of the summary ORs had been fluctuating at around 1 in the past years (in random effect model, summary OR for Gln/Gln versus Arg/Arg: 1.03 at the end of 2001, 0.98 at the end of 2002, 1.04 at the end of 2003, and 1.03 till now). In the funnel plot analysis of publication bias (contrast of homozygous genotype plotted against the precision), the shape of the funnel plot seems asymmetrical, suggesting that the larger estimate of the association belongs to the smaller study, whereas the larger study shows little effect (Fig. 5). Furthermore, an Egger's test was used to provide statistical evidence for funnel plot symmetry (62). In the linear regression analysis, the intercept value provides a measure of asymmetry—the larger its deviation from zero, the more pronounced the asymmetry. We observed the intercept values of 1.57, which significantly deviated from zero (t = 3.85, P < 0.001).
Funnel plot analysis to detect publication bias. Each point represents a separate study for the indicated association. For each study, the OR is plotted on a logarithmic scale against the precision (the reciprocal of the SE). If bias is absent, small studies will have ORs that are widely scattered but still centered around the OR estimates provided by large, more precise studies.
XRCC1 Arg194Trp. The magnitude of the summary ORs had been undergoing a trend toward a more and more stable effect as postulated (summary ORs for Trp/Trp + Arg/Trp versus Arg/Arg: 0.75 at the end of 2001, 0.86 at the end of 2002, 0.88 at the end of 2003, and 0.89 till now). The change of the summary OR was in dispersion in the last 2 years. Although the funnel plot analyses and Egger's test suggest that publication bias may influence the results (intercept values 1.31, P = 0.004), indeed it reinforces the postulated effect because the smaller studies suggest a risk effect on cancer with the variant genotypes (Fig. 5).
XRCC1 Arg280His. In cumulative meta-analysis and recursive cumulative meta-analysis, the summary ORs changed considerably in the year 2004 with a sharp reduction of the postulated effect (summary ORs for His/His + Arg/His versus Arg/Arg: 1.36 at the end of 2001, 1.38 at the end of 2002, 1.40 at the end of 2003, and 1.19 till now). This change mainly resulted from the study of esophageal cancer by Hao et al. (52), which is the only study that presented a protective effect of the variant genotypes. Furthermore, the funnel plot and Egger's test revealed that publication bias might have influenced the estimates (intercept values 2.15, P = 0.002).
Discussion
This meta-analysis, including a total of 11,957 cancer cases and 14,174 controls from 38 case-control studies, examined the association of three well-characterized polymorphisms of the DNA repair gene XRCC1 (Arg194Trp, Arg280His, and Arg399Gln) with cancer risk. There was no overall effect either in recessive or dominant modeling for the Arg399Gln polymorphism, and compared with the wild Arg/Arg genotype, variant Gln/Gln homozygote was not associated with overall cancer risk. For the Arg194Trp, the variant genotypes (Trp/Trp + Arg/Trp), compared with the wild-type homozygote (Arg/Arg), were associated with a significantly decreased cancer risk (OR, 0.89; 95% CI, 0.81-0.98) for all tumor types without between-study heterogeneity. Similarly, we observed a fixed overall 19% increased risk of cancer for the variant genotypes (His/His + Arg/His) of the Arg280His polymorphism, compared with the wild homozygote (OR, 1.19; 95% CI, 1.00-1.42). However, considering the relatively small sample size and marginal statistical evidence for Arg280His, our result in relation to this polymorphism should always be treated as preliminary. Nevertheless, our analysis shows that even if a common variant in the functional region of a definitively meaningful gene had an effect on human disease, such as cancer, it may play only a small role in the disease causation, which is consistent with the characteristics of low-penetrance genes (63).
Given the multiplicity of possible comparisons and the unavoidable flexibility of choosing and defining the correlates, associations may have been detected by chance alone. Several criteria have been proposed for assessing associations between genetic polymorphisms and disease (64); the claim was that studies “ideally should have large sample sizes, small P values, report associations that make biological sense, and alleles that affect the gene product in a physiologically meaningful way.” The sample size and scientific hypotheses of the study are obviously important to know the proportion of false-positive findings of meta-analysis that are attributable to constituent studies with poor study design, nondifferential misclassification errors, and selection bias from publication (65). Although considerable effort and resources have been put into testing possible associations between DNA repair gene XRCC1 polymorphisms and cancer risk, there are still serious limitations inherited from the published studies. First, selection bias could have played a role because the genotype distributions of at least one polymorphism among control subjects disobeyed the law of Hardy-Weinberg equilibrium in five studies (20, 21, 27, 32, 44). Second, demographic parameters are not well matched and are statistically adjusted in a few studies (26, 32, 36, 39, 44, 49). Third, misclassifications on disease status and genotypes may also influence the results because cases in several studies were not confirmed by pathology or other gold standard methods, and the quality control of genotyping was also not well-documented in some studies. Finally, the Egger's test revealed that there are considerable unpublished negative studies not included in this meta-analysis although they are likely to be small in terms of sample size. However, because the Trp194 allele was associated with a decreased cancer risk and the unpublished studies on this polymorphism are likely to show a protective effect, the association between the Arg194Trp variant and cancer risk may be slightly underestimated. To further evaluate the results derived from Egger's test, we searched the abstracts from annual meetings of AACR between 2000 and 2004 and found 24 eligible case-control studies. Among the 24 abstracts, 11 were fully or partially published and had been included in the current meta-analysis and only 2 unpublished studies were with a sample size of >300 cases. One study conducted in Finland with 483 breast cancer patients and 482 controls reported that no association between the Arg399Gln polymorphism and breast cancer risk (66). The other multicentric study in central and eastern Europe consisted of a large sample size (2,073 cases and 1,953 controls) showed that Arg280His, Arg399Gln did not confer an effect on lung cancer, whereas Arg194Trp showed a protective effect among current smokers (OR, 0.76; 95% CI, 0.42-0.99; ref. 67). These two meeting abstracts with relatively large sample size also support the current meta-analysis that XRCC1 Arg194Trp rather than Arg280His and Arg399Gln may play a role in individual susceptibility to cancers.
Several studies suggested a possible interaction between XRCC1 Arg399Gln and family history on breast cancer risk (48), indicating that family history, particularly in first-degree relatives, broadly represents shared genes and environmental factors, and that the weak effect of a single polymorphism on the individual's phenotype may not be measurable except in the context of these risk factors. In addition, cigarette smoking (36), alcohol consumption (35), and antioxidant vitamins intake (43) are associated with the production of free radical intermediates, including hydroxyethyl free radicals and reactive oxygen species, which are corrected in part by the involvement of XRCC1. Therefore, gene-to-environment interactions have been of great interest to evaluate the exact roles of genetic polymorphisms. However, lacking of the original data of the meta-analysis limited our further evaluation of potential gene-to-gene and gene-to-environment interactions.
The heterogeneity test revealed that there was no significant between-study heterogeneity in terms of the XRCC1 Arg194Trp polymorphism for all tumor types. In addition, different tumor types also did not significantly contribute to the overall heterogeneity in relation to the Arg399Gln polymorphism, indicating that our current combined analyses were unbiased, regardless of tumor types. Because we did not find any significant heterogeneity in the distribution of the Arg399Gln genotypes by ethnicity and the study sample size, it is possible that other unmeasured characteristics in different study populations and/or inherited limitations of the recruited studies may partially contribute to the observed overall heterogeneity. Although heterogeneity among different ethnic groups was not statistically significant in our combined analysis, other population stratification factors may have played a role in the heterogeneity when the sample size of a combined analysis become sufficient large (68) and, therefore, we suggest a careful matching on ethnicity for future larger genetic association studies.
In conclusion, our current study support that XRCC1 Arg194Trp polymorphism may contribute to individual susceptibility of cancers. To further evaluate gene-to-gene and gene-to-environment interactions on XRCC1 polymorphisms and cancer risk, a single larger study with thousands of subjects and tissue-specific biochemical and biological characterizations are required.
Footnotes
-
Grant support: National Key Basic Research Program grants 2002CB512900 and National Natural Science Foundation grant 30371240.
-
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
- Accepted April 18, 2005.
- Received October 28, 2004.
- Revision received March 24, 2005.