Molecular features of luminal breast cancer defined through spatial and single-cell transcriptomics.

BACKGROUND: Intratumour heterogeneity is a hallmark of most solid tumours, including breast cancers. We applied spatial transcriptomics and single-cell RNA-sequencing on patient-derived xenografts (PDXs) to profile spatially resolved cell populations within oestrogen receptor-positive (ER+ ) breast cancer and to elucidate their importance in oestrogen-dependent tumour growth. METHODS: Two PDXs of 'ER-high' breast cancers with opposite oestrogen-mediated growth responses were investigated: oestrogen-suppressed GS3 (80-100% ER) and oestrogen-dependent SC31 (40-90% ER) models. The observation was validated via single-cell analyses on an 'ER-low' PDX, GS1 (5% ER). The results from our spatial and single-cell analyses were further supported by a public ER+ breast cancer single-cell dataset and protein-based dual immunohistochemistry (IHC) of SC31 examining important luminal cancer markers (i.e., ER, progesterone receptor and Ki67). The translational implication of our findings was assessed by clinical outcome analyses on publicly available cohorts. RESULTS: Our space-gene-function study revealed four spatially distinct compartments within ER+ breast cancers. These compartments showed functional diversity (oestrogen-responsive, proliferative, hypoxia-induced and inflammation-related). The 'proliferative' population, rather than the 'oestrogen-responsive' compartment, was crucial for oestrogen-dependent tumour growth, leading to the acquisition of luminal B-like features. The cells expressing typical oestrogen-responsive genes like PGR were not directly linked to oestrogen-dependent proliferation. Dual IHC analyses demonstrated the distinct contribution of the Ki67+ proliferative cells toward oestrogen-mediated growth and their response to a CDK4/6 inhibitor. The gene signatures derived from the proliferative, hypoxia-induced and inflammation-related compartments were significantly correlated with worse clinical outcomes, while patients with the oestrogen-responsive signature showed better prognoses, suggesting that this compartment would not be directly associated with oestrogen-dependent tumour progression. CONCLUSIONS: Our study identified the gene signature in our 'proliferative' compartment as an important determinant of luminal cancer subtypes. This 'proliferative' cell population is a causative feature of luminal B breast cancer, contributing toward its aggressive behaviours.

Pan-tumour analysis of COX-2 expression in dogs.

Cyclooxgenase-2 (COX-2) is associated with inflammatory microenvironment and tumour progression. COX-2 expression was reported in canine tumours, and anti-COX treatment showed therapeutic effects in selected tumour types. Currently, direct comparisons between different tumour types or reports were impossible due to varying evaluation protocols. Additionally, COX-2 expression in relatively uncommon tumours were yet to be evaluated. Here, we analysed COX-2 expression across various tumour types in dogs in a consistent protocol, aiming to revisit accumulated evidence in the field and report novel candidate tumours for anti-COX therapy. COX-2 expression in 32 histological types of tumours, which consisted of 347 samples in total, was investigated using immunohistochemistry followed by the Belshaw's method scoring (range: 0-12). More than the half of the samples expressed COX-2 in mast cell tumours, transitional cell carcinoma in the urinary tract, squamous cell carcinoma, liposarcoma, and melanoma, with COX-2 median scores ranging from 1-8. On the other hand, <20% tissues expressed COX-2 in the half of tumour types investigated. Overall COX-2 positive rate was 27%. In conclusion, the results confirmed COX-2 expression in the well-known COX-2-expresing tumour types and suggested novel candidate tumours for anti-COX-2 therapy. At the same time, overall COX-2 expression was low, and inter- and intra-histology heterogeneity was apparent. This study will provide a foundation reference for future research in canine tumours.

Phenotypic and molecular characterization of novel pulmonary adenocarcinoma cell lines established from a dog.

Canine pulmonary adenocarcinoma (PAC) resembles human lung tumors in never-smokers, but it is rarer than human pulmonary adenocarcinoma. Therefore, research on canine PAC is challenging. In the present study, we successfully established various novel canine PAC cell lines from a single lesion in a dog, including two parent cell lines and fourteen cloned cell lines, and characterized their cellular properties in vitro. Several of these cell lines showed epithelial-mesenchymal transition (EMT)-like and/or cancer stem cell (CSCs)-like phenotypes. We additionally assessed the sensitivity of the cells to vinorelbine in vitro. Three clonal lines, two of which showed EMT- and CSC-like phenotypes, were resistant to vinorelbine. Furthermore, we evaluated the expression and activation status of EGFR, HER2, and Ras signaling factors. The findings indicated that the cell lines we established preserved the expression and activation of these factors to varying extents. These novel canine PAC cell lines can be utilized in future research for understanding the pathogenesis and development of treatments for canine PAC.

Identification and characterization of a proliferative cell population in estrogen receptor-positive metastatic breast cancer through spatial and single-cell transcriptomics.

Background: Intratumor heterogeneity is a hallmark of most solid tumors, including breast cancers. We applied spatial transcriptomics and single-cell RNA-sequencing technologies to profile spatially resolved cell populations within estrogen receptor-positive (ER + ) metastatic breast cancers and elucidate their importance in estrogen-dependent tumor growth. Methods: Spatial transcriptomics and single-cell RNA-sequencing were performed on two patient-derived xenografts (PDXs) of "ER-high" metastatic breast cancers with opposite estrogen-mediated growth responses: estrogen-suppressed GS3 (80-100% ER) and estrogen-stimulated SC31 (30-75% ER) models. The analyses included samples treated with and without 17ß-estradiol. The findings were validated via scRNA-seq analyses on "ER-low" estrogen-accelerating PDX, GS1 (5% ER). The results from our spatial and single-cell analyses were further supported by the analysis of a publicly available single cell dataset and a protein-based dual immunohistochemical (IHC) evaluation using three important clinical markers [i.e., ER, progesterone receptor (PR), and Ki67]. The translational implication of these results was assessed by clinical outcome analyses on public breast cancer cohorts. Results: Our novel space-gene-function study revealed a "proliferative" cell population in addition to three major spatially distinct compartments within ER + metastatic breast cancers. These compartments showed functional diversity (i.e., estrogen-responsive, proliferative, hypoxia-induced, and inflammation-related). The "proliferative ( MKI67 + )" population, not "estrogen-responsive" compartment, was crucial for estrogen-dependent tumor growth, leading to the acquisition of luminal B features. The cells with induction of typical estrogen-responsive genes such as PGR were not directly linked to estrogen-dependent proliferation. Additionally, the dual IHC analyses demonstrated the distinct contribution of the Ki67 + proliferative cells toward estrogen-mediated growth and their response to palbociclib, a CDK4/6 inhibitor. The gene signatures developed from the proliferative, hypoxia-induced, and inflammation-related compartments were significantly correlated with worse clinical outcomes, while patients with the high estrogen-responsive scores showed better prognosis, confirming that the estrogen-responsive compartment would not be directly associated with estrogen-dependent tumor progression. Conclusions: For the first time, our study elucidated a "proliferative" cell population distinctly distributed in ER + metastatic breast cancers. They contribute differently toward progression of these cancers, and the gene signature in the "proliferative" compartment is an important determinant of luminal cancer subtypes.

Expression profile of immunoregulatory factors in canine tumors.

Cancers utilize a variety of molecules to escape host immune responses. Better understanding the immune environment surrounding cancer may facilitate application of innovative cancer immunotherapies, such as immune checkpoint inhibitors, to dogs as well as humans. In this study, we screened the expression of 20 immune regulatory molecules in diverse canine tumors (n = 59). Quantitative RT-PCR (qPCR) analysis revealed that some immune regulatory molecules, such as LGALS9 (coding Galectin-9) and CD48, were expressed in most canine tumors, but other molecules, such as CD274 (coding PD-L1), IL4I1, PVR, TNFSF18, ICOSLG, and TNFSF4, were rarely expressed. NECTIN2 was highly expressed in epithelial tumors but was low in non-epithelial tumors. In contrast, VSIR and CD200 expressions were low in epithelial tumors but high in non-epithelial tumors. Interestingly, several tumors expressed distinctive immunoregulatory factors. Hepatocellular carcinomas expressed FGL1, mast cell tumors expressed PDCD1LG2 (coding PD-L2), transitional cell carcinomas expressed VTCN1 (coding B7x), and lymphomas and squamous cell carcinomas expressed CD70. Consistent with qPCR results, immunofluorescence staining confirmed that hepatocellular carcinomas expressed FGL-1 protein. Thus, this study reveals the expression profile of immunoregulatory molecules in canine tumors and opens the door to better understanding the relationship between canine tumors and host immunity.

