- Research
- Open access
- Published:
High FAAP24 expression reveals poor prognosis and an immunosuppressive microenvironment shaping in AML
Cancer Cell International volume 23, Article number: 117 (2023)
Abstract
Background
As a core member of the FA complex, in the Fanconi anemia pathway, FAAP24 plays an important role in DNA damage repair. However, the association between FAAP24 and patient prognosis in AML and immune infiltration remains unclear. The purpose of this study was to explore its expression characteristics, immune infiltration pattern, prognostic value and biological function using TCGA-AML and to verify it in the Beat AML cohort.
Methods
In this study, we examined the expression and prognostic value of FAAP24 across cancers using data from TCGA, TARGET, GTEx, and GEPIA2. To further investigate the prognosis in AML, development and validation of a nomogram containing FAAP24 were performed. GO/KEGG, ssGSEA, GSVA and xCell were utilized to explore the functional enrichment and immunological features of FAAP24 in AML. Drug sensitivity analysis used data from the CellMiner website, and the results were confirmed in vitro.
Results
Integrated analysis of the TCGA, TARGET and GTEx databases showed that FAAP24 is upregulated in AML; meanwhile, high FAAP24 expression was associated with poor prognosis according to GEPIA2. Gene set enrichment analysis revealed that FAAP24 is implicated in pathways involved in DNA damage repair, the cell cycle and cancer. Components of the immune microenvironment using xCell indicate that FAAP24 shapes an immunosuppressive tumor microenvironment (TME) in AML, which helps to promote AML progression. Drug sensitivity analysis showed a significant correlation between high FAAP24 expression and chelerythrine resistance. In conclusion, FAAP24 could serve as a novel prognostic biomarker and play an immunomodulatory role in AML.
Conclusions
In summary, FAAP24 is a promising prognostic biomarker in AML that requires further exploration and confirmation.
Introduction
AML is a heterogeneous hematological malignancy caused by malignant transformation of hematopoietic stem cells that is characterized by an increase in clonal myeloid primordial cells [1]. AML patients have a very variable prognosis and a high mortality rate: the 5-year survival rate is approximately 27% for people older than 20 years and approximately 69% for people younger than 20 years [2]. The pathogenesis and mechanism of AML are extremely complex and mainly include environmental factors and genetic factors.
Genetic mutations in cell proliferation, differentiation and apoptosis pathways are the basis for the pathogenesis of AML [3, 4]. Mutation also results in abnormal gene expression of oncogenes or suppressor genes. For example, mutation of the transcriptional regulatory factor TP53 upregulates the expression of genes in cancers [5, 6]. With the advent of personalized medicine, many diagnostic and prognostic biomarkers with disordered expression levels have been identified and contribute to a better understanding of the molecular pathology of the disease [7].
Forming a heterodimer complex with its partner FANCM, FAAP24 (FA-associated protein 24) can recognize and bind to damaged DNA and recruit other members of the FA pathway to initiate downstream repair [8]. Deficiency of FAAP24 causes a cancer-prone recessive genetic disorder characterized by congenital abnormalities, bone marrow failure, and cancer susceptibility [9]. According to structural analysis, FAAP24 has a high affinity with single-stranded DNA (ssDNA), which is required for optimized checkpoint activation. Its ssDNA-binding activity and FANCM-interacting functions enable FAAP24 to respond to DNA damage against crosslinking lesions [10]. The Fanconi anemia DNA repair pathway, in which FAAP24 involves, has recently been intensively studied for its contribution to anticancer drug resistance [11]. Therefore, functional research on FAAP24 has also been extended to other cancers, including human clone cancer [12]. However, its prognosis-impacting function in AML is rarely discussed.
Through an analysis of the public databases, we investigated FAAP24 expression differences across cancer types. We evaluated the independent prognostic value of FAAP24 in AML using multidimensional analysis. To further explore the functional network and pathogenetic mechanism related to FAAP24 in AML, we investigated its regulatory function in tumor immunity and its interaction with m6A RNA methylation and cuproptosis in AML. Our study fills the gap in discovering the biological function of FAAP24 and uncovers a promising prognostic biomarker for AML.
Methods and materials
The flow chart of this study is presented in Figure S1.
Sources of data and processing methods
To explore the expression of FAAP24 across tumors and normal controls, the gene expression levels of RNA sequencing (RNA-seq), including 33 types of cancer and 24 paracancerous tissues from TCGA (The Cancer Genome Atlas), were downloaded from the UCSC XENA website (https://xena.ucsc.edu/) [13]. A cohort combined with TCGA AML (n = 173), GTEx (Genotype-Tissue Expression, n = 337) and TARGET (Therapeutically Applicable Research To Generate Effective Treatments) AML (n = 196) was also collected to evaluate the differential expression of FAAP24 between AML and normal whole blood controls. Subsequently, GSE65409, consisting of 8 normal controls and 30 AML samples, was used to verify the difference in FAAP24 expression. The cohort of GDC TCGA AML (n = 151) was chosen to explore the prognostic value of FAAP24 in AML. Moreover, we downloaded gene counts of the Beat AML cohort from the data viewer (www.vizome.org), and the summary of clinical data was acquired from a previous study [14, 15]. The toil method was performed to process the expression profile data, and all RNA-seq data were normalized as log2 (TPM + 1) before analysis. Then, a total of 446 Beat AML patients were investigated to determine the relationship of FAAP24 expression with clinical features and validate the prognostic value of FAAP24 in AML. Processed datasets of composite expression and compound activity (DTP NCI-60) were also downloaded to perform drug sensitivity analysis (https://discover.nci.nih.gov/cellminer/home.do) [16].
Expression and survival analysis of FAAP24 across cancers
To explore FAAP24 expression across cancers, the Mann-Whitney U test was executed in each group. Box plots were used to exhibit the differential expression between tumors and normal tissues.
GEPIA2 (http://gepia2.cancer-pku.cn/) [17] was exploited to investigate the prognostic values of FAAP24 in different types of tumors. Based on the median expression of FAAP24 in different tumors, patients with each cancer were divided into two groups (FAAP24high and FAAP24low). The Kaplan -Meier method with the log-rank test, accompanied by the Cox proportional hazards model, was applied to explore and visualize the prognostic values of FAAP24 expression for overall survival (OS) across cancers.
Development and validation of a nomogram containing FAAP24 for AML prognosis
In the GDC TGCA AML cohort (n = 151), we performed univariate and multivariable survival analyses to investigate the prognostic role of the FAAP24 gene. Subsequently, a nomogram was established for predicting AML prognosis using significant variables in the multivariable Cox survival model (P < 0.05). Receiver operator characteristic (ROC) and calibration curves at different time points were sketched to estimate the discrimination and calibration ability. Decision curve analysis (DCA) was also plotted to evaluate the clinical utility. Meanwhile, we used Beat AML to prove FAAP24’s impact on AML prognosis and the applicability of the nomogram. According to the median risk score of the nomogram, both TCGA and Beat AML were divided into two risk groups (high risk and low risk). Kaplan-Meier curves with hazard ratios (HRs) were then generated to demonstrate the nomogram’s efficiency. Finally, a Sankey diagram was created to display the ability of the nomogram to reclassify risk stratification from cytogenetics. All analyses were carried out by using the R packages “survival”, “survminer”, “rms”, “timeROC”, “dcurves”, “ggalluvial” and “ggplot2”.
Analysis of FAAP24-related expressed genes
According to the median expression of FAAP24, AML patients from TCGA were segregated into two groups (FAAP24high and FAAP24low). The R package limma was adopted to search the differentially expressed genes (DEGs) between FAAP24high and FAAP24low patients. Pearson correlation analysis was used to screen the coexpressed genes of FAAP24 in AML by using the LinkedOmics database (http://www.linkedomics.org/) [18]. Both DEGs and coexpressed genes were visualized using a volcano map, and then a Venn diagram was drawn to explore the intersecting genes between them for functional enrichment analysis using the R packages VennDiagram and ggplot2.
Functional enrichment and immunological features in AML
Enrichment analyses of intersected genes, including GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, were carried out using the “clusteProfiler” R package. The hallmarks of tumors are of great significance for cancer research, and the enrichment in hallmark gene sets to a special phenotype may help to explain the mechanism related to FAAP24 in AML. “H.all.v2022.1.Hs.symbols.gmt” was downloaded from the MSigDB website (http://software.broadinstitute.org/gsea/msigdb) as a reference gene set. Gene set enrichment analysis (ssGSEA) was executed based on a single sample gene in TCGA AML by means of the “GSVA” R package. Pearson coefficients were calculated to explore the potential phenotype of cancer linked to FAAP24 expression.
The immune cycle in cancer reflects the anticancer immune response and consists of 7 steps. When these steps are activated to varying degrees, the destiny of cancer cells is determined at the same time [19]. We also performed ssGSEA to evaluate the anticancer immune cycle in TCGA AML based on individual gene expression. Then, Pearson correlation analyses were performed to scan the relationship between FAAP24 and the enrichment of these steps of the anticancer immunity cycle.
Afterwards, we further calculated the immunological features of TCGA AML patients in the tumor microenvironment (TME) from individual gene expression using the R package “xCell” [20]. Spearman coefficients were determined first to identify FAAP24-related immune cells. Thereafter, patients were also split into two groups (FAAP24high and FAAP24low) on account of the median expression, and the different abundances of immune cells were subsequently compared between these two groups through the Wilcoxon test. To explore whether FAAP24 expression is modified by m6A and interacts with cuproptosis in AML, Spearman coefficients were also used to investigate the relationship of FAAP24 with 23 m6A- and 10 cuproptosis-related genes [21, 22]. We visualized the results by using the R package “ggplot2”.
Patient samples and Western blotting
We obtained whole blood samples from 5 healthy controls and bone marrow samples from 4 AML patients. The test of each sample was carried out in triplicate. Western blotting was performed as described previously [8]. The FAAP24 primary antibody (GTX117277, United States) was purchased from GeneTex, and the GAPDH primary antibody (ab8245, United States) was obtained from Abcam.
Cell Culture
AML line cells were purchased from the American Type Culture Collection (ATCC). The cells were cultured in 90% RPMI 1640 medium (Gibco, Thermo Fisher Scientific) with 10% fetal bovine serum (FBS), 100 U/mL penicillin and 100 mg/mL streptomycin at 37 °C in a humidified 5% CO2 incubator. The cells were harvested, washed, and resuspended in phosphate buffered saline (PBS).
Real-time PCR
Four AML cell lines, HL-60, MV4-11, MOLM13 and U937, were harvested for expression analysis by quantitative PCR (qPCR) as previously described [23]. The following primers for FAAP24 were used: FAAP24forwards, 5’- ATGGCTTGACACCAGACTTTT-3’, and reverse, 5’-TACTGGGCTCTTTGGTTTGC-3’. Relative levels of expression were computed using the 2−ΔΔCt method with FAAP24 expression in HL60 for normalization. The samples were run in triplicate parallel reactions.
Drug sensitivity analysis and cell viability assay
Both expression levels and compound activity (DTP NCI-60, Z score) were downloaded in a processed format from the CellMiner website, and Pearson correlation analysis with a condition of |cor|>0.4 and P < 0.01 was performed to screen drug sensitivity associated with FAAP24 expression.
Then, these four AML cell lines were used to confirm drug sensitivity in vitro. When cells were in the logarithmic growth period, AML cells were seeded in 96-well plates at a density of 6000 cells per well. After treatment with chelerythrine (MCE company, 0, 1, 10 µM) for 72 h, the cells were incubated with 10 µL of CCK-8 reagent (Dojindo) for 2 h at 37 °C. In each well, the optical density (OD) was measured using a microplate reader at 450 nm, and the OD values were recorded as the means ± SDs. Each assay was performed in triplicate. IC50 values were analysed in GraphPad Prism 5 software.
Immunofluorescence staining
The distribution of FAAP24 in HL-60, MV4-11, MOLM13 and U937 cells was determined as previously described [24]. Cells were fixed first with 2% paraformaldehyde at room temperature for 10 min. After blocking with 5% normal goat serum, cells were incubated with FAAP24 antibody (Gene Tex) at 4 °C for 1 h and then incubated with secondary FITC-conjugated anti-rabbit Abs at 4 °C for 30 min. After washing with PBS, 100 ng/mL DAPI (Invitrogen) was added for nuclear staining at room temperature. Imaging was presented under a Zeiss LSM900 laser scanning confocal microscope (Zeiss, Germany).
Statistical analysis
R software (v 4.2.1) was used to conduct statistical analysis. χ2 or fish tests were performed to analyse the differences in clinical features between the FAAP24high and FAAP24low groups, including age (> 60 vs. ≤60), sex, cytogenetic risk, classification of French-American-British (FAB), and gene mutations. The Mann-Whitney U test was performed to compare differences in continuous variables, such as white blood cell (WBC) count, percentage of blasts in peripheral blood (PB) or bone marrow (BM), hemoglobin (HB) and platelets (PLT). P values less than 0.05 were considered statistically significant.
Results
Expression and prognostic analysis of FAAP24 across cancers
In pancancer analysis, the mRNA levels of FAAP24 in 33 types of tumors were explored and are presented in Fig. 1A and Supplementary Table S1. The results showed that almost all cancers expressed higher FAAP24 than normal tissues, including BLCA (bladder urothelial carcinoma), BRAC (breast cancer), and CESC (cervical cancer). Similarly, both TCGA AML (n = 173) and TARGET AML (n = 196) expressed higher FAAP24 than normal controls (whole blood, n = 337) from the GTEx database (Fig. 1B). This differential expression was also confirmed in the GSE65409 cohort (Fig. 1C) and in our western blot results (Fig. 1D). These expressive features suggest that FAAP24 may play a unique role in AML.
To analyse the effects of FAAP24 on survival across cancers, patients with 33 types of tumors were split into two different risk groups based on the median value of FAAP24 expression and analysed using the GEPIA2 database (Fig. 2A). The results showed that upregulation of FAAP24 had poor prognostic value in KICH (kidney chromophore, Fig. 2B), KIRC (kidney renal clear cell carcinoma, Fig. 2C), AML (Fig. 2D), LGG (brain lower grade glioma, Fig. 2E), LUAD (lung adenocarcinoma, Fig. 2F), and PAAD (pancreatic adenocarcinoma Fig. 2G) (P < 0.05).
Independent prognostic value of FAAP24 in AML
Survival curves with log-rank tests and Cox proportional analyses inferred that high FAAP24 expression presented a significantly inferior OS than low FAAP24 expression in AML from the GEPIA2 database. Subsequently, GDC TCGA AML (n = 151) was further used to explore the relationship of FAAP24 with clinical features and confirm the prognostic role of FAAP24 in AML.
The clinical characteristics between the FAAP24high and FAAP24low groups are summarized in Table 1. In the results, more M5 and male patients were assigned to the FAAP24high group. Univariate Cox proportional analysis showed that FAAP24 expression (high vs. low, P = 0.002), age (> 60 vs. ≤60, P < 0.001) and cytogenetic risk (intermediate vs. favorable, P = 0.002; poor vs. favorable, P < 0.001) significantly affected the OS of AML (Supplementary Table S2). An independent prognostic role for FAAP24 expression in AML was confirmed by a multivariate model (HR = 1.75, P = 0.013, Fig. 3A and Supplementary Table S2).
Development of a predicted nomogram containing FAAP24 in TCGA AML
A nomogram (Fig. 3B) was developed for predicting AML OS using a multivariate Cox regression model. Evidence of age older than 60 was assigned 100 points, while the score was 0 points if age ≤ 60. The favorable risk of cytogenetics was regarded as the reference and assigned a 0 point, while 76 and 88 points were given for intermediate and poor risk of cytogenetics, respectively. Similarly, zero points was assigned to the low expression of FAAP24, while fifty-three points was given for the evidence of high FAAP24 expression. The prognostic score for each patient was determined by the sum of the weight integral of each variable. According to the median prognostic risk score of the nomogram, patients with high prognostic scores had inferior OS than patients with low risk scores (HR = 3.35, P = 0.002, Fig. 3C). Calibration plots showed good performances at 1, 3 and 5 years (Fig. 3D), and the ROC curves also demonstrated good discrimination (Fig. 3E). At 1, 3 and 5 years, these values are 0.783, 0.803 and 0.906, respectively.
When compared with the clinical model combined with age and cytogenetic risk only, DCA showed that the nomogram presented various benefits at 5 years (Fig. 3F). Moreover, the prognosis of AML predicted by the nomogram exhibited more benefits than either the intervention-all or intervention-none strategy at the different threshold points. In addition, the cytogenetic risk of AML could be further reclassified on the basis of the nomogram (Fig. 3G).
Validation of the nomogram containing FAAP24 in beat AML
To confirm the poor prognostic value of FAAP24 in AML, an independent AML cohort (n = 446) was downloaded and analysed, and its clinical features are presented in Supplementary Table S3. Univariate Cox proportional analysis showed that FAAP24 expression (high vs. low, P = 0.003), age (> 60 vs. ≤60, P < 0.001), cytogenetic risk (intermediate vs. favorable/poor vs. favorable, P < 0.001), and risk recommendation of ELN (European LeukemiaNet) 2017 (Intermediate or Adverse vs. Favorable, P < 0.001) were significantly related to OS in AML (Supplementary Table S4). Subgroup survival analysis (Fig. 4A) implied that FAAP24 may play an important role in patients with age > 60 (P = 0.002), cytogenetic risk of intermediate (P = 0.006), Intermediate or Adverse of ELN 2017 risk recommendation (P < 0.001), patients with NPM1neg mutation (P = 0.002), and patients with TP53 mutation (P = 0.031). An independent prognostic value of FAAP24 was also confirmed in Beat AML (HR = 1.67, P < 0.001, Fig. 3B and Supplementary Table S5).
In Beat AML, calibration plots of the external validation verified the practicability of the predicted nomogram containing FAAP24 at different times (Fig. 4C, D and E). An inferior OS was presented in high-risk patients (HR = 2.48, P < 0.001, Fig. 4F). The ROC curves of external validation also confirmed good discrimination (Fig. 4G), with values of 0.697, 0.713 and 0.868 at 1, 3 and 5 years, respectively. Similar to TCGA AML results, reclassification of cytogenetic risk could be performed based on the predicted nomogram (Fig. 3H).
Analysis of FAAP24-related genes in AML
To better explore the underlying pathogenic mechanism of FAAP24 in AML, the coexpressed genes of FAAP24 were screened out first in TCGA-AML patients using the LinkedOmics database (Fig. 5A). A total of 586 genes were positively coexpressed with FAAP24 in AML, while 40 genes were negatively coexpressed with FAAP24 (Supplementary Table S6, |cor|>0.3, FDR < 0.05). Then, we also identified 2826 DEGs between FAAP24high and FAAP24low in AML (Fig. 5B, Supplementary Table S7), including 1229 upregulated (log2 FC > 0.3, P < 0.05) and 1597 downregulated protein-coding genes (log2 FC<-0.3, P < 0.05). Subsequently, a total of 217 FAAP24-related genes overlapped between coexpressed genes and DEGs of FAAP24 (Fig. 5C), consisting of 207 upregulated genes and 10 downregulated genes, which were included for further functional enrichment analysis.
Biological function of FAAP24 in AML
To investigate the potential biological function of FAAP24 in AML, we performed GO/KEGG analyses on the 217 overlapping genes (Supplementary Table S8) using bubble charts to present the top enriched items. The main sets of BP (biological process) were organelle fission, nuclear division and chromosome segregation (Fig. 5D), and chromosomal region was the major list of CC (cellular component) (Fig. 5E). The molecular functions (MF) of FAAP24 and its partner genes were primarily involved in catalytic activity that acts to modify DNA (Fig. 5F).
KEGG analyses (Fig. 5G) suggested that FAAP24 and its related genes were mainly involved in DNA damage and repair pathways. In addition to the cell cycle, DNA replication and cellular senescence, these genes were also involved in regulating the metabolic processes of AML, such as pyruvate, carbon, pyrimidine, glyoxylate and dicarboxylate metabolism. These results implied the potential biological role of FAAP24 in leukemogenesis.
Moreover, we also explored the correlation between FAAP24 and the enrichment scores of the tumor pathway from the hallmark gene set (Fig. 6A, Supplementary Table S9). FAAP24 was positively associated with proliferation and the cell cycle, including E2F TARGETS (r = 0.452), FATTY ACID METABOLISM (r = 0.436), PEROXISOME (r = 0.424), G2M CHECKPOINT (r = 0.376), DNA REPAIR (r = 0.363), MTORC1 SIGNALING (r = 0.358) and MYC TARGETS V1 (r = 0.346), but negatively associated with TGF BETA SIGNALING (r=-0.285).
FAAP24 shapes an immunosuppressive TME in AML
To investigate the role of FAAP24 in the immune regulation of AML, the relationship between FAAP24 and the antitumour immune response is presented in Fig. 6B and Supplementary Table S10. Different steps of the anticancer immune cycle showed that higher FAAP24 expression in AML strongly correlated with a stronger recruitment of MDSCs (r = 0.314), TH2 cells (r = 0.239) and Treg cells (r = 0.175) but was negatively related to the processes of killing cancer cells (r=-0.226) and priming of antitumour immune response (r=-0.200). All these discoveries supported why high FAAP24 expression is often related to poor prognosis in AML.
Because of the adverse antitumour immune response, xCell was also used to further evaluate the TME induced by FAAP24 in AML. The results demonstrated that FAAP24 was positively (Fig. 7A) linked to Th2 cells (r = 0.299) and M2 macrophages (r = 0.171) but negatively (Fig. 7B) related to CD4 + naive T cells (r=-0.220), CD8 + Tcm cells (r=-0.220), CD4 + memory T cells (r=-0.213), CD4 + Tcm cells (r=-0.194), CD4 + Tem cells (r=-0.178) and CD8 + T cells (r=-0.171). When compared with the FAAP24low group, the FAAP24high group displayed low-level infiltration of immune cells, including many subtypes of CD8 + and CD4 + cells (Fig. 7C), but more Th2 cells and macrophages (P < 0.05). Although not statistically significant, the FAAP24high group still presented a higher tendency of M2 macrophage cells than the FAAP24low group (P = 0.08). These findings suggested that FAAP24 shapes an immunosuppressive TME in AML, which may help to promote leukemia progression.
FAAP24 interacted with m6A RNA methylation and cuproptosis in AML
To determine whether FAAP24 expression is mediated by m6A modification, the relationship of FAAP24 with m6A-related genes was probed in AML. In Fig. 8A, the expression level of FAAP24 was positively related to 13 m6A-related genes, including METTL14 (r = 0.254), METTL16 (r = 0.192), VIRMA (r = 0.207), ZC3H13 (r = 0.290), RBM15 (r = 0.202), RBM15B (r = 0.291), YTHDF1 (r = 0.163), YTHDF2 (r = 0.316), HNRNPC (r = 0.207), LRPPRC (r = 0.328), HNRNPA2B1 (r = 0.291), RBMX (r = 0.274), and ALKBH5 (r = 0.360). Moreover, 8 of 10 cuproptosis-related genes were significantly associated with FAAP24 (Fig. 8B): LIAS (r = 0.343), LIPT1 (r = 0.193), DLD (r = 0.368), DLAT (r = 0.307), PDHA1 (r = 0.387), PDHB (r = 0.368), MTF1 (r = 0.291), and CDKN2A (r = 0.185). These conclusions implied that FAAP24 is associated with m6A modification and cuproptosis in AML.
High FAAP24 expression is associated with chelerythrine resistance
Using the CellMiner database, drug sensitivity analysis showed that chelerythrine was meaningfully correlated with FAAP24 expression (Fig. 9A, r = 0.443, P < 0.001). Then, 4 AML cell lines were used to confirm the correlation of FAAP24 with chelerythrine sensitivity in vitro. In these 4 cell lines, FAAP24 mainly localized to the nucleus and was visualized using fluorescent labelling (Fig. 9B). The qPCR results presented a gradual increasing trend in FAAP24 expression in HL-60, MV4-11, MOLM13 and U937 cells (Fig. 9C). Drug sensitivity showed that high FAAP24 expression was associated with chelerythrine resistance, with IC50 values of 0.679 µM, 2.490 µM, 4.332 µM, and 7.024 µM, respectively (Fig. 9D).
Discussion
AML is a highly heterogeneous stem/progenitor cell disease that results from chromosome translocation and mutation of hematopoietic proliferation- and differentiation-related genes. Recent studies have indicated that some genetic abnormalities and specific chromosomal translocations in AML are associated with deficiencies in DNA repair pathways [25]. Among the canonical DNA damage repair pathways, the Fanconi anemia pathway is a new player that intersects with many other repair processes to respond to interstrand crosslink DNA lesions and thus to maintain chromosome stability [26]. To date, at least 22 FA complementation groups and many other interacting proteins have been identified in the FA pathway. Although it has been discovered that genetic mutations or abnormal expression in components of the FA pathway play important roles in the development of AML, treatment strategies targeting this pathway have made little progress.
As a component of the FA core complex, FAAP24 functions in recognizing and binding to damaged DNA and stalled replication folks. According to the GTEx and TCGA databases, high expression of FAAP24 in AML may promote the development of AML and thus has the potential to be a novel prognostic biomarker. It is yet to be determined whether FAAP24 acts as a tumor suppressor gene and how it works.
We further analysed the association between FAAP24 and OS across cancers using GEPIA2, and FAAP24 expression was significantly associated with survival in patients with six cancers. To verify the prognostic value of FAAP24, we analysed its relationship with prognosis in the TCGA AML and Beat AML cohorts. The results revealed significant prognostic values of FAAP24 expression in subgroups, including age, cytogenetic risk, NPM1neg mutation and TP53 mutation. All of these results confirmed the independent prognostic value of FAAP24 expression for OS in AML. It is worth noting that, as we mentioned before, FAAP24 is expressed at high levels in AML, and the prognosis is poor in the FAAP24high group. This phenomenon might be caused by the DNA damage repair function of this gene. DNA damage caused by anticancer drugs is quickly recognized and repaired by the overexpression of FAAP24, leading to drug resistance.
The enrichment results confirmed that FAAP24 was closely linked to the DNA damage repair process and cell cycle-related process in AML. The mechanisms of DNA repair proteins are often unclear and can be altered during the development of tumors. FANCM-FAAP24 is important for ATR-mediated checkpoint signaling in response to replication stress [27]. These results suggest that abnormal expression of FAAP24 may impact the normal proliferation, differentiation, cell apoptosis and hematopoietic function of bone marrow. The role of genetic and epigenetic changes in regulating cancer development and immune tolerance must be further explored, as well as the protein-protein interactions of FAAP24.
Many immunosuppressive factors are present in the TME, which is also an important factor hindering the benefits of cancer treatment. This immunosuppressive effect is attributed to many substances produced by tumors and immune cells. Additionally, the TME is packed with Tregs, MDSCs, TAMs and TANs and promotes tumor survival through the secretion of TGF-β, IL-10, nitric acid, and IDO (indoleamine 2,3 dioxygenase) [28]. The relationship between FAAP24 and antitumour immune cells in AML was explored. The results showed that higher FAAP24 expression in AML strongly correlated with a stronger recruitment of MDSCs, TH2 cells and Treg cells but was negatively related to the process of killing cancer cells and priming of immune response. The proliferation and differentiation of MDSCs into TAMs is induced by AML blasts, thus further inhibiting their immunogenicity [29]. On the other hand, increased Treg cells mediate T-cell suppression. Meanwhile, the high expression of FAAP24 is positively related to M2 macrophages, which are immunosuppressive and promote tumor progression by facilitating angiogenesis and tissue phagocytosis [30]. All these factors lead to the promotion of leukemia progression, which explains the poor prognosis of high FAAP24 expression in AML in our study. Further exploration of the strategies of tumor immune escape and the regulatory factors of M2 macrophages that are related to FAAP24 in AML would provide immunologic insights into novel treatments.
Considering the close correlation between FAAP24 and m6A methylation and cuproptosis-related genes in our results, further experiments could be performed to verify their interactions and their mechanism of promoting the development of AML. For example, RBM15, which modifies m6A targets in the Notch signaling pathway, is highly expressed in AML [31]. Other m6A genes that have been found to play oncogenic roles in AML include WTAP, METTL3, METTL14, FTO, ALKBH5, IGF2BP1, and PRMT6 [32,33,34]. Among them, METTL14, RBM15 and ALKBH5 were significantly associated with FAAP24, which enhanced our understanding of the mechanism of FAAP24 in AML. The role of cuproptosis in AML is rarely studied. It was reported that 3 cuproptosis-related genes, GCSH, LIPT1, and DLAT, have independent and significant prognostic value in AML [35]. In our research, LIPT1 and DLAT were significantly associated with FAAP24. Our research provides evidence that m6A methylation and cuproptosis are related to AML and that FAAP24 might play an important role in these processes. Further investigation should be performed to explore the possible molecular pathways.
By retrieving data from the NCI-60 cell line and 4 AML cell lines, our study illustrated that increased FAAP24 expression levels were associated with chelerythrine resistance. Chelerythrine is reported to inhibit proliferation and induce apoptosis through several signaling pathways [36]. Although it is currently the focal point in antitumour research and its effect has been approved in prostate cancer, cervical cancer, gastric cancer cells, non-small cell lung carcinoma, renal cell carcinoma, etc. [37], little research has been performed on AML. The findings of our study provide new strategies for treating AML, and further studies are needed to confirm the pharmacological mechanism of chelerythrine in this disease.
Our study aimed to investigate the prognostic value of FAAP24 in AML, its relationship with tumor-infiltrating immune cells, and its associated pathways. It was also confirmed that FAAP24 expression was associated with poor prognosis and suppressed TME in AML. Finally, drug sensitivity analyses suggested a correlation between FAAP24 expression and chelerythrine resistance, which could be a novel treatment strategy. The results validated the value of FAAP24 as a prognostic biomarker in AML and identified significant areas for further research.
Conclusions
In summary, we identified significant areas for further exploration and confirmation of the value of FAAP24 as a prognostic biomarker in AML.
Data Availability
Related datasets of TCGA, TARGET, GTEx, GSE65409, Beat AML and DTP NCI-60 could be downloaded from the UCSC XENA website, cBioPortal website and NIH websites. The data that support the findings of this study are available from the corresponding author upon reasonable request.
Abbreviations
- AML:
-
Acute myeloid leukemia
- FAAP24:
-
FA-associated protein 24
- TP53:
-
Tumor protein p53
- FANCM:
-
Fanconi anemia complementation group M
- FA:
-
Fanconi Anemia
- ssDNA:
-
single-stranded DNA
- TCGA:
-
The Cancer Genome Atlas
- GTEx:
-
Genotype-Tissue Expression
- GEPIA:
-
Gene Expression Profiling Interactive Analysis
- GDC:
-
Genomic Data Commons
- ROC:
-
Receiver operator characteristic
- DCA:
-
Decision curve analysis
- HRs:
-
Hazard ratios
- DEGs:
-
expressed genes
- GO:
-
Gene Ontology
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- TME:
-
Tumor microenvironment
- ssGSEA:
-
Single sample gene set enrichment analysis
- FITC:
-
Fluorescein Isothiocyanate
- DAPI:
-
4’,6-Diamidino-2-phenylindole
- ACC:
-
adrenocortical cancer
- BLCA:
-
bladder urothelial carcinoma
- BRAC:
-
breast cancer
- CESC:
-
cervical cancer
- KICH:
-
Kidney Chromophobe
- KIRC:
-
Kidney renal clear cell carcinoma
- LGG:
-
Brain Lower Grade Glioma
- LUAD:
-
Lung adenocarcinoma
- PAAD:
-
Pancreatic adenocarcinoma
- OS:
-
Overall survival
- ELN:
-
European Leukemia Net
- BP:
-
biological process
- CC:
-
cellular component
- MF:
-
molecular functions
- Tcm:
-
Central memory T-cell
- Th2:
-
T helper 2 cells
- Tem:
-
Effective Memory T-Cell
- RT‒qPCR:
-
Quantitative real-time polymerase chain reaction
- IC50:
-
Half maximal inhibitory concentration
References
Rubnitz JE, Gibson B, Smith FO. Acute myeloid leukemia. Hematol Oncol Clin N Am. 2010;24(1):35–63.
Welsh J. Acute Myeloid Leukemia Survival Rates and Outlook 2021. Available from: https://www.verywellhealth.com/acute-myeloid-leukemia-survival-rates-5208865.
Steffen B, Müller-Tidow C, Schwäble J, Berdel WE, Serve H. The molecular pathogenesis of acute myeloid leukemia. Crit Rev Oncol/Hematol. 2005;56(2):195–221.
Tan Y, Zheng L, Du Y, Zhong Q, Zhu Y, Liu ZML, Shuang, et al. Identification of the hub genes and pathways involved in acute myeloid leukemia using bioinformatics analysis. Medicine. 2020;99(35):e22047.
Donehower L, Soussi T, Korkut A, Liu Y, Schultz A, Cardenas M et al. Integrated Analysis of TP53 Gene and Pathway Alterations in The Cancer Genome Atlas. Cell Reports. 2019;28:1370-84.e5.
Wang C, Zhang S, Ma B, Fu Y, Luo Y. TP53 mutations upregulate RCP expression via Sp1/3 to drive lung cancer progression. Oncogene. 2022;41.
Prada J, Arroyave Ospina J, Rothlisberger S. Molecular Biomarkers in Acute Myeloid Leukemia. Blood Rev. 2016;31.
Ciccia A, Ling C, Coulthard R, Yan Z, Xue Y, Meetei AR, et al. Identification of FAAP24, a Fanconi Anemia Core Complex protein that interacts with FANCM. Mol Cell. 2007;25(3):331–43.
de Winter JP, Joenje H. Fanconi Anemia. In: Maloy S, Hughes K, editors. Brenner’s Encyclopedia of Genetics (Second Edition). San Diego: Academic Press; 2013. p. 17–20.
Wang Y, Han X, Wu F, Leung JW, Lowery MG, Do H, et al. Structure analysis of FAAP24 reveals single-stranded DNA-binding activity and domain functions in DNA damage response. Cell Res. 2013;23(10):1215–28.
Ghosal K, Agatemor C, Han RI, Ku ATT, Sabu MS. Fanconi Anemia DNA Repair Pathway as a New Mechanism to Exploit Cancer Drug Resistance. Mini Rev Med Chem. 2020;20(9):pp. 779 – 87(9).
Wang Y, Leung Justin W, Jiang Y, Lowery Megan G, Do H, Vasquez Karen M, et al. FANCM and FAAP24 maintain Genome Stability via Cooperative as Well as Unique Functions. Mol Cell. 2013;49(5):997–1009.
Goldman MJ, Craft B, Hastie M, Repečka K, McDade F, Kamath A, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38(6):675–8.
Tyner JW, Tognon CE, Bottomly Dea. Functional genomic landscape of acute myeloid leukaemia. Nature. 2018;562(7728):526–31.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of Complex Cancer Genomics and Clinical Profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1–pl.
Reinhold WC, Sunshine M, Liu H, Varma S, Kohn KW, Morris J, et al. CellMiner: a web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 Cell Line Set. Cancer Res. 2012;72(14):3499–511.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47:W556 - W60.
Vasaikar S, Straub P, Wang J, Zhang B, LinkedOmics. Analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2017;46.
Chen DS, Mellman I. Oncology meets immunology: the cancer-immunity cycle. Immunity. 2013;39 1:1–10.
Aran D, Hu Z, Butte AJ. xCell: Digitally portraying the tissue cellular heterogeneity landscape. bioRxiv. 2017:114165.
Li Y, Xiao J, Bai J, et al. Molecular characterization and clinical relevance of m6A regulators across 33 cancer types. Mol Cancer. 2019;18:137. https://doi.org/10.1186/s12943-019-1066-322.
Tsvetkov P, Coy S, Petrova B, Dreishpoon M, Verma A, Abdusamad M, et al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Sci (New York NY). 2022;375(6586):1254–61.
Xu X, Liu R, Zhou X, Zhang Z, Zhu T, Huang Y et al. Characterization of exosomes derived from IPEC-J2 treated with probiotic Bacillus amyloliquefaciens SC06 and its regulation of macrophage functions. Front Immunol. 2022;13.
Xu X-g, Hu J-f, Ma J-x, Nie L, Shao T, Xiang L-x, et al. Essential roles of TIM-1 and TIM-4 homologs in adaptive humoral immunity in a zebrafish model. J Immunol. 2016;196(4):1686–99.
Esposito MT, So CWE. DNA damage accumulation and repair defects in acute myeloid leukemia: implications for pathogenesis, disease progression, and chemotherapy resistance. Chromosoma. 2014;123(6):545–61.
Ceccaldi R, Sarangi P, D’Andrea AD. The fanconi anaemia pathway: new players and new functions. Nat Rev Mol Cell Biol. 2016;17(6):337–49.
Hořejší Z, Collis SJ, Boulton SJ. FANCM-FAAP24 and HCLK2: roles in ATR signalling and the Fanconi Anemia pathway. Cell Cycle. 2009;8(8):1133–7.
Sendker S, Reinhardt D, Niktoreh N. Redirecting the Immune Microenvironment in Acute Myeloid Leukemia. Cancers [Internet]. 2021; 13(6).
Hyun SY, Na EJ, Jang JE, Chung H, Kim SJ, Kim JS, et al. Immunosuppressive role of CD11b + CD33 + HLA-DR – myeloid-derived suppressor cells-like blast subpopulation in acute myeloid leukemia. Cancer Med. 2020;9(19):7007–17.
Sumitomo R, Hirai T, Fujita M, Murakami H, Otake Y, Huang CL. M2 tumor–associated macrophages promote tumor progression in non–small–cell lung cancer. Exp Ther Med. 2019;18(6):4490–8.
Junhua Lv YZ, Suwei Gao C, Zhang Y, Chen W, Li Y-G, Yang Q, Zhou. Feng Liu endothelial-specific m6A modulates mouse hematopoietic stem and progenitor cell development via notch signaling. Cell Res. 2018;28(2):249–52.
Li Z, Weng H, Su R, Weng X, Zuo Z, Li C, et al. FTO Plays an Oncogenic Role in Acute Myeloid Leukemia as a N6-Methyladenosine RNA Demethylase. Cancer Cell. 2017;31(1):127 – 41.
Xiulin Jiang BL, Zhi Nie L, Duan QX, Jin Z, Yang C. Yongbin Chen. The role of m6A modification in the biological functions and diseases. Signal Transduction and Targeted Therapy 2021;6(74).
Cheng Y, Gao Z, Zhang T, Wang Y, Xie X, Han G et al. Decoding m6A RNA methylome identifies PRMT6-regulated lipid transport promoting AML stem cell maintenance. Cell Stem Cell. 2022;30.
Xuan Zhou ZX, Hehua M et al. A novel cuproptosis-related gene signature can predict prognosis in acute myeloid leukemia. Res Square. 2022.
Yang T, Xu R, Su Q, Wang H, Liu F, Dai B et al. Chelerythrine hydrochloride inhibits proliferation and induces mitochondrial apoptosis in cervical cancer cells via PI3K/BAD signaling pathway. Toxicol In Vitro. 2020;68(104965).
Kang K, Jiang H, Zhang S, Cheng B. Antitumor Effects of Chelerythrine: A literature review. Nat Prod Commun. 2022;17(6):1934578X221103028.
Acknowledgements
We acknowledge the TCGA, TARGET, GTEx, GSE65409, Beat AML and DTP NCI-60 databases for providing their platforms and contributors for uploading their meaningful datasets. We thank all the authors and funding for supporting this research.
Funding
This work was supported by grants from Science and technology innovation team of the Suzhou vocational health college (Grant No. SZWZYTD202201), Key Technology Program of Suzhou People’s Livelihood Technology Projects (Grant No. SKY2021029), Medical and Health Projects in Zhejiang Province (2021KY400), The Open Project of Jiangsu Biobank of Clinical Resources (TC2021B009), Application Study Project of Public Welfare of Zhejiang Province (GF22H085719), and Qing-Lan Project of Jiangsu Province in China (2021, 2022). Zhejiang Provincial Natural Science Foundation (LGF21H250002), Chinese Traditional Medicine Science and Technology Projects of Zhejiang Province (2021ZB002, 2022ZB002).
Author information
Authors and Affiliations
Contributions
X.B., J.C., Y.Z., M.Y., J.D., Z.T., X.X. , G.M., Z.W., J.C. performed the research, analysed the data, and reviewed the manuscript; SB.L., T.X., J.H. designed the research, analysed and interpreted the data, and wrote the manuscript. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
All patient and healthy control person care and experimental procedures were approved by the Research Ethics Committee of the First Affiliated Hospital of Soochow University. We confirmed that all methods were performed in accordance with the relevant guidelines and regulations. All patients signed the informed consent form.
Consent for publication
The data used in this study came from the TCGA, TARGET, GTEx, GSE65409, Beat AML and DTP NCI-60 databases, and this study was in full compliance with the published guidelines of public databases. Accordingly, no additional informed consent was provided.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Bao, X., Chi, J., Zhu, Y. et al. High FAAP24 expression reveals poor prognosis and an immunosuppressive microenvironment shaping in AML. Cancer Cell Int 23, 117 (2023). https://doi.org/10.1186/s12935-023-02937-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12935-023-02937-3