Epithelial-mesenchymal plasticity (EMP) underlies embryonic development, wound healing, and cancer metastasis and fibrosis. Cancer cells exhibiting EMP often have more aggressive behavior, characterized by drug resistance, and tumor-initiating and immuno-evasive traits. Thus, the EMP status of cancer cells can be a critical indicator of patient prognosis. Here, we compare three distinct transcriptomic-based metrics—each derived using a different gene list and algorithm—that quantify the EMP spectrum. Our results for over 80 cancer-related RNA-seq datasets reveal a high degree of concordance among these metrics in quantifying the extent of EMP. Moreover, each metric, despite being trained on cancer expression profiles, recapitulates the expected changes in EMP scores for non-cancer contexts such as lung fibrosis and cellular reprogramming into induced pluripotent stem cells. Thus, we offer a scoring platform to quantify the extent of EMP in vitro and in vivo for diverse biological applications including cancer.

Epithelial to mesenchymal transition (EMT) is a well-studied hallmark of epithelial-like cancers that is characterized by loss of epithelial markers and gain of mesenchymal markers. Melanoma, which is derived from melanocytes of the skin, also undergo phenotypic plasticity toward mesenchymal-like phenotypes under the influence of various micro-environmental cues. Our study connects EMT to the phenomenon of de-differentiation (i.e., transition from proliferative to more invasive phenotypes) observed in melanoma cells during drug treatment. By analyzing 78 publicly available transcriptomic melanoma datasets, we found that de-differentiation in melanoma is accompanied by upregulation of mesenchymal genes, but not necessarily a concomitant loss of an epithelial program, suggesting a more “one-dimensional” EMT that leads to a hybrid epithelial/mesenchymal phenotype. Samples lying in the hybrid epithelial/mesenchymal phenotype also correspond to the intermediate phenotypes in melanoma along the proliferative-invasive axis - neural crest and transitory ones. As melanoma cells progress along the invasive axis, the mesenchymal signature does not increase monotonically. Instead, we observe a peak in mesenchymal scores followed by a decline, as cells further de-differentiate. This biphasic response recapitulates the dynamics of melanocyte development, suggesting close interactions among genes controlling differentiation and mesenchymal programs in melanocytes. Similar trends were noted for metabolic changes often associated with EMT in carcinomas in which progression along mesenchymal axis correlates with the downregulation of oxidative phosphorylation, while largely maintaining glycolytic capacity. Overall, these results provide an explanation for how EMT and de-differentiation axes overlap with respect to their transcriptional and metabolic programs in melanoma.

Epithelial-mesenchymal plasticity (EMP) underlies embryonic development, wound healing, and cancer metastasis and fibrosis. Cancer cells exhibiting EMP often have more aggressive behavior, characterized by drug resistance, and tumor-initiating and immuno-evasive traits. Thus, the EMP status of cancer cells can be a critical indicator of patient prognosis. Here, we compare three distinct transcriptomic-based metrics - each derived using a different gene list and algorithm - that quantify the EMP spectrum. Our results for 96 cancer-related RNA-seq datasets reveal a high degree of concordance among these metrics in quantifying the extent of EMP. Moreover, each metric, despite being trained on cancer expression profiles, recapitulates the expected changes in EMP scores for non-cancer contexts such as lung fibrosis and cellular reprogramming into induced pluripotent stem cells. Thus, we offer a scoring platform to quantify the extent of EMP in vitro and in vivo for diverse biological applications including cancer.

Epithelial to mesenchymal transition (EMT) is a well-studied hallmark of epithelial-like cancers that is characterized by loss of epithelial markers and gain of mesenchymal markers. Interestingly, melanoma, which is derived from melanocytes of the skin, also undergo phenotypic plasticity toward mesenchymal-like phenotypes under the influence of various micro-environmental cues. Our study connects EMT to the phenomenon of de-differentiation (i.e., transition from proliferative to more invasive phenotypes) observed in melanoma cells during drug treatment. By analyzing 78 publicly available transcriptomic melanoma datasets, we found that de-differentiation in melanoma is accompanied by upregulation of mesenchymal genes, but not necessarily a concomitant loss of an epithelial program, suggesting a more 'one-dimensional' EMT that leads to a hybrid epithelial/ mesenchymal phenotype. Samples lying in the hybrid epithelial/mesenchymal phenotype also correspond to the intermediate phenotypes in melanoma along the proliferative-invasive axis - neural crest and transitory ones. Interestingly, as melanoma cells progress along the invasive axis, the mesenchymal signature does not increase monotonically. Instead, we observe a peak in mesenchymal scores followed by a decline, as cells further de-differentiate. This biphasic response recapitulates the dynamics of melanocyte development, suggesting close interactions among genes controlling differentiation and mesenchymal programs in melanocytes. Similar trends were noted for metabolic changes often associated with EMT in carcinomas in which progression along mesenchymal axis correlates with the downregulation of oxidative phosphorylation, while largely maintaining glycolytic capacity. Overall, these results provide an explanation for how EMT and de-differentiation axes overlap with respect to their transcriptional and metabolic programs in melanoma.