
Cancer Research 68, 1707, March 15, 2008. doi: 10.1158/0008-5472.CAN-07-6325
© 2008 American Association for Cancer Research
Molecular Biology, Pathobiology, and Genetics |
Promoter Methylation of Genes in and around the Candidate Lung Cancer Susceptibility Locus 6q23-25
Mathewos Tessema1,
Randy Willink1,
Kieu Do1,
Yang Y. Yu1,
Wayne Yu3,
Emi O. Machida3,
Malcolm Brock3,
Leander Van Neste4,
Christine A. Stidley2,
Stephen B. Baylin3 and
Steven A. Belinsky1
1 Lung Cancer Program, Lovelace Respiratory Research Institute; 2 Department of Internal Medicine, University of New Mexico, Albuquerque, New Mexico; 3 Cancer Biology Division, The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Baltimore, Maryland; and 4 Department of Molecular Biotechnology, Faculty of Biosciences Engineering, Ghent University, Ghent, Belgium
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
|
|---|
Chromosomal aberrations associated with lung cancer are frequently observed in the long arm of chromosome 6. A candidate susceptibility locus at 6q23-25 for lung cancer was recently identified; however, no tumor suppressor genes inactivated by mutation have been identified in this locus. Genetic, epigenetic, gene expression, and in silico screening approaches were used to select 43 genes located in 6q12-27 for characterization of methylation status. Twelve (28%) genes were methylated in at least one lung cancer cell line, and methylation of 8 genes was specific to lung cancer cell lines. Five of the 8 genes with the highest prevalence for methylation in cell lines (TCF21, SYNE1, AKAP12, IL20RA, and ACAT2) were examined in primary lung adenocarcinoma samples from smokers (n = 100) and never smokers (n = 75). The prevalence for methylation of these genes was 81%, 50%, 39%, 26%, and 14%, respectively, and did not differ by smoking status or age at diagnosis. Transcription of SYNE1, AKAP12, and IL20RA was completely silenced by hypermethylation and could be restored after treatment with 5-aza-2-deoxycytidine. Significant associations were found between methylation of SYNE1 and TCF21, SYNE1 and AKAP12, and AKAP12 and IL20RA, indicating a coordinated inactivation of these genes in tumors. A higher prevalence for methylation of these genes was not associated with early-onset lung cancer cases, most likely precluding their involvement in familial susceptibility to this disease. Together, our results indicate that frequent inactivation of multiple candidate tumor suppressor genes within chromosome 6q likely contributes to development of sporadic lung cancer. [Cancer Res 2008;68(6):1707–14]
 |
Introduction
|
|---|
Lung cancer remains the leading cause of cancer death among both men and women, with an estimated 160,440 deaths and 173,000 new cases annually in the United States alone (1). It kills more people than breast, colon, and prostate cancer combined. Lung tumors develop through mutation or epigenetic silencing via promoter hypermethylation of numerous genes involved in cell cycle control, growth regulation, signaling, apoptosis, and adhesion (2). Some of the earliest changes detected in the bronchial epithelium of smokers include loss of heterozygosity (LOH) at chromosomes 3p21, 9p21, and 17p13 (3). Within these areas of allelic loss, inactivation of the remaining allele by promoter hypermethylation of RASSF1A and p16 genes and by mutation of the p53 gene is commonly seen in non–small cell lung cancers (NSCLC) (4). Methylation of the p16 gene is one of the earliest changes in lung cancer development, occurring in the field of epithelial damage induced by carcinogens within tobacco and increasing in prevalence during histologic progression of adenocarcinoma and squamous cell carcinoma (5–7).
Another chromosome region where LOH is frequently observed in adenocarcinomas from both smokers and never smokers is 6q. LOH within 6q22-27 occurs in 30% to 55% of tumors (8), whereas deletions between 6q14 and 6q24 have been detected in 60% of primary lung tumors (9). Moreover, this locus has been associated with lung cancer susceptibility. A recently conducted genome-wide linkage analysis of 52 extended pedigrees with a minimum of three family members with aerodigestive cancer identified a lung cancer susceptibility locus at chromosome 6q23-25 (10). Smoking at any level also was found to increase risk in carriers of the inherited susceptibility locus on chromosome 6q. The high frequency for LOH and the potential existence of a susceptibility locus within 6q23-25 supports the existence of tumor suppressor genes inactivated through the classic Knudson's two-hit model in which complete loss of gene function arises through loss of one allele and mutation of the second allele (11). However, only one candidate tumor suppressor gene, p34, localized to 6q25 has been identified, but this gene was not found to be associated with familial lung cancer susceptibility (12). This scenario is reminiscent of the chromosome 3p14-25 in which LOH is commonly seen in lung tumors, although no major tumor suppressor gene inactivated by mutation has been identified in this locus. Rather, genes inactivated by promoter hypermethylation at prevalences ranging from 30% to 58% have been identified within this locus. These include RASSF1A, BLU, SEMA38, and retinoic acid receptor β (4).
Chromosome 6 is one of the gene-rich and CpG island–rich chromosomes that contains
1,557 genes and 1,070 CpG islands (13, 14). The fact that chromosome 6q accumulates genetic aberration in the form of LOH could also make this region a hotspot for silencing of genes by promoter hypermethylation. Support for this supposition is growing as genes silenced by methylation in lung tumors are now being identified within 6q. Estrogen receptor
, which maps to 6q25, has been shown by our laboratory to be silenced by promoter hypermethylation in 20% and 36% of lung tumors from smokers and never smokers, respectively (15). Recently, Smith et al. (16) identified a tumor suppressor gene, TCF21, within the 6q23-24 locus that is normally expressed in lung airway epithelial cells but silenced in aerodigestive tumors. The purpose of this investigation was to identify novel genes in and around the candidate lung cancer susceptibility locus 6q23-25 that are inactivated by promoter hypermethylation and to compare their prevalence in adenocarcinomas from smokers and never smokers. A multifaceted strategy was undertaken that used comparative genomic hybridization (CGH), a transcriptome microarray, in silico screening, and a candidate approach to select and evaluate the methylation patterns of 43 genes.
 |
