
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Departments of 1 Genetics and Complex Diseases and 2 Biostatistics, Harvard School of Public Health and 3 Department of Environmental Health, Boston University School of Public Health, Boston, Massachusetts
Requests for reprints: Karl T. Kelsey, Department of Genetics and Complex Diseases, Harvard School of Public Health, 665 Huntington Avenue, Boston, MA 02115. Phone: 617-432-3313; Fax: 617-432-0107. E-mail: kelsey{at}hsph.harvard.edu.
| Abstract |
|---|
|
|
|---|
Methods: Global DNA methylation levels were assessed using a modified version of the combined bisulfite restriction analysis of the LRE1 sequence in a population-based case-control study of HNSCC from the Boston area.
Results: Hypomethylation lead to a significant 1.6-fold increased risk for disease (95% confidence interval, 1.1-2.4), in models controlled for other HNSCC risk factors. Smoking showed a significant differential effect (P < 0.03) on blood relative methylation between cases and controls. Furthermore, in cases, variant genotype in the MTHFR gene and low folate intake showed relationships with decreased global methylation, whereas in controls, antibody response to human papillomavirus 16 was associated with an increased global methylation level.
Discussion: DNA hypomethylation in nontarget tissue was independently associated with HNSCC and had a complex relationship with the known risk factors associated with the genesis of HNSCC. (Cancer Epidemiol Biomarkers Prev 2007;16(1):10814)
| Introduction |
|---|
|
|
|---|
11,000 deaths (1). Although survival rates can be as high as 75% to 80% if HNSCC is detected in its early stages, the majority of patients are not diagnosed until the disease has progressed to an advanced stage, reducing the chances of survival after 5 years to 35%, a figure which has improved little in the past 25 years (2-4). Alcohol and tobacco use are long-established risk factors for head and neck cancers and are known to act synergistically (5-10). Infection with high-risk human papillomavirus (HPV) types, particularly HPV16, has been associated with increased risk for HNSCC, and DNA from this virus has been detected in tumors of HNSCC patients (11-16). Antibodies to the L1 protein of HPV16 reflect the presence of viral DNA and have also been associated with risk for HNSCC.4 Additionally, dietary factors, including deficiencies in dietary folate, are also hypothesized risk factors for HNSCC (17, 18).
Changes in methylation patterns, particularly promoter-specific hypermethylation and global (genome-wide) hypomethylation, are thought to contribute to neoplasia and tumor growth (19, 20). Gene promoter hypermethylation in tumor tissues is a common event in the development of many types of cancer, including HNSCC (21-23). This is because neoplastic growth is frequently preceded by aberrant promoter methylation of tumor suppressor genes, which leads to a loss of function that promotes cell proliferation (24). Cancer-linked global genomic hypomethylation in tumor tissue is a common characteristic in a wide variety of malignancies, ranging from solid tumors, such as breast, colon, oral, and lung cancers, to cancers of the blood (25-29). Whereas hypermethylation occurs chiefly in gene promoter regions, global hypomethylation occurs not only in transcription control regions, such as promoters, but also in repetitive DNA sequences, such as heterochromatic regions and retrotransposons (25). Hypomethylation is thought to contribute to carcinogenesis by inducing genomic instability (30, 31), thereby causing the formation of abnormal chromosomal structures (25, 32). Evidence also points to the possibility of this alteration activating or enhancing oncogene expression because (a) hypomethylation of promoter regions of certain genes increases the target gene expression (33) and (b) promoter hypomethylation is linked to global methylation levels (34). Furthermore, the causal role of hypomethylation in carcinogenesis has been established using mouse models with decreased methyltransferase activity (31).
The role of folates in carcinogenesis is also strongly linked to genomic methylation levels. Early studies showed that global hypomethylation occurs in rats with low-folate diets (35). A study comparing global DNA methylation levels of tumor tissues in squamous cell lung cancer with those in adjacent, unaffected tissues concluded that hypomethylation is associated with folate deficiency in both diseased and unaffected tissues, supporting the assertion that folate is necessary for proper DNA methylation (36). This study also suggests that although methylation levels are more pronounced in tumor tissues and are tissue specific, proxy tissues may be indicators of local methylation levels.
The family of LINE1 (long interspersed nuclear elements) retrotransposons is reportedly hypomethylated in many cancers and reflects global methylation status in the genome (28, 37), thus examination of methylation at LINE1 regions has served as a proxy for measuring global methylation levels. One long interspersed nuclear element repeat region, LRE1, located on 22q11-q12, is a consistent indicator of global methylation status (28, 37, 38).
It is increasingly clear that epigenetics plays a causal role in cancer development. HNSCC is a useful disease for studying global hypomethylation because the mechanisms of epigenetic maintenance are related to the risk factors linked to this disease. Examining hypomethylation in this context may shed light on the means by which these factors contribute to carcinogenesis, deepening the general understanding of the complex interaction between epigenetics and environmental exposures in cancer development. Therefore, we have examined, in a population-based case-control study of HNSCC, the association between global genomic methylation, measured as LRE1 methylation status, and HNSCC as well as have examined the associations between risk factors associated with HNSCC and LRE1 methylation status.
| Materials and Methods |
|---|
|
|
|---|
3.5 million people in 249 cities and towns within a 1-h drive of Boston. The institutional review boards at all participating institutions approved this study, and all volunteer participants provided informed consent. Briefly, incident cases of HNSCC were identified through multidisciplinary head and neck clinics, otolaryngology, and radiation oncology departments at nine medical facilities located in Boston, Massachusetts. We defined HNSCC as including International Classification of Disease Ninth Revision codes 141, 143-6, 148, 149, and 161. All patients with carcinoma in situ, lip, salivary gland, or nasopharyngeal cancer or recurrent cancer of the head and neck region were excluded. Histologic classification of malignancy was based on that reported by pathology at the participating hospitals. Population-based controls were drawn from the specified greater Boston population. The controls were frequency matched (1:1) to cases by age (±3 years), gender, and town of residence. These controls were identified through random selection from the resident lists for the 249 cities and towns within the study area using the address of the cases as reference. Participating cases and controls were given a self-administered questionnaire to collect medical history, demographic information, as well as information on tobacco and alcohol consumption. Each questionnaire was reviewed with each participant by a trained research coordinator. Smoking history was ascertained with a standardized instrument that assesses the number of years smoked, the number of cigarettes smoked daily, age at which an individual started smoking, number of years since quitting, and the duration of smoking in a decade-specific manner. Similar information was obtained about lifetime consumption of beer, wine, and liquor. The Willett food frequency questionnaire, a standardized and validated instrument, was used to assess diet history (40, 41). Subjects are asked to recall their usual diet over a 1-year time period 5 years previous to their diagnosis (in a fashion consistent with the validation study of diet recall using this instrument). Because many tumors of the oral cavity and pharynx may affect food intake before being diagnosed, moving the recall period well in advance of tumor development eliminates this possible bias. Questionnaires were given to case participants during an initial clinical visit and subsequently retrieved in-person. Control participants received their questionnaires in the mail and returned them in person to the research assistant.
Eight hundred and twenty-three eligible cases were invited to participate, of these 57 refused to participate. Among the 766 consented subjects, another 44 did not complete the questionnaire. Of these participants, complete questionnaire data, including food frequency questionaire and blood samples were available on 278 cases. Similarly, 1,623 subjects were identified and eligible for participation as controls. Eight hundred twenty-eight refused to participate, 815 subjects were consented, and 771 finally enrolled in the study. Six of the controls were withdrawn as they were matched to a case that became ineligible, such that 765 controls were enrolled and completed. Complete questionnaire data and blood samples were available on 526 people. The characteristics of the final study population are described in Table 1 .
|
LRE1 Relative Methylation Assay
To assess the relative methylation status of LRE1 in DNA derived from whole blood samples, we used a modified version of the previously described combined bisulfite restriction analysis LRE1 assay, allowing for fluorescence-based relative quantification of differentially methylated product (38). Sodium bisulfite-modified DNA was subjected to PCR using the fluorescently labeled primers Line1-F (5'-Hex-CCGTAAGGGGTTAGGGAGTTTTT-3') and Line1-R (5'-Hex-RTAAAACCCTCCRAACCAAATATAAA-3'). Amplification conditions were as follows: bisulfite-modified DNA, 1 nmol/L each of forward and reverse primers, 2 mmol/L of each deoxynucleotide triphosphate, 1.5 µmol/L MgCl2, 1x PCR buffer (Applied Biosystems, Foster City, CA), and 1.25 units AmpliTaq polymerase were combined in a final volume of 50 µL. Reactions were incubated for 5 min at 95°C and then for 40 cycles at 94°C denaturing for 30 s, 50°C annealing for 30 s, and 72°C extension for 30 s followed by a final 7-min extension at 72°C. To confirm amplification of the 160-bp product, 10 µL of the completed PCRs were resolved using a 3% agarose gel in 1x Tris-borate EDTA and visualized with ethidium bromide.
The remaining reaction volumes were then used for double digestion with 2 units TaqI and 8 units of TasI enzyme in 1x TaqI buffer (Fermentas, Hanover, MD) at 65°C for 18 to 24 h in the dark. Methylated amplicons are TaqI positive and yield two 80-bp fragments, whereas unmethylated amplicons, in which cytosines are converted to uracils during bisulfite modification, are TasI positive and yield 63- and 97-bp fragments. Fragments were prepared for analysis by mixing 1.5 µL of the digested product with 11.5 µL DI formamide and 0.5 µL Genescan 350 TAMRA size standard (Applied Biosystems) and then denatured at 95°C for 5 min. Samples were resolved using capillary electrophoresis in the ABI Prism 310 Genetic Analyzer. The accompanying Genescan software was used to determine peak heights for restriction products generated by TasI (63 and 97 bp) from unmethylated amplicons and TaqI (80 bp) from methylated amplicons. The degree of relative methylation was assessed by taking the ratio of the methylated peak height (80 bp) to the sum of peak heights from all digestions (63, 80, and 97 bp).
Controls were generated by using unmodified primers on unmodified DNA derived from whole blood to amplify a 390-bp region surrounding the 160-bp region of interest in LRE1, thus generating an unmethylated control, as the PCR product does not preserve methylation from the original template. Positive control DNA was generated by treating whole-blood DNA with DNA methylase, effectively methylating all available CpG islands. The methylated and unmethylated products were then modified with sodium bisulfite, subjected to PCR under the conditions described above, and the subsequent products were mixed in varying proportions to generate a standard curve under which samples were expected to decrease. Other controls used to validate this assay were DNA from the K-562 cell line (highly unmethylated) and sperm DNA (highly methylated).
Statistical Analysis
Data were analyzed using SAS statistical software. The Wilcoxon rank-sum test was used as a nonparametric comparison of median methylation levels in cases versus controls. To evaluate the effects of individual variables, including LRE1 methylation levels on case-control status, while controlling for intervariable confounding, unconditional logistic regression was used to determine odds ratios and their associated 95% confidence intervals. In this analysis, the LRE1 relative methylation level was broken into terciles based on the distribution in controls. Lifetime smoking history was similarly broken into terciles based on the distribution in controls, and lifetime average number of alcoholic drinks per week was broken into quartiles based on the distribution in controls. HPV16 serology was considered positive in all subjects whose HPV16 titer was greater than the limit of detection (>12 milliMerck units). Dietary folate was first broken into terciles based on the distribution in controls and then dichotomized to examine the effect of dietary folate among the lowest 33rd percentile of consumption.
To examine those factors that may be etiologically related to LRE1 relative methylation level, we did stratified analyses in cases and controls, as there may be fundamental differences in the biology of this alteration between these groups. For these analyses, we modeled LRE1 relative methylation as a linear function of covariates, using generalized linear models,with settings appropriate to a ß response: an identity link function, the binomial variance Var(µ) = µ(1 µ), and a scale variable estimated with Pearson residuals. The fit of the model was assessed using the cumulative residual method of Lin et al. (42), with 10,000 resampled residual processes, to assess goodness-of-fit. We also reparameterized the stratified model as an interaction model, including interaction terms for each of the covariates with case status, to assess whether differences in the effects of the covariates between cases and controls was statistically significant.
| Results |
|---|
|
|
|---|
T polymorphism. In assessing the association between HNSCC and dietary folate intake, we compared the lowest tercile of intake, <418 µg/d (similar to the U.S. recommended daily allowance level of 400 µg/d) witho the upper two terciles. Figure 1 depicts the distribution of the LRE1 relative methylation level in cases and controls as box plots of the data. The median methylation level of 0.753 in controls was slightly but significantly higher than the median level of 0.747 in cases (P < 0.03, Wilcoxon rank-sum test).
|
|
|
1% reduction in LRE1 relative methylation. At the same time, each one pack-year increase in lifetime pack-years of smoking was significantly associated (P < 0.04) with an
0.02% increase in the LRE1 relative methylation level among cases. In a nonstratified model examining the interaction of each of the covariates, only smoking history, measured as pack-years, showed a statistically significant interaction with case status (P < 0.03).
Controls, on the other hand, showed a different group of significant predictors for LRE1 relative methylation (Table 3). In a model controlled for age, lifetime average alcoholic drinks per week, lifetime pack-years smoked, dietary folate intake, and MTHFR codon 677 genotype, we observed a significant (P < 0.002) 1% increase in the relative LRE1 methylation level in males compared with females and significant
1% increases in LRE1 relative methylation for subjects with positive HPV16 antibody serology and for subjects of non-Caucasian race compared with Caucasians (P < 0.02 and P < 0.03, respectively).
| Discussion |
|---|
|
|
|---|
In stratified models examining predictors of LRE1 relative methylation, we observed that gender was significantly related to methylation status, with women being more likely to have reduced methylation level. Although the precise explanation for this is not clear, differences between male and female nutrient intake and loss might explain this distribution. Women may tend to have diets lower in protein or lower caloric intake in general, resulting in comparatively lower levels of nutrients, such as choline, methionine, homocysteine, and folate; these macromolecules are all metabolically related in the process of DNA methylation (43). In addition, folate is essential for erythrocyte formation and development. Because menstruation regularly depletes their supply of erythrocytes, females may tend to have a higher folate requirement, or conversely, lower levels of circulating folate. This may contribute to the difference in global methylation between males and females. In postmenopausal females, who do not have a regular depletion of folate-requiring erythrocytes, other dietary factors and behavioral factors may contribute to low methylation.
The significantly higher levels of relative methylation observed in non-Caucasian controls may also reflect effects of different environmental factors, such as dietary intake, which can contribute to higher global methylation levels, or equally likely, genetic predisposition to constitutively higher global methylation.
Both the association of an increasing LRE1 methylation level among non-Caucasian controls as well as among men is anomalous when compared with the relationship of these factors with disease risk; non-Caucasians and men are at higher risk for HNSCC. This suggests that these factors, in particular, but others as well, do not work in a linear pathway through global hypomethylation to elicit disease. Instead, these data might imply that global methylation functions in a complex pathway to influence individual cancer risk, perhaps interacting with other known and yet-to-be-known factors. Additional studies are needed to more carefully examine how gender and race may influence LRE1 methylation levels and, within these subgroups, how additional factors, such as diet and lifestyle, influence global methylation.
It is also of interest that, in our data, the relative LRE1 methylation level was significantly associated with detectable HPV16 antibodies in controls, with those subjects with positive serology having increased relative methylation. It has been observed that the HPV genome itself may become hypermethylated on infection or integration into the host (44, 45) and that the carcinogenicity of HPV16 is then related to hypomethylation of its genome during neoplastic progression (46) There is also evidence, in vitro, that immortalization by HPV leads to specific gene hypermethylation (47). We also observed previously an association between gene-specific hypermethylation of the SFRP4 gene promoter and HPV16 presence, particularly among nonsmokers. This may reflect a mechanism whereby aberrant promoter methylation results from mistargeted host defense methylation during viral integration or the genomic instability attributed to HPV16 presence (48). Our observation of greater relative methylation in controls with HPV16 antibodies, which are associated with the presence of viral DNA,4 also may explain, in part, why these individuals were not susceptible to the carcinogenic effects of the HPV16 virus; increased or more stable global methylation levels may reduce the carcinogenic potential of the HPV16 genome, thus preventing viral neoplastic progression in these individuals. Additionally, higher LRE1 methylation levels, particularly in lymphocytes, may be related to the ability of the individual to mount an immune response to HPV (measured by HPV serology), thus marking individuals who successfully cleared the HPV infection before carcinogenesis. Therefore, these results suggest that susceptibility to HPV16 carcinogenesis may be an epigenetically modified process.
In contrast to the data in controls, the data from cases show that low daily dietary folate intake and possessing the variant T allele at MTHFR codon 677 are correlated with reduced LRE1 relative methylation. Given the relationship between MTHFR function and levels of circulating folate (in the form of 5-methyl-THF), this is to be expected, particularly because various MTHFR polymorphisms have been implicated as risk factors for HNSCC (49). Preliminary analysis of our data has suggested that MTHFR 677genotype is not an independent risk factor for HNSCC, but its influence on hypomethylation may suggest the possibility for effect modification by this or other factors. Although low folate has been shown as a risk factor for HNSCC (50, 51), the mechanism of this association has not been fully elucidated. One possible explanation is that a change in the physiology of HNSCC patients, compared with the baseline physiologic characteristics of controls (or a general population), causes methylation levels to become highly folate dependent. A rapidly dividing population of cells, such as a tumor, has a high nutrient requirement. With elevated levels of protein synthesis and DNA replication, requirements for methionine, folate, and homocysteine will increase, upsetting the ratio of cofactors necessary to maintain balance and proper functioning of the methylation process throughout the body and particularly in dividing cells, such as lymphocytes. Essentially, because the folate requirement is increased in cancer patients, but their dietary patterns will not increase to fill that requirement (and may in fact lessen due to the disease), more patients are likely to decrease below that elevated threshold. The methylation status of these individuals becomes even more perturbed as cells struggle to function with lower folate levels. Research has shown that whereas folate supplementation in healthy cells is generally protective against the development of tumors, in epigenetically disrupted cancerous cells, supplementation in fact may increase tumor growth (43, 52, 53), supportive of the results obtained here.
To more closely examine if these differences in the effects of specific covariates on LRE1 relative methylation truly differ in cases and controls, we did a nonstratified analysis, using the same generalized linear model, but included terms for the interactions of each of the covariates with case status. According to this approach, only lifetime smoking history showed a significant interaction at the P < 0.05 level, suggesting that it is truly acting differentially to predict LRE1 relative methylation in cases and controls. The nonsignificant results among the other predictors which seem to significantly affect the LRE1 relative methylation differentially in the stratified models suggest that we are underpowered to examine these interactions. We also found that race and HPV16 serology were differentially associated with relative LRE1 methylation, having opposite effects on LRE1 relative methylation in cases and controls. These variables predict reduced LRE1 relative methylation in cases but increased LRE1 relative methylation in controls. Male gender, on the other hand, seems to have a stronger positive relationship to LRE1 relative methylation in controls (estimate of effect, 0.012) than in cases (estimate of effect, <0.0001). Similarly, both low dietary folate intake and variant MTHFR C677T genotype show stronger negative relationships with LRE1 relative methylation in cases compared with controls.
Therefore, we believe that our data from the stratified analyses in cases and controls, as well as our overall model of HNSCC risk, suggest that LRE1 relative methylation is an independent epigenetic biomarker of HNSCC. Feinberg, et al. (54) have proposed that cancer may derive from epigenetic progenitors and that epigenetic status, unlike genetic susceptibility, which is defined by inherited polymorphisms, may be more plastic. Our results build on this model, suggesting that epigenetic biomarkers, such as LRE1 relative methylation, can potentially be influenced by a variety of factors and may thus aid in explaining the modification of effect of these factors on risk for disease.
| Acknowledgments |
|---|
| 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.
4 C.S. Furniss, M.D. McClean, J.F. Smith, et al. Human papillomavirus 16 and head and neck squamous cell carcinoma, International Journal of Cancer, in press. ![]()
Received 7/28/06; revised 9/21/06; accepted 10/13/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. J. Marsit, C. C. Black, M. R. Posner, and K. T. Kelsey A Genotype-Phenotype Examination of Cyclin D1 on Risk and Outcome of Squamous Cell Carcinoma of the Head and Neck Clin. Cancer Res., April 15, 2008; 14(8): 2371 - 2377. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. S. Furniss, C. J. Marsit, E. A. Houseman, K. Eddy, and K. T. Kelsey Line Region Hypomethylation Is Associated with Lifestyle and Differs by Human Papillomavirus Status in Head and Neck Squamous Cell Carcinomas Cancer Epidemiol. Biomarkers Prev., April 1, 2008; 17(4): 966 - 971. [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 | Cell Growth & Differentiation |