Comprehensive analysis of peroxiredoxins expression profiles and prognostic values in breast cancer

Abstract

Background

The peroxiredoxins (PRDXs) gene family has been demonstrated to participate in carcinogenesis and development of numerous cancers and the prognostic values in several cancers have been evaluated already. Purpose of our research is to explore the expression profiles and prognostic values of PRDXs in breast cancer (BrCa).

Methods

The transcriptional levels of PDRX family members in primary BrCa tissues and their association with intrinsic subclasses were analyzed using UALCAN database. Then, the genetic alterations of PDRXs were examined by cBioPortal database. Moreover, the prognostic values of PRDXs in BrCa patients were investigated via the Kaplan-Meier plotter.

Results

The transcriptional levels of most PRDXs family members in BrCa tissues were significantly elevated compared with normal breast tissues. Meanwhile, dysregulated PRDXs expression was associated with intrinsic subclasses of BrCa. Besides, copy number alterations (CNA) of PRDXs positively regulated their mRNA expressions. Furthermore, high mRNA expression of PRDX4/6 was significantly associated with poor overall survival (OS) in BrCa patients, while high mRNA expression of PRDX3 was notably related to favorable OS. Simultaneously, high mRNA expression of PRDX1/2/4/5/6 was significantly associated with shorter relapse-free survival (RFS) in BrCa patients, while high mRNA expression of PRDX3 was notably related to favorable RFS. In addition, the prognostic value of PRDXs in the different clinicopathological features based on intrinsic subclasses and chemotherapeutic treatment of BrCa patients was further assessed in the KM plotter database.

Conclusion

Our findings systematically elucidate the expression profiles and distinct prognostic values of PRDXs in BrCa, which might provide novel therapeutic targets and potential prognostic biomarkers for BrCa patients.

Background

Breast cancer (BrCa) has the highest morbidity among all female’s cancers worldwide, which may cause 41,760 cancer-related deaths in the United States (US) in 2019 according to the prediction by the American Cancer Society (ACS) [1]. Being a multifaceted disease, BrCa can be classified into various subclasses based on the expression status of estrogen receptor (ER), progesterone receptor (PR), epidermal growth factor receptor 2 (HER2) and antigen ki-67 (Ki-67), which suggest different therapeutic guidance and prognostic implications for BrCa patients [2]. Although the constant amelioration of comprehensive therapies for BrCa has significantly decreased the mortality of BrCa in recent years, it is still necessary to further explore the potential mechanism of oncogenesis and progression of BrCa.

Peroxiredoxins (PRDXs), a family of antioxidant enzymes in eukaryotes, containing six isoforms (PRDX1, PRDX2, PRDX3, PRDX4, PRDX5, and PRDX6), which catalyze the reduction reaction of peroxide and maintain the balance of intracellular hydrogen peroxide (H₂O₂) levels [3]. As momentous regulators in diverse signaling pathways, PRDXs are of great significance to the signal transduction and cells metabolism [4]. Expression of PRDXs will be significantly upregulated when cells are under oxidative stress conditions. Several researches have suggested that overexpression of PRDXs may play dichotomous role in oncogenesis of tumors, where they could either stimulate the progression of cancers or suppress the development of cancers [5]. An increasing number of studies have observed the preliminary functions and ambiguous prognostic values of PRDXs in cancerous diseases [6,7,8,9]. However, the expression profiles and prognostic values of PRDXs in BrCa samples are still elusive up to now. Our research aims to explore the differential expression and potential roles of PRDXs in BrCa.

The Cancer Genome Atlas (TCGA) program, which launched by the US National Cancer Institute (NCI) and the National Human Genome Research Institute (NHGRI), attempts to sequence the entire genome of more than 10,000 tumor samples and to distinguish the genetic changes specific for each cancer [10,11,12]. Along with the successful implementation of the TCGA project, massive genomic information is accumulating exponentially. Over the past few years, many interactive and user-friendly online platforms based on the TCGA database greatly elevate the efficiency of TCGA database analysis and increasing amounts of tumor biomarkers have been identified on the strength of these websites [13, 14].