Assessment of changes in blood pancreatic lipase activities using FDC-v-LIP in dogs that underwent various surgical procedures.

OBJECTIVE: To describe the perioperative changes in blood pancreatic lipase activity and explore the contributing clinical factors associated with these changes. DESIGN: Prospective observational study. SETTING: University teaching hospital. ANIMALS: One hundred and four dogs underwent various surgical procedures under general anesthesia. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: Blood pancreatic lipase activities, which were measured using FUJI DRI-CHEM v-Lip-P (FDC-v-Lip), significantly increased postoperatively compared to preoperative measurements (premedian 58.5 U/L [range, 23-157] vs. postmedian 80 U/L [range, 22-1000], P < 0.0001). The patient with a postoperative increase in FDC-v-Lip over the normal range (35 dogs [33.6%]) had significantly higher preoperative FDC-v-Lip values. CONCLUSIONS: In this study, dogs had significantly increased pancreas-specific lipase activities after surgical procedures under general anesthesia. Direct contributors to the increase and its relevance to clinical and histological pancreatitis should be determined in the future.

Single-Cell Transcriptomics Identifies Heterogeneity of Mouse Mammary Gland Fibroblasts With Distinct Functions, Estrogen Responses, Differentiation Processes, and Crosstalks With Epithelium.

Fibroblasts have been shown to be one of the essential players for mammary gland organization. Here, we identify two major types of mouse mammary gland fibroblasts through single-cell RNA sequencing analysis: Dpp4 + fibroblasts and Dpp4 - fibroblasts. Each population exhibits unique functional characteristics as well as discrete localization in normal mouse mammary glands. Remarkably, estrogen, a crucial mediator of mammary gland organization, alters the gene expression profiles of fibroblasts in a population-specific manner, without distinct activation of estrogen receptor signaling. Further integrative analysis with the inclusion of five other publicly available datasets reveals a directional differentiation among the mammary gland fibroblast populations. Moreover, the combination with the mouse mammary epithelium atlas allows us to infer multiple potential interactions between epithelial cells and fibroblasts in mammary glands. This study provides a comprehensive view of mouse mammary gland fibroblasts at the single-cell level.

Influence of Estrogen Treatment on ESR1+ and ESR1- Cells in ER+ Breast Cancer: Insights from Single-Cell Analysis of Patient-Derived Xenograft Models.

A 100% ER positivity is not required for an endocrine therapy response. Furthermore, while estrogen typically promotes the progression of hormone-dependent breast cancer via the activation of estrogen receptor (ER)-α, estrogen-induced tumor suppression in ER+ breast cancer has been clinically observed. With the success in establishing estrogen-stimulated (SC31) and estrogen-suppressed (GS3) patient-derived xenograft (PDX) models, single-cell RNA sequencing analysis was performed to determine the impact of estrogen on ESR1+ and ESR1- tumor cells. We found that 17ß-estradiol (E2)-induced suppression of GS3 transpired through wild-type and unamplified ERα. E2 upregulated the expression of estrogen-dependent genes in both SC31 and GS3; however, E2 induced cell cycle advance in SC31, while it resulted in cell cycle arrest in GS3. Importantly, these gene expression changes occurred in both ESR1+ and ESR1- cells within the same breast tumors, demonstrating for the first time a differential effect of estrogen on ESR1- cells. E2 also upregulated a tumor-suppressor gene, IL-24, in GS3. The apoptosis gene set was upregulated and the G2M checkpoint gene set was downregulated in most IL-24+ cells after E2 treatment. In summary, estrogen affected pathologically defined ER+ tumors differently, influencing both ESR1+ and ESR1- cells. Our results also suggest IL-24 to be a potential marker of estrogen-suppressed tumors.

Optimization of SPIO Injection for Sentinel Lymph Node Dissection in a Rat Model.

