Clinico-genomic features predict distinct metastatic phenotypes in cutaneous melanoma

Metastasis drives mortality and morbidity in cancer. While some patients develop broad metastatic disease across multiple organs, others exhibit organ-specific spread. To identify mechanisms underlying metastatic organotropism, we analyzed clinico-genomic data from over 7,000 patients with metastatic cutaneous melanoma in three independent cohorts (one primary discovery and two validation cohorts including a nationwide electronic health record-derived deidentified database), leveraging machine learning approaches to clinical data. We found that female sex and increased tumor mutational burden associate with decreased metastatic potential, while older age associates with more lung and adrenal metastases. Using unsupervised analyses, patients clustered into five metastatic patterns: a "highly metastatic" cluster characterized by involvement of many organs, a "low metastatic" cluster characterized by few metastatic sites (mostly lymph node metastases), and three additional clusters each characterized by metastasis to specific sites (brain, lung, liver). Mutations in B2M and PTEN associated with increased overall metastatic potential. PTEN mutations were also associated with brain metastases but were enriched only in the "highly metastatic" cluster and not the brain-specific cluster. Mutations in GNAQ or GNA11 (GNA) associated with increased liver metastasis. To validate this association, we tested and demonstrated liver tropism in two GNA-mutant genetically engineered cutaneous melanoma mouse models of metastasis. Overall, our study elucidates distinct phenotypes of metastasis in patients with melanoma and identifies novel clinical and genomic associations that illuminate the clinical drivers of metastatic organotropism.