Abstract

Intra-tumor heterogeneity is characterized by a diverse population of tumor clones and subclones which are important drivers of tumor evolution and therapeutic response. However, accurate subclonal reconstruction at scale remains challenging. We developed a machine learning tool, CliPP, and surveyed 9,972 tumors from 32 cancer types. We found that high subclonal mutation load (sML), the fraction of subclonal single nucleotide variants (SNVs) to all SNVs in the coding region, was prognostic of survival (progression free survival or overall survival) in 18 cancer types. In 14 cancers with low to moderate tumor mutation burden (TMB), high sML was associated with better prognosis. In immunotherapy trials for 42 metastatic prostate cancer (mCRPC), high sML was predictive of favorable response to ipilimumab and associated with increased CD8+ T-cell infiltration and decreased macrophage population. A validation using 613 whole-genomes of esophageal adenocarcinoma confirms the favorable effect of high sML and the observed tumor-associated macrophage. Our study identifies sML as a key feature of cancer, suggesting a biphasic relationship between evolutionary dynamics and differential immune environments. Finally, sML may serve as an orthogonal approach to identify likely responders of immune checkpoint blockade in low to moderate TMB tumors.