Dynamic altruistic cooperation within breast tumors.

BACKGROUND: Social behaviors such as altruism, where one self-sacrifices for collective benefits, critically influence an organism's survival and responses to the environment. Such behaviors are widely exemplified in nature but have been underexplored in cancer cells which are conventionally seen as selfish competitive players. This multidisciplinary study explores altruism and its mechanism in breast cancer cells and its contribution to chemoresistance. METHODS: MicroRNA profiling was performed on circulating tumor cells collected from the blood of treated breast cancer patients. Cancer cell lines ectopically expressing candidate miRNA were used in co-culture experiments and treated with docetaxel. Ecological parameters like relative survival and relative fitness were assessed using flow cytometry. Functional studies and characterization performed in vitro and in vivo include proliferation, iTRAQ-mass spectrometry, RNA sequencing, inhibition by small molecules and antibodies, siRNA knockdown, CRISPR/dCas9 inhibition and fluorescence imaging of promoter reporter-expressing cells. Mathematical modeling based on evolutionary game theory was performed to simulate spatial organization of cancer cells. RESULTS: Opposing cancer processes underlie altruism: an oncogenic process involving secretion of IGFBP2 and CCL28 by the altruists to induce survival benefits in neighboring cells under taxane exposure, and a self-sacrificial tumor suppressive process impeding proliferation of altruists via cell cycle arrest. Both processes are regulated concurrently in the altruists by miR-125b, via differential NF-κB signaling specifically through IKKß. Altruistic cells persist in the tumor despite their self-sacrifice, as they can regenerate epigenetically from non-altruists via a KLF2/PCAF-mediated mechanism. The altruists maintain a sparse spatial organization by inhibiting surrounding cells from adopting the altruistic fate via a lateral inhibition mechanism involving a GAB1-PI3K-AKT-miR-125b signaling circuit. CONCLUSIONS: Our data reveal molecular mechanisms underlying manifestation, persistence and spatial spread of cancer cell altruism. A minor population behave altruistically at a cost to itself producing a collective benefit for the tumor, suggesting tumors to be dynamic social systems governed by the same rules of cooperation in social organisms. Understanding cancer cell altruism may lead to more holistic models of tumor evolution and drug response, as well as therapeutic paradigms that account for social interactions. Cancer cells constitute tractable experimental models for fields beyond oncology, like evolutionary ecology and game theory.

CXCR4 expression is elevated in TNBC patient derived samples and Z-guggulsterone abrogates tumor progression by targeting CXCL12/CXCR4 signaling axis in preclinical breast cancer model.

Environmental factors such as exposure to ionizing radiations, certain environmental pollutants, and toxic chemicals are considered as risk factors in the development of breast cancer. Triple-negative breast cancer (TNBC) is a molecular variant of breast cancer that lacks therapeutic targets such as progesterone receptor, estrogen receptor, and human epidermal growth factor receptor-2 which makes the targeted therapy ineffective in TNBC patients. Therefore, identification of new therapeutic targets for the treatment of TNBC and the discovery of new therapeutic agents is the need of the hour. In this study, CXCR4 was found to be highly expressed in majority of breast cancer tissues and metastatic lymph nodes derived from TNBC patients. CXCR4 expression is positively correlated with breast cancer metastasis and poor prognosis of TNBC patients suggesting that suppression of CXCR4 expression could be a good strategy in the treatment of TNBC patients. Therefore, the effect of Z-guggulsterone (ZGA) on the expression of CXCR4 in TNBC cells was examined. ZGA downregulated protein and mRNA expression of CXCR4 in TNBC cells and proteasome inhibition or lysosomal stabilization had no effect on the ZGA-induced CXCR4 reduction. CXCR4 is under the transcriptional control of NF-κB, whereas ZGA was found to downregulate transcriptional activity of NF-κB. Functionally, ZGA downmodulated the CXCL12-driven migration/invasion in TNBC cells. Additionally, the effect of ZGA on growth of tumor was investigated in the orthotopic TNBC mice model. ZGA presented good inhibition of tumor growth and liver/lung metastasis in this model. Western blotting and immunohistochemical analysis indicated a reduction of CXCR4, NF-κB, and Ki67 in tumor tissues. Computational analysis suggested PXR agonism and FXR antagonism as targets of ZGA. In conclusion, CXCR4 was found to be overexpressed in majority of patient-derived TNBC tissues and ZGA abrogated the growth of TNBC tumors by partly targeting the CXCL12/CXCR4 signaling axis.

