Abstract
The role of autophagy in colorectal cancer (CRC) pathogenesis appears to be crucial. Autophagy acts both as a tumor suppressor, by removing redundant cellular material, and a tumor-promoting factor, by providing access to components necessary for growth, metabolism, and proliferation. To date, little is known about the expression of genes that play a basal role in the autophagy in CRC. In this study, we aimed to compare the expression levels of 46 genes involved in the autophagy pathway between tumor-adjacent and tumor tissue, employing large RNA sequencing (RNA-seq) and microarray datasets. Additionally, we verified our results using data on 38 CRC cell lines. Gene set enrichment analysis revealed a significant deregulation of autophagy-related gene sets in CRC. The unsupervised clustering of tumors using the mRNA levels of autophagy-related genes revealed the existence of two major clusters: microsatellite instability (MSI)-enriched and -depleted. In cluster 1 (MSI-depleted), ATG9B and LAMP1 genes were the most prominently expressed, whereas cluster 2 (MSI-enriched) was characterized by DRAM1 upregulation. CRC cell lines were also clustered according to MSI-enriched/-depleted subgroups. The moderate deregulation of autophagy-related genes in cancer tissue, as compared to adjacent tissue, suggests a prominent field cancerization or early disruption of autophagy. Genes differentiating these clusters are promising candidates for CRC targeting therapy worthy of further investigation.