Novel Oxidative Stress-related Prognostic Biomarkers for Melanoma Associated With Tumor Metastasis

doi:10.21203/rs.3.rs-54156/v1

Download PDF

Primary research

Novel Oxidative Stress-related Prognostic Biomarkers for Melanoma Associated With Tumor Metastasis

https://doi.org/10.21203/rs.3.rs-54156/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 26 Feb, 2021

Read the published version in Medicine →

Version 1

posted

You are reading this latest preprint version

Background

Skin cutaneous melanoma (SKCM) is a prevalent skin cancer whose metastatic form is dangerous due to its high morbidity and mortality. Previous studies have systematically established the vital role of oxidative stress (OS) in melanoma progression. This study aimed to identify prognostic OS genes closely associated with SKCM and illustrate their potential mechanisms.

Methods

Transcriptome data and corresponding clinical traits of patients with SKCM were retrieved from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. A weighted gene co-expression network analysis was conducted to identify relationships between clinical features and OS genes in specific modules. Subsequently, Cox regression analysis was performed on candidate OS genes; four hub prognosis-associated OS genes (AKAP9, VPS13C, ACSL4, and HMOX2) were identified to construct a prognostic model.

Results

After a series of bioinformatics analysis, our prognostic model was identified significantly associated with the overall survival of patients with SKCM and metastatic ability of the cancer. Furthermore, our risk model demonstrated improved diagnostic accuracy in TCGA and GEO cohorts. In addition, we established two nomograms based on either risk score or hub genes, which displayed favorable discriminating ability for SKCM.

Conclusions

Together, our results provide novel insight into the potential applications of OS-associated genes in SKCM.

Cancer Biology

Skin cutaneous melanoma

Oxidative stress

Prognosis

Metastatic

Integrated bioinformatics analysis

Skin cutaneous melanoma (SKCM), a highly aggressive malignant tumor, is a serious threat to human health [1]. According to the GLOBOCAN database, 287,723 new patients with melanoma were diagnosed worldwide in 2018, of which 21.1% died [2]. Occurrence and progression of SKCM are significantly associated with the degree of skin pigmentation. Indeed, patients with low skin pigmentation content exposed to ultraviolet radiation comprise the majority of SKCM cases [3, 4]. Metastatic SKCM, with high rates of metastasis [5], is the most common form and is regarded as the leading cause of skin cancer deaths [6]. Globally, the 10-year overall survival of patients with stage I and II SKCM ranges from 98–75% [7], but one third of patients will progress to metastatic melanoma [8], making early diagnosis and individualized treatments vital. Currently, the optimal management of metastatic melanoma comprises targeted therapy and checkpoint suppression [9, 10]. However, therapeutic effectiveness remains limited due to lack of knowledge regarding the intricate mechanism underlying SKCM metastasis. In recent years, many studies have focused on identifying potential biomarkers for metastatic melanoma [11, 12]. Unfortunately, these developments have been inadequate for improving early diagnosis and prognosis of SKCM. Thus, development of more effective biomarkers is urgently needed for accurate identification and prognostication of SKCM.

The unequivocal mechanism leading to SKCM metastasis remains poorly understood; however, oxidative stress (OS) has been proposed as an important factor driving tumorigenesis and cancer progression through excessive production of reactive oxygen species (ROS) [13–15]. Imbalances between antioxidants and oxidants favor the oxidants, which potentially leads to biological damage, termed OS. As a characteristic of OS, ROS comprise a number of free radicals or reactive nonradical species, including singlet oxygen, hydrogen peroxide, superoxide anion, and so on [16], which are dramatically elevated in patients with SKCM [17]. Excessive ROS presence leads to genotoxicity and DNA damage when scavenging potential is outrun [18, 19], eventually inducing various genomic mutations and initiating tumorigenesis [20, 21]. In melanocytes, ROS are mainly derived from melanosomes, NADPH oxidase (NOX) family of enzymes, mitochondria, nitric oxide synthase (NOS) activity, and several arachidonic acid oxygenase activities [22]. Interestingly, the relationship between melanin and ROS appears to be a double-edged sword. On the one hand, melanin can protect melanocytes from UV radiation, but its synthesis in epidermal melanocytes results in additional intracellular ROS generation, which is also one of the reasons that melanocytes are particularly susceptible to OS [23, 24]. Together, these studies indicate that OS is closely associated with the progression of SKCM. Nevertheless, the value of OS genes in SKCM prognosis prediction remains largely unclarified and its underlying mechanisms require further validation.