Profiling Microbial Communities in Idiopathic Granulomatous Mastitis.

Idiopathic granulomatous mastitis (IGM) is a rare and benign inflammatory breast disease with ambiguous aetiology. Contrastingly, lactational mastitis (LM) is commonly diagnosed in breastfeeding women. To investigate IGM aetiology, we profiled the microbial flora of pus and skin in patients with IGM and LM. A total of 26 patients with IGM and 6 patients with LM were included in the study. The 16S rRNA sequencing libraries were constructed from 16S rRNA gene amplified from total DNA extracted from pus and skin swabs in patients with IGM and LM controls. Constructed libraries were multiplexed and paired-end sequenced on HiSeq4000. Metagenomic analysis was conducted using modified microbiome abundance analysis suite customised R-resource for paired pus and skin samples. Microbiome multivariable association analyses were performed using linear models. A total of 21 IGM and 3 LM paired pus and skin samples underwent metagenomic analysis. Bray−Curtis ecological dissimilarity distance showed dissimilarity across four sample types (IGM pus, IGM skin, LM pus, and LM skin; PERMANOVA, p < 0.001). No characteristic dominant genus was observed across the IGM samples. The IGM pus samples were more diverse than corresponding IGM skin samples (Shannon and Simpson index; Wilcoxon paired signed-rank tests, p = 0.022 and p = 0.07). Corynebacterium kroppenstedtii, reportedly associated with IGM in the literature, was higher in IGM pus samples than paired skin samples (Wilcoxon, p = 0.022). Three other species and nineteen genera were statistically significant in paired IGM pus−skin comparison after antibiotic treatment adjustment and multiple comparisons correction. Microbial profiles are unique between patients with IGM and LM. Inter-patient variability and polymicrobial IGM pus samples cannot implicate specific genus or species as an infectious cause for IGM.

Photoacoustic Tomography Appearance of Fat Necrosis: A First-in-Human Demonstration of Biochemical Signatures along with Histological Correlation.

A 50-year-old woman with no past medical history presented with a left anterior chest wall mass that was clinically soft, mobile, and non-tender. A targeted ultrasound (US) showed findings suggestive of a lipoma. However, focal "mass-like" nodules seen within the inferior portion suggested malignant transformation of a lipomatous lesion called for cross sectional imaging, such as MRI or invasive biopsy or excision for histological confirmation. A T1-weighted image demonstrated a large lipoma that has a central fat-containing region surrounded by an irregular hypointense rim in the inferior portion, confirming the benignity of the lipoma. An ultrasound-guided photoacoustic imaging (PA) of the excised specimen to derive the biochemical distribution demonstrated the "mass-like" hypoechoic regions on US as fat-containing, suggestive of benignity of lesion, rather than fat-replacing suggestive of malignancy. The case showed the potential of PA as an adjunct to US in improving the diagnostic confidence in lesion characterization.

αvß1 integrin is enriched in extracellular vesicles of metastatic breast cancer cells: A mechanism mediated by galectin-3.

Breast cancer cells release a large quantity of biocargo-bearing extracellular vesicles (EVs), which mediate intercellular communication within the tumour microenvironment and promote metastasis. To identify EV-bound proteins related to metastasis, we used mass spectrometry to profile EVs from highly and poorly metastatic breast cancer lines of human and mouse origins. Comparative mass spectrometry indicated that integrins, including αv and ß1 subunits, are preferentially enriched in EVs of highly metastatic origin over those of poorly metastatic origin. These results are consistent with our histopathological findings, which show that integrin αv is associated with disease progression in breast cancer patients. Integrin αv colocalizes with the multivesicular-body marker CD63 at a higher frequency in the tumour and is enriched in circulating EVs of breast cancer patients at late stages when compared with circulating EVs from early-stage patients. With a magnetic bead-based flow cytometry assay, we confirmed that integrins αv and ß1 are enriched in the CD63+ subsets of EVs from both human and mouse highly metastatic cells. By analysing the level of integrin αv on circulating EVs, this assay could predict the metastatic potential of a xenografted mouse model. To explore the export mechanism of integrins into EVs, we performed immunoprecipitation mass spectrometry and identified members of the galectin family as potential shuttlers of integrin αvß1 into EVs. In particular, knockdown of galectin-3, but not galectin-1, causes a reduction in the levels of cell surface integrins ß1 and αv, and decreases the colocalization of these integrins with CD63. Importantly, knockdown of galectin-3 leads to a decrease of integrin αvß1 export into the EVs concomitant with a decrease in the metastatic potential of breast cancer cells. Moreover, inhibition of the integrin αvß1 complex leads to a reduction in the binding of EVs to fibronectin, suggesting that integrin αvß1 is important for EV retention in the extracellular matrix. EVs retained in the extracellular matrix are taken up by fibroblasts, which differentiate into cancer associated fibroblasts. In summary, our data indicate an important link between EV-bound integrin αvß1 with breast cancer metastasis and provide additional insights into the export of integrin αvß1 into EVs in the context of metastasis.

