Identification of Fanconi anemia pathway genes as novel prognostic biomarkers and therapeutic targets for breast cancer
Highlight box
Key findings
• This study suggests that Fanconi anemia (FA) pathway genes have potential as prognostic biomarkers and therapeutic targets in breast cancer treatment.
What is known and what is new?
• FA pathway genes play a role in DNA repair, and mutations in these genes are linked to certain cancers, including breast cancer.
• This manuscript provides evidence that the expression and mutations of FA pathway genes could serve as indicators for breast cancer prognosis. Additionally, the study identifies these genes as possible targets for therapeutic intervention, suggesting that their modulation could have clinical benefits.
What is the implication, and what should change now?
• The findings imply that FA pathway genes could be utilized to improve prognostic accuracy in breast cancer patients, aiding in personalized treatment strategies. They may also offer a new avenue for targeted therapy.
Introduction
Breast cancer is the most common type of cancer among women. It accounts for 30% of all new cancer diagnoses, and is the second leading cause of cancer associated mortalities among women (1). It is comprised of a heterogeneous group of diseases with different histopathological characteristics and high genetic variability, and is therefore, characterized by different prognostic outcomes. Specific breast cancer subtypes are defined by their histopathological appearance and expression of hormone receptors, including oestrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) (2). Among the genetic risk factors, pathogenic mutations in high- and moderate-risk cancer predisposition genes such as BRCA1 and BRCA2, have an important impact on breast cancer development (3). Due to tumor heterogeneity, the current breast cancer biomarkers for predicting prognosis have some limitations, therefore, there is a need to establish new biomarkers as prognostic indicators to effectively enhance prognosis and individualize breast cancer treatment.
Fanconi anemia (FA) is a rare, cancer prone disease with mutations in at least 22 genes (4). In addition, protein products of these 22 FA genes along with the FA-associated proteins, interact in a common cellular pathway, known as the FA pathway, to repair DNA interstrand cross-links (ICLs). The FA pathway plays a major role in responses to replication stress by facilitating the resolution of DNA lesions arising from DNA replication (5). Moreover, amplification and gain-of-function mutations in FA genes are advantageous in cancer cells by alleviating replication stress and mitigating chemotherapeutics induced DNA damage (6). Studies have documented the key functions of FA genes in different kinds of cancers, including prostate cancer (7), colorectal cancer (8), hepatocellular carcinoma (9), bladder cancer (10), and breast cancer (11). Breast cancer susceptibility genes, BRCA1 and BRCA2, also known as FANCS and FANCD1, respectively, are involved in the FA pathway. Furthermore, some of the FA pathway genes are associated with clinicopathological features in breast cancer and could serve as cancer diagnostic or prognostic biomarkers (12). However, studies regarding the expression patterns and prognostic values of all the FA pathway genes are few. Therefore, on the basis of the analysis of thousands of gene expressions or variations in copy numbers published online, this study explored expressions and mutations in FA pathway genes in patients with breast cancer to determine their expression patterns, distinct prognostic values, and potential functions of these genes in breast cancer. We present this article in accordance with the REMARK reporting checklist (available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-772/rc).
Methods
The Oncomine Cancer Microarray (ONCOMINE) database
ONCOMINE (www.oncomine.org) is a cancer microarray database that allows genome-wide expression analysis (13). In this study, the FA pathway gene transcriptional levels in different cancers were analyzed using ONCOMINE in this study. Datasets were screened with thresholds of P value (1E−4), fold change (2) and gene rank (top 10%). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).
University of ALabama at Birmingham Cancer (UALCAN) database
UALCAN (http://ualcan.path.uab.edu) is a web-portal used for performing in-depth analyses on The Cancer Genome Atlas (TCGA) transcriptomic data (14). It was used to assess the messenger RNA (mRNA) expression levels of FA pathway genes in breast cancer tissues and in the corresponding normal tissues.
Kaplan-Meier plotter
The prognostic values of FA pathway genes in breast cancer were evaluated using the Kaplan-Meier Plotter (www.kmplot.com). This online platform can be used to assess the significance of the expression levels of various genes on clinical outcomes in cancer patients (15). Furthermore, the platform was used to analyze the associations between the expression levels of FA pathway genes and clinicopathological features in breast cancer.
cBio Cancer Genomics Portal (cBioPortal) database
The cBioPortal for Cancer Genomics (http://cbioportal.org) is a comprehensive web resource that can visualize and analyze multidimensional cancer genomic data (16). Therefore, data from this platform was used to analyze changes in the frequency of FA pathway genes in breast cancer.
Gene Expression Profiling Interactive Analysis (GEPIA) database
GEPIA (http://gepia.cancer-pku.cn/) is an interactive web application for gene expression analysis. It is based on 9,736 tumor and 8,587 normal tissue samples from TCGA and the Genotype-Tissue Expression (GTEx) databases (17). Thus, it was used to assess the correlations between the expression levels of FA pathway genes in breast cancer.
Gene Multi-Association Network Integration Algorithm (GeneMANIA) database
GeneMANIA (http://www.genemania.org) is a web interface for the identification of related genes from many large, publicly available biological datasets (18). In this study, the relationship between FA pathway genes and their interactive genes was analyzed using this database.
The Database for Annotation, Visualization and Integrated Discovery (DAVID)
DAVID (http://david.abcc.ncifcrf.gov) is an annotation, visualization and integrated discovery database that is able to extract biological features associated with specific genes (19). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of FA pathway genes and their closely associated genes in breast cancer were performed by DAVID and visualized using R language.
Tumor Immune Estimation Resource (TIMER) database
TIMER (cistrome.shinyapps.io/timer) is a web interface for comprehensive molecular characterization of tumor-immune interactions (20). The expression level of FA pathway genes in breast cancer and their correlation with tumor purity and infiltrating immune cells such as B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells were assessed using TIMER.
Statistical analysis
The statistical analysis was performed using R (version 4.2.3) and SPSS (version 25.0). Descriptive statistics for continuous data were presented as mean ± standard deviation (SD). For comparing two groups or multiple groups, independent t-tests or one-way analysis of variance (ANOVA) were applied for normally distributed data. For data that did not meet the assumption of normality, the Mann-Whitney U test or Kruskal-Wallis test was used. Survival analysis was performed using the Kaplan-Meier method, and differences between survival curves were evaluated using the log-rank test. Correlation analysis was conducted to assess the strength and direction of relationships between variables using Pearson correlation coefficients, depending on the data distribution. Hazard ratios and their 95% confidence intervals were calculated through Cox proportional hazards regression analysis to identify independent risk factors. All statistical tests were two-tailed, with a significance level set at P<0.05.
Results
mRNA expression levels of FA pathway genes in breast cancer
Based on the ONCOMINE data, the transcriptional levels of various FA pathway genes were found to be elevated in various types of cancers, such as colorectal cancer, cervical cancer and breast cancer (Figure 1). We further compared the mRNA expression levels of FA pathway genes in breast cancer and their corresponding normal tissues. Results from the UALCAN database revealed that the genes were significantly upregulated in breast cancer compared to their corresponding normal tissues, except for FANCE and FANCM which were downregulated in tumors (Figure 2).


The prognostic value of FA pathway genes in breast cancer
The Kaplan-Meier plotter was used to determine the potential prognostic value of FA pathway genes in breast cancer (Figure 3). About half of the genes that were highly expressed were shown to exhibit a significant positive correlation with worse recurrence-free survival (RFS) of breast cancer patients. These genes include: FANCB, FANCG, FANCL, UBE2T (FANCT), FANCI, BRIP1 (FANCJ), BRCC5 (FANCR), BRCA1 (FANCS), MAD2L2 (FANCV) and RFWD3 (FANCW). In contrast, elevated mRNA expression levels of FANCC, SLX4 (FANCP), PALB2 (FANCN), XRCC2 (FANCU) and ERCC4 (FANCQ) were associated with favorable RFS. The other genes did not exhibit a significant correlation with the RFS of breast cancer patients.

We further used the Cox proportional hazards regression analysis to investigate the correlations between the expression levels of FA pathway genes and patients’ overall survival (OS), distant metastasis-free survival (DMFS) and post progression survival (PPS). The result forest plots are presented in Figure 4.

Association between the prognostic value of FA pathway genes in breast cancer with different clinicopathological features
The relationship between FA pathway genes and different clinicopathological features was investigated to elucidate on the roles of these genes in breast cancer prognosis. The clinicopathological features included cancer grade, ER status, PR status, and HER2 status (Table S1). It was revealed that elevated mRNA expression levels of FANCI were associated with poor RFS in grade 1 breast cancer. Additionally, elevated mRNA expression levels of FANCA, FANCE, FANCL, FANCI, BRCC5 and RFWD3 were correlated to worse RFS in grade 2 breast cancer. However, FANCF and PALB2 were associated with better RFS in grade 2 breast cancer. Moreover, low mRNA expression levels of FANCI and BRCA1 were associated with better RFS in grade 3 breast cancer, while BRIP1 and RAD51C were identified as good prognostic factors for grade 3 breast cancer. These results are presented in Table 1.
Table 1
Gene | Grade | Cases | HR (95% CI) | P value |
---|---|---|---|---|
FANCA | 1 | 345 | 1.62 (0.95–2.77) | 0.07 |
2 | 901 | 1.41 (1.11–1.79) | 0.005** | |
3 | 903 | 1.04 (0.84–1.29) | 0.72 | |
FANCB | 1 | 108 | 1.13 (0.39–3.22) | 0.82 |
2 | 227 | 1.17 (0.71–1.95) | 0.54 | |
3 | 444 | 0.9 (0.66–1.22) | 0.49 | |
FANCC | 1 | 345 | 1.16 (0.69–1.95) | 0.58 |
2 | 901 | 0.93 (0.73–1.19) | 0.56 | |
3 | 903 | 0.94 (0.75–1.17) | 0.57 | |
FANCE | 1 | 345 | 1.32 (0.78–2.23) | 0.30 |
2 | 901 | 1.3 (1.02–1.65) | 0.03* | |
3 | 903 | 0.98 (0.79–1.22) | 0.84 | |
FANCF | 1 | 345 | 0.73 (0.43–1.22) | 0.23 |
2 | 901 | 0.75 (0.59–0.96) | 0.02* | |
3 | 903 | 0.93 (0.75–1.16) | 0.52 | |
FANCG | 1 | 345 | 1.66 (0.97–2.82) | 0.059 |
2 | 901 | 1.19 (0.93–1.51) | 0.16 | |
3 | 903 | 1.19 (0.96–1.48) | 0.11 | |
FANCL | 1 | 345 | 1.2 (0.71–2.03) | 0.49 |
2 | 901 | 1.35 (1.06–1.72) | 0.01* | |
3 | 903 | 1.16 (0.93–1.44) | 0.19 | |
FANCM | 1 | 108 | 1.34 (0.46–3.87) | 0.59 |
2 | 227 | 1.49 (0.9–2.49) | 0.12 | |
3 | 444 | 1.1 (0.81–1.5) | 0.54 | |
UBE2T (FANCT) | 1 | 108 | 1.33 (0.46–3.85) | 0.59 |
2 | 227 | 1.23 (0.74–2.04) | 0.43 | |
3 | 444 | 0.83 (0.61–1.14) | 0.24 | |
SLX4 (FANCP) | 1 | 108 | 1.29 (0.45–3.71) | 0.64 |
2 | 227 | 1.19 (0.71–1.99) | 0.51 | |
3 | 444 | 1.01 (0.74–1.37) | 0.97 | |
FANCD2 | 1 | 108 | 1.09 (0.38–3.11) | 0.87 |
2 | 227 | 0.72 (0.43–1.2) | 0.21 | |
3 | 444 | 0.89 (0.65–1.22) | 0.47 | |
FANCI | 1 | 345 | 1.81 (1.06–3.1) | 0.02* |
2 | 901 | 1.72 (1.35–2.2) | 1e−05*** | |
3 | 903 | 1.31 (1.05–1.62) | 0.01* | |
FANCD1 | 1 | 345 | 1.15 (0.68–1.93) | 0.61 |
2 | 901 | 1.22 (0.96–1.55) | 0.11 | |
3 | 903 | 0.96 (0.77–1.2) | 0.74 | |
BRIP1 (FANCJ) | 1 | 345 | 0.94 (0.33–2.69) | 0.91 |
2 | 901 | 1.16 (0.7–1.92) | 0.57 | |
3 | 903 | 0.72 (0.53–0.99) | 0.041* | |
PALB2 (FANCN) | 1 | 345 | 1.21 (0.72–2.03) | 0.47 |
2 | 901 | 0.78 (0.61–0.99) | 0.041* | |
3 | 903 | 1.13 (0.91–1.4) | 0.28 | |
RAD51C (FANCO) | 1 | 345 | 0.95 (0.57–1.59) | 0.84 |
2 | 901 | 1.04 (0.82–1.32) | 0.74 | |
3 | 903 | 0.8 (0.64–0.99) | 0.042* | |
BRCC5 (FANCR) | 1 | 345 | 1.1 (0.65–1.84) | 0.72 |
2 | 901 | 1.36 (1.07–1.73) | 0.01* | |
3 | 903 | 1.1 (0.89–1.37) | 0.37 | |
BRCA1 (FANCS) | 1 | 345 | 1.46 (0.86–2.47) | 0.16 |
2 | 901 | 1.21 (0.95–1.54) | 0.11 | |
3 | 903 | 1.5 (1.2–1.87) | 3e−04*** | |
XRCC2 (FANCU) | 1 | 345 | 1.5 (0.88–2.56) | 0.13 |
2 | 901 | 1.1 (0.87–1.4) | 0.41 | |
3 | 903 | 1.1 (0.89–1.37) | 0.38 | |
MAD2L2 (FANCV) | 1 | 108 | 0.78 (0.27–2.26) | 0.65 |
2 | 227 | 1.47 (0.88–2.45) | 0.14 | |
3 | 444 | 0.85 (0.62–1.16) | 0.30 | |
RFWD3 (FANCW) | 1 | 345 | 1.67 (0.97–2.85) | 0.06 |
2 | 901 | 1.5 (1.18–1.91) | 1e−04*** | |
3 | 903 | 1.23 (0.99–1.53) | 0.06 | |
ERCC4 (FANCQ) | 1 | 345 | 0.66 (0.23–1.9) | 0.44 |
2 | 901 | 0.85 (0.51–1.42) | 0.54 | |
3 | 903 | 1.05 (0.77–1.43) | 0.77 |
*, P<0.05; **, P<0.01; ***, P<0.001. FA, Fanconi anemia; HR, hazard ratio; CI, confidence interval.
Moreover, BRCA1 was found to be a promising marker for unfavorable prognosis in both ER positive and negative patients, while FANCA, FANCB, FANCE, FANCG, FANCL, UBE2T, FANCI, BRIP1, RAD51C, BRCC5 and RFWD3 were significantly associated with unfavorable RFS in ER positive patients. However, FANCC and FANCD1 were associated with favorable RFS in ER negative patients (Table 2, Table S2).
Table 2
Gene | ER status | Cases | HR (95% CI) | P value |
---|---|---|---|---|
FANCA | Positive | 2,061 | 1.41 (1.19–1.66) | 4.3e–05*** |
Negative | 801 | 0.95 (0.76–1.19) | 0.67 | |
FANCB | Positive | 762 | 1.61 (1.2–2.15) | 0.001** |
Negative | 347 | 0.82 (0.59–1.14) | 0.23 | |
FANCC | Positive | 2,061 | 1.08 (0.92–1.28) | 0.33 |
Negative | 801 | 0.75 (0.6–0.94) | 0.01* | |
FANCE | Positive | 2,061 | 1.19 (1.01–1.4) | 0.04* |
Negative | 801 | 1.02 (0.81–1.27) | 0.88 | |
FANCG | Positive | 2,061 | 1.39 (1.18–1.64) | 8.5e−05*** |
Negative | 801 | 1.02 (0.81–1.28) | 0.86 | |
FANCL | Positive | 2,061 | 1.18 (1–1.39) | 0.049* |
Negative | 801 | 1.16 (0.93–1.46) | 0.19 | |
UBE2T (FANCT) | Positive | 762 | 1.89 (1.4–2.55) | 2.4e−05*** |
Negative | 347 | 0.89 (0.64–1.24) | 0.50 | |
FANCI | Positive | 2061 | 1.86 (1.57–2.2) | 1.5e−13*** |
Negative | 801 | 1.12 (0.89–1.4) | 0.32 | |
FANCD1 | Positive | 2,061 | 1.16 (0.98–1.36) | 0.08 |
Negative | 801 | 0.75 (0.6–0.95) | 0.01* | |
BRIP1 (FANCJ) | Positive | 2,061 | 1.67 (1.24–2.24) | 6.3e−04*** |
Negative | 801 | 0.91 (0.65–1.26) | 0.56 | |
RAD51C(FANCO) | Positive | 2,061 | 1.32 (1.12–1.56) | 8.3e−04*** |
Negative | 801 | 0.85 (0.68–1.07) | 0.17 | |
BRCC5 (FANCR) | Positive | 2,061 | 1.33 (1.13–1.57) | 6.6e−04*** |
Negative | 801 | 1.06 (0.85–1.33) | 0.62 | |
BRCA1 (FANCS) | Positive | 2,061 | 1.6 (1.36–1.89) | 1.9e−08*** |
Negative | 801 | 1.53 (1.22–1.92) | 2.4e−04*** | |
RFWD3 (FANCW) | Positive | 2,061 | 1.32 (1.12–1.56) | 8.2e−04*** |
Negative | 801 | 1.18 (0.94–1.47) | 0.16 |
*, P<0.05; **, P<0.01; ***, P<0.001. FA, Fanconi anemia; ER, estrogen receptor; HR, hazard ratio; CI, confidence interval.
Furthermore, elevated mRNA expression levels of FANCA, FANCG, UBE2T, FANCI, BRIP1, RAD51C, BRCC5 and BRCA1 were shown to contribute to unfavorable RFS in PR positive patients while elevated expression levels of SLX4 and ERCC4 were associated with shorter RFS in PR negative patients (Table 3, Table S2).
Table 3
Gene | PR status | Cases | HR (95% CI) | P value |
---|---|---|---|---|
FANCA | Positive | 589 | 1.5 (1.06–2.13) | 0.02* |
Negative | 549 | 1.08 (0.81–1.45) | 0.59 | |
FANCG | Positive | 589 | 1.78 (1.24–2.54) | 0.001** |
Negative | 549 | 1.19 (0.89–1.59) | 0.24 | |
UBE2T (FANCT) | Positive | 489 | 1.71 (1.16–2.52) | 0.006** |
Negative | 372 | 1.05 (0.74–1.51) | 0.77 | |
SLX4 (FANCP) | Positive | 489 | 1.13 (0.77–1.65) | 0.52 |
Negative | 372 | 1.77 (1.23–2.55) | 0.002** | |
FANCI | Positive | 589 | 1.83 (1.28–2.61) | 7.4e−04*** |
Negative | 549 | 1.13 (0.85–1.51) | 0.41 | |
BRIP1 (FANCJ) | Positive | 589 | 1.97 (1.33–2.93) | 5.5e−04*** |
Negative | 549 | 0.94 (0.66–1.35) | 0.75 | |
RAD51C (FANCO) | Positive | 589 | 1.49 (1.05–2.11) | 0.02* |
Negative | 549 | 1.02 (0.76–1.37) | 0.88 | |
BRCC5 (FANCR) | Positive | 589 | 1.44 (1.01–2.04) | 0.041* |
Negative | 549 | 1.32 (0.98–1.76) | 0.06 | |
BRCA1 (FANCS) | Positive | 589 | 1.7 (1.19–2.43) | 0.003** |
Negative | 549 | 1.27 (0.95–1.7) | 0.11 | |
ERCC4 (FANCQ) | Positive | 589 | 1.02 (0.7–1.49) | 0.90 |
Negative | 549 | 1.58 (1.1–2.27) | 0.01* |
*, P<0.05; **, P<0.01; ***, P<0.001. FA, Fanconi anemia; PR, progesterone receptor; HR, hazard ratio; CI, confidence interval.
PALB2 exhibited a significant association with unfavorable RFS in HER2 positive patients. However, FANCA, FANCB, FANCC, FANCF, FANCG, UBE2T, SLX4, FANCI, BRIP1, RAD51C, BRCA1 and ERCC4 were associated with poor RFS while FANCF was associated with better RFS in HER2 negative patients. Nevertheless, FANCD2 was a good prognostic factor in both HER2 positive and negative patients (Table 4, Table S2).
Table 4
Gene | HER2 status | Cases | HR (95% CI) | P value |
---|---|---|---|---|
FANCA | Positive | 252 | 1.02 (0.66–1.57) | 0.95 |
Negative | 800 | 1.35 (1.04–1.75) | 0.02* | |
FANCB | Positive | 150 | 1.01 (0.59–1.73) | 0.97 |
Negative | 635 | 1.41 (1.05–1.9) | 0.02* | |
FANCC | Positive | 252 | 0.83 (0.54–1.29) | 0.41 |
Negative | 800 | 1.37 (1.05–1.78) | 0.02* | |
FANCF | Positive | 252 | 0.89 (0.58–1.38) | 0.60 |
Negative | 800 | 0.76 (0.58–0.99) | 0.04* | |
FANCG | Positive | 252 | 1.42 (0.91–2.19) | 0.12 |
Negative | 800 | 1.32 (1.01–1.71) | 0.04* | |
UBE2T (FANCT) | Positive | 150 | 0.61 (0.35–1.06) | 0.08 |
Negative | 635 | 2.28 (1.67–3.12) | 9.7e−08*** | |
SLX4 (FANCP) | Positive | 150 | 1.45 (0.84–2.51) | 0.18 |
Negative | 635 | 1.45 (1.07–1.96) | 0.02* | |
FANCD2 | Positive | 150 | 0.51 (0.29–0.9) | 0.02* |
Negative | 635 | 1.39 (1.03–1.87) | 0.03* | |
FANCI | Positive | 252 | 1.11 (0.72–1.71) | 0.64 |
Negative | 800 | 1.79 (1.37–2.34) | 1.5e−05*** | |
BRIP1 (FANCJ) | Positive | 252 | 0.61 (0.35–1.05) | 0.07 |
Negative | 800 | 1.98 (1.46–2.69) | 8.4e−06*** | |
PALB2 (FANCN) | Positive | 252 | 1.63 (1.05–2.53) | 0.03* |
Negative | 800 | 1.2 (0.92–1.56) | 0.17 | |
RAD51C (FANCO) | Positive | 252 | 0.74 (0.48–1.14) | 0.17 |
Negative | 800 | 1.44 (1.1–1.87) | 0.007** | |
BRCA1 (FANCS) | Positive | 252 | 1.13 (0.73–1.75) | 0.58 |
Negative | 800 | 1.73 (1.32–2.26) | 5.9e−05*** | |
ERCC4 (FANCQ) | Positive | 252 | 0.81 (0.47–1.4) | 0.46 |
Negative | 800 | 1.36 (1.01–1.83) | 0.044* |
*, P<0.05; **, P<0.01; ***, P<0.001. FA, Fanconi anemia; HER2, human epidermal growth receptor 2; HR, hazard ratio; CI, confidence interval.
Based on the TP53 status, elevated expression levels of UBE2T were strongly associated with worse RFS while FANCL, FAND2 and RAD51C were correlated with better RFS in TP53-mutated patients. In patients with wild type TP53, FANCA, FANCI, BRCA1 and RFWD3 were associated with unfavorable RFS (Table 5, Table S2).
Table 5
Gene | TP53 status | Cases | HR (95% CI) | P value |
---|---|---|---|---|
FANCA | Mutated | 188 | 0.87 (0.54–1.41) | 0.58 |
Wild type | 273 | 1.84 (1.19–2.83) | 0.005** | |
FANCL | Mutated | 188 | 0.62 (0.38–1) | 0.049* |
Wild type | 273 | 1.09 (0.72–1.66) | 0.68 | |
UBE2T (FANCT) | Mutated | 132 | 2.04 (1.12–3.73) | 0.02* |
Wild type | 82 | 1.45 (0.62–3.4) | 0.39 | |
FANCD2 | Mutated | 132 | 0.54 (0.3–0.99) | 0.043* |
Wild type | 82 | 1.58 (0.67–3.7) | 0.29 | |
FANCI | Mutated | 188 | 0.88 (0.55–1.41) | 0.60 |
Wild type | 273 | 1.75 (1.14–2.69) | 0.01* | |
RAD51C (FANCO) | Mutated | 188 | 0.58 (0.36–0.95) | 0.03* |
Wild type | 273 | 1.1 (0.72–1.68) | 0.64 | |
BRCA1 (FANCS) | Mutated | 188 | 1.09 (0.68–1.75) | 0.71 |
Wild type | 273 | 1.69 (1.1–2.59) | 0.02* | |
RFWD3 (FANCW) | Mutated | 188 | 0.98 (0.61–1.57) | 0.94 |
Wild type | 273 | 1.72 (1.12–2.65) | 0.01* |
*, P<0.05; **, P<0.01. FA, Fanconi anemia; HR, hazard ratio; CI, confidence interval.
Genetic alterations and interaction analysis of FA pathway genes in breast cancer
We further performed a comprehensive analysis of the molecular characteristics of genes in the FA pathway. The frequency of genetic alterations in these genes among breast cancer patients was determined using the cBioPortal database. It was found that mRNA deregulation was one of the most important single factors for genetic alterations in different kinds of breast cancers (Figure 5A). Mutation and amplification were the most common alterations in these samples. In addition, OncoPrint was used to show a visual summary of alterations in the FA pathway genes across a set of breast cancer samples (Figure 5B). Moreover, expression correlations were determined using GEPIA to further define the relationships among the FA pathway genes. There was a low to high positive correlation among most FA pathway genes (Figure S1A). Additionally, a network of FA pathway genes and their functionally related genes was constructed using GeneMANIA (Figure S1B). Twenty genes were found to be closely associated with the regulatory functions of differentially expressed FA pathway genes. These genes were FAAP24, FAAP100, RAD51B, XRCC3, RAD52, BARD1, BLM, RMI1, TOP3A, ATR, FAN1, APITD1, TOPBP1, USP1, G2E3, RAD51D, RPA1, RPA2, ERCC1 and ATRIP.

Functional enrichment analysis of FA-related genes in breast cancer
We used UALCAN to isolate the top 50 genes that were positively and negatively correlated with individual genes of the FA pathway in breast cancer. This was done to explore the underlying mechanisms of FA pathway genes in cancer. In addition, DAVID was used to perform GO and KEGG pathway enrichment analysis of the FA-associated genes in breast cancer. Figure S2 shows the most highly enriched GO items that were positively and negatively correlated with FA pathway genes in breast cancer. Among the positively correlated GO items, the most enriched biological process (BP) term was cell division, the most enriched cellular component (CC) term was nucleoplasm, and the most enriched molecular function (MF) term was protein binding (Figure S2A). In the negatively correlated GO items, the FA-associated genes were shown to participate in various functions, especially signal recognition particle (SRP)-dependent co-translational protein targeting to membrane as well as ribosome and structural constituent of ribosome (Figure S2B). Moreover, KEGG pathway enrichment analysis revealed that the positively correlated pathways were mainly involved in cell cycle, oocyte meiosis and viral carcinogenesis (Figure S2C). Besides, ribosomes, along with their metabolic pathways, were the most enriched pathways that were found to be negatively correlated with FA pathway genes in breast cancer (Figure S2D).
Correlation between immune infiltration and FA pathway genes in patients with breast cancer
Given that inflammatory responses and infiltrating immune cells can affect breast cancer prognosis, we evaluated the association between differentially expressed FA pathway genes and immune cell infiltration using the TIMER database. It was found that the mRNA expression levels of FA pathway genes were positively associated with tumor purity, while FANCE, FANCM, BRCA2 (FANCD1) and MAD2L2 (FANCV) had no significant correlation with the tumor purity of patients with breast cancer (Figure 6). Furthermore, most FA pathway genes were positively associated with immune infiltration levels of B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils and dendritic cells, except for FANCE of macrophages, FANCF of dendritic cells, FANCG of macrophages, FANCL of macrophages, UBE2T of macrophages, RAD51C of CD4+ T cells, MAD2L2 of CD8+ T cells and macrophages, which showed significant negative associations.

Discussion
Chemotherapy is one of the most important treatments for breast cancer after surgery. Approximately one-third of patients with breast cancer present metastases, which are the main cause of death in these patients (21). It has been documented that tumor responses to chemotherapeutic drugs are closely associated with the regulation of the DNA repair system (22). Some tumor cells can resist DNA damage drugs by activating self-DNA repair mechanisms (23). Moreover, deficiency in the proteins involved in DNA damage repair is considered a major determinant of the responses to chemotherapy in cancer cells (24). Previous study reported that the FA pathway, also referred to as the FA-BRCA pathway, can modulate tumor progression and immunotherapeutic effects (6). However, the prognostic values and biological functions of FA pathway genes in breast cancer have not been well elucidated.
DNA repair involves multiple enzymes and genes. Inactivating mutations in DNA repair components are common and often lead to certain DNA repair deficiencies. Therefore, cancer cells become hyperdependent on the remaining repair pathways for survival and proliferation (25). The FA pathway is a stepwise multiprotein complex pathway that confers cellular hypersensitivity to DNA intercalating substances, such as cisplatin, that trigger DNA ICLs (26). FA pathway activation status may serve as a clinical biomarker for cancer patients at different treatment stages. Herein, the expression levels of FA pathway genes in breast cancer were determined before evaluating their association with survival outcomes in breast cancer patients. Expression levels of 20 genes were shown to be significantly higher in breast cancer tissues than in the corresponding normal tissues, except for FANCE and FANCM, which were downregulated in tumors. Moreover, elevated expression levels of FANCB, FANCG, FANCL, UBE2T, FANCI, BRIP1, BRCC5, BRCA1, MAD2L2 and RFWD3 in breast cancer were associated with worse RFS. However, elevated mRNA levels of FANCC, SLX4, PALB2, XRCC2 and ERCC4 were correlated with a favorable RFS. Furthermore, elevated expression levels of FANCA, FANCG, UBE2T, FANCI, FANCD1, BRCC5, BRCA1 and RFWD3 were associated with worse OS. These findings imply that FANCG, UBE2T, FANCI, BRCC5, BRCA1 and RFWD3 exhibited better prospects for utilization as prognostic biomarkers in breast cancer patients. It has been documented that the prognosis of breast cancer patients is associated with tumor pathological tissue type, such as ER, PR and HER2 status, which have played a role in the identification of which patients are likely to benefit from endocrine therapy or targeted therapy (27). TP53 is the most frequent mutational target in human cancers. Mutations in TP53 are associated with different types of malignancies and adverse prognoses, including during breast cancer (28). In this study, most FA pathway genes showed a close relationship with worse RFS of breast cancer patients with different clinicopathological features including cancer grade, ER status, PR status, HER2 status and TP53 status. Collectively, the FA pathway genes are potential therapeutic targets and prognostic biomarkers for breast cancer.
Next-generation sequencing has uncovered the frequency of mutations and copy number alterations across different cancer types and demonstrated that alterations in DNA repair mechanisms are common events in carcinogenesis. Mutations with high variant allele frequencies (VAFs) indicate early appearance of tumorigenesis or tremendous contribution to the later expansion of tumor cells (29). Moreover, compensatory mutations in BRCA1 and BRCA2 that restore homologous recombination (HR) functionality in initially cisplatin sensitive tumors are able to develop cisplatin resistance (30). Comprehensively revealing mutation characteristics in breast cancer elucidates on the mutational diversity among different molecular subtypes, enables the identification of potential treatment biomarkers, and provides a basis for genomic targeting strategies and clinical trials (31). Given the significant differential expression of genes in the FA pathway, we further explored their molecular characteristics. It was found that mutation, fusion, amplification, deep deletion and multiple alterations were the main mutational signatures of FA pathway genes in breast cancer. Notably, amplification was the main characteristic of gene mutations in the FA pathway genes, which meant that the FA pathway could be a significant compensatory DNA repair pathway for cancer cells. More importantly, these mutational signatures may be new therapeutic targets for precision medicine, providing opportunities for personalized treatment strategies based on the imperfection of patient’s DNA repair networks.
Cell response to DNA damage is a complex mechanism involving multiple protein networks with interconnected functions that are responsible for damage detection, cell cycle regulation and DNA repair. Establishing the underlying mechanisms involved in the association between FA pathway and breast cancer, besides DNA damage repair, will have significant implications in clinical practice (32). In this study, a low to high expression correlation among FA pathway genes in breast cancer was obtained, suggesting that they played a synergistic role in tumorigenesis and cancer progression. Then, we determined the core genes that were potentially associated with FA pathway gene functions. Some of them were identified as important gene regulators. For instance, it has been shown that FANCM and its binding partner, FAAP24, suppress the formation of DNA double-stranded breaks and mitotic recombination in a manner that is dependent on FANCM translocase activity (33). Moreover, BRCA1-BARD1 are required for fork protection and are associated with cancer development (34). Functional enrichment analysis was then performed to elucidate on the biological functions of FA pathway genes in breast cancer. The FA-related genes were found to be primarily positively associated with the cell cycle and nucleoplasm. However, the FA pathway genes were also negatively correlated with genes involved in SRP-dependent cotranslational protein targeting to membrane and ribosome. The ribosome plays a critical role in normal cellular physiology, in cellular responses to internal and external environmental stimuli, and in the pathogenesis of human diseases (35). Under stress situations, a decreased ribosomal activity and reduced protein synthesis are shown, subsequently leading to nuclear mobilization and DNA repair activation to minimize the negative impact to cell growth (36). This may be the mechanism through which ribosome-related genes are down-regulated. Recent finding has begun to elucidate the interplay between ribosomal biogenesis, which means ribosomal synthesis and DNA repair (37). FANCI is required for ribosomal biogenesis, and may function by coordinating rDNA replication and transcription (38). Nonetheless, the exact functions of FA pathway genes in breast cancer should be investigated further.
The cancer immune microenvironment plays an important role in tumor progression (39). In recent years, immunotherapy has been found to be a promising therapy for cancer, and the development of immunological biomarkers has been of increasing importance (40). In this study, there was a significant positive correlation between the mRNA expressions of FA pathway genes and tumor purity in breast cancer. Tumor purity is highly associated with genomic patterns and immune phenotypes, which is substantially inversely correlated with tumor heterogeneity (41). Targeting DNA repair processes may influence the adaptive immune system by leading to an increased number of mutations, and subsequently increased burden of neoantigens, which in turn increases tumor heterogeneity, resulting in a higher probability of recognition by the immune system, and this has the potential to be exploited in therapeutic approaches (42). Pan-cancer analysis suggested that increasing mutation load is linearly correlated with increasing immune activity in the tumor microenvironment of a tumor and is likely to influence immune recognition (43). Therefore, FA pathway genes were found to be potential therapeutic targets for breast cancer and could be combined with immunotherapy. In addition, the mRNA expression levels of FA pathway genes were also highly correlated with various immune cell infiltrations. These findings imply that FA pathway genes are not only prognostic indicators but also reflect the “immune-hot” status in breast cancer. DNA repair can also influence how the innate immune system initially responds to a tumor and recruits the adaptive immune system to the malignancy site (44). Alterations in DNA repair can influence how the adaptive, innate, or both parts of the immune system respond to the underlying malignancy. However, further studies are required to verify the potential role of FA pathway genes in breast cancer as predictive biomarkers of immunotherapeutic responses.
There are some limitations in this study. Analysis on the transcriptional level can reflect some immune status aspects, but not wholesome changes. Independent cohort and in vitro or in vivo studies should be performed to validate our results.
In summary, there is a significant correlation between the mRNA expression levels of FA pathway genes and tumor prognosis as well as the cancer immune infiltration. This implies that the genes may mediate tumor progression and exert immunotherapeutic effects in breast cancer. Therefore, elucidation of how these genes is regulated during tumor progression may highlight their potential prognostic and therapeutic role in breast cancer.
Conclusions
This study elucidates on the expression, mutations and prognostic values of FA pathway genes in breast cancer. Analysis of the relationship between FA pathway gene expression and clinicopathological characteristics in breast cancer indicated that FA pathway genes could be promising prognostic biomarkers in patients with breast cancer, and may be novel targets for breast cancer therapy. More studies are needed to explore the exact mechanisms and therapeutic roles of FA pathway genes in breast cancer. It is possible that the FA pathway genes will be effective prognostic markers of breast cancer in future.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the REMARK reporting checklist. Available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-772/rc
Peer Review File: Available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-772/prf
Funding: This research was funded by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tcr.amegroups.com/article/view/10.21037/tcr-24-772/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013).
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Bray F, Laversanne M, Sung H, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2024;74:229-63. [Crossref] [PubMed]
- Li Z, Wei H, Li S, et al. The Role of Progesterone Receptors in Breast Cancer. Drug Des Devel Ther 2022;16:305-14. [Crossref] [PubMed]
- Marcotte R, Sayad A, Brown KR, et al. Functional Genomic Landscape of Human Breast Cancer Drivers, Vulnerabilities, and Resistance. Cell 2016;164:293-309. [Crossref] [PubMed]
- Degan P, Cappelli E, Regis S, et al. New Insights and Perspectives in Fanconi Anemia Research. Trends Mol Med 2019;25:167-70. [Crossref] [PubMed]
- Webster ALH, Sanders MA, Patel K, et al. Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer. Nature 2022;612:495-502. [Crossref] [PubMed]
- Nalepa G, Clapp DW. Fanconi anaemia and cancer: an intricate relationship. Nat Rev Cancer 2018;18:168-85. [Crossref] [PubMed]
- Kumar A, Coleman I, Morrissey C, et al. Substantial interindividual and limited intraindividual genomic diversity among tumors from men with metastatic prostate cancer. Nat Med 2016;22:369-78. [Crossref] [PubMed]
- Parsa FG, Nobili S, Karimpour M, et al. Fanconi Anemia Pathway in Colorectal Cancer: A Novel Opportunity for Diagnosis, Prognosis and Therapy. J Pers Med 2022;12:396. [Crossref] [PubMed]
- Takahashi J, Masuda T, Kitagawa A, et al. Fanconi Anemia Complementation Group E, a DNA Repair-Related Gene, Is a Potential Marker of Poor Prognosis in Hepatocellular Carcinoma. Oncology 2022;100:101-13. [Crossref] [PubMed]
- Roh YG, Mun JY, Kim SK, et al. Fanconi Anemia Pathway Activation by FOXM1 Is Critical to Bladder Cancer Recurrence and Anticancer Drug Resistance. Cancers (Basel) 2020;12:1417. [Crossref] [PubMed]
- Wang M, Pan B, Hu Y, et al. STIL facilitates the development and malignant progression of triple-negative breast cancer through activation of Fanconi anemia pathway via interacting with KLF16. Transl Oncol 2024;46:102010. [Crossref] [PubMed]
- Feng L, Jin F. Expression and prognostic significance of Fanconi anemia group D2 protein and breast cancer type 1 susceptibility protein in familial and sporadic breast cancer. Oncol Lett 2019;17:3687-700. [Crossref] [PubMed]
- Rhodes DR, Yu J, Shanker K, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia 2004;6:1-6. [Crossref] [PubMed]
- Chandrashekar DS, Karthikeyan SK, Korla PK, et al. UALCAN: An update to the integrated cancer data analysis platform. Neoplasia 2022;25:18-27. [Crossref] [PubMed]
- Győrffy B. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput Struct Biotechnol J 2021;19:4101-9. [Crossref] [PubMed]
- Cerami E, Gao J, Dogrusoz U, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov 2012;2:401-4. [Crossref] [PubMed]
- Li C, Tang Z, Zhang W, et al. GEPIA2021: integrating multiple deconvolution-based analysis into GEPIA. Nucleic Acids Res 2021;49:W242-6. [Crossref] [PubMed]
- Warde-Farley D, Donaldson SL, Comes O, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res 2010;38:W214-20. [Crossref] [PubMed]
- Sherman BT, Hao M, Qiu J, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res 2022;50:W216-21. [Crossref] [PubMed]
- Li T, Fu J, Zeng Z, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res 2020;48:W509-14. [Crossref] [PubMed]
- Ibragimova MK, Tsyganov MM, Kravtsova EA, et al. Organ-Specificity of Breast Cancer Metastasis. Int J Mol Sci 2023;24:15625. [Crossref] [PubMed]
- Lu X, Liu R, Wang M, et al. MicroRNA-140 impedes DNA repair by targeting FEN1 and enhances chemotherapeutic response in breast cancer. Oncogene 2020;39:234-47. [Crossref] [PubMed]
- Hopkins JL, Lan L, Zou L. DNA repair defects in cancer and therapeutic opportunities. Genes Dev 2022;36:278-93. [Crossref] [PubMed]
- Pitroda SP, Bao R, Andrade J, et al. Low Recombination Proficiency Score (RPS) Predicts Heightened Sensitivity to DNA-Damaging Chemotherapy in Breast Cancer. Clin Cancer Res 2017;23:4493-500. [Crossref] [PubMed]
- Klinakis A, Karagiannis D, Rampias T. Targeting DNA repair in cancer: current state and novel approaches. Cell Mol Life Sci 2020;77:677-703. [Crossref] [PubMed]
- Dauden MI, López-Perrote A, Llorca O. RUVBL1-RUVBL2 AAA-ATPase: a versatile scaffold for multiple complexes and functions. Curr Opin Struct Biol 2021;67:78-85. [Crossref] [PubMed]
- Zhang X. Molecular Classification of Breast Cancer: Relevance and Challenges. Arch Pathol Lab Med 2023;147:46-51. [Crossref] [PubMed]
- Marvalim C, Datta A, Lee SC. Role of p53 in breast cancer progression: An insight into p53 targeted therapy. Theranostics 2023;13:1421-42. [Crossref] [PubMed]
- Flynn J, Samant N, Schneider-Nachum G, et al. Mutational fitness landscape and drug resistance. Curr Opin Struct Biol 2023;78:102525. [Crossref] [PubMed]
- Tung NM, Garber JE. BRCA1/2 testing: therapeutic implications for breast cancer management. Br J Cancer 2018;119:141-52. [Crossref] [PubMed]
- Lang GT, Jiang YZ, Shi JX, et al. Characterization of the genomic landscape and actionable mutations in Chinese breast cancers by clinical sequencing. Nat Commun 2020;11:5679. [Crossref] [PubMed]
- Fang CB, Wu HT, Zhang ML, et al. Fanconi Anemia Pathway: Mechanisms of Breast Cancer Predisposition Development and Potential Therapeutic Targets. Front Cell Dev Biol 2020;8:160. [Crossref] [PubMed]
- Wang H, Li S, Oaks J, et al. The concerted roles of FANCM and Rad52 in the protection of common fragile sites. Nat Commun 2018;9:2791. [Crossref] [PubMed]
- Daza-Martin M, Starowicz K, Jamshad M, et al. Isomerization of BRCA1-BARD1 promotes replication fork protection. Nature 2019;571:521-7. [Crossref] [PubMed]
- Kang J, Brajanovski N, Chan KT, et al. Ribosomal proteins and human diseases: molecular mechanisms and targeted therapy. Signal Transduct Target Ther 2021;6:323. [Crossref] [PubMed]
- Korsholm LM, Gál Z, Lin L, et al. Double-strand breaks in ribosomal RNA genes activate a distinct signaling and chromatin response to facilitate nucleolar restructuring and repair. Nucleic Acids Res 2019;47:8019-35. [Crossref] [PubMed]
- Sakthivel D, Brown-Suedel A, Bouchier-Hayes L. The role of the nucleolus in regulating the cell cycle and the DNA damage response. Adv Protein Chem Struct Biol 2023;135:203-41. [Crossref] [PubMed]
- Sondalle SB, Longerich S, Ogawa LM, et al. Fanconi anemia protein FANCI functions in ribosome biogenesis. Proc Natl Acad Sci U S A 2019;116:2561-70. [Crossref] [PubMed]
- Bianchini G, De Angelis C, Licata L, et al. Treatment landscape of triple-negative breast cancer - expanded options, evolving needs. Nat Rev Clin Oncol 2022;19:91-113. [Crossref] [PubMed]
- Kumar N, Gann PH, McGregor SM, et al. Quantification of subtype purity in Luminal A breast cancer predicts clinical characteristics and survival. Breast Cancer Res Treat 2023;200:225-35. [Crossref] [PubMed]
- Revkov E, Kulshrestha T, Sung KW, et al. PUREE: accurate pan-cancer tumor purity estimation from gene expression data. Commun Biol 2023;6:394. [Crossref] [PubMed]
- Amodio V, Lamba S, Chilà R, et al. Genetic and pharmacological modulation of DNA mismatch repair heterogeneous tumors promotes immune surveillance. Cancer Cell 2023;41:196-209.e5. [Crossref] [PubMed]
- Qing T, Jun T, Lindblad KE, et al. Diverse immune response of DNA damage repair-deficient tumors. Cell Rep Med 2021;2:100276. [Crossref] [PubMed]
- Ding JH, Xiao Y, Yang F, et al. Guanosine diphosphate-mannose suppresses homologous recombination repair and potentiates antitumor immunity in triple-negative breast cancer. Sci Transl Med 2024;16:eadg7740. [Crossref] [PubMed]