Approximately 21% of patients with renal cell cancer (RCC) present with synchronous metastatic disease at the time of diagnosis, and metachronous metastatic disease occurs in 20-50% of cases within 5 years. Recent advances in adjuvant treatment of aggressive RCC following surgery suggest that biomarker-based prediction of risk for distant metastasis could improve patient selection. Biometrical analysis of TCGA-KIRC data identified candidate loci in the NK6 homeobox 2 gene (NKX6-2) that are hypermethylated in primary metastatic RCC. Analyses of NKX6-2 DNA methylation in three gene regions including a total of 16 CpG sites in 154 tumor-adjacent normal tissue, 189 RCC, and 194 metastatic tissue samples from 95 metastasized RCC patients revealed highly significant tumor-specific, primary metastatic-specific, and metastatic tissue-specific hypermethylation of NKX6-2. Combined CpG site methylation data for NKX6-2 and metastasis-associated genes (INA, NHLH2, and THBS4) demonstrated similarity between metastatic tissues and metastatic primary RCC tissues. The random forest method and evaluation of an unknown test cohort of tissues using receiver operator characteristic curve analysis revealed that metastatic tissues can be differentiated by a median area under the curve of 0.86 (p = 1.7 × 10-8-7.5 × 10-3) in 1000 random runs. Analysis of variable importance demonstrated an above median contribution for decision-making of at least one CpG site in each of the genes, suggesting superior informativity for sites annotated to NHLH2 and NKX6-2. Thus, DNA methylation of NKX6-2 is associated with the metastatic state of RCC tissues and contributes to a four-gene-based statistical predictor of tumoral and metastatic renal tissues.
Keywords: CpG-methylation; DNA-methylation; INA; NHLH2; NKX6-2; THBS4; hypermethylation; metastasis; prognosis; renal cell carcinoma; signature.