Abstract
Endothelin receptor type B (EDNRB) and kinesin family member 1A (KIF1A) are candidate tumor suppressor genes that are inactivated in cancers. In this study, we evaluated the promoter hypermethylation of EDNRB and KIF1A and their potential use for risk classification in prospectively collected salivary rinses from patients with premalignant/malignant oral cavity lesions. Quantitative methylation-specific PCR was performed to analyze the methylation status of EDNRB and KIF1A in salivary rinses of 191 patients. We proceeded to determine the association of methylation status with histologic diagnosis and estimate classification accuracy. On univariate analysis, diagnosis of dysplasia/cancer was associated with age and KIF1A or EDNRB methylation. Methylation of EDNRB highly correlated with that of KIF1A (P < 0.0001). On multivariable modeling, histologic diagnosis was independently associated with EDNRB (P = 0.0003) or KIF1A (P = 0.027) methylation. A subset of patients analyzed (n = 161) without prior biopsy-proven malignancy received clinical risk classification based on examination. On univariate analysis, EDNRB and risk classification were associated with diagnosis of dysplasia/cancer and remained significant on multivariate analysis (EDNRB: P = 0.047, risk classification: P = 0.008). Clinical risk classification identified dysplasia/cancer with a sensitivity of 71% and a specificity of 58%. The sensitivity of clinical risk classification combined with EDNRB methylation improved to 75%. EDNRB methylation in salivary rinses was independently associated with histologic diagnosis of premalignancy and malignancy and may have potential in classifying patients at risk for oral premalignant and malignant lesions in settings without access to a skilled dental practitioner. This may also potentially identify patients with premalignant and malignant lesions that do not meet the criteria for high clinical risk based on skilled dental examination. Cancer Prev Res; 3(9); 1093–103. ©2010 AACR.
Read the Perspective on this article by Lingen, p. 1056
Introduction
Head and neck squamous cell carcinoma (HNSCC) accounts for more than 37,000 new cases in the United States each year. Oral cavity cancers have an incidence of approximately 35,310 new cases per year and are the cause of 7,590 estimated deaths annually (1). Over the years, improvements have been made in the diagnosis, management, and targeted therapies for these cancers. However, despite these advances, a considerable number of patients continue to present with advanced stage disease. Advanced stage cancers have traditionally been associated with higher rates of mortality and decreased locoregional control rates. Intuitively, early detection of oral cancers would lead to improved quality of life and survival for the patients. The cost-effectiveness of screening for oral cavity cancers in high-risk patients has been previously shown (2).
The application of salivary rinses from high-risk patients has been explored as a potential for molecular screening in HNSCC (3–7). The inactivation of tumor suppressor genes caused by epigenetic changes such as promoter region CpG island hypermethylation has been well established in the literature (8–10). The use of real-time quantitative methylation-specific PCR (Q-MSP) provides a high-throughput mechanism for detecting promoter hypermethylation in patient samples. The ability to quantify the methylation through Q-MSP allows for the potential of identifying high-risk patients with premalignant lesions. This has been previously shown in salivary rinse samples obtained in lung cancer (11), oral cavity cancer (12), and oropharynx/hypopharynx cancer (3) patients.
We have previously published results of salivary rinse screening using promoter hypermethylation–based markers in patients with previously diagnosed HNSCC (12). To date, however, the effectiveness of this strategy, which includes salivary rinse samples collected in a prospective cohort who are at risk for oral cancer and oral premalignancy, has not yet been evaluated. In this study, we evaluated the utility of detection of methylation of two gene promoters, KIF1A and EDNRB, and the association of methylation of these promoter regions with the presence of oral cancer and premalignancy.
Materials and Methods
Tissue samples
Salivary rinse samples from 191 patients were prospectively obtained from the dental and oral medicine clinics at New York University College of Dentistry (140 patients, 73.30%), University of Puerto Rico (39 patients, 20.42%), St. Vincent's Cancer Center (1 patient, 0.52%), and Moffitt Cancer Center (11 patients, 5.76%). Institutional review board approval was obtained before the collection of the samples. A written informed consent was obtained from each subject. Enrollment included collection of demographic information and risk factor history (tobacco and alcohol). The inclusion criteria for enrollment were as follows: (a) English and/or Spanish speaking, (b) age >18 years, (c) the presence of a candidate oral epithelial lesion, and (d) the absence of a medical condition that would preclude a scalpel biopsy. Those lesions with an obvious etiology such as trauma, apthous ulceration, infection, or lichen planus were excluded. All patients were enrolled and assigned lesion risk classifications of low or high risk; this included an additional group of patients with histopathologically confirmed oral cancer who were also enrolled. Low-risk groups were defined as having leukoplakia without the presence of erythroplakia, ulceration, erosion, or submucosal extension/induration. High-risk patients exhibited one or more of the following: leukoplakia with ulceration, erosion, or submucosal/extension; erythroplakia; erythroleukoplakia; or ulceration. If a subject presented with more than one lesion, each lesion was separately classified and the overall risk classification was designated based on the worst lesion.
Salivary rinses were obtained from all subjects as previously described (12). In brief, salivary rinses were then obtained by rinsing and gargling with 25 mL of normal saline solution for 15 seconds. Three strokes were done using cotton-tipped applicators to collect exfoliated cells from the buccal mucosa, alveolar ridge, lateral tongue, floor of mouth, and pharyngeal inlet. The cellular materials from the applicators were also agitated and released into the salivary rinse specimen. This was then centrifuged to obtain a cell pellet after discarding the supernatant. Pellets were then immediately frozen and stored at −80°C. All low- and high-risk lesions underwent an incisional scalpel biopsy and a histologic diagnosis was obtained, which categorized the patients into one of six histologic diagnoses. Each specimen was examined by two calibrated oral pathologists blinded to the results of any other clinical data. Each specimen was examined at multiple levels and, in addition to a standard description and diagnosis, each pathologist recorded a possible primary histopathology outcome classification: benign (with or without atypia), dysplasia (mild, moderate, and severe grades), and invasive carcinoma (squamous cell carcinoma). Differences between oral pathologist classifications were resolved by joint review resulting in a consensus classification.
The tissue samples and clinical data were obtained from the New York University College of Dentistry in a collaborative effort by the Department of Otolaryngology-Head and Neck Surgery at Johns Hopkins Medical Institutions, Baltimore, Maryland. These patients constituted a subset of patients analyzed as part of a study previously conducted under a U54 mechanism at the New York University College of Dentistry (U54DE014257, Dr. David Sirois). The following procedures, including DNA extraction, bisulfite treatment, and Q-MSP, were all performed by two individuals blinded to the clinical data pertaining to the clinical risk classification and the histologic diagnoses.
DNA extraction
DNA obtained from the salivary rinse samples was extracted by the tissue bank by digestion with 50 μg/mL proteinase K (Boehringer) in the presence of 1% SDS at 48°C overnight followed by phenol/chloroform extraction and ethanol precipitation.
Bisulfite treatment
The DNA obtained from the salivary rinse samples was subjected to bisulfite treatment as has been previously described (13). Briefly, the EpiTect Bisulfite Kit was used to for the conversion of 2 μg of genomic DNA. The included Qiagen protocol was followed. After thermal denaturation and sodium bisulfite DNA conversion, the DNA was applied to an EpiTect spin plate. Optimized buffers and a vacuum manifold were used to wash and remove all traces of sodium bisulfate. The DNA was eluted. The eluted DNA was then ready for use for Q-MSP.
Quantitative methylation-specific PCR
The bisulfite-treated DNA was used as a template for fluorescence-based real-time Q-MSP as described previously (14). The EDNRB and KIF1A genes had been previously detected on a prior screen of salivary rinses in HNSCC patients.11
11S. Demokan, X. Chang, A. Chuang, W.K. Mydlarz, J. Kaur, P. Huang, Z. Khan, et al. KIF1A and EDNRB are differentially methylated in primary HNSCC and salivary rinses. Int J Cancer. Epub 2010 Feb 16.
Statistical analysis
The proportions of EDNRB or KIF1A gene methylation were compared between patient salivary rinse samples. The initial analyses were done with a cohort including 191 patients, where the subcategories of histologic outcome included benign, premalignant, and malignant. Two prespecified candidate genes (e.g., EDNRB and KIF1A) were evaluated. Gene hypermethylation was dichotomized at zero (i.e., no methylation versus any methylation). Predictors associated with head and neck cancers were also evaluated, including age, gender, race, smoking status, and alcohol consumption. Age was analyzed as a continuous variable, whereas all the other variables were considered as categorical variables. Univariate and multivariable proportional odds models were constructed sequentially to first explore the association of the variables of interest with the histologic outcome. Variables of significance based on the univariate models (P < 0.20) along with those deemed to be biologically/clinically important were retained for further analysis. Simultaneous effects expressed by these variables were studied using the multivariable proportional odds model. Odds ratios (OR) were reported with 95% confidence intervals (95% CI), which indicated the strength of the association and its uncertainty. Throughout the analyses, proportional odds assumption of common slope for all of the cumulative logits was checked by a score test.
A secondary analysis explored the association of methylation that is independent of other predictors of histology, in which patients with known cancer before coming to the clinic were excluded (n = 30). This subset of 161 patients (i.e., the patient cohort not including the known cancer patients) was categorized as benign versus dysplasia/cancer, exploring the association of clinical risk classification with histopathology. Univariate and multivariable logistic regression analyses were done using the same biologically/clinically important covariates as described above.
Receiver operating characteristic (ROC) analyses were conducted to estimate the classification accuracy, sensitivity (true-positive rate), and specificity (false-positive rate) of the predictor along with 95% CIs. The area under the curve (AUC), an index of predictive power, was also provided. A logistic prediction model using the clinical risk classification combined with methylation of EDNRB was developed and internally validated. Statistical analyses were done using SAS (v 9.2) and STATA software (v 8.2, SAS Institute); all statistical tests were two-sided with P < 0.05 considered statistically significant.
Results
Population characteristics
In analyzing the population characteristics of patients with oral cavity lesions (n = 191), we found that the majority of patients were males and Caucasian. There were 68.1%, 67.4%, and 74.3% males in the benign, dysplasia, and cancer categories, respectively. We also noted that there were 68.1%, 74.4%, and 68.6% Caucasians in the benign, dysplasia, and cancer categories, respectively. Median age was 54 years and ranged from 18 to 90 years in the entire cohort. The baseline characteristics of the groups were similar (Table 1). Furthermore, when assessing the exposure to clinically relevant risk factors such as tobacco and alcohol consumption, again we noted that the groups were inherently well matched. Tobacco consumption (current or past) was analyzed as a categorical variable and was noted in 67.3% of patients in the benign category and in 69.8% and 74.3% of patients in the dysplasia and cancer categories, respectively. Alcohol consumption was analyzed as a binary variable and was noted in 70.8%, 74.4%, and 77.1% of patients in the benign, dysplasia, and cancer categories, respectively. An exact χ2 test, however, revealed a statistically significant association between risk classification and histologic diagnosis in our 161-patient cohort (P < 0.008).
Univariate analysis of association between predictors and histology (n = 191)
Variable . | Benign (n = 113) . | Premalignant (n = 43) . | Malignant (n = 35) . | Total . |
---|---|---|---|---|
Age (y) | ||||
Mean (±SD) | 52 (±13) | 56 (±14) | 58 (±14) | 54 (±13) |
Median (range) | 52 (18-78) | 57 (26-82) | 56 (24-90) | 54 (18-90) |
Gender, n (%) | ||||
Female | 36 (31.9) | 14 (32.6) | 9 (25.7) | 59 (30.9) |
Male | 77 (68.1) | 29 (67.4) | 26 (74.3) | 132 (69.1) |
Race, n (%) | ||||
African American | 28 (24.8) | 9 (20.9) | 7 (20.0) | 44 (23.0) |
Caucasian | 77 (68.1) | 32 (74.4) | 24 (68.6) | 133 (69.6) |
Other | 8 (7.1) | 2 (4.7) | 4 (11.4) | 14 (7.3) |
Tobacco, n (%) | ||||
Never user | 37 (32.7) | 13 (30.2) | 9 (25.7) | 59 (30.9) |
Former user | 20 (17.7) | 7 (16.3) | 11 (31.4) | 38 (19.9) |
Current user | 56 (49.6) | 23 (53.5) | 15 (42.9) | 94 (49.2) |
Ethanol, n (%) | ||||
Never used | 33 (29.2) | 11 (25.6) | 8 (22.9) | 52 (27.2) |
Used | 80 (70.8) | 32 (74.4) | 27 (77.1) | 139 (72.8) |
Risk classification, n (%) | ||||
Low risk | 75 (66.4) | 21 (48.8) | 0 (0) | 96 (50.3) |
High risk | 38 (33.6) | 21 (48.8) | 6 (17.1) | 65 (34.0) |
Cancer | 0 (0) | 1 (2.3) | 29 (82.9) | 30 (15.7) |
EDNRB, n (%) | ||||
Unmethylated | 88 (77.9) | 26 (60.5) | 14 (40.0) | 128 (67.0) |
Methylated | 25 (22.1) | 17 (39.5) | 21 (60.0) | 63 (33.0) |
KIF1A, n (%) | ||||
Unmethylated | 84 (74.3) | 31 (72.1) | 17 (48.6) | 132 (69.1) |
Methylated | 28 (24.8) | 12 (27.9) | 18 (51.4) | 58 (30.4) |
Missing | 1 (0.9) | 0 (0) | 0 (0) | 1 (0.5) |
Variable . | Benign (n = 113) . | Premalignant (n = 43) . | Malignant (n = 35) . | Total . |
---|---|---|---|---|
Age (y) | ||||
Mean (±SD) | 52 (±13) | 56 (±14) | 58 (±14) | 54 (±13) |
Median (range) | 52 (18-78) | 57 (26-82) | 56 (24-90) | 54 (18-90) |
Gender, n (%) | ||||
Female | 36 (31.9) | 14 (32.6) | 9 (25.7) | 59 (30.9) |
Male | 77 (68.1) | 29 (67.4) | 26 (74.3) | 132 (69.1) |
Race, n (%) | ||||
African American | 28 (24.8) | 9 (20.9) | 7 (20.0) | 44 (23.0) |
Caucasian | 77 (68.1) | 32 (74.4) | 24 (68.6) | 133 (69.6) |
Other | 8 (7.1) | 2 (4.7) | 4 (11.4) | 14 (7.3) |
Tobacco, n (%) | ||||
Never user | 37 (32.7) | 13 (30.2) | 9 (25.7) | 59 (30.9) |
Former user | 20 (17.7) | 7 (16.3) | 11 (31.4) | 38 (19.9) |
Current user | 56 (49.6) | 23 (53.5) | 15 (42.9) | 94 (49.2) |
Ethanol, n (%) | ||||
Never used | 33 (29.2) | 11 (25.6) | 8 (22.9) | 52 (27.2) |
Used | 80 (70.8) | 32 (74.4) | 27 (77.1) | 139 (72.8) |
Risk classification, n (%) | ||||
Low risk | 75 (66.4) | 21 (48.8) | 0 (0) | 96 (50.3) |
High risk | 38 (33.6) | 21 (48.8) | 6 (17.1) | 65 (34.0) |
Cancer | 0 (0) | 1 (2.3) | 29 (82.9) | 30 (15.7) |
EDNRB, n (%) | ||||
Unmethylated | 88 (77.9) | 26 (60.5) | 14 (40.0) | 128 (67.0) |
Methylated | 25 (22.1) | 17 (39.5) | 21 (60.0) | 63 (33.0) |
KIF1A, n (%) | ||||
Unmethylated | 84 (74.3) | 31 (72.1) | 17 (48.6) | 132 (69.1) |
Methylated | 28 (24.8) | 12 (27.9) | 18 (51.4) | 58 (30.4) |
Missing | 1 (0.9) | 0 (0) | 0 (0) | 1 (0.5) |
The presence of methylated EDNRB promoter is associated with oral premalignancy and malignancy
A univariate analysis was then conducted to evaluate the association between histologic diagnosis and variables including age, gender, race, tobacco consumption, and alcohol consumption. In the analysis, we also included methylation status of EDNRB and KIF1A as potential predictors of histology (Table 2). Based on the univariate modeling, we observed that age was significantly associated with the diagnosis of benign, dysplasia, or cancer (OR, 1.3; 95% CI, 1.1-1.6; P = 0.014). Associations between histologic diagnosis and EDNRB and KIF1A methylation status were also found to be statistically significant [OR, 3.6 (95% CI, 2.0-6.4; P < 0.0001) and 2.2 (95% CI, 1.2-3.9; P = 0.011), respectively]. We did not observe a significant association with gender (P = 0.618) or race (P = 0.730) in the analysis. Although tobacco consumption (P = 0.372) and alcohol consumption (P = 0.435) did not suggest a significant association with histologic diagnosis in our cohort of patients, these variables have a well-established role in the progression to oral cancer.
Multivariable analyses of association between predictors and histology (n = 191)—EDNRB
Variable . | Univariate analysis . | Multivariable analysis . | ||
---|---|---|---|---|
OR* (95% CI) . | P* . | Adjusted OR† (95% CI) . | Adjusted P† . | |
Age (y) | 1.3‡ (1.1-1.6) | 0.014 | 1.2‡ (1.0-1.6) | 0.088 |
Gender | ||||
Female | Reference | Reference | ||
Male | 1.2 (0.6-2.1) | 0.618 | 0.9 (0.5-1.9) | 0.896 |
Race | ||||
African American | Reference | Reference | ||
Caucasian | 1.2 (0.6-2.5) | 0.730 | 1.1 (0.6-2.4) | 0.745 |
Other | 1.5 (0.5-4.9) | 1.6 (0.5-5.5) | ||
Tobacco | ||||
Never user | Reference | Reference | ||
Former user | 1.7 (0.8-3.8) | 0.372 | 1.7 (0.7-3.9) | 0.454 |
Current user | 1.1 (0.6-2.2) | 1.4 (0.7-3.0) | ||
Ethanol | ||||
Never used | Reference | Reference | ||
Used | 1.3 (0.7-2.4) | 0.435 | 1.2 (0.6-2.3) | 0.663 |
EDNRB | ||||
Unmethylated | Reference | Reference | ||
Methylated | 3.6 (2.0-6.4) | <0.0001 | 3.1 (1.7-5.8) | 0.0003 |
Variable . | Univariate analysis . | Multivariable analysis . | ||
---|---|---|---|---|
OR* (95% CI) . | P* . | Adjusted OR† (95% CI) . | Adjusted P† . | |
Age (y) | 1.3‡ (1.1-1.6) | 0.014 | 1.2‡ (1.0-1.6) | 0.088 |
Gender | ||||
Female | Reference | Reference | ||
Male | 1.2 (0.6-2.1) | 0.618 | 0.9 (0.5-1.9) | 0.896 |
Race | ||||
African American | Reference | Reference | ||
Caucasian | 1.2 (0.6-2.5) | 0.730 | 1.1 (0.6-2.4) | 0.745 |
Other | 1.5 (0.5-4.9) | 1.6 (0.5-5.5) | ||
Tobacco | ||||
Never user | Reference | Reference | ||
Former user | 1.7 (0.8-3.8) | 0.372 | 1.7 (0.7-3.9) | 0.454 |
Current user | 1.1 (0.6-2.2) | 1.4 (0.7-3.0) | ||
Ethanol | ||||
Never used | Reference | Reference | ||
Used | 1.3 (0.7-2.4) | 0.435 | 1.2 (0.6-2.3) | 0.663 |
EDNRB | ||||
Unmethylated | Reference | Reference | ||
Methylated | 3.6 (2.0-6.4) | <0.0001 | 3.1 (1.7-5.8) | 0.0003 |
*Univariate proportional odds model.
†Multivariate proportional odds model.
‡Unit of 10 y.
EDNRB and KIF1A were highly correlated (P < 0.0001) so that their individual effects on the observed histology could not be clearly separated if they were included simultaneously in the model. We then chose EDNRB along with other variables of significance suggested by the univariate model (e.g., age) as well as those deemed to play a clinically and biologically significant role in the development of oral cavity cancer (tobacco and alcohol consumption) to perform a multivariable proportional odds logistic regression analysis (Table 2). Once again, tobacco (P = 0.454) and alcohol (P = 0.663) consumption did not reach statistical significance. Age on multivariable analysis had an OR of 1.2 (1.0-1.6) but only reached borderline significance (P = 0.088). EDNRB methylation, however, remained significantly associated with histologic diagnosis (OR, 3.1; 95% CI, 1.7-5.8; P = 0.0003). These data indicate that EDNRB hypermethylation is an independent predictor of histologic diagnosis. Patients with methylated EDNRB were about three times as likely to have a diagnosis of premalignancy or malignancy as those with unmethylated EDNRB, after adjusting for age, tobacco use, and alcohol consumption.
EDNRB salivary rinse methylation analysis can improve risk classification when combined with clinical risk classification
A subset of our patient cohort initially presented with a known diagnosis of cancer. To obtain insight into the performance of clinical risk classification and salivary rinse methylation, we excluded those patients who had a diagnosis of known cancer on presentation (n = 30) from the cohort. The remaining patients (n = 161) were classified, using consensus risk classification methods based on the clinical appearance of the oral lesion, into high-risk (n = 21) and low-risk (n = 140) lesions (15).
We used the same predictors as used in our initial analysis and added risk classification designated during the initial examination to construct a univariate logistic regression model as shown in Table 3. Age, gender, race, tobacco and alcohol consumption, and KIF1A methylation status were not significantly associated with diagnosis. Risk classification, described as low risk or high risk, based on clinical examination was significantly associated with histologic diagnosis (OR, 2.5; 95% CI, 1.3-5.1; P = 0.008). EDNRB methylation status on univariate modeling was significant (OR, 2.1; 95% CI, 1.0-4.4; P = 0.046). We expanded the analysis to determine if there was any correlation between EDNRB and KIF1A. A Spearman correlation coefficient was calculated as R = 0.39 (P < 0.0001) and a scatter plot was constructed to better visualize if there appeared to be a correlation between the two markers (Fig. 1). The correlation between EDNRB and KIF1A was further evaluated by calculating the φ coefficient to measure the degree of correlation between two categorical variables, which was 0.225. Although these results may be limited because of the use of EDNRB and KIF1A as binary variables (methylated versus unmethylated) and the sample size, our findings suggested only a small relationship between EDNRB and KIF1A.
Univariate analysis of association between predictors and histology, excluding patients with known cancer at presentation (n = 161)
Variable . | Benign (n = 113) . | Dysplasia/cancer (n = 48) . | OR* (95% CI) . | P* . |
---|---|---|---|---|
Age (y) | ||||
Mean (±SD) | 52 (±13) | 56 (±15) | 1.2† (0.9-1.6) | 0.110 |
Median (range) | 52 (18-78) | 57 (24-86) | ||
Gender, n (%) | ||||
Female | 36 (31.9) | 15 (31.3) | Reference | |
Male | 77 (68.1) | 33 (68.7) | 1.0 (0.5-2.1) | 0.940 |
Race, n (%) | ||||
African American | 28 (24.8) | 10 (20.8) | Reference | |
White | 77 (68.1) | 36 (75.0) | 1.3 (0.6-3.0) | 0.639 |
Other | 8 (7.1) | 2 (4.2) | 0.7 (0.1-3.9) | |
Tobacco, n (%) | ||||
Never user | 37 (32.7) | 14 (29.2) | Reference | |
Former user | 20 (17.7) | 10 (20.8) | 1.3 (0.5-3.5) | 0.855 |
Current user | 56 (49.6) | 24 (50.0) | 1.1 (0.5-2.5) | |
Ethanol, n (%) | ||||
Never used | 33 (29.2) | 12 (25.0) | Reference | |
Used | 80 (70.8) | 36 (75.0) | 1.2 (0.6-2.7) | 0.587 |
Risk classification, n (%) | ||||
Low risk | 75 (66.4) | 21 (43.8) | Reference | |
High risk | 38 (33.6) | 27 (56.2) | 2.5 (1.3-5.1) | 0.008 |
EDNRB, n (%) | ||||
Unmethylated | 88 (77.9) | 30 (62.5) | Reference | |
Methylated | 25 (22.1) | 18 (37.5) | 2.1 (1.0-4.4) | 0.046 |
KIF1A, n (%) | ||||
Unmethylated | 84 (74.3) | 35 (72.9) | Reference | |
Methylated | 28 (24.8) | 13 (27.1) | 1.1 (0.5-2.4) | 0.782 |
Missing‡ | 1 (0.9) | 0 (0) |
Variable . | Benign (n = 113) . | Dysplasia/cancer (n = 48) . | OR* (95% CI) . | P* . |
---|---|---|---|---|
Age (y) | ||||
Mean (±SD) | 52 (±13) | 56 (±15) | 1.2† (0.9-1.6) | 0.110 |
Median (range) | 52 (18-78) | 57 (24-86) | ||
Gender, n (%) | ||||
Female | 36 (31.9) | 15 (31.3) | Reference | |
Male | 77 (68.1) | 33 (68.7) | 1.0 (0.5-2.1) | 0.940 |
Race, n (%) | ||||
African American | 28 (24.8) | 10 (20.8) | Reference | |
White | 77 (68.1) | 36 (75.0) | 1.3 (0.6-3.0) | 0.639 |
Other | 8 (7.1) | 2 (4.2) | 0.7 (0.1-3.9) | |
Tobacco, n (%) | ||||
Never user | 37 (32.7) | 14 (29.2) | Reference | |
Former user | 20 (17.7) | 10 (20.8) | 1.3 (0.5-3.5) | 0.855 |
Current user | 56 (49.6) | 24 (50.0) | 1.1 (0.5-2.5) | |
Ethanol, n (%) | ||||
Never used | 33 (29.2) | 12 (25.0) | Reference | |
Used | 80 (70.8) | 36 (75.0) | 1.2 (0.6-2.7) | 0.587 |
Risk classification, n (%) | ||||
Low risk | 75 (66.4) | 21 (43.8) | Reference | |
High risk | 38 (33.6) | 27 (56.2) | 2.5 (1.3-5.1) | 0.008 |
EDNRB, n (%) | ||||
Unmethylated | 88 (77.9) | 30 (62.5) | Reference | |
Methylated | 25 (22.1) | 18 (37.5) | 2.1 (1.0-4.4) | 0.046 |
KIF1A, n (%) | ||||
Unmethylated | 84 (74.3) | 35 (72.9) | Reference | |
Methylated | 28 (24.8) | 13 (27.1) | 1.1 (0.5-2.4) | 0.782 |
Missing‡ | 1 (0.9) | 0 (0) |
*Univariate logistic regression model.
†Unit of 10 y.
‡Excluded from the analysis of association.
Correlation between EDNRB and KIF1A. The above scatter plot provides a visual observation of correlation between numerical EDNRB and KIF1A (log transformed).
Correlation between EDNRB and KIF1A. The above scatter plot provides a visual observation of correlation between numerical EDNRB and KIF1A (log transformed).
We then proceeded to construct a multivariable model where we analyzed risk classification, age, gender, race, tobacco, ethanol, and EDNRB methylation status (Table 4). Once again, risk classification remained independently associated with histologic diagnosis (OR, 2.6; 95% CI, 1.3-5.2; P = 0.008). And similarly, EDNRB methylation showed statistical significance in being independently associated with diagnosis of dysplasia or cancer (OR, 2.1; 95% CI, 1.0-4.6; P = 0.047).
Risk factors for histologic diagnosis by multivariable analysis (n = 161)
Variable . | Adjusted OR* (95% CI) . | Adjusted P* . |
---|---|---|
Age (y) | 1.1† (0.8-1.5) | 0.455 |
Gender | ||
Male vs female | 0.9 (0.4-2.0) | 0.763 |
Race | ||
Caucasian vs African American | 1.3 (0.5-3.2) | 0.547 |
Other vs African American | 0.7 (0.1-4.0) | 0.659 |
Tobacco | ||
Former user vs never user | 1.5 (0.5-4.4) | 0.457 |
Current user vs never user | 1.6 (0.7-3.8) | 0.296 |
Ethanol | ||
Used vs never used | 1.3 (0.6-3.0) | 0.561 |
Risk classification | ||
High risk vs low risk | 2.6 (1.3-5.2) | 0.008 |
EDNRB | ||
Methylated vs unmethylated | 2.1 (1.0-4.6) | 0.047 |
Variable . | Adjusted OR* (95% CI) . | Adjusted P* . |
---|---|---|
Age (y) | 1.1† (0.8-1.5) | 0.455 |
Gender | ||
Male vs female | 0.9 (0.4-2.0) | 0.763 |
Race | ||
Caucasian vs African American | 1.3 (0.5-3.2) | 0.547 |
Other vs African American | 0.7 (0.1-4.0) | 0.659 |
Tobacco | ||
Former user vs never user | 1.5 (0.5-4.4) | 0.457 |
Current user vs never user | 1.6 (0.7-3.8) | 0.296 |
Ethanol | ||
Used vs never used | 1.3 (0.6-3.0) | 0.561 |
Risk classification | ||
High risk vs low risk | 2.6 (1.3-5.2) | 0.008 |
EDNRB | ||
Methylated vs unmethylated | 2.1 (1.0-4.6) | 0.047 |
*Multivariate logistic regression model.
†Unit of 10 y.
Through logistic regression modeling and ROC analyses, we then calculated the sensitivities and specificities of risk classification strategy and EDNRB methylation status (Table 5). Quantitative EDNRB was also considered to assess whether a cutoff other than zero may result in better predictive accuracy. The AUC and 95% CI were also calculated. Risk classification as a sole predictor of histology outcome had a sensitivity of 71% (95% CI, 56-83%) and a specificity of 58% (95% CI, 48-67%). The AUC for risk classification was 0.65 (95% CI, 0.56-0.075). EDNRB methylation in salivary rinse samples as a predictor of histologic diagnosis had a sensitivity of 65% (95% CI, 49-78%), a specificity of 51% (95% CI, 42-61%), and an AUC of 0.61 (95% CI, 0.51-0.71) when treated as a binary value. As described in Table 5, the ROC curves with adjustment for other covariates, illustrated in Fig. 2, show the AUC for EDNRB alone (0.61), risk classification alone (0.65), and the combination of risk classification and EDNRB (0.68).
Predictive accuracy of risk classification and EDNRB and their combination after adjusting for other predictors* associated with head and neck cancer (n = 161)
Predictor . | Cutoff† . | Sensitivity (%, 95% CI) . | Specificity (%, 95% CI) . | PPV‡ (%, 95% CI) . | NPV§ (%, 95% CI) . | AUC (95% CI) . |
---|---|---|---|---|---|---|
Risk classification | 0.267590 | 71 (56-83) | 58 (48-67) | 41 (31-53) | 77 (66-86) | 0.65 (0.56-0.75) |
EDNRB | 0.274459 | 65 (49-78) | 51 (42-61) | 36 (26-47) | 82 (72-90) | 0.61 (0.51-0.71) |
Risk and EDNRB combined | 0.246792 | 75 (60-86) | 50 (41-60) | 39 (29-50) | 83 (72-91) | 0.68 (0.58-0.77) |
0.450188 | 21 (10-35) | 92 (85-96) | 53 (29-76) | 73 (65-80) |
Predictor . | Cutoff† . | Sensitivity (%, 95% CI) . | Specificity (%, 95% CI) . | PPV‡ (%, 95% CI) . | NPV§ (%, 95% CI) . | AUC (95% CI) . |
---|---|---|---|---|---|---|
Risk classification | 0.267590 | 71 (56-83) | 58 (48-67) | 41 (31-53) | 77 (66-86) | 0.65 (0.56-0.75) |
EDNRB | 0.274459 | 65 (49-78) | 51 (42-61) | 36 (26-47) | 82 (72-90) | 0.61 (0.51-0.71) |
Risk and EDNRB combined | 0.246792 | 75 (60-86) | 50 (41-60) | 39 (29-50) | 83 (72-91) | 0.68 (0.58-0.77) |
0.450188 | 21 (10-35) | 92 (85-96) | 53 (29-76) | 73 (65-80) |
NOTE: The cutoffs are not EDNRB methylation values but a predicted probability from the logistic regression model that simultaneously includes risk classification and EDNRB.
*Other predictors include age, sex, race, tobacco, and ethanol use.
†Based on predicted probability of high-grade dysplasia/cancer using the multivariable logistic regression model.
‡Positive predictive value, depending on the prevalence of the disease (high-grade dysplasia/cancer), which was 13% for this study population.
§Negative predictive value, depending on the prevalence of the disease (high-grade dysplasia/cancer), which was 13% for this study population.
ROC curves with adjustment for other covariates as described in the legend to Table 5.
ROC curves with adjustment for other covariates as described in the legend to Table 5.
We also included the positive and negative predictive values in our analysis. Individuals with a positive test had a 41% (95% CI, 31-53%) chance of having high-grade dysplasia/cancer based solely on risk classification. The positive predictive value increased to 53% (95% CI, 29-76%) with the combination of risk classification and EDNRB when optimizing for specificity. We saw a similar increase in the negative predictive value with the combination of risk classification and EDNRB from 77% (95% CI, 66-86%) to 83% (95% CI, 72-91%) when optimizing for sensitivity.
We then performed logistic regression analysis and analyzed the variables by combining risk classification and EDNRB methylation to determine if there was any improvement in the predictive capability. Optimal cutoffs to maximize either sensitivity or specificity were obtained based on the predicted probability of high-grade dysplasia/cancer from the multivariable logistic regression model. Using the combination of risk classification and EDNRB at the selected cutoff to maximize specificity resulted in a specificity of 92% (95% CI, 85-96%) and a sensitivity of 21% (95% CI, 10-35%) with an AUC of 0.68 (95% CI, 0.58-0.77). At a cutoff threshold to maximize sensitivity, addition of EDNRB methylation improved sensitivity to 75% (95% CI, 60-86%) but decreased specificity to 50% (95% CI, 41-60%). In practical terms, the application of EDNRB salivary rinse methylation to define a high-risk category of patients within a cohort of 161 patients without prior biopsy would have changed the clinical low-risk assessment of 5 patients with mild dysplasia and 3 patients with severe dysplasia to high risk based on methylation analysis of salivary rinses.
Discussion
Aberrant promoter hypermethylation has been recently proposed as a means for detection of HNSCC in salivary rinses (3–7). We have studied a large cohort of prospectively collected salivary rinses obtained from patients with benign, dysplastic, and cancer diagnoses to determine the ability of Q-MSP to detect EDNRB or KIF1A promoter hypermethylation in high-risk patients. The sensitivity of the Q-MSP technique to detect the presence of EDNRB gene or KIF1A gene methylation enabled us to accurately correlate the risk classification strategy to the methylation status of the samples. This is the first report to show the use of molecular markers in salivary rinse samples for detection of premalignant oral disease.
Our group previously published the utility of evaluating the promoter region methylation status of various genes as a tool for detection of HNSCC (12). EDNRB hypermethylation has been studied extensively in prostate cancers with potential diagnostic as well as prognostic values (16–22). In addition, EDNRB methylation status has been studied in a variety of other cancers including lung cancer (23), bladder cancer (24, 25), hepatocellular carcinoma (26), nasopharyngeal carcinoma (27), and others (28, 29). More recently, our laboratory showed the presence of KIF1A promoter hypermethylation in breast cancer (30). We have also been successful in discovering that EDNRB and KIF1A are preferentially methylated in salivary rinse samples of HNSCC (see Supplementary Data 1)],12
12S. Demokan, X. Chang, A. Chuang, W.K. Mydlarz, J. Kaur, P. Huang, Z. Khan, et al. KIF1A and EDNRB are differentially methylated in primary HNSCC and salivary rinses. Int J Cancer 2009; submitted for publication.
We observed that the presence of EDNRB promoter methylation in salivary rinses was associated with the presence of dysplasia or invasive cancer, and that this was independent of clinical covariates including age and exposure history. This confirms that epigenetic alterations specific to dysplasia or invasive cancer can be detected in salivary rinses in the context of a dental clinic designed to assess patients at risk for oral cancer.
In our analysis, we found that risk classification by a specialist resulted in a 71% sensitivity in screening individuals with oral cavity lesions for dysplasia or cancer. This underscores the significance of clinical examination and risk classification as a gold standard of initial screening performed by highly trained individuals such as otolaryngologists/dentists. In addition, the combination of EDNRB with risk classification when using a selected cutoff threshold from the logistic regression analysis allowed an increase in sensitivity of 75% with a moderate decrease in specificity from 58% to 50%. In our analysis, we observed considerable variablility in the sensitivity and specificity based on selected cutoff values. Ideally, we would like to observe an increase in sensitivity without a substantial drop in specificity with the addition of EDNRB as a predictive variable. However, we noted that even in clinical situations where a skilled professional has designated a risk classification based on clinical examination, the addition of EDNRB methylation status may change the risk assessment of a lesion, prompting biopsy of a dysplastic or malignant oral lesion that would have otherwise not met the clinical criteria for biopsy. Our examined cohort has several demographic and exposure characteristics associated with oral cancer, including tobacco and ethanol exposure and advanced age. Therefore, the EDNRB promoter region methylation status in salivary rinse may be useful as a risk assessment tool in patients evaluated for potential oral malignancy based on exposure history and age but presenting with a clinically low-risk lesion.
The fundamental necessity of early detection in oral cancer can be confounded by numerous barriers. This can be due to the lack of a trained dentist/head and neck surgeon/otolaryngologist in the community as well as a basic lack of education/awareness in the public and health professionals (31, 32). In addition, lack of access to health care can also prevent patients from seeking care to facilitate earlier detection (33). The ideal test for oral premalignancy and cancer must be available for administration to a high-risk population and may be administered by health care workers without specialized training, yet still provide predictive outcome results. Q-MSP provides a cost-effective and easy-to-carry-out method that allows high-throughput and rapid analysis. This would indicate the potential use of this technique as a means for early detection of dysplastic oral cavity lesions and reinforces the potential usefulness of obtaining salivary rinses as a screening and surveillance strategy.
In a clinical setting, highly trained professionals examine patients with oral cavity lesions through conventional techniques of physical examination. Other established adjuncts include oral cytology, toluidine blue, and light-based detection systems (34). Based on history, clinical/environmental risk factors, and oral examination parameters, the clinician sets a threshold whereby those patients that are deemed to be at a higher risk of having oral cancers undergo the gold standard scalpel biopsy. This method often times results in a population of patients that are categorized as low-risk patients clinically, but may in fact harbor dysplastic lesions (35). Our study unveils the presence of false negatives by using a clinically based risk classification system alone by identifying patients who had a low-risk clinical lesion but had a histologic diagnosis of dysplasia. We established that the increase in sensitivity by combining both risk classification and EDNRB methylation status resulted from recognizing those clinically low-risk patients that in fact had dysplastic lesions. Due to the lack of long-term clinical outcomes for our patient cohort at this time, the true clinical implications of identifying more dysplastic lesions than would be discernable by a clinical risk stratification strategy alone need to be interpreted with caution. Although it is not apparent which patients with dysplastic lesions will proceed to develop cancer, we acknowledge that these patients most certainly warrant a more vigilant follow-up.
The use of molecular markers in salivary rinses for the detection of cancer or those harboring occult cancers has been explored with the intent to improve screening accuracy and cost-effectiveness. Salivary rinse samples potentially carry whole cells, DNA, RNA, and proteins, which allows for the capability of detecting alterations leading to cancer. Increasingly, saliva has been used to diagnose infectious diseases, hereditary disorders, autoimmune diseases, and endocrine disorders. Rosas et al. (6) published the first study demonstrating the detection of aberrant promoter hypermethylation in saliva from HNSCC patients. Carvalho et al. (12) identified differential hypermethylation patterns in salivary rinses and sera of patients with HNSCC in a panel of eight genes by Q-MSP. Similarly, Lallemant et al. (4) analyzed the expression levels of nine genes in HNSCC and control salivary rinse samples by quantitative reverse transcription-PCR. Zhao et al. (3) explored the feasibility of DNA PCR to screen for human papillomavirus in salivary rinse samples of head and neck cancer patients. The novelty of this article, however, is its prospective nature and the inclusion of salivary rinse samples from patients with premalignant lesions in the oral cavity.
Further investigations of the functionality of EDNRB and downstream pathways may yield additional insight into its role in oral cavity lesions. Future studies using CpG island microarray technologies may be useful in creating helpful panels of genes with increased sensitivities and retained specificities. Future studies to investigate the progression from premalignant changes to malignant transformation and the timing of the aberrant hypermethylation will also be of great value. Nonetheless, the use of salivary rinse molecular analysis may offer a feasible, rapid, and cost-effective tool for stratification of high-risk patients and early detection of premalignant lesions.
Disclosure of Potential Conflicts of Interest
D. Sidransky owns Oncomethylome Sciences, SA stock, which is subject to certain restrictions under University policy. D. Sidransky is a paid consultant to Oncomethylome Sciences, SA, and is a paid member of the company's Scientific Advisory Board. J.A. Califano is the Director of Research of the Milton J. Dance Head and Neck Endowment. The terms of this arrangement are being managed by the Johns Hopkins University in accordance with its conflict of interest policies.
Acknowledgments
Grant Support: National Institute of Dental and Craniofacial Research (NIDCR) and NIH Specialized Program of Research Excellence grant 5P50DE019032, NIDCR/NIH grant U54 DE14257, and Early Detection Research Network grant U01-CA084986. The funding agencies had no role in the design of the study, data collection or analysis, the interpretation of the results, the preparation of the manuscript, or the decision to submit the manuscript for publication.