Abstract
Poly (ADP-ribose) polymerase 1 (PARP1), which plays a critical role in the base excision DNA repair mechanism, and matrix metallopeptidase 2 (MMP2), a member of the matrix metalloprotease family, are involved in tumor formation and metastasis, respectively. In the present study, the possible association of single nucleotide polymorphisms (SNPs) and gene-gene interaction between PARP1 and MMP2 with the increased incidence of gastric cancer (GC) development and lymph node metastasis (LNM) was investigated in a Korean population. Samples were obtained from 326 patients with chronic gastritis and 153 patients with GC and genotyped using the GoldenGate® method. The PARP1 rs1136410 genotype showed a significant association with the frequency of LNM of GC (odds ratio [OR] = 2.19, p = 0.02), LNM stage (p = 0.035), and tumor invasion (p = 0.035). The allele frequency of MMP2 rs243865 was not associated with the development of GC or with the development of LNM of GC. Epistasis between the PARP1 SNP and the MMP2 SNP was associated with the development of LNM of GC. The combination of the MMP2 rs243865 CC genotype and the PARP1 rs1136410 CC or CC+CT genotypes showed a high risk of LNM of GC (OR = 2.47, p = 0.01; OR = 2.28, p = 0.01, respectively). In summary, PARP1 is associated with the risk of LNM of GC and the stage of LNM and tumor invasion. Epistasis between PARP1 rs1136410 and MMP2 rs243865 increased the risk of LNM of GC.
Introduction
The development of lymph node metastasis (LNM) in GC is a critical factor in the determination of the prognosis and survival after treatment (Coburn, 2009). The most important step in metastasis and invasion is the break down and penetration of the basement membrane (Schwartz, 1996). The matrix metallopeptidase (MMP) family plays a critical role in cancer invasion and metastasis through the degradation of extracellular matrix or basement membrane barriers (Curran and Murray, 1999; Stamenkovic, 2000). Among the MMP family, MMP2 hydrolyzes type IV collagen (Stamenkovic, 2000) and its expression is associated with LNM stage in GC (Monig, 2001; Wu et al., 2006). In particular, the frequency of the MMP2 rs243865 genotype is associated with the incidence of nasopharyngeal cancer (Zhou et al., 2007), GC (Miao et al., 2003), and LNM from GC (Wu et al., 2007). The MMP2 rs243865 single nucleotide polymorphism (SNP) is located in the promoter region, and the C to T transition is associated with reduced MMP2 gene expression levels (Price et al., 2001).
Considering interaction of polymorphisms in more than two genes is important to explain the inheritance of polygenic diseases. Recent studies revealed that gene-gene interactions show a stronger association with increased cancer risk than individual SNPs (Zhang et al., 2005; Chen et al., 2009). The present study investigated the potential association of PARP1, which is involved in the BER pathway, and MMP2, which is involved in metastasis, with GC development and metastasis. In addition, the interactions between these two proteins and their association with GC and LNM were assessed through a case-control study.
Materials and Methods
Subjects for the case-control study and sample preparation
A total of 326 patients with chronic gastritis (CG) and 153 patients with GC were recruited from the division of gastroenterology of the Ajou University hospital. The patients with CG were diagnosed by endoscopic examination between January 2002 and December 2006; and only patients with H. pylori-positive CG participated in this study. Patients with GC who were all operative cases for the eradication of a gastric tumor after being diagnosed by endoscopic examination and biopsy between January 2001 and December 2008 and who were H. pylori-positive were included in the study. H. pylori infection was detected using the Campylobacter-like organism test and the 13C-urea breath test.
TNM stage of patients with GC was determined using isolated tumor mass and lymph nodes according to the 6th edition of the GC staging system, which was devised by an American Joint Committee on Cancer (AJCC) (Greene et al., 2002). The TNM classification by AJCC sixth was decided by a surgeon in the division of gastroenterology of the Ajou University Hospital where patients with GC were treated.
Genomic DNA (gDNA) was extracted from blood samples using G-DEX™ blood gDNA purification kits (Intron Biotechnology, Seungnam, Korea). Purified DNA was quantified by using the picogreen dsDNA quantification reagent according to a standard protocol (Molecular Probes, Eugene, OR). The gDNA was then stored at −20°C until genotypic analysis was performed. Genotypic analysis was performed by using the GoldenGate® Assay Kit (Illumina, San Diego, CA). Oligos was amplified by allele-specific primer extension. After hybridization to a sentrix array matrix, signal intensities were read by BeadArray Reader (Illumina). BeadStudio III software (Illumina) was used for the clustering of genotype data.
Statistical analysis
The associations between genotype and risk of GC were estimated by odds ratio (OR), and the 95% confidence interval was calculated by logistic regression adjusting for age and sex. Genotype distributions were compared between the CG and LNM of GC groups using the same statistical analysis. The genetic models were divided into additive (AA vs. AB vs. BB), dominant (AA vs. AB plus BB), and recessive (AA plus AB vs. BB). Associations between the PARP1 rs1136410 genotype and clinicopathological features were evaluated by using Fisher's exact test. Logistic regression and Fisher's exact test were performed by using the SAS (SAS Institute, Cary, NC; version 9.1.3). Epistasis analyses were performed by logistic regression adjusting for age and gender using the PLINK program. (
Results
The characteristics of the 326 CG and 153 GC patients are summarized in Table 1. The age and sex distribution among patients with GC and CG was statistically different. The mean ages of the patients with GC and CG were 57.8 years and 55.2 years, respectively. Among the patients with GC, 61.4%, 20.9%, 14.4%, and 3.3% were classified as N0, N1, N2, and N3, respectively, based on LNM stage; and 45.1%, 24.2%, 28.8%, and 2.0% were classified as T1, T2, T3, and T4, respectively, based on tumor invasion stage.
The p-value was calculated by Chi-square test.
LNM and tumor invasion were determined according to the sixth edition of the gastric cancer staging system, which was devised by an American joint committee on cancer.
LNM, lymph node metastasis.
To study the association of PARP1 and MMP2 with the risk of GC, the genotype distribution of well-characterized SNPs (PARP1 rs1136410 and MMP2 rs243865) from both genes was determined (Table 2). The distribution of genotypes satisfied the Hardy-Weinberg equilibrium (p > 0.05) by χ 2 test (data not shown). The frequencies of the PARP1 rs1136410 TT, TC, and CC genotypes were 31.9%, 50.3%, and 17.8% among patients with CG, and 27.8%, 46.4%, and 25.8% among patients with GC, respectively. The frequencies of the MMP2 rs243865 CC, CT, and TT genotypes were 83.1%, 16.6%, and 0.3%, respectively, among patients with CG, and 78.4%, 20.9%, and 0.7% in patients with GC, respectively. The distribution of the rs1136410 and rs243865 SNPs was not significantly associated with the risk of GC development.
Six patients with chronic gastritis and two patients with gastric cancer were excluded from the analysis, because the PARP1 genotype could not be defined.
The p-value was calculated by logistic regression.
Power was calculated by Power and Precision V4 software.
OR, odds ratio; CI, confidence interval; PARP1, poly (ADP-ribose) polymerase 1; MMP2, matrix metallopeptidase.
The 153 patients with GC were subdivided into 59 patients with LNM and 94 patients without LNM. The genotype distribution of the two polymorphisms (PARP1 rs1136410 and MMP2 rs243865) was compared between patients with CG and LNM of patients with GC to assess for the risk of the development of LNM (Table 3). The frequencies of the PARP1 rs1136410 TT, TC, and CC genotypes were 28.8%, 39.0%, and 32.2% in patients with LNM GC. The PARP1 rs1136410 CC genotype was significantly associated with an increased risk of LNM in patients with GC (OR = 2.19, p = 0.02). The frequencies of the MMP2 rs243865 CC, CT, and TT genotypes were 83.1%, 15.3%, and 1.7% in patients with LNM of GC. For the MMP2 rs243865 SNP, no significant deviation was observed in the genotype distribution of the polymorphism between CG and LNM of GC.
The p-value was calculated by logistic regression.
Power was calculated by Power and Precision V4 software.
Bold face indicates significant p-value (p<0.05).
The association between the clinicopathological features of GC and the distribution of the PARP1 rs1136410 polymorphism was assessed (Table 4). A significant correlation between the PARP1 rs1136410 genotype and tumor invasion stage (p = 0.035) and LNM stage (p = 0.035) was found in the GC patient group. The risk allele CC genotype frequency was significantly higher in the T3 and T4 subgroups and in the N1 and N3 subgroups.
The p-value was calculated by Fisher's exact test.
Bold face indicates significant p-value (p<0.05).
Gene-gene interactions between PARP1 and MMP2 were also investigated (Table 5). A significant association between LNM GC risk and the combined effects of the rs243865 and rs1136410 SNPs was found. The MMP2 rs243865 CC and PARP1 rs1136410 CC polymorphism combination contributed to an increased risk of LNM among patients with GC (OR = 2.47, p = 0.010). The combination variant rs243865 CC+CT and rs1136410 CC exhibited a 2.28-fold increase in the risk for LNM development (OR = 2.28, p = 0.012). The p-value after correction by Bonferroni's multiple testing showed a significant association between gene-gene interaction and LNM of GC (p-value <0.0167). Power values of the interaction in a dominant manner and in a recessive manner were 0.76 and 0.68, respectively.
The p-value less than 0.0167 was considered significant after Bonferroni correction.
Power was calculated by Power and Precision V4 software.
GC, gastric cancer; CG, chronic gastritis.
Bold face indicates significant p-value (p<0.05).
Discussion
The present study investigated the interaction between PARP1 rs1136410 and MMP2 rs243865 and the association between these two SNPs and the incidence of GC and LNM of GC using logistic regression analyses adjusted for age and gender. The PARP1 rs1136410 genotype was significantly associated with the development of LNM of GC, LNM stage, and tumor invasion, but not with GC development. The MMP2 rs243865 polymorphism was not associated with GC susceptibility or LNM of GC. Interestingly, the interaction between these two SNPs was found to increase the risk of LNM of GC more than each SNP alone.
H. pylori infection is one of the most important environmental factors involved in the pathogenesis of GC. Patients with H. pylori-infected CG and GC were included in the present study to rule out the effect of H. pylori infection and to emphasize the effect of genetic factors. Uemura et al. (2001) reported an association between H. pylori infection and the incidence of GC. H. pylori induces cell proliferation, ROS formation, and causes DNA damage in gastric epithelial cells (Obst et al., 2000); and infection with H. pylori results in the down-regulation of the BER pathway in vitro and in vivo (Machado et al., 2009).
The BER pathway is involved in the repair of DNA damaged by ROS generated by H. pylori infection, ionizing radiation, or alkylating agents (Hoeijmakers, 2001). Various genes such as APEX nuclease (multifunctional DNA repair enzyme) 1, polymerase (DNA directed) beta, XRCC1, NEIL1, and PARP1 are involved in the BER pathway. These genes play a role in tumor suppression mechanisms (Sweasy et al., 2006). Polymorphisms of genes that are part of the BER pathway have frequently been associated with various cancers (Mitra et al., 2008; Sliwinski et al., 2009; Goto, 2010). In particular, PARP1 plays a critical role in the BER pathway by forming a functional complex involved in DNA repair. PARP1 is located on chromosome 1 and consists of a DNA-binding domain, an automodification domain, and a catalytic domain. PARP1 enzymatic activity is affected by the T to C substitution in the rs1136410 polymorphism, which results in an amino acid change from Val to Ala. The V to A substitution increases the Km, which results in a decrease in enzyme activity (Mitra et al., 2008). PARP1 rs1136410 CC type was found to significantly decrease enzyme activity in prostate cancer in an allele dosage-dependent manner (Lockett et al., 2004). Zhang et al. (2005, 2009) reported an association between PARP1 rs1136410 and GC risk and lung cancer risk in the Chinese population.
In the present study, PARP1 rs1136410 was significantly associated with patients with CG and LNM of patients with GC (Table 2, OR = 2.19, p = 0.02). The PARP1 rs1136410 CC genotype was significantly associated with the LNM stage (p = 0.035) and tumor invasion stage (p = 0.035) in the GC patient group by Fisher's exact test, which is appropriate for small sample sizes (Table 4). The PARP1 rs1136410 CC genotype correlated with the development of LNM of GC and cancer progression, but did not influence GC development.
In addition, the PARP1 rs1136410 SNP showed almost perfect linkage disequilibrium (LD) with rs1805415 (r2 = 0.998), as revealed by Korean LD data obtained from the Korea association resource dataset (Cho et al., 2009) merged with imputed data by using the IMPUTE program (Marchini et al., 2007). Cao et al. (2007) reported a perfect LD between two SNPs (rs1136410 and rs1805415) in a French population. The minor A allele of PARP1 rs1805415 is associated with increased interleukin-6 (IL-6) levels (Walston et al., 2009), and high IL-6 levels were associated with LNM stage (Kim et al., 2009). Among the two SNPs (rs1136410 and rs1805415) that have the two haplotypes (Hts) TG and CA, the CA Ht has a lower DNA repair capacity and higher IL-6 level than the TG Ht. Therefore, the PARP1 rs1136410 and rs1805415 CA Ht might be associated with the development of LNM of GC.
MMP2 is a critical protein for cancer progression due to its role in metastasis and invasion. The MMP2 rs243865 SNP is located on the promoter sequence, and its C to T transition disrupts the SP1-type promoter site (CCACC box), resulting in lower promoter activity (Price et al., 2001). The MMP2 rs243865 polymorphism is associated with various diseases, including nasopharyngeal carcinoma (Zhou et al., 2007), esophageal cancer (Yu et al., 2004), lung cancer (Yu et al., 2002), oral carcinoma (Lin et al., 2004), GC (Miao et al., 2003), and lymphatic invasion of GC (Wu et al., 2007). The MMP2 rs243865 polymorphism is also associated with stroke in the Portuguese population (Manso et al., 2010) and rheumatoid arthritis among the Spanish population (Rodriguez-Lopez et al., 2006). However, certain studies reported that the MMP2 rs243865 polymorphism is not associated with intestinal metaplasia (Zhu et al., 2009), breast cancer, (Beeghly-Fadiel et al., 2009), or melanoma (Cotignola et al., 2007).
In the present study, MMP2 rs243865 genotype frequency was not significantly associated with GC development among patients with CG and GC, and it was not associated with LNM in patients with GC. The MMP2 rs243865 T minor allele was found at a frequency of 9% in CG, 11% of GC according to the present data. MMP2 polymorphisms might not be associated with the development of GC or LNM of GC in the Korean population. Further studies, including a larger number of samples, are necessary to clarify the involvement of MMP2 polymorphisms in the development of GC and LNM in GC. Additionally, a more detailed LNM staging system such as lymphatic invasion designation might be useful.
Gene-gene interactions are important in polygenic diseases such as cancers. The interaction between the PARP1 and XRCC1 genes, which are involved in the BER pathway, is associated with an increased risk of lung cancer (Zhang et al., 2005) and LNM in thyroid carcinoma (Chiang et al., 2008). Interaction between the MMP9 and TIMP3 genes, which play important roles in metastasis, is associated with breast cancer risk (Lei et al., 2007). However, there are no studies investigating the interactions between two different but functionally related pathways such as the BER pathway and the LNM process. In the present study, the interaction between PARP1 rs1136410 and MMP2 rs243865 was significantly associated with LNM of GC. In the MMP2 rs243865 CC group, the PARP1 CC genotype was significantly associated with an increase in the risk of LNM from GC. In the MMP2 rs243865 CC+CT group, the PARP1 CC genotype significantly increased the risk of LNM from GC. Notably, the combination of polymorphisms from two different genes affected the incidence of GC formation and metastasis of GC. The statistical powers of this interaction in a dominant and recessive manner were 0.76 and 0.68, respectively (Table 5). Recent association reports based on similar sample sizes to that of this study were found to have similar or lower power than that in this study (Supplementary Table S1; Supplementary Data are available online at
Although functional validation of PARP1 gene and MMP2 gene interaction will require a bigger sample size, speculation on the functional role of SNPs in the genes involved in BER and metastasis pathways might be helpful to design future experiments. The PARP1 rs1136410 CC genotype might be characterized by reduced DNA repair capacity and enhanced oncogenic capacity. The PARP1 rs1805415 A allele, which is associated with elevated levels of IL-6, might promote the LNM of GC. The MMP2 rs243865 C allele, which could enhance the activity of the promoter resulting in increased MMP2 gene expression, could result in an elevated level of extracellular matrix degradation. Koreans who have been infected with H. pylori and have the PARP1 rs1136410 CC genotype and the rs1805415 A and MMP2 rs243865 C alleles might have an elevated risk of LNM of GC. Therefore, the interaction between PARP1 and MMP2 may be a valuable cancer prognosis marker.
Conclusion
In conclusion, the present case-control study of H. pylori infection-associated GC showed a significantly increased risk of LNM of GC, LNM stage, and tumor invasion stage in patients with PARP1 CC type SNPs. In addition, the present results indicated that the risk of LNM of GC was associated with gene-gene interactions between MMP2 rs243865 CC+CT type and PARP1 CC type SNPs in the Korean population.
Footnotes
Acknowledgments
We are thankful to every individual who gave us informed consent for this study. Informed consent was acquired from each patient, and the sampling protocols and informed consent forms used were approved by the IRB of Ajou University. Protocols for handling the human genomic materials and all the experimental procedures used were approved by the IRB of CHA University. This work was funded by grants from the Ministry of Health and Welfare, Republic of Korea (A010383, A080734).
Disclosure Statement
No competing financial interests exist.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