Biochemical "decoding" of breast ultrasound images with optoacoustic tomography fusion: First-in-human display of lipid and collagen signals on breast ultrasound.

To date, studies which utilized ultrasound (US) and optoacoustic tomography (OT) fusion (US-OT) in biochemical differentiation of malignant and benign breast conditions have relied on limited biochemical data such as oxyhaemoglobin (OH) and deoxyhaemoglobin (DH) only. There has been no data of the largest biochemical components of breast fibroglandular tissue: lipid and collagen. Here, the authors believe the ability to image collagen and lipids within the breast tissue could serve as an important milestone in breast US-OT imaging with many potential downstream clinical applications. Hence, we would like to present the first-in-human US-OT demonstration of lipid and collagen differentiation in an excised breast tissue from a 38-year-old female.

Infiltrative pattern of metastatic invasive lobular breast carcinoma in the abdomen: a pictorial review.

Invasive lobular carcinoma (ILC) has a greater tendency to metastasize to the peritoneum, retroperitoneum, and gastrointestinal (GI) tract as compared to invasive carcinoma of no special type (NST). Like primary ILC in the breast, ILC metastases are frequently infiltrative and hypometabolic, rather than mass forming and hypermetabolic in nature. This renders them difficult to detect on conventional and metabolic imaging studies. As a result, intra-abdominal ILC metastases are often detected late, with patients presenting with clinical complications such as liver failure, hydronephrosis, or bowel obstruction. In patients with known history of ILC, certain imaging features are very suggestive of infiltrative metastatic ILC. These include retroperitoneal or peritoneal nodularity and linitis plastica appearance of the bowel. Recognition of linitis plastica on imaging should prompt deep or repeat biopsies. In this pictorial review, the authors aim to familiarize readers with imaging features and pitfalls for evaluation of intra-abdominal metastatic ILC. Awareness of these will allow the radiologist to assess these patients with a high index of suspicion and aid detection of metastatic disease. Also, this can direct histopathology and immunohistochemical staining to obtain the correct diagnosis in suspected metastatic disease.

Optoacoustic characterization of breast conserving surgery specimens - A pilot study.

In this pilot study, we tested an ultrasound-guided optoacoustic tomography (US-OT) two-dimensional (2D) array scanner to understand the optoacoustic patterns of excised breastconserving surgery (BCS) specimens. We imaged 14 BCS specimens containing malignant tumors at eight wavelengths spanning 700-1100 nm. Spectral unmixing across multiple wavelengths allowed for visualizing major intrinsic chromophores in the breast tissue including hemoglobin and lipid up to a depth of 7 mm. We identified less/no lipid signals within the tumor and intense deoxy-hemoglobin (Hb) signals on the rim of the tumor as unique characteristics of malignant tumors in comparison to no tumor region. We also observed continuous broad lipid signals as features of negative margins and compromised lipid signals interrupted by vasculature as features of positive margins. These differentiating patterns can form the basis of US-OT to be explored as an alternate, fast and efficient intraoperative method for evaluation of tumor resection margins.

Primary Pleomorphic Lipoma of the Parotid Gland with Prominent Myxoid Change: Report of a Rare Case Mimicking Carcinoma Ex Pleomorphic Adenoma on Fine Needle Aspiration Cytology.

We present a case of a primary 2.5 cm pleomorphic lipoma of the right parotid gland with prominent myxoid change which on FNA displayed features that mimicked a carcinoma or sarcoma ex pleomorphic adenoma. The patient was a 79 year old man with no oncological history or tumor elsewhere. On immunohistochemistry the neoplastic cells strongly expressed CD34. There was no expression of retinoblastoma protein, smooth muscle actin, S100-protein or cytokeratins (AE1/3 and CAM5.2). The Ki-67 proliferation index was low (< 2%). Fluorescence in situ hybridization was negative for MDM2 gene amplification.