Materials and Methods
|
|---|
Cell lines, tumors, and control samples. Twelve lung cancer cell lines (H23, H358, H1568, H1993, H2023, H2085, H2228, H2009, Calu-3, Calu-6, SKLU-1, and SW-900) and 4 breast cancer cell lines (MCF-7, T47D, MDA-MB-231, and MDA-MB-435) obtained from the American Type Culture Collection (ATCC) and maintained in ATCC-recommended medium were used in this study. A total of 175 frozen lung adenocarcinomas (100 from smokers and 75 from never smokers) were obtained from frozen tumor banks at Johns Hopkins and the Mayo Clinic (Supplementary Table S1). Each tumor was histologically classified as adenocarcinoma according to the current WHO criteria and considered suitable for study based on the presence of >80% tumor cells (17). Peripheral blood mononuclear cells (PBMC) isolated from never smokers using Histopaque-1077 (n = 20) and human bronchial epithelial cells (HBEC) isolated from cancer-free smokers who underwent bronchoscopy (n = 20) were used as controls.
CGH analysis. Regional copy number alterations were examined using array CGH analysis as described (18). Briefly, DNA from four lung cancer cell lines (H1435, H1568, H2009, and Calu-6) was isolated using the Puregene DNA isolation kit (Gentra Systems). Genomic DNA was fragmented, random prime labeled, and hybridized to an oligonucleotide microarray that contains 22,500 elements designed for expression profiling (Human 1A V2, Agilent Technologies), for which 16,097 unique map positions were defined (Build 35). The median interval between mapped elements is 54.8 kb; 96.7% of intervals are <1 Mb and 99.5% are <3 Mb. Fluorescence ratios of scanned images of the arrays were calculated as the average of two paired arrays (dye swap), and the raw array CGH profiles were processed to identify statistically significant transitions in copy number using a segmentation algorithm as described (18). Significant copy number changes are determined based on segmented profiles only.
Treatment and genome-wide transcriptome array. Six lung adenocarcinoma cell lines, three (H23, H1568, and H1993) originated from smokers and three (H2023, H2085, and H2228) from never smokers, were used for the transcriptome array. Cells at log phase of growth were grouped in duplicate into control, trichostatin A (TSA) alone, and 5-aza-2-deoxycytidine (DAC) alone treatment groups and treated as follows: control (culture medium), TSA [culture medium containing 300 nmol/L TSA (stock solution 5 µmol/L in ethanol; Sigma) for 18 h], or DAC [culture medium containing 5 µmol/L DAC (stock solution 10 mmol/L in PBS; Sigma) for 96 h with fresh medium containing the drug changed every 24 h]. Cells in the control and TSA groups underwent fresh medium changes in parallel with the DAC group, and TSA treatment was done 18 h before all treatment groups were harvested in TRI reagent (Sigma). Genome-wide transcriptome array analysis was conducted essentially as described (19). Briefly, total RNA was isolated and quantified, and its quality was assessed with the 2100 Bioanalyzer (Agilent Technologies). Sample amplification and labeling procedures were carried out using the Low RNA Input Fluorescent Linear Amplification kit (Agilent Technologies), and the labeled cRNA was purified and quantified. RNA spike-in controls (Agilent Technologies) were added to RNA samples before amplification. Samples (0.75 µg) labeled with Cy3 or Cy5 were mixed with control targets (Agilent Technologies), assembled on the oligo microarray, hybridized, and processed according to the Agilent microarray protocol. Scanning was performed with the Agilent G2565BA microarray scanner using settings recommended by Agilent Technologies. All calculations were performed using the R statistical computing platform and packages from the Bioconductor bioinformatics software project (20–22).
DNA extraction and modification. High molecular weight DNA was isolated from cell lines, PBMCs, HBECs, and tumor samples using a standard phenol-chloroform extraction method. Genomic DNA was modified using the EZ DNA Methylation-Gold kit (Zymo Research) as described by the manufacturer, and 40 ng of modified DNA were used per PCR for combined bisulfite modification and restriction analysis (COBRA) or methylation-specific PCR (MSP).
Combined bisulfite modification and restriction analysis. Forty-three CpG islands located in gene promoter regions within 6q were selected using genetic (array CGH), epigenetic (gene expression array after demethylation), and in silico screening (presence of a promoter CpG island and predicted candidate tumor suppressor function). The selection criteria used for each gene are shown in Table 1
. The DNA methylation profile of these 43 genes was analyzed in 12 lung cancer cell lines using COBRA. The COBRA primers were designed to amplify gene promoter regions where multiple methylation-sensitive restriction enzyme recognition sites were present. Primer sequences, amplification conditions, and the number of restriction sites within the amplified region of each gene are described in Table 1 and Supplementary Table S2. PCRs were conducted in a total volume of 50 µL, containing 1x PCR buffer with 2.5 mmol/L MgCl2, 62.5 µmol/L of each deoxynucleotide triphosphate, 0.25 pmol each of the forward and reverse primers, and 2.5 units of Taq Gold polymerase (Perkin-Elmer). The reaction mixes were first heat denatured for 10 min at 94°C followed by 40 amplification cycles and 5 min of extension at 72°C. Each PCR cycle consisted of equal time of denaturation at 94°C, annealing (see Supplementary Table S2), and extension at 72°C. The resulting PCR products were digested with the specific enzyme (Supplementary Table S2), mixed with 5 µL of loading buffer, resolved by electrophoresis through 3% agarose gel, stained with ethidium bromide, and visualized under UV illumination.
MSP and bisulfite sequencing. Promoter regions that showed tumor-specific methylation in the lung cancer cell lines (i.e., no methylation in normal HBECs and PBMCs) from COBRA were further analyzed in primary lung adenocarcinomas from smokers and never smokers using MSP (for primers and amplification conditions, see Supplementary Table S3). For genes with a high prevalence of methylation (>20%), the density of methylation at each CpG within the COBRA-assayed region was further characterized by bisulfite sequencing. Bisulfite-modified DNA from a HBEC line and a primary adenocarcinoma, and lung cancer cell lines was amplified (using the COBRA primers that do not discriminate between methylated and unmethylated DNA) and cloned into pCR II cloning vector (Invitrogen), and five clones were sequenced per sample. Within the promoter regions of SYNE1, AKAP12, and IL20RA, 41, 33, and 30 CpGs, respectively, were analyzed by bisulfite sequencing, and the density of promoter methylation of each gene was described as percent of methylated CpGs from the total number of CpGs analyzed in five clones.
Gene expression analysis. Sham-, TSA-, and DAC-treated cell lines were harvested in TRI reagent, and RNA was isolated following TRI reagent instructions (Molecular Research Center, Inc.). Total RNA (3 µg) was reverse transcribed using the SuperScript First-Strand Synthesis System for reverse transcription-PCR (RT-PCR; Invitrogen) according to the cDNA synthesis protocol from Invitrogen. Expression of SYNE1, AKAP12, and IL20RA was evaluated using RT-PCR and electrophoresis in 3% agarose gels. In total, five pairs of RT-PCR primers were used. SYNE1-F1/R1 was used to amplify the longest (L1) and longer (L2) SYNE1 transcripts and SYNE1-F2/R2 to amplify all four SYNE1 transcripts: L1, L2,
(A), and β (B). AKAP12-F1/R1 and AKAP12-F2/R2 amplified the A and B transcripts of AKAP12, respectively, whereas IL20RA-F1/R1 amplified the IL20RA transcript. Large introns were included in the RT-PCR amplification regions and isolated RNA was treated with DNase to exclude PCR products from contaminating DNA. RT-PCR primers and amplification conditions are described in Supplementary Table S4.
Data analysis. Patient characteristics, including age, gender, smoking status, and race, and gene methylation are summarized with mean and SD for continuous variables and proportions for categorical variables. Fisher's exact test was used to assess the association between methylation and patient characteristics and the association between pairs of genes. Methylation of genes and patient characteristics were assessed for association with survival. Survival time was calculated from time of diagnosis until death or last follow-up visit. Individual genes and variables were assessed with Kaplan-Meier plots and the log-rank test. Proportional hazards models were used to simultaneously examine multiple genes while adjusting for patient characteristics. All analyses were conducted using Statistical Analysis System version 9.1.3.
 |