The magnetic technique, consisting of a magnetic tracer and a handheld magnetometer, is a promising alternative technique for sentinel lymph node dissection (SLND) and was shown to be non-inferior to the standard technique in terms of identification rates. In this study, injection characteristics (iron dose, dilution, time course and massaging) were evaluated to optimize magnetic tracer uptake in the sentinel lymph nodes (SLN) in a rat hindleg model. 202 successful SLNDs were performed. Iron uptake in the SLN is proportional (10% utilization rate) to the injection dose between 20 and 200 µg, showing a plateau uptake of 80 µg in the SLN around 1000 µg injection. Linear regression showed that time had a higher impact than dilution, on the SLN iron uptake. Massaging showed no significant change in iron uptake. The amount of residual iron at the injection site was also proportional to the injection dose without any plateau. Time was a significant factor for wash-out of residual iron. From these results, preoperative injection may be advantageous for SLN detection as well as reduction in residual iron at the injection site by potential decrease in required injection dose.

Mammary cell gene expression atlas links epithelial cell remodeling events to breast carcinogenesis.

The female mammary epithelium undergoes reorganization during development, pregnancy, and menopause, linking higher risk with breast cancer development. To characterize these periods of complex remodeling, here we report integrated 50 K mouse and 24 K human mammary epithelial cell atlases obtained by single-cell RNA sequencing, which covers most lifetime stages. Our results indicate a putative trajectory that originates from embryonic mammary stem cells which differentiates into three epithelial lineages (basal, luminal hormone-sensing, and luminal alveolar), presumably arising from unipotent progenitors in postnatal glands. The lineage-specific genes infer cells of origin of breast cancer using The Cancer Genome Atlas data and single-cell RNA sequencing of human breast cancer, as well as the association of gland reorganization to different breast cancer subtypes. This comprehensive mammary cell gene expression atlas ( https://mouse-mammary-epithelium-integrated.cells.ucsc.edu ) presents insights into the impact of the internal and external stimuli on the mammary epithelium at an advanced resolution.

Functional characterization of androgen receptor in two patient-derived xenograft models of triple negative breast cancer.

Extensive efforts, through cell line-based models, have been made to characterize the androgen receptor (AR) signaling pathway in triple-negative breast cancer (TNBC). However, these efforts have not yet reached a consensus with regards to the mechanism of AR in TNBC. Considering that patient-derived xenografts (PDXs) are more appropriate than cell line-based models for recapitulating the structural and molecular features of a patient's tumor, we have identified and molecularly characterized two new AR-positive TNBC PDX models and assessed the impacts of AR agonist [dihydrotestosterone (DHT)] and antagonist (enzalutamide) on tumor growth and gene expression profiles by utilizing immunohistochemistry, western blots, and RNA-Seq analyses. Two PDX models, termed TN1 and TN2, were derived from two grade-3 TNBC tumors, each harboring 1∼5% of AR nuclear positive cancer cells. DHT activated AR in both PDX tumors by increasing nuclear localization and AR protein levels. However, the endpoint tumor volume of DHT-treated TN1 was 3-folds smaller than that of non-treated TN1 tumors. Conversely, the endpoint tumor volume of DHT-treated TN2 was 2-folds larger than that of non-treated TN2. Moreover, enzalutamide failed to antagonize DHT-induced tumor growth in TN2. The RNA-Seq analyses revealed that DHT mainly suppressed gene expression in TN1 (961 down-regulated genes versus 149 up-regulated genes), while DHT promoted gene expression in TN2 (673 up-regulated genes versus 192 down-regulated genes). RNA-Seq data predicted distinct TNBC molecular subtypes for TN1 and TN2. TN1 correlated to a basal-like 1 (BL1) subtype, and TN2 correlated to a basal-like 2 (BL2) subtype. These analyses suggest that TN1 and TN2, which both express functional AR, are two molecularly distinct PDX models. The molecular characterization of these PDX models expands our current knowledge on AR-positive TNBC. Our results do not support that AR is a suitable therapeutic target in TNBC. To our best knowledge, the molecular mechanisms of AR in TNBC are equivocal and should be evaluated using clinically relevant models, considering both the heterogeneous expression of AR in TNBC and the general complexities of AR signaling.

Author Correction: Aberrant expression of the COX2/PGE2 axis is induced by activation of the RAF/MEK/ERK pathway in BRAFV595E canine urothelial carcinoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Spindle assembly checkpoint competence in aneuploid canine malignant melanoma cell lines.

The spindle assembly checkpoint (SAC) is a surveillance mechanism that prevents unequal segregation of chromosomes during mitosis. Abnormalities in the SAC are associated with chromosome instability and resultant aneuploidy. This study was performed to evaluate the SAC competence in canine malignant melanoma (CMM) using four aneuploid cell lines (CMeC1, CMeC2, KMeC, and LMeC). After treatment with nocodazole, a microtubule disrupting agent, CMeC1, KMeC, and LMeC cells were arrested in M phase, whereas CMeC2 cells were not arrested, and progressed into the next cell cycle phase without cytokinesis. Chromosome spread analysis revealed a significantly increased rate of premature sister chromatid separation in CMeC2 cells. Expression of the phosphorylated form of the SAC regulator, monopolar spindle 1 (Mps1), was lower in CMeC2 cells than in the other CMM cell lines. These results indicate that the SAC is defective in CMeC2 cells, which may partially explain aneuploidy in CMM. Thus, CMeC2 cells may be useful for further studies of the SAC mechanism in CMM and in determining the relationship between SAC incompetence and aneuploidy.

PDPN Is Expressed in Various Types of Canine Tumors and Its Silencing Induces Apoptosis and Cell Cycle Arrest in Canine Malignant Melanoma.

Podoplanin (PDPN), a small transmembrane mucin-like glycoprotein, is ectopically expressed. It is also known to be linked with several aspects of tumor malignancy in some types of human tumors, including invasion, metastasis, and cancer stemness. However, there are few reports on the expression of dog PDPN (dPDPN) in canine tumors, and the association between dPDPN and tumor malignancy has not been elucidated. We identified that 11 out of 18 types of canine tumors expressed dPDPN. Furthermore, 80% of canine malignant melanoma (MM), squamous cell carcinoma, and meningioma expressed dPDPN. Moreover, the expression density of dPDPN was positively associated with the expression of the Ki67 proliferation marker. The silencing of dPDPN by siRNAs resulted in the suppression of cell migration, invasion, stem cell-like characteristics, and cell viability in canine MM cell lines. The suppression of cell viability was caused by the induction of apoptosis and G2/M phase cell cycle arrest. Overall, this study demonstrates that dPDPN is expressed in various types of canine tumors and that dPDPN silencing suppresses cell viability through apoptosis and cell cycle arrest, thus providing a novel biological role for PDPN in tumor progression.

Aberrant expression of the COX2/PGE2 axis is induced by activation of the RAF/MEK/ERK pathway in BRAFV595E canine urothelial carcinoma.

Cancer-promoting inflammation is an important event in cancer development. Canine urothelial carcinoma (cUC) overexpresses prostaglandin E2 (PGE2) and has a unique sensitivity to cyclooxygenase 2 (COX2)-inhibiting therapy. In addition, majority of cUC harbour BRAFV595E mutation. However, mechanisms underlying aberrant PGE2 production in BRAFV595E cUC patients remain unclear. Drug screening revealed that inhibition of RAF/MEK/ERK pathway, p38 and JNK pathway reduced PGE2 production in cUC cells. By pharmacological inhibition of the multiple components in the pathway, activation of the ERK MAPK pathway was shown to mediate overexpression of COX2 and production of PGE2 in BRAFV595E cUC cells. In silico gain-of-function analysis of the BRAF mutation also implicated involvement of mutation in the process. The positive association between ERK activation and COX2 expression was further validated in the clinical patients. Moreover, it was also suggested that p38 and JNK regulates PGE2 production independently of ERK pathway, possibly through COX2-dependent and COX1-/COX2- independent manner, respectively. In conclusion, this study demonstrated that activation of ERK induces production of PGE2 in BRAFV595E cUC cells, which is also independently regulated by p38 and JNK. With its unique vulnerability to COX-targeted therapy, BRAFV595E cUC may serve as a valuable model to study the tumour-promoting inflammation.

Overexpression of human epidermal growth factor receptor 2 in canine primary lung cancer.

Human epidermal growth factor receptor 2 (HER2) overexpression has been reported in various human cancers. HER2-targeted therapies showed clinical responses in humans with HER2-positive tumors. The incidence of canine primary lung cancer (cPLC) is increasing, but there are no effective systemic therapies for dogs with late-stage cPLC. HER2-targeted therapy could be an option for cPLC, but HER2 expression in cPLC remains unknown. We evaluated HER2 expression in cPLC. Immunohistochemical analysis revealed that 3 samples (19%) scored 3+; 8 (50%), 2+; 5 (31%); and 1+ and 0 (0%), 0. Of the cPLC tissues, 69% were HER2 positive (scored ≥2+). These data would lead to further evaluation of the role of HER2 in cPLC as a mechanism of malignancy and therapeutic target.

Environmental Carcinogenesis at the Single-Cell Level.

Elucidating the mechanisms behind how exposure to environmental chemicals can lead to cancer is not easy due to the complex natures of these compounds and the challenges to establish biologically relevant experimental models to study them. Environmental chemicals often present selective mechanisms of action on different cell types and can be involved in the modulation of targeted cells and their microenvironment, including immune cells. Currently, the limitations of traditional epidemiologic correlation analyses, in vitro cell-based assays, and animal models are that they are unable to comprehensively examine cellular heterogeneity and the tissue-selective influences. To this end, we propose utilizing single-cell RNA-sequencing (scRNA-seq) to more effectively capture the subtle and complex effects of environmental chemicals and how their exposure could lead to cancer. scRNA-seq's capabilities for studying gene expression level data at a significantly higher resolution relative to bulk RNA-sequencing (RNA-seq) enable studies to evaluate how environmental chemicals regulate gene transcription on different cell types as well as how these compounds impact signaling pathways and interactions between cells in the tissue microenvironment. These studies will be valuable for evaluating environmental chemicals' carcinogenic properties at the individual cell level.See all articles in this CEBP Focus section, "Environmental Carcinogenesis: Pathways to Prevention."

ErbB2 Copy Number Aberration in Canine Urothelial Carcinoma Detected by a Digital Polymerase Chain Reaction Assay.

Urothelial carcinoma (UC) is the most common tumor affecting the urinary bladder of dogs. Protein overexpression of ErbB2 (the canine homolog of HER2) has been observed in dogs with UC. However, no study regarding ErbB2 copy number aberration (CNA) is reported in dogs with UC. In this study, a digital PCR assay for detecting CNA of canine ErbB2 was developed. DNA samples were isolated from 83 formalin-fixed, paraffin-embedded urinary bladder tissues (36 UC, 8 polypoid cystitis, and 39 normal) and 94 urinary sediments (54 UC, 30 nonneoplastic, and 10 normal). The copy number of canine chromosome 8 (CFA8) was used as a control. In the urinary bladder tissues, ErbB2 CNA was detected in 12 of 36 (33%) UC, 2 of 8 (25%) polypoid cystitis, and 0 of 39 (0%) normal controls. In the urinary sediments, ErbB2 CNA was also detected in 19 of 54 (35%) UC; however, no ErbB2 CNA was detected in nonneoplastic diseases or normal controls. The sensitivity and specificity of ErbB2 CNA in urinary sediment for the detection of UC were 35% and 100%, respectively. There was a positive correlation between the copy number ratios of ErbB2 to CFA8 in the urinary bladder tissues and urinary sediments. Our findings indicate that the digital PCR assay of urinary sediments may be a useful, noninvasive method for detecting ErbB2 CNA in dogs with UC.

Single-cell RNA-sequencing analysis of estrogen- and endocrine-disrupting chemical-induced reorganization of mouse mammary gland.

Menopause is a critical window of susceptibility for its sensitivity to endocrine disrupting chemicals due to the decline of endogenous estrogen. Using a surgical menopausal (ovariectomized) mouse model, we assessed how mammary tissue was affected by both 17ß-estradiol (E2) and polybrominated diphenyl ethers (PBDEs). As flame retardants in household products, PBDEs are widely detected in human serum. During physiologically-relevant exposure to E2, PBDEs enhanced E2-mediated regrowth of mammary glands with terminal end bud-like structures. Analysis of mammary gland RNA revealed that PBDEs both augmented E2-facilitated gene expression and modulated immune regulation. Through single-cell RNA sequencing (scRNAseq) analysis, E2 was found to induce Pgr expression in both Esr1+ and Esr1- luminal epithelial cells and Ccl2 expression in Esr1+ fibroblasts. PBDEs promote the E2-AREG-EGFR-M2 macrophage pathway. Our findings support that E2 + PBDE increases the risk of developing breast cancer through the expansion of estrogen-responsive luminal epithelial cells and immune modulation.