Please cite this paper as:

Chuffa, L., Carvalho, R., Camargo, V.L., Cury, S., Domeniconi, R., Zuccari, D.A. and Seiva, F.R. 2023. Melatonin and cancer: Exploring gene networks and functional categories. Melatonin Research. 6, 4 (Dec. 2023), 431-451. DOI:https://doi.org/https://doi.org/10.32794/mr112500161.


Research Article

Melatonin and cancer: Exploring gene networks and functional categories

Luiz Gustavo de Almeida Chuffa1*, Robson Francisco Carvalho1, Victória Larissa Schimidt Camargo1, Sarah Santiloni Cury1, Raquel Fantin Domeniconi1, Debora Aparecida Pires de Campos Zuccari2, Fábio Rodrigues Ferreira Seiva3

1Department of Structural and Functional Biology, UNESP - São Paulo State University, Institute of Biosciences, Botucatu, 18618–689, São Paulo, Brazil.

2Faculdade de Medicina de São José do Rio Preto, São José do Rio Preto, SP, 15090–000, Brazil.

3Department of Chemical and Biological Sciences, São Paulo State University (UNESP), Institute of Bioscience, Botucatu, 18618–689, São Paulo, Brazil.

*Correspondence: guchuffa@yahoo.com.br, Tel: +55 (14) 3880-0027

Running title: Gene networks in melatonin-treated cancer 

Received: June 2, 2023; Accepted: December 7, 2023

 

ABSTRACT

     While melatonin is known for its multifaceted properties and its potential to combat cancer, there has been limited exploration of the cancer-melatonin interaction at the gene network level. One of the ways to better understand the molecular mechanisms of melatonin’s anti-cancer effects is to use text-mining strategies to extract relevant information that creates knowledge networks of entities and their associations. In this study, we mined gene-publication associations to search for genes most relevant to the terms of “melatonin” and “cancer”. A total of 152 genes were identified and ranked, among which 15 were kinase-related and three G-protein coupled receptor genes. The hub genes (STAT3JUNTP53MAPK3EP300SRCHSP90AA1AKT1ESR1, and IL6) were involved with several pathways in cancer. After examining the melatonin-treated cancers, we mapped 25 upregulated and 51 downregulated genes; these were strongly associated with cancer hallmarks such as resisting cell death, sustaining proliferative signaling, and inducing invasion and metastasis. Upregulated genes showed molecular functions including apoptotic protease activator, caspase activator, enzyme regulator, and protein binding, whereas the downregulated genes affected protein kinase activities, transcription factor binding, protein, enzyme, DNA, and promoter bindings. By connecting gene subsets, we detected a closer relationship among breast, hepatocellular, prostate, and oral cancers, in addition to neuroblastoma and osteosarcoma in terms of changes in melatonin-related signaling pathways. TCGA data were analyzed to understand the impact of gene signatures on survival of patients, and melatonin-downregulated genes were associated with longer survival of patients with glioblastoma, bladder, breast, colon, stomach, liver, lung, and ovarian carcinomas. These results provide a global view of gene interaction networks in melatonin-treated cancers and their functional value, opening new opportunities to consider melatonin for cancer therapy.

Key words: Melatonin, cancer, target genes, regulatory network, functional enrichment, gene mining

___________________________________________________________________________


1.      INTRODUCTION

     Cancer is an enormous public health problem, with an estimated 19.3 million new cases and 10.0 million cancer deaths in 2020 globally (1). Based on demographic projections, global cancer burden is expected to reach 28.4 million cases in 2040. Cancer ranks as the most important barrier to increasing life expectancy worldwide in the 21st century (2). The multifactorial etiology, genetic factors, and lifestyle conditions may further exacerbate the epidemiology of cancers (3). There are more than 150 different types of cancer and the lack of effective early stage diagnosis hinders the cure of patients (4, 5). Cancer is considered a disease of the genes since various cellular dynamics rely on specific gene activation or repression (5). These cells are usually accompanied by evasion of apoptosis, uncontrolled cellular proliferation, cellular migration, invasion, and metastasis, and changes in energy metabolism (6). The cancer-related signaling networks are heterogeneous and provide informative data for dynamics given that the networks contain both directed and undirected interaction types (7). The identification of functional agents capable of regulating a wide range of molecules in the pathological cells would be an important therapeutic ally.    

     Melatonin (N-acetyl-5-methoxytryptamine), a highly pleiotropic molecule, is an anti-tumor regulator synthesized in a circadian and in a non-circadian manner (8). At the gene expression level, melatonin acts directly or indirectly on RNAs and proteins, including microRNAs, long-non-coding RNAs, Piwi-interacting RNA, DNA methyltransferases, and histone-modifying enzymes, in diverse conditions such as in cancer (8, 9). Using a meta-analysis approach, we have reported that melatonin regulates various microRNAs in breast, gastric, oral, colorectal, and prostate cancers, and glioma, and these microRNAs were potential regulators of target genes implicated in cell proliferation, differentiation, apoptosis, senescence, and autophagy (10).   

     Although the relationship between melatonin and cancer has been known for decades (11), the knowledge of how melatonin modulates molecular networks in different cancer types remains sparingly investigated. The rather slow progress in this area has been due to the numerous convoluted means of action of melatonin in cancer inhibition; because of this, the indolamine has been described as a selective smart killer of cancer cells (12). Melatonin has well-known actions on cancer cells such as pro-oxidative effects, regulation of estrogenic receptor expression, depletion of telomerase activity, epigenetic alterations, and disruption of energy metabolism (13-15). Additional information on the close crosstalk between genetic and epigenetic regulation of cancer by melatonin would be helpful to clarify proper molecular mechanisms that qualify the use of this molecule, along with other agents, for clinical application.

     It is well documented that melatonin exerts pro-apoptotic, immunomodulatory, antiproliferative, and anti-angiogenic actions (16-18). Melatonin has long been considered a multifaceted molecule interfering with cancer hallmarks, and epidemiological studies have documented shiftwork and light-at-night as contributors associated with increased incidence of cancers (19-21). A particular interest is the fact that melatonin levels are reduced in cancer cells (22, 23), prompting us to examine whether its synthetic rate in mitochondria of cancer cells are impaired or whether the circulating melatonin is not being sufficiently taken up by these cells; the prediction of melatonin synthesis using gene expression data of cancer patients has revealed a significant decline of melatonin levels (24-26).

     Considering that melatonin has functional effects on regulatory RNAs (27), we mapped gene interactions in multiple cancers and uncovered relevant gene sets that are regulated by melatonin in each cancer type. We also showed the connections of these genes with cancer hallmarks, functional enrichment, and their clinical significance.

 

2. MATERIALS AND METHODS

