Metastatic colorectal adenocarcinoma tumor purity assessment from whole exome sequencing data

Tumors rich in stroma are associated with advanced stage and poor prognosis in colorectal adenocarcinoma (CRC). Abundance of stromal cells also has implications for genomic analysis of patient tumors as it may prevent detection of somatic mutations. As part of our efforts to interrogate stroma-cancer cell interactions and to identify actionable therapeutic targets in metastatic CRC, we aimed to determine the proportion of stroma embedded in hepatic CRC metastases by performing computational tumor purity analysis based on whole exome sequencing data (WES). Unlike previous studies focusing on histopathologically prescreened samples, we used an unbiased in-house collection of tumor specimens. WES from CRC liver metastasis samples were utilized to evaluate stromal content and to assess the performance of three in silico tumor purity tools, ABSOLUTE, Sequenza and PureCN. Matching tumor derived organoids were analyzed as a high purity control as they are enriched in cancer cells. Computational purity estimates were compared to those from a histopathological assessment conducted by a board-certified pathologist. According to all computational methods, metastatic specimens had a median tumor purity of 30% whereas the organoids were enriched for cancer cells with a median purity estimate of 94%. In line with this, variant allele frequencies (VAFs) of oncogenes and tumor suppressor genes were undetectable or low in most patient tumors, but higher in matching organoid cultures. Positive correlation was observed between VAFs and in silico tumor purity estimates. Sequenza and PureCN produced concordant results whereas ABSOLUTE yielded lower purity estimates for all samples. Our data shows that unbiased sample selection combined with molecular, computational, and histopathological tumor purity assessment is critical to determine the level of stroma embedded in metastatic colorectal adenocarcinoma.