| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Epidemiology |
1 Lung Cancer Program, Lovelace Respiratory Research Institute; 2 Department of Internal Medicine, University of New Mexico; 3 New Mexico VA Health Care System, Albuquerque, New Mexico and 4 Keck School of Medicine, University of Southern California, Los Angeles, California
Requests for reprints: Steven A. Belinsky, Lung Cancer Program, Lovelace Respiratory Research Institute, 2425 Ridgecrest Drive Southeast, Albuquerque, NM 87108. Phone: 505-348-9465; Fax: 505-348-4990; E-mail: sbelinsk{at}LRRI.org.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
30% of all deaths from cancer (1). The lack of a validated screening approach for early detection and the resistance of advanced-stage tumors to therapy are largely responsible for the 5-year survival rate of 14% for lung cancer patients (2). The discovery of field cancerization in the respiratory tract of smokers prompted studies leading to the discovery that inactivation of genes, such as p16, by promoter hypermethylation occurs in precursor lesions to non–small cell lung cancer (3). This finding suggested that methylation, when detected in exfoliated cells within sputum, could serve as a biomarker for the early stages of lung carcinogenesis (4). To test this hypothesis, our group examined a large panel of genes for their ability to predict lung cancer in a nested case-control study. A combination of six genes was identified whose methylation in sputum predicted lung cancer before clinical diagnosis with both a sensitivity and specificity of 65% (5). A better understanding of the susceptibility factors that predispose smokers to the acquisition of multiple epigenetic alterations in key cellular regulatory genes within the respiratory epithelium could improve prediction of lung cancer risk and affect strategies for early detection and chemoprevention.
The precise mechanisms by which carcinogens disrupt the capacity of the cells to maintain the normal epigenetic code during DNA replication and repair are largely unknown. Smoking accounts for >90% of lung cancer. Carcinogens within tobacco induce single- and double-strand breaks (DSB) in DNA, and reduced capacity for repair of DNA damage has been associated with lung cancer (6). Accumulating evidence suggests that extensive DNA damage, manifested through DSBs, could in part be responsible for the acquisition of aberrant gene promoter methylation during lung carcinogenesis. For example, the prevalence of promoter methylation of the p16 gene was significantly greater in adenocarcinomas from workers occupationally exposed to plutonium, an exposure that predominantly produces DSBs, than in cancer from unexposed smokers (7). The prevalence of p16 methylation increased with increasing plutonium exposure. In a second study, the prevalence of methylation of the estrogen receptor
gene promoter was greater in plutonium-induced adenocarcinomas in rodent lung tumors compared with tumors induced by 4-(methylnitrosamino)-I-(3-pyridyl)-1-butanone, diesel exhaust, or carbon black exposures, which mainly induce single-strand breaks of DNA (8). These studies support the hypothesis that DSBs may play a key role in the development of aberrant gene promoter hypermethylation. The purpose of this study was to test the hypothesis that a high methylation index (defined as the methylation of three or more gene-specific promoters detected in sputum) is associated with a reduced capacity to repair DSBs. We also hypothesize that sequence variation in genes from the DSB repair pathway will predict for high methylation index.
| Materials and Methods |
|---|
|
|
|---|
Cytologically adequate sputum samples from 824 cohort subjects were evaluated for gene promoter methylation of eight genes as described below. High methylation index was defined as the methylation of three or more gene-specific promoters in sputum. We selected persons from our cohort that exhibited a high (cases) or low [controls (zero of eight genes)] methylation index. To increase the stringency for case selection, GATA4, which was most commonly methylated in sputum, was excluded as one of the three methylated genes needed for case classification and 131 of 824 cohort subjects met this criteria. Cases were frequency matched by gender to controls. Cases (n = 131) and controls (n = 130) were selected for the genetic association study. Among the 131 cases, 77 had adequate number of cryopreserved lymphocytes for the mutagen sensitivity assay. Seventy-eight controls were selected from the 130 controls, with frequency matching by gender maintained, for the mutagen sensitivity assay.
Sputum cytology and nested methylation-specific PCR. Sputum samples were stored in Saccomanno's fixative. Three slides were made for each sputum sample to check for adequacy defined as the presence of deep lung macrophages or Curschmann's spiral (9). The methylation-specific PCR (MSP) assay was only performed on cytologically adequate sputum samples. Eight genes [p16, O6-methylguanine-DNA methyltransferase (MGMT), death-associated protein kinase (DAPK), ras effector homologue 1 (RASSF1A), PAX5
, PAX5β, GATA4, and GATA5] were selected for analysis of methylation in sputum based on our previous studies establishing their association with risk for lung cancer (5, 10–12). Nested MSP was used to detect methylated alleles in DNA recovered from the sputum samples as described (5, 10–12).
Evaluation of DSB repair capacity in peripheral lymphocytes. Phytohemagglutinin (PHA)-stimulated lymphocytes were treated with bleomycin to evaluate the generation of chromosome aberrations as an index of DSB repair capacity (DSBRC; ref. 13). Briefly, cryopreserved lymphocytes were thawed and cultured in RPMI 1640 supplemented with fetal bovine serum (20%) and PHA (1.5%) at a cell density of <0.5 x 106/mL. Sixty-seven hours after PHA stimulation, the culture was split into two T25 flasks and treated with bleomycin or vehicle for 5 h. The final concentration for bleomycin in culture medium was 3 units/L, a concentration defined through dose-response studies using isolated lymphocytes from cohort subjects and two lymphoblastoid cell lines: GM02782 [mutant ataxia-telangiectasia mutated (ATM)] and GM00131 (wild-type ATM; data not shown). The dose selected was within the linear dose-response range and caused obvious genotoxicity but minimal cytotoxicity. One hour before harvest, colcemid was added to the cultures at a final concentration of 0.06 mg/mL. Slides were prepared according to conventional procedures and 100 well-spread metaphases were examined for chromatid breaks. Samples were assayed as a batch, and slides were scored by a person blinded to case-control status. The criteria of Hsu et al. (13) were used to record the aberrations: a chromatid break was scored as one break and each isochromatid break set was scored as two breaks. Chromosome/chromatid gaps, chromosome-type aberrations (dicentrics, ring, and acentric fragments), or chromatid exchanges were recorded but not added to the frequencies of chromatid breaks. On rare occasions, a metaphase with >12 breaks was observed on a slide with bleomycin treatment. When this occurred, the number of breaks was recorded as 12. The DSBRC was expressed as the mean number of chromatid breaks per cell.
The means of spontaneous chromatid breaks per cell derived from 100 metaphases of untreated cells were 0.013 in cases and 0.021 in controls, which were similar to the spontaneous frequency reported in the literature (13), and less than one fifteenth the mean number of breaks seen in bleomycin-treated cells (0.32). Therefore, for statistical comparisons, the spontaneous breaks were not subtracted from the breaks observed following treatment with bleomycin.
Single nucleotide polymorphism selection and genotyping by Illumina platform for 16 genes in the DSB repair pathway and related cell cycle control genes. A total of 294 single nucleotide polymorphisms (SNP) were selected for 16 candidate genes from DSBR and cell cycle control pathways. Tag SNPs (n = 245) were derived from Latino and White data from the University of Southern California plus phase 1 HapMap for whites for 15 genes. Tag SNPs were selected using r2
0.8, with nonsynonymous SNPs retained as the tag SNPs (14). One additional SNP for bins with at least six or more SNPs was selected as a redundant SNP in case of genotyping failure. For the remaining gene, NBN, 49 SNPs were selected using dbSNPs based on a SNP density of one to three SNPs per kb depending on the haplotype block structure, validation status, Illumina design score, and functional potential of the SNPs. The number of SNPs selected for each of these 16 genes is shown in Supplementary Table S1, and a SNP list is available on request. These SNPs were genotyped by the Illumina Golden Gate Assay for 261 DNA samples isolated from lymphocytes of cases and controls.
Selection of subjects and construction of MRE11A promoter constructs. Five common haplotypes (6–34%) were constructed based on the 14 tag SNPs assayed for MRE11A in the population (15). A Bayesian statistical method implemented in the program PHASE (version 2.1) was used to reconstruct the haplotypes from the SNPs in the MRE11A gene for the 261 subjects. Two subjects homozygous for the haplotype that contained the rs7117042 SNP associated with high methylation index were selected. The other four people selected were each homozygous for one of the other four haplotypes. The MRE11A promoter fragment (–2,541 to –5 with +1 being the translational start site) was amplified from lymphocyte DNA from these six subjects. The promoter fragment was directionally subcloned into the pGL2-basic luciferase reporter vector (Promega) upstream of the luciferase coding sequence. Five clones from each person were commercially sequenced to identify variants within the promoter region (Sequetech).
Transient transfection and reporter gene assays. The Calu6 lung tumor–derived cell line was used for transient transfections. Cells (1.5 x 105) were plated into six-well dishes and transfected the following day. Plasmid DNA (1 µg) and the pSV-β-galactosidase control vector (0.5 µg; Promega) were cotransfected into cells with Fugene 6 transfection reagent (Roche Diagnostics) at a Fugene 6 to DNA ratio of 3:1. A promoterless pGL2-basic vector and the pGL2-control vector that contains the SV40 promoter were used as negative and positive controls, respectively. Forty-eight hours after transfection, cells were harvested and lysed. Immediately after lysing, cell extracts were assayed in a luminometer for luciferase activity using the Luminoskan Ascent luminometer (Thermo Electron) for luciferase activity using the Luciferase Assay System (Promega). β-Galactosidase activity in cell lysates was measured using the Galacto-Star Reporter Gene Assay System (Tropix). Promoter activity was calculated as the ratio of activities of luciferase and β-galactosidase. Transfections were done in duplicate in four independent experiments.
Statistical analysis. The two-sample t test, Wilcoxon rank sum test, and
2 test were used to compare the mean or distribution of several demographic variables and DSBRC results between cases and controls as appropriate. Because the DSBRC data and the number of spontaneous breaks were not normally distributed, analysis was also performed on log-transformed data. The results based on log-transformed data were similar to those based on untransformed data, so only results based on untransformed data are shown. Analysis of covariance and logistic regression were used to assess the association between selected variables, such as SNPs and case-control status, and the outcome variable, DSBRC, with adjustment of covariates selected a priori (age at sputum collection, sex, race, current smoking status, and pack-years). DSBRC was dichotomized for logistic regression models using the upper quartile of DSBRC in control participants. The selection of the upper quartile of DSBRC in controls as the cutoff value was based on the distribution of DSBRC in cases and controls. Analysis of covariance and logistic regression models, stratified by status, were also examined for different associations between SNPs and DSBRC by case-control status. A receiver operator characteristic (ROC) curve was also drawn to compare the sensitivity and specificity of DSBRC induced by bleomycin for classifying cases (16). Multivariate unconditional logistic regression assessed the association between SNPs and the outcome of case-control status, with the same covariates outlined above. Model results are presented as odds ratios (OR) with 95% confidence intervals (95% CI) for having three or more methylated genes. Logistic regression modeling was extended to generalized logit models to more precisely examine the high methylation index. ORs and 95% CIs for the risk of having three, four, or five or more methylated genes with zero methylated genes as the reference group were obtained with adjustment for the same covariates.
The call rate for each SNP was assessed before data analysis. For the 294 SNPs assayed, 42 were deemed unsuitable because they were monomorphic, had a mean allele frequency of <0.05, had low yield (<80%), or showed a highly significant distortion from Hardy-Weinberg equilibrium (P < 0.0001). These SNPs were removed from analysis. Four models were tested: codominant, dominant, additive, and recessive. Because of power limitations, only results for the additive model are presented for each SNP, and common homozygote, heterozygote, and rare homozygote were coded as 0, 1, and 2, respectively. A logistic regression model was used to calculate the ORs and 95% CIs for each individual SNP with adjustment for age, sex, ethnicity, and smoking selected a priori. A ROC curve was drawn to evaluate the classification accuracy of this panel of variables for promoter methylation. An analysis excluding the 23% of study subjects that were not of non-Hispanic white origin had no effect on the identified associations. Therefore, all 261 subjects were included in the data analysis.
Monte Carlo estimates of exact P values were calculated by permuting the case-control status for all subjects 10,000 times. False-positive report probability (FPRP) was also calculated to address the robustness of our findings for individual SNPs (17). In assigning a prior probability for these genes, we considered the strong association between DSBRC and risk for promoter methylation and the stringent r2 value (0.8) for selecting tag SNPs. Based on the evidence for associations between SNPs in CHEK2, XRCC3, DNA-PKc, NBN, LIG4, and XRCC2 and several cancers (18–25), we assigned a relatively high prior probability range (0.1–0.25) for SNPs of these six genes. In contrast, for MRE11A, Ku80, RAD50, and CHEK1, a relatively low prior probability range (0.01–0.1) was assigned because there are no studies that have addressed the association of variants within these genes to cancer. All data analyses were performed with SAS/STAT and SAS/GENETICS 9.1.3.
| Results |
|---|
|
|
|---|
, and PAX5β was evaluated. Methylation of these genes has been associated with increased risk for lung cancer (5, 10–12). The prevalence of methylation ranged from 1.2% for RASSF1A to 31% for GATA4 and was not associated with family history for lung cancer (Supplementary Table S2). Nineteen percent of cohort members were methylated for three or more genes (Supplementary Table S2). Our previous nested case-control study within the Colorado Cohort revealed that methylation of three or more genes from a six-gene panel (excluding GATA4 and PAX5
) was associated with a 6.5-fold increased risk for lung cancer (5).
|
|
The OR and 95% CI were calculated to further characterize the association between methylation index and DSBRC. The OR associated with an increment of 0.1 chromatid breaks per cell for having three or more methylated genes was 6.6 (95% CI, 3.7–13.2) for bleomycin treatment after adjustment for selected covariates. When the number of methylated genes was used as a multinomial response variable, the OR associated with an increment of 0.1 chromatid breaks per cell for methylation of three, four, and five or more genes was 5.3 (95% CI, 2.7–11.4), 7.1 (95% CI, 3.8–14.8), and 8.5 (95% CI, 4.2–19.3), respectively. A histogram detailing the distribution of chromatid breaks per cell by case-control status revealed that 75% of controls compared with 18% of cases accumulated <0.38 breaks per cell (Fig. 1D). We chose to dichotomize the number of DSBs per cell at the 75th percentile for controls and this gave an adequate overlap in the distribution of breaks per cell in cases and controls. The results did not differ significantly with other cut points. When methylation index was compared with chromatid breaks per cell, the ORs for detecting methylation increased monotonically from 10 to 15 to 26 (Table 2
). Overall, chromatid breaks
0.38 per cell were associated with a 14.5-fold increased risk of having three or more methylated genes in sputum (Table 2).
|
Forty-four SNPs were associated with risk for promoter methylation (P < 0.15) with adjustment for covariates (Supplementary Table S4). Because of the relatively high correlation between SNPs in these genes, we tested which SNP, or set of SNPs, was most significantly associated with risk for promoter methylation by using a stepwise logistic regression model. The underlined SNPs in Supplementary Table S4 were selected from each gene (P < 0.15) to represent the allelic status for those genes. These 16 SNPs were then included with the covariates in one model and stepwise selection was used to identify the SNPs with the lowest P value. Ten SNPs from different genes were identified with four SNPs associated with increased risk for promoter methylation (OR, 1.6–4.0) and 6 SNPs with reduced risk for promoter methylation (OR, 0.4–0.7; Table 3
). Monte Carlo estimates of exact P values were calculated by permuting the case-control status for all subjects 10,000 times. The exact P value for five SNPs (rs7117042, rs5762763, rs2295146, rs7830743, and rs6998169) was <0.05 (Table 3). This result indicates that if a similar study were repeated under a null distribution (i.e., no SNPs associated with risk for promoter methylation), an association similar to that observed with any of these five SNPs would occur by chance <5% of the time. The calculated FPRP was <0.2 for four SNPs (rs7117042, rs5762763, rs2295146, and rs7830743) under the assigned prior probability range (Table 3). Findings with a FPRP of
0.2 are considered to be noteworthy.
|
|
|
Two subjects homozygous for the haplotype containing rs7117042 and four subjects, each homozygous for one of the other four common haplotypes, were selected for assessment of promoter activity (described in Materials and Methods; Supplementary Table S5). Sequencing of the 2,500-bp promoter construct revealed three haplotypes (ACGACTG, GCACTAT, and AGGCTTG), with each haplotype present in two subjects. The most distinct sequence difference was the G to C change at –590 bp. We genotyped 100 subjects selected randomly from our study population for this SNP and found that the G allele was in complete linkage disequilibrium with the T allele of the risk SNP, rs7117042, identified to be most strongly associated with high methylation index. The highest promoter activity was seen in constructs containing the ACGACTG haplotype. With this haplotype as the reference, a 23% and 38% reduction in promoter activity was seen for the GCACTAT and AGGCTTG haplotypes, respectively (Fig. 2B). These results show that the risk tag SNP is associated with a marked reduction in transcription of the MRE11A gene. MRE11A has a critical role in recognition of DSB damage. It complexes with Rad50 and Nbs1 to directly sense the DSBs, binds to the DNA, modifies the ends via 3' to 5' exonuclease activity, recruits ATM to the damaged DNA template, and dissociates the ATM dimer (32). Therefore, a reduction in level of the MRE11A protein could have a major effect on DSBRC.
| Discussion |
|---|
|
|
|---|
Another key contributing factor for aberrant de novo methylation during DNA damage is the rapid recruitment of DNMT1 to sites of DNA damage (36). Le Gac et al. (37) found that in cells treated with doxorubicin, which induces DSBs, DNMT1 is recruited by activated p53 and binds to functional Sp1 sites within promoters of the survivin, cdc2, and cdc25 genes. Moreover, the transcriptional repressor HDAC1 and the repressive chromatin mark H3K9me2 were also found at these promoters following DNA damage (37, 38). Subsequent in vitro studies showed that, following DSB DNA damage induced by doxorubicin, DNMT1 complexed with p53 was recruited to the survivin gene promoter followed by de novo methylation and gene silencing (39). Cuozzo et al. (40) provides even stronger support for a mechanistic link between DNA damage and methylation. In that study, a recombinant plasmid containing a 1-SCE1 restriction site within one copy of two inactivated tandem repeated green fluorescent protein (GFP) genes was introduced into HeLa or mouse embryonic stem cells. The restriction endonuclease 1-SCE1 was added to the cell to induce a DSB in the GFP gene at this site. Rapid gene silencing associated with homologous recombination and DNA methylation of the recombinant gene was seen and could be blocked by treatment with the demethylating agent 5-aza-deoxycytidine. Chromatin immunoprecipitation revealed that DNMT1 was bound specifically to the homologous GFP DNA. Together, these in vitro studies strongly support a direct mechanistic link between DNA damage and induction of de novo methylation by DNMT1. Our population-based studies now provide for the first time an in vivo association between DNA repair capacity and gene promoter methylation, both through a functional assay and genetic variants in genes within the DSB repair pathway. Thus, in the absence of efficient repair, the recruitment of p53, DNMT1, and transcriptional repressors to many genes, such as p16, which also contain Sp1 sites within its promoter, could lead to de novo methylation and gene silencing.
The identification of DSBRC and specific genes within this pathway as a critical determinant for gene promoter hypermethylation has important implications for basic and translational science. Our study substantiates that DNA damage that has long been recognized as an initiating event for mutagenesis is also likely a major factor in initiating aberrant promoter hypermethylation. Other DNA damage response pathways, such as apoptosis, nucleotide, and base excision repair, may also contribute to the induction of aberrant promoter hypermethylation. A major priority for our research is to replicate the provocative findings in this study along with our emerging methylation gene panel in a prospective population-based study. Genetic variants associated with promoter hypermethylation could be used to identify young smokers who would be most susceptible to induction of preneoplasia and, thus, should receive chemoprevention. In addition, the integration of these genetic variants with detection of gene promoter hypermethylation in sputum in long-term heavy smokers could lead to the first diagnostic test for incident lung cancer and affect long-term survival from this fatal disease.
| Acknowledgments |
|---|
Grant support: U01 CA097356 and the State of New Mexico as a direct appropriation from the Tobacco Settlement Fund.
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.
We thank Trisstin Maroney for her technical assistance.
| Footnotes |
|---|
Received 11/21/07. Revised 1/24/08. Accepted 2/15/08.
| References |
|---|
|
|
|---|
and β in human cancers. Cancer Res 2003;63:4620–5.
-H2AX induction. EMBO J 2006;25:3986–97.[CrossRef][Medline]This article has been cited by other articles:
![]() |
M. Tessema, Y. Y. Yu, C. A. Stidley, E. O. Machida, K. E. Schuebel, S. B. Baylin, and S. A. Belinsky Concomitant promoter methylation of multiple genes in lung adenocarcinomas from current, former and never smokers Carcinogenesis, July 1, 2009; 30(7): 1132 - 1138. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Cheng, S. Leng, H. Li, C. Huang, Y. Niu, L. Zhang, X. Liang, H. Lin, and Y. Zheng Suboptimal DNA Repair Capacity Predisposes Coke-Oven Workers To Accumulate More Chromosomal Damages in Peripheral Lymphocytes Cancer Epidemiol. Biomarkers Prev., March 1, 2009; 18(3): 987 - 993. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Verri, L. Roz, D. Conte, T. Liloglou, A. Livio, A. Vesin, A. Fabbri, F. Andriani, C. Brambilla, L. Tavecchio, et al. Fragile Histidine Triad Gene Inactivation in Lung Cancer: The European Early Lung Cancer Project Am. J. Respir. Crit. Care Med., March 1, 2009; 179(5): 396 - 401. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Leng, A. Bernauer, C. A. Stidley, M. A. Picchi, X. Sheng, M. A. Frasco, D. Van Den Berg, F. D. Gilliland, R. E. Crowell, and S. A. Belinsky Association between Common Genetic Variation in Cockayne Syndrome A and B Genes and Nucleotide Excision Repair Capacity among Smokers Cancer Epidemiol. Biomarkers Prev., August 1, 2008; 17(8): 2062 - 2069. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Leng, C. A. Stidley, A. M. Bernauer, M. A. Picchi, X. Sheng, M. A. Frasco, D. Van Den Berg, F. D. Gilliland, R. E. Crowell, and S. A. Belinsky Haplotypes of DNMT1 and DNMT3B are associated with mutagen sensitivity induced by benzo[a]pyrene diol epoxide among smokers Carcinogenesis, July 1, 2008; 29(7): 1380 - 1385. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |