In this study, we measured and analyzed the expression of 750 tumor-related genes in 34 breast fibroepithelial lesions, as a potential aid in their classification and to explore their underlying biology. Although fibroadenomas and benign phyllodes tumors are recognized as distinct entities, their common overlapping histologic features (particularly between cellular FA and benign PT) makes the distinction challenging in daily practice. There are several comprehensive genetic and proteomic studies of FELs suggests that FAs and benign PTs share similar gene and protein expression profile(9, 10, 18, 23). Our study revealed that fibroadenomas (FAs), cellular fibroadenomas, and benign phyllodes tumors (PTs) manifest similar gene expression profiles. Genes that are involved in different cancer-related pathways like MAPK, epigenetic regulation, matrix remodeling and metastasis, angiogenesis, TGF-beta signaling, Hedgehog signaling, and PI3K-Akt showed no difference when performing gene set analysis between these three groups. Similar patterns were also recognized in these three groups of FELs using cell type analysis. While the current treatment protocol for asymptomatic fibroadenomas is conservative management (such as annual breast examination and ultrasound as clinically indicated(24)), benign phyllodes tumors are treated with surgical intervention in many medical centers(25, 26). The overlapping histologic and gene expression profiles of FAs, cellular FAs, and benign phyllodes tumors could imply a similarly indolent clinical behavior and may suggest that a conservative management strategy may be appropriate for each of these entities.
A recent study revealed that borderline phyllodes tumor could be separated into either benign or malignant PTs using molecular assays(27). Interestingly, in our study, the hierarchical clustering heatmap show several borderline phyllodes cases clustered with the malignant PTs group, and some clustered closely with the benign phyllodes tumor group (supplement Fig. 2). This suggests that a subset of borderline PTs exhibits a gene expression profile which overlaps with that of malignant PTs. Various gene alterations were found both in malignant PTs and borderline PTs (supplement Fig. 3), which could imply that some malignant PTs derive from a subset of borderline PTs through clonal evolution and/or acquisition of somatic mutations leading to malignant transformation Prospective studies utilizing gene expression analysis with larger sample size could aid in optimizing the diagnostic and prognostic stratification of borderline PTs. Additional studies are also required to determine whether the remaining subset of borderline PTs lacking a gene expression profile comparable to that of malignant PTs may also progress to malignancy via de novo pathways.
Extracellular matrix-related gene alterations are known to promote metastases in various cancers(28–30). The ECM molecules are major components in the tumor microenvironment integral for cellular adhesion and harboring growth-factors(31). Fibroblasts secrete collagen which is the main structural protein of the ECM responsible for intercellular adhesion, differentiation, and integrity(32). Previous studies report ECM gene dysregulation such as, but not limited to, COL1A1, FN1, TIMP1, and MMP9, to play an important role in tumorigenesis in various malignancies and are associated with poor prognosis(31, 33–35). In this study, we also found higher expressions of collagen type I alpha 1 chain (COL1A1) as well as matrix metallopeptidases 7 and 9 (MMP7, MMP9) genes in the borderline and malignant PTs groups when compared with benign PTs. Further studies are required to construct and validate a comprehensive molecular panel composed of ECM-related genes might aid with proper risk stratification of borderline PTs.
Our gene expression profiling results identified many important cancer-related biological processes such as angiogenesis, matrix remodeling & metastasis, PI3K-Akt pathways, TGF-beta signaling, MAPK pathways, and hypoxia which discriminate malignant PTs from other FELs. This is consistent with our current understanding of the pathogenesis and biologic behavior of malignant PTs, where stromal overgrowth, increased angiogenesis, and extracellular matrix factors are profoundly involved in tumor progression(9). Pathway score and gene set analysis revealed that angiogenesis-related genes including VEGFA, ITGAV, NFIL3, FDFR1, and CCND2 are highly upregulated in malignant PTs. Recent studies have demonstrated that VEGFRs are expressed in several cancer cell types and may dictate tumor invasion(36). This finding could potentially predict response to the targeted therapy (angiogenesis inhibitor). Overexpression of EGFR and PDGFA are also present in malignant phyllodes tumors, which suggests the potential use of dual VEGFR/EGFR inhibitors (Vandetanib), a combination of anti-VEGFR and anti-EGFR and multikinase inhibitor (Panzopnib).
The growth and survival of malignant cells are often driven by constitutive activation in the mitogen-activated protein kinase (MAPK) and phosphor-inositide 3-kinase (PI3K/AKT) signaling pathways(37–39). VEGFR-2 mediated activation of PI3K/Alk cascade is also important for tumor survival. It is well established that no pathways exist independently, the crosstalk between the pathways promotes tumor cell proliferation, survival, and invasion. Gene set analysis (GAS) revealed that 28 genes (PIK3R1, PRLR, EIF4EBP1, NRAS, HRAS, MAP3K12, KIT), involved in PI3K/AKT and MAPK pathways, were differentially expressed in malignant PTs group. Moreover, genes in NK-kappa B (IKBKB, KIBKG, NFKB1) and NOTCH (NOTCH1, HDAC11, APH1B) signaling pathways are also differentially expressed. PIK3CA/KRAS and HRAS alterations have been already reported in borderline and malignant phyllodes tumor(40). In contrast to the previous study, we identified high frequency and more gene alterations in comparison with the previous data which implies a potential therapeutic benefit from IKK inhibitors and mTOR inhibitors.
An increased macrophage abundance score in malignant phyllodes tumors was found in this cohort by cell type analysis. The score is calculated by measuring the expression level of four macrophage-related genes: CD163, CD68, CD84, and MS4A4A. We found CD68, CD84, and CD163 to be upregulated in malignant PTs. Many studies have found that breast cancer with high tumor-associated macrophage (TAM) infiltration was significantly correlated with aggressive biological behavior and could serve as an independent predictor of overall survival and recurrence-free survival(41, 42). Previous studies also proposed that TAMs could stimulate myofibroblast differentiation and promote proliferation as well as invasion of phyllodes tumor by CCL18 driven NF-kB/PTEN/AKT axis(43). The interactions of the tumor microenvironment and tumor cells have been shown to drive the progression of various cancers, and TAMs are the most abundant inflammatory cell type within the microenvironment of malignant PTs(44). Additional TAM-associated markers are also identified in our study by differential expression analysis including CCL18, CCL5, MMP-9, and SPP-1(45). Consistent with prior signaling pathway alteration results(46–48), our data also suggest TAMs play a crucial role in the tumorigenesis of malignant phyllodes tumors. However, the specific pathways of TAMs interacting with myofibroblasts and other inflammatory cells leading to malignant progression would require further investigation.
Our analysis was limited by a relatively small cohort size. Future studies with larger sample size are needed to validate our findings. In addition, correlating the molecular result with long term clinical follow-up data might provide relevant information to guide clinical management.
In conclusion, this study provided a comprehensive profile of gene alterations across various subtypes of breast fibroepithelial lesions and might bring new insights in their classification and improve understanding of their pathogenesis. Our results raise the possibility that fibroadenoma, cellular fibroadenoma and benign phyllodes tumor might be merged in the same subgroup for the purpose of clinical management. Further studies would be of particular interest to predict the biological behavior of borderline phyllodes tumors.