Regulation of gene expression by DNA methylation with cytotoxic T lymphocytes evaluation in consensus molecular subtypes of colorectal cancer.

Abstract

25 Background: Low cytotoxic T lymphocyte (CTLs) infiltration in colorectal cancer (CRC) tumors is a challenge to treatment with immune checkpoint inhibitors. Consensus molecular subtypes (CMS) classify patients based on tumor attributes, and CMS1 patients include the majority of patients with high CTL infiltration and “inflamed” tumors. Epigenetic modification plays a critical role in gene expression and therapy resistance. Therefore, in this study we compared DNA methylation, gene expression, and CTL infiltration of CMS1 patients to other CMS groups to determine targets for improving immunotherapy in CRC. Methods: RNA-seq (n = 511) and DNA methylation (n = 316) from The Cancer Genome Atlas databases were used to determine gene expression and methylation profiles based on CMSs. CMS1 was used as a reference and compared to other subtypes (CMS2-4). Microenvironment Cell Populations- counter (MCPcounter) was used to determine tumor CTL infiltration. Genes with significantly different expression (p < 0.01, LogFC≥|1.5|) and difference of mean methylation β value ≥|0.25| were integrated for Pearson correlation coefficient analysis with MCPcounter score (r > |0.7|). Results: Comparing CMS1 and CMS2, ARHGAP9, TBX21, and LAG3 were differentially methylated and correlated with CTL scores. ARHGAP9 and TBX21 were decreased and hypomethylated in CMS2. Comparing CMS1 and CMS3, ARHGAP9, TBX21, FMNL1, HLA-DPB1, and STX11 were downregulated in CMS3 and highly correlated with CTL scores. ARHGAP9, FMNL1, HLA-DPB1, and STX11 were hypomethylated in CMS3 and TBX21 was methylated in both, but had a higher methylation ratio in CMS1. Comparing CMS1 and CMS4, TBX21 was the only gene downregulated, hypomethylated, and highly correlated with CTL scores in CMS4 patients. Conclusions: We found six genes differentially expressed, differentially methylated, and highly correlated with CTL infiltration when comparing CMS1 to other CMS groups. Specifically, TBX21 was the only gene highly correlated with CTL scores with differential gene expression and methylation in CMS2-4 when compared to CMS1. Thus, T-bet may be a critical regulator of T cell responses in CRC.