|Home | About | Journals | Submit | Contact Us | Français|
The detection of tumor suppressor gene promoter methylation in sputum-derived exfoliated cells predicts early lung cancer. Here we identified genetic determinants for this epigenetic process and examined their biological effects on gene regulation. A two-stage approach involving discovery and replication was employed to assess the association between promoter hypermethylation of a 12-gene panel and common variation in 40 genes involved in carcinogen metabolism, regulation of methylation, and DNA damage response in members of the Lovelace Smokers Cohort (n=1434). Molecular validation of three identified variants was conducted using primary bronchial epithelial cells. Association of study-wide significance (P<8.2×10−5) was identified for rs1641511, rs3730859, and rs1883264 in TP53, LIG1, and BIK, respectively. These SNPs were significantly associated with altered expression of the corresponding genes in primary bronchial epithelial cells. In addition, rs3730859 in LIG1 was also moderately associated with increased risk for lung cancer among Caucasian smokers. Together, our findings suggest that genetic variation in DNA replication and apoptosis pathways impacts the propensity for gene promoter hypermethylation in the aerodigestive tract of smokers. The incorporation of genetic biomarkers for gene promoter hypermethylation with clinical and somatic markers may improve risk assessment models for lung cancer.
Lung cancer is the leading cause of cancer-related mortality in both men and women in the United States and occurs largely from chronic exposure to tobacco carcinogens (1). The development of this disease over 30 to 40 years involves field cancerization, characterized as the acquisition of genetic and epigenetic changes in oncogenes and tumor suppressor genes (TSG) throughout the airway epithelium (2–4). The silencing of TSGs through promoter methylation is now recognized as a major and causal epigenetic event that occurs during lung cancer initiation and progression (4). Genes involved in all aspects of normal cell function that include cell cycle regulation, differentiation, adhesion, and death, are silenced by promoter methylation in lung tumors (4). Based on the silencing of key TSGs in the lungs of smokers, we hypothesized that the detection of promoter methylation of TSGs in exfoliated cells in sputum would provide an assessment of the extent of field cancerization that in turn may predict early lung cancer. This hypothesis has been validated in several studies (5–7), suggesting that gene methylation in sputum could be used as a molecular marker for identifying people at high risk for cancer incidence.
The fact that promoter methylation of TSGs is a promising marker for early detection of lung cancer makes understanding factors that influence the individual susceptibility for this epigenetic process throughout the respiratory epithelium a high priority. The precise mechanism by which carcinogens disrupt the capacity of cells to maintain the epigenetic code during DNA replication and repair is largely unknown. Carcinogens within tobacco induce single- and double-strand DNA breaks which, if left unrepaired or if mis-repaired, lead to mutagenic and/or carcinogenic changes in cells (8). Accumulating evidence from our group and others suggests that extensive DNA damage could also be responsible for acquisition of promoter methylation of TSGs during lung carcinogenesis. Support for this supposition was provided through several recent studies from our group (9–12). A highly significant association was observed between DNA repair capacity (DRC) for double strand breaks measured in lymphocytes and the propensity for gene methylation detected in sputum from cancer-free smokers from the Lovelace Smokers Cohort (LSC) (9). A subsequent study in the same cohort identified the dietary factors including folate, leafy green vegetables, and multivitamin use as protective against the acquisition of gene methylation possibly through the modulation of DRC and/or the reduction of DNA damage induced by tobacco-derived carcinogens due to their anti-oxidative effect (10). Chronic exposure of human bronchial epithelial cells to non-cytotoxic levels of tobacco-derived carcinogens induced de novo methylation of TSGs, epithelial-to-mesenchymal transition, and eventually transformation (11). Cuozzo et al. (12) provides a mechanistic link between DNA damage and methylation by demonstrating activation of homologous recombination following introduction of a double strand break and subsequent methylation of the recombinant gene. Together, these studies suggest that chronic DNA damage and reduced DRC could be important determinants for inducing gene methylation.
Several sequence patterns within gene promoters that contain CpG islands and embryonic targets of polycomb-repressive complex 2 are predictive for gene predisposition for methylation in cancer, but cannot discriminate the inter-individual susceptibility for gene silencing (13–17). Sequence variants in promoters associated with reduced gene transcription lead to allele-specific methylation (ASM) and silencing in glutathione S-transferase pi (GSTP1) and O6-methylguanine-DNA methyltransferase (MGMT) in tumors and premalignant tissues (18,19). Mechanisms independent of effects on gene transcription were also identified for ASM of the reversion-induced LIM gene (20). Several studies conducting chromosome-wide or genome-wide surveys on non-imprinted, autosomal regions in human lymphocytes suggest that the majority of TSGs are not silenced by sequence variant dependent ASM (21,22). Based on the likelihood that DNA damage induced by tobacco carcinogens is an important step in the acquisition of de novo methylation and that reduced carcinogen detoxification and DRC have been associated with lung cancer (9–12,23), we tested the hypothesis that genetic variation in some genes involved in these pathways are associated with susceptibility for smokers to acquire gene-specific promoter methylation detected in sputum that contains exfoliated lung cells. A two-stage approach involving discovery and replication was employed to assess the association between promoter methylation of a 12-gene panel in members of the LSC and common variation in 40 genes involved in carcinogen metabolism, regulation of methylation, and DNA damage response, the latter including DNA damage repair, cell cycle regulation, and apoptosis. Molecular validation of significant variants was conducted using primary bronchial epithelial cell cultures.
The LSC was established in 2001 to conduct longitudinal studies on molecular markers of respiratory carcinogenesis in biological fluids such as sputum from people at risk for lung cancer (9). The enrollment initially focused on female smokers and was expanded to include male smokers in 2004. Enrollment was restricted to current and former smokers age 40 to 74 y with a minimum of 20 pack-years of smoking. Detailed information regarding sample collection was described in Supplementary Materials and Methods. All participants signed a consent form, and the Western Institutional Review Board approved this project.
Methylation of a 12-gene panel was successfully assessed in cytological adequate sputum samples from 1434 cohort members who are either Caucasian or Hispanic and for whom the genotyping call rate was ≥ 75%. Members with other ethnicities were not included in this study because of their low representation in the LSC (overall < 6%). Cohort members were split into two populations for the discovery (n=713) and replication (n=721) based on their methylation index and several non-genetic risk factors for gene methylation including gender, ethnicity, current smoking status, and age at enrollment (10). The demographic characteristics for the cohort members in the discovery and replication stages are shown in Table 1.
A detailed procedure for sputum collection and processing has been described (10). Briefly, sputum samples were stored in Saccomanno’s fixative. Sputum adequacy, defined as the presence of deep lung macrophages or Curschmann’s spiral (7), was assessed by a pathologist. Twelve genes were selected for analysis of methylation in sputum based on our previous studies establishing their association with risk for lung cancer and their specificity to methylation in epithelial cells (5–7). These twelve genes include p16, MGMT, DAPK, RASSF1A, PAX5α, PAX5β, GATA4, GATA5, SULF2, PCDH20, DAL1, and JPH3. Our ongoing study to identify the best gene panel for early lung cancer detection has screened over 40 genes and suggests that a panel of approximately 12 genes (inclusive of those described above) may provide the best sensitivity and specificity with additional genes not conferring additional risk (Leng unpublished). Genomic DNA isolated from the sputum samples was bisulfite-modified. Given the low percentage (<3%) of lung epithelial cells in sputum samples that also varied significantly between individuals, a two-stage nested MSP was used to detect methylated alleles (5–7). Our assay can reproducibly detect one methylated allele in a background of 20,000 unmethylated alleles (4). Genomic DNA isolated from cell lines that were obtained from and authenticated by the American Type Culture Collection (Manasssas, VA) with methylated or unmethylated genes from the 12-gene panel was included in each batch of the MSP assays for quality control. Experiments were conducted in cell lines passes for a maximum of 6-months post resuscitation.
Genes (n = 40) were selected based on the literature derived citations establishing the relevance of genes to the regulation of gene methylation, lung cancer, and pulmonary function (Supplementary Table 1, available online). Tag SNPs (n = 718) were selected using pairwise r2 algorithm based on the phase 2 HapMap CEU database, linkage disequilibrium (LD) database for DNA repair and carcinogen metabolism genes from the University of Southern California (24; Gilliland unpublished), or the LD database for methylation regulatory and apoptotic genes from LRRI (Leng unpublished) for Caucasian and Hispanic populations. Fifty ancestry informative markers (AIMs) were chosen to assess ancestry admixture (25). The description of the two LD databases, strategy for tag SNP selection, efficiency of selected SNPs for tagging unmeasured SNPs, and the selection of AIMs were described in detail in Supplementary Materials and Methods.
A 768-plex oligo pool assay was designed to genotype lymphocyte DNA from the LSC using the Illumina GoldenGate technology. Detailed information about genotyping and quality control was provided in Supplementary Materials and Methods.
Normal human bronchial epithelial cells (NHBECs) that were obtained by bronchoscopy from current or former smokers (n = 58) were used in this study. Expression of genes was assessed in cells at passage 1 or 2. Cells were harvested in TRI reagent (Sigma-Aldrich, St. Louis, MO) at approximately 80% confluence and mRNA was isolated following TRI reagent instructions. TaqMan real time PCR was conducted to measure the expression of candidate genes in cDNA using the ΔCT method with both PCNA and β-actin as the endogenous controls. PCNA and β-actin were selected as the endogenous controls in NHBECs because they show minimal variation in expression across the NHBEC samples and are highly correlated (Pearson correlation coefficient = 0.83, P<0.0001).
Population admixture was analyzed using the model-based, Bayesian Markov Monte Carlo algorithm implemented in STRUCTURE 2.3.3 (26). A detailed description of ascertaining population admixture in the LSC members was described in the Supplementary Materials and Methods and in reference 27.
Genetic association analysis was first conducted for all 600 SNPs that successfully passed the quality check in the 713 members in the discovery stage (Supplementary Table 2, available online). Only SNPs associated with methylation index with P values ≤ 0.10 in the discovery stage were then analyzed in the replication stage (n=721). SNPs had to pass the replication stage with P < 0.05. The SNPs that satisfied the criteria for both the discovery and replication stages also were analyzed with the combined data set. A Bonferroni adjustment over all SNPs included in the analysis (n=600) was used (P=8.2×10−5).
The methylation index based on a 12- gene panel is a discrete count variable defined as the number of genes in the panel that are methylated and ranges from 0 to 10 in the LSC (Figure 1). Poisson regression was conducted to examine the association between the independent variables, the SNPs/haplotypes, and the outcome variable, the methylation index, with adjustment for covariates selected a priori that include age, gender, ethnicity, current smoking status, and packyears. The additive model was tested for each SNP; and common homozygotes, heterozygotes, and rare homozygotes were coded as 0, 1, and 2, respectively. The Poisson regression model is a generalized linear model that can be used to estimate the mean frequency ratio (FR) and its 95% confidence interval. The FR was calculated from the regression coefficient for an explanatory variable in the Poisson regression model, β, through the equation FR=eβ. A FR value > 1 suggests that a larger value of a covariate is associated with a higher methylation index, while a FR value < 1 suggests that a larger value of a covariate is associated with a lower methylation index.
The PHASE program was used to reconstruct the haplotypes and to calculate their estimated probabilities from the htSNPs data in a block in Caucasians and Hispanics separately (28,29). The probabilities of the common haplotypes for each individual were used as explanatory variables in the Poisson regression model with adjustment for non-genetic factors to assess the association between the haplotypes and methylation index. Haplotypes with frequency <5% were combined into one group.
The genotype – expression correlation in NHBECs was examined for the SNPs associated with the methylation index. The expression data was presented by setting the level of gene expression with wild homozygote genotypes at 100%. The expression level with heterozygote and variant homozygote was calculated as the ratio relative to that of the wild homozygote. The genotype expression correlation was then analyzed using logarithm-transformed values of the relative quantification of gene expression because this transformation satisfies the normality and homoscedasticity assumption. Logarithm transformed gene expression with high MAF was compared across genotypes using the one way ANOVA. For SNPS with low MAF, a dominant model that combines variant homozygotes and heterozygotes was applied and the difference between the two genotype groups was compared using the student t test. The statistical analyses were conducted using SAS software, version 9.2 (SAS Institute, Cary, NC).
Individual SNP association analysis identified 109 SNPs associated with the methylation index with P value < 0.10 in the discovery stage (Table 2; not shown). These 109 SNPs were then analyzed in the replication stage. Eleven SNPs in LIG1, TP53, BIK, BMF, and BAK1 were associated with methylation index with P ≤ 0.05 in the replication stage (Table 2). Seven of these 11 SNPs were tag SNPs in LIG1. Combined analysis showed that the overall P value for three SNPs including rs1641511, rs3730859, and rs1883264 reached study-wide significance (Bonferroni corrected P = 8.2×10−5). All other SNPs genotyped in the replication stage failed to replicate the associations observed in the discovery stage at the significance level of 0.05. Relaxing the P value from 0.10 to 0.15 for identifying SNPs to be carried over into the replication stage did not identify any additional SNPs reaching study-wide significance.
Because multiple tag SNPs were associated with methylation index in LIG1, a haplotype-based analysis was conducted to integrate the genetic variation in each individual genomic locus. LIG1 resides in an 84-kb haplotype block (chr19: 53296285– 53381205, NCBI B36) that contains three genes: PLA2G4C, LIG1, and LOC374920 (Supplementary Figure 1A, available online). Six haplotype tagging SNPs (htSNPs) including rs11564509, rs251692, rs3730859, rs3730895, rs4802436, and rs972800 allowed the construction of the four common haplotype alleles with minimal Rh2 over 0.82 in both Caucasians and Hispanics (24). These six htSNPs accounted for > 84% of the haplotype variation in the LIG1 block in Caucasians and Hispanics in the LSC (Table 3). The likelihood ratio global tests of association indicate that the genetic variation in the haplotype block where LIG1 resides as a locus is strongly associated with methylation index in the entire LSC (P=0.00029; not shown). Hap1 strongly tagged by rs3730859 is the only risk haplotype allele associated with increased methylation index compared to all other haplotype alleles (FR=1.09, P=0.00036, Table 3). Hap3 solely tagged by rs4802436 was associated with a reduction of methylation index compared with all other haplotype alleles (FR=0.93, P=0.023, Table 3). However, when Hap2, Hap4, and others of LIG1 were used as the reference for comparison in a Poisson regression model, only Hap1 was significantly associated with increased methylation index (P=0.0032, Table 3), supporting Hap1 as the most influential haplotype allele defining the global association between this locus and methylation index.
Characterization of the LD for SNPs located 500kb surrounding each of the top three SNPs using the 1000 Genome Pilot1 CEU population identified 60 SNPs that have high LD (r2 > 0.8) with these three SNPs. Interestingly, all 60 SNPs locate within the expanded regions of the three candidate genes (Supplementary Figure1B – 1D, available online), suggesting that it is unlikely that the association seen for these three SNPs is due to long-range LD with functional SNPs located further away from these candidate regions. The functional potential for these SNPs was assessed by searching the FuncPred: Functional SNP Prediction module of SNPinfo Web Server (30). No known nonsynonymous SNPs from the genes in the four candidate regions were identified to be in high LD with these SNPs. However, 13 of the 60 SNPs were predicted either to affect the binding of transcription factors or microRNAs or to locate in an exonic splicing enhancer or silencer element (Supplementary Table 3, available online). Specifically, seven SNPs surrounding the transcription start site of LIG1 that are in high LD with rs3730859 (r2 > 0.93) were predicted to affect the binding of transcription factors.
Thus, we examined whether the top three SNPs could in cis affect the expression of LIG1, TP53, and BIK in NHBECs. The expression of LIG1, TP53, and BIK was highly correlated in NHBECs with Spearman correlation coefficients from 0.60 to 0.75 (P values <0.0001; not shown). As a control, the expression of PLA2G4C which resides in the same haplotype block and directly upstream of LIG1 had no or weak correlation with the expression of LIG1 and the other genes (Spearman correlation coefficients from 0.16 to 0.42, not shown). Genotype – expression correlation analysis found strong associations between rs3730859 and LIG1, between rs1641511 and TP53, and between rs1883264 and BIK (Table 4). The most striking difference was observed for rs3730859 with the variant homozygotes associated with a 74% reduction of LIG1 expression compared with the wild homozygotes (P=0.002). In addition, as expected, the genotype expression correlation is more significant for LIG1 using PCNA as the endogenous control compared with β-actin because LIG1 expression is cell cycle regulated (31). The genotype expression correlation for TP53 has comparable results between PCNA and β-actin, a gene whose expression during the cell cycle is influenced by the presence of DNA damage. BIK expression is not cell cycle related and as such the significant genotype expression correlation was identified using β-actin as the endogenous control. More importantly, the repressive effect of current smoking status on BIK expression (32) can only be confirmed with β-actin as the endogenous control (36% reduction in current smokers, P=0.048, Leng unpublished), supporting the use of β-actin for the relative quantification of BIK expression. No association between SNPs and PLA2G4C was identified (not shown).These results suggest that SNPs associated with methylation index could in cis affect the expression of LIG1, TP53, and BIK.
Hispanics in the LSC are highly admixed between Native American ancestry and European ancestry. The component of Native American ancestry in self-reported Hispanics ranges from 18.2 to 56.0% with an average of 33.7%. In contrast, the average component of Native American ancestry in self-reported Caucasians is < 1.0% with only 10 subjects having Native American ancestry >10%. The African ancestry is very low in both Hispanics (2.1%) and Caucasians (0.3%) enrolled in the LSC. In a comparison of methods for adjusting for ancestry, no substantial differences were observed between the results obtained from Poisson regression models that included self-reported ethnicity and those that included the components of the Native American ancestry and African ancestry (not shown). Furthermore, limiting the analysis to Caucasians with European ancestry > 90% identified similar SNP associations. Thus, population admixture is an unlikely explanation for the genetic association identified.
This is the first study to comprehensively evaluate the relationship between methylation index in a large population of lung-cancer free smokers and common variation in 40 candidate genes involved in carcinogen metabolism, regulation of methylation, and DNA damage response pathways. Our findings support a role for genetic variation in LIG1, TP53, and BIK as predictors for the acquisition of gene promoter methylation in exfoliated cells in smokers’ sputum. Specifically, molecular validation through a genotype – expression correlation analysis indicates that the three most significant SNPs in LIG1, TP53, and BIK affect the expression of these three genes in cis in NHBECs. Thus, genetic variation in genes affecting DNA replication and apoptosis impacts the propensity for gene promoter methylation in the aerodigestive tract of smokers and this effect may be driven by the change in the endogenous expression of LIG1, TP53, and BIK.
Epigenetic silencing of TSGs in lung cells is an intermediate biomarker for lung cancer. Therefore, genetic variation that is strongly associated with epigenetic silencing of TSGs should also impact the risk for lung cancer. Indeed, this premise is supported by studies that systematically assessed the association between tag SNPs in LIG1 and risk for lung cancer (33–36). The criteria for selecting these studies to be summarized include large samples size (>500 pairs cases and controls), Caucasian ethnicity, pathway or genome wide association studies, and adequate SNP coverage in LIG1. These criteria were set to minimize the confounding from the publication bias and to ensure inclusion of the sufficiently powered studies. Only association with lung cancer in moderate and heavy smokers (current or former) was presented because over 90% of cohort members enrolled in the LSC are smokers with over 20 pack-year smoking history. A case-control study conducted in the Pittsburgh Lung Screening Study cohort included 722 lung cancer patients and 929 controls who were Caucasian and current or ex-cigarette smokers with > 10 pack-year smoking history and found that subjects carrying the variant homozygote for rs10500298 that is in high LD with rs3730859 (r2=0.89) in the 1k genome project CEU population have a 96% increased risk for lung cancer (P<0.0001) (33). A second case-control study conducted in the Beta-Carotene and Retinol Efficacy Trial cohort also identified rs156640 and rs156641 that is in high LD with rs3730859 (r2=0.89) in CEU as a strong risk factor for lung cancer in Caucasian smokers with >20 pack-year of smoking history (ORs=1.20, Ps<0.003, 746 cases and 1477 controls) (34). Two lung cancer GWAS studies conducted in United States populations were queried as well and subjects with homozygous rs156641 have a 22% increased risk for lung cancer in subjects with over 20 pack-year smoking history in MD Anderson GWAS study (P<0.14, 1004 cases and 926 controls) (35). However, no association was identified between LIG1 SNPs and risk for lung cancer in two US populations used in the discovery stage in the NCI GWAS study (1652 cases and 1212 controls) (36). The SNPs in TP53 and BIK associated with methylation index were also queried in the MD Anderson and NCI GWAS with no significant association to lung cancer identified (not shown). Together, these studies support a moderate association for rs3730859 in LIG1 with risk for the smoking-induced lung cancer and this effect may be mediated by its affect on the predisposition for epigenetic silencing of TSGs.
Among the three DNA ligase families (LIG1, LIG3, and LIG4) in vertebrates, LIG1 plays an essential role in the joining of Okazaki fragments during lagging strand synthesis and is also implicated in multiple DNA repair pathways, although the role of LIG1 in DNA repair can be substituted by other DNA ligases (37–40). Deficiency in LIG1 expression results in the delayed joining of Okazaki fragments, genome instability, and increased incidence of epithelial cancers and lymphoma (37–40). Expression of LIG1 in lung tumors (n=23) is reduced by 34% compared to that in paired distant normal tissues providing further support for a functional deficiency of this gene in lung carcinogenesis (P=0.028, Leng unpublished). The strong association between genetic variation in LIG1 and methylation index found in this study provides in vivo evidence that incomplete loss of LIG1 function due to SNPs could promote the epigenetic silencing of TSGs which in turn contributes to the epithelial carcinogenesis. The mechanism for the effect of LIG1 deficiency on gene methylation is unclear. A recent study found that the replication defect in 46BR.1G1 cells which are homozygote for the Arg771Trp mutation and retain only 3 to 5% of normal LIG1 activity results in the generation of endogenous single- and double-stranded DNA breaks during S phase (40). These DNA breaks behind the replication fork induce a persistent activation of the ATM/CHK2 pathway without triggering the S-phase-specific DNA damage checkpoint and lead to an altered chromatin structure needed for efficient processing of the DNA damage throughout the cell cycle. Thus, the replication timing could be disrupted in cells with LIG1 deficiency because of the inefficient joining of Okazaki fragments in newly synthesized DNA and the extensively altered chromatin structure. This may be a mechanism by which the guidance of the cytosine-DNA methytransferase to regions that are to be methylated is disrupted during DNA replication and the normally unmethylated regions of DNA become aberrantly methylated (41). Although such a scenario is clearly speculative, an in vitro malignant transformation model or in vivo tumorigenesis model using LIG1 knock out human epithelial cells or animals could be used to directly test this hypothesis.
Evading apoptosis is a hallmark of cancer and occurs during the progression from premalignant lesions to carcinoma in the lung (42). Apoptosis occurs through two mechanisms: the death receptor FAS/FAS ligand (also called extrinsic) pathway and the mitochondrial (DNA damage–induced and p53-mediated, also called intrinsic) pathway (43). The intrinsic apoptosis pathway plays a key role in eliminating lung epithelial cells with extensive DNA damage upon carcinogen exposure (43). Thus, the association between functional SNPs in TP53 and BIK and methylation index support the concept that the TP53 mediated apoptosis pathway could prevent the acquisition of heritable genetic and epigenetic changes through eliminating lung cells with extensive DNA damage. The loss of TP53 function by gene mutation is a common genetic alteration found in human cancers and over 50% of lung cancer cases carry somatic mutation of TP53 (44). TP53 mutation or reduced level of wild-type TP53 protein was associated with apoptosis suppression during progression from premalignant lesions to carcinoma in the lung (42). Thus, rs1641511 that was associated with a 40% reduction of TP53 transcription could render lung cells more resistant to apoptosis upon extensive DNA damage in smokers, which in turn could result in an increased methylation index. Although rs1641511 is geographically closer to ATP1B2, the gene downstream of TP53, LD analysis using the Phase 2 HapMap CEU data and the 1000 genome project identify one SNP (rs72829457) around 3kb downstream of TP53 in moderate LD with rs1641511 (D'=1, r2=0.5). Moreover, the deletion/insertion polymorphisms that usually account for 15–20% of the entire genetic variation, are not available from the above two databases. Thus, it is possible that the missing deletion/insertion polymorphisms in high LD with rs1641511 are causal for the reduced expression of TP53.
BIK (Bcl2-interacting killer) is the founding member of the pro-apoptotic BH3-only proteins (45). BIK functions as a death sensor to mediate activation of the mitochondrial apoptosis pathway as shown in response to oncogenic stress signals or DNA damage in epithelial cancer cell lines (45). BIK’s role in cancer development may vary by tumor type. Inactivation of BIK in some primary tissues occurs frequently by chromosomal deletion in renal cell carcinomas (40%), colorectal cancers (22%), and gliomas (42%), but only occasionally in head and neck cancers (5%; 45). Somatic mutations that lead to amino acid changes were also detected in 8% of patients with B cell lymphomas (45). Epigenetic silencing was also proposed as a potential mechanism to silence BIK expression in specific tumor types, although this mechanism was only indirectly tested based on BIK re-expression in tumor cell lines treated with DNA demethylating agent and/or histone deacetylase inhibitor (45). Rs1883264 associated with reduced BIK expression in NHBECs is actually protective for gene methylation in smokers. This seems to be in conflict with the pro-apoptotic role of BIK identified in epithelial cancer cell lines. However, BIK is over expressed in lung tumors (n=23) versus paired distant normal tissues (P=0.015, Leng unpublished). Furthermore, CpG sites in the BIK promoter and exon1 are not methylated in lung adenocarcinomas interrogated on the Illumina HumanMethylation 27K BeadChip (Leng unpublished). Consistent with these findings in lung cancer, BIK over expression was also reported in breast tumors (46). In addition, poor prognosis of non-small cell lung cancers correlated with high expression of BIK (47). These results support different mechanisms for BIK regulation in cancer biology by tumor site and carcinogen exposure pattern. Thus, the protective effect of rs1883264 on gene methylation in smokers identified in our studies is consistent with the increased expression of BIK in lung tumors. Additional mechanistic studies are required to understand the regulation of BIK biology in cigarette smoke-induced lung cancer.
These studies identified the DNA replication and apoptosis pathways as new determinants for acquiring gene methylation in lung cells from smokers. Lung cancer GWAS studies successfully identified several loci in the human genome that are associated with risk for lung cancer at 10−8 significance level (35,36). However, these validated genetic variants only explain a small fraction of lung cancer heritability. Some of the missing heritability likely exists in the genetic variants with less significant P values. The detection of promoter methylation of TSGs in exfoliated cells in sputum as an assessment of the extent of field cancerization in the lung can be used as a functional readout to help identify the true association from the candidate SNPs with GWAS P values > 10−8. Thus, the genetic variants in LIG1 provide a proof-of-concept that SNPs strongly associated with risk for gene methylation in lung cells from smokers can be moderately associated with risk for lung cancer and should be incorporated into risk assessment models for lung cancer.
We thank Kieu C. Do, Amanda M. Bernauer, and Cynthia L. Thomas at Lovelace Respiratory Research Institute and Autumn Gaitherdavis at University of Pittsburgh School of Medicine for their technical assistance in conducting methylation-specific PCR assay, gene expression assay, and culture of human bronchial epithelial cells.
This work was primarily supported by National Cancer Institute (NCI) R01 CA097356 and the State of New Mexico as a direct appropriation from the Tobacco Settlement Fund to S.A.B. Support for previous lung cancer case control studies was through NCI U19CA148127 and R01CA121197 to C.I.A., NCI P50-CA090440 to J.M.S., individual contracts from the NCI to the University of Colorado Denver (NO1-CN-25514), Georgetown University (NO1-CN-25522), the Pacific Health Research Institute (NO1-CN-25515), the Henry Ford Health System (NO1-CN-25512), the University of Minnesota, (NO1-CN-25513), Washington University (NO1-CN-25516), the University of Pittsburgh (NO1-CN-25511), the University of Utah (NO1-CN-25524), the Marshfield Clinic Research Foundation (NO1-CN-25518), the University of Alabama at Birmingham (NO1-CN-75022), Westat, Inc. (NO1-CN-25476), and the University of California, Los Angeles (NO1-CN-25404), American Cancer Society, and the Intramural Research Program of the National Institutes of Health, NCI, Division of Cancer Epidemiology and Genetics.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org)