Skip to main content

Classification of PR-positive and PR-negative subtypes in ER-positive and HER2-negative breast cancers based on pathway scores



PR loss in ER+/HER2- breast cancer indicates worse prognosis and insensitivity to anti-estrogen therapy, while the mechanisms of PR loss in ER+/HER2- breast cancer remain unrevealed.


In this study, ER+/PR+/HER2- and ER+/PR-/HER2- breast cancer cases from TCGA were used. 1387 pathways were analyzed and used as variables for classifying the two groups with LASSO regression.


ER+/PR+/HER2- and ER+/PR-/HER2- breast cancer groups can be classified by a combination of 13 pathways using their activity score. Among the 13 pathways, those involving growth factors and ion-channel transporters were most significant in the distinction, followed by pathways involving immune modulation and cell metabolism. Two growth factor pathways, EGF and IGF-1, were deferentially regulated in ER+/PR+/HER2- and ER+/PR-/HER2- groups.


In conclusion, this study indicated in ER+/HER2- breast cancers the various status of PR expression can be an indication of molecular variation, particularly for the growth factor pathway activation.

Peer Review reports


About 80 % of breast cancers were hormonal receptor-positive which means they express at least one of the two hormonal receptors, estrogen receptor (ER) and progesterone receptor (PR) [1, 2]. Both ER and PR are ligand-activated transcription factors that promote the expression of specific gene sets by binding to their promoters [3]. Although the expression of ER and PR were often closely correlated and highly consistent, there is still discordance in some breast cancers. It was reported that 15 % of ER-positive breast cancer were PR negative while in PR-positive breast cancer, only 2 % were ER-negative [4], which suggests that ER expressed more widely than PR. Indeed, about 12 % of all breast cancer patients have the hormonal receptor status as ER+/PR- [4]. One of the possible mechanisms for PR loss could be the copy number loss of the PGR gene which encodes for PR.

The expression of PR was mostly controlled by activated ER [5], while it can also be regulated by growth factor pathways [6] and cyclin D1 [7]. Analysis of both the Surveillance, Epidemiology and End Results (SEER) [4] and the National Cancer Database (NCD) [8] datasets has confirmed that ER+/PR- breast cancer has worse survival than ER+/PR + breast cancers. Moreover, it was shown that the ER+/PR- group was more resistant to selective estrogen receptor modulator (SERM) therapies [9, 10]. The exact mechanism of the association between PR loss and worse prognosis was remained to be elucidated although several clues existed [11]. One possible mechanism was that the relative overexpression of HER2 in the ER+/PR- group compared with the ER+/PR + group made ER+/PR- breast cancer resistant to tamoxifen [12, 13]. However, in ER+/HER2- breast cancer, the effect of PR expression on prognosis still exists, indicating the presence of other mechanisms [14]. PR has long been considered as the downstream gene of ER since its expression was shown to be induced by estradiol [15], while recent researches proved that PR has strong negative regulatory effect on ER activity [16]. The complex correlation between ER and PR makes it essential to elucidate the molecular characteristiscs and clinical significance of PR loss in ER+/HER2- breast cancer [16].

In this study, we analyzed and compared the pathway activities in ER+/PR+/HER2- and ER+/PR-/HER2- breast cancer using transcriptomic and gene amplification data from TCGA. The two groups were abbreviated as ER+/PR + and ER+/PR-, thus all the studied patients were HER2- unless otherwise specified.


Clinicopathological characteristics and survival analysis of the ER+/PR + and ER+/PR- breast cancer patients

In the TCGA cohort, 592 patients have information on ER and PR expression based on immunohistochemical staining. The expression of ER and PR were reported as the percentage of cells with positive expression. ER + and PR + were defined as expressions of more than 1 %, according to the ASCO/CAP guideline [17].

Cases with positive expression of ER and PR were further subgrouped into 10 categories with a 10 % interval. By this definition, 123 patients belong to the ER-/PR-, while 60 cases were ER+/PR- and 409 cases were ER+/PR+ (Fig. 1).

Fig. 1
figure 1

