Zhihao Yang
Yaoshu Song1,3†
Guobo Du
Hongpan Zhang- 1Department of Oncology, Affiliated Hospital of North Sichuan Medical College, Nanchong, China
- 2Tianjin Key Laboratory of Medical Epigenetics, Key Laboratory of Breast Cancer Prevention and Therapy (Ministry of Education), Department of Biochemistry and Molecular Biology, Tianjin Medical University, Tianjin, China
- 3North Sichuan Medical College, Nanchong, China
- 4Department of Pathology and Medical Research Center, Beijing Chaoyang Hospital, Capital Medical University, Beijing, China
- 5Suzhou Kowloon Hospital, Shanghai Jiao Tong University School of Medicine, Suzhou, China
Extensive evidence has revealed that ferroptosis plays a vital role in HCC development and progression. Fanconi anemia complementation group D2 (FANCD2) has been reported to serve as a ferroptosis-associated gene and has a close relationship with tumorigenesis and drug resistance. However, the impact of the FANCD2-related immune response and its mechanisms in HCC remains incompletely understood. In the current research, we evaluated the prognostic significance and immune-associated mechanism of FANCD2 based on multiple bioinformatics methods and databases. The results demonstrated that FANCD2 was commonly upregulated in 15/33 tumors, and only the high expression of FANCD2 in HCC was closely correlated with worse clinical outcomes by OS and DFS analyses. Moreover, ncRNAs, including two major types, miRNAs and lncRNAs, were closely involved in mediating FANCD2 upregulation in HCC and were established in a ceRNA network by performing various in silico analyses. The DUXAP8-miR-29c-FANCD2 and LINC00511-miR-29c-FANCD2 axes were identified as the most likely ncRNA-associated upstream regulatory axis of FANCD2 in HCC. Finally, FANCD2 expression was confirmed to be positively related to HCC immune cell infiltration, immune checkpoints, and IPS analysis, and GSEA results also revealed that this ferroptosis-associated gene was primarily involved in cancer-associated pathways in HCC. In conclusion, our investigations indicate that ncRNA-related modulatory overexpression of FANCD2 might act as a promising prognostic and immunotherapeutic target against HCC.
Introduction
Hepatocellular carcinoma (HCC) remains the sixth and third leading cause of cancer-related morbidity and death globally among men and women (Sung et al., 2021; Siegel et al., 2022). HCC accounts for 75–85% of primary liver cancer cases, with an increasing incidence in China and the United States (Islami et al., 2017; Tapper and Parikh, 2018; Zhou et al., 2019). Despite advances in traditional treatment methods, including surgical resection, orthotopic liver transplantation, and locoregional therapies, the survival of patients with HCC is still particularly poor (Forner et al., 2018; Llovet et al., 2021). Currently, four forms of immune checkpoint inhibitors have proven to be promising and effective immunotherapy strategies for HCC (Killock, 2017; Vogel and Saborowski, 2022). However, only a subset of patients can benefit from immunotherapy, with a generally efficient response rate of no more than 20% (Cheng et al., 2020; Sangro et al., 2020). Our main aim is to identify immune-modulatory biomarkers for diagnosis, prognosis, and treatment of HCC.
Fanconi anemia complementation group D2 (FANCD2) belongs to one of the DNA repair markers, which is a class of Fanconi anemia (FA) proteins that regulate DNA damage responses and maintain genomic integrity (Liu et al., 2010; Semlow and Walter, 2021). Recent studies have demonstrated that FANCD2 in HCC cells maintains resistance to DNA-damaging tumor suppression by increasing DNA repair activity, FANCD2 deletion restricts chemoresistance of HCC cells in vitro (Palagyi et al., 2010), and elevated FANCD2 is closely linked to adverse prognosis in HCC (Komatsu et al., 2017), glioma cells (Chen et al., 2021), breast cancer (Fagerholm et al., 2013), colorectal cancer (Ozawa et al., 2010), and non-small-cell lung cancer (Ferrer et al., 2005). FANCD2 has been reported to elicit an active role in the negative regulation of ferroptosis (Song et al., 2016), a recently new form of programmed cell death characterized by iron-dependent peroxide lipid accumulation, leading to many forms of tumors that are sensitive to ferroptosis induction treatment (Dixon et al., 2012). However, the immune-regulating effect on this ferroptosis-associated gene in HCC onset and progression remains unexplored.
In our current study, we deciphered the role of ferroptosis-associated gene FANCD2 expression and prognosis in pancancers anchored by several databases. Subsequently, ncRNAs, including two major types, miRNAs and lncRNAs, were found to be closely involved in mediating FANCD2 in HCC and were analyzed through ceRNA. Finally, a high level of FANCD2 in HCC was also positively associated with immune cell infiltration, immune checkpoints, and IPS analysis, suggesting that ncRNA modulatory upregulation of FANCD2 is an immune-associated biomarker with an adverse influence on HCC patient clinical survival.
Methods
Data collection and expression
The expression data of FANCD2 and corresponding clinical information in 33 types of tumors were processed from UCSC Xena (Navarro Gonzalez et al., 2021) (https://xena.ucsc.edu/) using R software (Version 4.0.3; https://www.R-project.org) and Strawberry Perl (5.3.0.1). The expression of FANCD2 in HCC was also verified by the Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB), which contains 15 public HCC expression datasets and provides candidate-targeted gene expression prognosis and coexpression analysis (http://lifeome.net/database/hccdb/search.html) (Lian et al., 2018). Furthermore, the Cancer Cell Line Encyclopedia (CCLE) database was used to determine the diverse expression of FANCD2 in HCC cell lines (https://portals.broadinstitute.org/ccle/about) (Barretina et al., 2019).
starBase online database and lncRNA subcellular locations
The starBase database (http://starbase.sysu.edu.cn/) (Li et al., 2014) was applied to verify FANCD2 expression and prognostic status in many tumors. Then, this website, including PITA, RNA22, miRmap, microT, miRanda, PicTar, and TargetScan, was adopted to predict the upstream miRNA of target genes (FANCD2) and the expression and survival status of screened miRNAs in HCC. Subsequently, to ascertain candidate lncRNAs via a competing endogenous RNA (ceRNA) mechanism, the upstream lncRNAs of the miR-29c-3p-FANCD2 axis were predicted, and the web was conducted to examine lncRNA expression and its correlation with FANCD2 together with miR-29c-3p. Finally, Cytoscape software (v3.6.0) was used to visualize the FANCD2-regulated ceRNA network.
Arched by binomial distribution approaches, the iLoc-lncRNA database (Su et al., 2018) and lncLocator database (Lin et al., 2021) were applied to accurately predict the subcellular locations of lncRNAs, encompassing four locations of lncRNAs in a cell (nucleus, cytoplasm, ribosome, and exosome).
Overall survival and disease-free survival analyses
Kaplan–Meier (KM) survival curves, including overall survival (OS) and disease-free survival (DFS) analyses, combined with a log-rank test mapped differences of FANCD2 for all types of TCGA tumors employing the GEPIA tool (http://gepia.cancer-pku.cn) (Tang et al., 2019), which was validated by the HCCDB (Lian et al., 2018). p < 0.05 was treated as the significance threshold for all statistical analyses.
Analysis of tumor immunity aimed at Fanconi anemia complementation group D2
The association between FANCD2 expression and infiltrating immune cells, such as B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells, was analyzed across the HCC patients captured from the TIMER2 web server (https://cistrome.shinyapps.io/timer/) (Li et al., 2020), which was incorporated into 10,897 samples from TCGA.
Immunotherapy has proven to be an effective and promising strategy against tumors, especially involving the use of CTLA-4, PD-1, and PD-L1 antibodies to treat several types of tumors (Galluzzi et al., 2020). Therefore, we used the TIMER2 database to explore the interplay between FANCD2 and these immune checkpoints. The immunophenoscore (IPS), which can be obtained from the TCIA website (https://tcia.at/home), has been verified to predict the patients’ response to immune checkpoint inhibitors (ICIs). A higher IPS has a better outcome with ICI treatment (Charoentong et al., 2017).
Gene set enrichment analysis
Gene set enrichment analysis (GSEA) (Subramanian et al., 2005), including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses, was used to determine the potential mechanisms and functions of FANCD2 in HCC, with the following parameters: nPerm = 1,000, minGSSize = 10, maxGSSize = 1,000, and nominal p-value < 0.05.
Results
Overview of the workflow
Our goal was to detect the potential function of FANCD2 in various TCGA cancer types and to investigate the correlation of immunotherapy between FANCD2 and HCC (Figure 1). To this end, we collected various public data sources and conducted a series of bioinformatics and functional genomic analyses.
FIGURE 1. Research workflow of construction and analysis of the FANCD2 ceRNA network.
Fanconi anemia complementation group D2 mRNA expression across cancers
Aberrantly expressed analysis from TCGA data revealed that FANCD2 was differentially expressed in many human tumor tissues, and that FANCD2 mRNA overexpression was found in BLCA, BRCA, CHOL, COAD, ESCA, GBM, HNSC, KIRC, KIRP, HCC, LUSC, PRAD, THCA, and UCEC, whereas FANCD2 loss was found in KICH (Figure 2A). We further evaluated whether FANCD2 expression is high in most human cancers compared with matched normal samples utilizing the starBase database. This result indicated that FANCD2 mRNA expression in BLCA, BRCA, CHOL, COAD, GBM, HCC, LUSC, READ, STAD, and UCEC was also high (Figures 2B–K), suggesting that FANCD2 acts as an oncogene in multiple tumors, which was also verified by the HCCDB (Supplementary Figure S1A and Table 1).
FIGURE 2. Pancancer expression profiling analysis of FANCD2. (A) Interactive body map of FANCD2 mRNA expression constructed by GEPIA. (B–K) Expression of FANCD2 was analyzed between TCGA cancer tissues and GTEx normal tissues (*p < 0.05, **p < 0.01, ***p < 0.001).
TABLE 1. Differential expression analysis of FANCD2 from clinical cohorts impinged on the HCCDB.
Prognostic implication of Fanconi anemia complementation group D2 in screened cancers
Next, we used the GEPIA database within the RNA sequencing data from TCGA and GTEx data to ascertain the prognostic implication of FANCD2, including OS and DFS, in the aforementioned 10 cancers (BLCA, BRCA, CHOL, COAD, GBM, HCC, LUSC, READ, STAD, and UCEC). The results indicated that only the high level of FANCD2 showed worse OS in HCC (p = 0.0047, HR = 1.7, Figure 3), which was consistent with DFS in HCC (p = 0.00017, HR = 1.8) (Figure 4). Intriguingly, we also used the CCLE database and found that FANCD2 in most HCC cell lines had at least partly high levels (Supplementary Figure S1B). The finding that high FANCD2 activity in HCC was associated with unfavorable OS (p = 0.00767 and 0.00325, Supplementary Figures S1C,D) indicates that overexpressed FANCD2 might be a crucial prognostic factor in HCC.
FIGURE 3. Correlation of FANCD2 expression with patient OS. (A–J) K-M method showed that high expression of FANCD2A was closely correlated with shorter lifetimes in BLCA, BRCA, CHOL, COAD, GBM, LIHC, LUSC, READ, STAD, and UCEC.
FIGURE 4. Correlation of FANCD2 expression with patient DFS. (A–J) K-M method showed that increased expression of FANCD2A was closely related to disease-free survival lifetimes in BLCA, BRCA, CHOL, COAD, GBM, LIHC, LUSC, READ, STAD, and UCEC.
Prediction and analysis of upstream miRNAs of Fanconi anemia complementation group D2
Non-coding RNAs, including miRNAs and lncRNAs, that account for nearly 98% of transcriptomes have long been thought of as pivotal players in the biological regulation of cancer cells, such as ceRNA and gene expression (Beermann et al., 2016). Therefore, our first aim was to identify possible miRNAs involved in modulating FANCD2. Notably, miRNAs directly bind to the 3′-untranslated region (UTR) of the corresponding target mRNA, leading to their repression or posttranscriptional gene silencing (Goodall and Wickramasinghe, 2021). These findings revealed that eight screened miRNAs directly interact with FANCD2 upstream (Figure 5A) and are negatively associated with hsa-let-7c-5p, hsa-miR-29a-3p, hsa-miR-34a-5p, and hsa-miR-29c-3p in HCC via correlation analysis (p < 0.05, Figures 5B–I). Among these miRNAs, hsa-miR-29c-3p was markedly distinct from FANCD2 expression (Figure 6A). Intriguingly, low hsa-miR-29c-3p in HCC revealed unfavorable patient OS through the K-M plotter database, and it was significantly low in the aforementioned tumors (Figures 6B–D). The 3’ UTR arising from FANCD2 that targeted miR-29c-3p was also successfully predicted by the starBase database (Figure 5E). These outcomes suggested that hsa-miR-29c-3p is the most likely regulatory miRNA of FANCD2 in HCC.
FIGURE 5. Identification of the miRNA-FANCD2 regulatory network in HCC. (A) FANCD2-associated miRNA network was visualized using Cytoscape. (B–I) Expression correlation between predicted miRNAs and FANCD2 in HCC analyzed by the starBase database.
FIGURE 6. Identification of miR-29c-3p as the most potential regulatory miRNA of FANCD2 in HCC. (A) Association between FANCD2 and miR-29c-3p expression is shown in the dot plot. (B) Expression distribution of miR-29c-3p in HCC. (C,D) Kaplan‒Meier survival analysis of miR-29c-3p. (F) Base pairing of miR-29c-3p of the 3′ UTR in FANCD2-binding sites.
Prediction and analysis of upstream lncRNAs of hsa-miR-29c-3p
In addition to screened miRNAs, we predicted hsa-miR-29c-3p-modulatory lncRNAs in HCC using the starBase database and observed eight possible lncRNAs (Figures 7A–J). Anchored by the ceRNA hypothesis, lncRNAs have been viewed as sponges of miRNAs to restrict the reduction in the role of target mRNAs (Salmena et al., 2011). Accordingly, DUXAP8 and LINC00511 in HCC were elevated compared with those in normal samples, and their high levels were reciprocally relevant to better patient OS among all eight lncRNAs (Figures 8A–D). More importantly, DUXAP8 and LINC00511 were strongly associated with FANCD2 expression but were antithetical to hsa-miR-29c-3p, according to the starBase database (Figures 8E,F).
FIGURE 7. Expression of upstream lncRNAs predicted by miR-29c-3p. (A) AC098828.2, (B) H19, (C) LINC0068, (D) LINC0051, (E) SNHG17, (F) LINC01270, (G) DUXAP8, (H) TUGI, (I) LINC01521, (J) AL031587.5.
FIGURE 8. DUXAP8 and LINC00511 are the most likely upstream lncRNAs of the hsa-miR-29c-3p-FANCD2 axis in HCC. (A–D) Kaplan‒Meier survival analysis of DUXAP8 and LINC00511 through OS and DFS analyses. (E,F) Comparison of DUXAP8 and LINC00511 and FANCD2 expression in HCC. (G) Base pairing of DUXAP8 and LINC00511 of the 3′ UTR in miR-29c-5-binding sites. (H) Prediction of DUXAP8 and LINC00511 subcellular localization using lncLocator. (I) FANCD2-associated ceRNA network was visualized using Cytoscape.
To determine whether these two lncRNAs were involved in regulating FANCD2-based ceRNAs, we used the lncLocator database and the iLoc-lncRNA database to analyze their subcellular localization. As seen in Figures 8H,I and Supplementary Figures 2A,B, DUXAP8 was mainly distributed in the cytoplasm, while LINC00511 was mostly located in the ribosome, indicating that the former in HCC are more likely to participate in the FANCD2-specific ceRNA network. The target sequence of the 3’ UTR of DUXAP8 and LINC00511 paired with miR-29c-3p was predicted via the starBase online database (Figure 8G). To this end, both lncRNAs most likely are upstream lncRNAs of the hsa-miR-29c-3p-FANCD2 axis in HCC, especially DUXAP8 (Figure 8J).
Fanconi anemia complementation group D2 induces immune infiltration in the hepatocellular carcinoma tumor microenvironment
It has been well documented that infiltrating immune cells exert a crucial role in the tumor microenvironment (TME), which can regulate the initiation and development of cancers (Su et al., 2018; Cózar et al., 2021; Jhunjhunwala et al., 2021). When we examined the tumor immunity of FANCD2, we found that the FANCD2 expression level in HCC showed a significant interaction with immune infiltration (Figure 9A). As shown in Figure 9A, the FANCD2 expression level in HCC had an enormously positive correlation with neutrophils, purity, macrophages, CD4+ T cells, CD8+ T cells, dendritic cells, and B cells, indicating that FANCD2 could be of great importance in the regulation of the immune infiltration function of HCC.
FIGURE 9. Correlation of FANCD2 expression with immune infiltration level (A) and immune checkpoints (B–D) in HCC.
Interaction between Fanconi anemia complementation group D2 and immune checkpoints and IPS analysis
Immune checkpoints, including PD-1/PD-L1 and CTLA-4, have gained attention in cancer-associated therapy via immune surveillance and escape (Dall'Olio et al., 2022; Doroshow et al., 2021). We ascertained the role of the immune response of FANCD2 in HCC with PD-1, PD-L1, or CTLA-4 using TIMER data and found that FANCD2 expression was strongly statistically associated with PD-1, PD-L1, and CTLA-4 in HCC (Figures 9B–D), demonstrating that FANCD2 might mediate HCC carcinogenesis of immune escape.
The IPS values (ips_ctla4_neg_pd1_pos, ips_ctla4_pos_pd1_neg, ips_ctla4_pos_pd1_neg, and ips_ctla4_pos_pd1_pos) have been validated in predicting cancer patients’ responses to anti-CTLA-4 treatment. Here, we noted a marked rise in scores in the low-FANCD2 group (Figures 10A–D; p < 0.05), implying that the low-FANCD2 group is more suitable for anti-CTLA-4 and anti-PD-1/PD-L1 immunotherapy.
FIGURE 10. Correlation of FANCD2 expression with the immunophenoscore (IPS) in HCC. (A) ips_ctla4_neg_pd1_neg, (B) ips_ctla4_neg_pd1_pos, (C) ips_ctla4_pos_pd1_neg, (D) ips_ctla4_pos_pd1_pos.
Enrichment analysis of Fanconi anemia complementation group D2
To predict the functions or pathways that FANCD2 in HCC may modulate, GSEA was carried out. The results showed that activated FANCD2 was positively enriched in the cell cycle, p53 signaling pathway, mTOR signaling pathway, DNA replication, TGF-beta signaling pathway, Wnt signaling pathway, and RNA degradation (p < 0.05; Figure 11), revealing the tumor-associated and even tumor-promoting roles of FANCD2 and shedding light on its potential mechanisms in HCC.
FIGURE 11. GSEA functional enrichment of FANCD2 in HCC.
Discussion
HCC is a malignant tumor that threatens human health (Llovet et al., 2021). Despite advancements in immunotherapy against HCC, only a minority of patients are eventually able to prolong their survival time (Cheng et al., 2020; Faivre et al., 2020; Pinter et al., 2021). The identification of better immune-associated prognostic biomarkers for timely detection, early diagnosis, and treatment of HCC is urgently needed. Accumulating evidence has shown that FANCD2 can significantly resist ferroptosis-related concentrations of iron and lipid ROS in various diseases (Song et al., 2016; Li et al., 2021), acting as a ferroptosis-related regulator. Ferroptosis can also regulate immune activity within the HCC microenvironment (Zhao et al., 2020). Nevertheless, little progress has been achieved in studying the roles and mechanisms of this ferroptosis-related gene against HCC.
Our current study investigated whether the ferroptosis-related gene FANCD2 in the pancancer analysis was more activated than in counterpart tissues, as validated in the starBase database, other databases, and many published articles. Survival analysis further revealed that only FANCD2 at a high expression level had shorter OS and DFS survival in HCC, indicating that this gene could serve as a potential and prognostic target in this deadly cancer.
It has been well documented that lncRNAs can act as ceRNA networks by relieving the suppressive effects of targeted genes through miRNAs to facilitate their expression (Salmena et al., 2011; Tay et al., 2014). Our study found that there were eight possible miRNAs that contained binding sites for FANCD2 identified by bioinformatics analysis. Further analysis of only hsa-miR-29c was incorporated in our following study, as this miRNA deficiency was closely associated with a more favorable prognosis and reversely correlated with FANCD2 mRNA expression, implying that hsa-miR-29c is more likely to modulate FANCD2 mRNA stability. This hypothesis is supported by previous studies, reporting that miR-29c plays a tumor-suppressive role in various cancers, such as hepatocellular carcinoma (Wu et al., 2019), non-small-cell lung cancer (NSCLC) (Sun et al., 2018), endometrial cancer (EC) (Li et al., 2019), colon cancer (Chen et al., 2017), prostate cancer (Li et al., 2018), and pancreatic cancer (Lu et al., 2021). For instance, overexpression of miR-29c-3p in HCC can activate large tumor suppressor gene 1 (LATS1) expression and demethylation of LATS1 via DNMT3B, thereby impeding HCC cell tumor growth and migration both in vitro and in vivo (Wu et al., 2019). Intriguingly, this is in agreement with other results reported by SW Nam’s group (Bae et al., 2014). These data suggest that FANCD2 is possibly a direct target of miR-29c-3p in HCC.
To date, extensive evidence has shown that lncRNAs could act as a molecular sponge for miRNAs in cancers via the ceRNA mechanism to modulate ferroptosis. For example, lncRNAs HULC in HCC induced ATF4 stability by competitively binding to miR-3200-5p, resulting in the inhibition of ferroptosis (Guan et al., 2022). Among these lncRNAs, we found that nine possible lncRNAs in HCC were screened through the starBase database, and the high activity of both DUXAP8 and LINC00511 within HCC had poorer survival outcomes. Interestingly, we observed that both lncRNAs displayed a positive correlation with FANCD2 expression but were opposite to miR-29c-3p. Our result further supports previous observations that DUXAP8 (Liu et al., 2021) and LINC00511 (Hu et al., 2020) act as ceRNA for miR584-5p together with miR195, further contributing to HCC cell proliferation, invasion, and migration, and DUXAP8 has also been reported to be ferroptosis-associated lncRNA (Xing et al., 2021). Overall, these data uncover that this lncRNA functions as the molecular sponge in HCC impinged on ceRNA machinery.
Tumor-infiltrating immune cells have generally emerged as critical controllers of the immunological response, as they are surrounded by various infiltrating inflammatory cells that can be exploited to modulate ferroptosis-inducing resistance in tumor cells or mediate antitumor immunity (Whiteside, 2002; Wattenberg and Beatty, 2020; Zhao et al., 2020). For example, CD4+ T cells represent a major component of the adaptive immune system, establishing effective antitumor properties (Togashi et al., 2019; Oh and Fong, 2021). Indeed, we also found that FANCD2 in HCC was strongly associated with the high infiltration of CD4+ T cells. Apart from this, it harbors the same trend in the infiltration of neutrophils, purity, macrophages, CD8+ T cells, dendritic cells, and B cells. These data demonstrate that FANCD2 is an essential component of the adaptive immune system against HCC cells. Immunotherapy has been shown to be the backbone of effective treatment against cancer cells. Its blockade of the interaction between ICIs, such as CTLA-4, PD-1, and PD-L1 antibodies, and their ligands with IC inhibitors has been approved to treat several types of tumors, depending on the results of the trials (Sharma and Allison, 2020). Importantly, FANCD2 in HCC revealed a similar trend to the aforementioned three checkpoints, which has been successfully validated by IPS analysis, indicating that it might exert a pivotal role in restoring immune cytotoxic activity.
Although GSEA was explored to decipher the tumor-associated and even tumor-promoting roles of FANCD2, encompassing the cell cycle, p53 signaling pathway, mTOR signaling pathway, DNA replication, and TGF-beta signaling pathway, there are some restrictions of the current study that need to be highlighted. First, gathering more robust data or samples of HCC is needed to verify the tumor-promoting role of FANCD2. Furthermore, despite using comprehensive and deep in silico analysis, we provided few experiments with meaningful insights into the DUXAP8-miR-29c-FANCD2 and LINC00511-miR-29c-FANCD2 axes via the ceRNA network and the effect of FANCD2 on antitumor immunity.
In summary, this is the first study to report that ncRNA-regulated FANCD2 in HCC plays a tumor-promoting role based on the starBase database, HCCDB, TIMER database, and GEPIA database, indicating that it can be employed as a potential prognostic marker. In addition, our studies reveal that this ferroptosis-associated regulator harbors a strong association with HCC immunity, including immune infiltration cells, immune infiltration markers, and immune checkpoints, suggesting its significant immunological antitumor value. Thus, FANCD2-mediated ceRNA offers a promising strategy for novel targeted HCC therapies.
Data availability statement
Publicly available datasets were analyzed in this study. These data can be found at: (https://xena.ucsc.edu/) and (http://starbase.sysu.edu.cn/).
Author contributions
ZY designed the experiments and wrote the manuscript; YS and YL modified the language of the manuscript. YM and GD collected the data and provided references; HZ analyzed the data, conducted the experiments, and revised the manuscript; and BT provided financial support.
Acknowledgments
We sincerely acknowledge the contributions from TCGA project and starBase.
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.
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.955225/full#supplementary-material
Supplementary Figure S1 | Prediction of DUXAP8 and LINC00511 subcellular localization using iLoc-lncRNA.
Supplementary Figure S2 | Comprehensive analysis of FANCD2 in HCC. (A) Differential expression of FANCD2 across cancers according to the HCCDB. (B) Differential expression of FANCD2 in multiple HCC cell lines according to the CCLE database. (C,D) Kaplan‒Meier survival analysis of FANCD2 in HCC through the HCCDB.
References
Bae, H. J., Noh, J. H., Kim, J. K., Eun, J. W., Jung, K. H., Kim, M. G., et al. (2014). MicroRNA-29c functions as a tumor suppressor by direct targeting oncogenic SIRT1 in hepatocellular carcinoma. Oncogene 33 (20), 2557–2567. doi:10.1038/onc.2013.216
Barretina, J., Caponigro, G., Stransky, N., Venkatesan, K., Margolin, A. A., Kim, S., et al. (2019). Addendum: The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 565 (7738), E5–E6. doi:10.1038/s41586-018-0722-x
Beermann, J., Piccoli, M. T., Viereck, J., and Thum, T. (2016). Non-coding RNAs in development and disease: Background, mechanisms, and therapeutic approaches. Physiol. Rev. 96 (4), 1297–1325. doi:10.1152/physrev.00041.2015
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
Chen, G., Zhou, T., Li, Y., Yu, Z., and Sun, L. (2017). p53 target miR-29c-3p suppresses colon cancer cell invasion and migration through inhibition of PHLDB2. Biochem. Biophys. Res. Commun. 487 (1), 90–95. doi:10.1016/j.bbrc.2017.04.023
Chen, Q., Wang, W., Wu, Z., Chen, S., Chen, X., Zhuang, S., et al. (2021). Over-expression of lncRNA TMEM161B-AS1 promotes the malignant biological behavior of glioma cells and the resistance to temozolomide via up-regulating the expression of multiple ferroptosis-related genes by sponging hsa-miR-27a-3p. Cell. Death Discov. 7 (1), 311. doi:10.1038/s41420-021-00709-4
Cheng, A. L., Hsu, C., Chan, S. L., Choo, S. P., and Kudo, M. (2020). Challenges of combination therapy with immune checkpoint inhibitors for hepatocellular carcinoma. J. Hepatol. 72 (2), 307–319. doi:10.1016/j.jhep.2019.09.025
Cózar, B., Greppi, M., Carpentier, S., Narni-Mancinelli, E., Chiossone, L., and Vivier, E. (2021). Tumor-infiltrating natural killer cells. Cancer Discov. 11 (1), 34–44. doi:10.1158/2159-8290.Cd-20-0655
Dall'Olio, F. G., Marabelle, A., Caramella, C., Garcia, C., Aldea, M., Chaput, N., et al. (2022). Tumour burden and efficacy of immune-checkpoint inhibitors. Nat. Rev. Clin. Oncol. 19 (2), 75–90. doi:10.1038/s41571-021-00564-3
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
Doroshow, D. B., Bhalla, S., Beasley, M. B., Sholl, L. M., Kerr, K. M., Gnjatic, S., et al. (2021). PD-L1 as a biomarker of response to immune-checkpoint inhibitors. Nat. Rev. Clin. Oncol. 18 (6), 345–362. doi:10.1038/s41571-021-00473-5
Fagerholm, R., Sprott, K., Heikkinen, T., Bartkova, J., Heikkilä, P., Aittomäki, K., et al. (2013). Overabundant FANCD2, alone and combined with NQO1, is a sensitive marker of adverse prognosis in breast cancer. Ann. Oncol. 24 (11), 2780–2785. doi:10.1093/annonc/mdt290
Faivre, S., Rimassa, L., and Finn, R. S. (2020). Molecular therapies for HCC: Looking outside the box. J. Hepatol. 72 (2), 342–352. doi:10.1016/j.jhep.2019.09.010
Ferrer, M., Span, S. W., Vischioni, B., Oudejans, J. J., van Diest, P. J., de Winter, J. P., et al. (2005). FANCD2 expression in advanced non-small-cell lung cancer and response to platinum-based chemotherapy. Clin. Lung Cancer 6 (4), 250–254. doi:10.3816/CLC.2005.n.005
Forner, A., Reig, M., and Bruix, J. (2018). Hepatocellular carcinoma. Lancet 391 (10127), 1301–1314. doi:10.1016/s0140-6736(18)30010-2
Galluzzi, L., Humeau, J., Buqué, A., Zitvogel, L., and Kroemer, G. (2020). Immunostimulation with chemotherapy in the era of immune checkpoint inhibitors. Nat. Rev. Clin. Oncol. 17 (12), 725–741. doi:10.1038/s41571-020-0413-z
Goodall, G. J., and Wickramasinghe, V. O. (2021). RNA in cancer. Nat. Rev. Cancer 21 (1), 22–36. doi:10.1038/s41568-020-00306-0
Guan, L., Wang, F., Wang, M., Han, S., Cui, Z., Xi, S., et al. (2022). Downregulation of HULC induces ferroptosis in hepatocellular carcinoma via targeting of the miR-3200-5p/ATF4 Axis. Oxid. Med. Cell. Longev. 2022, 9613095. doi:10.1155/2022/9613095
Hu, W. Y., Wei, H. Y., Li, K. M., Wang, R. B., Xu, X. Q., and Feng, R. (2020). LINC00511 as a ceRNA promotes cell malignant behaviors and correlates with prognosis of hepatocellular carcinoma patients by modulating miR-195/EYA1 axis. Biomed. Pharmacother. 121, 109642. doi:10.1016/j.biopha.2019.109642
Islami, F., Miller, K. D., Siegel, R. L., Fedewa, S. A., Ward, E. M., and Jemal, A. (2017). Disparities in liver cancer occurrence in the United States by race/ethnicity and state. Ca. Cancer J. Clin. 67 (4), 273–289. doi:10.3322/caac.21402
Jhunjhunwala, S., Hammer, C., and Delamarre, L. (2021). Antigen presentation in cancer: Insights into tumour immunogenicity and immune evasion. Nat. Rev. Cancer 21 (5), 298–312. doi:10.1038/s41568-021-00339-z
Killock, D. (2017). Immunotherapy: Nivolumab keeps HCC in check and opens avenues for checkmate. Nat. Rev. Clin. Oncol. 14 (7), 392. doi:10.1038/nrclinonc.2017.70
Komatsu, H., Masuda, T., Iguchi, T., Nambara, S., Sato, K., Hu, Q., et al. (2017). Clinical significance of FANCD2 gene expression and its association with tumor progression in hepatocellular carcinoma. Anticancer Res. 37 (3), 1083–1090. doi:10.21873/anticanres.11420
Li, C., Zhang, Y., Liu, J., Kang, R., Klionsky, D. J., and Tang, D. (2021). Mitochondrial DNA stress triggers autophagy-dependent ferroptotic death. Autophagy 17 (4), 948–960. doi:10.1080/15548627.2020.1739447
Li, J., Fu, F., Wan, X., Huang, S., Wu, D., and Li, Y. (2018). Up-regulated miR-29c inhibits cell proliferation and glycolysis by inhibiting SLC2A3 expression in prostate cancer. Gene 665, 26–34. doi:10.1016/j.gene.2018.04.086
Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 42, D92–D97. doi:10.1093/nar/gkt1248
Li, L., Shou, H., Wang, Q., and Liu, S. (2019). Investigation of the potential theranostic role of KDM5B/miR-29c signaling axis in paclitaxel resistant endometrial carcinoma. Gene 694, 76–82. doi:10.1016/j.gene.2018.12.076
Li, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020). TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 48 (W1), W509–W514. doi:10.1093/nar/gkaa407
Lian, Q., Wang, S., Zhang, G., Wang, D., Luo, G., Tang, J., et al. (2018). Hccdb: A database of hepatocellular carcinoma expression Atlas. Genomics Proteomics Bioinforma. 16 (4), 269–275. doi:10.1016/j.gpb.2018.07.003
Lin, Y., Pan, X., and Shen, H. B. (2021). lncLocator 2.0: a cell-line-specific subcellular localization predictor for long non-coding RNAs with interpretable deep learning. Bioinformatics 37, 2308–2316. doi:10.1093/bioinformatics/btab127
Liu, T., Ghosal, G., Yuan, J., Chen, J., and Huang, J. (2010). FAN1 acts with FANCI-FANCD2 to promote DNA interstrand cross-link repair. Science 329 (5992), 693–696. doi:10.1126/science.1192656
Liu, Z., Lu, J., Fang, H., Sheng, J., Cui, M., Yang, Y., et al. (2021). m6A modification-mediated DUXAP8 regulation of malignant phenotype and chemotherapy resistance of hepatocellular carcinoma through miR-584-5p/MAPK1/ERK pathway Axis. Front. Cell. Dev. Biol. 9, 783385. doi:10.3389/fcell.2021.783385
Llovet, J. M., Kelley, R. K., Villanueva, A., Singal, A. G., Pikarsky, E., Roayaie, S., et al. (2021). Hepatocellular carcinoma. Nat. Rev. Dis. Prim. 7 (1), 6. doi:10.1038/s41572-020-00240-3
Lu, Y., Tang, L., Zhang, Z., Li, S., Liang, S., Ji, L., et al. (2021). Corrigendum to "long noncoding RNA TUG1/miR-29c Axis Affects cell proliferation, invasion, and migration in human pancreatic cancer. Dis. Markers 2021, 5150272. doi:10.1155/2021/5150272
Navarro Gonzalez, J., Zweig, A. S., Speir, M. L., Schmelter, D., Rosenbloom, K. R., Raney, B. J., et al. (2021). The UCSC genome browser database: 2021 update. Nucleic Acids Res. 49 (D1), D1046–d1057. doi:10.1093/nar/gkaa1070
Oh, D. Y., and Fong, L. (2021). Cytotoxic CD4(+) T cells in cancer: Expanding the immune effector toolbox. Immunity 54 (12), 2701–2711. doi:10.1016/j.immuni.2021.11.015
Ozawa, H., Iwatsuki, M., Mimori, K., Sato, T., Johansson, F., Toh, H., et al. (2010). FANCD2 mRNA overexpression is a bona fide indicator of lymph node metastasis in human colorectal cancer. Ann. Surg. Oncol. 17 (9), 2341–2348. doi:10.1245/s10434-010-1002-7
Palagyi, A., Neveling, K., Plinninger, U., Ziesch, A., Targosz, B. S., Denk, G. U., et al. (2010). Genetic inactivation of the Fanconi anemia gene FANCC identified in the hepatocellular carcinoma cell line HuH-7 confers sensitivity towards DNA-interstrand crosslinking agents. Mol. Cancer 9, 127. doi:10.1186/1476-4598-9-127
Pinter, M., Jain, R. K., and Duda, D. G. (2021). The current landscape of immune checkpoint blockade in hepatocellular carcinoma: A review. JAMA Oncol. 7 (1), 113–123. doi:10.1001/jamaoncol.2020.3381
Salmena, L., Poliseno, L., Tay, Y., Kats, L., and Pandolfi, P. P. (2011). A ceRNA hypothesis: The rosetta stone of a hidden RNA language? Cell. 146 (3), 353–358. doi:10.1016/j.cell.2011.07.014
Sangro, B., Chan, S. L., Meyer, T., Reig, M., El-Khoueiry, A., and Galle, P. R. (2020). Diagnosis and management of toxicities of immune checkpoint inhibitors in hepatocellular carcinoma. J. Hepatol. 72 (2), 320–341. doi:10.1016/j.jhep.2019.10.021
Semlow, D. R., and Walter, J. C. (2021). Mechanisms of vertebrate DNA interstrand cross-link repair. Annu. Rev. Biochem. 90, 107–135. doi:10.1146/annurev-biochem-080320-112510
Sharma, P., and Allison, J. P. (2020). Dissecting the mechanisms of immune checkpoint therapy. Nat. Rev. Immunol. 20 (2), 75–76. doi:10.1038/s41577-020-0275-8
Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2022). Cancer statistics, 2022. Ca. Cancer J. Clin. 72 (1), 7–33. doi:10.3322/caac.21708
Song, X., Xie, Y., Kang, R., Hou, W., Sun, X., Epperly, M. W., et al. (2016). FANCD2 protects against bone marrow injury from ferroptosis. Biochem. Biophys. Res. Commun. 480 (3), 443–449. doi:10.1016/j.bbrc.2016.10.068
Su, Z. D., Huang, Y., Zhang, Z. Y., Zhao, Y. W., Wang, D., Chen, W., et al. (2018). iLoc-lncRNA: predict the subcellular location of lncRNAs by incorporating octamer composition into general PseKNC. Bioinformatics 34 (24), 4196–4204. doi:10.1093/bioinformatics/bty508
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
Sun, D. M., Tang, B. F., Li, Z. X., Guo, H. B., Cheng, J. L., Song, P. P., et al. (2018). MiR-29c reduces the cisplatin resistance of non-small cell lung cancer cells by negatively regulating the PI3K/Akt pathway. Sci. Rep. 8 (1), 8007. doi:10.1038/s41598-018-26381-w
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Cagdoi:10.3322/caac.21660
Tang, Z., Kang, B., Li, C., Chen, T., and 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
Tapper, E. B., and Parikh, N. D. (2018). Mortality due to cirrhosis and liver cancer in the United States, 1999-2016: Observational study. Bmj 362, k2817. doi:10.1136/bmj.k2817
Tay, Y., Rinn, J., and Pandolfi, P. P. (2014). The multilayered complexity of ceRNA crosstalk and competition. Nature 505 (7483), 344–352. doi:10.1038/nature12986
Togashi, Y., Shitara, K., and Nishikawa, H. (2019). Regulatory T cells in cancer immunosuppression - implications for anticancer therapy. Nat. Rev. Clin. Oncol. 16 (6), 356–371. doi:10.1038/s41571-019-0175-7
Vogel, A., and Saborowski, A. (2022). Medical therapy of HCC. J. Hepatol. 76 (1), 208–210. doi:10.1016/j.jhep.2021.05.017
Wattenberg, M. M., and Beatty, G. L. (2020). Overcoming immunotherapeutic resistance by targeting the cancer inflammation cycle. Semin. Cancer Biol. 65, 38–50. doi:10.1016/j.semcancer.2020.01.002
Whiteside, T. L. (2002). Tumor-induced death of immune cells: Its mechanisms and consequences. Semin. Cancer Biol. 12 (1), 43–50. doi:10.1006/scbi.2001.0402
Wu, H., Zhang, W., Wu, Z., Liu, Y., Shi, Y., Gong, J., et al. (2019). miR-29c-3p regulates DNMT3B and LATS1 methylation to inhibit tumor progression in hepatocellular carcinoma. Cell. Death Dis. 10 (2), 48. doi:10.1038/s41419-018-1281-7
Xing, X. L., Yao, Z. Y., Ou, J., Xing, C., and Li, F. (2021). Development and validation of ferroptosis-related lncRNAs prognosis signatures in kidney renal clear cell carcinoma. Cancer Cell. Int. 21 (1), 591. doi:10.1186/s12935-021-02284-1
Zhao, Y., Li, M., Yao, X., Fei, Y., Lin, Z., Li, Z., et al. (2020). HCAR1/MCT1 regulates tumor ferroptosis through the lactate-mediated AMPK-SCD1 activity and its therapeutic implications. Cell. Rep. 33 (10), 108487. doi:10.1016/j.celrep.2020.108487
Keywords: HCC, ferroptosis, FANCD2, prognosis, ceRNA, immunotherapy
Citation: Yang Z, Song Y, Li Y, Mao Y, Du G, Tan B and Zhang H (2022) Integrative analyses of prognosis, tumor immunity, and ceRNA network of the ferroptosis-associated gene FANCD2 in hepatocellular carcinoma. Front. Genet. 13:955225. doi: 10.3389/fgene.2022.955225
Received: 28 May 2022; Accepted: 05 September 2022;
Published: 29 September 2022.
Edited by:
Feng Xu, Shantou University, ChinaReviewed by:
Fahui Liu, University of Gothenburg, SwedenGui Chengpeng, The First Affiliated Hospital of Sun Yat-sen University, China
Copyright © 2022 Yang, Song, Li, Mao, Du, Tan and Zhang. 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.
*Correspondence: Bangxian Tan, tbx_nsmc@126.com; Hongpan Zhang, hongpanzhangsccn@163.com
†These authors have contributed equally to this work