Abstract
Colorectal cancer (CRC) is the third-leading cause of cancer-related deaths in the United States and worldwide. Obesity-a worldwide public health concern-is a known risk factor for cancer including CRC. However, the mechanisms underlying the link between CRC and obesity have yet to be fully elucidated in part because of the molecular heterogeneity of CRC. We hypothesized that obesity modulates CRC in a consensus molecular subtype (CMS)-dependent manner. RNA-seq data and associated tumor and patient characteristics including body weight and height data for 232 patients were obtained from The Cancer Genomic Atlas-Colon Adenocarcinoma (TCGA-COAD) database. Tumor samples were classified into the four CMSs with the CMScaller R package; body mass index (BMI) was calculated and categorized as normal, overweight, and obese. We observed a significant difference in CMS categorization between BMI categories. Differentially expressed genes (DEGs) between obese and overweight samples and normal samples differed across the CMSs, and associated prognostic analyses indicated that the DEGs had differing associations on survival. Using Gene Set Enrichment Analysis, we found differences in Hallmark gene set enrichment between obese and overweight samples and normal samples across the CMSs. We constructed Protein-Protein Interaction networks and observed differences in obesity-regulated hub genes for each CMS. Finally, we analyzed and found differences in predicted drug sensitivity between obese and overweight samples and normal samples across the CMSs. Our findings support that obesity impacts the CRC tumor transcriptome in a CMS-specific manner. The possible associations reported here are preliminary and will require validation using in vitro and animal models to examine the CMS-dependence of the genes and pathways. Once validated the obesity-linked genes and pathways may represent new therapeutic targets to treat colon cancer in a CMS-dependent manner.