Ultrasound Guided Optoacoustic Tomography in Assessment of Tumor Margins for Lumpectomies.

PURPOSE: To determine the accuracy of a handheld ultrasound-guided optoacoustic tomography (US-OT) probe developed for human deep-tissue imaging in ex vivo assessment of tumor margins postlumpectomy. METHODS: A custom-built two-dimensional (2D) US-OT-handheld probe was used to scan 15 lumpectomy breast specimens. Optoacoustic signals acquired at multiple wavelengths between 700 and 1100 nm were reconstructed using model linear algorithm, followed by spectral unmixing for lipid and deoxyhemoglobin (Hb). Distribution maps of lipid and Hb on the anterior, posterior, superior, inferior, medial, and lateral margins of the specimens were inspected for margin involvement, and results were correlated with histopathologic findings. The agreement in tumor margin assessment between US-OT and histopathology was determined using the Bland-Altman plot. Accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) of margin assessment using US-OT were calculated. RESULTS: Ninety margins (6 × 15 specimens) were assessed. The US-OT probe resolved blood vessels and lipid up to a depth of 6 mm. Negative and positive margins were discriminated by marked differences in the distribution patterns of lipid and Hb. US-OT assessments were concordant with histopathologic findings in 87 of 89 margins assessed (one margin was uninterpretable and excluded), with diagnostic accuracy of 97.9% (kappa = 0.79). The sensitivity, specificity, PPV, and NPV were 100% (4/4), 97.6% (83/85), 66.7% (4/6), and 100% (83/83), respectively. CONCLUSION: US-OT was capable of providing distribution maps of lipid and Hb in lumpectomy specimens that predicted tumor margins with high sensitivity and specificity, making it a potential tool for intraoperative tumor margin assessment.

Using next-generation sequencing (NGS) platform to diagnose pathogenic germline BRCA1/2 mutations from archival tumor specimens.

OBJECTIVE: Clinical genetic testing to diagnose germline mutations often requires blood sample or saliva smear from a cancer-affected individual. This rules out testing in families when cancer-affected individuals are deceased. We explored the use of a next-generation sequencing (NGS) platform to diagnose germline pathogenic mutations from tumors. METHODS: Archival tumors (ovarian = 26, breast = 25, others = 9) were retrieved from 60 cancer patients who have undergone multi-gene panel blood testing. Genomic DNA was extracted and sequenced for BRCA1/2 using a NGS platform. 41/60 specimens were sequenced for 5 other genes (APC, ATM, PALB2, PTEN, TP53). Tumor testing and results interpretation were performed blinded to the blood test result. RESULTS: All 38 patients with no BRCA1/2 mutations on blood testing were correctly tested negative on tumor. Tumor testing correctly diagnosed BRCA1/2 pathogenic mutations in 15/22 (68%) patients while in 7/22 (32%) patients, the mutation was either detected but incorrectly classified as VUS (n = 3) or not detected at all (n = 4). Overall concordance rate for tumor and blood testing for BRCA1/2 mutations was 88%, with 0% false positive and 32% false negative rate for pathogenic mutations. Tumor testing correctly diagnosed 1/2 pathogenic germline ATM mutation, 1/1 pathogenic germline PALB2 mutation and 2/2 pathogenic germline TP53 mutations. False positive germline mutations were diagnosed in 4 genes at a rate of 2.4%-10.3% (APC = 2.4%, PALB2 = 2.4%, PTEN = 4.9%, TP53 = 10.3%). CONCLUSION: Tumor testing for BRCA1/2 germline mutations using an NGS platform is fairly reliable with no false positive findings, and correctly diagnosed more than two-thirds of pathogenic germline BRCA1/2 mutations. However, it is not reliable to diagnose pathogenic germline mutations in genes frequently mutated in sporadic cancers, such as PTEN and TP53.

Elevated WBP2 Expression in HER2-positive Breast Cancers Correlates with Sensitivity to Trastuzumab-based Neoadjuvant Therapy: A Retrospective and Multicentric Study.

