Spectral Clustering of Microarray Data Elucidates the Roles of Microenvironment Remodeling and Immune Responses in Survival of Head and Neck Squamous Cell Carcinoma

TY - JOUR

T1 - Spectral Clustering of Microarray Data Elucidates the Roles of Microenvironment Remodeling and Immune Responses in Survival of Head and Neck Squamous Cell Carcinoma

AU - Thurlow, Johanna K.

AU - Murillo, Claudia L. Pena

AU - Hunter, Keith D.

AU - Buffa, Francesca M.

AU - Patiar, Shalini

AU - Betts, Guy

AU - West, Catharine M. L.

AU - Harris, Adrian L.

AU - Parkinson, Eric K.

AU - Harrison, Paul R.

AU - Ozanne, Bradford W.

AU - Partridge, Max

AU - Kalna, Gabriela

PY - 2010/6/10

Y1 - 2010/6/10

N2 - Purpose To identify functionally related prognostic gene sets for head and neck squamous cell carcinoma (HNSCC) by unsupervised statistical analysis of microarray data. Patients and Methods Microarray analysis was performed on 14 normal oral epithelium and 71 HNSCCs from patients with outcome data. Spectral clustering (SC) analysis of the data set identified multiple vectors representing distinct aspects of gene expression heterogeneity between samples. Gene ontology (GO) analysis of vector gene lists identified gene sets significantly enriched within defined biologic pathways. The prognostic significance of these was established by Cox survival analysis. Results The most influential SC vectors were V2 and V3. V2 separated normal from tumor samples. GO analysis of V2 gene lists identified pathways with heterogeneous expression between HNSCCs, notably focal adhesion (FA)/extracellular matrix remodeling and cytokine-cytokine receptor (CR) interactions. Similar analysis of V3 gene lists identified further heterogeneity in CR pathways. V2CR genes represent an innate immune response, whereas high expression of V3CR genes represented an adaptive immune response that was not dependent on human papillomavirus status. Survival analysis demonstrated that the FA gene set was prognostic of poor outcome, whereas classification for adaptive immune response by the CR gene set was prognostic of good outcome. A combined FA&CR model dramatically exceeded the performance of current clinical classifiers (P <.001 in our cohort and, importantly, P = .007 in an independent cohort of 60 HNSCCs). Conclusion The application of SC and GO algorithms to HNSCC microarray data identified gene sets highly significant for predicting patient outcome. Further large-scale studies will establish the usefulness of these gene sets in the clinical management of HNSCC.

AB - Purpose To identify functionally related prognostic gene sets for head and neck squamous cell carcinoma (HNSCC) by unsupervised statistical analysis of microarray data. Patients and Methods Microarray analysis was performed on 14 normal oral epithelium and 71 HNSCCs from patients with outcome data. Spectral clustering (SC) analysis of the data set identified multiple vectors representing distinct aspects of gene expression heterogeneity between samples. Gene ontology (GO) analysis of vector gene lists identified gene sets significantly enriched within defined biologic pathways. The prognostic significance of these was established by Cox survival analysis. Results The most influential SC vectors were V2 and V3. V2 separated normal from tumor samples. GO analysis of V2 gene lists identified pathways with heterogeneous expression between HNSCCs, notably focal adhesion (FA)/extracellular matrix remodeling and cytokine-cytokine receptor (CR) interactions. Similar analysis of V3 gene lists identified further heterogeneity in CR pathways. V2CR genes represent an innate immune response, whereas high expression of V3CR genes represented an adaptive immune response that was not dependent on human papillomavirus status. Survival analysis demonstrated that the FA gene set was prognostic of poor outcome, whereas classification for adaptive immune response by the CR gene set was prognostic of good outcome. A combined FA&CR model dramatically exceeded the performance of current clinical classifiers (P <.001 in our cohort and, importantly, P = .007 in an independent cohort of 60 HNSCCs). Conclusion The application of SC and GO algorithms to HNSCC microarray data identified gene sets highly significant for predicting patient outcome. Further large-scale studies will establish the usefulness of these gene sets in the clinical management of HNSCC.

U2 - 10.1200/JCO.2009.24.8724

DO - 10.1200/JCO.2009.24.8724

M3 - Article

SN - 1527-7755

VL - 28

SP - 2881

EP - 2888

JO - Journal of Clinical Oncology

JF - Journal of Clinical Oncology

IS - 17

ER -