Abstract
Melatonin is a pleiotropic anti-cancer molecule that controls cancer growth by multiple mechanisms. RNA-Seq can potentially evaluate therapeutic response and its use in xenograft tumor models can differentiate the changes that occur specifically in tumor cells or in the tumor microenvironment (TME). Melatonin actions were evaluated in a xenograft model of triple-negative breast cancer. Balb/c nude mice bearing MDA-MB-231 tumors were treated with melatonin or vehicle. RNA-Seq was performed on the Illumina HiSeq. 2500 and data were mapped against human and mouse genomes separately to differentiate species-specific expression. Differentially expressed (DE) genes were identified and Weighted Gene Co-expression Network Analysis (WGCNA) was used to detect clusters of highly co-expressed genes. Melatonin treatment reduced tumor growth (p < 0.01). 57 DE genes were identified in murine cells, which represented the TME, and were mainly involved in immune response. The WGCNA detected co-expressed genes in tumor cells and TME, which were related to the immune system among other biological processes. The upregulation of two genes (Tnfaip8l2 and Il1f6) by melatonin was validated in the TME, these genes play important roles in the immune system. Taken together, the transcriptomic data suggests that melatonin anti-tumor actions occur through modulation of TME in this xenograft tumor model.