Prediction of survival of HPV16-negative, p16-negative oral cavity cancer patients using a 13-gene signature: A multicenter study using FFPE samples.

OBJECTIVES: To test the performance of an oral cancer prognostic 13-gene signature for the prediction of survival of patients diagnosed with HPV-negative and p16-negative oral cavity cancer. MATERIALS AND METHODS: Diagnostic formalin-fixed paraffin-embedded oral cavity cancer tumor samples were obtained from the Fred Hutchinson Cancer Research Center/University of Washington, University of Calgary, University of Michigan, University of Utah, and seven ARCAGE study centers coordinated by the International Agency of Research on Cancer. RNA from 638 Human Papillomavirus (HPV)-negative and p16-negative samples was analyzed for the 13 genes using a NanoString assay. Ridge-penalized Cox regressions were applied to samples randomly split into discovery and validation sets to build models and evaluate the performance of the 13-gene signature in predicting 2-year oral cavity cancer-specific survival overall and separately for patients with early and late stage disease. RESULTS: Among AJCC stage I/II patients, including the 13-gene signature in the model resulted in substantial improvement in the prediction of 2-year oral cavity cancer-specific survival. For models containing age and sex with and without the 13-gene signature score, the areas under the Receiver Operating Characteristic Curve (AUC) and partial AUC were 0.700 vs. 0.537 (p < 0.001), and 0.046 vs. 0.018 (p < 0.001), respectively. Improvement in predicting prognosis for AJCC stage III/IV disease also was observed, but to a lesser extent. CONCLUSIONS: If confirmed using tumor samples from a larger number of early stage oral cavity cancer patients, the 13-gene signature may inform personalized treatment of early stage HPV-negative and p16-negative oral cavity cancer patients.

Telomere Length and Lung Cancer Mortality among Heavy Smokers.

Background: Accumulating evidence suggests that short telomere length is associated with increased overall mortality, but the relationship with cancer mortality is less clear. We examined whether telomere length (global, and chromosome arm 5p- and 13q-specific) is associated with lung cancer mortality among cases from the ß-Carotene and Retinol Efficacy Trial of heavy smokers.Methods: Telomere length was measured on average 6 years before diagnosis for 788 lung cancer cases. Adjusted Cox proportional hazards models of all-cause and lung cancer-specific mortality were assessed for lung cancer overall and by histotype.Results: Short telomere length was associated with increased mortality for small cell lung cancer (SCLC), particularly stage III/IV SCLC [HR and 95% confidence interval for shortest vs. longest telomere length tertile: 3.32 (1.78-6.21)]. Associations were strongest for those randomized to the active intervention and when telomere length was measured ≤5 years before diagnosis. All-cause mortality patterns were similar. Short chromosome 5p telomere length was suggestively associated with lung cancer mortality, but there was no association with chromosome 13q telomere length.Conclusions: Our large prospective study suggests that among heavy smokers who developed lung cancer, short prediagnosis telomere length is associated with increased risk of death from SCLC.Impact: This is the first study to examine telomere length and mortality in lung cancer cases by histotype. If the association between short telomere length and SCLC mortality is replicated, elucidation of mechanisms through which telomere length influences survival for this highly aggressive cancer may inform more effective use of telomere-targeted therapeutics. Cancer Epidemiol Biomarkers Prev; 27(7); 829-37. ©2018 AACR.

Nested case-control study of telomere length and lung cancer risk among heavy smokers in the ß-Carotene and Retinol Efficacy Trial.

BACKGROUND: Telomeres protect cells from genomic instability. We examined telomere length and lung cancer risk prospectively in heavy smokers. METHODS: In a nested case-control study with 709 cases and 1313 controls, conditional logistic regression was used to evaluate associations between telomere length (global, chromosome 5p, and 13q) and lung cancer risk by histotype, controlling for detailed smoking history. RESULTS: Risks of overall lung cancer and adenocarcinoma were suggestively elevated among individuals with telomere length in the longest tertile. No clear patterns were observed for other histotypes, or for chromosome 5p or 13q telomere length. Associations with adenocarcinoma were strongest among (OR, 95% CI for longest versus shortest tertile): former smokers (2.26, 1.03-4.96), individuals <65 years (2.22, 1.13-4.35), and women (2.21, 0.99-4.93). CONCLUSIONS: Our large study of heavy smokers adds additional evidence that long telomere length prior to diagnosis is associated with risk of lung adenocarcinoma, but not other histotypes.

