Prostate cancer radiogenomics reveals proliferative gene expression programs associated with distinct MRI-based hypoxia levels

The biology behind individual hypoxia levels in patient tumors is poorly understood. Here, we used radiogenomics to identify associations between magnetic resonance imaging (MRI)-based hypoxia levels and biological processes derived from gene expression data in prostate cancer.For 85 prostate cancer patients, MRI-based hypoxia images were constructed by combining diffusion-weighted images reflecting oxygen consumption and supply. The ability to differentiate hypoxia levels in these images was verified by comparison with matched biopsy sections stained for the hypoxia marker pimonidazole. For MRI-defined hypoxia levels, corresponding hypoxic fractions were calculated and correlated with biopsy gene expression profiles. Biological processes were predicted by gene set enrichment analysis (GSEA) and validated by immunohistochemistry (Ki67 proliferation marker, reactive stroma grade) and RT-PCR (MYC).Genes with correlation between expression level and hypoxic fraction were identified for 56 MRI-based hypoxia levels. At all levels, GSEA identified proliferation as the predominant biological process enriched among the correlating genes. Two independent proliferative gene signatures were developed. The Peak1 signature, upregulated at moderate/severe hypoxia, reflected MYC upregulation and high Ki67-proliferation index of cancer cells in pimonidazole-positive regions. The Peak2 signature, upregulated at mild to non-hypoxic levels, was associated with fibroblast gene signature and reactive stroma grade. High scores of both Peak1 and Peak2 indicated elevated risk of biochemical recurrence in multiple cohorts.Radiogenomics identified two gene expression programs activated at different hypoxia levels, reflecting proliferation of cancer cells and stroma cells. Genes involved in these programs could be candidate targets for intervention.