AUTHOR=Shan Guangyao , Zhang Huan , Bi Guoshu , Bian Yunyi , Liang Jiaqi , Valeria Besskaya , Zeng Dejun , Yao Guangyu , Zhan Cheng , Fan Hong 

TITLE=Multi-Omics Analysis of Cancer Cell Lines with High/Low Ferroptosis Scores and Development of a Ferroptosis-Related Model for Multiple Cancer Types

JOURNAL=Frontiers in Cell and Developmental Biology

VOLUME=9

YEAR=2021

URL=https://www.frontiersin.org/journals/cell-and-developmental-biology/articles/10.3389/fcell.2021.794475

DOI=10.3389/fcell.2021.794475

ISSN=2296-634X

ABSTRACT=<p><bold>Background:</bold> Ferroptosis is a newly identified regulated cell death characterized by iron-dependent lipid peroxidation and subsequent membrane oxidative damage, which has been implicated in multiple types of cancers. The multi-omics differences between cancer cell lines with high/low ferroptosis scores remain to be elucidated.</p><p><bold>Methods and Materials:</bold> We used RNA-seq gene expression, gene mutation, miRNA expression, metabolites, copy number variation, and drug sensitivity data of cancer cell lines from DEPMAP to detect multi-omics differences associated with ferroptosis. Based on the gene expression data of cancer cell lines, we performed LASSO-Logistic regression analysis to build a ferroptosis-related model. Lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), esophageal cancer (ESCA), bladder cancer (BLCA), cervical cancer (CESC), and head and neck cancer (HNSC) patients from the TCGA database were used as validation cohorts to test the efficacy of this model.</p><p><bold>Results:</bold> After stratifying the cancer cell lines into high score (HS) and low score (LS) groups according to the median of ferroptosis scores generated by gene set variation analysis, we found that IC50 of 66 agents such as oxaliplatin (<italic>p &lt; 0.001</italic>) were significantly different, among which 65 were higher in the HS group. 851 genes such as KEAP1 and NRAS were differentially muted between the two groups. Differentially expressed genes, miRNAs and metabolites were also detected—multiple items such as IL17F (logFC = 6.58, <italic>p &lt; 0.001</italic>) differed between the two groups. Unlike the TCGA data generated by bulk RNA-seq, the gene expression data in DEPMAP are from pure cancer cells, so it could better reflect the traits of tumors in cancer patients. Thus, we built a 15-signature model (AUC = 0.878) based on the gene expression data of cancer cell lines. The validation cohorts demonstrated a higher mutational rate of NFE2L2 and higher expression levels of 12 ferroptosis-related genes in HS groups.</p><p><bold>Conclusion:</bold> This article systemically analyzed multi-omics differences between cancer cell lines with high/low ferroptosis scores and a ferroptosis-related model was developed for multiple cancer types. Our findings could improve our understanding of the role of ferroptosis in cancer and provide new insight into treatment for malignant tumors.</p>