Transcriptome analysis reveals differentially expressed lncRNAs between oral squamous cell carcinoma and healthy oral mucosa.

Oral cavity and oropharyngeal squamous cell carcinoma (OSCC) is a major cancer type in the head and neck region. To better understand the roles long non-coding RNA (lncRNA) play in OSCC carcinogenesis, we compared the expression levels of 3,054 probe sets for lncRNAs between 167 OSCCs and 45 healthy oral mucosa using an Affymetrix HG U133 plus 2.0 array dataset. We found 658 lncRNA transcripts (790 probe sets) to be significantly differentially expressed using a criteria of FDR < 0.01, with 36 of them (39 probe sets) showing more than a 2-fold change. We further validated the top differentially expressed lncRNAs in three independent datasets from Gene Expression Omnibus (GEO) repository: GSE42743, GSE9844, and GSE6791. Fourteen lncRNAs (15 probe sets) were validated in all three datasets using the criteria FDR < 0.01: LOC441178, C5orf66-AS1, HCG22, FLG-AS1, CCL14/CCL15-CCL14, LOC100506990, TRIP10, PCBP1-AS1, LINC01315, LINC00478, COX10-AS1/LOC100506974, MLLT4-AS1, MIR31HG, and DUXAP10/LINC01296. Three lncRNAs in the validated list which showed the highest fold change (LOC441178, HCG22 and C5orf66-AS1) were verified by quantitative RT-PCR in a subset of 20 OSCCs and 10 control samples. In silico prediction of their functional role has given us directions for further investigation.

Genome-Wide Loss of Heterozygosity and DNA Copy Number Aberration in HPV-Negative Oral Squamous Cell Carcinoma and Their Associations with Disease-Specific Survival.

Oral squamous cell cancer of the oral cavity and oropharynx (OSCC) is associated with high case-fatality. For reasons that are largely unknown, patients with the same clinical and pathologic staging have heterogeneous response to treatment and different probability of recurrence and survival, with patients with Human Papillomavirus (HPV)-positive oropharyngeal tumors having the most favorable survival. To gain insight into the complexity of OSCC and to identify potential chromosomal changes that may be associated with OSCC mortality, we used Affymtrix 6.0 SNP arrays to examine paired DNA from peripheral blood and tumor cell populations isolated by laser capture microdissection to assess genome-wide loss of heterozygosity (LOH) and DNA copy number aberration (CNA) and their associations with risk factors, tumor characteristics, and oral cancer-specific mortality among 75 patients with HPV-negative OSCC. We found a highly heterogeneous and complex genomic landscape of HPV-negative tumors, and identified regions in 4q, 8p, 9p and 11q that seem to play an important role in oral cancer biology and survival from this disease. If confirmed, these findings could assist in designing personalized treatment or in the creation of models to predict survival in patients with HPV-negative OSCC.

Integrative genomics in combination with RNA interference identifies prognostic and functionally relevant gene targets for oral squamous cell carcinoma.

In oral squamous cell carcinoma (OSCC), metastasis to lymph nodes is associated with a 50% reduction in 5-year survival. To identify a metastatic gene set based on DNA copy number abnormalities (CNAs) of differentially expressed genes, we compared DNA and RNA of OSCC cells laser-microdissected from non-metastatic primary tumors (n = 17) with those from lymph node metastases (n = 20), using Affymetrix 250K Nsp single-nucleotide polymorphism (SNP) arrays and U133 Plus 2.0 arrays, respectively. With a false discovery rate (FDR)<5%, 1988 transcripts were found to be differentially expressed between primary and metastatic OSCC. Of these, 114 were found to have a significant correlation between DNA copy number and gene expression (FDR<0.01). Among these 114 correlated transcripts, the corresponding genomic regions of each of 95 transcripts had CNAs differences between primary and metastatic OSCC (FDR<0.01). Using an independent dataset of 133 patients, multivariable analysis showed that the OSCC-specific and overall mortality hazards ratio (HR) for patients carrying the 95-transcript signature were 4.75 (95% CI: 2.03-11.11) and 3.45 (95% CI: 1.84-6.50), respectively. To determine the degree by which these genes impact cell survival, we compared the growth of five OSCC cell lines before and after knockdown of over-amplified transcripts via a high-throughput siRNA-mediated screen. The expression-knockdown of 18 of the 26 genes tested showed a growth suppression ≥ 30% in at least one cell line (P<0.01). In particular, cell lines derived from late-stage OSCC were more sensitive to the knockdown of G3BP1 than cell lines derived from early-stage OSCC, and the growth suppression was likely caused by increase in apoptosis. Further investigation is warranted to examine the biological role of these genes in OSCC progression and their therapeutic potentials.