Cell percentage categories of ER and PR expression in HER2- breast cancers from TCGA. The case number of each specific ER and PR expression status was summarized and presented. ER and PR expression were reported either as negative or as the percentage of cells showed positive expression which was subgrouped into 10 categories with a 10 % interval

Clinicopathological characteristics of the ER+/PR- and ER+/PR + groups were compared (Table 1). In concordance with previous studies, the two groups showed a significant difference in clinicopathological characteristics including clinical stage and nodal status. Specifically, ER+/PR- group patients have both more advanced clinical stage and nodal status. In terms of PAM50 intrinsic subtypes distribution of the two groups, it was found that the PR- group was more enriched in luminal B and basal-like type than the PR + group. In addition, comprehensive molecular portraits of the two groups were analyzed including DNA methylation, copy number variation, and miRNA profile which have all been tested and clustered in the previous study [18]. Distributions in methylation cluster and copy number cluster were found to be different between the two groups while no difference was detected in the miRNA cluster distribution. For the DNA methylation cluster, the PR- group was enriched in cluster 5 which was a cluster that overlapped with basal-like mRNA subtype while the PR + group was more enriched in cluster 2 which seems to be a mixture of mRNA subtypes. For copy number cluster, PR- group was enriched in cluster 2 with PR + enriched in cluster 1. Copy number clusters 1 and 2 were previously found to be correlated with luminal A and basal-like subtypes respectively. All the above results indicated that molecularly PR- group was more similar with basal-like subtype and has a more advanced clinical stage than the PR + group.

Table 1 Patient Clinicopathological features of the two selected groups and the multivariate and univariate analysis in the TCGA dataset

The survival analyses were performed by comparing the overall survival of the ER+/PR-/HER2- and ER+/PR+/HER2- group (Fig. 2). In accordance with previous studies, the overall survival of ER+/PR-/HER2- was significantly poorer than the ER+/PR+/HER2- group. To further validate the result, the same analysis was performed in the Surveillance, Epidemiology, and End Results (SEER) database (Fig. 2). A total of 110,930 ER+/HER2- breast cancer patients were included in the SEER dataset analysis with 97,397 of them being PR + and 13,533 of them being PR-. The clinicopathological characteristics of the included population were presented in Table S1. Li et al. have previously studied the ER+/PR + and ER+/PR- group in the SEER dataset while not considering the status of the HER2 receptor [4].

Fig. 2
figure 2

Survival analysis of the ER+/PR + and ER+/PR- groups. Kaplan-Meier survival curves were plotted and the log-rank test was performed to compare the overall survival of ER+/PR + and ER+/PR- groups in TCGA a and SEER b database

Both multivariate and univariate analyses were conducted to determine the prognostic value of different clinicopathological features including PR status in all the selected ER+/HER2- patients (Table 2). In univariate analysis, it was found that age, disease stage, tumor size and nodal status were significant prognostic factors for overall survival while factors including age, disease stage, tumor size and copy number cluster were statistically significant in multivariate analysis.

Table 2 Univariate and multivariate analysis of the prognostic value of clinicopathological features in ER+/HER2- patients using Cox proportional hazard modeling in the TCGA dataset

Classification of ER+/PR + and ER+/PR- breast cancer with pathway activities using LASSO methodology

A strategy named least absolute shrinkage and selection operator (LASSO) was used to select the most accurate and compact set of pathways that can differentiate the ER+/PR + and ER+/PR- group. Logistic-LASSO is a regression method that minimizes the usual sum of squared errors which penalizes the regression coefficients. From the 1387 variables/pathways pool, LASSO was able to reduce the number of pathways needed down to thirteen in the final model (Fig. 3). The final model with all the selected pathways and statistics is shown in Table 3. In the final model, the ER+/PR+/HER2- and ER+/PR-/HER2- groups can be distinguished by pathway activities of the 13 pathways, achieving a 0.8625 Area Under the receiver operation characteristic Curve (AUC) (Fig. 3). The most distinct pathways between PR + and PR- groups were growth factors and ion-channel transporter. Also, the FOXA1 pathway which represents the luminal features, and the NOTCH1 pathway which controls cell proliferation were included in the selected pathways.