High-throughput technologies have been widely applied to illustrate connections between hub genes and SKCM [25, 26]. A weighted gene co-expression network (WGCN) is frequently constructed to describe gene association patterns between different phenotypic traits [27]. SKCM data can thus be used to construct a powerful scale-free network to identify hub biomarkers or therapeutic targets for prognosis evaluation and clinical diagnosis [28, 29], with virtual applications in pathological or biological process identification (such as tumor progression), as well as brain imaging data and genetics analyses [30]. In the present study, we identified key modules and hub OS genes associated with the prognosis of and metastasis in patients with SKCM using WGCN analysis (WGCNA) to construct a risk model. We also applied a series of systematic analyses to identify the function and clinical significance of each identified gene in patients with SKCM, which may provide potential prognostic and diagnostic biomarkers for SKCM.

Data acquisition

An RNA-sequencing dataset containing 471 SKCM samples and 1 normal skin tissue sample was downloaded from The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov); meanwhile, its corresponding clinical and prognostic parameters were acquired from the University of California Santa Cruz Xena (UCSC Xena; http://xena.ucsc.edu/) [31]. Transcriptome data of 812 whole skin samples were retrieved from the Genome Tissue Expression (GTEx) database (https://gtexportal.org/home/datasets) [32, 33]. Additionally, gene expression profiles and clinical information of 214 patients with SKCM from the Gene Expression Omnibus (GEO) GSE65904 corhort (https://www.ncbi.nlm.nih.gov/geo/) were downloaded as a validation cohort [34]. All gene expression files were log2-transformed, and profiles from different platforms were normalized using the R package “sva” to remove batch effects reported by previous studies [35, 36]. Furthermore, 1399 protein domains of OS were downloaded from the GeneCards database (https://www.genecards.org) with relevance score ≥ 7, and subsequently applied for further exploration.

Construction of WGCN and hub gene identification

In this study, a WGCN was constructed to explore the relationship between expression modules and numerous clinicopathological features (overall survival time, overall survival state, age, gender, tumor stage, and metastasis). Expression profiles of OS genes in TCGA cohort were utilized to create a co-expression network with R package “WGCNA” [37, 38], and pairwise Pearson’s correlation coefficients were calculated between genes. Two parameters [module eigengenes (MEs) and gene significance (GS)] were identified to reveal modules most relevant to SKCM overall survival and metastasis. Hub genes comprising highly interconnected nodes within the module were regarded as functionally significant [39]. Thus, after choosing a significant module, genes with high module membership (MM > 0.6 and GS > 0.2) was selected. Transcriptional expression levels were then compared in normal tissues and SKCM samples; significantly differentially expressed genes were defined as candidate OS genes.

Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses

GO enrichment and KEGG pathway analyses were applied to systematically investigate the biological functions of OS genes in key modules using Database for Annotation, Visualization, and Integrated Discovery (DAVID) version 6.8 [40]. GO analysis comprehensively comprised three terms: biological process (BP), cellular component (CC), and molecular function (MF). Both p and FDR values < 0.05 were considered statistically significant.

Prognostic model construction and efficacy evaluation

All candidate OS genes were subjected to univariate Cox regression analysis with R package “survival” to explore the relationship between each gene and patient overall survival; genes with P < 0.05 were identified as prognosis-related OS genes. Thereafter, these genes were integrated to construct a multivariate Cox proportional hazards regression model, categorizing patients with SKCM into high or low risk groups. The formula of each sample’s risk score was calculated as follows:

where expgenei represents the relative expression value of OS gene i, and β represents the regression coefficient.

The Kaplan–Meier method and log-rank test using R/Bioconductor package “survival” were further conducted to compare outcomes between two risk subgroups. Additionally, the R package “survivalROC” was used to validate the predictive accuracy of the gene signature [41]. Univariate and multivariate Cox regression analyses were also conducted to evaluate the relationship between clinical characteristics and risk score. Finally, a nomogram incorporated with calibration plots was constructed to forecast the clinical outcome of patients with SKCM using R package “rms” [42].

Hub genes evaluation

To clarify the differential expression of the hub OS genes at a protein level, data in the Human Protein Atlas (HPA) online database (http://www.proteinatlas.org/) was used to elucidate the differences between normal and SKCM tumor tissues [43]. Expression of these OS genes in SKCM was also verified in TCGA dataset. Furthermore, to calculate the prognostic value of each gene for SKCM, the Kaplan-Meier method was employed on TCGA cohort, and the relationship between hub gene expression and different clinicopathological features of SKCM was systematically compared, especially connections with the cancer’s metastatic abilities. Nomogram and calibration plots based on the expression of hub genes were also constructed to forecast the clinical outcome of patients with SKCM.

Gene set enrichment analysis (GSEA)

To explore the potential functional roles that were significantly enriched in our constructed prognostic gene signature, GSEA was applied to a curated gene set “c2.cp.kegg.v7.1.symbols.gmt” and a GO gene set “c5.all.v7.1.symbols.gmt.” A normalized p value < 0.05 and false discovery rate (FDR) q value < 0.05 were considered statistically significant.

Identification of hub module by constructing WGCN

Bioinformatics analysis of publicly available datasets was performed according to the flow chart shown in Fig. 1. First, WGCNA of 471 SKCM samples with corresponding clinical data from TCGA cohort was carried out on 1399 extracted OS genes. As shown in Fig. 2A, six clinical characteristics of patients with SKCM, including overall survival state, overall survival time, age, gender, tumor stage, and metastasis, were included for analysis. To construct a scale-free network, β = 3 (scale-free R2 = 0.89) was selected as the soft threshold (Fig. 2B), and a total of four co‑expressed modules were identified (Fig. 3A–B). Subsequently, each module was assigned a different color to identify the module most associated with overall survival of SKCM. Among all modules, the turquoise module was specifically positively associated with tumor survival (p < 0.05). Interestingly, the 695 genes in the turquoise module also had the closest associations with SKCM metastasis (p < 0.05, Fig. 3C–D). Thus, the turquoise module was identified as the key module of interest in TCGA cohort.

Functional enrichment analysis

GO enrichment analysis suggested that the 695 OS genes were mainly enriched in the BP category, related to response to oxidative stress, oxygen levels, nutrient levels, hypoxia, and neuron death (Fig. 4A). In terms of MF, the key module was notably enriched in protein serine/threonine kinase activity, coenzyme binding, and oxidoreductase activity acting on NAD(P)H. With regard to the CC category, the key module was significantly enriched in mitochondrial matrix, neuronal cell body, and axon part. Gene symbols and their interaction with the enriched functions in GO are shown in Fig. 4B. In addition, the results of KEGG enrichment analysis indicated that the 695 OS genes in the turquoise module were mainly enriched in pathways of Hepatitis B, pancreatic cancer, fluid shear stress and atherosclerosis, chronic myeloid leukemia, and Kaposi sarcoma-associated herpesvirus infection (Fig. 4C). OS genes enriched in these pathways are displayed in Fig. 4D.

Screening of prognosis-related OS genes and construction of a genetic risk score model for patients with SKCM

Twenty-six highly relevant OS genes were screened from the turquoise module (Fig. 5A–B). Subsequently, transcription expression levels of these genes were compared between SKCM samples and normal tissues (Fig. 5C). Ultimately, 23 significantly differentially expressed OS genes were identified as candidate SKCM metastasis- and prognosis-associated genes. These genes further underwent univariate Cox regression analysis to identify nine OS genes with p < 0.05 (Fig. 6A). Thereafter, a multivariate Cox proportional hazards regression model was constructed (Fig. 6B), and four OS genes (AKAP9, VPS13C, ACSL4, and HMOX2) were ultimately selected for calculation of the genetic risk score. Meanwhile, all patients with SKCM in TCGA or GSE65904 cohorts were separated into low- and high-risk groups according to the median risk score (Fig. 6C–D).

Evaluating prognostic value of hub OS genes in patients with SKCM

First, we extracted the expression value of each key gene from TCGA cohort and drew a violin plot and heatmap. As shown in Fig. 7A–B, ACSL4 and HMOX2 were significantly overexpressed in SKCM samples, while expression of AKAP9 and VPS13C was decreased in comparison to that in normal tissues. Similar results were validated by analyzing the protein expression level of the key OS genes in accordance with the immunohistochemistry results from the HPA database (Fig. 7C–F). Subsequently, Kaplan-Meier analysis was also carried out to investigate the prognostic value of the four OS genes in SKCM. As shown in Fig. 8A–D, increased overall survival was significantly associated with elevated expression of AKAP9 (p = 0.004), VPS13C (p = 1.43e-04), and ACSL4 (p = 0.042); however, HMOX2 (p = 6.374e-05) expression had the opposite effect on the prognosis of patients with SKCM. These results indicated that our four identified OS genes were significantly correlated with the prognosis of SKCM.

Validating the relationship between OS genes and clinicopathological features of SKCM

While investigating the association between expression of hub OS genes and SKCM clinical features, we discovered that patients with metastatic melanoma had significantly increased expression of AKAP9, VPS13C, and ACSL4 (Fig. 8I–K), while HMOX2 demonstrated higher expression in patients with primary melanoma (Fig. 8L). Therefore, we further investigated connections between AKAP9, VPS13C, ACSL4, and HMOX2 expression and cancer metastasis. Results indicated that the four hub genes were all significantly related to metastasis (p < 0.05; Fig. 8E–H); AKAP9, VPS13C, and ACSL4 were positively associated with metastatic ability, while HMOX2 was negatively associated with cancer metastasis. Additionally, T stage of SKCM was also significantly associated with VPS13C, ACSL4, and HMOX2 transcription levels (p < 0.05; Fig. 8M–O), and patients with T3 or T4 stages had significantly enhanced HMOX2 expression and decreased VPS13C and ACSL4 expression. These results suggested that our identified hub genes played a vital role in metastasis and progression of SKCM. Therefore, we considered whether these was the reason that our identified genes controlled the overall survival of patients with SKCM and investigated the relationship between clinical features and patient survival. As shown in Fig. 8P–S, overall survival was significantly associated with patient age, T stage, TNM stage, and metastatic ability of the cancer (p < 0.01). Intriguingly, patients with primary melanoma or higher tumor stage had significantly worse outcomes, which suggested that AKAP9, VPS13C, ACSL4, and HMOX2 might play a vital role in the overall survival of SKCM through mediating cancer growth and metastasis.

Associations between prognostic risk score and clinical characteristics of SKCM

As per the survival analysis indicated in Fig. 9A, overall survival of patients with SKCM was significantly decreased with an increased risk score in TCGA cohort. In addition, ROC analysis indicated that our prediction model was more credible than the clinicopathological characteristics at both 1, 3, and 5 years in TCGA cohort (Fig. 9C). Similar conclusions were also validated in the GSE65904 cohort (Fig. 9B, 9D), demonstrating the moderate specificity and sensitivity of our prognostic model. Further, univariate and multivariate Cox regression analysis determined that risk score could be identified as an independent prognostic feature associated with SKCM prognosis (Fig. 10A–B). Moreover, we evaluated the connection between risk score and each clinicopathological characteristic, revealing that patients in T0-2 stage or those with metastatic melanoma were significantly associated with lower risk scores (Fig. 10E–F). Meanwhile, risk score was also significantly associated with the age of patients with SKCM (Fig. 10D). Additionally, we discovered that risk score was negatively associated with cancer metastasis (Fig. 10G), indicating that our constructed risk model was significantly associated with SKCM prognosis possibly by predicting the cancer’s metastatic ability. A heatmap revealed the expression of the four hub OS genes in the high- and low-risk groups (Fig. 10C), and their expression level in the high-risk group was remarkably consistent with their respective prognostic value in patients with SKCM. Meanwhile, the two risk groups differed significantly with respect to T stage, TNM stage, and metastasis in TCGA cohort. A nomogram plot based on risk score and other clinical characteristics was modeled to calculate the clinical outcomes of patients with SKCM (Fig. 11A); calibration plots at 3 and 5 years demonstrated good conformity with our prognostic model (Fig. 11B). A nomogram plot was also developed based on expression of the four hub OS genes, allowing us to calculate the survival probabilities of each patient with SKCM at 1, 3, and 5 years (Fig. 11C). Calibration plots also indicated good conformity between the predicted and observed outcomes at 3 and 5 years (Fig. 11D). These results indicated that our prognostic model showed great promise for predicting SKCM outcomes and clinical features.

GSEA analysis

GSEA analysis indicated that our selected four hub genes were significantly associated with several GO functions, especially death in response to oxidative stress, Golgi vesicle transport, proteasomal protein catabolic process, protein modification by small protein removal, regulation of cellular protein catabolic process, regulation of the ERAD pathway, regulation of intracellular protein transport, regulation of response to endoplasmic reticulum stress, ubiquitin like protein ligase binding, and vesicle targeting (Fig. 12A). In KEGG pathway analysis, these genes were also significantly associated with apoptosis, chronic myeloid leukemia, the ERBB signaling pathway, long term potentiation, the neurotrophin signaling pathway, pancreatic cancer, regulation of actin cytoskeleton, renal cell carcinoma, the Toll-like receptor signaling pathway, and ubiquitin mediated proteolysis (Fig. 12B).

Melanoma is one of the most common skin malignancies, and poor patient prognosis is significantly associated with its high metastatic capacity [44]. Although many novel diagnostic techniques and molecular biomarkers have been discovered in recent years, they have not sufficiently improved early diagnosis and prognosis of SKCM [11]. Therefore, it is imperative to identify SKCM prognosis-related molecules and the exact mechanism of cancer metastasis. In the present study, we constructed a co-expression network of OS genes to identify candidate gene clusters and hub OS genes involved in the prognosis of SKCM. Interestingly, we found that genes in the key module (turquoise) were not only significantly associated with SKCM overall survival but were also associated with the metastatic phenotype. Thus, we speculated that the key module could be employed to predict the outcomes of SKCM. What’s more, function analysis indicated that module turquoise was mainly enriched in oxidative stress and oxygen levels, while pathway enrichment analysis suggested this module could possibly play a role in the progression of several tumors, including pancreatic cancer, colorectal cancer, prostate cancer, hepatocellular carcinoma, glioma, endometrial cancer, and especially melanoma.

In the key module, 26 OS genes with MM > 0.6 and GS > 0.2 were screened from 695 module genes. Among them, four differentially expressed prognosis-related OS genes (AKAP9, VPS13C, ACSL4, and HMOX2) were identified as hub genes for further exploration. mRNA and protein expression patterns of these four genes using the expression data from TCGA and HPA cohorts revealed that ACSL4 and HMOX2 were overexpressed, while AKAP9 and VPS13C were downregulated in SKCM tissues. Meanwhile, Kaplan-Meier analysis revealed that HMOX2 was negatively associated with the overall survival of patients with SKCM, while AKAP9, VPS13C, and ACSL4 were positively related to patient outcome. These results might correspond with the modulation effects of these genes in cancer metastasis and growth. Omar et al. determined that AKAP9 acts as a metastasis-promoting gene in SKCM [45]. In addition, as an ROS scavenger in cellular compartments, downregulation of HMOX2 was shown to significantly inhibit tumor growth in lung cancer in vivo [46]. These findings concurred with our results that patients with SKCM had poor prognosis if their cancers were in stage III-IV or in primary sites. Meanwhile, downregulation of HMOX2 or overexpression of AKAP9, VPS13C, and ACSL4 was significantly associated with SKCM metastasis, growth inhibition, and higher overall survival, which was also associated with the prognostic effects of SKCM clinicopathological features. This suggested that our four specific OS genes might provide valuable biomarkers for adjusting treatment strategies of patients with SKCM.

Moreover, to identify whether these specific OS genes could be used as a prognostic factor, we constructed a novel prognostic prediction model based on the four hub genes. To our knowledge, this is the first OS-associated risk model for prognostication. Univariate and multivariate Cox regression analyses revealed that our risk model was an independent prognostic factor and demonstrated reliable prognostic value for SKCM. In addition, survival and ROC analyses confirmed the advantage of its biological implications for predicting SKCM prognosis, and its improved predictive accuracy compared with clinicopathological criteria. Our nomogram analysis confirmed the credibility of the genetic risk signature in predicting the overall survival of patients with SKCM. Taken together, our explorations demonstrate the prognostic value of an OS gene-associated risk model for patients with SKCM and suggest a novel direction for further research.

OS plays a critical role in various stages of carcinogenesis and cancer progression [47, 48]. SKCM research has also recently identified that OS may be closely associated with metastasis [23]. Considering the role of our four hub genes in cancer metastasis and growth, we assessed the relationship between risk score and SKCM clinical factors and found that our risk model was closely associated with SKCM metastasis and T stage of cancer. Therefore, our risk model was not only highly correlated with SKCM prognosis but may also provide potential biomarkers for predicting tumor progression.

Nonetheless, there are some limitations in this study. First, this study was designed as a retrospective analysis; more prospective research should be performed to verify our results. Second, our results lack in vitro or in vivo exploration to confirm the reliability of our mechanism analysis. Therefore, we need to conduct several experiments to prove the mechanistic connections between these genes and SKCM progression.

In conclusion, after constructing a co-expression network and performing bioinformatics analyses, we identified four prognosis-associated OS genes related to metastasis in patients with SKCM. We also successfully constructed a prognostic model with powerful predictive effects. To our knowledge, this is the first OS-associated model to predict the prognosis of malignancies. Our study provides a novel research approach for studying the pathogenesis of SKCM, and greatly contributes to elucidating the metastatic mechanism of SKCM.

Authors’ contributions

Conceptualization, JM-Z; Methodology, XP-W, ZY-W; Software, XP-W; Validation, ZY-W; Investigation, JM-Z, XP-W; Writing-Original Draft Preparation, XP-W; Writing-Review & Editing, ZY-W; Supervision, JM-Z; Project Administration, JM-Z. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge and appreciate our colleagues for their valuable efforts and comments on this paper.

Competing interests

All authors declare that they have no competing interests.

Availability of data and materials

The datasets analysed during the current study are available were sourced from the publicly available TCGA, GTEx, GEO, and HPA database (https://portal.gdc.cancer.gov, https://gtexportal.org/home/datasets, https://www.ncbi.nlm.nih.gov/geo/, http://www.proteinatlas.org/).

Animal ethics Committee

Not applicable.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

Not applicable.

Ekwueme DU, Guy GP Jr, Li C, Rim SH, Parelkar P, Chen SC. The health burden and economic costs of cutaneous melanoma mortality by race/ethnicity-United States, 2000 to 2006. J Am Acad Dermatol. 2011;65(5 Suppl 1):133–43.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Board PDQATE. Melanoma Treatment (PDQ®): Health Professional Version. In: PDQ Cancer Information Summaries. Bethesda (MD): National Cancer Institute (US); 2002.
Kanavy HE, Gerstenblith MR. Ultraviolet radiation and melanoma. Semin Cutan Med Surg. 2011;30(4):222–8.
Finn L, Markovic SN, Joseph RW. Therapy for metastatic melanoma: the past, present, and future. BMC Med. 2012;10:23.
Liu-Smith F, Jia J, Zheng Y. UV-Induced Molecular Signaling Differences in Melanoma and Non-melanoma Skin Cancer. Adv Exp Med Biol. 2017;996:27–40.
Gershenwald JE, Scolyer RA, Hess KR, Sondak VK, Long GV, Ross MI, Lazar AJ, Faries MB, Kirkwood JM, McArthur GA, et al. Melanoma staging: Evidence-based changes in the American Joint Committee on Cancer eighth edition cancer staging manual. CA Cancer J Clin. 2017;67(6):472–92.
Hamm C, Verma S, Petrella T, Bak K, Charette M. Biochemotherapy for the treatment of metastatic malignant melanoma: a systematic review. Cancer Treat Rev. 2008;34(2):145–56.
Tsai J, Lee JT, Wang W, Zhang J, Cho H, Mamo S, Bremer R, Gillette S, Kong J, Haass NK, et al. Discovery of a selective inhibitor of oncogenic B-Raf kinase with potent antimelanoma activity. Proc Natl Acad Sci U S A. 2008;105(8):3041–6.
Leach DR, Krummel MF, Allison JP. Enhancement of antitumor immunity by CTLA-4 blockade. Science. 1996;271(5256):1734–6.
Wang LX, Li Y, Chen GZ. Network-based co-expression analysis for exploring the potential diagnostic biomarkers of metastatic melanoma. PLoS One. 2018;13(1):e0190447.
Wang HZ, Wang F, Chen PF, Zhang M, Yu MX, Wang HL, Zhao Q, Liu J. Coexpression network analysis identified that plakophilin 1 is associated with the metastasis in human melanoma. Biomed Pharmacother. 2019;111:1234–42.
Zhou F, Pan Y, Wei Y, Zhang R, Bai G, Shen Q, Meng S, Le XF, Andreeff M, Claret FX. Jab1/Csn5-Thioredoxin Signaling in Relapsed Acute Monocytic Leukemia under Oxidative Stress. Clin Cancer Res. 2017;23(15):4450–61.
Brown JM, Wilson WR. Exploiting tumour hypoxia in cancer treatment. Nat Rev Cancer. 2004;4(6):437–47.
Kangari P, Zarnoosheh Farahany T, Golchin A, Ebadollahzadeh S, Salmaninejad A, Mahboob SA, Nourazarian A. Enzymatic Antioxidant and Lipid Peroxidation Evaluation in the Newly Diagnosed Breast Cancer Patients in Iran. Asian Pac J Cancer Prev. 2018;19(12):3511–5.
Lü JM, Lin PH, Yao Q, Chen C. Chemical and molecular mechanisms of antioxidants: experimental approaches and model systems. J Cell Mol Med. 2010;14(4):840–60.
Liu-Smith F, Dellinger R, Meyskens FL Jr. Updates of reactive oxygen species in melanoma etiology and progression. Arch Biochem Biophys. 2014;563:51–5.
Wang JY, Liu GZ, Wilmott JS, La T, Feng YC, Yari H, Yan XG, Thorne RF, Scolyer RA, Zhang XD, et al. Skp2-Mediated Stabilization of MTH1 Promotes Survival of Melanoma Cells upon Oxidative Stress. Cancer Res. 2017;77(22):6226–39.
Zhou F, Shen Q, Claret FX. Novel roles of reactive oxygen species in the pathogenesis of acute myeloid leukemia. J Leukoc Biol. 2013;94(3):423–9.
Oates JC, Gilkeson GS. The biology of nitric oxide and other reactive intermediates in systemic lupus erythematosus. Clin Immunol. 2006;121(3):243–50.
Smith J, Tho LM, Xu N, Gillespie DA. The ATM-Chk2 and ATR-Chk1 pathways in DNA damage signaling and cancer. Adv Cancer Res. 2010;108:73–112.
Denat L, Kadekaro AL, Marrot L, Leachman SA, Abdel-Malek ZA. Melanocytes as instigators and victims of oxidative stress. J Invest Dermatol. 2014;134(6):1512–8.
Obrador E, Liu-Smith F, Dellinger RW, Salvador R, Meyskens FL, Estrela JM. Oxidative stress and antioxidants in the pathophysiology of malignant melanoma. Biol Chem. 2019;400(5):589–612.
Jenkins NC, Grossman D. Role of melanin in melanocyte dysregulation of reactive oxygen species. Biomed Res Int. 2013;2013:908797.
Qiu T, Wang H, Wang Y, Zhang Y, Hui Q, Tao K. Identification of genes associated with melanoma metastasis. Kaohsiung J Med Sci. 2015;31(11):553–61.
Haqq C, Nosrati M, Sudilovsky D, Crothers J, Khodabakhsh D, Pulliam BL, Federman S, Miller JR 3rd, Allen RE, Singer MI, et al. The gene expression signatures of melanoma progression. Proc Natl Acad Sci U S A. 2005;102(17):6092–7.
Udyavar AR, Hoeksema MD, Clark JE, Zou Y, Tang Z, Li Z, Li M, Chen H, Statnikov A, Shyr Y, et al. Co-expression network analysis identifies Spleen Tyrosine Kinase (SYK) as a candidate oncogenic driver in a subset of small-cell lung cancer. BMC Syst Biol. 2013;7(Suppl 5):1. ., Suppl 5.
Yang L, Xu Y, Yan Y, Luo P, Chen S, Zheng B, Yan W, Chen Y, Wang C: Common Nevus and Skin Cutaneous Melanoma: Prognostic Genes Identified by Gene Co-Expression Network Analysis. Genes (Basel) 2019, 10(10).
Wan Q, Tang J, Han Y, Wang D. Co-expression modules construction by WGCNA and identify potential prognostic markers of uveal melanoma. Exp Eye Res. 2018;166:13–20.
Shi K, Bing ZT, Cao GQ, Guo L, Cao YN, Jiang HO, Zhang MX. Identify the signature genes for diagnose of uveal melanoma by weight gene co-expression network analysis. Int J Ophthalmol. 2015;8(2):269–74.
Liu J, Li R, Liao X, Hu B, Yu J. Comprehensive investigation of the clinical significance and molecular mechanisms of plasmacytoma variant translocation 1 in sarcoma using genome-wide RNA sequencing data. J Cancer. 2019;10(20):4961–77.
Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science 2015, 348(6235):648–660.
Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.
Huang R, Mao M, Lu Y, Yu Q, Liao L. A novel immune-related genes prognosis biomarker for melanoma: associated with tumor microenvironment. Aging. 2020;12(8):6966–80.
Zhang M, Wang X, Chen X, Zhang Q, Hong J. Novel Immune-Related Gene Signature for Risk Stratification and Prognosis of Survival in Lower-Grade Glioma. Front Genet. 2020;11:363.
Xiao Y, Zhu Z, Li J, Yao J, Jiang H, Ran R, Li X, Li Z. Expression and prognostic value of long non-coding RNA H19 in glioma via integrated bioinformatics analyses. Aging. 2020;12(4):3407–30.
Mason MJ, Fan G, Plath K, Zhou Q, Horvath S. Signed weighted gene co-expression network analysis of transcriptional regulation in murine embryonic stem cells. BMC Genom. 2009;10:327.
Horvath S, Dong J. Geometric interpretation of gene coexpression network analysis. PLoS Comput Biol. 2008;4(8):e1000117.
Albert R. Scale-free networks in cell biology. J Cell Sci. 2005;118(Pt 21):4947–57.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Heagerty PJ, Zheng Y. Survival model predictive accuracy and ROC curves. Biometrics. 2005;61(1):92–105.
Gu HY, Zhang C, Guo J, Yang M, Zhong HC, Jin W, Liu Y, Gao LP, Wei RX. Risk score based on expression of five novel genes predicts survival in soft tissue sarcoma. Aging. 2020;12(4):3807–27.
Thul PJ, Åkesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, Alm T, Asplund A, Björk L, Breckels LM, et al: A subcellular map of the human proteome. Science 2017, 356(6340).
Palmieri G, Capone M, Ascierto ML, Gentilcore G, Stroncek DF, Casula M, Sini MC, Palla M, Mozzillo N, Ascierto PA. Main roads to melanoma. J Transl Med. 2009;7:86.
Kabbarah O, Nogueira C, Feng B, Nazarian RM, Bosenberg M, Wu M, Scott KL, Kwong LN, Xiao Y, Cordon-Cardo C, et al. Integrative genome comparison of primary and metastatic melanomas. PLoS One. 2010;5(5):e10770.
Kim JJ, Lee YA, Su D, Lee J, Park SJ, Kim B, Jane Lee JH, Liu X, Kim SS, Bae MA, et al. A Near-Infrared Probe Tracks and Treats Lung Tumor Initiating Cells by Targeting HMOX2. J Am Chem Soc. 2019;141(37):14673–86.
Reuter S, Gupta SC, Chaturvedi MM, Aggarwal BB. Oxidative stress, inflammation, and cancer: how are they linked? Free Radic Biol Med. 2010;49(11):1603–16.
Hecht F, Pessoa CF, Gentile LB, Rosenthal D, Carvalho DP, Fortunato RS. The role of oxidative stress on breast cancer development and therapy. Tumour Biol. 2016;37(4):4281–91.

Download PDF

Journal Publication

published 26 Feb, 2021

Read the published version in Medicine →

Version 1

posted

You are reading this latest preprint version

Novel Oxidative Stress-related Prognostic Biomarkers for Melanoma Associated With Tumor Metastasis

Status:

Journal Publication

Version 1

Abstract

Background

Methods

Results

Conclusions

Figures

Background

Methods

Data acquisition

Construction of WGCN and hub gene identification

Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses

Prognostic model construction and efficacy evaluation

Hub genes evaluation

Gene set enrichment analysis (GSEA)

Results

Identification of hub module by constructing WGCN

Functional enrichment analysis

Evaluating prognostic value of hub OS genes in patients with SKCM

Validating the relationship between OS genes and clinicopathological features of SKCM

Associations between prognostic risk score and clinical characteristics of SKCM

GSEA analysis

Discussion

Conclusions

Declarations

References

Status:

Journal Publication

Version 1