Integrative analysis in oral squamous cell carcinoma reveals DNA copy number-associated miRNAs dysregulating target genes.

OBJECTIVE: To better understand possible mechanisms involved in the dysregulation of gene expression unique to oral squamous cell carcinoma (OSCC) metastasis, the investigators examined the differential expression of microRNAs (miRNAs) in OSCC metastasis and their functional impact on target gene expression. STUDY DESIGN: Observational assessment of DNA copy number, miRNA, and RNA expression in primary and metastatic OSCC. SETTING: University of Washington Medical Center and affiliated hospitals. SUBJECTS: Tumor samples were taken from patients with primary incident OSCC; cells were laser-capture microdissected from 17 nonmetastatic primary tumors and 20 metastatic lymph nodes. METHODS: DNA copy number aberrations and gene expression profiles were previously determined using Affymetrix 250K Nsp I SNP arrays and HU133 plus 2.0 expression arrays. miRNAs were interrogated with Exiqon's Ready-to-Use PCR Panels assessing the expression of 368 human miRNAs. RESULTS: Investigators found 31 miRNAs differentially expressed between metastatic and nonmetastatic samples (false discovery rate <0.4; 26 overexpressed and 5 underexpressed in metastatic samples). Expression of 7 of these miRNAs was significantly associated with their DNA copy numbers, and expressions of 8 of these miRNAs were significantly associated with their target genes. Among these unique miRNAs, miR-140-3p, miR-29c, and miR-29a were differentially expressed in metastasis versus nonmetastatic samples and had a strong positive correlation with their DNA copy numbers and a negative correlation with the expression of their target genes. CONCLUSION: Results suggest that DNA copy number aberration may play a role in the dysregulation of some differentially expressed miRNAs in OSCC metastasis, warranting further investigation.

Tumor and salivary matrix metalloproteinase levels are strong diagnostic markers of oral squamous cell carcinoma.

BACKGROUND: The matrix metalloproteinases (MMP) cause degradation of the extracellular matrix and basement membranes, and thus may play a key role in cancer development. METHODS: In our search for biomarkers for oral squamous cell carcinomas (OSCC), we compared primary OSCC, oral dysplasia and control subjects with respect to: (i) expression of MMP1, MMP3, MMP10, and MMP12 in oral epithelial tissue using Affymetrix U133 2.0 Plus GeneChip arrays, followed by quantitative reverse transcription-PCR (qRT-PCR) for MMP1, and (ii) determination of MMP1 and MMP3 concentrations in saliva. RESULTS: MMP1 expression in primary OSCC (n = 119) was >200-fold higher (P = 7.16 × 10(-40)) compared with expression levels in nonneoplastic oral epithelium from controls (n = 35). qRT-PCR results on 30 cases and 22 controls confirmed this substantial differential expression. The exceptional discriminatory power to separate OSCC from controls was validated in two independent testing sets (AUC% = 100; 95% CI: 100-100 and AUC% = 98.4; 95% CI: 95.6-100). Salivary concentrations of MMP1 and MMP3 in OSCC patients (33 stage I/II, 26 stage III/IV) were 6.2 times (95% CI: 3.32-11.73) and 14.8 times (95% CI: 6.75-32.56) higher, respectively, than in controls, and displayed an increasing trend with higher stage disease. CONCLUSION: Tumor and salivary MMPs are robust diagnostic biomarkers of OSCC. IMPACT: The capacity of MMP gene expression to identify OSCC provides support for further investigation into MMPs as potential markers for OSCC development. Detection of MMP proteins in saliva in particular may provide a promising means to detect and monitor OSCC noninvasively.

Can a metastatic gene expression profile outperform tumor size as a predictor of occult lymph node metastasis in oral cancer patients?