Fig. 3
figure 3

Classifying ER+/PR + and ER+/PR- using 1387 pathways with LASSO methodology a Plots for LASSO regression coefficients over different values of the penalty parameter. b Cross-validation plot for the penalty term. c Receiver operating characteristic curve of the final model

Table 3 Pathways selected in the final model by LASSO

EGFR and IGF-1 pathway were deferentially regulated in ER+/PR + and ER+/PR- breast cancer

In the 13 selected pathways used in the final model, three of them were associated with growth factors including epidermal growth factor receptor (EGFR) and insulin-like growth factor-1 (IGF-1), indicating the direct association between growth factor pathway activation and PR loss. Previous studies have shown that in breast cancer cell lines, activation of EGF or IGF-1 pathways can sharply lower the expression of PR [6]. However, in this study, it was found that the EGFR pathway was more activated in the PR- group while the IGF-1 pathway was more activated in the PR + group (Fig. 4).

Fig. 4
figure 4

Pathway scores of EGFR and IGF-1 pathways between the PR + and PR- groups. Box plot of the pathway scores in three selected pathways by LASSO between the two groups

The correlation between the pathway score with the mRNA expression value of the PGR gene was plotted (Fig. 5). The IGF-1 pathway was found to be positively correlated with the expression of PGR while the EGFR pathway was negatively correlated, consistent with our above findings. To further look at the impact of the three growth factor pathways on the prognosis, survival analysis was performed between the pathway-defined subtypes by classifying the studied TCGA population into high-activity and low-activity groups according to the value of pathway score (Fig. 6). In the three analyzed pathways, only the IGF1R pathway-defined two groups showed survival difference while the other two pathways did not show a significant difference which was not surprising due to the fact that the IGF1R pathway score showed the greatest difference between PR + and PR- group in Fig. 4. The insignificance of prognosis between the IGF-1 and EGFR defined groups could possibly be explained by the complexity of the PR loss mechanism which was contributed by a network of pathways.

Fig. 5
figure 5

Correlation between growth factor pathway activity with PR expression. Correlation plot of the pathway activity in EGFR and IGF-1 pathways with the mRNA expression value of PGR gene in all the samples. Regression line, correlation co-efficiency and p-value were also displayed

Fig. 6
figure 6

Survival analysis of the pathway activity defined subgroups. Kaplan-Meier survival curves were plotted and the log-rank test was performed to compare the overall survival of low and high activity subgroups as defined by the pathway activity of the IGF1R a, IGF-1 b and EGFR c pathway respectively


Although PR loss in ER+/HER2- breast cancer had a worse prognosis and showed insensitivity towards SERM therapy, the molecular mechanism related to PR loss remains to be equivocal. In this study, we analyzed ER+/PR + and ER+/PR- breast cancer in the TCGA cohort. The PR- group was found to be more enriched in basal-like mRNA subtype than the PR + group. Also, in terms of DNA methylation and copy number alteration features, the PR- group was also more similar to the basal-like subtype than the PR + group. Pathway activities which are calculated from transcriptomic and gene copy number data were comprehensively analyzed and compared between the two groups. The PR + and PR- groups could be classified using the pathway activities of 13 selected pathways. Growth factor pathways including EGFR and IGF-1 were found to be essential for the distinction of the two groups which agrees with previous studies. Although previous studies showed that both EGFR and IGF-1 pathway could suppress the PR expression in a way independent of ER [6, 19], our results found that while IGF-1 pathway activity was negatively correlated with the expression of PR, the EGFR pathway activity showed the opposite. Expression of PR was previously found to be repressed by activation of the IGF-1 pathway through PI3K/Akt/mTOR signaling which was independent of ER [6]. Also, genes in PI3K/Akt/mTOR signaling pathways were upregulated in ER+/PR- breast cancer patients compared with ER+/PR + patients [20]. For the EGF pathway, studies were showing that activation of EGFR by adding EGF could enhance the phosphorylation of PR and promote its transcription activity [21], which could explain the positive correlation identified in our result between the EGFR pathway and PR expression. Other important pathways contributing to the classifying were ion-channel transporting pathway, immune-modulating pathway, and cell metabolism pathway. The above findings suggest that in ER+/HER2- breast cancers the various status of PR expression can be an indication of molecular variation, particularly for the growth factor pathway activation. Although, current clinical practice treats ER+/PR+/HER2- and ER+/PR-/HER2- breast cancer patients in the same way, our study suggests that the PR- group should be more intensively studied to search for an effective therapy due to its poor prognosis. Our study, together with other researches, suggested that growth factor pathways such as IGF-1 pathway could be served as potential targets.