2.1. Identification of potential genes involved with melatonin and cancer.

     To comprehensively identify the genes related to the terms of “cancer” and “melatonin” based on published literature, we initially used Geneshot software (28) to identify and filter the relevant genes. For better data refinement, we searched for studies using the following terms: cancer and melatonin, cancers and melatonin, tumor and melatonin, tumors and melatonin, melatonin and neoplasia, and melatonin and neoplasms. After uncovering these genes, they were ranked according to the terms into publications over the time so that any correlation between cancer and melatonin was identified. Next, we manually curated the data identified in the studies through the PubMed electronic database (http://www.ncbi.nlm.nih.gov/pubmed), linking the cancer type with the appropriate response to melatonin treatment. The inclusion criteria were: studies validated in cancer samples (in vivo or in vitro), studies with untreated tissues for comparison (control groups), and data providing statistical significance. Studies involving non-cancerous tissues/cells, reviews, studies that were not case-controlled and without complete access, and studies involving in silico analyses were excluded.

2.2. Protein-protein interaction (PPI) networks based on melatonin and cancer-associated genes.

     Functional molecular network was generated from those genes associated with melatonin and overall cancers. The metasearch STRING database (Search Tool for Retrieval of Interacting Genes, v. 11.5) was used for creating PPI enrichment (29). The following interaction sources were set: database, text mining, experiments, and co-occurrence. We also applied the minimum interaction score of 0.900 (highest confidence), and all disconnected nodes in the network were hidden to facilitate visual analysis (P-values represented the statistical significance provided by STRING). Access in January 2023. The complex gene network interaction was used for annotation and visualization using the open-source Cytoscape platform v.3.9.1 (30) available at https://cytoscape.org/; a total of 152 nodes were integrated into a regulatory network consisting of 713 interactions (average node degree = 9.38, average clustering coefficient = 0.507). Using cytoHubba (31), the entire network was explored for important nodes and subnetworks through several topological algorithms to detect the top hub genes. These included degree, edge percolated component, maximum neighborhood component, density of maximum neighborhood component, maximal centrality and centralities based on shortest paths, such as bottleneck, EcCentricity, closeness, radiality, betweenness, and stress.

2.3. Identification of genes across melatonin-treated cancers.

     We initially examined all the genes together regardless of cancer type and discriminated between those genes, which were up and downregulated by melatonin based on each cancer modality. The CancerGeneNet (32) was used to link genes that are frequently mutated in cancer-related phenotypes. This software bridges two interaction layers by connecting molecules whose activities are affected by cancer gene products to proteins that impact cancer phenotypes; this is achieved by implementing graph algorithms linking any gene of interest to the hallmarks of cancer (the score of 0.9 was defined). Next, we used melatonin-related genes to identify different ontology categories. Biological processes (Gene Ontology) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were analyzed using the EnrichR tool (available at http://amp.pharm.mssm.edu/Enrichr/), and the top 10 enriched terms were generated according to the lowest p-value < 0.05 (Fisher exact test) provided by the EnrichR database (33, 34). To examine gene interactions that were over and under-expressed in melatonin-treated cancers, we used WebGestalt (35), available at https://www.webgestalt.org/, to translate the gene list into biological insights. First, the networks were mapped under topology-based analysis using PPI Biogrid (FDR < 0.05; Benjamini-Hochberg adjusted p-value). Next, over-representation analysis of the genes was applied to find different categories (biological process, cellular component, and molecular function).

2.4. Survival analysis.

     Based on the representative individual gene set affected by melatonin, we used GEPIA2 (http://gepia2.cancer-pku.cn/#survival) (36) to determine the risk assessment and to perform overall survival of patients according to TCGA datasets [bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), stomach adenocarcinoma (STAD), glioblastoma multiforme (GBM), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), ovarian serous cystadenocarcinoma (OV)]. Using different gene signatures, we assessed individual correlations with overall survival by Cox Proportional Hazard regression, which stratified the risk groups (high-risk and low-risk), estimated by an optimization algorithm. The prognostic value of each signature predicting survival was visualized as the Kaplan–Meier method and log-rank test. Through the survival maps, we examined the survival contribution of the gene signatures influenced by melatonin in multiple cancer types using the Mantel-Cox test (FDR < 0.01). 

2.5. Single-cell RNA-sequencing (scRNA-seq) analysis.

     We re-analyzed publicly available scRNA-seq data to determine the transcriptional profile of 15 melatonin-related genes which showed involvement with patients’ survival. Next, BRCA (GSE148673, GSE161529, and E-MTAB-8107), LUAD (GSE131907, GSE123904, HRA000154, and E-MTAB-6149), OV (E-MTAB-8107), GBM (EGAS00001002185, EGAS00001001900, and EGAS00001003845), and COAD (EGAS00001003779) scRNA-seq data were downloaded from Curated Cancer Cell Atlas (https://www.weizmann.ac.il/sites/3CA/). We selected untreated primary samples belonging to the 10X platform with count data available. The selected datasets and samples are described in the Supplementary Table 3. LIHC, BLCA, and STAD were excluded considering the lack of data availability and our inclusion criteria. The count data were processed using Seurat (v4.0.6) package (37) in RStudio software. BRCA and LUAD had more than one available dataset and were integrated using the IntegrateData function. The identities of cell types were based on expression of canonical markers.

2.6. Data representation and analysis.

     To provide a reliable representation of the molecular functions into a gene set based on the GO hierarchy (FDR<0.01), we used the Biological Networks Gene Ontology tool (BiNGO) (38), visualized in Cytoscape. A hypergeometric test with False Discovery Rate (FDR) correction was applied (significance level < 0.05) when searching for overrepresented categories. Bar plots were constructed using WebGestalt software. KEGG analysis of the total gene set enrichment was created using ShinyGo 0.76.3 (39). The Circos plots were generated by connecting individual up and downregulated genes with specific hallmarks of cancer (40), considering all the cancer types (http://circos.ca/) (41). Alluvial diagrams were constructed with SankeyMATIC linking either up or downregulated genes with their particular cancer (https://sankeymatic.com/). PCA plot and Heatmap were created using the web tools ClustVis (http://biit.cs.ut.ee/clustvis/) (42) and Morpheus (https://software.broadinstitute.org/morpheus) (43).


3.      RESULTS

3.1. Identification of the genes over different cancer types correlated with melatonin.

     We initially searched for genes most relevant to the publication terms of “melatonin” and “cancer” using the platform Geneshot. A total of 152 genes were identified and ranked according to the appearance and number of publications. Based on term counts, the top 10 genes were ESR1CASP3TP53BCL2AKT1IL2IL6INSTNF, and HIF1A (Table S1). Among all genes, 15 of them were kinases-related genes (AKT1ERBB2MAPK3EGFRMTOREIF2AK3SRCATMCDK4KDRPRKCACSNK1EERN1MAPK8, and ROCK1) and three belonged to G-protein coupled receptors (MTNR1AMTNR1B, and OPN4). Specifically, MTNR1A and MTNR1B genes encodes both melatonin receptors termed MT1 and MT2, respectively. After identifying the genes using the highest confidence of 0.9, we identified 141 connected genes (Figure 1A). This complex network was again examined and 10 hub genes (STAT3, score 46.0; JUN, score 41.0; TP53, score 40.0; MAPK3, score 39.0; EP300, score 39.0; SRC, score 39.0; HSP90AA1, score 35.0; AKT1, score 34.0; ESR1, score 31.0; and IL6, score 31.0) were detected with potential interest in the context of melatonin and cancer (Figure 1B); these genes are involved with several pathways in cancer, estrogen signaling, and immune response (false discovery rate (FDR) < 10-7, combined score > 3,000), and most of them are transcription factors.

     Using KEGG database, the gene network was mainly associated with pathways in cancer, lipid and atherosclerosis, PI3K-Akt signaling, viral infections, proteoglycans in cancer among others (-Log10 FDR > 30; Figure 1C). Furthermore, the most representative GO biological processes involved with the genes associated with melatonin and cancer were cellular responses to cytokine stimuli, positive regulation of transcription, gene expression, intracellular signal transduction, metabolic processes, and regulation of apoptosis (Figure 1D, Table S2). To deeper explore this network with a pathological perspective related to cancer phenotypes, we used CancerGeneNet to develop a rational framework for the genes that are frequently mutated in cancers. Figure 1E summarizes the cellular regulatory circuits displaying membrane receptors and ligands, cytoplasmic signaling molecules, and nuclear activators (physical interactions with score 0.9).

3.2. Upregulated and downregulated genes in melatonin-treated cancers.

     Among the molecules associated with melatonin and cancer, we sorted them into upregulated and downregulated genes (Table S1). A total of 25 genes were upregulated by melatonin in different cancers (P<0.05); CASP3TP53 (score 22.0), ANXA5 (score 20.0), CASP9, and CDH1 (score 17.0), which are related to apoptosis, p53 signaling, tumor suppression, were recognized as top five hub genes by degree. Figure 2A depicts the upregulated gene network connected with top ranking neighbors (FBXO5CASP10XIAPTRIM25APPRNF4BCL2L1ELAVL1FBXW11, and PPP1CA). These upregulated genes were associated with important cancer hallmarks suggesting that melatonin may reverse some of these cancer phenotypes. Of those, resisting cell death, sustaining proliferative signaling, and induction of invasion and metastasis are major hallmarks involved with these genes in overall tumors (Figure 2B). The bar charts show different categories of GO terms involved with upregulated genes (Figure 2C). The top biological processes were response to stimulus, biological regulation, cell communication, and metabolic processing. Regarding the cellular components and molecular functions, we identified associations with the nucleus, membrane-enclosed lumen cytosol, membrane, protein-containing complex, and the main functions included protein binding, hydrolase and enzyme regulator activities, nucleic acid and ion binding, and others.

Figure 1.jpg

Fig. 1. Overall gene interactions and term enrichment.  

     A) Network interaction obtained from Geneshot list was generated in STRING and visualized in Cytoscape (PPI enrichment p-value < 1.0e-16; minimum confidence score: 0.9). B) Top 10 hub genes from the gene network were generated in Cytoscape (score > 31). C) KEGG analysis of the gene set was created in ShinyGo 0.76.3 (FDR < 0.05). D) Biological processes were based on functional enrichment of the gene set and reflect the Log10 p-value of the terms. E) Cellular scheme exhibiting the main molecular regulators of the gene list according to the most frequent mutations and cancer phenotypes (blue line = upregulation, red line = downregulation, dashed line = molecular binding, continuous line = direct interaction; score 0.9).

