
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Departments of 1 Biochemistry and 2 Preventive Medicine, School of Medicine; 3 Cancer Research Institute; Departments of 4 Internal Medicine, 5 Thoracic Surgery, 6 Pathology, and 7 Obstetrics and Gynecology, Kyungpook National University Hospital; and 8 Department of Internal Medicine, School of Medicine, Keimyung University, Daegu, Korea
Requests for reprints: Jae Yong Park, Department of Internal Medicine, School of Medicine, Kyungpook National University, Samduk 2Ga 50, Daegu, 700-412, Korea. Phone: 82-53-420-5536; Fax: 82-53-426-2046. E-mail: jaeyong{at}kyungpook.ac.kr
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
DNA cytosine methylation in CpG dinucleotides is a major epigenetic mechanism that regulates chromosomal stability and gene expression (3, 4). Many human cancers, including lung cancer, have both global hypomethylation and regional hypermethylation of CpG islands (5-8). Such aberrant DNA methylation may contribute to carcinogenesis in several ways. Hypomethylation may lead to chromosomal instability, reactivation of transposable elements, and loss of imprinting (6, 9). Methylation of CpG sequences may facilitate C-to-T transition mutations in tumor suppressor genes and/or oncogenes through deamination of 5-methylcytosine to thymine (10). Methylated CpG sequences may also increase susceptibility to attack by some environmental carcinogens (11, 12). Finally, de novo hypermethylation of promoter CpG islands may lead to silencing of tumor suppressor genes and DNA repair genes (4, 6, 9).
Methylated CpG sites are recognized by a family of protein factors containing the methyl-CpG binding domain (MBD); to date, five family members have been identified in mammals: MeCP2, MBD1, MBD2, MBD3, and MBD4 (13-15). Four of these proteins (MeCP2, MBD1, MBD2, and MBD3) play important roles for methylation-mediated transcriptional silencing by recruiting chromatin-modifying factors, such as histone deacetylases, to the methylated promoters (14, 15). In contrast to the other family members, MBD4 protein has a thymine glycosylase activity and binds preferentially to 5mCpG-TpG mismatches, which are the primary products of deamination at methyl-CpG. Therefore, MBD4 protein is thought to function as a DNA repair enzyme to minimize mutation at 5-methylcytosine (16, 17).
MBD1 is known to act as a transcriptional repressor through the cooperation of MBD, cysteine-rich CXXC domains, and a COOH-terminal transcription repression domain (18, 19). Among the MBD family of proteins, MBD1 is characterized by two or three cysteine-rich CXXC domains that were originally found in DNA methyltransferase and human trithorax protein HRX (20). MBD1 has at least five isoforms that are the result of alternative splicing within the regions of the CXXC domains and the COOH terminus. These MBD1 isoforms preferentially repress transcription from the methylated gene promoters, but the MBD1 isoforms containing three CXXC domains can also repress transcription from the unmethylated promoters, suggesting that MBD1 plays an important role for the establishment and maintenance of local chromatin states to regulate gene activities (18, 21). In addition to transcriptional gene regulation, a recent study suggests that MBD1 is also involved in DNA repair through its interaction with methylpurine DNA glycosylase, which removes the damaged purines produced by methylating or oxidative agents (22).
Single nucleotide polymorphisms are the most common form of human genetic variation and may contribute to the individual susceptibility for lung cancer. We previously showed that some variants in the DNA repair and DNA methyltransferase (DNMT) genes affect either the expression or the activities of enzymes and are therefore associated with the risk of lung cancer (23-27). Several candidate single nucleotide polymorphisms in the MBD1 gene have been recently deposited in the public databases (http://www.ncbi.nlm.nih.gov/SNP). Although the functional effects of these polymorphisms have not been elucidated, we hypothesized that some of these variants, particularly their haplotypes, may influence MBD1 activity on gene regulation and genome stability, thereby modulating the susceptibility to lung cancer. To test this hypothesis, a case-control study was conducted to evaluate the association between MBD1 genotypes/haplotypes and the risk of lung cancer. Among the candidate single nucleotide polymorphisms in the MBD1 gene, we have focused on the amino acid substitution variants (Pro13Leu, 29295427G>A in exon 3 and Pro401Ala, 29289281C>G in exon 12; Genbank accession no. NT_010966) and the variants in promoter regions (870G>A, 750C>A, 634G>A, and 501delT; Genbank accession no. NT_010966), because these are the most likely to affect gene functioning. In the present study, we evaluated the association of 634G>A, 501delT (501 T/T, T/, /), and Pro401Ala genotypes and their haplotypes with lung cancer, because Pro13Leu, 870G>A, and 750C>A polymorphisms were not detected in a preliminary study that consisted of 27 lung cancer patients and 27 healthy controls.
| Materials and Methods |
|---|
|
|
|---|
MBD1 Genotyping
Genomic DNA was extracted from peripheral blood lymphocytes by proteinase K digestion and phenol/chloroform extraction. The MBD1 634G>A, 501delT, and Pro401Ala genotypes were determined by using a PCR-RFLP assay. The PCR primers were designed based on the Genbank reference sequence (accession no. NT_010966). The PCR primers for the MBD1 634G>A, 501delT, and Pro401Ala polymorphisms were 5'-CCTCTGCCCGTGGGAGGCT-3' (forward) and 5'-CCCTGCAGGAGGCGGAG(mutated T
G)CC-3' (reverse); 5'-TTCCAGCCTCAACCTGAAGG-3' (forward) and 5'-CTCAGTTTACCTGCGGTGTG-3' (reverse); and 5'-AAGCAGATTCCGTGCTGGGA-3' (forward) and 5'-CTTTCCACGACGGTAAGGTC(mutated G
C)-3' (reverse), respectively. The PCR reactions were done in a total volume of 20 µL containing 100 ng genomic DNA, 25 pmol/L of each primer, 0.2 mmol/L deoxynucleotide triphosphates, 75 mmol/L Tris-HCl (pH 9.0), 15 mmol/L ammonium sulfate, 0.1 g/µL bovine serum albumin, 2.5 mmol/L MgCl2, and 1 unit of Taq polymerase (Takara Shuzo Co., Otsu, Shiga, Japan). The PCR cycle conditions consisted of an initial denaturation step at 95°C for 5 minutes followed by 35 cycles of 30 seconds at 94°C; 30 seconds at 56°C for 634G>A, 59°C for 501delT, and 53°C for Pro401Ala; 30 seconds at 72°C; and a final elongation at 72°C for 10 minutes. The PCR products were digested overnight at 37°C with the appropriate restriction enzymes (New England BioLabs, Beverly, MA). The restriction enzymes for the 634G>A, 501delT, and Pro401Ala genotypes were BanII, AvaI, and AvaII, respectively. The digested PCR products were resolved on 6% acrylamide gels and stained with ethidium bromide for visualization under UV light. For quality control, the genotyping analysis was done "blind" with respect to the case/control status. About 10% of the samples were randomly selected to be genotyped again by a different author, and the results were 100% concordant. To confirm the genotyping results, selected PCR-amplified DNA samples (n = 2, respectively, for each genotype) were examined by DNA sequencing, and the results were also 100% concordant.
Promoter-Luciferase Constructs
To examine the potential effects of the 634G>A and 501delT polymorphisms on the MBD1 transcription activity, the fragments, including the 634G>A (872 to 576, transcription start site of exon 1 counted as +1) and the 501delT (601 to 290), were synthesized, respectively, by PCR using the genomic DNA from donors that were carrying either the wild-type allele or the polymorphic allele of each polymorphism. The PCR primers used for the 634G>A and the 501delT promoter regions were 5'-GAAGCTGTC TCCACATTGCT-3' (forward) and 5'-CACACCGCAGGTAAACTGA-3' (reverse) and 5'-CGCGTGCCTCAGTTTACCT-3' (forward) and 5'-CGCTTTCCAGCCTCAACCT-3' (reverse), respectively. The PCR products were inserted upstream of the SV40 promoter in the pGL3-promoter plasmid (Promega, Madison, WI). For comparison of the haplotypes' promoter activities (haplotype 634G/501T versus haplotype 634A/501), the fragments of the MBD1 promoter region (994 to +113) were synthesized by PCR using the genomic DNA from donors carrying either the 634G/501T haplotype or the 634A/501 haplotype. The PCR primers for the MBD1 promoter were 5'-GCCCATCTGCCTTTATCAGA-3' (forward) and 5'-CAGCCCCAAGGCTGTCT-3' (reverse). The PCR products were inserted upstream of the luciferase gene in the pGL3-basic plasmid (Promega). The correct sequence of all the clones was verified by DNA sequencing.
Transient Transfection and Luciferase Assay
The promoter activity was measured using the Dual-Luciferase Reporter Assay System (Promega) in Chinese hamster ovary cells and A549 lung cancer cells. The Chinese hamster ovary cells were cultured in minimal essential medium with 10% heat-inactivated fetal bovine serum, and the A549 cells were grown in RPMI 1640 that was supplemented with 10% heat-inactivated fetal bovine serum. Cells (1 x 105) were plated in a six-well plate the day before transfection so that the cells were
60% confluent by the next day. The pRL-SV40 plasmid and the pGL3-basic plasmid with the synthesized fragments of the MBD1 promoter region were cotransfected using Lipofectine reagent (Invitrogen, Carlsbad, CA). The pRL-SV40 vector that provided the constitutive expression of Renilla luciferase was used as an internal control to correct for the differences in transfection and harvesting efficiency. The cells were collected 48 hours after transfection, and the cell lysates were prepared according to Promega's instruction manual. Luciferase activity was measured using a Lumat LB953 luminometer (EG&G Berthhold, Bad Wildbad, Germany), and the results were normalized using the activity of Renilla luciferase. Independent triplicate experiments were done four times, and the results were reported as mean ± SD.
Statistical Analysis
The cases and controls were compared using Student's t test for continuous variables and the
2 test for categorical variables. Hardy-Weinberg equilibrium was tested for with a goodness-of-fit
2 test with one degree of freedom to compare the observed genotype frequencies among the subjects with the expected genotype frequencies. We examined the widely used measure of linkage disequilibrium between pairs of biallelic loci, Lewontin's D' (|D'|; ref. 30). The haplotypes and their frequencies were estimated based on a Bayesian algorithm using the Phase program (31), which is available at http://www.stat.washington.edu/stephens/phase.html. Unconditional logistic regression analysis was used to calculate the odds ratios (ORs) and 95% confidence intervals (95% CI), with adjustment for possible confounders (gender and family history of lung cancer as a nominal variable and age and pack-years as continuous variables). In addition to the overall association analysis, we did a stratified analysis by age, gender, smoking status, and tumor histology to further explore the association between MBD1 genotypes/haplotypes and the risk of lung cancer in each stratum. To test which one of the three polymorphisms is more likely to be the main cause of the observed association, we compared seven different logistic regression models (each polymorphism alone, any two of the three polymorphisms, and then all the three polymorphisms together). When multiple comparisons are made, the corrected Ps (Pc) were also calculated for multiple testing by using Bonferroni's inequality method. All analyses were done using Statistical Analysis Software for Windows, version 8.12 (SAS Institute, Cary, NC).
| Results |
|---|
|
|
|---|
|
|
|
|
The association between the MBD1 genotypes/haplotypes and the risk of adenocarcinoma was further examined after stratifying for age, gender, smoking status, and histologic subtypes of adenocarcinoma [adenocarcinoma with/without bronchioloalveolar histology (BAC)]. There was no clear evidence that age, gender, or smoking status modified the effect of the MBD1 genotypes/haplotypes on the risk of adenocarcinoma in the stratified analysis. Moreover, the stratified analysis on the histologic subtypes of adenocarcinoma did not show any difference between adenocarcinomas with and without BAC histology (data not shown).
We investigated the effects of the 634G>A and 501delT polymorphisms on the promoter activity of MBD1 by luciferase assay. The promoter activity of the 634A allele was significantly higher (1.5-fold) compared with the 634G allele in the Chinese hamster ovary cells (P < 0.001; Fig. 1A), but the 501delT polymorphism did not have an effect on the promoter activity (Fig. 1B). Because the 634G>A and 501delT polymorphisms were in linkage disequilibrium, we also compared the transcription activity of the haplotypes (haplotype 634G/501T versus haplotype 634A/501). In the Chinese hamster ovary cells, the 634A/501 haplotype increased transcription activity by 2.1-fold compared with the 634G/501T haplotype (P < 0.001). Similar result was observed in the A549 cells (P < 0.01; Fig. 1C).
|
| Discussion |
|---|
|
|
|---|
In the present study, carriers with the 634G/501T/401Pro haplotype were at increased risk of lung cancer compared with individuals having the 634A/501/401Ala haplotype. To determine whether the association between the MBD1 polymorphisms and the risk of lung cancer is due to differences in the transcriptional activity of MBD1 promoter, we compared the promoter activity of the wild-type allele or the polymorphic allele of these two polymorphisms by luciferase assay. In vitro promoter assay revealed that the 634A allele had a significantly higher transcriptional activity than the 634G allele, and the 501delT polymorphism did not have an effect on the transcriptional activity of the MBD1 promoter. When we compared the promoter activity of the MBD1 haplotypes, the 634A/501 haplotype had significantly higher transcriptional activity than the 634G/501T haplotype. These results suggest that the genetic effects of MBD1 polymorphisms on the risk of lung cancer can be mainly attributed to the 634G>A polymorphism, and these findings also suggest that the 634G>A change influences MBD1 expression, thus contributing to the genetic susceptibility to lung cancer. The mechanism by which the MBD1 634A/501 haplotype leads to higher promoter activity is currently unknown. Analysis of the potential transcription factor binding sites by the Alibaba2 program (32) showed that the 634G>A change leads to the creation of an cytoplasmic polyadenylation element binding site and eliminates an activator protein 2
binding site and thyroid hormone receptor-
binding site, whereas the 501delT polymorphism has no effects on the transcription binding site of the MBD1 promoter. Therefore, it is possible that the predicted changes in the putative transcription factor binding sites due to the 634G>A change may lead to enhanced promoter activity. However, this hypothesis has to be verified in future studies.
MBD1 has been thought to be involved in mediating methylation-associated gene silencing in various human cancers (18, 19). Therefore, it is possible that up-regulation of MBD1 expression may lead to a predisposition towards silencing of the methylated genes, thereby increasing the susceptibility to lung cancer. In contrast to this hypothesis, however, we observed in the present study that the MBD1 genotype/haplotype with lower transcriptional activity was associated with an increased risk of lung cancer. In fact, this finding is supported by several previous reports (33-36) showing that MBD proteins were down-regulated in a variety of human cancers. Although the mechanism by which down-regulation of MBD1 expression leads to the enhanced susceptibility to lung cancer remains to be elucidated, this may happen because a decrease of MBD1 expression leads to a disturbance in the normal cellular homeostasis of the gene expression profiles; thus, it allows tumorigenesis via epigenetic instability (33-36). Another possible explanation is that MBD1 may function as a DNA repair system that is associated with methyl-CpG dinucleotides (22), as is the case of MBD4, thereby influencing the susceptibility to lung cancer.
Another interesting finding of the present study is that the MBD1 polymorphisms, and their haplotypes had a more pronounced association with adenocarcinoma. Although the reason for the observed histology-specific difference in the risk conferred by the MBD1 polymorphisms remains to be elucidated, this difference may be attributable to the differences in the carcinogenic pathways among the histologic types of lung cancer. Therefore, certain genotypes could confer a greater susceptibility to a particular histologic type of lung cancer (37-39). In our previous study (26), we have observed a similar finding that DNMT3B polymorphisms and their haplotypes were significantly associated with only the risk of adenocarcinoma. These findings of our previous and present studies suggest that alterations of the DNA methylation machinery, including DNMT3B and MBD1, may have a pronounced association with development of adenocarcinoma.
Several recent studies have shown that mutations in the kinase domain of epidermal growth factor receptor gene, like the K-ras mutations, frequently target adenocarcinoma but are more frequent in never smokers, females, adenocarcinomas with nonmucinous BAC histology, and the East Asian populations, whereas the K-ras mutations are more frequent in smokers, males, adenocarcinomas with mucinous BAC histology, and the Western populations (40, 41). These observations suggest that adenocarcinomas arising in the never smokers and smokers may be caused by different etiologies, not only in relation to environmental risk factors but also in relation to genetic susceptibility factors (40-42). Therefore, we did a stratification analysis to examine if the MBD1 genotypes/haplotypes may have differential effects on the risk of adenocarcinoma according to age, gender, and smoking status and the histologic subtypes of adenocarcinoma (adenocarcinoma with/without BAC histology). In the current study, no risk modification was found with regard to age, gender, and smoking status. Moreover, the stratified analysis on the histologic subtypes of adenocarcinoma did not show any difference between adenocarcinomas with and without BAC histology (data not shown). However, because the number of subjects in the subgroups was small, our findings from the stratified analyses should be interpreted with caution before being confirmed in larger studies.
In the current study, the MBD1 Pro13Leu, 870G>A, and 750C>A polymorphisms were not detected in the preliminary study that included 27 healthy controls. These samples included 54 chromosomes, which provides at least a 95% confident level to detect alleles with frequencies >5%. Thus, it is very likely that if these polymorphisms exist, they may not play a major role in the genetic susceptibility to lung cancer in the Korean population (43, 44).
In conclusion, we found that the MBD1 634G>A, 501delT, and Pro401Ala polymorphisms and their haplotypes were significantly associated with the risk of lung cancer and particularly adenocarcinoma. These results suggest that the MBD1 gene may be involved in the development of lung cancer, although additional studies having larger sample sizes are required to confirm our findings. Future studies on the other MBD1 sequence variants and their biological function are also needed to understand the role of the MBD1 polymorphisms in determining the risk of lung cancer. Moreover, because genetic polymorphisms often vary between different ethnic groups, further studies are needed to clarify the association of the MBD1 polymorphism with lung cancer in diverse ethnic populations.
| Footnotes |
|---|
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.
Note: J-S. Jang and S.J. Lee contributed equally to this work.
Received 6/ 9/05; revised 8/16/05; accepted 9/ 1/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. M. Allan, X. Liang, Y. Luo, C. Pak, X. Li, K. E. Szulwach, D. Chen, P. Jin, and X. Zhao The loss of methyl-CpG binding protein 1 leads to autism-like behavioral deficits Hum. Mol. Genet., July 1, 2008; 17(13): 2047 - 2057. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Miremadi, M. Z. Oestergaard, P. D.P. Pharoah, and C. Caldas Cancer genetics of epigenetic genes Hum. Mol. Genet., April 15, 2007; 16(R1): R28 - R49. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. C. Shin, S. J. Lee, J. E. Choi, S. I. Cha, C. H. Kim, W. K. Lee, S. Kam, Y. M. Kang, T. H. Jung, and J. Y. Park Glu346Lys Polymorphism in the Methyl-CpG Binding Domain 4 Gene and the Risk of Primary Lung Cancer Jpn. J. Clin. Oncol., August 1, 2006; 36(8): 483 - 488. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Y. Park, J. M. Park, J. S. Jang, J. E. Choi, K. M. Kim, S. I. Cha, C. H. Kim, Y. M. Kang, W. K. Lee, S. Kam, et al. Caspase 9 promoter polymorphisms and risk of primary lung cancer Hum. Mol. Genet., June 15, 2006; 15(12): 1963 - 1971. [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 |