Data acquisition and pre-processing

UCSC Xena, an online exploration tool for public and private, multi-omic, and clinical/phenotype data, was used to download data of selected samples [22]. The ‘TCGA Breast Cancer (BRCA)’ cohort in the UCSC Xena was selected. All raw data used to generate Fig. 1; Table 1 was downloaded from the ‘Phenotypes TCGA Hub’, pathway scores were downloaded from the ‘z score of 1387 constituent PARADIGM pathways TCGA Hub’. The pathway scores were generated by the TCGA [23], using the PARADIGM-inferred activation of pathway features [24]. Surveillance, Epidemiology, and End Results (SEER) 18 registries research database (Nov 2018 Submission) was used for the analysis which includes cases diagnosed from 1975 to 2016 [25].

Feature selection and LASSO regression

Regularized regression is often applied in genetic studies of molecular phenotypes to select the most promising set of variants associated with a phenotype of interest [26]. A widely applied regularized regression method is LASSO, which adds a penalty term for the shrinkage of the parameter estimates to the least-squares loss function [27].

Feature selection and LASSO regression were performed using the “glmnet” R package [28]. The ‘cv.glmnet’ function was used to do the k-fold cross-validation for glmnet, and the lambda.1se of the result was used as the optimal lambda. The glmnet function was used to fit a generalized linear model via penalized maximum likelihood. The alpha penalty was set to be 1 and the family was set to be ‘binomial’ in the glmnet function. The ‘coef.glmnet’ function was used to extract coefficients from the object generated by the glmnet function, with the ‘s’ argument being the optimal lambda generated by the cv.glmnet function.

Statistical analysis

Survival analyses were performed using the ‘survival’ (version 2.41) package [14]. The Kaplan-Meier method was used to estimate the survival outcomes of all patients by different categories; groups were compared using the log-rank statistic [29]. For univariate and multivariate analysis, the prognostic values of different clinicopathological factors were calculated by Cox proportional hazard modeling. The optimum threshold for distinguishing low and high activity groups was in similar method as described in our previous study [30]. P values were calculated as two-sided, with statistical significance declared for P less than 0.05.

Availability of data and materials

The dataset supporting the conclusions of this article is included within the manuscript.



Area Under the receiver operation characteristic Curve.


Epidermal Growth Factor Receptor.


Estrogen Receptor.


Human Epidermal growth factor Receptor 2.


Infiltrating Ductal Carcinoma.


Insulin-like Growth Factor-1.


Infiltrating Lobular Carcinoma.


Least Absolute Shrinkage and Selection Operator.


National Cancer Database.


progesterone receptor.


Surveillance, Epidemiology, and End Results.


Selective Estrogen Receptor Modulator.