Results
|
|---|
Gross copy number alteration of chromosome 6 in lung cancer. Four lung cancer cell lines (H1435, H1568, H2009, and Calu-6) were examined by array CGH analysis to identify specific loci with copy number alterations. The X chromosome of each cell line was evaluated to verify the quality of array performance and for data normalization (Fig. 1A
). The array CGH profile of H1435 (adenocarcinoma cell line from a never smoker) revealed a copy number alteration of an
70 Mbp segment in the region between 6p12.1 and 6q22.3 (Fig. 1B). Similarly, a single deletion of
15.41 Mbp segment in the region between 6p22.1 and 6p12.2 and 6q25-27 was observed in the Calu-6 cell line (data not shown). Chromosome 6 of the other two cell lines, H1568 and H2009, showed no significant copy number alteration (data not shown). Genes within 6q12-22.3 and 6q25-27 that showed significant copy number alteration and harbor CpG islands in their promoter regions were selected for assessment of methylation (Table 1).

View larger version (15K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 1. Array CGH analysis of H1435, a lung adenocarcinoma cell line from a never smoker. A, results from the X chromosome were used to verify that the quality of array performance and data normalization were good (Window = 2.0 Mb; Z-score = 1.5). B, a large single-copy deletion between 6p12.1 and 6q22.33 was found within chromosome 6 (Window = 2.0 Mb; Z-score = 2.5).
|
|
Promoter hypermethylation of candidate genes within chromosome 6q. A total of 43 genes with CpG-rich promoters that met the CpG island criteria of Takai and Jones (23) were selected from 6q for methylation analysis (Table 1). Thirteen genes were identified through array CGH analysis of H1435 and Calu-6 that localized to the region of LOH in 6q. Twenty-one genes were obtained from the transcriptome array, and 4 of these were detected by both assays. The remaining 13 genes were selected in silico based on biological function. The COBRA analysis revealed that 31 genes (72.1%) were devoid of promoter methylation (0 of 12) in lung cancer cell lines (Fig. 2
). Two genes, PNLDC1 and PLAGL1, that showed a high prevalence for methylation in the lung cancer cell lines (100% and 83.3%, respectively) also were methylated in 100% of normal HBECs and PBMCs. PNLDC1 was completely methylated (no unmethylated signal) in all cell lines and control samples (Fig. 3A
). In contrast, PLAGL1 was partially methylated (approximately 1:1 methylated and unmethylated alleles) in 75% (9 of 12) of the lung cancer cell lines and all normal HBECs and PBMCs (Fig. 3B). The PLAGL1 promoter CpG island analyzed in this study is located within the differentially methylated region (DMR) observed in various human tissue (24). However, 25% (3 of 12) of lung cancer cell lines showed loss of differential (allele specific) methylation of the DMR (Fig. 3C). Two of the cell lines (SKLU-1 and SW-900) were completely unmethylated, whereas the third cell line (H358) was completely methylated (biallelic methylation; Fig. 3C). Two genes, NR2E1 and DACT2, were methylated in 67% and 83% of the lung cancer cell lines, respectively, and 100% of the PBMCs but not in normal HBEC (Fig. 2; Table 2
). The remaining eight genes (TCF21, IL20RA, SYNE1, AKAP12, ACAT2, LATS1, AIM1, and SASH1) showed tumor-specific methylation [i.e., no methylation in normal HBECs (0 of 20) and PBMCs (0 of 20); Fig. 2; Table 2]. Five of these eight genes were methylated at prevalence of >20% in cell lines and thereby selected for further analysis in primary adenocarcinomas.

View larger version (24K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 2. DNA methylation profile of selected chromosome 6 genes in lung cancer cell lines. Promoter regions of 43 genes located within the long arm of chromosome 6 (6q) were selected based on the presence of a major copy number alteration, significant increase in gene expression after DAC treatment, and presence of CpG island and candidate tumor suppressor function. DNA methylation profile of these CpG islands was screened in 12 lung cancer cell lines using COBRA, and the percentage of cell lines with methylation of a specific promoter region is shown for each of the 43 CpG islands.
|
|

View larger version (67K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 3. Hypermethylation of PNLDC1 promoter and PLAGL1 DMR in normal and malignant cells. A, COBRA using the BstUI restriction enzyme that recognizes only methylated CGCG sites revealed complete methylation (no unmethylated DNA) of PNLDC1 promoter in lung cancer cell lines as well as normal HBECs and PBMCs. B, similar analysis of the PLAGL1 DMR with BstUI showed partial digestion [a roughly 1:1 mixture of methylated (cut) and unmethylated (uncut) DNA] in all normal HBECs and PBMCs as well as 75% of the lung cancer cell lines (two examples of each sample type are shown). C, three lung cancer cell lines (25%) have lost the normal differential (partial) methylation. Whereas SKLU-1 and SW-900 were completely unmethylated, H358 was completely methylated. D, loss of differential methylation (complete methylation) of the PLAGL1 DMR was seen in three of four (75%) breast cancer cell lines tested.
|
|
View this table:
[in this window]
[in a new window]
|
Table 2. Prevalence of promoter hypermethylation of chromosome 6q genes in primary lung adenocarcinoma samples from smokers and never smokers
|
|
Hypermethylation of 6q genes in primary lung adenocarcinoma. Primary lung adenocarcinomas from 100 smokers (37 current, 59 former, and 4 unclassified) and 75 never smokers were screened for methylation. The prevalence for promoter hypermethylation of TCF21, SYNE1, AKAP12, IL20RA, and ACAT2 was 81%, 50%, 39%, 26%, and 14%, respectively (Table 2). None of the genes showed significant difference in prevalence of methylation between tumors from smokers versus never smokers and did not vary by age of cancer diagnosis. The methylation status of each CpG located in the promoter regions of SYNE1, AKAP12, and IL20RA was characterized further by bisulfite sequencing. The TCF21 promoter was not sequenced because dense methylation of this gene has already been described in NSCLC and head and neck tumors (16). SYNE1 and AKAP12 gene promoters showed virtually no methylation (0 of 205 and 1 of 165 CpGs, respectively) in normal HBECs, whereas IL20RA showed minimal methylation (12 of 150 CpGs; Fig. 4A
). In contrast, the promoter CpG islands of SYNE1, AKAP12, and IL20RA showed dense methylation in the primary adenocarcinomas (80%, 68%, and 43% of the CpGs) and lung cancer cell lines (88%, 95%, and 89% of the CpGs), respectively (Fig. 4A). The slightly lower density of methylation in the primary tumor for SYNE1 and AKAP12 compared with the cancer cell lines is likely due to the presence of normal inflammatory and stromal cells, whereas some heterogeneity for methylation is likely present in the promoter regions of IL20RA in the sequenced tumor specimen.

View larger version (32K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 4. Dense methylation of SYNE1, AKAP12, and IL20RA promoters silences expression of these genes in lung adenocarcinoma. A, a total of 41, 33, and 30 CpG sites within the promoter CpG islands of SYNE1, AKAP12, and IL20RA, respectively, were amplified with the COBRA primers that do not discriminate between methylated and unmethylated DNA, and five clones were sequenced per sample. As shown by the key, each circle is divided into five equal parts, representing the methylation status of one clone for the specified CpG. Filled regions, methylated clones; open regions, unmethylated clones. B, in the Calu-6 cell line where the promoter region of SYNE1 was methylated, expression of the longer and longest SYNE1 transcripts (SYNE1 L) as well as all four SYNE1 transcripts (SYNE1 all) was silenced. Expression could be restored after DAC treatment. Similarly, AKAP12-A and AKAP12-B transcripts (C) as well as the IL20RA transcript (D) were silenced in Calu-6 cell lines where the promoter regions of the two genes were methylated. Expression of these genes in Calu-6 was restored after DAC treatment. In contrast, robust expression of each gene was observed in the H2228 cell line where the promoter regions of all three genes were devoid of methylation, and this expression pattern in H2228 was not altered after DAC treatment.
|
|
Promoter methylation is associated with silencing of gene expression. The relationship between promoter hypermethylation and gene expression of SYNE1, AKAP12, and IL20RA was further characterized by determining the levels of transcription of these genes in cell lines. For SYNE1 and AKAP12, two different primer pairs (SYNE1-F1/R1 and SYNE1-F2/R2, and AKAP12-F1/R1 and AKAP12-F2/R2) were used to measure the expression of the different transcripts. Robust expression of the different transcripts of the three genes was seen in H2228, consistent with the lack of methylation seen for these genes in this cell line (Fig. 4B–D). In contrast, cell lines such as Calu-6, where the promoter regions of the three genes were hypermethylated, all four transcripts of SYNE1, both A and B transcripts of AKAP12, as well as the IL20RA transcript were not expressed (Fig. 4B–D).
DAC treatment restores expression of silenced genes. Cell lines with methylated and unmethylated promoters for each of the three genes were treated with DAC to verify the cause and effect relationship between promoter CpG island hypermethylation and loss of gene expression. In Calu-6 cells, the expression of SYNE1-
, SYNE1-β, AKAP12-A, AKAP12-B, and IL20RA transcripts was completely silenced by promoter hypermethylation but was restored to normal levels (levels seen in cell lines with unmethylated promoter) after DAC treatment (Fig. 4B–D). The longer and longest transcripts of SYNE1 (SYNE1 L) were only partially reexpressed after DAC treatment (Fig. 4B). In contrast, the expression of these genes in the H2228 cell line where the promoter CpG islands were unmethylated was not altered by the DAC treatment (Fig. 4B–D). Treatment with TSA did not affect the expression of these genes.
Survival and the association of methylation between genes. As expected, patients with advanced-stage lung adenocarcinoma (stages II–IV) had significantly shorter survival than those with early-stage (stage I) disease (P = 0.001; hazard ratio, 2.24). No significant association was seen between survival and promoter hypermethylation for any of the five genes individually or in combination. The association of methylation of the five genes was assessed in tumors independent of smoking status. A positive association was found in early-stage adenocarcinomas between methylation of SYNE1 and TCF21 (P = 0.02), SYNE1 and AKAP12 (P < 0.001), and AKAP12 and IL20RA (P = 0.045). Strikingly, 88% of tumors methylated for SYNE1 were methylated for TCF21, whereas 74% of tumors positive for methylation of AKAP12 were also methylated for SYNE1. Furthermore, the methylation status of SYNE1 was associated with the methylation of three or more of the other four genes (P < 0.001). A similar scenario was also seen for AKAP12 (P < 0.002).
 |
Discussion
|
|---|
Chromosomal regions with recurrent LOH in malignancies, such as the long arm of chromosome 6 in lung cancer, have been predicted to contain one or more tumor suppressor genes. Genetic, epigenetic, and gene expression assays were used in this study to identify frequent promoter hypermethylation–mediated inactivation of five genes (TCF21, SYNE1, AKAP12, IL20RA, and ACAT2) in primary lung adenocarcinomas from smokers and never smokers. All five genes are candidate tumor suppressor genes located in and around the candidate 6q23-25 lung cancer susceptibility locus and, with the exception of TCF21, have not been shown to be epigenetically silenced in lung cancer. Promoter hypermethylation of TCF21 has been described in a limited number of aerodigestive tumors, primarily with squamous cell histology (16). Overall, 12 of the 43 promoter CpG islands screened (28%) showed methylation in one or more lung cancer cell lines. Four (9%) were also methylated in normal HBEC and PBMC, corroborating the recently conducted comprehensive CpG island methylation screening (256 CpG islands were analyzed from chromosome 6 alone) where
12% of the promoter CpG islands were methylated in normal human tissue (13).
The genes selectively methylated in lung adenocarcinomas have diverse functions that when abrogated could contribute to neoplastic transformation. TCF21 is a basic helix-loop-helix transcription factor that is critical for lung development. TCF21-null mice are born alive but die shortly after birth due to severely underdeveloped lungs and kidneys that lack alveoli and mature glomeruli, indicating a critical role of TCF21 in lung development (25). Promoter hypermethylation of TCF21 silences gene expression, and its reexpression through retroviral transfection as well as DAC treatment reverses neoplastic properties (16). The high prevalence (81%) for inactivation of this gene in primary lung adenocarcinomas shown in this study was also seen in squamous cell carcinomas (16), making it one of the most common targets in NSCLC.
The SYNE1 gene spans 0.5 Mb of genomic DNA, comprising 147 exons and encoding four different transcripts through alternative splicing (26). It is a multifunctional gene involved in cytokinesis, nuclear organization and the structural integrity, and function of the Golgi apparatus (27, 28). Hypermethylation in the promoter CpG island and mutation in the coding region of SYNE1 were recently reported in approximately 95% and 10% of colorectal cancer cases, respectively, indicating a predominant role of epigenetic regulation of SYNE1 in this cancer (19, 29). Hypermethylation of SYNE1 associated with the silencing of all four transcripts was seen in 50% of lung adenocarcinoma samples from smokers and never smokers. The size of this gene prohibited sequencing of unmethylated cell lines and tumors; however, the low prevalence for mutation in colon tumors suggests that methylation is the dominant mechanism for altering function of this gene.
AKAP12/Gravin is one of the A-kinase anchoring proteins (AKAP) that regulates mitogenesis by anchoring key signaling proteins (such as protein kinase A and protein kinase C) and modulating the expression of genes involved in cell cycle and apoptosis (30). It suppresses tumor cell viability and growth by inducing apoptosis via caspase-3, up-regulation of Bax, and down-regulation of Bcl-2 expression (30). AKAP12 also induces cell cycle arrest by increasing the expression of cell cycle checkpoint proteins Cip1/p21 and Kip1/p27 and reducing the cyclin D1 expression (30). AKAP12 also is implicated in decreasing vascular endothelial growth factor expression and vascular growth that is critical for tumor growth and metastasis (31). Overexpression of SSeCKS (rodent orthologue of human AKAP12) in rats reduces metastasis of prostate cancer cells to the lungs (31). The tumor suppressor function of AKAP12 is strongly suppressed in various cancers, including gastric, prostate, ovarian, and breast cancer, and promoter hypermethylation–mediated inactivation of AKAP12 has been shown in gastric carcinoma and childhood myeloid malignancies (32, 33). These studies, coupled with our findings of methylation of AKAP12 in primary lung adenocarcinoma, support the premise that this gene represents another common target for silencing in cancer.
IL20RA encodes a receptor for interleukin (IL)-20 and IL-24 (34). IL-24 (MDA-7) induces G2-M cell cycle arrest and apoptotic cell death through up-regulation of proapoptotic proteins (Bax and Bak) and down-regulation of antiapoptotic proteins (Bcl-2 and BCL-xL; refs. 35, 36). IL-24 gene transfer kills various cancer cells, including lung and breast cancer, by inhibition of Wnt/phosphatidylinositol 3-kinase signaling pathways (37). This cytotoxicity of IL-24 could be abolished by anti–IL-24 as well as anti-IL20R antibodies, indicating that the bystander tumor-specific cytotoxicity of IL-24 is mediated via IL-24 binding to IL-20R1/IL-20R2 heterodimeric receptor complex (38). Therefore, silencing of IL20RA by promoter hypermethylation represents another survival strategy of lung cancer cells by blocking cell cycle arrest and apoptotic signals from IL-24. Although IL20RA methylation was observed in a portion (26%) of lung adenocarcinomas, it was significantly associated with methylation of AKAP12. Interestingly, methylation of the SYNE1 promoter was also significantly associated with methylation of TCF21 and AKAP12 genes. Taken together, these findings indicate a coordinated silencing of multiple growth-regulatory genes within chromosome 6q in lung cancer.
Our screen of 6q for methylation also identified a partial (allele specific) methylation of PLAGL1 in normal HBECs and PBMCs, indicating epigenetic-mediated regulation of this gene in adult human lung and peripheral blood cells. PLAGL1 is a maternally imprinted zinc finger protein with functional similarities to p53. It interacts physically and functionally with cyclic AMP-responsive element binding protein–binding protein and p300 and induces cell cycle arrest and apoptosis (39–41). PLAGL1 is differentially methylated and imprinted in normal human peripheral blood lymphocytes, skin fibroblasts, ovary, placenta, and some human fetal tissues including fetal lung (24, 42–44). Down-regulation of PLAGL1 has been reported in lung, breast, and ovarian cancers, and in ovarian cancers, it correlates with hypermethylation of the DMR (43, 45). In this study, 25% of lung cancer cell lines showed loss of differential methylation that could be indicative of loss of imprinting. However, PCR-based assays could not be used to assess PLAGL1 methylation in primary lung tumors due to the presence of lymphocytic infiltration that would contain normal monoallelic methylation. Interestingly, three of four breast cancer cell lines showed complete methylation, suggesting that abrogation of function of this gene could be a contributing factor in breast cancer development (Fig. 3D). In addition to PLAGL1 (6q24), the region between 6q24 and 6q27 also harbors several imprinted genes that include HYMAI (6q24.1-2), IGF2R (6q25.3), SLC22A2 (6q26), and SLC22A3 (6q26-27; reviewed in ref. 46). Normal allele-specific expression of imprinted genes that are located in close proximity as well as their deregulation in cancer is often coordinately regulated, making this an interesting locus for future investigation.
This study also identified cell- and tumor-specific methylation of the NR2E1 and DACT2 promoters. Methylation of these genes was seen in all normal PBMCs and the majority of lung tumor–derived cell lines that are free of PBMC but not in HBECs. In contrast, PNLDC1 was methylated in all normal PBMCs and HBECs evaluated. The function of these genes is not well described; however, the differential methylation patterns seen in the cell types suggest that promoter methylation–mediated regulation may be occurring at the embryonic stage and/or in response to differentiation.
A final key question from this study is the contribution of the identified genes silenced in adenocarcinoma to familial susceptibility for lung cancer. It is clear for early-onset cancers (such as breast cancers associated with germ-line BRCA1/2 mutations) that the second hit in the affected locus plays a major role in the genesis of the disease. However, the three methylated genes (IL20RA, SYNE1, and AKAP12) within the 6q23-25 locus where LOH is frequent in susceptible families did not show an increased prevalence for silencing in early-onset (<50 years) lung cancer patients. Thus, although the genes discovered in this study probably do not predispose one for lung cancer, their high prevalence for silencing likely contributes to the development of sporadic lung cancer in smokers and never smokers.
 |
Acknowledgments
|
|---|
Conflict of interest: S.A. Belinsky is a consultant to Oncomethylome Sciences. Under a licensing agreement between Lovelace Respiratory Research Institute and Oncomethylome Sciences, nested MSP was licensed to Oncomethylome Sciences, and the author is entitled to a share of the royalties received by the Institute from sales of the licensed technology. The Institute, in accordance with its conflict-of-interest policies, is managing the terms of these arrangements. The commercial rights to standard MSP also belong to Oncomethylome Sciences. S.B. Baylin is a consultant to Oncomethylome Sciences and is entitled to royalties from any commercial use of this procedure.
Grant support: NIH grant R01 ES008801 (S.A. Belinsky).
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.
 |
Footnotes
|
|---|
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).
Received 11/19/07.
Revised 1/ 9/08.
Accepted 1/17/08.
 |
References
|
|---|
- Jemal A, Siegel R, Ward E, Murray T, Xu J, Thun MJ. Cancer statistics, 2007. CA Cancer J Clin 2007;57:43–66.[Abstract/Free Full Text]
- Minna JD, Fong K, Zochbauer-Muller S, Gazdar AF. Molecular pathogenesis of lung cancer and potential translational applications. Cancer J 2002;8 Suppl 1:S41–6.[Medline]
- Wistuba II, Lam S, Behrens C, et al. Molecular damage in the bronchial epithelium of current and former smokers. J Natl Cancer Inst 1997;89:1366–73.[Abstract/Free Full Text]
- Ito M, Ito G, Kondo M, et al. Frequent inactivation of RASSF1A, BLU, and SEMA3B on 3p21.3 by promoter hypermethylation and allele loss in non-small cell lung cancer. Cancer Lett 2005;225:131–9.[CrossRef][Medline]
- Belinsky SA. Gene-promoter hypermethylation as a biomarker in lung cancer. Nat Rev Cancer 2004;4:707–17.[CrossRef][Medline]
- Belinsky SA, Nikula KJ, Palmisano WA, et al. Aberrant methylation of p16(INK4a) is an early event in lung cancer and a potential biomarker for early diagnosis. Proc Natl Acad Sci U S A 1998;95:11891–6.[Abstract/Free Full Text]
- Kim DH, Nelson HH, Wiencke JK, et al. p16(INK4a) and histology-specific methylation of CpG islands by exposure to tobacco smoke in non-small cell lung cancer. Cancer Res 2001;61:3419–24.[Abstract/Free Full Text]
- Girard L, Zochbauer-Muller S, Virmani AK, Gazdar AF, Minna JD. Genome-wide allelotyping of lung cancer identifies new regions of allelic loss, differences between small cell lung cancer and non-small cell lung cancer, and loci clustering. Cancer Res 2000;60:4894–906.[Abstract/Free Full Text]
- Goeze A, Schluns K, Wolf G, Thasler Z, Petersen S, Petersen I. Chromosomal imbalances of primary and metastatic lung adenocarcinomas. J Pathol 2002;196:8–16.[CrossRef][Medline]
- Bailey-Wilson JE, Amos CI, Pinney SM, et al. A major lung cancer susceptibility locus maps to chromosome 6q23-25. Am J Hum Genet 2004;75:460–74.[CrossRef][Medline]
- Knudson AG. Two genetic hits (more or less) to cancer. Nat Rev Cancer 2001;1:157–62.[CrossRef][Medline]
- Wang M, Vikis HG, Wang Y, et al. Identification of a novel tumor suppressor gene p34 on human chromosome 6q25.1. Cancer Res 2007;67:93–9.[Abstract/Free Full Text]
- Eckhardt F, Lewin J, Cortese R, et al. DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet 2006;38:1378–85.[CrossRef][Medline]
- Mungall AJ, Palmer SA, Sims SK, et al. The DNA sequence and analysis of human chromosome 6. Nature 2003;425:805–11.[CrossRef][Medline]
- Issa JP, Baylin SB, Belinsky SA. Methylation of the estrogen receptor CpG island in lung tumors is related to the specific type of carcinogen exposure. Cancer Res 1996;56:3655–8.[Abstract/Free Full Text]
- Smith LT, Lin M, Brena RM, et al. Epigenetic regulation of the tumor suppressor gene TCF21 on 6q23-q24 in lung and head and neck cancer. Proc Natl Acad Sci U S A 2006;103:982–7.[Abstract/Free Full Text]
- Brambilla E, Travis WD, Colby TV, Corrin B, Shimosato Y. The new World Health Organization classification of lung tumours. Eur Respir J 2001;18:1059–68.[Abstract/Free Full Text]
- Tonon G, Wong KK, Maulik G, et al. High-resolution genomic profiles of human lung cancer. Proc Natl Acad Sci U S A 2005;102:9625–30.[Abstract/Free Full Text]
- Schuebel KE, Chen W, Cope L, et al. Comparing the DNA hypermethylome with gene mutations in human colorectal cancer. PLoS Genet 2007;3:e157.[CrossRef]
- Gentleman RC, Carey VJ, Bates DM, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004;5:R80.[CrossRef][Medline]
- Smyth GK, Speed T. Normalization of cDNA microarray data. Methods 2003;31:265–73.[CrossRef][Medline]
- Smyth GK, Yang YH, Speed T. Statistical issues in cDNA microarray data analysis. Methods Mol Biol 2003;224:111–36.[Medline]
- Takai D, Jones PA. Comprehensive analysis of CpG islands in human chromosomes 21 and 22. Proc Natl Acad Sci U S A 2002;99:3740–5.[Abstract/Free Full Text]
- Valleley EM, Cordery SF, Bonthron DT. Tissue-specific imprinting of the ZAC/PLAGL1 tumour suppressor gene results from variable utilization of monoallelic and biallelic promoters. Hum Mol Genet 2007;16:972–81.[Abstract/Free Full Text]
- Quaggin SE, Schwartz L, Cui S, et al. The basic-helix-loop-helix protein pod1 is critically important for kidney and lung organogenesis. Development 1999;126:5771–83.[Abstract]
- Gros-Louis F, Dupre N, Dion P, et al. Mutations in SYNE1 lead to a newly discovered form of autosomal recessive cerebellar ataxia. Nat Genet 2007;39:80–5.[CrossRef][Medline]
- Fan J, Beck KA. A role for the spectrin superfamily member Syne-1 and kinesin II in cytokinesis. J Cell Sci 2004;117:619–29.[Abstract/Free Full Text]
- Gough LL, Fan J, Chu S, Winnick S, Beck KA. Golgi localization of Syne-1. Mol Biol Cell 2003;14:2410–24.[Abstract/Free Full Text]
- Sjoblom T, Jones S, Wood LD, et al. The consensus coding sequences of human breast and colorectal cancers. Science 2006;314:268–74.[Abstract/Free Full Text]
- Yoon DK, Jeong CH, Jun HO, et al. AKAP12 induces apoptotic cell death in human fibrosarcoma cells by regulating CDKI-cyclin D1 and caspase-3 activity. Cancer Lett 2007;254:111–8.[CrossRef][Medline]
- Su B, Zheng Q, Vaughan MM, Bu Y, Gelman IH. SSeCKS metastasis-suppressing activity in MatLyLu prostate cancer cells correlates with vascular endothelial growth factor inhibition. Cancer Res 2006;66:5599–607.[Abstract/Free Full Text]
- Choi MC, Jong HS, Kim TY, et al. AKAP12/Gravin is inactivated by epigenetic mechanism in human gastric carcinoma and shows growth suppressor activity. Oncogene 2004;23:7095–103.[CrossRef][Medline]
- Flotho C, Paulun A, Batz C, Niemeyer CM. AKAP12, a gene with tumour suppressor properties, is a target of promoter DNA methylation in childhood myeloid malignancies. Br J Haematol 2007;138:644–50.[CrossRef][Medline]
- Chada S, Mhashilkar AM, Ramesh R, et al. Bystander activity of Ad-mda7: human MDA-7 protein kills melanoma cells via an IL-20 receptor-dependent but STAT3-independent mechanism. Mol Ther 2004;10:1085–95.[CrossRef][Medline]
- Lebedeva IV, Sarkar D, Su ZZ, et al. Bcl-2 and Bcl-x(L) differentially protect human prostate cancer cells from induction of apoptosis by melanoma differentiation associated gene-7, mda-7/IL-24. Oncogene 2003;22:8758–73.[CrossRef][Medline]
- Lebedeva IV, Su ZZ, Chang Y, Kitada S, Reed JC, Fisher PB. The cancer growth suppressing gene mda-7 induces apoptosis selectively in human melanoma cells. Oncogene 2002;21:708–18.[CrossRef][Medline]
- Chada S, Bocangel D, Ramesh R, et al. mda-7/IL24 kills pancreatic cancer cells by inhibition of the Wnt/PI3K signaling pathways: identification of IL-20 receptor-mediated bystander activity against pancreatic cancer. Mol Ther 2005;11:724–33.[CrossRef][Medline]
- Zheng M, Bocangel D, Doneske B, et al. Human interleukin 24 (MDA-7/IL-24) protein kills breast cancer cells via the IL-20 receptor and is antagonized by IL-10. Cancer Immunol Immunother 2007;56:205–15.[CrossRef][Medline]
- Huang SM, Stallcup MR. Mouse Zac1, a transcriptional coactivator and repressor for nuclear receptors. Mol Cell Biol 2000;20:1855–67.[Abstract/Free Full Text]
- May P, May E. Twenty years of p53 research: structural and functional aspects of the p53 protein. Oncogene 1999;18:7621–36.[CrossRef][Medline]
- Spengler D, Villalba M, Hoffmann A, et al. Regulation of apoptosis and cell cycle arrest by Zac1, a novel zinc finger protein expressed in the pituitary gland and the brain. EMBO J 1997;16:2814–25.[CrossRef][Medline]
- Abdollahi A, Pisarcik D, Roberts D, Weinstein J, Cairns P, Hamilton TC. LOT1 (PLAGL1/ZAC1), the candidate tumor suppressor gene at chromosome 6q24-25, is epigenetically regulated in cancer. J Biol Chem 2003;278:6041–9.[Abstract/Free Full Text]
- Kamikihara T, Arima T, Kato K, et al. Epigenetic silencing of the imprinted gene ZAC by DNA methylation is an early event in the progression of human ovarian cancer. Int J Cancer 2005;115:690–700.[CrossRef][Medline]
- Kamiya M, Judson H, Okazaki Y, et al. The cell cycle control gene ZAC/PLAGL1 is imprinted—a strong candidate gene for transient neonatal diabetes. Hum Mol Genet 2000;9:453–60.[Abstract/Free Full Text]
- Singhal S, Amin KM, Kruklitis R, et al. Alterations in cell cycle genes in early stage lung adenocarcinoma identified by expression profiling. Cancer Biol Ther 2003;2:291–8.[Medline]
- Morison IM, Ramsay JP, Spencer HG. A census of mammalian imprinting. Trends Genet 2005;21:457–65.[CrossRef][Medline]
This article has been cited by other articles:

|
 |

|
 |
 
S. Ocak, M. L. Sos, R. K. Thomas, and P. P. Massion
High-throughput molecular analysis in lung cancer: insights into biology and potential clinical applications
Eur. Respir. J.,
August 1, 2009;
34(2):
489 - 506.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
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]
|
 |
|

|
 |

|
 |
 
S.P. Barros and S. Offenbacher
Epigenetics: Connecting Environment and Genotype to Phenotype and Disease
Journal of Dental Research,
May 1, 2009;
88(5):
400 - 408.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
M. Tessema and S. A. Belinsky
Mining the Epigenome for Methylated Genes in Lung Cancer
Proceedings of the ATS,
December 1, 2008;
5(8):
806 - 810.
[Abstract]
[Full Text]
[PDF]
|
 |
|