frontiers
Check for updates
OPEN ACCESS
EDITED BY William C Cho, QEH, Hong Kong SAR, China
REVIEWED BY Md Khadem Ali, Stanford University, United States Emil Bulatov, Kazan Federal University, Russia
*CORRESPONDENCE Cong Ye, sphyc1989@163.com
“These authors have contributed equally to this work and share first authorship
SPECIALTY SECTION This article was submitted to Cancer Genetics and Oncogenomics, a section of the journal Frontiers in Genetics
RECEIVED 21 May 2022 ACCEPTED 31 October 2022 PUBLISHED 15 November 2022
CITATION Wang W, Chang G, Zhuo R and Ye C (2022), Identification of a ferroptosis- related long noncoding RNA signature with a prognostic value in adrenocortical carcinoma. Front. Genet. 13:949457. doi: 10.3389/fgene.2022.949457
COPYRIGHT
@ 2022 Wang, Chang, Zhuo and Ye. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Identification of a ferroptosis-related long noncoding RNA signature with a prognostic value in adrenocortical carcinoma
Weixi Wang1t, Guilin Chang11, Ran Zhuo2 and Cong Ye3*
1Department of Geriatrics, Zhongshan Hospital, Fudan University, Shanghai, China, 2Department of Urology, Zhongshan Hospital, Fudan University, Shanghai, China, 3Department of Thoracic Surgery, Shanghai Pulmonary Hospital, Tongji University, Shanghai, China
Background: Adrenocortical carcinoma (ACC) is an uncommon endocrine malignancy associated with poor clinical outcome. As a novel form of cell death, ferroptosis is reliant on the accumulation of iron and reactive oxygen species and is involved in the pathogenesis of various tumors, including ACC. Our study aimed to identify and characterize the prognostic ferroptosis-related lncRNA signature (FerRLSig) in ACC.
Methods: A regulatory network of ferroptosis-related lncRNAs (FerRLs) and mRNAs was constructed based on The Cancer Genome Atlas (TCGA). Univariate and multivariate Cox regression assays were performed to construct the FerRLSig.
Results: Twenty-four FerRLs were identified in the prognostic model, and the high-risk FerRLSig was related to the worse overall survival (OS) in ACC [hazard ratio (HR): 1.936 (1.484-2.526), p < 0.001]. The area under the curve (AUC) value of the FerRLSig was 0.936 according to the receiver operating characteristic (ROC) analyses, superior to other traditional clinicopathological features, further supported the utility in prognosis prediction of ACC. We further established a prognostic nomogram combining clinical factors with the FerRLSig, which showed favorable efficacy for survival prediction. Next, gene set enrichment analysis (GSEA) revealed that gene sets were involved in many immune regulatory biological processes related to malignancies. T-cell function of type II INF response and the immune checkpoints, including CD40, CD276, IDO2, NRP1, and CD80, were expressed with a significant difference between the low- and high-risk groups.
Conclusion: This study offered new insights into the pathogenesis of ACC. The novel FerRLSig could be useful in predicting survival and may provide information of immunological research and treatment for ACC patients.
KEYWORDS
adrenocortical carcinoma (ACC), ferroptosis, lncRNA, immune infiltration, data mining
Introduction
Adrenocortical carcinoma (ACC) is an aggressive endocrine malignancy with a rare incidence but high mortality (Else et al., 2014; Zheng et al., 2016). Although patients could benefit from complete surgical resection or mitotane treatment, the prognosis of ACC is dismal with the 5-year survival of 20%-59% (Fassnacht et al., 2011; Zheng et al., 2016). In addition, the prognosis of ACC varies with the patient’s age, surgical scope, mitotic intensity, and secretion of hormones. Due to the tumor heterogeneity of ACC, the current tumor, lymph node, and metastasis (TNM) staging system is not reliable in clinical outcome prediction. It is challenging to make a prediction of the survival in ACC patients owing to the low incidence rate, diverse pathogenic determinants, high heterogeneity, and dismal prognosis. Therefore, it is particularly critical to uncover new and reliable prognostic biological markers for patients with ACC.
Past decades have seen a growing number of research studies to investigate the function of ferroptosis in malignancies. As a novel form of cell death, ferroptosis is reliant on the accumulation of intracellular iron and reactive oxygen species, distinguished from other forms of cell death such as apoptosis and autophagy (Dixon et al., 2012; Huang et al., 2020). Imbalance of iron metabolism could facilitate tumor growth and also act as a risk factor for tumorigenesis. Compared with normal cells, malignant cells excessively depend on iron for cell proliferation, that is, iron addiction (Manz et al., 2016). One recent study uncovered that ACC is remarkably sensitive to ferroptosis in human, and compared to mitotane, the ferroptosis inducer may be more specific, more effective, and less toxic for the treatment of ACC patients (Belavgeni et al., 2019). Weigand et al. (2020) demonstrated that adrenocortical cells are extraordinary sensitive to ferroptosis owing to the active steroid synthesis pathway. In addition, Chen and colleagues identified a ferroptosis-related signature which could be useful in prognosis prediction and immunotherapy or targeted therapy screening (Chen et al., 2021). Therefore, it is of vital importance to further investigate the changing profiles of ferroptosis-related genes and the underlying mechanisms in ACC.
As a subset of RNA molecules, a long noncoding RNA (lncRNA) has the length of over 200 nt and is involved in gene regulation (Zhang et al., 2021a). In addition, lncRNA is involved in various biological process and plays an important role in the development and progression of malignancies (Gupta et al., 2010), including ACC (Buishand et al., 2020). At present, few studies have focused on the functions of ferroptosis-related lncRNAs (FerRLs) in tumors. Recently, Chao et al. reported that cytosolic lncRNA P53RRA suppressed tumor progression by promoting ferroptosis. The P53RRA-G3BP1 interaction in the cytoplasm leads to the
| Variable | Number of samples |
|---|---|
| Gender | |
| Male/Female | 32/60 |
| Age at diagnosis | |
| ≤65/>65 | 81/11 |
| Stage | |
| I/II/III/IV/NA | 9/44/19/18/2 |
| T | |
| T1/T2/T3/T4/NA | 9/49/11/21/2 |
| M | |
| M0/M1/NA | 69/21/2 |
| N | |
| N0/N1/NA | 80/10/2 |
accumulation of p53 in the nucleus, conferring to ferroptosis and cell-cycle arrest. In patients with breast cancers or lung cancers harboring the wild-type p53, the decreased expression of P53RRA was significantly related to a poor prognosis (Mao et al., 2018), and p53 is considered as one of the most promising molecular targets for cancer therapy (Chasov et al., 2020). In another related study, Zhang et al. (2021b) demonstrated that the chronic cadmium exposure could significantly increase the expression of lncRNA OIP5-AS1 in prostate cancer. OIP5-AS1 regulates the expression of SLC7A11 (an alternative marker of ferroptosis) by sponging miR-128-3p, leading to enhanced ferroptosis and decreased cell viability. Similarly, LINC00336 serves as a competing endogenous sponge of miR6852 to inhibit ferroptosis in lung cancer (Wang et al., 2019a). However, no research has been performed to identify the function of FerRLs in predicting the overall survival (OS) for ACC patients. Thus, we constructed a prognostic ferroptosis-related lncRNAs signature (FerRLSig) according to The Cancer Genome Atlas (TCGA) database. Additionally, functions and biological pathways of ferroptosis-related genes and the N6-methyladenosine (m6A) mRNA status and immune checkpoint in the survival prediction of ACC patients were further explored.
Materials and methods
Data collection
On 5 February 2022, we extracted the RNA-sequence data of 92 patients from the TCGA-ACC data set (https://portal.gdc. cancer.gov/). Table 1 lists the clinical features, and the identified lncRNAs were annotated using the GTF annotation files of human lncRNAs retrieved from the GENCODE (https://www. gencodegenes.org/). We obtained 382 ferroptosis-related genes
(n = 150 for the drivers, n = 109 for the suppressors, and n = 123 for the markers) from the FerrDb database, which is a data network with comprehensive and timely updates, covering the latest progress of ferroptosis-related genes along with their regulatory molecules and related diseases (Zhou and FerrDb, 2020) (http://www.zhounan.org/ferrdb, Supplementary Table S1). Information of the tumor grade in ACC was missing; thus, we collected clinicopathological data including gender, age, TNM stage, and survival information. Since the study was based on a public database, ethical approval is not required.
Identification of a prognostic FerRLSig
The relationship between the FerRLs and ferroptosis-related genes was assessed by Pearson correlation. Statistical significance was achieved at a correlation coefficient |R2|>0.8 and p < 0.001. Next, a FerRLSig was constructed based on the univariate and multivariate Cox regression analyses. Information of age, gender, stage, and risk score was included in the multivariate Cox regression analysis. The risk score was calculated as the sum of (coefficient of IncRNA1 x IncRNA1 expression) + (coefficient of IncRNA2 x IncRNA2 expression) + + (coefficient of IncRNAn x IncRNAn expression). Taking the median value as the cut-off point, the lncRNAs were divided into two groups: high-risk group and low-risk group. Differentially expressed genes (DEGs) between the low- and high-risk groups were identified and statistically significant was considered at a false discovery rate (FDR) < 0.05 and |log2FC|≥1. A Kaplan-Meier survival analysis was then performed to evaluate the clinical correlation of risk stratification and patient prognosis. The FerRLSig was present as a risk plot, including the distribution of risk scores and survival status in ACC patients. We further calculated the area under the curve (AUC) value in the receiver operating characteristic (ROC) curve to evaluate the prediction accuracy of OS in ACC patients. Sensitivity and specificity of the FerRLSig were compared with other clinical and pathological factors via ROC curves and decision curve analysis (DCA) (Vickers and Elkin, 2006). Cytoscape 3.9.0 was used to visualize the co-expression network of ferroptosis-related genes, lncRNAs, and signaling pathways.
Functional analyses of DEGs
Functions of DEGs were compared between the low-risk group and the high-risk group. We evaluated the associated biological functions using gene ontology (GO), including biological processes (BPs), molecular functions (MFs), and cellular components (CCs). Based on the data in the Kyoto Encyclopedia of Genes and Genomes (KEGG), we further analyzed the biological pathways using the ggplot2 package in R 4.1.2.
Gene set enrichment analysis and the predictive nomogram
To explore potential molecular mechanisms of the FerRLSig, we performed a gene set enrichment analysis (GSEA) to define the enriched terms in the KEGG. Significant gene sets were achieved at p-value <0.05 (GSEA 4.2.2). Statistical significance was expressed by the normalized enrichment score (NES) and FDR (Subramanian et al., 2005). A nomogram was developed by integrating the FerRLSig and clinical characteristics to predict the 1-, 3-, and 5-year OS in ACC.
Immune infiltration and gene expression analyses
Potential correlation of the FerRLSig and immune cells infiltration was analyzed using the CIBERSORT algorithms (Newman et al., 2015; Charoentong et al., 2017) and TIMER website (http://timer.comp-genomics.org/) (Li et al., 2017). The differences of immune responses between the two groups were illustrated on a cluster heatmap. Next, a single-sample gene set enrichment analysis (ssGSEA) was performed to analyze the scores of tumor-infiltrating immune cells along with their functions in TCGA. Furthermore, we retrieved the previous literature to extract potential immune checkpoints (Charoentong et al., 2017).
Statistical analysis
R software (version 4.1.2), along with its appropriate packages, was used to perform the statistical analysis. The Wilcoxon test and unpaired Student’s t-test were used to analyze the non-normally and normally distributed variables, respectively. Differentially expressed genes between the low- and high-risk groups were identified using the Benjamini-Hochberg method. The association between FerRLs and clinicopathological characteristics of ACC patients was assessed by logistic regression analysis, presented in a cluster heatmap. Statistical significance was achieved at p < 0.05.
Results
Identification of a prognostic FerRLSig in adrenocortical carcinoma
According to the Pearson correlation analysis, a total of 85 FerRLs were identified from the 382 ferroptosis-related genes in TCGA-ACC dataset (Supplementary Table S2). Subsequently, 32 significant FerRLs were screened by the univariate analysis and then included in the multivariate Cox analysis. Overall, 24 FerRLs (AC046143.1, AC099850.3,
LINC02268
RGS4
Leishmaniasis
MIR222HG
AC012073.1
NCF2
NOS2
SAT1
AC046143.1
LINC01503
PTGS2
PRKAR1B-AS1
NNMT
RRM2
U62317.1
Amoebiasis
PRDX1
CAV1
IL-17 signaling pathway
IL6
SLC7A5
BX470102
AC008105.2
AC099850.3
DGUOK-AS1
AC004816.2
CAPG
CXCL2
MUC1
ASTR
ZFP36
AL513477.2
STMN1
CDKN2A
Kaposi sarcoma-associated herpesvirus infection
TNF signaling pathway
AC091057.1
LINC01843
LINC00460
DPP4
GAPLINC
HNF4A
SLC2A3
AL391422.4
HIC1
LINC01614
AC123912.4
MYOSLID
DPP4-DT
AC112715.1
FIGURE 1 Regulatory network of the ferroptosis-related lncRNAs, mRNAs, and signaling pathways. Diamonds represent lncRNAs in red, rectangles indicate signaling pathways in blue, and ellipses represent mRNAs in green.
LINC01614, LASTR, DPP4-DT, AL391422.4, LINC01503, LINC01843, AC008105.2, AC004816.2, AC112715.1, AC012073.1, LINC00460, U62317.1, BX470102.1, LINC02268, GAPLINC, MYOSLID, AC123912.4, AC091057.1, DGUOK- AS1, PRKAR1B-AS1, MIR222HG, and AL513477.2) were identified to be independent prognostic factors for ACC (Supplementary Table S3). Annotations and potential mechanisms of these identified FerRLs are presented in Supplementary Table S4. Next, a prognostic FerRLSig was constructed based on the calculated risk scores.
The regulatory network of FerRLs, mRNAs, and pathways is shown in Figure 1, which indicated that these FerRLs regulate multiple ferroptosis-related pathways by targeting the mRNAs. For example, lncRNA BX470102.1 is involved in amoebiasis,
Kaposi sarcoma-associated herpesvirus infection, TNF signaling pathway, and IL-17 signaling pathway by targeting multiple target genes such as IL-6, CXCL2, and PTGS2. The target genes of LINC01503, U62317.1, AC008105.2, LASTR, and AC004816.2 are also involved in the aforementioned pathways, which have been proved to be related to ferroptosis in multiple tumors (Liu et al., 2021).
Enrichment analysis of the identified ferroptosis-related genes
Overall, 540 genes were identified to be differentially expressed between the low- and high-risk groups (n = 103 for
A
organelle fission nuclear division
B
chromosome segregation
Cell cycle
mitotic nuclear division
nuclear chromosome segregation
mitotic cell cycle phase transition
8
Oocyte meiosis
sister chromatid segregation -
mitotic sister chromatid segregation
Cellular senescence
spindle organization -
Count
Count
regulation of nuclear division
Progesterone-mediated oocyte maturation
10
chromosomal region
20
spindle
40
15
20
condensed chromosome
60
Hippo signaling pathway
chromosome, centromeric region
25
microtubule
condensed chromosome, centromeric region
8
qvalue
p53 signaling pathway
kinetochore
0.0125
qvalue
0.0100
0.04
mitotic spindle
spindle pole
Protein digestion and absorption-
0.0075
0.03
DNA replication preinitiation complex
0.0050
0.02
tubulin binding
0.0025
Fanconi anemia pathway
0.01
microtubule binding
extracellular matrix structural constituent
catalytic activity, acting on DNA
ECM-receptor interaction
extracellular matrix structural constituent conferring tensile strength
microtubule motor activity
₹
cyclin-dependent protein serine/threonine kinase regulator activity
Small cell lung cancer
single-stranded DNA helicase activity
proteoglycan binding
DNA replication
platelet-derived growth factor binding
.
0.04
0.08
0.12
0.16
0.03
0.06
0.09
0.12
GeneRatio
GeneRatio
A
Risk
High risk
Low risk
B
C
D
1.00
0.10-
40
0
· High- risk
Survival probability
0.75
Low-Risk
O.
Risk score
0.50
9
Net Benefit
0.05
Risk
0.8
Age
0.25
Sensitivity
Gender
p<0.001
-80
Stage
All
0.00
0.4
None
0
20
40
60
80
0
1
2
3
4
5
6
7
9
9
10
11
12
0.00
Time(years)
2
Patients (increasing risk socre)
Risk, AUC=0.936
High risk Low risk
Age, AUC=0.707
39
35
22
16
6
2
0
0
0
0
D
0
0
Gender, AUC=0.438
Dead
40
40
36
28
24
22
16
11
8
7
4
2
2
0
Stage, AUC=0.587
Survival time (years)
2
Alive.
0
1
2
3
4
5
8
7
8
9
10
11
12
Time(years)
0.0
0.2
0.4
0.6
0.8
1.0
0.00
0.05
0.10
0.15
Risk Threshold
3
1-Specificity
00
៛
CV
o
E
F
G
0
20
40
60
80
9
Patients (increasing risk socre)
pealue
Hazard ratio
pvalua
Hazard ratio
U.S
0000143.1
2
8
C293850.3
LINGOT 614
:
/qu
0.328
1.012(0.588-1/458)
Age
0.320
1.0140 656-1.0421
LASTR
Sensitivity
0.6
OPP4-OT
NL291422.4
LINGO:533
#
Gender
0.972
0.983(0.451-2.154)
Gender
C.578
077700 322-1.800]
LINCO: 843
A08581052
0.4
NB:43
-
**
0 CÔI
2 914(1.860-4.555)
Stage
0 450
1.2240 030-2 1751
2
Meyen
AUC at 1 years: 0.936
MAYA
062317.1
AUC at 2 years: 0.983
@X470102.1
makScore
00.001
1 875(1 527-2 553)
AUC at 3 years: 0.989
LING02266
H
15K50000
4001
1.939(1,494-2.526)
GAPLING
Q
1
2
3
3
0.0
0.5
MYOBUID
1.0
Hazard ratio
1.5
2.0
2.6
AG123312.4
Hazard ratio
0.0
0.2
0.4
0.6
0.8
1.0
0091057.1
DGJON-AS1
1-Specificity
FRKARIB-ASI
MIR222HG
AL513477.2
A
B
Points
0
10
20
30
40
50
60
70
80
90
100
N
N
-
NO 41 unknow
M
MO
Stage **
Gender
5
M1
Age
M **
risk
MO M. unknow
Age
AC046143.1
10
AC099850.3
Q
Gender
MALE
LINC01614
T1 12
·
LASTR
13 T4
-5
DPP4-DT
unknow
N
FEMALE
AL391422 4
Stage”
Stage I
N1
LINCO1503
Stage II
Stage II
T ***
T2,
LINCO1843
Stage I unknow
AC008106.2
Gendar
AC004816.2
Stage ***
FEMALE MALE
Stage I
Stage HI
AC112715.1
Ape
low
i
SyngeMl
AC012073.1
< 65
risk
LINC00460
>65
risk
·
·
U62317.1
high
Total points
high
BX470102.1
LINC02268
380
400
420
440
460
490
GAPLINC
480
500
MYOSLID
Pr( futime < 5)
0.731
AC123912.4
8.0020.08 0.998
AC091057.1
Pr( futime < 3
0.286
DGUOK-AS1
0.000.04 0.96
PRKAR18-AS1
Prí futime
0.0549
MIR222HG
6.006
0.4
AL513477 2
FIGURE 4 (A) A nomogram for the identified FerRLs and other clinicopathological factors. (B) Cluster heatmap for the association of clinicopathological manifestations and the prognostic FerRLSig.
the low-risk group and n = 437 for the high-risk group, Supplementary Table S5). BP participated in organelle fission, mitotic nuclear division, chromosome segregation, and mitotic cell cycle phase transition. CC were mainly involved in the chromosomal region, mitotic spindle, and DNA replication preinitiation complex pathways. MF mainly regulated the single-stranded DNA helicase activity, cyclin-dependent protein serine/threonine kinase regulator activity, tubulin binding, and extracellular matrix structural constituent (Figure 2A). The KEGG pathway analyses revealed that the genes in the high-risk group were particularly involved in oocyte meiosis, cell cycle, DNA replication, protein digestion and absorption, Hippo signaling pathway, small-cell lung cancer, cellular senescence, and p53 signaling pathway (Figure 2B; Supplementary Table S6).
Prognostic value of the FerRLSig
High-risk FerRLs were associated with a worse OS according to the Kaplan-Meier analyses (Figure 3A, p < 0.001). In addition, the AUC value for the FerRLSig was 0.936 according to the ROC curve, significantly higher than that for traditional clinicopathological factors such as gender, age, and tumor stage (Figures 3B,C), further supporting the potential superior performance of the FerRLSig in predicting the survival in ACC. As shown in the prognostic curve and scatter plot, the risk scores of ACC patients were negatively correlated with their OS.
Meanwhile, the cluster heatmap revealed the elevated expressions of candidate FerRLs in the high-risk group (Figure 3D), which needs to be further explored.
To explore whether the generated FerRLSig could independently predict the OS for patients with ACC, we then performed univariate and multivariate Cox regression analyses. The hazard ratios (HRs) of the risk scores in the univariate analysis and multivariate Cox regression analysis were 1.975 (1.527-2.533) and 1.936 (1.484-2.526), respectively (p < 0.001, Figures 3E,F). Additionally, the predictive accuracy of this novel signature was evaluated by a ROC analysis, with the AUC values of 0.936, 0.983, and 0.989 for the 1-, 2-, and 3-year survival rates, respectively (Figure 2G). To make the model more applicable in clinical management, a nomogram incorporating the FerRLSig and clinicopathological features was constructed to predict the 1- ,3-, and 5-year OS of ACC patients based on TCGA database (Figure 4A). The cluster heatmap for the association of clinicopathological manifestations and FerRLSig is illustrated in Figure 4B.
Identification of biological pathways in the prognostic FerRLSig
A total of 178 enriched KEGG pathways were identified in the GSEA (GSEA 4.2.2). Most of the identified FerRLs were involved in tumor-related and immunological pathways including cell cycle (NES = 2.07; p < 0.001), p53 signaling
A
B
C
D
E
Enrichment plot: KEGG_CELL_CYCLE
Enrichment plot: KEGG_P53_SIGNALING_PATHWAY
Enrichment plot: KEGG_SMALL_CELL_LUNG_CANCER
Enrichment plot: KEGG_PANCREATIC_CANCER
Enrichment score (ES)
Enrichment score (28]
Enrichment score (BS)
Enrichment plot: KEGG_HOMOLOGOUS_RECOMBINATION
0.7
Enrichment score (ES)
0.6
Enrichment score (ES)
Q.B
0.7
F
A
.3
G
7
N
A
V
A
0
1.0
0
Ranked list mette (Bigmac2Naise)
Ranked list metric (Bigmal2Naise)
Ranked list metric (Signal2Naise)
Ranked list metric (Bignal”Naise)
Ranked list methe (5 gnal_Tara*)
a
AS
·
Y [segitivity caselate-d)
0
10,00€
20.000
23 000
4 300
50 200
0
10.000
20 000
23.000
C 300
50.000
2
29:#30
30 300
50:000
0
10.000
20,000
60.000
a
10.000
20.000
30.000
Rark in Ordered Dataset
Rank in Ordered Dareset
Rank in Codered Dwasel
Ranik in Orderes Dataset
Rank in Orderad Delaset
Enrichment profile - Hits
Ranking rietris soomes
Enrichment profile -His
Ranking riti: stores
Enrichment profile - Hits
Ranking retnie soomes
Enrichment pictle -Hr’s
Ranking mettic scores
Enricinment profile - Hits
Ranking metric scores
F
G
H
I
J
Enrichment plot: KEGG_PYRIMIDINE_METABOLISM
Enrichment plot:
Enrichment score (ES)
Enrichment plot: KEGG_TOF_BETA_SIGNALING_PATHWAY
Enrichment plot: KEGG_PRIMARY_BILE_ACID_BIOSYNTHESIS
KEGG_PATHOGENIC_ESCHERICHIA_COLI_INFECTION
Enrichment plot: KEGG_DRUG_METABOLISM_CYTOCHROME_P450
-0.6
2
Drrichment score (CS)
Enrichment score (ES)
115
Enrichment score (ES)
Enrichment score (ES)
9
0.4
0.1
0.1
0
2
1
8
5.3
17
1.4
3.6
Ranseo int meine (99ignacio se)
Rankad list metric [Bigna/2Noise)
Farkad lidt mat tic [Signal)No ls)
Ranked list metric [Signal2Naine]
Ranked list metic [8 gnar2Noise]
E
.
.
6
27.000
13.000
10.000
20.000
4 300
0
20.000
30.000
40.000
00000
.
16.000
20.000
30.070
40,000
ø
1.000
20.000
40.000
00.030
Rank In Oreered Datset
Bank If Offered Dansstt
Rank in OMeres CMaset
Hank In Offerte Datesel
-Envichmant profita - Hits
Ranking ristris soonas
Envichman prodia - His
Ranking risti scores
-Ewichmint profile - Hits
Ranking matric soonas
-Enichmert prefle -Hi’s
Ranking mettic scores
Erwichmont prodia - Hits
Ranking matic scores
FIGURE 5 GSEA of the ferroptosis-related signaling pathways according to TCGA database. (A) Cell cycle, (B) p53 signaling pathway, (C) Small cell lung cancer, (D) Pancreatic cancer, (E) Homologous recombination, (F) Pyrimidine metabolism, (G) Pathogenic Escherichia coli infection, (H) TGF beta signaling pathway, (I) Primary bile acid biosynthesis, (J) Drug metabolism cytochrome P450 signaling pathway.
pathway (NES = 1.92; p < 0.001), small-cell lung cancer (NES = 1.75; p = 0.002), pancreatic cancer (NES = 1.74; p < 0.001), homologous recombination (NES = 1.91; p < 0.001), pyrimidine metabolism (NES = 1.89; p < 0.001), pathogenic Escherichia coli infection (NES = 1.78; p = 0.002), TGF-beta signaling pathway (NES = 1.67; p = 0.002), primary bile acid biosynthesis (NES = - 1.66; p = 0.013), and drug metabolism cytochrome P450 (NES = - 1.49; p = 0.027) signaling pathways (Figure 5 and Supplementary Table S7). Taken together, the novel prognostic signature may be involved in the microenvironment of tumor immune.
Ferroptosis-related genes and immunity
Figure 6 illustrates the immune responses clustering heatmap based on TIMER and CIBERSORT algorithms, which showed that immune infiltration was increased in cases with high-risk scores. Based on the ssGSEA of TCGA-ACC data, correlation analyses were performed to evaluate the association of immune cell subsets and corresponding functions, which revealed that T-cell function of type II INF response decreased significantly in the high-risk group (Figure 7A), suggesting that high-risk patients might be related to a suppressive immune microenvironment. Given the clinical importance of the checkpoint blocking therapy in tumors, we further analyzed the expression differences of immune checkpoints between the
two groups, and the expressions of CD40, CD276, IDO2, NRP1, and CD80 were of significant differences (Figure 7B). The comparison of m6A-related mRNA expression between the two risk groups showed that RBM15, WTAP, YTHDF1, and HNRNPC were differentially expressed (Figure 7C). The aforementioned results indicate that the generated FerRLSig might be a candidate marker for the immune checkpoint inhibitor therapy in ACC patients.
Discussion
ACC is an uncommon endocrine malignancy associated with high mortality (Wang et al., 2019b). Early diagnosis and more accurate therapies are of great importance to improve the prognosis. Some nomograms and risk models have been constructed to predict the survival in ACC (Wang et al., 2021a; Chen et al., 2021); however, no prognostic FerRLSig has been identified. Ferroptosis is associated with the development and progression of malignancies and might be a potentially novel approach for tumor therapy; additionally, its pivotal role in ACC has recently been revealed (Belavgeni et al., 2019). In this study, a FerRLSig that can evaluate the prognosis of ACC was constructed based on TCGA-ACC dataset. We further explored the association of survival and immune infiltrating cells and immune checkpoint inhibitors, which could provide potential biomarkers and immune targets for ACC treatment.
riskScore risk
riskScore
20
NK cell resting_CIBERSORT
5
-80
Macrophage MO_CIBERSORT
risk
low
NK cell resting_CIBERSORT-ABS
high
0
Methods
Macrophage MO_CIBERSORT-ABS
CIBERSORT
CIBERSORT-ABS
QUANTISEQ
B cell_QUANTISEQ
MCPCOUNTER
-5
XCELL
EPIC
Monocyte_QUANTISEQ
T cell CD4+ (non-regulatory)_QUANTISEQ
T cell regulatory (Tregs)_QUANTISEQ
Myeloid dendritic cell_QUANTISEQ
Cancer associated fibroblast_MCPCOUNTER
T cell CD4+ effector memory_XCELL
Monocyte_XCELL
T cell CD4+ Th1_XCELL
T cell CD4+ Th2_XCELL
B cell_EPIC
Cancer associated fibroblast_EPIC
T cell CD4+_EPIC
Methods
FIGURE 6 Immune response correlated to the ferroptosis-related lncRNA signature, illustrated in a cluster heatmap.
A
B
Risk
low
high
C
Risk
low
high
Risk
#
low
E
high
6
ns
ns
ns
ns
ns
ns
ns
ns
1.00
ns
ns
ns
ns
ns
ns
ns
ns
ns
ns
ns
ns
6
Gene expression
0.75
4
Gene expression
Score
4.
.50
2
0.25
S
2
0.00
0
:
APC_co_inhibition
APC_co_stimulation
CCR
Check-point
Cytolytic_activity
HLA
Inflammation-promoting
MHC_class_J
Parainflammation
T_cell_co-inhibition
T_cell_co-stimulation
Type_I_IFN_Reponse
Type_Il_IFN_Reponse
CD40
CD276
IDO2
NRP1
CD80
RBM15
WTAP
FTO
YTHDC1
ZC3H13
METTL3
YTHDC2
METTL14
YTHDF2
ALKBH5
YTHDF1
HNRNPC
FIGURE 7 (A) Correlation analysis for immune cell subsets and the related functions by ssGSEA. (B) Comparisons of immune checkpoints expressions between the low- and high-risk groups. (C) Expressions of m6A-related genes based on the FerRLSig.
Overall, 24 FerRLs were identified in the prognostic FerRLSig. In review of the literature, cell types that predominantly express these FerRLs along with the potential mechanisms are presented in
Supplementary Table S4. The AUC value of the ROC analysis indicated that the FerRLSig might have a superior predictive performance for OS in ACC. Additionally, a nomogram
incorporating the FerRLSig and clinicopathological features of ACC was constructed to facilitate the development of treatment strategies and might guide clinical decision-making. One recent study revealed that overexpression of lncRNA AC099850.3 could promote cell proliferation and invasion via the PRR11/PI3K/AKT pathway, predicting a poor prognosis in hepatocellular carcinoma (Zhong et al., 2022). In a related study, LINC01614 was significantly upregulated in osteosarcoma and could facilitate tumor progression through the miR-520a-3p/SNX3 axis (Cai et al., 2021). In addition, the overexpression of LINC01614 in bladder cancer could promote tumor proliferation, migration, and invasion through the miR-217/ RUNX2 and Wnt/B-Catenin pathway (Wang et al., 2021b). Upregulated expression of lncRNA LASTR was found in several epithelial tumors and could facilitate cancer cell fitness in hypoxic breast cancer through SART3 (De Troyer et al., 2020). LINC01503 is involved in a variety of malignancies, including hepatocellular carcinoma, cholangio-carcinoma, non-small-cell lung cancer, cervical cancer, and gastric cancer (Qu et al., 2019; Shen et al., 2020; Wang et al., 2020; Zhang et al., 2020; Feng et al., 2021; Ma et al., 2021). Wu et al. (2021) reported that the overexpression of LINC00460 in pancreatic ductal adenocarcinoma could promote tumor aggressiveness by targeting miR-491-5p, conferring a poor clinical outcome. Similarly, LINC00460 in head and neck squamous cell carcinoma facilitates peroxiredoxin-1 entry into the nucleus, thereby promoting cell proliferation and metastasis (Jiang et al., 2019). Knockdown of lncRNA GAPLINC repressed the tumor growth in renal cell cancer, while overexpression of GAPLINC facilitated tumorigenesis via the miR-135b-5p/CSF1 axis, indicating poor survival (Wang et al., 2021c). On the other hand, Yang and colleagues revealed that lncRNA MYOSLID played a critical role in the progression of osteosarcoma through the miR-1286/RAB13 axis (Yang et al., 2019a). In cervical cancer, upregulated expression of LncRNA DGUOK-AS1 could sponge miR-653-5p to repress DNA repair and promote cell proliferation (Wu et al., 2020). However, to the best of our knowledge, there has been no research conducted to explore whether FerRLs participate in the pathogenesis and prognosis of ACC.
Next, the present study identified 540 ferroptosis-related DEGs and stratified them by the novel FerRLSig. KEGG-based analyses further revealed that the genes of the high-risk group mainly participated in the Hippo signaling pathway, p53 signaling pathway, cancer-related pathway, ECM-receptor interaction, cellular senescence, and cell cycle. Recently, Lin and colleagues revealed that recurrent breast tumors are highly sensitive to ferroptosis, and upregulated expression of DDR2 could promote ferroptosis through the Hippo signaling pathway (Lin et al., 2021). Wen-Hsuan et al. reported that TAZ (a Hippo pathway effector) mediates the cell density-regulated ferroptosis by regulating the EMP1-NOX4 axis in renal cell carcinoma (Yang et al., 2019b). Another related study identified an ACSL4-independent and ALOX12-mediated pathway of ferroptosis that is essential for p53-dependent tumor suppression (Chu et al., 2019). Taken together, the
aforementioned signaling pathways are closely related to ferroptosis in tumors.
Another significant contribution of the current study is the uncovering of the association between the FerRLSig and tumor immune microenvironment. Increasing evidence suggests that T-cells could promote tumor ferroptosis to enhance antitumor activity, which might be a potential therapeutic target for cancers combining with an immune checkpoint blockade (Wang et al., 2019c). Functional enrichment analyses revealed that FerRLs mainly participated in tumor-related and immune pathways. Type II INF response in T-cell function was significantly decreased in the high-risk group, suggesting that FerRLs may be involved in the regulation of tumor immune infiltration.
As a novel form of cell death, ferroptosis provides a new way for tumor development and treatment (Hassannia et al., 2019; Liang et al., 2019). However, the exact mechanism of ferroptosis in ACC remains unclear. To date, we are the first to construct a FerRLsig for ACC, which could shed light on novel treatment modalities for this rare tumor. Nonetheless, several limitations exist in the present study. First and most importantly, our findings are based on a single database and have not been validated by a clinical cohort. We explored the Gene Expression Omnibus (GEO) and the International Cancer Genome Consortium (ICGC); however, ACC is extremely rare and only 26 ACCs were found in the GEO, and no case of ACC was found in the ICGC. Thus, we could not compare or validate our findings with other data sets. In addition, the retrospective nature of this study made it less applicable. In addition, potential predictors such as tumor markers and tumor grade were not included in the research due to the incomplete data. Furthermore, functional experimental validation of the molecular mechanisms of the generated FerRLSig is currently lacking. In general, the role of the prognostic signature established in this study needs further exploration, and international multicenter studies with a larger sample size are expected.
In conclusion, a FerRLSig was established to predict the survival of ACC. Additionally, this model might be related to immune infiltration, providing potential immune targets for the control of ACC.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.
Author contributions
WW and CY conceived and designed the study. WW and GC wrote and revised the manuscript. RZ collected and interpreted
the data. All authors have contributed to the article and approved the submitted version.
Funding
The present study was supported by the National Natural Science Foundation of China (Grant Number: 82000084 and 82002828).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Belavgeni, A., Bornstein, S. R., von Massenhausen, A., Tonnus, W., Stumpf, J., Meyer, C., et al. (2019). Exquisite sensitivity of adrenocortical carcinomas to induction of ferroptosis. Proc. Natl. Acad. Sci. U. S. A. 116 (44), 22269-22274. doi:10.1073/pnas.1912700116
Buishand, F. O., Liu-Chittenden, Y., Fan, Y., Tirosh, A., Gara, S. K., Patel, D., et al. (2020). Adrenocortical tumors have a distinct, long, non-coding RNA expression profile and LINC00271 is downregulated in malignancy. Surgery 167 (1), 224-232. doi:10.1016/j.surg.2019.04.067
Cai, Q., Zhao, X., Wang, Y., Li, S., Wang, J., Xin, Z., et al. (2021). LINC01614 promotes osteosarcoma progression via miR-520a-3p/SNX3 axis. Cell. Signal. 83, 109985. doi:10.1016/j.cellsig.2021.109985
Charoentong, P., Finotello, F., Angelova, M., Mayer, C., Efremova, M., Rieder, D., et al. (2017). Pan-cancer immunogenomic analyses reveal genotype- immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 18 (1), 248-262. doi:10.1016/j.celrep.2016.12.019
Chasov, V., Mirgayazova, R., Zmievskaya, E., Khadiullina, R., Valiullina, A., Stephenson Clarke, J., et al. (2020). Key players in the mutant p53 team: Small molecules, gene editing, immunotherapy. Front. Oncol. 10, 1460. doi:10.3389/fonc.2020.01460
Chen, X., Yan, L., Jiang, F., Lu, Y., Zeng, N., Yang, S., et al. (2021). Identification of a ferroptosis-related signature associated with prognosis and immune infiltration in adrenocortical carcinoma. Int. J. Endocrinol. 2021, 4654302. doi:10.1155/2021/4654302
Chu, B., Kon, N., Chen, D., Li, T., Liu, T., Jiang, L., et al. (2019). ALOX12 is required for p53-mediated tumour suppression through a distinct ferroptosis pathway. Nat. Cell Biol. 21 (5), 579-591. doi:10.1038/s41556-019-0305-6
De Troyer, L., Zhao, P., Pastor, T., Baietti, M. F., Barra, J., Vendramin, R., et al. (2020). Stress-induced IncRNA LASTR fosters cancer cell fitness by regulating the activity of the U4/U6 recycling factor SART3. Nucleic Acids Res. 48 (5), 5198-5199. doi:10.1093/nar/gkaa280
Dixon, S. J., Lemberg, K. M., Lamprecht, M. R., Skouta, R., Zaitsev, E. M., Gleason, C. E., et al. (2012). Ferroptosis: An iron-dependent form of nonapoptotic cell death. Cell 149 (5), 1060-1072. doi:10.1016/j.cell.2012.03.042
Else, T., Kim, A. C., Sabolch, A., Raymond, V. M., Kandathil, A., Caoili, E. M., et al. (2014). Adrenocortical carcinoma. Endocr. Rev. 35 (2), 282-326. doi:10.1210/ er.2013-1029
Fassnacht, M., Libe, R., Kroiss, M., and Allolio, B. (2011). Adrenocortical carcinoma: A clinician’s update. Nat. Rev. Endocrinol. 7 (6), 323-335. doi:10. 1038/nrendo.2010.235
Feng, J., Gao, F. Y., Li, Y. Y., Xiao, X. Y., and Xi, J. (2021). Upregulation of LINC01503 promotes cervical cancer progression by targeting the miR-615-3p/ CCND1 axis. J. Cancer 12 (15), 4552-4560. doi:10.7150/jca.54148
Gupta, R. A., Shah, N., Wang, K. C., Kim, J., Horlings, H. M., Wong, D. J., et al. (2010). Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature 464 (7291), 1071-1076. doi:10.1038/nature08975
Hassannia, B., Vandenabeele, P., and Vanden Berghe, T. (2019). Targeting ferroptosis to iron out cancer. Cancer Cell 35 (6), 830-849. doi:10.1016/j.ccell.2019.04.002
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene. 2022.949457/full#supplementary-material
Huang, L., McClatchy, D. B., Maher, P., Liang, Z., Diedrich, J. K., Soriano-Castell, D., et al. (2020). Intracellular amyloid toxicity induces oxytosis/ferroptosis regulated cell death. Cell Death Dis. 11 (10), 828. doi:10.1038/s41419-020-03020-9
Jiang, Y., Cao, W., Wu, K., Qin, X., Wang, X., Li, Y., et al. (2019). LncRNA LINC00460 promotes EMT in head and neck squamous cell carcinoma by facilitating peroxiredoxin-1 into the nucleus. J. Exp. Clin. Cancer Res. 38 (1), 365. doi:10.1186/s13046-019-1364-z
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). Timer: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77 (21), e108-e110. doi:10.1158/0008-5472.CAN-17-0307
Liang, C., Zhang, X., Yang, M., and Dong, X. (2019). Recent progress in ferroptosis inducers for cancer therapy. Adv. Mat. 31 (51), 1904197. doi:10. 1002/adma.201904197
Lin, C. C., Yang, W. H., Lin, Y. T., Tang, X., Chen, P. H., Ding, C. C., et al. (2021). DDR2 upregulation confers ferroptosis susceptibility of recurrent breast tumors through the Hippo pathway. Oncogene 40 (11), 2018-2034. doi:10.1038/s41388- 021-01676-x
Liu, T., Li, X., Cui, Y., Meng, P., Zeng, G., Wang, Y., et al. (2021). Bioinformatics analysis identifies potential ferroptosis key genes in the pathogenesis of intracerebral hemorrhage. Front. Neurosci. 15, 661663. doi:10.3389/fnins.2021. 661663
Ma, Z., Gao, X., Shuai, Y., Wu, X., Yan, Y., Xing, X., et al. (2021). EGR1-mediated linc01503 promotes cell cycle progression and tumorigenesis in gastric cancer. Cell Prolif. 54 (1), 12922. doi:10.1111/cpr.12922
Manz, D. H., Blanchette, N. L., Paul, B. T., Torti, F. M., and Torti, S. V. (2016). Iron and cancer: Recent insights. Ann. N. Y. Acad. Sci. 1368 (1), 149-161. doi:10. 1111/nyas.13008
Mao, C., Wang, X., Liu, Y., Wang, M., Yan, B., Jiang, Y., et al. (2018). A G3BP1- interacting lncRNA promotes ferroptosis and apoptosis in cancer via nuclear sequestration of p53. Cancer Res. 78 (13), 3484-3496. doi:10.1158/0008-5472. CAN-17-3454
Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 12 (5), 453-457. doi:10.1038/nmeth.3337
Qu, Y. K., Qu, X. S., Chen, G., Feng, Y., Teng, X. L., Liu, W. X., et al. (2019). LINC01503 promotes cell proliferation, invasion and EMT process in cholangio- carcinoma. Eur. Rev. Med. Pharmacol. Sci. 23 (15), 6445-6452. doi:10.26355/ eurrev_201908_18526
Shen, Q., Sun, Y., and Xu, S. (2020). LINC01503/miR-342-3p facilittes malignancy in non-small-cell lung cancer cells via regulating LASP1. Respir. Res. 21 (1), 235. doi:10.1186/s12931-020-01464-3
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102 (43), 15545-15550. doi:10.1073/pnas.0506580102
Vickers, A. J., and Elkin, E. B. (2006). Decision curve analysis: A novel method for evaluating prediction models. Med. Decis. Mak. 26 (6), 565-574. doi:10.1177/ 0272989X06295361
Wang, M., Mao, C., Ouyang, L., Liu, Y., Lai, W., Liu, N., et al. (2019). Long noncoding RNA LINC00336 inhibits ferroptosis in lung cancer by functioning as a competing endogenous RNA. Cell Death Differ. 26 (11), 2329-2343. doi:10.1038/s41418-019-0304-y
Wang, M. R., Fang, D., Di, M. P., Guan, J. L., Wang, G., Liu, L., et al. (2020). Long non- coding RNA LINC01503 promotes the progression of hepatocellular carcinoma via activating MAPK/ERK pathway. Int. J. Med. Sci. 17 (9), 1224-1234. doi:10.7150/ijms.45256
Wang, S., Yang, X., Xie, W., Fu, S., Chen, Q., Li, Z., et al. (2021). LncRNA GAPLINC promotes renal cell cancer tumorigenesis by targeting the miR-135b-5p/ CSF1 Axis. Front. Oncol. 11, 718532. doi:10.3389/fonc.2021.718532
Wang, W., Chang, G., Sun, Y., Zhuo, R., Li, H., Hu, Y., et al. (2021). Nomograms for individualized evaluation of prognosis in adrenocortical carcinomas for the elderly: A population-based analysis. J. Invest. Surg. 35, 1153-1160. doi:10.1080/ 08941939.2021.1968981
Wang, W., Green, M., Choi, J. E., Gijon, M., Kennedy, P. D., Johnson, J. K., et al. (2019). CD8(+) T cells regulate tumour ferroptosis during cancer immunotherapy. Nature 569 (7755), 270-274. doi:10.1038/s41586-019-1170-y
Wang, W., Han, R., Ye, L., Xie, J., Tao, B., Sun, F., et al. (2019). Adrenocortical carcinoma in patients with MEN1: A kindred report and review of the literature. Endocr. Connect. 8 (3), 230-238. doi:10.1530/EC-18-0526
Wang, Z., Yan, H., Cheng, D., Xu, L., Shen, T., Chen, Y., et al. (2021). Novel lncRNA LINC01614 facilitates bladder cancer proliferation, migration and invasion through the miR-217/RUNX2/Wnt/B-Catenin Axis. Cancer Manag. Res. 13, 8387-8397. doi:10.2147/CMAR.S330019
Weigand, I., Schreiner, J., Rohrig, F., Sun, N., Landwehr, L. S., Urlaub, H., et al. (2020). Active steroid hormone synthesis renders adrenocortical cells highly susceptible to type II ferroptosis induction. Cell Death Dis. 11 (3), 192. doi:10.1038/s41419-020-2385-4
Wu, J., Sun, S., Liao, W., Chen, E., Wang, X., Song, Y., et al. (2021). LINC00460 promotes pancreatic cancer progression by sponging miR-491-5p. J. Gene Med. 23 (6), 3333. doi:10.1002/jgm.3333
Wu, N., Song, H., Ren, Y., Tao, S., and Li, S. (2020). DGUOK-AS1 promotes cell proliferation in cervical cancer via acting as a ceRNA of miR-653-5p. Cell biochem. Funct. 38 (7), 870-879. doi:10.1002/cbf.3506
Yang, S., Chen, M., and Lin, C. (2019). A novel lncRNA MYOSLID/miR-1286/ RAB13 Axis plays a critical role in osteosarcoma progression. Cancer Manag. Res. 11, 10345-10351. doi:10.2147/CMAR.S231376
Yang, W. H., Ding, C. C., Sun, T., Rupprecht, G., Lin, C. C., Hsu, D., et al. (2019). The hippo pathway effector TAZ regulates ferroptosis in renal cell carcinoma. Cell Rep. 28 (10), 2501-2508. doi:10.1016/j.celrep.2019.07.107
Zhang, M. L., Zhao, T. T., Du, W. W., Yang, Z. F., Peng, W., and Cui, Z.J. (2020). C-MYC-induced upregulation of LINC01503 promotes progression of non-small cell lung cancer. Eur. Rev. Med. Pharmacol. Sci. 24 (21), 11120-11127. doi:10.26355/ eurrev_202011_23599
Zhang, T., Li, K., Zhang, Z. L., Gao, K., and Lv, C. L. (2021). LncRNA Airsci increases the inflammatory response after spinal cord injury in rats through the nuclear factor kappa B signaling pathway. Neural Regen. Res. 16 (4), 772-777. doi:10.4103/1673-5374.295335
Zhang, Y., Guo, S., Wang, S., Li, X., Hou, D., Li, H., et al. (2021). LncRNA OIP5- AS1 inhibits ferroptosis in prostate cancer with long-term cadmium exposure through miR-128-3p/SLC7A11 signaling. Ecotoxicol. Environ. Saf. 220, 112376. doi:10.1016/j.ecoenv.2021.112376
Zheng, S., Cherniack, A. D., Dewal, N., Moffitt, R. A., Danilova, L., Murray, B. A., et al. (2016). Comprehensive pan-genomic characterization of adrenocortical carcinoma. Cancer Cell 29 (5), 723-736. doi:10.1016/j.ccell. 2016.04.002
Zhong, F., Liu, S., Hu, D., and Chen, L. (2022). LncRNA AC099850.3 promotes hepatocellular carcinoma proliferation and invasion through PRR11/PI3K/AKT axis and is associated with patients prognosis. J. Cancer 13 (3), 1048-1060. doi:10. 7150/jca.66092
Zhou, N., and FerrDb, Bao J. (2020). A manually curated resource for regulators and markers of ferroptosis and ferroptosis-disease associations. Database J. Biol. databases curation, 2020. doi:10.1093/database/baaa021