PURPOSE: To determine the differential gene expression between oral squamous cell carcinoma (OSCC) with and without metastasis to cervical lymph nodes and to assess prediction of nodal metastasis by using molecular features. EXPERIMENTAL DESIGN: We used Affymetrix U133 2.0 plus arrays to compare the tumor genome-wide gene expression of 73 node-positive OSCCs with 40 node-negative OSCCs (≥ 18 months). Multivariate linear regression was used to estimate the association between gene expression and nodal metastasis. Stepwise logistic regression and receiver operating characteristics (ROC) analysis were used to generate predictive models and to compare these with models by using tumor size alone. RESULTS: We identified five genes differentially expressed between node-positive and node-negative OSCCs after adjusting for tumor size and human papillomavirus status: REEP1, RNF145, CTONG2002744, MYO5A, and FBXO32. Stepwise regression identified a four-gene model (MYO5A, RFN145, FBXO32, and CTONG2002744) as the most predictive of nodal metastasis. A leave-one-out ROC analysis revealed that our model had a higher area under the curve (AUC) for identifying occult nodal metastasis compared with that of a model by tumor size alone (respective AUC: 0.85 and 0.61; P = 0.011). A model combining tumor size and gene expression did not further improve the prediction of occult metastasis. Independent validation using 31 metastatic and 13 nonmetastatic cases revealed a significant underexpression of CTONG2002744 (P = 0.0004). CONCLUSIONS: These results suggest that our gene expression markers of OSCC metastasis hold promise for improving current clinical practice. Confirmation by others and functional studies of CTONG2002744 is warranted.

Integrative analysis of DNA copy number and gene expression in metastatic oral squamous cell carcinoma identifies genes associated with poor survival.

BACKGROUND: Lymphotropism in oral squamous cell carcinoma (OSCC) is one of the most important prognostic factors of 5-year survival. In an effort to identify genes that may be responsible for the initiation of OSCC lymphotropism, we examined DNA copy number gains and losses and corresponding gene expression changes from tumor cells in metastatic lymph nodes of patients with OSCC. RESULTS: We performed integrative analysis of DNA copy number alterations (CNA) and corresponding mRNA expression from OSCC cells isolated from metastatic lymph nodes of 20 patients using Affymetrix 250 K Nsp I SNP and U133 Plus 2.0 arrays, respectively. Overall, genome CNA accounted for expression changes in 31% of the transcripts studied. Genome region 11q13.2-11q13.3 shows the highest correlation between DNA CNA and expression. With a false discovery rate < 1%, 530 transcripts (461 genes) demonstrated a correlation between CNA and expression. Among these, we found two subsets that were significantly associated with OSCC (n = 122) when compared to controls, and with survival (n = 27), as tested using an independent dataset with genome-wide expression profiles for 148 primary OSCC and 45 normal oral mucosa. We fit Cox models to calculate a principal component analysis-derived risk-score for these two gene sets ('122-' or '27-transcript PC'). The models combining the 122- or 27-transcript PC with stage outperformed the model using stage alone in terms of the Area Under the Curve (AUC = 0.82 or 0.86 vs. 0.72, with p = 0.044 or 0.011, respectively). CONCLUSIONS: Genes exhibiting CNA-correlated expression may have biological impact on carcinogenesis and cancer progression in OSCC. Determination of copy number-associated transcripts associated with clinical outcomes in tumor cells with an aggressive phenotype (i.e., cells metastasized to the lymph nodes) can help prioritize candidate transcripts from high-throughput data for further studies.

The effect of acid suppression on upper airway anatomy and obstruction in patients with sleep apnea and gastroesophageal reflux disease.

STUDY OBJECTIVES: This study was designed to assess the effect of acid suppression on upper airway structure and function in patients with obstructive sleep apnea syndrome (OSAS) and gastroesophageal reflux disease (GERD). METHODS: This is a single-site within-subjects design. Twenty five patients with documented mild OSAS and objectively documented GERD via 24-hour pH measurement were included in the study. Patients were studied before and after 8 weeks of treatment with rabeprazole, 20 mg, twice a day. Subjects underwent laryngoscopy, polysomnography, and 24-hour pH monitoring. Subjective assessments of sleep obtained included the Pittsburg Sleep Quality Index and the Epworth Sleepiness Scale. RESULTS: Posterior commissure edema was significantly reduced (p < 0.05), and the Reflux Finding Score was improved (p < 0.07). Objective and subjective sleep parameters were significantly improved, sleep-onset latency was significantly reduced (26.2 vs 11.2, p < 0.05), and sleep-related acid contact time was significantly reduced (8.0% vs 1.7% p < 0.001). There was no significant change in the apnea-hypopnea index. CONCLUSIONS: In patients with mild OSAS and documented GERD, acid suppression improves upper airway abnormalities, as well as objective and subjective measures of sleep quality. Aggressive treatment of GERD in patients with OSAS may be helpful in the overall treatment of this select patient population.

A genetic expression profile associated with oral cancer identifies a group of patients at high risk of poor survival.