The Cancer Genome Atlas.


  1. Ignatiadis M, Sotiriou C. Luminal breast cancer: from biology to treatment. Nat Rev Clin Oncol. 2013;10(9):494.

    Article  CAS  Google Scholar 

  2. Perou CM, Sørlie T, Eisen MB, Van De Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.

    Article  CAS  Google Scholar 

  3. Jacobsen BM, Horwitz KB. Progesterone receptors, their isoforms and progesterone regulated transcription. Mol Cell Endocrinol. 2012;357(1–2):18–29.

    Article  CAS  Google Scholar 

  4. Li Y, Yang D, Yin X, Zhang X, Huang J, Wu Y, Wang M, Yi Z, Li H, Li H. Clinicopathological characteristics and breast cancer–specific survival of patients with single hormone receptor–positive breast cancer. JAMA Netw Open. 2020;3(1):e1918160.

    Article  Google Scholar 

  5. Vienonen A, Syvala H, Miettinen S, Tuohimaa P, Ylikomi T. Expression of progesterone receptor isoforms A and B is differentially regulated by estrogen in different breast cancer cell lines. J Steroid Biochem Mol Biol. 2002;80(3):307–13.

    Article  CAS  Google Scholar 

  6. Cui X, Zhang P, Deng W, Oesterreich S, Lu Y, Mills GB, Lee AV. Insulin-like growth factor-I inhibits progesterone receptor expression in breast cancer cells via the phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin pathway: progesterone receptor as a potential indicator of growth factor activity in breast cancer. Mol Endocrinol. 2003;17(4):575–88.

    Article  CAS  Google Scholar 

  7. Yang C, Chen L, Li C, Lynch MC, Brisken C, Schmidt EV. Cyclin D1 enhances the response to estrogen and progesterone by regulating progesterone receptor expression. Mol Cell Biol. 2010;30(12):3111–25.

    Article  Google Scholar 

  8. Dauphine C, Moazzez A, Neal JC, Chlebowski RT, Ozao-Choy J. Single Hormone Receptor-Positive Breast Cancers Have Distinct Characteristics and Survival. Ann Surg Oncol. 2020;27(12):4687–94.

    Article  Google Scholar 

  9. Boland M, Ryan É, Dunne E, Aherne T, Bhatt N, Lowery A. Meta-analysis of the impact of progesterone receptor status on oncological outcomes in oestrogen receptor‐positive breast cancer. Br J Surg. 2020;107(1):33–43.

    Article  CAS  Google Scholar 

  10. Cui X, Schiff R, Arpino G, Osborne CK, Lee AV. Biology of progesterone receptor loss in breast cancer and its implications for endocrine therapy. J Clin Oncol. 2005;23(30):7721–35.

    Article  CAS  Google Scholar 

  11. Daniel AR, Gaviglio AL, Knutson TP, Ostrander JH, D’Assoro AB, Ravindranathan P, Peng Y, Raj GV, Yee D, Lange CA. Progesterone receptor-B enhances estrogen responsiveness of breast cancer cells via scaffolding PELP1-and estrogen receptor-containing transcription complexes. Oncogene. 2015;34(4):506–15.

    Article  CAS  Google Scholar 

  12. Balleine R, Earl M, Greenberg M, Clarke C. Absence of progesterone receptor associated with secondary breast cancer in postmenopausal women. Br J Cancer. 1999;79(9):1564–71.

    Article  CAS  Google Scholar 

  13. Bamberger A-M, Milde-Langosch K, Schulte HM, Löning T. Progesterone receptor isoforms, PR-B and PR-A, in breast cancer: correlations with clinicopathologic tumor parameters and expression of AP-1 factors. Hormone Res Paediatr. 2000;54(1):32–7.

    Article  CAS  Google Scholar 

  14. Liu X-Y, Ma D, Xu X-E, Jin X, Yu K-D, Jiang Y-Z, Shao Z-M. Genomic landscape and endocrine-resistant subgroup in estrogen receptor-positive, progesterone receptor-negative, and HER2-negative breast cancer. Theranostics. 2018;8(22):6386.

    Article  CAS  Google Scholar 

  15. Horowitz K, McGuire WL. Predicting response to endocrine therapy in human breast cancer: a hypothesis. Science. 1975;189(4204):726–7.

    Article  Google Scholar 

  16. Mohammed H, Russell IA, Stark R, Rueda OM, Hickey TE, Tarulli GA, Serandour AA, Birrell SN, Bruna A, Saadi A. Progesterone receptor modulates ERα action in breast cancer. Nature. 2015;523(7560):313–7.

    Article  CAS  Google Scholar 

  17. Allison KH, Hammond MEH, Dowsett M, McKernin SE, Carey LA, Fitzgibbons PL, Hayes DF, Lakhani SR, Chavez-MacGregor M, Perlmutter J. Estrogen and progesterone receptor testing in breast cancer: American society of clinical oncology/college of American pathologists guideline update. Arch Pathol Lab Med. 2020;144(5):545–63.

    Article  CAS  Google Scholar 

  18. Network CGA. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61.

    Article  Google Scholar 

  19. Thakkar JP, Mehta DG. A review of an unfavorable subset of breast cancer: estrogen receptor positive progesterone receptor negative. Oncologist. 2011;16(3):276.

    Article  CAS  Google Scholar 

  20. Creighton CJ, Osborne CK, van de Vijver MJ, Foekens JA, Klijn JG, Horlings HM, Nuyten D, Wang Y, Zhang Y, Chamness GC. Molecular profiles of progesterone receptor loss in human breast tumors. Breast Cancer Res Treat. 2009;114(2):287–99.

    Article  CAS  Google Scholar 

  21. Daniel AR, Qiu M, Faivre EJ, Ostrander JH, Skildum A, Lange CA. Linkage of progestin and epidermal growth factor signaling: phosphorylation of progesterone receptors mediates transcriptional hypersensitivity and increased ligand-independent breast cancer cell growth. Steroids. 2007;72(2):188–201.

    Article  CAS  Google Scholar 

  22. Goldman MJ, Craft B, Hastie M, Repečka K, McDade F, Kamath A, Banerjee A, Luo Y, Rogers D, Brooks AN: Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol 2020:1–4.

  23. Hoadley KA, Yau C, Hinoue T, Wolf DM, Lazar AJ, Drill E, Shen R, Taylor AM, Cherniack AD, Thorsson V. Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of cancer. Cell. 2018;173(2):291–304 e296.

    Article  CAS  Google Scholar 

  24. Vaske CJ, Benz SC, Sanborn JZ, Earl D, Szeto C, Zhu J, Haussler D, Stuart JM. Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM. Bioinformatics. 2010;26(12):i237–45.

  25. Howlader N, Noone AM, Krapcho M, Miller D, Brest A, Yu M, Ruhl J, Tatalovich Z, Mariotto A, Lewis D. SEER cancer statistics review, 1975–2016. Natl Cancer Institut. 2019. Online:

  26. Deutelmoser H, Scherer D, Brenner H, Waldenberger M, Suhre K, Kastenmüller G, Lorenzo Bermejo J. Robust Huber-LASSO for improved prediction of protein, metabolite and gene expression levels relying on individual genotype data. Brief Bioinform. 2020;bbaa230.

  27. Tibshirani R. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series B (Methodological). 1996;58(1):267–88.

    Google Scholar 

  28. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1.

    Article  Google Scholar 

  29. Bland JM, Altman DG. The logrank test. BMJ. 2004;328(7447):1073.

    Article  Google Scholar 

  30. Long M, Hou W, Liu Y, Hu T. A Histone Acetylation Modulator Gene Signature for Classification and Prognosis of Breast Cancer. Curr Oncol. 2021;28(1):928–39.

    Article  Google Scholar 