PURPOSE: Trastuzumab-based chemotherapy has shown remarkable clinical benefits for patients with HER2-positive breast cancer. However, treatment regimens involving trastuzumab had little or no effect for a subset of patients. Preliminary studies revealed WW-binding protein 2 (WBP2), an oncogenic transcription coactivator, to be coamplified with HER2 in 36% of HER2-positive breast cancers. We hypothesize that WBP2 regulates and correlates with the response of HER2-positive breast cancer to trastuzumab. EXPERIMENTAL DESIGN: The coexpression of WBP2 and HER2 in breast tumors was validated using IHC. The role and mechanism of WBP2 in regulating breast cancer response to trastuzumab was elucidated using in vitro, patient-derived xenograft and murine xenograft models. A multicenter retrospective study involving 143 patients given neoadjuvant trastuzumab-based chemotherapy was conducted to determine whether WBP2 expression correlates with pathologic complete response (pCR). RESULTS: Elevated expression of WBP2 significantly enhanced breast cancer's response to trastuzumab by augmenting trastuzumab-induced HER2 downregulation and cell-cycle arrest via inhibition of cyclin D expression. High level of WBP2 correlated with better pCR (67.19%) compared with low WBP2 level (26.58%). The highest response was observed in subgroups of patients with high WBP2-expressing tumors also aged below 50 years (77.78%) or were premenopausal in status (73.33%). Retrospectively, WBP2 demonstrated sensitivity of 80% to 81% and specificity of 76.5% to 80% in discriminating between patients showing pCR and non-pCR. CONCLUSIONS: WBP2 expression correlates with the response of HER2-positive breast cancer to trastuzumab-based neoadjuvant chemotherapy.

Phosphorylation of E-box binding USF-1 by PI3K/AKT enhances its transcriptional activation of the WBP2 oncogene in breast cancer cells.

WW domain binding protein 2 (WBP2), a transcriptional coactivator, plays a vital role in breast tumorigenesis. It positively regulates estrogen receptor, Hippo, and Wnt pathways, which subsequently enhance the transcription of downstream target genes contributing to cancer. Understanding the regulation of the expression and activity of WBP2 oncoprotein has implication in cancer therapy. We have previously reported that WBP2 is regulated at the post-translational and post-transcriptional levels. However, its regulation at the transcriptional level is not known. In this study, the minimal promoter region of WBP2 that is critical for its transcription was identified. The E-box motif in the WBP2 promoter was demonstrated to be essential for its transcription. The E-box binding protein upstream stimulatory factor 1 (USF-1) was discovered to be a key transcription factor for WBP2 by yeast one-hybrid analysis and was validated through reporter and chromatin immunoprecipitation assays and tandem mass spectrometry, which also suggested that USF-1 acts by regulating a network of genes, in addition to WBP2, associated with cell movement, proliferation, cell-cycle, and survival cellular processes. USF-1 is overexpressed in majority of the breast cancer cell lines and tissues tested, and has profound effects on cancer cell proliferation. USF-1-mediated transcription of WBP2 was demonstrated to be inducible by insulin, which led to AKT-mediated phosphorylation of USF-1 that modulated its ability to bind to the WBP2 promoter and activate its transcription. This study sheds new light onto the regulation of the WBP2 oncogene at the transcriptional level by a novel oncogenic transcription factor, USF-1. USF-1 is a potential drug target for treatment of WBP2-positive breast cancer.-Ramos, A., Miow, Q. H., Liang, X., Lin, Q. S., Putti, T. C., Lim, Y. P. Phosphorylation of E-box binding USF-1 by PI3K/AKT enhances its transcriptional activation of the WBP2 oncogene in breast cancer cells.

High Densities of Tumor-Associated Plasma Cells Predict Improved Prognosis in Triple Negative Breast Cancer.