PURPOSE: To determine if gene expression signature of invasive oral squamous cell carcinoma (OSCC) can subclassify OSCC based on survival. EXPERIMENTAL DESIGN: We analyzed the expression of 131 genes in 119 OSCC, 35 normal, and 17 dysplastic mucosa to identify cluster-defined subgroups. Multivariate Cox regression was used to estimate the association between gene expression and survival. By stepwise Cox regression, the top predictive models of OSCC-specific survival were determined and compared by receiver operating characteristic analysis. RESULTS: The 3-year overall mean+/-SE survival for a cluster of 45 OSCC patients was 38.7+/-0.09% compared with 69.1+/-0.08% for the remaining patients. Multivariate analysis adjusted for age, sex, and stage showed that the 45 OSCC patient cluster had worse overall and OSCC-specific survival (hazard ratio, 3.31; 95% confidence interval, 1.66-6.58 and hazard ratio, 5.43; 95% confidence interval, 2.32-12.73, respectively). Stepwise Cox regression on the 131 probe sets revealed that a model with a term for LAMC2 (laminin gamma2) gene expression best identified patients with worst OSCC-specific survival. We fit a Cox model with a term for a principal component analysis-derived risk score marker and two other models that combined stage with either LAMC2 or PCA. The area under the curve for models combining stage with either LAMC2 or PCA was 0.80 or 0.82, respectively, compared with 0.70 for stage alone (P=0.013 and 0.008, respectively). CONCLUSIONS: Gene expression and stage combined predict survival of OSCC patients better than stage alone.

Simultaneous isolation of DNA and RNA from the same cell population obtained by laser capture microdissection for genome and transcriptome profiling.

Laser capture microdissection (LCM) is used extensively for genome and transcriptome profiling. Traditionally, however, DNA and RNA are purified from separate populations of LCM-harvested cells, limiting the strength of inferences about the relationship between gene expression and gene sequence variation. There have been no published protocols for the simultaneous isolation of DNA and RNA from the same cells that are obtained by LCM of patient tissue specimens. Here we report an adaptation of the Qiagen AllPrep method that allows the purification of DNA and RNA from the same LCM-harvested cells. We compared DNA and RNA purified by the QIAamp DNA Micro kit and the PicoPure RNA Isolation kit, respectively, from LCM-collected cells from adjacent tissue sections of the same specimen. The adapted method yields 90% of DNA and 38% of RNA compared with the individual methods. When tested with the GeneChip 250K Nsp Array, the concordance rate of the single nucleotide polymorphism heterozygosity calls was 98%. When tested with the GeneChip U133 Plus 2.0 Array, the correlation coefficient of the raw gene expression was 97%. Thus, we developed a method to obtain both DNA and RNA material from a single population of LCM-harvested cells and herein discuss the strengths and limitations of this methodology.

Genetic factors in catechol estrogen metabolism in relation to the risk of endometrial cancer.

2-Hydroxylated metabolites of estrogen have been shown to have antiangiogenic effects and inhibit tumor cell proliferation, whereas 4-hydroxylated metabolites have been implicated in carcinogenesis. We examined whether polymorphisms in certain genes involved in estrogen metabolism are associated with endometrial cancer risk in a population-based case-control study with 371 cases and 420 controls. Based on previously published genotype-phenotype correlation studies, we defined variant alleles thought to increase estrogen 2-hydroxylation as presumptively low-risk (CYP1A1 m1 T6235C and m2 Ile(462)Val) and those thought to increase estrogen 4-hydroxylation as high-risk (CYP1A1 m4 Thr(461)Asn, CYP1A2 A734C, and CYP1B1 Leu(432)Val). Odds ratios (OR) and 95% confidence intervals (95% CI) were calculated using unconditional logistic regression. Carrying at least one CYP1A1 m1 or m2 variant allele was associated with a decreased risk of endometrial cancer [ORs (95% CIs), 0.64 (0.44-0.93) and 0.54 (0.30-0.99), respectively]. No strong alteration in risk was observed among women with any of the putative high-risk alleles. When CYP1A1, CYP1A2, and CYP1B1 genotypes were combined and ranked by the number of putative low-risk genotypes carried, women with four or five low-risk genotypes had a reduced risk of endometrial cancer (OR, 0.29; 95% CI, 0.15-0.56) compared with women with one or none. No appreciable alteration in risk was observed among women carrying two or three low-risk genotypes. Some of our findings are consistent with the hypothesis that increased estrogen 2-hydroxylation is associated with decreased endometrial cancer risk, but replication of these results is required before any firm conclusions can be reached.