Therefore, in the current research, we first compared the transcriptional levels of PRDXs in BrCa and adjacent breast tissues using UALCAN database. In addition, the cBioPortal database was used to analyze the genetic alterations of PDRXs and the correlation with transcriptional levels. Moreover, the Kaplan-Meier plotter database was used to assess the prognostic effects of PRDXs mRNA expression in patients with BrCa. Overall, our research preliminarily but systematically characterizes the expression profiles of PRDXs in BrCa and reveals that the detection of the PRDXs expression status of BrCa patients may be valuable and potential biomarkers for prognostic assessment.

Materials and methods

Gene expression analysis via UALCAN

UALCAN (http://ualcan.path.uab.edu/) is an online open-access platform based on level 3 RNA-seq and clinical information from TCGA database [15]. It can be used to analyze relative transcriptional levels of potential genes of interest between cancerous and paired normal tissues and association of the transcriptional levels with clinicopathologic features. In the current study, UALCAN was applied to analyze the transcriptional levels of PDRXs family members in primary BrCa tissues and their association with intrinsic subclasses. All the BrCa cases available on UALCAN were included in our research.

Data-mining analysis based on cBioPortal

cBioPortal (www.cbioportal.org/) is a user-friendly, interactive website resource and provides visualization, analysis, and download of large-scale cancer genomics datasets [16, 17]. In the current study, we analyzed the genetic alterations of PDRXs family members, which contained mutations and putative copy-number alterations (CNA) from GISTIC. Furthermore, we download the data of putative copy-number alterations and mRNA expression z-Scores to evaluate the association between various CNAs and transcriptional levels of PRDXs. Tumor samples with RNA-seq and CAN data on cBioPortal were included in our research which contains total 1076 BrCa samples.

Survival analysis by Kaplan-Meier plotter

Kaplan-Meier Plotter (KM Plotter, http://kmplot.com/analysis/) is an online database containing gene expression profiles and survival information of cancer patients [18]. The prognostic values of PRDXs (PRDX1, PRDX2, PRDX3, PRDX4, PRDX5, and PRDX6) at mRNA level in BrCa was analyzed using all BrCa samples available on KM Plotter. The patients’ cohorts were split at the median expression of each PRDXs mRNA level. The subgroup analysis of the prognostic value of PRDXs in BrCa patients was further performed according to intrinsic subclasses and different regimens of chemotherapy. All cohorts were compared with Kaplan-Meier survival plots. Hazard ratio (HR), 95% confidence interval (95% CI), and log-rank P value were calculated and displayed online.

Statistical analysis

All statistical analyses were performed on the bioinformatics database online or using SPSS 25.0 software (Chicago, IL). The differential mRNA expression of PRDXs in BrCa tissues was analyzed by Student’s t-test. Kaplan-Meier survival plots were generated with survival curves compared by log-rank test. For all analyses, Differences were considered statistically significant if P values were less than 0.05.

Results

Transcriptional levels of PRDXs in BrCa samples

In order to evaluate the exact expression profiles of PRDXs members in BrCa patients, the differential transcriptional levels of PRDX family members between BrCa and paired normal breast tissue was evaluated by UALCAN database. As shown in Fig. 1, the transcriptional level of PRDX1 (Fig. 1a, P < 0.001), PRDX2 (Fig. 1b, P < 0.001), PRDX4 (Fig. 1d, P < 0.001), and PRDX5 (Fig. 1e, P < 0.001) was significantly upregulated in BrCa tissues compared with paracancerous tissues. However, the transcriptional level of PRDX3 (Fig. 1c, P = 0.251) showed a non-significant difference in BrCa tissues compared with paracancerous tissues. Besides, the transcriptional level of PRDX6 was significantly downregulated in BrCa tissues compared with paracancerous tissues (Fig. 1f, P < 0.001).

Fig. 1

Transcriptional levels of PRDXs in paracancerous and BrCa tissues. Comparison of PRDX1, PRDX2, PRDX3, PRDX4, PRDX5, PRDX6 mRNA expression in paracancerous (n = 114) and BrCa (n = 1097) tissues in TVGA dataset based on data mining via UALCAN. a, b, d, e) The transcriptional level of PRDX1, PRDX2, PRDX4, and PRDX5 was significantly upregulated in BrCa tissues compared with paracancerous tissues. c The transcriptional level of PRDX3 showed a non-significant difference in BrCa tissues compared with paracancerous tissues. f The transcriptional level of PRDX6 was significantly downregulated in BrCa tissues compared with paracancerous tissues

Transcriptional levels of PRDXs in different BrCa subclasses

Classification of intrinsic subclasses is helpful in the prediction of therapeutic response and prognosis of BrCa [19]. So, we next compared the differential transcriptional levels of PRDX family members according to different intrinsic subclasses of BrCa. As shown in Fig. 2, mRNA expressions of PRDXs family members were significantly correlated with intrinsic subclasses of BrCa. Patients who were with HER2-positive and triple-negative subclasses BrCa tended to express higher PRDXs (exclude PRDX2, PRDX3) mRNA, while express lower PRDX2 and PRDX3 mRNA. The highest mRNA expressions of PRDX1/5/6 were found in HER2-positive tissues (Fig. 2a, e, f), and the highest mRNA expressions of PRDX4 were found in triple-negative tissues (Fig. 2d). Besides, the lowest mRNA expressions of PRDX2/3 were found in triple-negative tissues (Fig. 2b, c). Taken together, these findings above revealed that transcriptional levels of PRDXs family members were significantly correlated with intrinsic subclasses in BrCa patients.

Fig. 2

Transcriptional levels of PRDXs in various BrCa subclasses. The transcriptional level of PRDXs in BrCa patients with different subclasses, PRDXs mRNA was significantly downregulated (PRDX3) or upregulated (other PRDXs) in HER2-positive and triple-negative BrCa tissues compared with luminal BrCa tissues. a PRDX1. b PRDX2. c PRDX3. d PRDX4. e PRDX5. f PRDX6

Genetic alterations of PRDXs in BrCa samples

DNA copy number alterations (CNA) are most common genetic alterations which participates in oncogenesis of cancers via regulating cancer-related gene expression [20,21,22]. In the fact that most PRDX family members was dysregulated in BrCa tissues, we speculated that DNA CNA may regulate the transcriptional levels of PRDXs. Next, we analyzed genetic alteration in PRDXs and correlations with their mRNA expressions based on cBioPortal website. As shown in Fig. 3a and Table 1, low amplification rate of PRDXs was found in BrCa patients. However, although copy gain (gain and amplification) of PRDXs was not frequent, it was still associated with notably upregulated PRDXs mRNA levels compared with the copy-neutral (diploid) and copy-loss (shallow deletion and deep deletion) cases (Fig. 3b-g). To conclude, the results suggested that PRDX mRNA expressions were regulated by their DNA copy number alterations.

Fig. 3

Correlation between the genetic alterations of PRDXs and mRNA levels in BrCa tissues. a Oncoprint in cBioPortal database exhibited the proportion and distribution of specimens with genetic alterations in PRDXs. The Figure was cropped on the right to exclude specimens without any alterations. b-g Copy gain (gain and amplification) of PRDXs was associated with notably upregulated PRDXs mRNA levels compared with the copy-neutral (diploid) and copy-loss (shallow deletion and deep deletion) cases. b PRDX1. c PRDX2. d PRDX3. e PRDX4. f PRDX5. g PRDX6

Table 1 Frequency and proportion of genetic alterations of PRDXs in BrCa

Prognostic values of PRDXs mRNA expression in all BrCa samples

Further, we employed KM plotter to evaluate the prognostic values of PRDX family members. As shown in Fig. 4 and Fig. 5, high mRNA expression of PRDX4 (HR = 1.54, 95% CI: 1.24–1.91, P < 0.001), and PRDX6 (HR = 1.26, 95% CI: 1.02–1.56, P = 0.033) were significantly associated with poor overall survival (OS) of BrCa patients, while high mRNA expression of PRDX3 was notably related to favorable OS of BrCa patients (HR = 0.74, 95% CI: 0.59–0.91, P = 0.005). However, other PRDXs mRNA expression showed a null association with prognosis of BrCa patients.

Fig. 4

Prognostic value of PRDXs mRNA in all BrCa patients (OS). OS curves were plotted to evaluate the prognostic value of PRDXs mRNA expression. a PRDX1 (208680_at). b PRDX2 (39729_at). c PRDX3 (201619_at). d PRDX4 (201923_at). e PRDX5 (1560587_s_at). f PRDX6 (200845_s_at)

Fig. 5

Prognostic value of PRDXs mRNA in BrCa patients (RFS). RFS curves were plotted to evaluate the prognostic value of PRDXs mRNA expression. a PRDX1 (208680_at). b PRDX2 (39729_at). c PRDX3 (201619_at). d PRDX4 (201923_at). e PRDX5 (1560587_s_at). f PRDX6 (200845_s_at)

We next analyzed the associations between PRDXs mRNA expression and RFS of BrCa patients and found that high mRNA expression of PRDX1 (HR = 1.31, 95% CI: 1.17–1.46, P < 0.001), PRDX2 (HR = 1.30, 95% CI: 1.17–1.45, P < 0.001), PRDX4 (HR = 1.54, 95% CI: 1.38–1.72, P < 0.001), PRDX5 (HR = 1.21, 95% CI: 1.04–1.42, P = 0.015) and PRDX6 (HR = 1.20, 95% CI: 1.08–1.34, P < 0.001) were significantly associated with shorter relapse-free survival (RFS) of BrCa patients, while high mRNA expression of PRDX3 was notably related to favorable RFS of BrCa patients (HR = 0.81, 95% CI: 0.72–0.90, P < 0.001). Overall, the findings above implied that mRNA expressions of PRDX3/4/6 were remarkably correlated with BrCa patients’ both OS and RFS, which might be identified as promising biomarkers to predict the survival of BrCa patients.

Prognostic values of PRDXs mRNA in different BrCa subclasses

To further analyze the association of PRDXs mRNA expression with various BrCa subclasses, we detected the survival effects of PRDX family members in 4 subclasses of BrCa patients, including basal-like, luminal A, luminal B, and HER2 positive BrCa. As shown in Table 2, for PRDX1 (HR = 1.52, 95% CI: 1.06–2.17, P = 0.022), PRDX4 (HR = 1.58, 95% CI: 1.10–2.27, P = 0.012), and PRDX6 (HR = 1.49, 95% CI: 1.04–2.13, P = 0.027), high mRNA expression was associated with unfavorable OS in luminal A BrCa patients, respectively. Besides, high mRNA expression of PRDX4 (HR = 1.55, 95% CI: 1.06–2.27, P = 0.023) predicted poor OS in luminal B BrCa patients as well.

Table 2 Association between prognostic value of PRDXs mRNA expression and different subclasses in BrCa

Moreover, we also analyzed the associations between PRDXs mRNA expression and RFS of various subclasses BrCa patients and the results indicated that high expressions of PRDX1 (luminal A: HR = 1.23, 95% CI: 1.06–1.49, P = 0.009), PRDX2 (luminal A: HR = 1.23, 95% CI: 1.04–1.46, P = 0.016; luminal B: HR = 1.67, 95% CI: 1.37–2.03, P < 0.001), PRDX4 (luminal A: HR = 1.56, 95% CI: 1.31–1.85, P < 0.001; luminal B: HR = 1.56, 95% CI: 1.28–1.89, P < 0.001), PRDX6 (luminal A: HR = 1.29, 95% CI: 1.08–1.53, P = 0.004; luminal B: HR = 1.27, 95% CI: 1.05–1.54, P = 0.015) were correlated with worse RFS in luminal BrCa patients. In addition, low expression of PRDX3 (HR = 0.77, 95% CI: 0.60–1.00, P = 0.046) and high expression of PRDX5 (HR = 1.71, 95% CI: 1.23–2.37, P = 0.001) predicted unfavorable RFS in basal-like BrCa patients. Thus, these results suggested the roles of PRDXs as potential prognostic predictors in BrCa patients with different subclasses.

Prognostic values of PRDXs mRNA in BrCa patients with diverse regimens of chemotherapy

Next, we also checked the prognostic effects of PRDX family members in BrCa patients with different chemotherapies, including adjuvant chemotherapy, neoadjuvant chemotherapy and non-chemotherapy. As shown in Table 3, high expression of PRDX1, PRDX2, PRDX3, and PRDX4 were significantly correlated with poor OS in BrCa patients with adjuvant chemotherapy. In addition, in BrCa patients who didn’t receive any chemotherapies, high expression PRDX1 (HR = 1.22, 95% CI: 1.03–1.44, P = 0.019), PRDX4 (HR = 1.37, 95% CI: 1.16–1.62, P < 0.001), PRDX6 (HR = 1.19, 95% CI: 1.00–1.40, P = 0.044) and low expression of PRDX3 (HR = 0.84, 95% CI: 0.71–1.00, P = 0.043) were associated with worse RFS. However, contrary to the prognostic effect of PRDX5 in total BrCa patients, high expression of PRDX5 (HR = 0.50, 95% CI: 0.30–0.81, P = 0.005) predicted better RFS in BrCa patients with adjuvant chemotherapy. Thus, these results suggested the roles of PRDXs as potential prognostic predictors in BrCa patients with different regimens of chemotherapy.

Table 3 Association between prognostic value of PRDXs mRNA expression and various chemotherapies in BrCa

Discussion

Reactive oxygen species (ROS), including the superoxide radical, the hydroxyl radical, H₂O₂ and etc., are the most important type of free radicals which produces secondary toxic metabolic products, such as peroxynitrites and nitrogen oxides, posing a lethal threat to cells by damaging DNA [23, 24]. Peroxiredoxins, one of the most significant antioxidant enzyme systems that include SOD, CAT, and GPx, were significnatly upregulated under oxidative stress conditions and mainly participate in the defense against oxidative [25, 26]. Several studies have observed that imbalances between the generation of ROS and PRDXs in tumor cells could lead to oxidative stress and the induction of cell apoptosis [27].

It has been demonstrated that PRDXs expression was significantly dysregulated during carcinogenesis of cancers and played dichotomous roles in oncogenesis. Overexpression of PRDX1 in BrCa has been observed to be positively associated with tumor grade and acted as dominant role in management of exogeneous oxidative stress [28, 29]. PRDX2, has been reported to specifically regulate the oxidative and metabolic stress response of metastatic breast cancer cells in lungs. Besides, overexpressed PRDX2 participates in chemo-resistant in BrCa cells [30, 31]. The function of PRDX3 and PRDX4 in BrCa is largely ambiguous, Liu et al. reveals that downregulation of PRDX3 potentiates PP2-induced apoptosis in MCF-7 cells, which suggests the tumor suppressor role of PRDX3 [32]. PRDX4 has been demonstrated to mediate osteoclast activation by human BrCa cells and enhance the aggressive phenotype [33]. Last but not least, there is no available research about the exact function of PRDX5 and PRDX6 in BrCa, but overexpression of PRDX5 and prognostic values have been observed in numerous cancers, including ovarian cancer and endometrial cancer [34, 35]. As well, the tumor promoter role of PRDX6 in cancers has also been suggested in colorectal cancer, lung cancer and so on [36, 37]. Although several studies that investigate the role of PRDXs in BrCa have been published, little is known about individual PRDXs expression and their effects on survival of BrCa patients.

In the present study, we first investigated the differential transcriptional levels of PRDX family members between BrCa and adjacent tissues and the results showed that the transcriptional levels of PRDX1, PRDX2, PRDX4, and PRDX5 were significantly upregulated in BrCa tissues. Besides, the transcriptional level of PRDX6 was significantly downregulated in BrCa tissues. However, the transcriptional level of PRDX3 showed a non-significant difference in BrCa tissues compared with adjacent tissues. we also compared the differential transcriptional levels of PRDX family members according to different intrinsic subclasses of BrCa and found mRNA expressions of PRDXs family members were significantly correlated with intrinsic subclasses of BrCa.

DNA copy number alterations (CNA) are most common genetic alterations which affect carcinogenesis and development of cancers by regulating cancer-related gene expression [20,21,22]. When we used cBioPortal to inspect genetic alteration in PRDXs and correlations with their mRNA expressions, the results showed that copy gain (gain and amplification) of most PRDXs was not frequent in BrCa, but it was still associated with notable upregulated PRDXs mRNA levels. Amplification is a positive factor to upregulate gene expression [38, 39]. However, to our largely surprise, the copy-gain frequency of PRDX6 accounts for a large proportion in total BrCa samples, but a significant decrease expression of PRDX6 was exhibited in BrCa tissues. We speculated that amplification of PRDX6 gene may upregulate PRDX6 expression in BrCa tissues compared with paired breast tissues. However, limited amounts of normal tissues expression data were included in TCGA dataset, thus, the opposite phenomenon that PRDX6 amplification decreased transcriptional level was exhibited.

The Kaplan-Meier plotter is an online database which is available to assess the prognostic effect of genes expression on survival in designate cancers. The primary purpose of the tool is a meta-analysis-based biomarker assessment and a lot of prognostic biomarkers have been identified based on this platform [40,41,42,43]. It has been reported that overexpression PRDX6 participates in cisplatin resistance in ovarian cancer and predicts poor OS and PFS [8, 44]. Several researches also observed the promising prognostic values of PRDXs in lung cancer and endometrial cancer [7, 35, 45]. However, the prognostic values of PRDXs in BrCa patients are largely unknown. Here, we found that high mRNA expression of PRDX4/6 were significantly associated with poor OS of BrCa patients and high mRNA expression of PRDX1/2/4/5/6 were significantly associated with shorter RFS of BrCa patients, while high mRNA expression of PRDX3 was notably related to favorable OS and DFS, which suggests the tumor suppressor role of PRDX3 in BrCa. Besides, the prognostic values of PRDXs mRNA in different BrCa subclasses and in BrCa patients with diverse regimens of chemotherapy were also assessed and results suggested the potential roles of PRDXs in predicting prognosis of BrCa patients with various subclasses and different regimens of chemotherapy.

Although this study systematically demonstrates the prognostic value of PRDXs in breast cancer, this research has several limitations as well. The major limitation is that online database only provides the expression of PRDXs mRNA level, which may not fully represent the expression of PRDXs at the protein level. In further study, western blotting, immumohistochemical staining and other protein detection techniques will be applied to determinate the protein level of PRDXs in breast cancer. Furthermore, the possible mechanisms that PRDXs is involved in the tumorigenesis and progression of breast cancer need to be further studied. Besides, although the large sample analyses based on meta-analysis have some advantages, but some essential information form one single center may be missing, such as some therapeutic information.

Conclusion

In summary, we systemically analyzed the expression profiles and prognostic values of PRDXs in BrCa. Our results revealed that PRDX1/2/4/5/6 might be the potential therapeutic targets for BrCa therapy, whereas PRDX3/4/6 were promising prognostic biomarkers for predicting OS and RFS of BrCa patients. Overall, our research provided a systematic insight into the heterogeneous and complex roles of PRDXs in the carcinogenesis of BrCa.