Download references




The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (Grant No. 82002979), Scientific Research and Development Funds of Peking University People’s Hospital (Grant No. RDY2020-16) and Peking University Medicine Fund of Fostering Young Scholars’ Scientific & Technological Innovation supported by “the Fundamental Research Funds for the Central Universities” (Grant No. BMU2020PYB022, BMU2021PYB013).

Author information

Authors and Affiliations



TH and ML contributed to the conception of the study; TH, YC and JP performed the data analyses; YL and DZ contributed to the methodology; YL, DZ and JP helped the manuscript writing with constructive discussion; ML wrote the final manuscript. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Jiankang Pan or Mengping Long.

Ethics declarations

Ethics approval and consent to participate

As this study relied on publicly available data, no ethics approval was required.

Consent for publication

Not applicable.

Competing interests

None. All data used in this publication were generated by the TCGA project. The authors declare no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Table S1. Clinicopathological features of the two selected groups in the SEER dataset. Table S2. Univariate and multivariate analysis of the prognostic value of clinicopathological features in ER+/HER2- patients in the SEER dataset.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hu, T., Chen, Y., Liu, Y. et al. Classification of PR-positive and PR-negative subtypes in ER-positive and HER2-negative breast cancers based on pathway scores. BMC Med Res Methodol 21, 108 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: