Analysis of NanoString data from validation sets using Non-Negative Matrix Factorization (NMF) analysis (26 genes up regulated in hypoxia) and dual clustering (samples and genes) further revealed that the increased level of BNIP3 (Bcl-2 adenovirus E1B 19 kDa-interacting protein 3)/GBE1 (glycogen branching enzyme1) differential pair correlates with the lack of response of melanoma patients to anti-PD1 (pembrolizumab) immunotherapy. Here, GBE1 is linked to melanoma.