Breast cancer is the most common malignancy affecting women, but the heterogeneity of the condition is a significant obstacle to effective treatment. Triple negative breast cancers (TNBCs) do not express HER2 or the receptors for estrogen or progesterone, and so often have a poor prognosis. Tumor-infiltrating T cells have been well-characterized in TNBC, and increased numbers are associated with better outcomes; however, the potential roles of B cells and plasma cells have been large. Here, we conducted a retrospective correlative study on the expression of B cell/plasma cell-related genes, and the abundance and localization of B cells and plasma cells within TNBCs, and clinical outcome. We analyzed 269 TNBC samples and used immunohistochemistry to quantify tumor-infiltrating B cells and plasma cells, coupled with NanoString measurement of expression of immunoglobulin metagenes. Multivariate analysis revealed that patients bearing TNBCs with above-median densities of CD38+ plasma cells had significantly better disease-free survival (DFS) (HR = 0.44; 95% CI 0.26-0.77; p = 0.004) but not overall survival (OS), after adjusting for the effects of known prognostic factors. In contrast, TNBCs with higher immunoglobulin gene expression exhibited improved prognosis (OS p = 0.029 and DFS p = 0.005). The presence of B cells and plasma cells was positively correlated (p < 0.0001, R = 0.558), while immunoglobulin gene IGKC, IGHM, and IGHG1 mRNA expression correlated specifically with the density of CD38+ plasma cells (IGKC p < 0.0001, R = 0.647; IGHM p < 0.0001, R = 0.580; IGHG1 p < 0.0001, R = 0.655). Interestingly, after adjusting the multivariate analysis for the effect of intratumoral CD38+ plasma cell density, the expression levels of all three genes lost significant prognostic value, suggesting a biologically important role of plasma cells. Last but not least, the addition of intratumoral CD38+ plasma cell density to clinicopathological features significantly increased the prognostic value for both DFS (ΔLRχ2 = 17.28, p = 1.71E-08) and OS (ΔLRχ2 = 10.03, p = 6.32E-08), compared to clinicopathological features alone. The best combination was achieved by integrating intratumoral CD38+ plasma cell density and IGHG1 which conferred the best added prognostic value for DFS (ΔLRχ2 = 27.38, p = 5.22E-10) and OS (ΔLRχ2 = 21.29, p = 1.03E-08). Our results demonstrate that the role of plasma cells in TNBC warrants further study to elucidate the relationship between their infiltration of tumors and disease recurrence.

Low Levels of NDRG1 in Nerve Tissue Are Predictive of Severe Paclitaxel-Induced Neuropathy.

INTRODUCTION: Sensory peripheral neuropathy caused by paclitaxel is a common and dose limiting toxicity, for which there are currently no validated predictive biomarkers. We investigated the relationship between the Charcot-Marie-Tooth protein NDRG1 and paclitaxel-induced neuropathy. METHODS/MATERIALS: Archived mammary tissue specimen blocks of breast cancer patients who received weekly paclitaxel in a single centre were retrieved and NDRG1 immunohistochemistry was performed on normal nerve tissue found within the sample. The mean nerve NDRG1 score was defined by an algorithm based on intensity of staining and percentage of stained nerve bundles. NDRG1 scores were correlated with paclitaxel induced neuropathy. RESULTS: 111 patients were studied. 17 of 111 (15%) developed severe paclitaxel-induced neuropathy. The mean nerve NDRG1 expression score was 5.4 in patients with severe neuropathy versus 7.7 in those without severe neuropathy (p = 0.0019). A Receiver operating characteristic (ROC) curve analysis of the mean nerve NDRG1 score revealed an area under the curve of 0.74 (p = 0.0013) for the identification of severe neuropathy, with a score of 7 being most discriminative. 13/54 (24%) subjects with an NDRG1 score < = 7 developed severe neuropathy, compared to only 4/57 (7%) in those with a score >7 (p = 0.017). CONCLUSION: Low NDRG1 expression in nerve tissue present within samples of surgical resection may identify subjects at risk for severe paclitaxel-induced neuropathy. Since nerve biopsies are not routinely feasible for patients undergoing chemotherapy for early breast cancer, this promising biomarker strategy is compatible with current clinical workflow.

Genomic landscapes of breast fibroepithelial tumors.

Breast fibroepithelial tumors comprise a heterogeneous spectrum of pathological entities, from benign fibroadenomas to malignant phyllodes tumors. Although MED12 mutations have been frequently found in fibroadenomas and phyllodes tumors, the landscapes of genetic alterations across the fibroepithelial tumor spectrum remain unclear. Here, by performing exome sequencing of 22 phyllodes tumors followed by targeted sequencing of 100 breast fibroepithelial tumors, we observed three distinct somatic mutation patterns. First, we frequently observed MED12 and RARA mutations in both fibroadenomas and phyllodes tumors, emphasizing the importance of these mutations in fibroepithelial tumorigenesis. Second, phyllodes tumors exhibited mutations in FLNA, SETD2 and KMT2D, suggesting a role in driving phyllodes tumor development. Third, borderline and malignant phyllodes tumors harbored additional mutations in cancer-associated genes. RARA mutations exhibited clustering in the portion of the gene encoding the ligand-binding domain, functionally suppressed RARA-mediated transcriptional activation and enhanced RARA interactions with transcriptional co-repressors. This study provides insights into the molecular pathogenesis of breast fibroepithelial tumors, with potential clinical implications.