Figure 2.jpg

Fig. 2Upregulated target genes in melatonin-treated cancer and related functional enrichment.

     A) Network connecting the gene set (major circles = seed genes, minor circles = top ranking neighbors). B) Circos plot showing the associations of upregulated genes with cancer hallmarks; each color ribbon represents a specific target gene. C) Bar charts depicting GO categories (biological process, cellular components, and molecular functions) of the genes.

     Melatonin also showed an involvement with other genes, which were downregulated in several cancer types. We detected 51 downregulated genes (P<0.05) and the top five hub genes were STAT3HIF1A (score 49.0), CTNNB1 (score 48.0), AKT1 (score 47.0), and ESR1 (score 46.0); these genes are mostly involved with cell survival and differentiation, angiogenesis, epithelial-mesenchymal transition, and estrogen receptor binding. As shown in Figure 3A, the downregulated genes associated with different ranking neighbors (HRH1IL6RSH3GL2ZBTB16TP53ELAVLNTRK1BAG6COPS6, and UBQLN4). The most relevant functions related to the expression of these genes were sustained proliferative signaling, invasion and metastasis, cell death resistance, and replicative immortality (Figure 3B), all of which were potentially attenuated by the melatonin-linked downregulated genes. The terms enriched by the downregulated genes were similar to those of upregulated genes and biological processes were over-represented by the response to stimuli, biological regulation, metabolic processes, and multicellular organismal processes (Figure 3C). The cellular components included associations with membrane-enclosed vesicles, nucleus, protein-containing complex, and membrane, and the top categories of molecular functions were protein binding, ion and nucleic acid binding, transferase activity, and others.

 Figure 3.jpg 

Fig. 3Downregulated target genes in melatonin-treated cancer and related functional enrichment.

     A) Network connecting the gene set (major circles = seed genes, minor circles = top ranking neighbors). B) Circos plot showing the associations of downregulated genes with cancer hallmarks; each color ribbon represents a specific target gene. C) Bar charts depicting GO categories (biological process, cellular components, and molecular functions) of the genes.

     In addition to identifying GO terms linked to up- and downregulated genes altered by melatonin in tumors, we used BiNGO to deeply assess the overrepresentation of gene ontology categories in the biological networks (hypergeometric test, FDR < 0.01). The gene set was referred to as the molecular functions. For instance, while upregulated genes showed representative categories such as apoptotic protease activator, caspase activator activity, enzyme regulator and activator activities, protein binding, and TRAIL binding (Figure 4A), the downregulated genes were involved with protein kinase activities, transcription factor binding, protein, enzyme, DNA, and promoter bindings, transcription regulator and activator activities (Figure 4B).  

3.3.   Gene profile and functional enrichment are specifically associated with cancer type in response to melatonin.

     We investigated the diversity of targets by which melatonin interferes with each cancer type by examining the literature; only those genes involved in at least two different cancer types were computed. From the upregulated genes, 15 out of 25 were directly associated with several cancer types. The most abundant genes regulated by melatonin in the tumors were involved with apoptosis, cell adhesion, and oncogenic signaling (CASP3CDH1CYCS, CASP9, and DDIT3), whereas the top cancer types affected by the upregulated targets were breast cancer, colorectal cancer, and hepatocellular carcinoma (Figure 5 A). Conversely, 18 genes downregulated by melatonin interacted to a greater or lesser extent with 27 distinct cancer types. Specifically, BCL2HIF1AMMP9STAT3SNAI1CDH2CCND1, and CDK4, which are associated with kinase signaling, negative regulation of apoptosis, extracellular matrix remodeling, transcription regulator complex, cell migration and differentiation, were the most common genes attenuated by melatonin treatment, and breast cancer, colorectal cancer, gastric cancer, hepatocellular carcinoma, ovarian cancer, and osteosarcoma were the cancer types most influenced by these downregulated genes (Figure 5 B).

Figure 4.jpg

Fig. 4. The interactome of the molecular functions in response to melatonin in cancers.

     BiNGO was used to overrepresent the gene set functions. A) Network generated with upregulated genes. B) Network generated with downregulated genes. The hypergeometric test corrected by Benjamini-Hochberg (FDR < 0.01).

Figure 5.jpg

