
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Imaging, Diagnosis, Prognosis |
Authors' Affiliations: Departments of 1 Pathology, 2 Leukemia, 3 Head and Neck Surgery, and 4 Thoracic/Head and Neck Medical Oncology, The University of Texas M. D. Anderson Cancer Center, Houston, Texas and 5 Department of Pathology, Ansan Hospital, Korea University, Ansan, Gyeonggi-Do, South Korea
Requests for reprints: Adel K. El-Naggar, Department of Pathology, Unit 85, The University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, TX 77030-4009. Phone: 713-792-3109; Fax: 713-792-5532; E-mail: anaggar{at}mdanderson.org.
| Abstract |
|---|
|
|
|---|
Experimental Design: The hypermethylation status of the retinoid acid receptor β2 (RARβ2), RAS association domain family 1A (RASSF1A), O6-methylguanine-DNA methyltransferase (MGMT), and E-cadherin genes was analyzed in five salivary carcinoma cell lines and 69 human salivary gland carcinoma specimens by pyrosequencing and MSP techniques. The two datasets were compared by linear regression. Correlations between methods and with clinicopathologic characteristics were assessed by Pearson's
2 test or the two-tailed Fisher exact test, as applicable, using cutoff points determined from the regression curves and empirical fitting. We also investigated the effect of demethylating agents on methylated genes in cell lines to assess their effect on the expression of these genes.
Results: Overall, regression analysis indicated high degrees of correlation of the two methods for measurement of methylation for the RARb2, RASSF1A, and MGMT genes (adjusted R2 = 0.319, 0.835, and 0.178; P < 0.001, <0.001, and 0.0002, respectively) among the 69 tumors tested. However, the pyrosequencing technique yielded four more instances of methylation above background levels than MSP for RARβ2 and three more for RASSF1. Methylation of either RARβ2 and RASSF1A alone or both by pyrosequencing were correlated with tumor type (P = 0.027, 0.014, and 0.012, respectively). Methylation of RARβ2 alone and in combination with RASSF1A by pyrosequencing were also significantly correlated with tumor grade (P = 0.014 and 0.011, respectively) and 3-year survival (P = 0.002 and 0.004, respectively). The survival curves of patients who had hypermethylation at both RARβ2 and RASSF1A were significantly lower than those of patients who had hypermethylation at neither or just for the RASSF1A (P = 0.008 and 0.007, respectively). 5-Azadeoxycytidine treatment of methylated cell lines led to the reactivation of RARβ2 expression in only one of the five cell lines.
Conclusions: (a) Although the methylation status of RARb2, RASSF1A, and MGMT genes by both techniques were significantly correlated, pyrosequencing is generally more sensitive and its results correlate better with the clinical variables than those of MSP. (b) The methylation level of the RARβ2 and/or RASSF1A by pyrosequencing is significantly associated with aggressive tumor phenotypes and patients survival.
10% of all head and neck cancers. Known for their marked phenotypic heterogeneity and wide range of biological behavior (1), these tumors can be categorized on the basis of their cytologic and architectural features into two broad classes, high and low to intermediate categories (2). Conventional primary therapy for SGCs is complete surgical resection, but for advanced surgically unresectable and recurrent tumors, the therapeutic options are limited and patient mortality is high (3). In search of objective variables for the biological stratification of patients with such advanced tumors, we and others have studied several biomarkers (4–6), but none of them has been validated for clinical utility. Hypermethylation of CpG-enriched sites, CpG islands, at the promoter and the first exon of genes is an important mechanism in the functional suppression of genes associated with mammalian development and cancer (7–9). This process entails the enzymatic addition of a methyl group to the fifth position of the deoxycytidine residue in each affected CpG dinucleotide (8–10). A number of reports have revealed an association between the aberrant hypermethylation of these CpG islands in several tumor suppressor genes and tumor development and thus implicated methylation as an alternative epigenetic mechanism of gene inactivation (11–15). In contrast to the integral nature of DNA-related alterations, epigenetic modification is a dynamic and reversible process of considerable therapeutic implications (16–20).
Techniques for assessing promoter methylation status vary considerably. Several methods that assess this feature have been reported (21–27), but the most widely used is the methylation-specific PCR (MSP) analysis of DNA after bisulfite treatment (24). However, the popularity of the MSP methods (whether nested or unnested) is tempered by notable shortcomings, including its qualitative nature, the limited CpG sites targeted, lack of internal control, and false-positive signals at high PCR cycles (25, 28). Recent attempts to remedy some of these deficiencies have led to the development of an alternative high-throughput technique for methylation analysis, known as pyrosequencing (27, 29, 30). Pyrosequencing, a sensitive real-time sequencing-by-synthesis method is based on the conversion of DNA by PPi enzymatic treatment and a stepwise luminescence incorporation of nucleotide for the quantitative measurement of CpG islands methylation (31–33).
Several studies have reported certain tumor suppressor genes to be frequently methylated in epithelial malignancies, including those originating in the salivary glands. Among these genes are retinoid acid receptor β2 (RARβ2), RAS association domain family 1 (RASSF1A), O6-methylguanine-DNA methyltransferase (MGMT), and E-cadherin genes (7–14). All of which are known for their ability to suppress vital cellular processes, including cell-cycle regulation, apoptosis, DNA repair, differentiation, and metastasis (34–38). In a previous study, we used the MSP method to analyze the methylation status of some of these genes in these tumors (39). However, because of the arbitrary and qualitative nature of the MSP technique, the need for more objective and quantitative assessment of methylation has recently been raised (26).
In the present study, we evaluated by means of MSP and pyrosequencing the methylation status of the promoters for RARβ2, RASSF1A, MGMT, and E-cadherin in a large number of human SGC tissue specimens and in several SGC tumor cell lines and then compared the two datasets to determine which method might have greater potential clinical utility.
| Materials and Methods |
|---|
|
|
|---|
Tissue specimens. A total of 69 fresh-frozen tissue specimens from human SGCs harvested by one head and neck pathologic and stored in the Head and Neck Pathology Section at University of Texas M.D. Anderson Cancer Center, in accordance with institutional policies, formed the materials for this study. All of the patients whose specimens were used had presented to M.D. Anderson Cancer Center for surgical treatment of their tumors between 1994 and 2006. In each case, tumor and matching nonneoplastic salivary gland tissues were harvested separately. A frozen section evaluation of each tumor specimen was done at the acquisition and before storage. Frozen section evaluation was used as a quality control step and insured that all specimens contained an excess of 80% neoplastic cells. Harvested tissues were immediately frozen in liquid nitrogen and stored at –80°C until use. Each specimen was classified histopathologically according to WHO histologic classification system for SGC and was staged according to the tumor-node metastasis (TNM) classification system for malignant tumors defined by the International Union Against Cancer. To establish the baseline methylation status of targeted genes, 12 samples of normal salivary gland tissue were also analyzed.
DNA extraction and bisulfite treatment. Genomic DNA was extracted from fresh-frozen tissue specimens and cultured cells using DNAzol (Molecular Research Center, Inc., Cincinnati, OH) in accordance with the manufacturer's instructions. Human genomic DNA (5 µg) was artificially methylated as a positive control using SssI (CpG) methylase (New England Biolabs). Extracted DNA samples were modified by treatment with bisulfite using a CpGenome DNA modification kit (Intergen) in accordance with the manufacturer's instructions. In brief, 2 µg of genomic DNA were denatured by treatment with 2 mol/LNaOH and modified by treatment with 3 mol/L sodium bisulfate for 16 h. DNA samples were purified using Wizard DNA purification resin (Promega), treated with 3 mol/L NaOH, precipitated with ethanol, and resuspended in 20 µL water. Finally, 2-µL aliquots were used as templates for PCRs.
MSP methylation analysis. Because our materials consisted of fresh-frozen specimens, where the quality of extracted DNA is not a limiting factor, we used the conventional MSP method instead of the two-step nested MSP. The latter is more efficient in yielding better DNA quality from archived materials (40, 41). The PCR primer sequences of each gene and the annealing temperatures for pyrosequencing, MSP, and reverse transcriptase–PCR (RT-PCR) are listed (see electronic files).
Pyrosequencing. Pyrosequencing analysis was done using PCR assays designed to amplify specific CpG islands in the promoter of each target gene. Pyrosequencing primers were designed to amplify five or six target CpG dinucleotides and to not amplify CpGs within the primer sequence. Methylation of target CpGs was assessed by determining the ratio of cytosine to thymine incorporated during pyrosequencing; cytosine incorporation indicated a methylated CpG and thymine incorporation, an unmethylated CpG. This guaranteed that the ratio of incorporated cytosine to thymine would be proportional to the degree of methylation at a particular CpG site in the template DNA. Analysis of non-CpG cytosines provided an internal control for the completeness of bisulfite treatment. Hot-start PCRs were carried out in a volume of 50 µL containing Taq Gold polymerase (Applied Biosystems). Each reaction mixture contained 2 µg bisulfite-treated DNA, 0.25 mmol/L deoxynucleoside triphosphate, 2.5 mmol/L MgCl2, 5 µL PCR buffer, and 2 units Taq Gold polymerase. For analysis of RASSF1A and MGMT, 2 pmol of forward primer and 2 pmol of reverse primer were used; for analysis of RARβ2 and E-cadherin, 2 pmol of forward primer, 0.8 pmol of reverse primer, and 1.8 pmol of universal primer were used. The quality and purity of the PCR products were established by their electrophoresis through 2.5% agarose gels and subsequent ethidium bromide staining. Pyrosequencing was carried out using single-strand binding proteins (PyroGold reagents) and the PSQ96MA System (Biotage) in accordance with the manufacturer's protocol. All pyrosequencing experiments were done twice.
Methylation-specific PCR. Bisulfite-modified DNAs were amplified using primer sets specific for unmethylated and methylated sequences in each gene. Each PCR was done in a final volume of 25 µL containing template DNA, 1x PCR buffer, 2.5 mmol/L MgCl2, 0.25 mmol/L deoxynucleotide triphosphate, 2 µmol/L each primer, and 1 unit AmpliTaq Gold (Applied Biosystems). PCR amplification was done as follows: at 94°C for 10 min (one cycle); at 94°C for 45 s (35 cycles), at the specific annealing temperature for each gene for 45 s (one cycle), at 72°C for 1 min (one cycle); and, in a final extension step, at 72°C for 5 min. The resulting PCR products were electrophoresed on 2% agarose gels. DNA isolated from normal salivary gland tissue and modified by SssI methylase was used as a positive control. To assess methylation semiqualitatively, SssI-treated control DNA diluted stepwise (90-5%) with untreated DNA was used as a reference. On the basis of the dilution results, a distinct methylated band whose signal was at least 10% of the ratio of the signals for the methylated and unmethylated bands was identified and used as a cutoff point for scoring methylation (data not shown). The extent of methylation (i.e., methylation status) was recorded as the percentage of the relative ratio of the signals for the methylated and unmethylated bands (Fig. 1 ).
|
|
RNA preparation and RT-PCR analysis. RNA was isolated from the eight groups of cultured cells by treating the cells with TRIzol reagent in accordance with the manufacturer's instructions. RT-PCR analysis of RARβ2 mRNA expression was done under the following conditions: at 45°C for 45 min (one cycle); at 94°C for 2 min (one cycle); at 94°C for 30 s, 60°C for 1 min, and 68°C for 2 min (40 cycles); and, in a final extension step, at 68°C for 7 min. The resulting RT-PCR products were electrophoresed through 2.5% agarose gel, stained with ethidium bromide, and examined under UV light. GADPH mRNA expression in each RNA sample was also analyzed and used as an internal loading control.
Statistical analysis. To compare the results of the pyrosequencing-based and MSP-based methylation assessments, regression coefficients were calculated via linear regression methods. Scatter plots were drawn with least-squares regression lines bounded by 95% confidence limits. Cutoff values above which samples were coded as being hypermethylated were established by comparison with values from normal tissues and inspection of the scatter plots of values obtained from the two techniques and two-by-two tables of correlations between the two techniques at various cutoff points. Cutoff points for pyrosequencing and MSP, respectively, were as follows: RARβ2,
30% and
20%; RassF1A,
30% and
20%; MGMT,
15% and
10%; E-cadherin, pyrosequencing, >15%.
Correlations between methylation variables and clinicopathologic and end point variables were assessed by Pearson's
2 test or, where there were <10 subjects in any cell of a 2 x 2 grid, by the two-tailed Fisher exact test. Comparison of Kaplan-Meier product limit survival curves was done by application of the log-rank test. Statistical analyses were done using a commercially available software program (Statistical, version 7.1; StatSoft, Inc.).
| Results |
|---|
|
|
|---|
2 cm in 55 of the 66 cases (83%) for which it could be determined. Metastasis to regional cervical lymph nodes was noted in 21 of 69 patients (30%). The disease stage was I or II in 31 of 68 patients (46%) and III or IV in 37 (54%). The mean follow up time, from the date of presentation at our institution for this tumor until last contact or death was 56.4 months (SD, 41.8 months; median, 47.2 months with a range of 1.3-174.3 months). The overall 3-year survival rate was 65% (41 of 63). Six of the patients were lost to follow up before 3 years had elapsed. The 5-year survival rate was 46% (24 of 52).
Methylation findings. In normal salivary tissue specimens, pyrosequencing analysis revealed low levels of target gene methylation (i.e.,
10% in RASSF1A and MGMT and
15% in RARβ2 and E-cadherin). Because these were considered to represent background levels of methylation, they were used as cutoff points for analyzing methylation in the cell lines (data not shown). By both pyrosequencing and MSP techniques, only the RARβ2 promoter methylation was observed in all five cell lines (100%); in contrast, methylation of the RASSF1A, MGMT, and E-cadherin promoters was detected in only one cell line each. The pyrosequencing analysis of human tumor tissue specimens revealed a widely distributed variation in methylation levels among methylated tumors (n = 69; Fig. 3
), but the great majority of assays gave a result of <50% methylation (248 of 256, 96.9%). The most frequently methylated genes were RARβ2 (46%) and RASSF1 (35%). MGMT and E-cadherin were infrequently methylated. Twenty-seven of the patients (39%) showed no hypermethylation at any of the four sites. Methylation of multiple target gene promoters was relatively rare; only one promoter was methylated in 38% (26 of 69) of tumors, two promoters in 22% (15 of 69), three promoters in 1% (1 of 69), and all four in none (0%).
|
|
Correlation of methylation status and clinicopathologic factors. Table 2 presents the correlations between target gene methylation status the cutoff values listed (as determined by pyrosequencing and MSP) and clinicopathologic factors. Data for E-cadherin are not shown because only one patient had hypermethylation at that site. There were no statistically significant correlations between methylation status and sites, gender, or stage of their disease. There was a significant correlation between RARβ2 methylation status by pyrosequencing and tumor categories (P = 0.027; only one patient with ACC or adenoid cystic carcinoma had hypermethylation at RARβ2), tumor grade (P = 0.014), tumor size (P = 0.008), 3-year survival (P = 0.002), and 5-year survival (P = 0.028). There was a significant correlation between hypermethylation at RASSF1A, as measured by pyrosequencing, and diagnosis (P = 0.014). It is interesting to note that over half of patients with ACC had high methylation of this gene. Concurrent hypermethylation at RARβ2 and RASSF1A, correlated with diagnosis, high tumor grade, lymph node metastasis, and death within 3 years, was observed. The correlation with lower survival was confirmed by log-rank analysis of Kaplan-Meier survival curves (Fig. 4 ).
|
|
| Discussion |
|---|
|
|
|---|
Moreover, a faint shadow band may not be infrequently observed on gel electrophoresis of MSP analysis, and this may potentially lead to interpretive difficulties if only present in tumor specimen. We also encountered these faint bands in some cases on MSP analysis. In all of these instances, no methylation was detected in the corresponding pyrosequencing analysis. These artifacts, together with the semiquantitative and the subjective nature of the MSP scoring, provide further support for the application of the pyrosequencing. Although our results show that, at a certain level of methylation, MSP is a simple qualitative technique for screening and detecting hypermethylated genes in tumors, it is less sensitive at low level methylation and is prone to false-positive results. Although the two-step MSP method may enhance sensitivity, the method fundamentally shares many of the shortcomings of the MSP and may render similar results on fresh-frozen materials; further comparative studies are recommended (40, 41).
Our results also indicate that certain genes may be selectively methylated during salivary gland carcinogenesis. RARβ2 and RASSF1A, in contrast to MGMT and E-cadherin, were found to be frequently methylated in cell lines, as well as in human tumor tissue specimens, suggesting an association between the functional status of both genes and the development and/or progression of some SGCs. Contrary to previous reports, however, we found the E-cadherin promoter to be infrequently methylated in our tumor samples, a discrepancy that can be attributed to differences in the nature of the materials analyzed (archival versus fresh frozen), the interpretation of results, or both (41). In our study, we used fresh-frozen tissue, which yields better results than those obtained from archival materials (41). We also noted that concurrent methylation of RARβ2 and RASSF1A genes correlated significantly with certain clinicopathologic factors and patients survival, suggesting a synergistic effect of both genes on salivary gland tumorigenesis. Other investigators have implicated methylation at both RARβ-2 and RASSF-1 in the tumorigenesis of different human malignancies (13, 14, 34–39, 43–47).
In our study, RARβ2 was the only gene found to be highly methylated in all five cell lines and in most of the SGC tissue specimens tested. Moreover, regardless of the method of detection, RARβ2 methylation correlated significantly with the adenocarcinoma phenotype, tumor size, high tumor grade, and 3-year and 5-year survival, thus implicating this gene in the progression of certain types of SGC. A similar association between gene methylation and aggressive tumor phenotypes has been reported for other malignancies as well (38, 39, 43–47). Located on chromosome 3p24, the RARβ gene encodes multiple retinoid receptor isoforms (β1, β2, β4) by two different promoters and alternative splicing (38). The RARβ2 isoform has specifically been shown to be involved in carcinogenesis of several malignant tumors (40, 41, 43–50). Altered expression of retinoid receptors including RARβ2 has been reported in various cancer cell types, premalignant oral lesions (7–10) and head and neck squamous cell carcinoma (11–16). These alterations have in turn been linked to resistance to retinoids (7).
We also noted significant hypermethylation of the RASSF1-A in salivary duct and acinic cell carcinomas compared with other tumor types. Also, methylation of the same gene was found more in high grade, large size, and lymph node–positive tumors, but this did not reach a statistical significance. The RASSF-1 gene is located on chromosome 3p21.3 and encodes two major isoforms of a Ras-associated oncoprotein (51). RASSF-1A and RASSF-1C via alternative splicing and promoter usage. RASSF-1 protein has been shown to heterodimerize with other family members to induce Ras gene activation and tumorigenesis (43). Thus, suppression of RASSF-1 activity, whether dependent or not on the Ras pathway, seems to be a prerequisite for the induction of cell cycle arrest and apoptosis during tumorigenesis (50). Interestingly, our data revealed frequent concomitant RASSF1A and RARβ2 methylation in both ACCs and salivary duct carcinomas phenotypes. This finding suggests that the two genes may synergistically involve in the carcinogenesis of these two entities.
Our in vitro analysis of salivary cell lines showed that DAC-treatment led selectively to the reactivation of RARβ2 expression in only one of the five cell lines. However, the addition of Trichostatin either alone or in combination with DAC failed to show any appreciable changes. This finding indicates a minor role for histone deacetylation in the functional suppression of the RARβ2 gene in this cell line (42). Moreover, the inability of DAC to reactivate RARβ2 expression in the remaining cell lines suggests variable response of tumor cells to demethylating agents (38). Further studies with other demethylating agents are required to determine the differential response to these agents. Nonetheless, attempts are under way to delineate the functional consequences of reactivating these genes as a novel therapeutic target for patients with advanced disease (19, 20). The potential therapeutic application of demethylating agents in salivary carcinomas patients may await validation for prospective clinical trials.
In conclusion, the results of our present study (i.e., the first comprehensive and quantitative comparison of the gene methylation profiles of four major types of SGC) have revealed that the pyrosequencing method is more sensitive than the MSP method and its results correlate significantly with clinicopathologic factors and patients survival. The study also shows that the RARβ2 and RASSF1A genes are frequently methylated in ACC, adenoid cystic carcinoma, and salivary duct carcinoma and thus may warrant further study as potential therapeutic targets in that setting. Our findings also imply that future methylation-based investigations of SGC may provide useful information in the assessment and management of salivary gland tumors.
| 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.
Received 5/18/07; revised 10/26/07; accepted 12/27/07.
| References |
|---|
|
|
|---|
protein and hypermethylation of 14-3-3
gene in salivary gland adenoid cystic carcinoma. Br J Cancer 2004;91:1131–8.[Medline]This article has been cited by other articles:
![]() |
S. N. Vasilatos, G. Broadwater, W. T. Barry, J. C. Baker Jr., S. Lem, E. C. Dietze, G. R. Bean, A. D. Bryson, P. G. Pilie, V. Goldenberg, et al. CpG Island Tumor Suppressor Promoter Methylation in Non-BRCA-Associated Early Mammary Carcinogenesis Cancer Epidemiol. Biomarkers Prev., March 1, 2009; 18(3): 901 - 914. [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 |