Fig. 5. Gene signatures and enrichment analysis of different cancer types treated with melatonin.

     Alluvial diagrams connecting the representative upregulated (A) and downregulated (B) genes by melatonin in respective cancer types were generated using the Sankeymatic online tool (http://sankeymatic.com/). Upregulated and downregulated plots showed respectively 80 and 108 flows linking genes and cancers. C) Heatmap top-ranked categories associated with different cancers in which genes were upregulated or downregulated by melatonin, computed by the gene set enrichment analysis using EnrichR. Signaling pathways were considered for over- or under-representation (Log2 combined score, minimum and maximum values). Rows and columns were clustered using Euclidean distance. D) Principal component analysis (PCA) based on functional enrichments by target genes in different cancer types.

     All over-represented genes upregulated or downregulated by melatonin (P<0.05) were selected for enrichment analysis using KEGG pathways (only tumors exhibiting > 5 genes affected by the indole were considered in this analysis; Figure 5 C). A total of 12 cancer types were more commonly associated with melatonin-related genes. Notably, a clear pattern of clusterization highlighted breast cancer, hepatocellular carcinoma, prostate cancer, neuroblastoma, oral cancer, and osteosarcoma as having more changes in functional enrichment-related genes (higher combined score) compared to glioma, gastric, colorectal, ovarian, gallbladder, and lung cancers (lower combined score). The main functions in the cancer subset with higher combined scores were the p53 signaling pathway, circadian rhythms, FOXO signaling, mitophagy, the hedgehog signaling pathway, apoptosis, HIF-1 signaling, JAK-STAT signaling, TNF signaling, cellular senescence, Wnt signaling, tight junction, and the Hippo signaling pathway. Conversely, the functions in the cancer types displaying a lower combined score were mainly associated with circadian rhythms, necroptosis, and sphingolipid signaling pathway. Based on these functional mechanisms, the principal component analysis (PCA) showed that neuroblastoma, lung cancer, glioma, and oral cancer were distinguishable from each other and among all cancer types regarding melatonin’s action (Figure 5 D). Next, we extrapolated these findings associating individual cancer signatures to patients’ overall survival.

3.4. Melatonin alters gene signature related to a shorter overall survival of patients.

     Considering that the vast majority of the molecular pathways and the gene sets regulated by the melatonin are related to cancer aggressiveness, we investigated whether the melatonin-downregulated genes indeed predict a worse prognosis for patients. To explore survival data of cancer patients, we used GEPIA2 to determine the prognostic index correlating survival curve with risk assessment. These patients were stratified into two groups (50% at high- and 50% at low-risk; median values) so that to individually evaluate gene signatures expressed in each cancer type available in the TCGA platform. The shorter survival of patients was linked to a high expression of specific genes in BLCA (SNAI1MTORCDH2, and FN1; hazard ratio (HR) = 1.4, log-rank p-value = 0.027; Figure 6A), BRCA (STAT3HIF1ASNAI1, and CDH2; HR = 1.4, log-rank p-value = 0.037; Figure 6B), COAD (STAT3MTORSNAI1POU5F1, and CDH2; HR = 1.7, log-rank p-value = 0.04; Figure 6C), STAD (BCL2HIF1ASNAI1EDN1, and FN1; HR = 1.4, log-rank p-value = 0.032; Figure 6D), GBM (STAT3HIF1AMMP2, and MMP9; HR = 2, log-rank p-value = 0.00033; Figure 6E), LIHC (BCL2IL6STAT3MMP9, CCND1, HIF1AMTOR, and CDK4; HR = 1.6, log-rank p-value = 0.0062; Figure 6F), LUAD (CTNNB1, POU5F1, and CDK4; HR = 1.3, log-rank p-value = 0.04; Figure 6G), and OV (HIF1ASNAI1, and MTOR; HR = 1.3, log-rank p-value = 0.026; Figure 6H). Collectively, the overall survival map revealed that FN1CDH2SNAI1HIF1AMMP2, and MMP9 genes were important contributors to shortening the survival of patients in most of the tumors (Figure 6 I). The gene signatures are also key candidates for disease-free survival of patients especially in BLCA, COAD, GBM, LIHC, and LUAD (Figure 6 J).

     To deeper understand the gene expression signatures on specific cancer cell type, we used a curated cancer cell atlas providing collected, annotated, and analyzed cancer scRNA-seq datasets. Based on heatmap clustering, there is a clear profile of gene expression in endothelial cells of BRCA, LUAD, OV, COAD, and in microglia of GBM patients. This expression pattern was closely associated with that found in fibroblasts of BRCA, LUAD, OV, and COAD, and also with genes majorly expressed in the malignant cells of LUAD, OV, and COAD (Figure 7). More generally, FN1 and MMP2 genes were highly expressed in tumor-associated fibroblasts of patients with BRCA, OV, LUAD, and COAD. The genes CTNNB1EDN1SNAI1MTORHIF1ASTAT3, and IL6 were mostly increased in endothelial cells. Furthermore, the genes CCND1, CDK4, and POU5F1 had increased expression in malignant cells of patients with LUAD, OV, and COAD. Specifically, BRCA patients also had high expression of CCND1 gene in malignant cells whereas MMP9 and HIF1A genes were increased in macrophages. COAD patients showed high expression of EDN1 gene in the malignant cells, and MMP9HIF1A, and IL6 in tumor-associated macrophages. OV patients showed increased expression of FN1 and HIF1A genes in macrophages whereas STAT3 was in mesothelial cells, BCL2 in T cells, and CTNNB1 and EDN1 genes in endothelial cells. Finally, GBM patients showed moderate expression of BCL2 in astrocytes, CDK4 and CDH2 in malignant cells, and CTNNB1HIF1A, and STAT3 in microglia (Figure 7). Exploring scRNA-seq datasets aid revealing upregulated targets in specific cancer cell types which are potentially affected by melatonin treatment.

Figure 6.jpg

Fig. 6. The impact of gene signatures regulated by melatonin on the prognosis of cancer patients.

     Overall survival curves comparing patients at high (red line) and low (blue line) risks based on overexpressed genes affected by melatonin in BLCA (A), BRCA (B), COAD (C), STAD (D), GBM (E), LIHC (F), LUAD (G), and OV (H). Kaplan–Meier curves determined by univariate Cox regression analysis show the hazard ratio (HR) corresponding to 95% confidence intervals, log-rank p-value, and the number of patients by groups. Heatmaps showing the main gene set belonging to each tumor signature relative to overall survival (I) and disease-free survival (J). Values are based on Log10 HR (+ 0.2 / - 0.2). BLCA, bladder urothelial carcinoma; BRCA, breast invasive carcinoma; COAD, colon adenocarcinoma; STAD, stomach adenocarcinoma; GBM, glioblastoma multiforme; LIHC, liver hepatocellular carcinoma, LUAD, lung adenocarcinoma; OV, ovarian serous cystadenocarcinoma.

 

Figure 7.jpg


Fig. 7. Single cell RNA-seq analysis reveals a specific cell type expression profile of melatonin-related genes in BRCA, COAD, GBM, LUAD, and OV.

     Unsupervised hierarchical clustering (Euclidean distance) of melatonin-related genes is presented in a heatmap according to average expression scaled (color) and percentage of cells expressing the gene (circle size). Heatmap was generated and analyzed using Morpheus online tool (https://software.broadinstitute.org/morpheus/). BRCA, breast invasive carcinoma; COAD, colon adenocarcinoma; GBM, glioblastoma multiforme; LUAD, lung adenocarcinoma; OV, ovarian serous cystadenocarcinoma.


4.      DISCUSSION

      The understanding of cancer molecular networks targeted by anticancer agents may help control its aggressiveness while favoring the development of novel personalized therapies (44). Here, we unveiled which genes whereby melatonin interacts with cancer and defined a core framework for specific cancer types. Through Geneshot, a search engine that ranks lists of relevant genes based on arbitrary terms, we first identified the genes involved with melatonin and cancer, and further individualized each gene set for each cancer type. A total of 141 highly connected genes (confidence score = 0.9) were found in the general cancer network, some of which only interact with melatonin at one level (MTNR1AMTNR1BASMT, and AANAT), while a vast majority was composed by driver genes undergoing frequent mutations in cancer studies. As expected, these genes enriched functions are commonly observed in many cancer processes (e.g., responses to cytokine stimuli, gene expression, intracellular signal transduction, apoptotic and metabolic processes, and others). Regarding melatonin production, ASMTL gene has been associated with higher capacity of melatonin synthesis in females, thus enhancing protection against stressful conditions (45). We previously showed that ASMTL gene positively correlated with metabolism-related genes (e.g., oxidative phosphorylation) in ovarian cancer (22).  The hub genes mostly affected by melatonin in cancers were STAT3, JUN, TP53MAPK3EP300SRCHSP90AA1AKT1ESR1, and IL6; these molecules are involved with IL-6/JAK/STAT3 signaling, TNF-alpha signaling via NF-kB, hypoxia, p53 pathway, and epithelial-mesenchymal transition (FDR < 0.01). Examination of melatonin’s main actions in multiple cancers from a gene expression standpoint allowed for reinforcing various anticancer properties, highlighting its role as an immune enhancer, anti-inflammatory, anti-angiogenic, anti-metastatic, and pro-apoptotic agent (46). Of note, the anticancer activity of melatonin is not limited to the aforementioned mechanisms and signaling pathways, but rather represents the central elements of its gene interaction.

     To better understand the involvement of melatonin-related genes with cancer hallmarks, we grouped them as upregulated and downregulated genes. Since hallmark capabilities are regulated by partially redundant signaling pathways, their value depends on the cancer's underlying molecular features (40). The genes upregulated by melatonin classified as tumor suppressor genes and transcription factors potentially alter the sustained proliferative signaling and evasion of growth suppressors. Also, genes participating in the TNF receptor pathway, apoptosis pathway, Fas signaling, caspase family, CARD family, TP53 activators and regulators involved with cell death resistance were upregulated by melatonin. Other important features counteracted by melatonin were invasion and metastasis (with the genes classified as cell-matrix adhesion and Rho signaling, proteolysis, cell motility, invasive projections, chemotaxis, cell-cell adhesion, and transcription factors), replicative immortality (genes classified as proteases and telomerase regulators), and genomic instability (genes classified as cell cycle checkpoint genes). The downregulated genes by melatonin were numerous and negatively impacted cancer cell activities. The hallmark and respective internal features associated with those genes were as follow: sustaining proliferative signaling (transcription factors, cell cycle, oncogenes, and multiple features), activating invasion and metastasis (cell-cell adhesion, Rho signaling, cell motility, Cdc42 signaling, matrix metalloproteinases, cell-matrix adhesion, invasive projections, receptor, filopodia, Rac signaling, proteolysis, and chemotaxis), resisting cell death (cell cycle pathway members, TP53 activators and cooperators, TNF superfamily, and Bcl2 family-regulated pathway), enabling replicative immortality (senescence pathway and telomerase regulation), inducing angiogenesis (growth factors and receptors, cytokines and chemokines), genome instability (DNA repair, chromatin modification molecules, and cell cycle checkpoint genes), and deregulation of cellular energetics (activator molecules). A comprehensive review documented the anticancer mechanisms of melatonin against cancer hallmarks associated with cell signaling pathways using in vivo and in vitro studies (12). These findings provide a foundation for the melatonin-mediated attenuation of cancer hallmarks at the gene expression level. This gene regulation is strongly reflected in the modulation of important molecular functions related to cancer progression as those described by BiNGO analysis.

     Upon melatonin treatment, most of the upregulated genes were detected in breast and colorectal cancers, and in hepatocellular carcinoma (shared genes included TP53CASP3, and CYCS), whereas the downregulated genes were primarily observed in breast cancer, colorectal cancer, gastric cancer, hepatocellular carcinoma, ovarian cancer, and osteosarcoma (shared genes were BCL2 and MMP9). This is attributed to the high number of melatonin-related studies and a larger number of genetic signatures per cancer compared to that observed in a smaller set of tumors. Through KEGG analysis, we explored the similarities and differences among melatonin-treated cancers considering a particular gene set and terms enrichment. For instance, breast cancer, hepatocellular carcinoma, prostate cancer, neuroblastoma, oral cancer, and osteosarcoma showed similar enriched signaling pathways in contrast to glioma, gastric, colorectal, ovarian, gallbladder, and lung cancers; this information may be utilized for the planning of personalized therapies using melatonin. Melatonin has been documented to directly affect cancer cells via apoptotic, antiangiogenic, antiproliferative, and metastasis-inhibitory signaling pathways (47). As a non-toxic, potent, and synergistic therapeutic agent, melatonin combined with conventional drugs improves cancer sensitivity and life quality of patients (48, 49). A comprehensive review involving cellular, animal, and human studies summarized the underlying molecular mechanisms, epidemiology, risk factors, and clinical significance of melatonin alone or in combination with other chemotherapies against breast, colorectal, liver, lung, prostate, skin, cervical, and ovarian cancers (47). Most of these regulatory events are dependent on gene networks explored in our current study.    

     Recently compiled information on the ability of melatonin to regulate diverse signaling pathways in different cancer types revealed its involvement with the PI3K-Akt-mTOR pathway, MAPK signaling, Wnt/β-catenin signaling, Notch and NF-kB signaling, JAK/STAT signaling, IGF, VEGF, FGF, and PDGF signaling pathway, and Hedgehog pathway (50). Many of these pathways and others were individually analyzed in each tumor type. The set of genes that enriched the most signaling pathways (e.g., FOXO signaling, p53 signaling, mitophagy, hedgehog pathway, HIF-1 signaling, proteoglycans in cancer, JAK/STAT signaling, PI3K-Akt and TNF signaling, cellular senescence, tight junctions, and Wnt signaling) were shared especially among breast cancer, hepatocellular carcinoma, prostate cancer, and neuroblastoma. These signaling pathways targeted by melatonin are differentially modulated in different cancer phenotypes and represent potentially important mechanisms for cancer biology and behavior. From this analysis, new therapeutic proposals for such cancers could include conventional chemotherapy or radiotherapy with melatonin. A recent study reported an important role of melatonin on the inhibition of VEGF and inactivation of HIF-1α showing that melatonin neutralizes pro-angiogenic while potentiating antiangiogenic effects induced by chemo or radiotherapeutic agents (51); this cooperative action of melatonin highlights its ability to turn cancer totally resistant to chemotherapeutics into a more sensitive tumor state. Additional studies are urgently needed to test if other regulatory factors or genes are involved with melatonin-induced chemosensitizing effects.    

     To verify the significance of these melatonin-related gene signatures in the context of cancer patients, we explored the melatonin-downregulated genes to uncover their impact on overall survival of patients based on the TCGA dataset. Clearly, functional gene signatures for each cancer type revealed shorter overall patient survival since the high-signature groups presented the worst prognosis; in this case, melatonin has the potential to improve patient survival by reversing the overexpression of specific cancer-associated genes. Together, we described a relevant set of genes altered by melatonin with a negative impact on overall survival of patients with bladder, breast, colorectal, gastric, lung, and ovarian cancers, in addition to glioblastoma multiforme, and hepatocellular carcinoma. Considering specific cancer types, we further identified the transcriptional profile of single cell signatures which are affected by melatonin; specific genes were similarly clustered in malignant cells, macrophages, fibroblasts, and endothelial cells. Particularly, CCND1 and CDK4 genes, which are cell cycle regulators, were increased in malignant cells, and FN1 and MMP2 genes, especially involved in cell invasion and metastasis, were highly expressed in cancer-associated fibroblasts; these molecules are targeted by melatonin and may serve as candidate biomarkers. Although mostly represented with few cells, the CTNNB1EDN1IL6SNAI1MTOR, and STAT3 genes showed a similar expression profile in cancer-associated endothelial cells.

     Mechanistically, up and downregulated genes mediated by melatonin may be a result of its interference in any element of gene expression machinery, potentially being modulated by non-coding RNAs and other epigenetic-related mechanisms (e.g., methylation, acetylation or LNCmRNA and miRNA regulation). The epigenetic regulation of different genes mediated by melatonin has already been documented in breast cancer (52), oral cancer (53), prostate cancer (54), malignant glioma (55), and others. Additional studies that focus on gene dynamics in the presence of melatonin will be helpful in identifying more precisely its biological mechanisms of action in different cancer phenotypes. The current data, however, can be used to guide treatment protocols. Melatonin exerts oncostatic activities in numerous malignancies-related events with important roles in cell cycle, apoptosis, and autophagy-related processes (56); these actions are especially mediated through receptors MT1 and MT2, cytosolic binding sites including calmodulin and quinone reductase II enzyme, and, although debatable, possibly involving activation of nuclear receptors (RZR/ROR). A Pan-cancer study, involving 33 cancer types from TCGA database, identified different melatonin regulators (57); in the genomic landscape, there was heterozygous amplification of AANAT and GPR50 genes, and heterozygous deletion of PER3CYP2C19, and MTNR1A as the dominant alterations.

     In addition to its intricate actions on molecular signatures, melatonin has numerous metabolic actions that affect the biology of tumors (58-61). The bioavailability of melatonin depends on its pharmacokinetic features, including absorption, metabolism, and elimination; these features may vary among individuals and should be kept in mind when establishing a dose protocol to support its use in cancer prevention and treatment. Because in vivo and in vitro studies have broadly detailed the mechanism of melatonin’s actions on different cancer types, it is imperative to maintain the use of melatonin as an adjuvant opportunity in clinical contexts.

     This study provides a new integrative view of the genes associated with melatonin into a cancer scenario. The general characterization of up and downregulated molecules in melatonin-treated cancers revealed important targets which are dysregulated in several cancer types, including rare cancers, and could be exploited to novel experiments for different cancer modalities. Melatonin is capable of influencing gene signatures associated with cancer aggressiveness, thus reflecting its anti-cancer regulatory role at both the molecular and functional levels. It should be disclosed that the number and interaction of genes for each cancer type is based on literature data and in no way does it restrict or prioritize some cancer types; as a result, we feel that melatonin has even a greater involvement in controlling the tumor gene expression. We also note that global transcriptomic analysis provides a more complete gene landscape to guide the molecular responses of melatonin in cancer studies.  

 

ACKNOWLEDGMENTS

     The following grants supported the development of this study: National Council for Scientific and Technological Development (CNPq, Process number 304108/2020-0 to LGAC) and São Paulo Research Foundation (FAPESP grant # 2021/12971-7 to LGAC).

 

AUTHORSHIP

     LGAC: conceptualization, data acquisition and curation, and drafted the manuscript. RFC, RFD, VLSC, SSC: data processing and collection and drafted the manuscript. FRFS, DACZ: collection, interpretation, and data analysis. All authors significantly contributed to the compilation of the literature and approved the final version of the manuscript.

 

CONFLICTS OF INTEREST

     Theauthors have no relevant financial or non-financial interests to disclose.

 

DATA AVAILABILITY

     All data supporting the findings of this study are available within the paper. Data will be made available on request.

 

REFERENCES

 

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F (2021) Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA. Cancer J. Clin. 71 (3): 209-249. DOI: 10.3322/caac.21660.

  2. Bray F, Laversanne M, Weiderpass E, Soerjomataram I (2021) The ever-increasing importance of cancer as a leading cause of premature death worldwide. Cancer 127 (16): 3029-3030. DOI: 10.1002/cncr.33587.

  3. Fitzmaurice C, Abate D, Abbasi N, Abbastabar H, Abd-Allah F, Abdel-Rahman O, Abdelalim A, Abdoli A, Abdollahpour I, Abdulle ASM, Abebe ND, et al. (2019) Global burden of disease cancer: Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: A systematic analysis for the global burden of disease study. JAMA Oncol. 5 (12): 1749-1768. DOI: 10.1001/jamaoncol.2019.2996).

  4. Nussinov R, Tsai CJ, Jang H (2021) Anticancer drug resistance: An update and perspective. Drug Resist. Updat. 59: 100796. DOI: 10.1016/j.drup.2021.100796.

  5. Tufail AB, Ma YK, Kaabar MKA, Martinez F, Junejo AR, Ullah I, Khan R (2021) Deep learning in cancer diagnosis and prognosis prediction: A minireview on challenges, recent trends, and future directions. Comput. Math. Methods Med. 2021: 9025470. DOI: 10.1155/2021/9025470.

  6. Hanahan D (2022) Hallmarks of cancer: New dimensions. Cancer Discov. 12 (1): 31-46. DOI: 10.1158/2159-8290.CD-21-1059.

  7. Creixell P, Schoof EM, Erler JT, Linding R (2012) Navigating cancer network attractors for tumor-specific therapy. Nat. Biotechnol. 30 (9): 842-848. DOI: 10.1038/nbt.2345.

  8. Davoodvandi A, Nikfar B, Reiter RJ, Asemi Z (2022) Melatonin and cancer suppression: insights into its effects on DNA methylation. Cell Mol. Biol. Lett. 27 (1): 73. DOI: 10.1186/s11658-022-00375-z.

  9. Hardeland R (2014) Melatonin, noncoding RNAs, messenger RNA stability and epigenetics--evidence, hints, gaps and perspectives. Int. J. Mol. Sci. 15 (10): 18221-18252. DOI: 10.3390/ijms151018221.

  10. Chuffa LGA, Carvalho RF, Justulin LA, Cury SS, Seiva FRF, Jardim-Perassi BV, Zuccari D, Reiter RJ (2020) A meta-analysis of microRNA networks regulated by melatonin in cancer: Portrait of potential candidates for breast cancer treatment. J. Pineal Res. 69 (4): e12693. DOI: 10.1111/jpi.12693.

  11. Cohen M, Lippman M, Chabner B (1978) Role of pineal gland in aetiology and treatment of breast cancer. Lancet 2 (8094): 814-816. DOI: 10.1016/s0140-6736(78)92591-6.

  12. Talib WH (2018) Melatonin and cancer hallmarks. Molecules 23 (3): 518. DOI: 10.3390/molecules23030518.

  13. Hill SM, Belancio VP, Dauchy RT, Xiang S, Brimer S, Mao L, Hauch A, Lundberg PW, Summers W, Yuan L, Frasch T, Blask DE (2015) Melatonin: an inhibitor of breast cancer. Endocr. Relat. Cancer. 22 (3): R183-204. DOI: 10.1530/ERC-15-0030.

  14. Mediavilla MD, Sanchez-Barcelo EJ, Tan DX, Manchester L, Reiter RJ (2010) Basic mechanisms involved in the anti-cancer effects of melatonin. Curr. Med. Chem. 17 (36): 4462-4481. DOI: 10.2174/092986710794183015.

  15. Cruz EMS, Concato VM, de Morais JMB, Silva TF, Inoue FSR, de Souza Cremer M, Bidoia DL, Machado RRB, de Almeida Chuffa LG, Mantovani MS, Panis C, Pavanelli WR, et al. (2023) Melatonin modulates the Warburg effect and alters the morphology of hepatocellular carcinoma cell line resulting in reduced viability and migratory potential. Life Sci. 319: 121530. DOI: 10.1016/j.lfs.2023.121530.

  16. Ashrafizadeh M, Ahmadi Z, Yaribeygi H, Sathyapalan T, Jamialahmadi T, Sahebkar A (2021) Antitumor and protective effects of melatonin: The potential roles of MicroRNAs. Adv. Exp. Med. Biol. 1328: 463-471. DOI: 10.1007/978-3-030-73234-9_31.

  17. Chuffa LG, Alves MS, Martinez M, Camargo IC, Pinheiro PF, Domeniconi RF, Junior LA, Martinez FE (2016) Apoptosis is triggered by melatonin in an in vivo model of ovarian carcinoma. Endocr. Relat. Cancer. 23 (2): 65-76. DOI: 10.1530/ERC-15-0463.

  18. Zonta YR, Martinez M, Camargo IC, Domeniconi RF, Lupi Junior LA, Pinheiro PF, Reiter RJ, Martinez FE, Chuffa LG (2017) Melatonin reduces angiogenesis in serous papillary ovarian carcinoma of ethanol-preferring rats. Int. J. Mol. Sci. 18 (4): 763. DOI: 10.3390/ijms18040763.

  19. Blask DE, Hill SM, Dauchy RT, Xiang S, Yuan L, Duplessis T, Mao L, Dauchy E, Sauer LA (2011) Circadian regulation of molecular, dietary, and metabolic signaling mechanisms of human breast cancer growth by the nocturnal melatonin signal and the consequences of its disruption by light at night. J. Pineal Res. 51 (3): 259-269. DOI: 10.1111/j.1600-079X.2011.00888.x.

  20. Bonmati-Carrion MA, Tomas-Loba A (2021) Melatonin and cancer: A polyhedral network where the source matters. Antioxidants 10 (2): 210. DOI: 10.3390/antiox10020210.

  21. Stevens RG, Brainard GC, Blask DE, Lockley SW, Motta ME (2014) Breast cancer and circadian disruption from electric lighting in the modern world. CA. Cancer J. Clin. 64 (3): 207-218. DOI: 10.3322/caac.21218.

  22. Cucielo MS, Cesario RC, Silveira HS, Gaiotte LB, Dos Santos SAA, de Campos Zuccari DAP, Seiva FRF, Reiter RJ, de Almeida Chuffa LG (2022) Melatonin reverses the warburg-type metabolism and reduces mitochondrial membrane potential of ovarian cancer cells independent of MT1 receptor activation. Molecules 27 (14): 4350. DOI: 10.3390/molecules27144350.

  23. Gaiotte LB, Cesário RC, Silveira HS, De Morais Oliveira DA, Cucielo MS, Romagnoli GG, Kaneno R, De Campos Zuccari DAP, Reiter RJ, Chuffa LGdA (2022) Combination of melatonin with paclitaxel reduces the TLR4-mediated inflammatory pathway, PD-L1 levels, and survival of ovarian carcinoma cells. Melatonin Res. 5 (1): 34-51. DOI: 10.32794/mr112500118.

  24. Chuffa LGDA, Carvalho RF, Seiva FR, Zuccari DAPdC, Reiter RJ (2021) Melatonergic index as a prognostic biomarker of reproductive organ cancers: correlations with metabolic parameters as well as clock genes PER1 and TIMELESS. Melatonin Res. 4 (2): 299-315. DOI: 10.32794/mr11250096.

  25. Kinker GS, Oba-Shinjo SM, Carvalho-Sousa CE, Muxel SM, Marie SK, Markus RP, Fernandes PA (2016) Melatonergic system-based two-gene index is prognostic in human gliomas. J. Pineal Res. 60 (1): 84-94. DOI: 10.1111/jpi.12293.

  26. Lv JW, Zheng ZQ, Wang ZX, Zhou GQ, Chen L, Mao YP, Lin AH, Reiter RJ, Ma J, Chen YP, Sun Y (2019) Pan-cancer genomic analyses reveal prognostic and immunogenic features of the tumor melatonergic microenvironment across 14 solid cancer types. J. Pineal Res. 66 (3): e12557. DOI: 10.1111/jpi.12557.

  27. Maleki M, Khelghati N, Alemi F, Younesi S, Asemi Z, Abolhasan R, Bazdar M, Samadi-Kafil H, Yousefi B (2021) Multiple interactions between melatonin and non-coding RNAs in cancer biology. Chem. Biol. Drug Des. 98 (3): 323-340. DOI: 10.1111/cbdd.13849.

  28. Lachmann A, Schilder BM, Wojciechowicz ML, Torre D, Kuleshov MV, Keenan AB, Ma'ayan A (2019) Geneshot: search engine for ranking genes from arbitrary text queries. Nucleic Acids Res. 47 (W1): W571-W577. DOI: 10.1093/nar/gkz393.

  29. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, Jensen LJ, Mering CV (2019) STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 (D1): D607-D613. DOI: 10.1093/nar/gky1131.

  30. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11): 2498-2504. DOI: 10.1101/gr.1239303.

  31. Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY (2014) cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. (4): S11. DOI: 10.1186/1752-0509-8-S4-S11.

  32. Iannuccelli M, Micarelli E, Surdo PL, Palma A, Perfetto L, Rozzo I, Castagnoli L, Licata L, Cesareni G (2020) CancerGeneNet: linking driver genes to cancer hallmarks. Nucleic Acids Res. 48 (D1): D416-D421. DOI: 10.1093/nar/gkz871.

  33. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, Clark NR, Ma'ayan A (2013) Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics 14: 128. DOI: 10.1186/1471-2105-14-128.

  34. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, McDermott MG, Monteiro CD, et al. (2016) Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44 (W1): W90-97. DOI: 10.1093/nar/gkw377.

  35. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B (2019) WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 47 (W1): W199-W205. DOI: 10.1093/nar/gkz401.

  36. Tang Z, Kang B, Li C, Chen T, Zhang Z (2019) GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 47 (W1): W556-W560. DOI: 10.1093/nar/gkz430.

  37. Satija R, Farrell JA, Gennert D, Schier AF, Regev A (2015) Spatial reconstruction of single-cell gene expression data. Nat. Biotechnol. 33 (5): 495-502. DOI: 10.1038/nbt.3192.

  38. Maere S, Heymans K, Kuiper M (2005) BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21 (16): 3448-3449. DOI: 10.1093/bioinformatics/bti551.

  39. Ge SX, Jung D, Yao R (2020) ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics. 36 (8): 2628-2629. DOI: 10.1093/bioinformatics/btz931.

  40. Menyhart O, Harami-Papp H, Sukumar S, Schafer R, Magnani L, de Barrios O, Gyorffy B (2016) Guidelines for the selection of functional assays to evaluate the hallmarks of cancer. Biochim. Biophys. Acta. 1866 (2): 300-319. DOI: 10.1016/j.bbcan.2016.10.002.

  41. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA (2009) Circos: an information aesthetic for comparative genomics. Genome Res. 19 (9): 1639-1645. DOI: 10.1101/gr.092759.109.

  42. Metsalu T, Vilo J (2015) ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res. 43 (W1): W566-570. DOI: 10.1093/nar/gkv468.

  43. Starruss J, de Back W, Brusch L, Deutsch A (2014) Morpheus: a user-friendly modeling environment for multiscale and multicellular systems biology. Bioinformatics 30 (9): 1331-1332. DOI: 10.1093/bioinformatics/btt772.

  44. Tran TD, Pham DT (2021) Identification of anticancer drug target genes using an outside competitive dynamics model on cancer signaling networks. Sci. Rep. 11 (1): 14095. DOI: 10.1038/s41598-021-93336-z.

  45. Tan D-X, Hardeland R (2021) The reserve/maximum capacity of melatonin's synthetic function for the potential dimorphism of melatonin production and its biological significance in mammals. Molecules 26 (23): 7302. DOI: 10.3390/molecules26237302.

  46. Bhattacharya S, Patel KK, Dehari D, Agrawal AK, Singh S (2019) Melatonin and its ubiquitous anticancer effects. Mol. Cell. Biochem. 462 (1-2): 133-155. DOI: 10.1007/s11010-019-03617-5.

  47. Gurunathan SQasim MKang M-H, Kim J-H (2021) Role and therapeutic potential of melatonin in various type of cancers. Onco Targets Ther14: 2019–2052. DOI: 10.2147/OTT.S298512.

  48. Zhao Y, Wang C, Goel A (2022) A combined treatment with melatonin and andrographis promotes autophagy and anticancer activity in colorectal cancer. Carcinogenesis 43 (3): 217-230. DOI: 10.1093/carcin/bgac008.

  49. Reiter RJ, Sharma R, Rosales-Corral S, Manucha W, Chuffa LGA, Zuccari DAPC (2021) Melatonin and pathological cell interactions: Mitochondrial glucose processing in cancer cells. Int. J. Mol. Sci22 (22): 12494. DOI: 10.3390/ijms222212494.

  50. Mihanfar A, Yousefi B, Azizzadeh B, Majidinia M (2022) Interactions of melatonin with various signaling pathways: implications for cancer therapy. Cancer Cell Int. 22 (1): 420. DOI: 10.1186/s12935-022-02825-2.

  51. González AAlonso-González CGonzález-González AMenéndez-Menéndez JCos SMartínez-Campa C (2021) Melatonin as an adjuvant to antiangiogenic cancer treatments. Cancers 13 (13): 3263. DOI: 10.3390/cancers13133263.

  52. Korkmaz A, Sanchez-Barcelo EJ, Tan DX, Reiter RJ (2009) Role of melatonin in the epigenetic regulation of breast cancer. Breast Cancer Res. Treat115 (1): 13-27. DOI: 10.1007/s10549-008-0103-5.

  53. Capote-Moreno A, Ramos E, Egea J, López-Muñoz F, Gil-Martín E, Romero A (2019)  Potential of melatonin as adjuvant therapy of oral cancer in the era of epigenomics. Cancers 11 (11): 1712. DOI: 10.3390/cancers11111712.

  54. Zhou L, Zhang CYang X, Liu L, Hu J, Hou Y, Tao HSugimura H, Chen Z, Wang L, Chen K (2021) Melatonin inhibits lipid accumulation to repress prostate cancer progression by mediating the epigenetic modification of CES1. Clin. Transl. Med. 11 (6): e449. DOI: 10.1002/ctm2.449.

  55. Davoodvandi A, Nikfar B, Reiter RJ, Asemi Z (2022) Melatonin and cancer suppression: insights into its effects on DNA methylation. Cell Mol. Biol. Lett27 (1): 73. DOI: 10.1186/s11658-022-00375-z.

  56. Targhazeh N, Reiter RJ, Rahimi M, Qujeq D, Yousefi T, Shahavi MH, Mir SM (2022) Oncostatic activities of melatonin: Roles in cell cycle, apoptosis, and autophagy. Biochimie200: 44-59. DOI: 10.1016/j.biochi.2022.05.008.

  57. Zhang JJiang H, Du K, Xie T, Wang B, Chen C, Reiter RJ, Cen B, Yuan Y (2021)  Pan-cancer analyses reveal genomics and clinical characteristics of the melatonergic regulators in cancer. J. Pineal Res71 (3): e12758. DOI: 10.1111/jpi.12758.

  58. Chuffa LGA, Reiter RJ, Lupi LA (2017) Melatonin as a promising agent to treat ovarian cancer: molecular mechanisms. Carcinogenesis 38 (10): 945-952. DOI: 10.1093/carcin/bgx054.

  59. de Almeida Chuffa LG, Seiva FRF, Cucielo MS, Silveira HS, Reiter RJ, Lupi LA (2019) Mitochondrial functions and melatonin: a tour of the reproductive cancers. Cell Mol. Life Sci. 76 (5): 837-863. DOI: 10.1007/s00018-018-2963-0.

  60. Li Y, Li S, Zhou Y, Meng X, Zhang JJ, Xu DP, Li HB (2017) Melatonin for the prevention and treatment of cancer. Oncotarget 8 (24): 39896-39921. DOI: 10.18632/oncotarget.16379.

  61. Reiter RJ, Sharma R, Ma Q, Rorsales-Corral S, de Almeida Chuffa LG (2020) Melatonin inhibits Warburg-dependent cancer by redirecting glucose oxidation to the mitochondria: a mechanistic hypothesis. Cell Mol. Life Sci. 77 (13): 2527-2542. DOI: 10.1007/s00018-019-03438-1.

CCBY.png

This work is licensed under a Creative Commons Attribution 4.0 International License