An adaptive Epithelial-Mesenchymal Transition Program Enables Basal Epithelial Cells to Bypass Stress-Induced Stasis and Contributes to Metaplastic Breast Cancer Progenitor State

An adaptive Epithelial-Mesenchymal Transition Program Enables Basal Epithelial Cells to Bypass Stress-Induced Stasis and Contributes to Metaplastic Breast Cancer Progenitor State | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article An adaptive Epithelial-Mesenchymal Transition Program Enables Basal Epithelial Cells to Bypass Stress-Induced Stasis and Contributes to Metaplastic Breast Cancer Progenitor State Joseph A. Caruso, Thea D. Tlsty This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4980285/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 18 Dec, 2024 Read the published version in Breast Cancer Research → Version 1 posted 10 You are reading this latest preprint version Abstract Background: Human mammary epithelial cell (HMEC) cultures encounter a stress-associated barrier termed stasis, during which most cells adopt a senescence-like phenotype. From these cultures, rare variants emerge from the basal epithelial population, re-initiating growth. Variants exhibit pre-malignant properties, including an aberrant epigenetic program that enables continued proliferation and acquisition of genetic changes. Following oncogenic transformation, variants produce tumors that recapitulate the histopathological characteristics of metaplastic breast cancer (MBC), a rare subtype characterized by squamous and mesenchymal differentiation. Methods: Using the conventional serum-free HMEC culture system, we probed the capacity for phenotypic plasticity inherent to basal epithelial cell populations from human breast tissue as they navigated stasis and emerged as variant populations. Results: We observed robust activation of a TGF-β-dependent epithelial-mesenchymal transition (EMT) program in basal epithelial cells during stasis, followed by subsequent attenuation of this program in emerging variants. Inhibiting the TGF-β pathway or depleting the EMT regulators Snail or Slug allowed basal epithelial cells to collectively bypass stasis, demonstrating that cellular dysfunction and arrest resulting from TGF-β and EMT activation are central to this in vitro barrier. The spontaneous emergence of variants from stasis cultures was associated with a restricted EMT trajectory, which diverted cells away from a complete mesenchymal state characterized by irreversible growth arrest, and instead limited variants to epithelial and intermediate EMT states associated with greater proliferative capacity and stemness. Epigenetic mechanisms, which contributed to the dysregulated growth control characteristic of the variant phenotype, also contributed to the constrained EMT program in variants. By overcoming the cellular dysfunction and growth arrest resulting from TGF-β and EMT activation, variants exhibited increased oncogenic transformation efficiency compared to pre-stasis basal epithelial cells. Inhibiting the TGF-β pathway prior to stasis significantly reduced EMT in the basal epithelial population, alleviated selective pressure driving variant emergence, and enhanced oncogenic transformation efficiency, resulting in tumors with markedly diminished metaplastic differentiation. Conclusions: This study reveals how adaptive EMT reprogramming governs basal epithelial cell fate decisions and contributes to the development of MBC progenitors by restricting access to terminal mesenchymal states that induce growth arrest and, instead, favoring intermediate states with enhanced tumorigenic potential. Human mammary epithelial cells Basal Myoepithelial TGF-β Pathway Epithelial-mesenchymal transition Epigenetic Regulation Metaplastic Breast Cancer Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Background The milk-producing lobules and ductal networks within a healthy human breast consist of an inner luminal layer of polarized columnar epithelial cells and an outer basal layer of contractile myoepithelial cells. Both luminal and basal epithelia contain hierarchies of progenitor cells [ 1 – 6 ], whose transformation significantly contributes to the histological and molecular heterogeneity observed among breast cancers [ 7 – 12 ]. These cancers are commonly classified into subtypes that can express markers of luminal or basal differentiation states [ 13 – 15 ]. However, most, including those with characteristically basal expression patterns, are believed to originate from progenitor cell populations in the luminal epithelium [ 7 , 11 , 16 ]. Only rare metaplastic breast cancers (MBCs), representing 0.2-5 percent of all breast cancers [ 17 ], are proposed to originate from basal progenitor populations [ 9 ] primarily responsible for the maintenance of the myoepithelial layer [ 1 , 2 ]. MBC is a histologically defined classification typified by differentiation into cell types that are not normally found in breast tissue, including squamous and mesenchymal (e.g., spindle cell, chondroid, and osseous metaplasia) phenotypes generating foci of bone, muscle and cartilage cells. These cancers predominantly fall within claudin-low or basal-like intrinsic molecular subtypes and are generally clinically classified as triple-negative breast cancers (TNBCs) [ 17 , 18 ]. Epithelial-mesenchymal transition (EMT) is closely associated with the pathogenesis of MBC and is believed to contribute to its comparatively poor prognosis, chemoresistance, and high metastatic potential. EMT is a developmental process in which epithelial cells acquire mesenchymal characteristics, including increased motility and invasiveness, which can be reactivated in response to injurious or pathological conditions, including carcinogenesis [ 19 , 20 ]. Notably, EMT is not a binary switch from the epithelial to mesenchymal state. Instead, it produces a spectrum of intermediate cellular states that exhibit both epithelial and mesenchymal characteristics. These intermediate states are associated with an increase in stem cell properties, tumorigenicity, drug resistance, and metastasis [ 21 – 28 ]. The human mammary epithelial cell (HMEC) culture system, developed more than 40 years ago [ 29 ] by the Stampfer group, remains one of the most widely used platforms for studying epithelial cell biology and carcinogenesis. Establishing HMEC cultures requires mincing healthy human breast tissue, enzymatically digesting matrix components, and removing stromal cell types [ 29 , 30 ]. In healthy tissues, cellular cooperation relies on complex communication networks involving a vast array of biochemical and biophysical interactions, including interactions between parenchymal and stromal cells, the extracellular matrix, and the systemic milieu [ 31 ]. Constantly integrating information from these networks defines a cell's operational state and sets limits on cellular behavior [ 32 ]. Indeed, after 5–20 population doublings in culture, HMECs enter stasis, a 'stress-associated' barrier that is dependent on the engagement of the tumor suppressor retinoblastoma (Rb) and characterized by most cells adopting a senescence-like phenotype [ 8 , 33 , 34 ]. Variants, exclusively derived from the basal compartment [ 9 ], evade this fate and initiate a second phase of logarithmic growth as clonal outgrowths from stasis cultures. These variant basal epithelial cells (vHMECs) continue to experience telomere erosion, allowing for an additional 20–70 population doublings before reaching a telomeric crisis state known as agonescence [ 29 , 34 ]. Because of their extended lifespans, variants are more broadly available from commercial sources and are often marketed simply as HMECs. These cells have been immortalized and oncogenically transformed into cell lines that are widely used in laboratories worldwide to model breast epithelial biology and carcinogenesis. Thus, vHMECs have profoundly influenced many fields, including the study of stromal-epithelial interactions, aging, EMT, and stem cell biology. Variant basal epithelial outgrowths that emerge from stasis and dominate post-selection HMEC cultures exhibit compromised growth control and acquire other characteristics associated with malignant phenotypes. Many of these alterations result from the activation of an epigenetic program, which is surprisingly consistent between donors [ 35 – 41 ]. Characteristic epigenetic alterations include hypermethylation of specific genes and loci targeted by polycomb repressor complex 2 (PRC2), including the promoters of genes regulating growth control and differentiation, such as cyclin-dependent kinase inhibitor 2A (CDKN2A) and homeobox A9 [ 35 , 37 , 41 ]. In addition to alterations in the chromatin landscape, cytogenetic analysis reveals an accumulation of chromosomal abnormalities in variant cultures, including translocations, deletions, telomeric associations, polyploidy, and aneuploidy [ 34 ]. Collectively, modifications in the variant HMEC state are sufficient to allow these cells to tolerate oncogenic events, such as the introduction of oncogenic forms of RAS, without undergoing oncogene-induced senescence [ 39 ]. Additional factors, such as the viral oncoprotein SV40 large T antigen and telomerase reverse transcriptase (TERT), are required to fully immortalize these cells, which is a prerequisite for the introduction of a transforming oncogene to achieve tumorigenicity in immunocompromised mice [ 42 ]. Notably, the malignant transformation of variant cultures characteristically results in the development of tumors that recapitulate the histopathological characteristics of the rare MBC subtype. Thus, many widely used breast cancer models derived from HMEC cultures, such as MCF-10AT, HMT-3522 T4-2, and HMLER, form tumors exhibiting metaplastic squamous and spindle cell differentiation [ 8 , 42 , 43 ]. In this study, we observed robust activation of the TGF-β pathway and EMT as HMEC cultures entered stasis, which was followed by the attenuation of these programs in emerging variant basal epithelial populations. In this system, activation of the TGF-β pathway and EMT was associated with cellular dysfunction, driving cells toward a senescent-like phenotype. Rare variants that spontaneously emerged from stasis cultures exhibited a restricted EMT trajectory, which prevented their complete transition to a mesenchymal state characterized by P-cadherin loss and irreversible growth arrest, instead confining them to dynamic transitions between epithelial and intermediate EMT states that maintained P-cadherin expression. The aberrant epigenetic program activated in variants not only compromised growth control and differentiation [ 35 – 41 ] but also contributed to the restricted EMT program observed in variants. The variant phenotype was associated with increased susceptibility to oncogenic transformation compared to unselected, pre-stasis basal epithelial cells and produced tumors with metaplastic elements. Pharmacological inhibition of the TGF-β pathway in cultured basal epithelial cells limited EMT program activation, reduced the selective pressures underlying variant emergence, and enhanced tumorigenic potential in pre-stasis basal epithelial populations, which produced tumors with dramatically reduced metaplastic differentiation. Our findings suggest that in a perturbed microenvironment, the adaptability of the EMT program, favoring intermediate states over complete mesenchymal differentiation, is crucial for directing basal epithelial cells toward a progenitor state that can progress to MBC. Methods Reagents were obtained from Thermo Fisher Scientific unless otherwise stated. Tissue Dissociation Disease-free breast tissue from women undergoing reduction mammoplasty was provided by the Cooperative Human Tissue Network (CHTN) Western Division, Nashville, TN, and Kaiser Foundation Research Institute, Oakland, CA. Human breast tissue was minced and enzymatically dissociated in RPMI with 2.1 mM L-glutamine and 25 mM HEPES supplemented with 10% (v/v) FBS (Atlanta Biologicals), 10 U/mL penicillin, 10 µg/mL streptomycin, 2.5 µg/mL Amphotericin B, 50 µg/mL gentamycin, Collagenase Type-2 (Worthington), and 100 U/mL hyaluronidase (Sigma-Aldrich) at 37°C for 12–16 hours. The digest was centrifuged at 300 g for 10 minutes and washed with RPMI-1640 supplemented with 10% FBS. Epithelial-enriched structures (referred to as tissue organoids) were recovered by filtration through a 150 µm nylon mesh filter and frozen for long-term storage in Dulbecco's Modified Eagle Medium (DMEM) containing 50% (v/v) FBS and 10% (v/v) DMSO. The age and race of the donors are listed in Table S2 . Fluorescence-Activated Cell Sorting (FACS) and Flow Cytometry : Tissue organoids were dissociated into single cells using 5 U/mL Dispase (Stem Cell Technologies) containing 1 µg/mL DNAse I (Stem Cell Technologies) for 5–10 minutes and then 0.05% trypsin with 0.14 g/L Ethylenediaminetetraacetic Acid (EDTA) for 5–10 minutes. Cultured cells were collected from the tissue culture plates by incubating them for 10–15 minutes with TrypLE. The cells were washed with Phosphate-Buffered Saline (PBS) + 1% Bovine Serum Albumin (BSA), filtered through a 40-µm cell strainer (Falcon), resuspended at a concentration of 1 × 10⁶ cells in 100 µL of PBS + 1% BSA, and stained for 15 minutes at room temperature with the appropriate antibody combination. When necessary, the cells were washed with PBS + 1% BSA and stained for 15 minutes at room temperature with 1 µL streptavidin conjugated with BV785 (BioLegend) and 0.5 µg/mL DAPI (Sigma). The cells were washed, resuspended at 1 × 10⁶ cells in 100 µL PBS + 1% BSA, filtered through a 40-µm cell strainer, and loaded onto a FACS Aria II Cell Sorter (BD) fitted with a 100 µm nozzle. To isolate basal epithelial cells from tissue organoids, we immunostained cells with 5 µL of anti-CD10 mouse monoclonal antibody (mAb) (clone: HI10a) conjugated with PE (Biolegend), 1.5 µL of anti-Epcam mouse mAb (clone VU1D9) conjugated with FITC (Stem Cell Technologies), 4 µL of anti-CD49f rat mAb (clone GoH3) conjugated with APC (Biolegend), 8 µL of anti-CD2 mouse mAb (clone RPA-2.10) conjugated with Biotin (Becton Dickinson [BD]), 8 µL of anti-CD3 mouse mAb (clone HIT3a) conjugated with Biotin (BD), 8 µL of anti-CD16 mouse mAb (clone 3G8) conjugated with Biotin (BD), 8 µL of anti-CD64 mouse mAb (clone 10.1) conjugated with Biotin (BD), 4 µL of anti-CD31 mouse mAb (clone MBC78.2) conjugated with Biotin (ThermoFisher Scientific), and 1 µL of anti-CD45 mouse mAb (clone HI30) conjugated with Biotin (Biolegend). To isolate EMT phenotypes from cultured basal epithelial cells, we used 15 µL of anti-P-cadherin mouse mAb (clone CSTEM29) conjugated with APC (Thermo Fisher Scientific) and 4 µL of anti-N-cadherin mouse mAb (clone 8C11) conjugated with PE (BioLegend). After gating out debris, doublets, dead cells (DAPI+), and, if applicable, unwanted cell types (BV785+), the remaining target populations expressing the desired immunoprofiles were collected into a 1.5 mL tube containing 0.5 mL cell culture medium. The HMECs were pulsed with 10 µM EdU for one hour before the proliferation index was determined. The Click-iT Plus EdU Alexa Fluor 488 Flow Cytometry Assay Kit and FxCycle PI/RNase Staining Solution (Thermo Fisher Scientific) were used according to the manufacturer’s instructions. The cells were analyzed using an LSRFortessa flow cytometer (BD Biosciences). Cell Culture Dissociated human breast tissue organoids were plated in T-75 flasks (Corning) in Mammary Epithelial Cell Growth Medium (MEGM, Lonza), which includes hydrocortisone, human epidermal growth factor, insulin, bovine pituitary extract (BPE), gentamicin, and amphotericin B. After the initial plating, differential trypsinization was performed to selectively detach and aspirate unwanted fibroblasts, followed by collection of epithelial cells for further passaging. Sorted populations were cultured in one well of a six-well plate, transferred to a p100 plate, and then to a T-75 flask. Subculturing was performed at 70–80% confluence using TrypLE Express for cell dissociation. Cells were incubated at 37°C in a humidified atmosphere containing 5% CO₂. The medium was changed every two to three days. Cell viability and counts were assessed using trypan blue exclusion and a Countess automated cell counter (Thermo Fisher Scientific) before reseeding or experimental use. The number of population doubling (PD) per passage was determined using the equation PD = log [A/B]/log2, where A is the number of collected cells and B is the number of plated cells. Attachment efficiency was consistently greater than 95% during the exponential growth phase. Cell populations were frozen in cryopreservation medium composed of 90% FBS and 10% DMSO (Sigma-Aldrich), aliquoted into cryovials, and gradually cooled to -80°C before being transferred to liquid nitrogen for long-term storage. The TGF-β inhibitor A83-01 (Tocris) was dissolved in DMSO at 1 mM and used at a concentration of 0.5 µM; an equivalent amount of DMSO was added to control cultures. Mammosphere Culture Single cells were directly sorted into 24-well ultra-low attachment plates (Corning) at a density of 10,000 cells per well in serum-free mammary epithelial basal medium (MEBM) (Lonza) supplemented with B27 (Thermo Fisher Scientific), 20 ng/mL Epidermal Growth Factor (EGF, Peprotech), 20 ng/mL Fibroblast Growth Factor 2 (FGF2, Peprotech), and 4 µg/mL heparin (Stem Cell Technologies). Mammospheres were counted and collected by centrifugation at 100 × g after 10 days of culture, dissociated enzymatically for 5–10 minutes in TrypLE Express, filtered through a 40-µm cell strainer, and cultured again in suspension [ 44 ]. Lentiviral Transduction : Basal epithelial cell populations, 1 × 10⁶ in a single well of a six-well plate, were incubated with concentrated lentiviral particles at a multiplicity of infection of three with 8 µg/mL polybrene overnight for 12–16 hours. For transformation, the cells were infected sequentially with (1) TERT with a hygromycin selection marker (Genecopoeia), (2) SV40 small T and large T antigens with a puromycin selection marker (Genecopoeia), and (3) KRAS G12V with a blasticidin selection marker (Amsbio). Cells were selected using 200 µg/mL hygromycin (Sigma-Aldrich), 1 µg/mL puromycin (InvivoGen), and 10 µg/mL blasticidin (InvivoGen). Transformed cells were grown in DMEM/F12 (1:1) supplemented with 5% Newborn Calf Serum (Hyclone), 10 µg/mL insulin, 10 ng/mL EGF, and 1 µg/mL hydrocortisone. CRISPR-mediated knockout of core EMT transcription factors, Snail (SNAI1) and Slug (SNAI2), Zeb (ZEB1), and Twist (TWIST1) was carried out using transEDIT CRISPR gRNA target gene sets from Transomics. Three constructs were tested for each gene. To suppress the EMT program, we chose to continue with: SNAI1 TEVH-1091413-pCLIP-All-EFS-ZsGreen, SNAI2 TEVH-1091323-pCLIP-All-EFS-ZsGreen, and non-targeting control TELA1013-CRISPR-NT#1-pCLIP-All-EFS-ZsGreen. CRISPR-mediated knockout of EZH2 was carried out using the Edit-R Human hEF1a-EGFP All-in-one Lentiviral sgRNA construct (Horizon Discovery Biosciences): EZH2 VSGH12180-256507490 and non-targeting VSGC12063. Soft-Agar Assay A solution of 1% agar (Sigma-Aldrich) in PBS was prepared, and 1.5 mL was added to each well of a six-well plate and allowed to gel at room temperature for 30 minutes as the bottom layer. Cells were prepared by suspending 1 × 10⁵ cells in 3 mL pre-warmed DMEM/F12 containing 20% (v/v) FBS, 10 U/mL penicillin, and 10 µg/mL streptomycin. The cell suspension was gently mixed with 2 mL of 1% agar in PBS to achieve a cell density of 20,000 cells/mL in 0.4% agar; 1 mL of this mixture was added to each well of the six-well plate and kept at 4°C for 10 minutes to allow quick gelling, followed by adding 0.5 mL culture medium on top of the agar gel. The cells were incubated at 37°C with 5% CO₂, and the medium was changed every four days for four weeks. Colonies were stained with 0.5 mL of 0.005% Crystal Violet (Sigma-Aldrich) in PBS containing 4% formaldehyde for one hour and then rinsed with distilled water before counting the visible colonies. Xenograft Study Exponentially growing cultures of oncogenically transformed cells (4 × 10⁶) were injected subcutaneously into 8–12-week-old NSG female mice in 50% Matrigel. Tumors were allowed to grow for 3–6 months. After the study, xenograft tumors were weighed, formalin-fixed, and paraffin-embedded at the UCSF Histology & Biomarker Core. The maximum tumor burden was not reached in any of the mice in this study. Sandwich ELISA Analysis of conditioned media from cultured HMECs from three donors at P2 (pre-stasis), P4 (stasis), and P8 (post-stasis) was initially performed using the Raybiotech Quantibody Human Kiloplex array, which detects 1,000 human biomarkers, as a service by the manufacturer. We followed up on interesting candidates using sandwich ELISA kits from Raybiotech to detect TGF-β, Activin A, Serpin E1, Wnt4, R-Spondin2, DKK1, and Follistatin according to the manufacturer’s protocols. For these assays, conditioned media were generated by culturing 300,000 cells per well in a six-well plate in 1 mL of serum-free basal media for 24 h. Each analyte was measured in eight biological and three technical replicates. The concentration of each factor in the conditioned media was extrapolated from the corresponding standard curves. Immunofluorescence : Cells were cultured on four-well glass chamber slides (Millipore). The cells were fixed with 4% paraformaldehyde in PBS for 15 minutes at room temperature and permeabilized for 10 minutes in PBS containing 0.25% Triton X-100. Non-specific binding was blocked with 10% donkey serum (Jackson Immunoresearch), 0.3M glycine, and 0.1% Tween-20 in PBS for one hour. Primary antibodies were applied overnight in 1% donkey serum and 0.1% Tween-20 in PBS. The following primary antibodies were used: p16 (E6H4, Roche) mouse mAb (pre-diluted), p21 Waf1/Cip1 (12D1, Cell Signaling Technologies) rabbit mAb (1:400), p63-α (D2K8X, Cell Signaling Technologies) XP rabbit mAb (1:200), α-tubulin (DM1A, Sigma-Aldrich) mouse mAb (1:4000), MUC1 (HMFG2, Abcam) mouse mAb (1:200), Lamin A + Lamin C (EPR4100, Abcam) rabbit mAb (1:500), cytokeratin 19 (EPR1579Y, Abcam), rabbit mAb (1:400), and cytokeratin 14 (LL002, Abcam) rabbit mAb (1:200). Alexa Fluor 555-conjugated phalloidin (1:400) was used to stain F-actin (Thermo Fisher Scientific). Secondary donkey anti-mouse 488 and anti-rabbit 555 antibodies (diluted 1:500) were added to 1% donkey serum and 0.1% Tween-20 in PBS for one hour. DAPI (0.2 µg/mL in PBS) was added for five minutes at room temperature. The slides were washed three times for three minutes between each step in PBS containing 0.1% Tween-20. The coverslips were mounted with Vectashield HardSet Mounting Medium (Vector Laboratories). Images were acquired using a Leica SP8 laser scanning confocal microscope (only Tubulin/F-actin and Lamin A/C) at a resolution of 1,024 × 1,024, processed using LASX software (Leica Microsystems) or a BZ-X800 fluorescence microscope (Keyence), and processed using ImageJ (version 2.14.0/1.54f). Quantitative Polymerase Chain Reaction (qPCR) : Total RNA was isolated from cells and treated with DNase I using the RNeasy Mini Kit (Qiagen). cDNA synthesis was performed using the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific). Quantitative PCR was performed on a CFX-96 (Bio-Rad Laboratories) thermocycler using 2x SsoFast Master Mix (Bio-Rad Laboratories) and analyzed using the standard curve method. A standard curve was prepared using the cDNA produced from Human Reference RNA (Agilent Technologies). Pre-made TaqMan Gene Expression Assays (Thermo Fisher Scientific) were used: ZEB1: Hs00232783_m1, SNAI1: Hs00195591_m1, SNAI2: Hs00161904_m1, TWIST1: Hs01675818_s1, CDH3: Hs00999915_m1, VIM: Hs00958111_m1, CDH1: Hs01023895_m1, CDH2: Hs00983056_m1, CDKN1A: Hs00355782_m1, MMP-2: Hs01548727_m1, FN1: Hs01549976_m1, CHRDL2: Hs01060234_m1, GREM1: Hs01879841_s1, and WNT5a: Hs00998537_m1. Custom primer-probe sets were used for GUSB: Forward Primer: CTCATTTGGAATTTTGCCGATT, Reverse Primer: CCGAGGAAGATCCCCTTTTTA, and Probe: FAM-TGAACAGTCACCGACGAGAGTGCTGGTA-TAM and CDKN2A (specific for p16 INK4A and not p15 ARF ): Forward Primer: CCAACGCACCGAATAGTTACG, Reverse Primer: GAGTGGCGGAGCTGCT, and Probe: FAM-CCGATCCAGGTCATGATG-TAM produced by Integrated DNA Technologies. GUSB expression was used to normalize variance in the input cDNA. Western Blotting : Cells were washed with PBS and lysed using radioimmunoprecipitation assay (RIPA) buffer containing the HALT protease and phosphatase inhibitor cocktail. For each sample, a cell extract (corresponding to 50 µg of protein as determined by a Micro BCA Protein Assay Kit) was prepared in NuPAGE LDS Sample Buffer (4X) with Reducing Agent (10X), heated at 95°C for 15 minutes, and electrophoresed in each lane of a NuPAGE 4–12% Bis-Tris gel along with a BenchMark Pre-stained Protein Ladder (Thermo Fisher Scientific). Proteins were transferred onto a nitrocellulose membrane (Biorad) overnight at 4°C at 36 mV. The membrane was washed in TBS-T (0.05 M Tris-HCl pH 7.5, 0.15 M NaCl, 0.05% Tween-20), blocked for one hour at room temperature in blotto (5% nonfat dry milk in TBS-T), incubated with primary antibodies overnight at 4°C, washed in TBS-T, incubated with goat anti-rabbit horseradish peroxidase conjugate secondary antibody (Jackson Immunoresearch) 1:5000 in Blotto for one hour, washed in TBS-T, and developed with ECL Western Blotting Substrate. Chemiluminescent signals were detected by exposing the CL-XPosure Film to the membranes or using the KwikQuant Digital Western Blot Detection System. Primary antibodies used at a 1:1000 dilution from Cell Signaling Technologies were Slug (C19G7) rabbit mAb, ZEB1 (D80D3) rabbit mAb, Slug (C19G7) rabbit mAb, E-cadherin (24E10) rabbit mAb, N-cadherin (D4R1H) XP rabbit mAb, Vimentin (D21H3) XP rabbit mAb, RAS rabbit pAb, EZH2 (D2C9) XP rabbit mAb, Histone H3 (D1H2) XP rabbit mAb, and β-Actin (13E5) rabbit mAb; and from Active Motif Histone H3K27me3 (MABI 0323) mouse mAb. Full images of the western blots are shown in Supplemental Figs. S4 and S5. Multiplex Immunohistochemistry : Sections (5 µm thick) were cut from FFPE tissue blocks and placed on positively charged Superfrost Plus microscopy slides. Slides were baked at 60°C overnight, deparaffinized in xylene, and rehydrated in graded ethanol (100%, 100%, 95%, 85%, and 70%) in distilled H₂O. Endogenous peroxidases were quenched with 3% H₂O₂ (Sigma-Aldrich) diluted in phosphate-buffered saline (PBS). Heat-induced antigen retrieval was performed in citrate buffer pH 6.0 (Sigma-Aldrich) at 95°C for 15 minutes. Non-specific antibody binding was blocked using a Background Sniper (Biocare Medical). The tissue sections were incubated for one hour at room temperature (RT) with each primary antibody in 1% BSA and 30 minutes in pre-diluted MACH 2 conjugated anti-mouse or anti-rabbit secondary antibodies (Biocare Medical). The slides were washed in TNT buffer (0.1 M TRIS-HCL pH 7.5, 0.15 M NaCl, and 0.05% Tween-20) following blocking and three times for five minutes each after applying both primary and secondary antibodies. The signal was developed using a Tyramide Signal Amplification (TSA) solution (Akoya): FITC (two minutes), Cy3 (three minutes), or Cy5 (seven minutes). The antibody complex was removed by heating at 95°C in citrate buffer pH 6.0 for five minutes to allow for multiplex staining. Nuclei were counterstained with 3 µM DAPI in PBS for five minutes, washed in distilled H₂O, and mounted using Vectashield HardSet Mounting Medium (Vector Laboratories). The primary antibodies used from Cell Signaling Technologies were Slug (C19G7) rabbit mAb (1:1000), Snail (C15D3) rabbit mAb (1:1000), p63-α (D2K8X) XP rabbit mAb (1:2000), and vimentin (D21H3) rabbit mAb (1:5000); from AbCam Cytokeratin 5 (SP27) rabbit mAb (1:4000), Lamin A + Lamin C (EPR4100) rabbit mAb (1:10000), Cytokeratin 13 (EPR3671) rabbit mAb (1:2000), Ki67 (SP6) rabbit mAb (1:6000), and Cytokeratin 14 (SP53) rabbit mAb (1:3000).The slides were imaged using a BZ-X800 analyzer. Images were prepared for analysis using the ImageJ software (version 2.14.0/1.54f). Analysis was performed using the Qupath software (v0.5.0)[ 45 ]. Using the pre-trained models, nuclear segmentation was performed using StarDist[ 46 ], normalizePercentiles = 1, 99, threshold = 0.5, pixel size = 0.5, and cell expansion = 5.0. Individual classifiers were trained for each marker and combined to create a composite classifier for scoring the cells within the annotated regions. Statistical Analysis For each graph, the mean and standard deviation are shown in red. Unless otherwise noted in the figure legend, each black circle indicates a single donor. Pairwise comparisons were performed using an unpaired Welch’s t-test (GraphPad Prism Software). Levels of significance used were *0.05 to 0.01, **0.01 to 0.001, ***0.001 to 0.0001, and ***<0.0001. Results Robust Activation of the TGF-β Pathway and Induction of EMT in Stasis HMECs Precede Variant Emergence When cultured in serum-free Mammary Epithelial Cell Growth Medium (MEGM), HMECs isolated from three tissue donors experienced growth cessation after 5–10 population doublings. Stasis lasted between 20 and 40 days before the emergence of clonal variant populations, which re-initiated logarithmic growth (Fig. 1 A). Nearly complete proliferative arrest was observed at stasis (Fig. 1 B), marked by the upregulation of p16 and development of a senescence-like morphology (Fig. 1 C). These results are consistent with previously published studies that have shown the clonal emergence of rare variants with pre-malignant properties within a field of growth-arrested, senescent-like epithelial cells [ 29 , 33 , 34 , 41 , 47 , 48 ]. Previous studies of this culture system have primarily focused on comparing the initial pre-selection HMEC population to post-selection variants. In contrast, we chose to characterize stasis cultures, specifically focusing on the factors that enforce growth arrest. Conditioned media (CM) from stasis cultures, diluted 1:1 with fresh MEGM, induced significant growth arrest in early passage HMECs compared to fresh MEGM or CM from pre-stasis HMECs at passage 2 (P2), mixed 1:1 with fresh MEGM (Fig. 1 D). Senescent and DNA double-strand breaks are associated with complex secretory programs, termed senescence-associated secretory phenotype (SASP) and stress-elicited extrinsic phenotype (SEEP), respectively, which can initiate or reinforce growth arrest in both autocrine and paracrine manners [ 49 – 51 ]. To characterize alterations in the cellular secretory program during HMEC culture, we collected CM and extracted cellular RNA and protein from cultures of pre-stasis HMECs at P2, growth-arrested stasis HMECs (S), and variant HMECs at two passages beyond stasis (S + 2). In stasis cultures, quantitative proteomic analysis using an antibody array (Table S1 ) indicated increased levels of transforming growth factor-beta (TGF-β) and Activin A, a member of the TGF-β superfamily. ELISA assays confirmed higher levels of TGF-β, Activin A, and Plasminogen Activator Inhibitor-1 (PAI-1), a biomarker of TGF-β pathway activation, in stasis cultures than in proliferating pre-stasis or variant HMEC cultures (Fig. 1 E). Additionally, we detected an increase in matrix metalloproteinase 2 (MMP-2) and fibronectin (FN) transcript levels during stasis (Fig. 1 F), which can be upregulated in response to the activation of the TGF-β pathway [ 52 , 53 ]. Several cytokines and chemokines typical of SASP (IL-6, IL-8) [ 51 ] remained unchanged in the stasis cultures (Fig. S1 A). In stasis HMECs, we observed significant changes in gene expression indicative of EMT activation, which is consistent with the fundamental role of the TGF-β pathway in this process [ 54 ]. At the transcript level, we observed decreased levels of the epithelial markers E-cadherin (CDH1, concentrated in luminal epithelial cells) and P-cadherin (CDH3, typically expressed by basal epithelial cells), and no change in N-cadherin (CDH2, typically expressed by mesenchymal cells) in stasis HMECs compared to pre-stasis HMECs (Fig. 1 G). We also observed increased levels of core EMT transcription factors Slug (SNAI2), Snail (SNAI1), Twist (TWIST1), and Zeb (ZEB1). Transcript levels of E-cadherin, Snail, Twist, and Zeb decreased in variant cultures compared to those in stasis cultures and were not significantly different from those observed in pre-stasis HMECs (Fig. 1 G). In contrast, Slug transcript levels were increased in stasis HMECs and maintained in the variants. The derivation of variants from the basal compartment [ 9 ] likely accounts for the increased levels of P-cadherin observed in variant cultures compared with both pre-stasis and stasis cultures (Fig. 1 G). Western blot analyses were consistent with the increased EMT activation, particularly during stasis. We observed decreased E-cadherin expression and increased Snail, Zeb, N-cadherin, and Vimentin expressions. Notably, protein expression was inconsistent with the changes in certain transcript levels (N-cadherin, Snail, and Slug), suggesting regulation by post-transcriptional processes commonly observed in EMT (Fig. 1 H). In addition to TGF-β (Fig. 1 D), several paracrine and autocrine factors involved in EMT [ 54 ] were robustly increased in stasis HMECs compared with those in either of the two active growth phases (Fig. S1 B and S1C). In summary, our initial observations demonstrated transient activation of the TGF-β pathway and engagement of EMT programs in cultured HMECs as they entered stasis. Remarkably, the emergence of variants from these cultures was associated with the attenuation of the TGF-β and EMT pathways. Basal Epithelial Cells Exhibit Enhanced EMT Activation Compared to Luminal Epithelial Cells The human breast epithelium is composed of two distinct layers: an inner luminal layer of polarized columnar epithelial cells and an outer basal layer of contractile myoepithelial cells. Both layers comprise a heterogeneous mixture of mature cells and cells at various stages of differentiation, including stem and progenitor cells [ 1 – 6 ]. For simplicity, we will refer to these complex populations as 'luminal' and 'basal' cells, acknowledging that these terms encompass a range of cellular states and phenotypes. To reduce the confounding variable of comparing a mixture of luminal and basal derivatives in pre-stasis cultures with exclusive basal derivatives in variant cultures [ 9 ], we used fluorescence-activated cell sorting (FACS) to separate basal epithelial cells (EPCAM − CD49f + CD10 + ) from luminal epithelial cells (EPCAM + ) (Fig. 2 A). When cultured in serum-free MEGM, the basal-enriched population from the three tissue donors entered stasis between 9 and 15 population doublings. These cultures consistently produced variant outgrowths 30 days after the cessation of growth at P5 (passage 5). In contrast, luminal cultures proliferated between 8 and 11 population doublings before growth arrest at P3, and did not produce variant colonies (Fig. 2 B), consistent with a previous report [ 9 ]. Basal epithelial cells contributed significantly greater levels of TGF-β but not Activin A compared to luminal epithelial cells when analyzed at stasis, P3 for isolated luminal cells, and P5 for isolated basal epithelial cells (Fig. 2 C). As early as P2, even before entering stasis, basal-enriched populations demonstrated a more dispersed and migratory phenotype and incorporated vimentin, a canonical marker of the mesenchymal phenotype, into their cytoskeleton, whereas luminal cells maintained tight colony morphologies and did not express significant levels of vimentin (Fig. 2 D). At the transcriptional level, the basal epithelial population demonstrated a robust EMT response, as evidenced by significant increases in Twist, Zeb, and N-cadherin, and decreases in E-cadherin and P-cadherin between P2 and P5. Luminal epithelial cells demonstrated a significant increase in Slug and N-cadherin expression, and a decrease in E-cadherin expression (Fig. 2 E). Thus, both luminal and basal epithelial cells demonstrated evidence of EMT; however, the magnitude of EMT-related changes in gene expression was markedly higher in the basal epithelial cells. Similar to our observations in bulk HMECs, basal epithelial cell cultures were characterized by high levels of TGF-β and Activin A at early passages (P3 and P5), followed by a significant reduction in expression of these cytokines at later passages (P8) following the emergence of variants (Fig. 2 F). Western blot analysis revealed high levels of N-cadherin, Snail, and Zeb in early passage (P3) and stasis (P5) basal epithelial cells, indicating activation of EMT. However, these components of the EMT program were downregulated in the variant cultures (P8). Notably, among the EMT transcription factors analyzed, Slug levels were consistently high across all passages, including the variants (Fig. 2 G). A role of Slug in maintaining basal epithelial gene expression has been observed previously [ 55 ]. In summary, only the basal epithelial population contains variants [ 9 ]. This population recapitulated the transient activation of TGF-β and EMT phenotypes observed in bulk HMEC cultures, leading us to focus further on the basal epithelial cell population Inhibition of TGF-β Signaling and EMT Transcription Factors Prevents Stasis in Basal Epithelial Cells The application of A83-01, a potent inhibitor of TGF-β type I receptor ALK5 kinase, type I activin/nodal receptor ALK4, and type I nodal receptor ALK7, was sufficient to prevent stasis in basal epithelial cells (Fig. 3 A). A83-01 inhibited EMT at the transcriptional level, as evidenced by the significantly increased expression of P-cadherin and the decreased expression of N-cadherin, Snail, Twist, and Zeb (Fig. 3 B). TGF-β inhibition was associated with a highly significant decrease in cell size (Fig. 3 C), reduction in nuclear area, and elimination of micronuclei (Fig. 3 D). A83-01 treatment also significantly reduced p21 and p16 levels (Fig. 3 E). We observed that treatment with A83-01 maintained the protein expression of p63, an important transcriptional regulator of basal cell differentiation (Fig. 3 F). Together, these results implicate the engagement of the TGF-β pathway in cellular dysfunction and the activation of tumor suppressor pathways during stasis. Similar to the treatment of basal epithelial cells with A83-01, we found that CRISPR-mediated knockout (KO) of Snail or Slug eliminated the stasis barrier compared with non-targeting controls (Fig. 4 A). Following selection, Snail and Slug KO basal epithelial cells generated tight colony morphologies similar to those produced by the variant basal epithelial cultures emerging from stasis, whereas non-targeting controls were more dispersed across the culture surface. In contrast, Zeb and Twist KO basal epithelial cells demonstrated a growth pattern similar to that of the non-targeting control, and were therefore not followed up in culture (Fig. 4 A). Targeting Snail and Slug using CRISPR vectors reduced their transcript levels. Consistent with their role in the EMT network, Snail and Slug knockout significantly increased the transcript levels of P-cadherin and decreased those of N-cadherin, Twist, and Zeb compared to non-targeting controls (Fig. 4 B). Basal epithelial cells with Snail or Slug KO demonstrated substantially lower transcription of p21 and p16 tumor suppressors than the controls (Fig. 4 C). We observed vimentin in the cytoskeleton of basal epithelial cells expressing p16, consistent with the idea that EMT is correlated with the senescence phenotype in this culture system (Fig. S2 A). Notably, the reduced growth arrest observed upon Snail and Slug KO was associated with decreased levels of TGF-β and the TGF-β superfamily member activin A (Fig. 4 D), which are part of the autocrine and paracrine loops that promote the EMT program[ 54 ] and perpetuate senescence[ 49 ]. In contrast, overexpression (OE) of the EMT-related transcription factors Snail and Slug immediately suppressed the growth of basal epithelial cells, preventing their expansion (Fig. S2 B). In this experiment, Snail and Slug were expressed at non-physiological levels, which were far higher than those observed in basal epithelial cells at stasis. Analysis of these cultures at the first passage showed a high degree of EMT, as evidenced by the downregulation of P-cadherin and upregulation of N-cadherin, Twist, and Zeb (Fig. S2 C). Further examination of this model was not feasible because of the rapid and complete cessation of cell growth. In summary, TGF-β pathway inhibition or Snail or Slug knockdown prevented EMT activation and stasis in basal epithelial cells. Activation of the TGF-β pathway induces the expression of EMT-related transcription factors, such as Snail, Slug, Zeb, and TWIST, which promote the transition from an epithelial to a mesenchymal phenotype. In turn, mesenchymal cells resulting from EMT secrete additional TGF-β and Activin A creating a positive feedback loop that sustains and amplifies the EMT process. The activation of TGF-β and EMT in basal epithelial cells leads to growth suppression and cellular dysfunction, as is evident in conventional culture systems (Fig. 4 E). EMT Generates a Spectrum of Basal Epithelial Cell States with Varying Proliferative and Self-Renewal Capacities Limiting Activation of the TGF-β pathway and EMT induction was sufficient to avoid stasis in the entire cultured basal epithelial cell population. However, rare basal epithelial variants were intrinsically capable of avoiding stasis. Variants can be visualized as epithelial colonies (tightly clustered and adhered to their neighbors) emerging among cells with spindle or senescent-like morphologies (Fig. 5 A). EMT is a dynamic process that enables epithelial cells to acquire mesenchymal properties, leading to increased plasticity and the ability to transition between multiple cellular states [ 21 – 26 ]. This process creates a spectrum of intermediate phenotypes that allow cells to exhibit varying degrees of epithelial and mesenchymal characteristics. Many studies have used CD24 and CD44 to distinguish between epithelial and mesenchymal states [ 25 ]. Indeed, almost all cells in stasis cultures expressed the CD44 hi CD24 lo phenotype, consistent with the high levels of EMT at this time point. In contrast, pre-stasis HMECs and post-stasis variants mostly lacked CD44 (Fig. 5 B). However, this method was unable to separate epithelial, intermediate, and more complete mesenchymal phenotypes. To better distinguish between these EMT states in basal epithelial cell cultures, we developed an alternative method that utilized antibodies against P-cadherin, an epithelial marker highly expressed by basal epithelial cells, and N-cadherin, a marker highly expressed in mesenchymal cell types. In early and stasis cultures, we observed complete EMT states expressing N-cadherin but lacking P-cadherin (P − N + ), an intermediate EMT state expressing both P- and N-cadherins (P + N + ), and an epithelial state expressing P-cadherin but lacking N-cadherin (P + N − ). In cultures dominated by variants, the complete EMT state, P − N + , was mostly absent (Fig. 5 C). We sorted basal epithelial cells in the P + N + and P − N + states from stasis cultures (P5). At the mRNA level, P − N + cells exhibited significantly lower levels of P-cadherin and higher levels of N-cadherin, Snail, Slug, Twist, and Zeb than did P + N + cells, confirming that P − N + populations represent a higher degree of EMT (Fig. 5 D). Notably, basal epithelial cells with a complete EMT phenotype were growth arrested and could not be passaged, although they remained adherent to the culture surface for several weeks (Fig. 5 E). Early passage basal epithelial cells were sorted based on the immunophenotypes P + N − , P + N + , and P − N + ) corresponding to different states along the EMT spectrum, and cultured for 14 days. When the resulting populations were re-analyzed, both P + N − (blue) and P + N + (purple) cells were able to recapitulate the entire spectrum of EMT states, whereas P − N + (red) cells were restricted to the ‘fully’ mesenchymal state (Fig. 5 F). When separated from stasis cultures before the emergence of variants, we demonstrated, using limiting dilution analysis, that basal epithelial cells in the P + N − and P + N + states could give rise to variant colonies, whereas cells in the P − N + state could not (Fig. 5 G). In addition to the ability to form variant colonies in monolayer culture, we observed that basal epithelial cells in the P + N − and P + N + states could form mammospheres under non-adherent conditions, which is suggestive of self-renewal and stem cell properties [ 44 ]. In contrast, cells in the P − N + state were unable to create mammospheres (Fig. 5 H). The P − N + state was completely absent in the variant cultures, and cells could only move between the P + N − and P + N + states, suggesting attenuation of the EMT program (Fig. 5 I). In summary, isolated basal epithelial cells initially activated an EMT program that could progress to a fully mesenchymal state, ultimately leading to growth arrest. However, over several passages, we observed stabilization of the epithelial and intermediate EMT states, restricting access to the fully mesenchymal state. This restriction was associated with the emergence of variants that retained plasticity and could transition between epithelial and intermediate EMT states but not the fully mesenchymal state. PRC2-Mediated Epigenetic Regulation Governs Access to EMT States in Basal Epithelial Variants PRC2 is known to play a critical role in epigenetic regulation of gene expression in variant basal epithelial cells. Indeed, by catalyzing the trimethylation of histone H3 on lysine 27 (H3K27me3), PRC2 silences the target genes involved in growth and differentiation. Over time in culture, silencing of many PRC2 targets can be reinforced by promoter hypermethylation [ 35 , 41 ]. EZH2, the catalytic subunit of PRC2, is central to this function, and its knockout disrupts the ability of this complex to regulate gene expression. We found that the knockout of EZH2 in basal epithelial cells at P3 reliably prevented the formation of variant colonies, which was expected given that PRC2 plays a well-established role in silencing CDKN2A and removing other obstacles to proliferation (Fig. 6 A). EZH2 knockout in the P8 variant basal epithelial cells also impaired cell proliferation (Fig. 6 B). Recent studies have revealed that PRC2 plays a crucial role in controlling the EMT trajectory by repressing specific genes involved in EMT [ 56 ]. In EZH2 knockout basal epithelial variants, we confirmed a decrease in EZH2 expression and global H3K27me3 levels at the protein level, and observed increased levels of Zeb and N-cadherin (Fig. 6 C). EZH2 knockout promoted EMT at the transcriptional level, as evidenced by decreased P-cadherin, but increased N-cadherin, Snail, Twist, and Zeb expression (Fig. 6 D). By reducing the activation of the EMT program, PRC2 and potentially other epigenetic regulators maintain variant basal epithelial cells in an intermediate EMT state that favors continued proliferation and stemness, and disfavors the path to complete EMT states associated with growth arrest (Fig. 6 E). TGF-β Inhibition Prior to Cellular Transformation Reduced Metaplastic Phenotypes in Basal Epithelial-Derived Tumors Basal epithelial progenitors have been proposed as the cells of origin for MBC [ 9 ], a rare histological subtype associated with elevated levels of EMT compared with other breast cancer subtypes [ 18 ] (Fig. S3A and S3B). In culture, basal epithelial cells strongly activated the TGF-β pathway and induced EMT. The selection of variants was associated with the stabilization of intermediate EMT states capable of continued proliferation, limiting access to a complete mesenchymal state subject to growth arrest. Pharmacological inhibition of the TGF-β pathway suppressed the activation of the EMT program and eliminated the selective pressure that promoted the outgrowth of the basal variants (Fig. 7 A). To better understand the role of TGF-β and EMT-dependent selection of the variant phenotype in MBC carcinogenesis, we employed an established transformation protocol [ 9 ] involving the sequential transduction of TERT, SV40 T antigen, and mutant K-Ras G12V . We compared the oncogenic transformation of P3 basal epithelial cell populations cultured in MEGM for three passages, poised to enter stasis, with variant basal epithelial cell populations cultured in MEGM for eight passages that had already overcome stasis. Additionally, we compared P3 basal epithelial cell populations cultured in MEGM + DMSO with those cultured in MEGM + A83-01 for three passages, which suppressed EMT and allowed them to avoid stasis and variant selection. Notably, both the oncogenically transformed variant and A83-01 basal epithelial populations expressed levels of the mutant K-Ras G12V equivalent to those of the P3 controls (Fig. 7 B). Proliferation of oncogenically transformed P3 basal epithelial cells was suppressed, whereas the oncogenically transformed variant and A83-01 basal epithelial populations exhibited exponential growth kinetics (Fig. 7 C). Additionally, oncogenically transformed P3 basal epithelial cells formed dramatically fewer colonies in soft agar, with most replicates showing no colony formation at all, indicating a reduced capacity for anchorage-independent growth compared to the variant and A83-01 basal epithelial populations, which showed no significant difference in their ability to form colonies in soft agar (Fig. 7 D). A greater level of EMT activation was observed in oncogenically transformed P3 basal epithelial cells compared to both the variant and A83-01 populations, as evidenced by the lower levels of P-cadherin and higher levels of N-cadherin, Snail, Twist, and Zeb compared to the variant and A83-01 basal epithelial cell populations. In contrast, transformed variant populations exhibited intermediate levels of EMT, characterized by significantly higher levels of P-cadherin and lower levels of N-cadherin, Snail, Twist, and Zeb, compared to A83-01-treated basal epithelial populations (Fig. 7 E). Flow cytometric analysis of P-cadherin and N-cadherin expression confirmed that oncogenically transformed variant populations predominantly exhibited intermediate EMT phenotypes (P + N + ). In contrast, oncogenically transformed basal epithelial cells treated with A83-01 mostly expressed P-cadherin, with relatively little N-cadherin expression, thus mostly preserving their epithelial state (Fig. 7 F). Only the oncogeneically transformed variant and A83-01 basal epithelial cell populations formed tumors when xenotransplanted subcutaneously into NSG mice (Fig. S3C). Tumors derived from both A83-01-treated and variant basal epithelial cells were highly proliferative based on the number of cells that stained positive for Ki67 (Fig. S3D). Both sets of tumors expressed p63, which is a typical phenotype of MBCs, and several even rarer histological subtypes likely arising from basal epithelial progenitors [ 57 – 59 ]. Notably, A83-01 tumors exhibited a lower frequency of cells expressing the squamous differentiation marker cytokeratin 13, the EMT markers Snail and Slug, and the mesenchymal differentiation marker vimentin when compared to variant tumors (Fig. 7 G). In summary, high levels of EMT in the malignantly transformed P3 basal epithelial group were associated with poor growth and tumorigenicity, whereas the variant basal epithelial group, which exhibited intermediate levels of EMT gene expression, demonstrated enhanced growth and tumorigenicity. Inhibition of the TGF-β pathway, which significantly reduced EMT gene expression, promoted growth and tumorigenicity in the P3 basal epithelial group but led to a loss of metaplastic squamous and metaplastic elements. Discussion In conventional monolayer cultures [ 29 ], HMECs are deprived of instructive biophysical and biochemical cues that maintain cellular homeostasis and cooperation within tissues. Instead, they form attachments with a non-compliant plastic substratum and proliferate in response to signals presented in an engineered medium. Consequently, HMECs enter stasis, the first, ‘stress-associated’ barrier to indefinite proliferation, during which most cells adopt a senescence-like phenotype [ 8 , 33 , 34 ]. To address this limitation, more recently developed culture systems optimized to support epithelial proliferation incorporate fibroblast feeders or hydrogels containing basement membrane components, mimicking key aspects of the in vivo microenvironment lost in conventional monolayer cultures. These systems also employ chemical inhibitors to reduce the ability of the cells to engage in stress response pathways. Collectively, these modifications allow cultured breast epithelial cells to avoid both stasis and agonescence, the second barrier to indefinite proliferation resulting from telomere attrition, without the need for genetic manipulation. Notably, these culture systems support genetic and epigenetic stability during long-term propagation [ 60 , 61 ]. In such conventional monolayer cultures, rare variants activate an aberrant yet remarkably consistent epigenetic program that enables them to adapt to unfavorable culture conditions and reinitiate proliferation after encountering the stasis barrier [ 35 – 37 ]. Many epigenetic alterations favored in variant cultures, including hypermethylation of the CDKN2A promoter, have also been observed during human tumor progression [ 35 ]. In intact tissues, etiological factors, such as aging, chronic inflammation, and environmental exposure, gradually disrupt homeostatic interactions. This progressive disruption contributes to a decline in the quality of the tissue microenvironments, which is hypothesized to favor the selection of clones (epigenetic progenitors) that are adaptive to altered microenvironments, potentially leading to the emergence of cancer over the course of many decades [ 62 , 63 ]. Arguably, conventional monolayer cultures represent an analogous process, albeit with a significantly accelerated time scale. Thus, variants within these cultures have proven useful in elucidating the mechanisms underlying the emergence of cancer progenitor states [ 25 , 26 ]. In this study, we report that robust activation of the TGF-β pathway and induction of the EMT program are drivers of stasis and, therefore, constitute an important selective pressure underlying the emergence of the variant phenotype. Inhibition of the TGFβ pathway allowed basal epithelial cells to bypass stasis and dramatically reduced cell size, chromosomal instability, and expression of both p21 and p16. TGF-β plays a direct role in enhancing cell size by activating the mammalian target of rapamycin (mTOR) [ 64 ]. Notably, increased cell size is not simply a consequence of cellular senescence, but also a contributing factor through mechanisms that include failure to scale macromolecule production with cell volume causing cytoplasmic dilution, altered nuclear organization, chromosomal missegregation, DNA damage, and mitochondrial dysfunction [ 65 – 69 ]. TGF-β induces cellular senescence through SMADs, which can upregulate p15 and p21 and contribute to the establishment of SASP, which can reinforce senescence in an autocrine manner or induce senescence in neighboring cells in a paracrine manner [ 49 – 51 ]. Knockout of the key EMT factors Snail and Slug also allowed basal epithelial cells to bypass stasis. The EMT program is associated with increased activation of the TGFβ pathway through autocrine and paracrine secretion [ 54 ] but can also have negative consequences through other mechanisms. Indeed, during EMT, if cells continue to proliferate, they exhibit mitotic abnormalities such as centrosome clustering, cytokinesis failure, and chromosome missegregation, which result in increased DNA damage, aneuploidy, binucleated cells, and micronuclei. These mitotic defects are associated with the suppression of nuclear envelope proteins like LaminB1, critical for mitotic regulation [ 70 , 71 ]. Thus, TGFβ pathway activation and EMT are key components of the stasis barrier, resulting in cellular arrest and dysfunction. Notably, EMT is not a binary switch between epithelial and mesenchymal states but generates a spectrum of intermediates with both epithelial and mesenchymal traits. When basal epithelial cells underwent a complete transition to the mesenchymal state, which peaked at stasis, they were subject to proliferative arrest. In contrast, the re-initiation of proliferation within stasis cultures by variants was associated with the stabilization of intermediate states. These hybrid states are characterized by enhanced cellular plasticity and adaptability, which can lead to increased stem cell properties and tumorigenic potential [ 21 – 26 ]. PRC2 is known to contribute to the epigenetic program underlying the emergence of the variant phenotype by silencing genes involved in growth control and differentiation [ 35 , 37 , 41 ]. PRC2 also suppresses the expression of EMT-inducing transcription factors, such as ZEB1 and ZEB2 [ 56 ]. Indeed, in variant basal epithelial cells, knockout of EZH2, the catalytic subunit of PRC2, increased EMT gene expression, influencing the trajectory of the EMT program. Significantly, epigenetic alterations in the TGF-β pathway have been previously observed in variants [ 35 , 38 , 72 ]. EZH2 knockout in basal epithelial cells and variants impaired PRC2 function and diminished basal epithelial proliferation, limiting our study. In a prior study, fully transformed HMLER cells, which were derived from variant cultures, were not subject to growth arrest in response to compromised PRC2 activity; instead, they transitioned into a quasi-mesenchymal state associated with high metastatic potential. In contrast, the loss of another epigenetic regulator, KMT2D, results in a highly mesenchymal state [ 56 ]. Further research is warranted to elucidate the interplay between epigenetic regulators and EMT during the emergence of the variant phenotype. Nevertheless, our findings suggest that the adaptation of EMT toward favorable intermediate EMT states is a consequence of the epigenetic program, which facilitates the emergence of variants by modulating gene expression involved in growth control, differentiation, and other hallmark programs characteristic of a cancer progenitor state. Consistent with the idea that variants are representative of a cancer progenitor state, the variant state confers tolerance to oncogenic insults such as RAS activation, thereby bypassing oncogene-induced senescence [ 39 ]. In this study, we also observed the increased vulnerability of variants to oncogenic transformation when compared to pre-stasis basal epithelial cells, which had not undergone the same epigenetic reprogramming characteristic of the variant population. TGF-β inhibition, which bypassed stasis in basal epithelial cultures, improved the transformation efficiency of pre-stasis basal epithelial cells, suggesting that overcoming the cellular arrest and dysfunction resulting from activation of the TGFβ pathway and EMT is necessary for achieving malignant transformation in these cells. This result was consistent with a previous report that found that either a dominant-negative TGFβ type II receptor or a TGFβ pathway inhibitor was able to bypass RAS-induced growth arrest variants and HMECs with compromised p16/Rb and p53. This study reported no difference in oncogenic transformation efficiency between pre-stasis HMEC and variant HMEC populations [ 73 ]. However, they compared pre-stasis HMEC (specimen 48R, batch T) with post-selection HMEC (specimen 48R, batch S), which were developed in different media systems—serum-free MCDB170 medium (a non-commercially available version of MEGM) versus serum-containing M87A supplemented with cholera toxin and oxytocin [ 74 ], which differed from our analysis of basal epithelial cells cultured in MEGM. Strikingly, the oncogenic transformation of basal epithelial cells resulted in metaplastic tumors with squamous and mesenchymal components in immunocompromised mice, which is consistent with previous reports [ 8 , 9 ]. However, treatment with TGF-β pathway inhibitors generated tumors with a marked reduction in metaplastic elements, even though the TGF-β pathway inhibitor was not maintained in culture after selection for KRAS G12V or in mice. Basal cells adapt to conventional HMEC culture systems by restricting the EMT program to a hybrid epithelial/mesenchymal state that maintains a high proliferation, stemness, and tumorigenicity. This cellular plasticity is required to achieve a basal progenitor state that gives rise to rare MBCs. Notably, the MBC-like tumor phenotypes observed in the variant tumors show similarities to squamous cell carcinomas originating from basal epithelial cells in stratified epithelia throughout the body, such as lung, head and neck, esophageal, and cutaneous squamous cell carcinomas [ 75 ]. A subset of these squamous tumors are associated with mesenchymal differentiation, akin to that sometimes observed in MBCs [ 17 , 18 ], often referred to as spindle cells or sarcomatous elements [ 76 ]. Both MBC and squamous cell carcinoma are associated with inactivation of CDKN2A [ 77 – 81 ]. Chronic inflammation is often a contributory factor in the development of squamous cell carcinoma, often linked to tobacco use in lung, head, neck, and esophageal cancers, as well as to ultraviolet radiation in cutaneous squamous cell carcinomas [ 75 ]. MBC has been observed in the clinical context of chronic inflammation, such as breast abscess associated with a recurrent or long-standing breast cyst or implant capsule [ 82 ]. Evidence of intermediate EMT states has been identified in single-cell RNA-sequencing and immunohistochemical studies of squamous cell carcinomas, which are hypothesized to play a role in promoting aggressive tumor properties, including the formation of circulating tumor cell clusters that are highly metastatic and exhibit cancer stem cell characteristics [ 83 ]. Variants and their derivatives have been widely used to model breast epithelial cell biology and carcinogenesis, facilitating many influential observations. However, it is essential to note that these cultures selectively favor basal epithelial cells, which are thought to be the origin of a small fraction of breast cancers, mainly MBCs [ 9 ]. Despite this limitation, findings from HMEC culture systems are often extrapolated to breast cancer as a whole, although they may not accurately represent the biology of the more common breast cancer subtypes. Conclusions This study illuminated dynamic processes operative in the HMEC model system, which is widely used to study epithelial cell biology and carcinogenesis. Following disruption of their native microenvironment, activation of the TGF-β pathway and induction of EMT drove stasis in HMEC cultures, exerting selective pressure that favored the emergence of rare variant basal epithelial cells. The epigenetic stabilization of intermediate EMT states in variants contributed to their enhanced proliferative capacity and avoided complete mesenchymal differentiation, which was associated with irreversible growth arrest. The ability of basal epithelial cells to navigate the spectrum of EMT states and stabilize favorable intermediate states contributed to the emergence of a cancer progenitor state that gives rise to tumors with histopathological features consistent with rare MBCs and may overlap with tumor types initiated from basal epithelial cells, particularly squamous cell carcinomas, in other tissues. Abbreviations HMEC: Human mammary epithelial cell; CDH1: cadherin1 (E-cadherin); CDH2: cadherin2 (N-cadherin); CDH3: cadherin 3 (P-cadherin); CDKN2A: Cyclin-dependent kinase inhibitor 2A; CM: Conditioned medium; EMT: Epithelial-mesenchymal transition; EZH2: Enhancer of zeste homolog 2; H3K27me3: Trimethylation of lysine 27 on histone H3; MBC: Metaplastic breast cancer; MEGM: Mammary epithelial cell growth medium; PRC2: Polycomb repressive complex 2; SASP: Senescence-associated secretory phenotype; SEEP: Stress-elicited extrinsic phenotype; SNAI1: Snail family transcriptional repressor 1(Snail); SNAI2: Snail family transcriptional repressor 2 (Slug); TGF-β: Transforming growth factor-beta; TNBC: Triple-negative breast cancer; TWIST1: Twist family BHLH transcription factor 1 (Twist); ZEB1: Zinc finger E-box binding homeobox 1 (Zeb); Declarations Acknowledgments: The authors would like to acknowledge Dr. Philippe Gascard for his scientific input and assistance in editing the manuscript, Dr. Deng Pan for his scientific input, Chira Chen-Tanyolac for processing donor material, Xianhong Wang for assistance with experiments, the staff within the Biological Imaging Development CoLab (BIDC) at UCSF Parnassus Heights for their training and support in using the Leica SP8 Confocal microscope, and the Parnassus Flow Cytometry Core (PFFC) for their training and support in using the FACS Aria II Cell Sorter and LSRFortessa flow cytometer. Funding: This work was supported by the California Breast Cancer Research Program (CBCRP) Translational Research Award 22OB-0032 to T.D.T., National Cancer Institute R35 Outstanding Investigator Award 1R35CA197694-01 to T.D.T, and Grand Challenge Cancer Research UK award A27145 to T.D.T. The UCSF Parnassus Flow Cytometry Core (RRID: SCR_018206 ) was supported in part by Grant NIH P30 DK063720 and by the NIH S10 Instrumentation Grant S10 1S10OD021822-01. Availability of data and materials: All data generated for this study are available in the paper and its supplementary information files. Author Contributions: J.A.C. and T.D.T Designed research; J.A.C. Performed research; J.A.C. and T.D.T. Analyzed data; J.A.C. and T.D.T. Wrote the paper Ethics approval and consent to participate: All specimens were obtained with patient consent under the Institutional Review Board-approved protocol #10-01532. The samples were identified using unlinked codes, following the Federal Health Insurance Portability and Accountability Act (HIPAA) guidelines. Animal experiments were conducted following protocols approved by the University of California, San Francisco Institutional Animal Care and Use Committee. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests. References Van Keymeulen A, Rocha AS, Ousset M, Beck B, Bouvencourt G, Rock J, Sharma N, Dekoninck S, Blanpain C: Distinct stem cells contribute to mammary gland development and maintenance . Nature 2011, 479 (7372):189-193. Davis FM, Lloyd-Lewis B, Harris OB, Kozar S, Winton DJ, Muresan L, Watson CJ: Single-cell lineage tracing in the mammary gland reveals stochastic clonal dispersion of stem/progenitor cell progeny . Nat Commun 2016, 7 :13053. Stingl J, Eaves CJ, Zandieh I, Emerman JT: Characterization of bipotent mammary epithelial progenitor cells in normal adult human breast tissue . Breast Cancer Res Treat 2001, 67 (2):93-109. Villadsen R, Fridriksdottir AJ, Ronnov-Jessen L, Gudjonsson T, Rank F, LaBarge MA, Bissell MJ, Petersen OW: Evidence for a stem cell hierarchy in the adult human breast . J Cell Biol 2007, 177 (1):87-101. Gudjonsson T, Villadsen R, Nielsen HL, Ronnov-Jessen L, Bissell MJ, Petersen OW: Isolation, immortalization, and characterization of a human breast epithelial cell line with stem cell properties . Genes Dev 2002, 16 (6):693-706. Eirew P, Stingl J, Raouf A, Turashvili G, Aparicio S, Emerman JT, Eaves CJ: A method for quantifying normal human mammary epithelial stem cells with in vivo regenerative ability . Nat Med 2008, 14 (12):1384-1389. Lim E, Vaillant F, Wu D, Forrest NC, Pal B, Hart AH, Asselin-Labat ML, Gyorki DE, Ward T, Partanen A et al : Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers . Nat Med 2009, 15 (8):907-913. Ince TA, Richardson AL, Bell GW, Saitoh M, Godar S, Karnoub AE, Iglehart JD, Weinberg RA: Transformation of different human breast epithelial cell types leads to distinct tumor phenotypes . Cancer Cell 2007, 12 (2):160-170. Keller PJ, Arendt LM, Skibinski A, Logvinenko T, Klebba I, Dong S, Smith AE, Prat A, Perou CM, Gilmore H et al : Defining the cellular precursors to human breast cancer . Proc Natl Acad Sci U S A 2012, 109 (8):2772-2777. Van Keymeulen A, Lee MY, Ousset M, Brohee S, Rorive S, Giraddi RR, Wuidart A, Bouvencourt G, Dubois C, Salmon I et al : Reactivation of multipotency by oncogenic PIK3CA induces breast tumour heterogeneity . Nature 2015, 525 (7567):119-123. Molyneux G, Geyer FC, Magnay FA, McCarthy A, Kendrick H, Natrajan R, Mackay A, Grigoriadis A, Tutt A, Ashworth A et al : BRCA1 basal-like breast cancers originate from luminal epithelial progenitors and not from basal stem cells . Cell Stem Cell 2010, 7 (3):403-417. Cancer Genome Atlas N: Comprehensive molecular portraits of human breast tumours . Nature 2012, 490 (7418):61-70. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS et al : Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications . Proc Natl Acad Sci U S A 2001, 98 (19):10869-10874. Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, He X, Perou CM: Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer . Breast Cancer Res 2010, 12 (5):R68. Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA et al : Molecular portraits of human breast tumours . Nature 2000, 406 (6797):747-752. Santagata S, Thakkar A, Ergonul A, Wang B, Woo T, Hu R, Harrell JC, McNamara G, Schwede M, Culhane AC et al : Taxonomy of breast cancer based on normal cell phenotype predicts outcome . J Clin Invest 2014, 124 (2):859-870. Reddy TP, Rosato RR, Li X, Moulder S, Piwnica-Worms H, Chang JC: A comprehensive overview of metaplastic breast cancer: clinical features and molecular aberrations . Breast Cancer Res 2020, 22 (1):121. Taube JH, Herschkowitz JI, Komurov K, Zhou AY, Gupta S, Yang J, Hartwell K, Onder TT, Gupta PB, Evans KW et al : Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes . Proc Natl Acad Sci U S A 2010, 107 (35):15449-15454. Thiery JP, Acloque H, Huang RY, Nieto MA: Epithelial-mesenchymal transitions in development and disease . Cell 2009, 139 (5):871-890. Kalluri R, Weinberg RA: The basics of epithelial-mesenchymal transition . J Clin Invest 2009, 119 (6):1420-1428. Kroger C, Afeyan A, Mraz J, Eaton EN, Reinhardt F, Khodor YL, Thiru P, Bierie B, Ye X, Burge CB et al : Acquisition of a hybrid E/M state is essential for tumorigenicity of basal breast cancer cells . Proc Natl Acad Sci U S A 2019, 116 (15):7353-7362. Jolly MK, Boareto M, Huang B, Jia D, Lu M, Ben-Jacob E, Onuchic JN, Levine H: Implications of the Hybrid Epithelial/Mesenchymal Phenotype in Metastasis . Front Oncol 2015, 5 :155. Grosse-Wilde A, Fouquier d'Herouel A, McIntosh E, Ertaylan G, Skupin A, Kuestner RE, del Sol A, Walters KA, Huang S: Stemness of the hybrid Epithelial/Mesenchymal State in Breast Cancer and Its Association with Poor Survival . PLoS One 2015, 10 (5):e0126522. Morel AP, Lievre M, Thomas C, Hinkal G, Ansieau S, Puisieux A: Generation of breast cancer stem cells through epithelial-mesenchymal transition . PLoS One 2008, 3 (8):e2888. Mani SA, Guo W, Liao MJ, Eaton EN, Ayyanan A, Zhou AY, Brooks M, Reinhard F, Zhang CC, Shipitsin M et al : The epithelial-mesenchymal transition generates cells with properties of stem cells . Cell 2008, 133 (4):704-715. Chaffer CL, Brueckmann I, Scheel C, Kaestli AJ, Wiggins PA, Rodrigues LO, Brooks M, Reinhardt F, Su Y, Polyak K et al : Normal and neoplastic nonstem cells can spontaneously convert to a stem-like state . Proc Natl Acad Sci U S A 2011, 108 (19):7950-7955. Qin S, Jiang J, Lu Y, Nice EC, Huang C, Zhang J, He W: Emerging role of tumor cell plasticity in modifying therapeutic response . Signal Transduct Target Ther 2020, 5 (1):228. Corgnac S, Damei I, Gros G, Caidi A, Terry S, Chouaib S, Deloger M, Mami-Chouaib F: Cancer stem-like cells evade CD8(+)CD103(+) tumor-resident memory T (T(RM)) lymphocytes by initiating an epithelial-to-mesenchymal transition program in a human lung tumor model . J Immunother Cancer 2022, 10 (4). Hammond SL, Ham RG, Stampfer MR: Serum-free growth of human mammary epithelial cells: rapid clonal growth in defined medium and extended serial passage with pituitary extract . Proc Natl Acad Sci U S A 1984, 81 (17):5435-5439. Labarge MA, Garbe JC, Stampfer MR: Processing of human reduction mammoplasty and mastectomy tissues for cell culture . J Vis Exp 2013(71). Bissell MJ, Hall HG, Parry G: How does the extracellular matrix direct gene expression? J Theor Biol 1982, 99 (1):31-68. Discher DE, Mooney DJ, Zandstra PW: Growth factors, matrices, and forces combine and control stem cells . Science 2009, 324 (5935):1673-1677. Brenner AJ, Stampfer MR, Aldaz CM: Increased p16 expression with first senescence arrest in human mammary epithelial cells and extended growth capacity with p16 inactivation . Oncogene 1998, 17 (2):199-205. Romanov SR, Kozakiewicz BK, Holst CR, Stampfer MR, Haupt LM, Tlsty TD: Normal human mammary epithelial cells spontaneously escape senescence and acquire genomic changes . Nature 2001, 409 (6820):633-637. Locke WJ, Zotenko E, Stirzaker C, Robinson MD, Hinshelwood RA, Stone A, Reddel RR, Huschtscha LI, Clark SJ: Coordinated epigenetic remodelling of transcriptional networks occurs during early breast carcinogenesis . Clin Epigenetics 2015, 7 (1):52. Novak P, Jensen TJ, Garbe JC, Stampfer MR, Futscher BW: Stepwise DNA methylation changes are linked to escape from defined proliferation barriers and mammary epithelial cell immortalization . Cancer Res 2009, 69 (12):5251-5258. Huschtscha LI, Noble JR, Neumann AA, Moy EL, Barry P, Melki JR, Clark SJ, Reddel RR: Loss of p16INK4 expression by methylation is associated with lifespan extension of human mammary epithelial cells . Cancer Res 1998, 58 (16):3508-3512. Hinshelwood RA, Huschtscha LI, Melki J, Stirzaker C, Abdipranoto A, Vissel B, Ravasi T, Wells CA, Hume DA, Reddel RR et al : Concordant epigenetic silencing of transforming growth factor-beta signaling pathway genes occurs early in breast carcinogenesis . Cancer Res 2007, 67 (24):11517-11527. Dumont N, Crawford YG, Sigaroudinia M, Nagrani SS, Wilson MB, Buehring GC, Turashvili G, Aparicio S, Gauthier ML, Fordyce CA et al : Human mammary cancer progression model recapitulates methylation events associated with breast premalignancy . Breast Cancer Res 2009, 11 (6):R87. Berman H, Zhang J, Crawford YG, Gauthier ML, Fordyce CA, McDermott KM, Sigaroudinia M, Kozakiewicz K, Tlsty TD: Genetic and epigenetic changes in mammary epithelial cells identify a subpopulation of cells involved in early carcinogenesis . Cold Spring Harb Symp Quant Biol 2005, 70 :317-327. Reynolds PA, Sigaroudinia M, Zardo G, Wilson MB, Benton GM, Miller CJ, Hong C, Fridlyand J, Costello JF, Tlsty TD: Tumor suppressor p16INK4A regulates polycomb-mediated DNA hypermethylation in human mammary epithelial cells . J Biol Chem 2006, 281 (34):24790-24802. Elenbaas B, Spirio L, Koerner F, Fleming MD, Zimonjic DB, Donaher JL, Popescu NC, Hahn WC, Weinberg RA: Human breast cancer cells generated by oncogenic transformation of primary mammary epithelial cells . Genes Dev 2001, 15 (1):50-65. Rizki A, Weaver VM, Lee SY, Rozenberg GI, Chin K, Myers CA, Bascom JL, Mott JD, Semeiks JR, Grate LR et al : A human breast cell model of preinvasive to invasive transition . Cancer Res 2008, 68 (5):1378-1387. Dontu G, Abdallah WM, Foley JM, Jackson KW, Clarke MF, Kawamura MJ, Wicha MS: In vitro propagation and transcriptional profiling of human mammary stem/progenitor cells . Genes Dev 2003, 17 (10):1253-1270. Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, McQuaid S, Gray RT, Murray LJ, Coleman HG et al : QuPath: Open source software for digital pathology image analysis . Sci Rep 2017, 7 (1):16878. Schmidt U, Weigert M, Broaddus C, Myers G: Cell Detection with Star-Convex Polygons . In : 2018; Cham : Springer International Publishing; 2018: 265-273. Holst CR, Nuovo GJ, Esteller M, Chew K, Baylin SB, Herman JG, Tlsty TD: Methylation of p16(INK4a) promoters occurs in vivo in histologically normal human mammary epithelia . Cancer Res 2003, 63 (7):1596-1601. Crawford YG, Gauthier ML, Joubel A, Mantei K, Kozakiewicz K, Afshari CA, Tlsty TD: Histologically normal human mammary epithelia with silenced p16(INK4a) overexpress COX-2, promoting a premalignant program . Cancer Cell 2004, 5 (3):263-273. Acosta JC, Banito A, Wuestefeld T, Georgilis A, Janich P, Morton JP, Athineos D, Kang TW, Lasitschka F, Andrulis M et al : A complex secretory program orchestrated by the inflammasome controls paracrine senescence . Nat Cell Biol 2013, 15 (8):978-990. Fordyce CA, Patten KT, Fessenden TB, DeFilippis R, Hwang ES, Zhao J, Tlsty TD: Cell-extrinsic consequences of epithelial stress: activation of protumorigenic tissue phenotypes . Breast Cancer Res 2012, 14 (6):R155. Coppe JP, Desprez PY, Krtolica A, Campisi J: The senescence-associated secretory phenotype: the dark side of tumor suppression . Annu Rev Pathol 2010, 5 :99-118. Miettinen PJ, Ebner R, Lopez AR, Derynck R: TGF-beta induced transdifferentiation of mammary epithelial cells to mesenchymal cells: involvement of type I receptors . J Cell Biol 1994, 127 (6 Pt 2):2021-2036. Kim ES, Kim MS, Moon A: TGF-beta-induced upregulation of MMP-2 and MMP-9 depends on p38 MAPK, but not ERK signaling in MCF10A human breast epithelial cells . Int J Oncol 2004, 25 (5):1375-1382. Scheel C, Eaton EN, Li SH, Chaffer CL, Reinhardt F, Kah KJ, Bell G, Guo W, Rubin J, Richardson AL et al : Paracrine and autocrine signals induce and maintain mesenchymal and stem cell states in the breast . Cell 2011, 145 (6):926-940. Phillips S, Prat A, Sedic M, Proia T, Wronski A, Mazumdar S, Skibinski A, Shirley SH, Perou CM, Gill G et al : Cell-state transitions regulated by SLUG are critical for tissue regeneration and tumor initiation . Stem Cell Reports 2014, 2 (5):633-647. Zhang Y, Donaher JL, Das S, Li X, Reinhardt F, Krall JA, Lambert AW, Thiru P, Keys HR, Khan M et al : Genome-wide CRISPR screen identifies PRC2 and KMT2D-COMPASS as regulators of distinct EMT trajectories that contribute differentially to metastasis . Nat Cell Biol 2022, 24 (4):554-564. Tavassoli FA: Myoepithelial lesions of the breast. Myoepitheliosis, adenomyoepithelioma, and myoepithelial carcinoma . Am J Surg Pathol 1991, 15 (6):554-568. Tse GM, Tan PH, Chaiwun B, Putti TC, Lui PC, Tsang AK, Wong FC, Lo AW: p63 is useful in the diagnosis of mammary metaplastic carcinomas . Pathology 2006, 38 (1):16-20. Koker MM, Kleer CG: p63 expression in breast cancer: a highly sensitive and specific marker of metaplastic carcinoma . Am J Surg Pathol 2004, 28 (11):1506-1512. Rosenbluth JM, Schackmann RCJ, Gray GK, Selfors LM, Li CM, Boedicker M, Kuiken HJ, Richardson A, Brock J, Garber J et al : Organoid cultures from normal and cancer-prone human breast tissues preserve complex epithelial lineages . Nat Commun 2020, 11 (1):1711. Liu X, Ory V, Chapman S, Yuan H, Albanese C, Kallakury B, Timofeeva OA, Nealon C, Dakic A, Simic V et al : ROCK inhibitor and feeder cells induce the conditional reprogramming of epithelial cells . Am J Pathol 2012, 180 (2):599-607. Feinberg AP, Koldobskiy MA, Gondor A: Epigenetic modulators, modifiers and mediators in cancer aetiology and progression . Nat Rev Genet 2016, 17 (5):284-299. Feinberg AP, Ohlsson R, Henikoff S: The epigenetic progenitor origin of human cancer . Nat Rev Genet 2006, 7 (1):21-33. Lamouille S, Derynck R: Cell size and invasion in TGF-beta-induced epithelial to mesenchymal transition is regulated by activation of the mTOR pathway . J Cell Biol 2007, 178 (3):437-451. Neurohr GE, Terry RL, Lengefeld J, Bonney M, Brittingham GP, Moretto F, Miettinen TP, Vaites LP, Soares LM, Paulo JA et al : Excessive Cell Growth Causes Cytoplasm Dilution And Contributes to Senescence . Cell 2019, 176 (5):1083-1097 e1018. Manohar S, Estrada ME, Uliana F, Vuina K, Alvarez PM, de Bruin RAM, Neurohr GE: Genome homeostasis defects drive enlarged cells into senescence . Mol Cell 2023, 83 (22):4032-4046 e4036. Lanz MC, Zatulovskiy E, Swaffer MP, Zhang L, Ilerten I, Zhang S, You DS, Marinov G, McAlpine P, Elias JE et al : Increasing cell size remodels the proteome and promotes senescence . Mol Cell 2022, 82 (17):3255-3269 e3258. Demidenko ZN, Blagosklonny MV: Growth stimulation leads to cellular senescence when the cell cycle is blocked . Cell Cycle 2008, 7 (21):3355-3361. Crozier L, Foy R, Adib R, Kar A, Holt JA, Pareri AU, Valverde JM, Rivera R, Weston WA, Wilson R et al : CDK4/6 inhibitor-mediated cell overgrowth triggers osmotic and replication stress to promote senescence . Mol Cell 2023, 83 (22):4062-4077 e4065. Comaills V, Kabeche L, Morris R, Buisson R, Yu M, Madden MW, LiCausi JA, Boukhali M, Tajima K, Pan S et al : Genomic Instability Is Induced by Persistent Proliferation of Cells Undergoing Epithelial-to-Mesenchymal Transition . Cell Rep 2016, 17 (10):2632-2647. Knouse KA, Lopez KE, Bachofner M, Amon A: Chromosome Segregation Fidelity in Epithelia Requires Tissue Architecture . Cell 2018, 175 (1):200-211 e213. Lopez-Diaz FJ, Gascard P, Balakrishnan SK, Zhao J, Del Rincon SV, Spruck C, Tlsty TD, Emerson BM: Coordinate transcriptional and translational repression of p53 by TGF-beta1 impairs the stress response . Mol Cell 2013, 50 (4):552-564. Garbe JC, Vrba L, Sputova K, Fuchs L, Novak P, Brothman AR, Jackson M, Chin K, LaBarge MA, Watts G et al : Immortalization of normal human mammary epithelial cells in two steps by direct targeting of senescence barriers does not require gross genomic alterations . Cell Cycle 2014, 13 (21):3423-3435. Cipriano R, Kan CE, Graham J, Danielpour D, Stampfer M, Jackson MW: TGF-beta signaling engages an ATM-CHK2-p53-independent RAS-induced senescence and prevents malignant transformation in human mammary epithelial cells . Proc Natl Acad Sci U S A 2011, 108 (21):8668-8673. Johnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR: Head and neck squamous cell carcinoma . Nat Rev Dis Primers 2020, 6 (1):92. Fletcher CD, Beham A, Schmid C: Spindle cell haemangioendothelioma: a clinicopathological and immunohistochemical study indicative of a non-neoplastic lesion . Histopathology 1991, 18 (4):291-301. Yarbrough WG, Aprelikova O, Pei H, Olshan AF, Liu ET: Familial tumor syndrome associated with a germline nonfunctional p16INK4a allele . J Natl Cancer Inst 1996, 88 (20):1489-1491. Reed AL, Califano J, Cairns P, Westra WH, Jones RM, Koch W, Ahrendt S, Eby Y, Sewell D, Nawroz H et al : High frequency of p16 (CDKN2/MTS-1/INK4A) inactivation in head and neck squamous cell carcinoma . Cancer Res 1996, 56 (16):3630-3633. Gonzalez-Zulueta M, Shibata A, Ohneseit PF, Spruck CH, 3rd, Busch C, Shamaa M, El-Baz M, Nichols PW, Gonzalgo ML, Elbaz M et al : High frequency of chromosome 9p allelic loss and CDKN2 tumor suppressor gene alterations in squamous cell carcinoma of the bladder . J Natl Cancer Inst 1995, 87 (18):1383-1393. Foulkes WD, Flanders TY, Pollock PM, Hayward NK: The CDKN2A (p16) gene and human cancer . Mol Med 1997, 3 (1):5-20. Bartels S, van Luttikhuizen JL, Christgen M, Magel L, Luft A, Hanzelmann S, Lehmann U, Schlegelberger B, Leo F, Steinemann D et al : CDKN2A loss and PIK3CA mutation in myoepithelial-like metaplastic breast cancer . J Pathol 2018, 245 (3):373-383. Behranwala KA, Nasiri N, Abdullah N, Trott PA, Gui GP: Squamous cell carcinoma of the breast: clinico-pathologic implications and outcome . Eur J Surg Oncol 2003, 29 (4):386-389. Pal A, Barrett TF, Paolini R, Parikh A, Puram SV: Partial EMT in head and neck cancer biology: a spectrum instead of a switch . Oncogene 2021, 40 (32):5049-5065. Hu Y, Smyth GK: ELDA: extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays . J Immunol Methods 2009, 347 (1-2):70-78. Additional Declarations No competing interests reported. Supplementary Files SupplementalFigs.docx SupplementaryTable1.xlsx Supplemental Table 1: Quantitative Proteomic Analysis of Conditioned Media from HMECs Using Antibody Array. Quantitative proteomic analysis was performed using the Quantibody array-based multiplex ELISA system to simultaneously quantify 1,000 protein targets, including cytokines, growth factors, proteases, soluble receptors, and other proteins. Conditioned media (CM) were collected from three different donors: passage 2 (P2) HMECs, stasis HMECs (S), and variant HMECs at two passages beyond stasis (S+2). The table presents the expression levels of the analyzed proteins in pg/mL. An unpaired t-test with Welch's correction was performed on each row of the dataset. Following this, a two-stage linear step-up procedure (Benjamini, Krieger, and Yekutieli) was applied to control the False Discovery Rate (FDR) and adjust for multiple comparisons. SupplementalTable2.xlsx Supplemental Table 2: Characteristics of Donor Tissues Used in the Study. The table summarizes the tissue source, age, and ethnicity of the 12 donor tissues used in this study. Cite Share Download PDF Status: Published Journal Publication published 18 Dec, 2024 Read the published version in Breast Cancer Research → Version 1 posted Editorial decision: Revision requested 22 Oct, 2024 Reviews received at journal 22 Oct, 2024 Reviews received at journal 02 Oct, 2024 Reviewers agreed at journal 24 Sep, 2024 Reviewers agreed at journal 23 Sep, 2024 Reviewers agreed at journal 13 Sep, 2024 Reviewers invited by journal 05 Sep, 2024 Editor assigned by journal 29 Aug, 2024 Submission checks completed at journal 29 Aug, 2024 First submitted to journal 26 Aug, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4980285","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":358577027,"identity":"d4f32284-f032-4a9c-ba5c-620fc08fc07a","order_by":0,"name":"Joseph A. Caruso","email":"","orcid":"","institution":"University of California San Francisco","correspondingAuthor":false,"prefix":"","firstName":"Joseph","middleName":"A.","lastName":"Caruso","suffix":""},{"id":358577028,"identity":"bd1b2dd1-4136-408d-ad59-4b47bf8b0a17","order_by":1,"name":"Thea D. Tlsty","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA2UlEQVRIiWNgGAWjYBAC+xtAgvGfBZBkPgAk5IA4Ab8WNgkQaSDBwMPABlJqTJIWHgMitUg3H/wA1GJvL5HzTeLHHwMGfvYcA/xaZI4lSwC1JPZI5G6T7G0zYJDseUNAi0SOAUhLAg9QizRjwx8GgxuEbJHI//wD5DAeiZxn0gxAh9kT1pLDBrKFsQfIkGZgMwCyCWpJM7NIAPnlzDNjS6BfeCTOPCsgoCX58Y0PBjb27O3JD28AQ0yOvz15A14tYJAAIgTAJDB+iAf8B0hQPApGwSgYBSMKAABnojmcE/AiAwAAAABJRU5ErkJggg==","orcid":"","institution":"University of California San Francisco","correspondingAuthor":true,"prefix":"","firstName":"Thea","middleName":"D.","lastName":"Tlsty","suffix":""}],"badges":[],"createdAt":"2024-08-26 21:54:11","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4980285/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4980285/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s13058-024-01920-8","type":"published","date":"2024-12-18T15:58:35+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":65453027,"identity":"3715dba2-3780-4179-ad71-6d99d94d19ce","added_by":"auto","created_at":"2024-09-27 15:26:04","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1077166,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eTransient Activation of EMT in Stasis HMECs Precedes Variant Emergence.\u003c/strong\u003e\u003cem\u003e\u003cstrong\u003e \u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e(A)\u003c/strong\u003e Human breast tissue samples from three donors were used to generate HMEC cultures in MEGM. The first passage after partial trypsinization was considered passage 1 (P1). At each passage, the cells were counted, and the cumulative population doubling was calculated. The resulting growth curves illustrate three distinct phases: the initial logarithmic growth phase (a), the stasis growth plateau (b), and the second logarithmic growth phase initiated by the expansion of variants (c). Protein, RNA, and conditioned media (CM) were collected from passage 2 (P2) HMECs, stasis HMECs (S), and variant HMECs at two passages beyond stasis (S+2). Crystal violet-stained flasks depict the pattern of variant emergence from the stasis cultures. \u003cstrong\u003e(B)\u003c/strong\u003e Cell proliferation was measured in HMEC cultures using the Click-iT EdU cell proliferation assay and quantified by flow cytometry using cells collected from P2, S, and S+2 cultures, which were re-plated overnight and pulsed with 10 µM EdU for one hour. \u003cstrong\u003e(C)\u003c/strong\u003e Gene expression analysis of p21 and p16 was performed on cDNA produced from cells collected from P2, S, and S+2 cultures by qPCR. Representative images of P2 and S HMECs immunostained with antibodies against p16 and p21, and counterstained with DAPI are shown. \u003cstrong\u003e(D)\u003c/strong\u003e CM collected directly from P2 and S HMEC cultures (incubated with cells for 24 hours) was mixed with an equal amount of fresh MEGM and used to culture freshly isolated HMECs (P1). \u003cstrong\u003e(E)\u003c/strong\u003eSandwich ELISA kits were used to determine concentrations of TGF-β, Activin A, and PAI-1. An equal number of cells was harvested from P2, S, and S+2 cultures and cultured in basal media (supplements excluded) for 24 hours. Concentrations were normalized to the final cell number. \u003cstrong\u003e(F)\u003c/strong\u003e MMP-2 and FN were analyzed by qPCR. \u003cstrong\u003e(G)\u003c/strong\u003e E-cadherin, N-cadherin, P-cadherin, Slug, Snail, Twist, and Zeb were analyzed by qPCR. \u003cstrong\u003e(H)\u003c/strong\u003e Western blot analysis of lysates from P2, S, and S+2 cells using antibodies against Snail, Zeb, Slug, E-cadherin (E-Cad.), N-cadherin (N-Cad.), Vimentin, and Actin (loading control).\u003c/p\u003e","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/6eef1bbd284ccdd17b907683.png"},{"id":65453973,"identity":"699ab662-151c-471c-abf6-694b45dbb724","added_by":"auto","created_at":"2024-09-27 15:42:04","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1991285,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eBasal Epithelial Cells Exhibit Enhanced TGF-β Signaling and EMT Activation Compared to Luminal Epithelial Cells. (A) \u003c/strong\u003eHuman breast tissues from the same three donors utilized in Fig. 1A were further processed into single cells. The cell suspension was stained with conjugated antibodies against CD10, EpCAM, CD49f, CD2, CD3, CD16, CD64, CD31, and CD45 then DAPI and conjugated streptavidin. Dead cells (DAPI\u003csup\u003e+\u003c/sup\u003e), doublets, and immune/stromal cell types (CD2\u003csup\u003e+\u003c/sup\u003e, CD3\u003csup\u003e+\u003c/sup\u003e, CD16\u003csup\u003e+\u003c/sup\u003e, CD64\u003csup\u003e+\u003c/sup\u003e, CD31\u003csup\u003e+\u003c/sup\u003e, and/or CD45\u003csup\u003e+\u003c/sup\u003e) were removed. The breast epithelium was separated into luminal (EpCAM\u003csup\u003e+\u003c/sup\u003e) and basal (EpCAM\u003csup\u003e-\u003c/sup\u003eCD49f\u003csup\u003e+\u003c/sup\u003eCD10\u003csup\u003e+\u003c/sup\u003e) enriched populations by FACS and expanded separately in MEGM in 6-well plates. Representative wells were immunostained prior to passaging (P0) with the luminal markers cytokeratin 19 (CK19) and mucin 1 (MUC1) or the basal markers p63 and cytokeratin 14 (CK14). \u003cstrong\u003e(B)\u003c/strong\u003e Basal and luminal epithelial cells were further expanded in MEGM, cells were counted at each passage, and the growth curve of the cumulative population doubling over time was plotted. Cells and CM were collected from basal epithelial cells at passages B(P2), B(P5), and B(P7) and from luminal cells at passages L(P2) and L(P3). Crystal violet-stained flasks depict the pattern of variant cell emergence from basal epithelial cell cultures after stasis. \u003cstrong\u003e(C)\u003c/strong\u003e Sandwich ELISA \u0026nbsp;kits were used to determine the TGF-β and Activin A concentrations. An equal number of cells was harvested from the B(P2), B(P5), L(P2), and L(P3) cultures and then cultured in basal media (supplements excluded) for 24 h. Concentrations were normalized to the final cell number. \u003cstrong\u003e(D)\u003c/strong\u003e Representative images of P2 basal and luminal epithelial cells immunostained with antibodies against pan-cytokeratin and vimentin and then counterstained with DAPI. \u003cstrong\u003e(E)\u003c/strong\u003e E-cadherin, P-cadherin, N-cadherin, Slug, Snail, Twist, and Zeb were analyzed by qPCR using cDNA produced from B(P2), B(P5), L(P2), and L(P3) cells. \u003cstrong\u003e(F)\u003c/strong\u003e Same as C, using basal epithelial cells from the B(P3), B(P5), and B(P8) cultures. \u003cstrong\u003e(G)\u003c/strong\u003e Western blot analysis of lysates from B(P2), B(P5), and B(P7) cultures, using antibodies against N-cadherin (N-Cad.), Snail, Zeb, Slug, and Actin (as loading controls).\u003c/p\u003e","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/34be82018308acb7df9a5f99.png"},{"id":65453031,"identity":"eaf078f9-35ba-4fee-bcd0-fa3e1581f9b6","added_by":"auto","created_at":"2024-09-27 15:26:04","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":2284022,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eInhibition of the TGF-β Pathway Prevents Stasis in Basal Epithelial Cells. (A)\u003c/strong\u003e Basal epithelial cells, generated as introduced in Fig. 2A, were expanded in MEGM with 500 nM A83-01 or DMSO, the cells were counted, and the growth curve of the cumulative population doubling over time was plotted. Cells and CM were collected at P3, as indicated by arrows.\u003cstrong\u003e (B)\u003c/strong\u003eP-cadherin, N-cadherin, Slug, Snail, Twist, and Zeb were analyzed by qPCR using cDNA produced from basal epithelial cells at P3 treated with A83-01 or DMSO. \u003cstrong\u003e(C)\u003c/strong\u003eRepresentative images of P3 basal epithelial cells treated with A83-01 or DMSO immunostained with an antibody against β-tubulin and phalloidin, which stain F-actin, and counterstained with DAPI. The cells were segmented using ImageJ software, and the cell size was calculated. Data from the four donors are combined in the graph; each black circle represents a single cell (DMSO, n=94; A83-01, n=78).\u003cstrong\u003e (D) \u003c/strong\u003eSame as in C for cells immunostained with an antibody against Lamin A/C and counterstained with DAPI. Nuclei were segmented using ImageJ software and the nucleus size was calculated (DMSO, n=82; A83-01, n=91). The percentage of cells with micronuclei was determined by staining fixed cells dropped on glass slides with Giemsa and determining the number of nuclei with associated micronuclei as a percentage of the total number of nuclei imaged (DMSO, n=549; A83-01, n=528). \u003cstrong\u003e(E)\u003c/strong\u003e p21 and p16 were analyzed by qPCR using cDNA produced from basal epithelial cells at P3 treated with A83-01 or DMSO. Representative images of basal epithelial cells treated with A83-01 or DMSO immunostained with antibodies against p16 and p21, and counterstained with DAPI. \u003cstrong\u003e(F)\u003c/strong\u003e Representative images of basal epithelial cells cultured in MEGM or MEGM containing 500 nM A83-01 were immunostained with antibodies against p63 and cytokeratin 14 (CK14), followed by counterstaining with DAPI. Western blot analysis of cell lysates from these populations was performed using antibodies against CK14, p63, SMA, and Actin (loading control).\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/16ff41572c9409692f5c1477.png"},{"id":65453385,"identity":"401d712f-9469-4f8a-8189-6f4367ccc15f","added_by":"auto","created_at":"2024-09-27 15:34:04","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":552158,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eDepletion of EMT Transcription Factors Prevents Stasis in Basal Epithelial Cells. (A)\u003c/strong\u003e Basal epithelial cells (P1) were transduced with lentiviral CRISPR constructs, containing CAS9 and sgRNA targeting Slug, Snail, Zeb, Twsit, or non-targeting controls (Cont.). These constructs also included an Internal Ribosome Entry Site (IRES)-GFP sequence, facilitating the selection of cells incorporating the vector using FACS. Snail knockout (KO), Slug KO, and control cells were expanded in MEGM, cells were counted at each passage, and a growth curve of the cumulative population doubling over time was plotted. Cells and CM were collected at P2+3 (P2 indicates that the transduced cells were at passage 2, and +3 indicates the number of passages following GFP\u003csup\u003e+\u003c/sup\u003e selection), as indicated by the arrows. \u003cstrong\u003e(B)\u003c/strong\u003e P-cadherin, N-cadherin, Slug, Snail, Twist and Zeb were analyzed by qPCR using cDNA produced from basal epithelial cells at P2+3, selected for integration of Snail KO, Slug KO, or the control vector. \u003cstrong\u003e(C)\u003c/strong\u003e p21 and p16 were analyzed as described in B. \u003cstrong\u003e(D) \u003c/strong\u003eSandwich ELISA \u0026nbsp;kits were used to determine the concentrations of TGF-β and activin A . An equal number of cells were harvested from basal epithelial cells at P2+3 selected for integration of Snail KO, Slug KO, or control vector, and cultured in basal media (supplements excluded) for 24 hours. Concentrations were normalized to the final cell number. \u003cstrong\u003e(E)\u003c/strong\u003e Graphic illustration of the key findings presented in Figs. 3 and 4.\u0026nbsp;In control basal epithelial cell cultures in serum-free\u0026nbsp;MEGM, positive feedback between activation of the\u0026nbsp;TGF-β\u0026nbsp;pathway and induction of\u0026nbsp;EMT\u0026nbsp;generates an environment that suppresses proliferation, causing stasis. In basal epithelial cell cultures treated with a\u0026nbsp;TGF-β\u0026nbsp;inhibitor or transduced with lentiviral\u0026nbsp;CRISPR\u0026nbsp;vectors targeting\u0026nbsp;Snail\u0026nbsp;or\u0026nbsp;Slug, the levels of TGF-β pathway activation and\u0026nbsp; EMT induction, as well as the positive feedback loop between them, were diminished, permitting proliferation and avoiding stasis.\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/c745b41b7179f1f9944b0149.png"},{"id":65453971,"identity":"9f17072a-18f6-4384-95bc-12f45bc06ce0","added_by":"auto","created_at":"2024-09-27 15:42:04","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":2201917,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eEMT Generates a Spectrum of Basal Epithelial Cell States with Varying Proliferative and Self-Renewal Capacities\u003c/strong\u003e. \u003cstrong\u003e(A) \u003c/strong\u003ePhase contrast images of P5 basal epithelial cells demonstrating different morphologies. \u003cstrong\u003e(B)\u003c/strong\u003e Dead cells (DAPI\u003csup\u003e+\u003c/sup\u003e) and doublets were removed, and CD44 and CD24 were visualized using flow cytometry. The CD44\u003csup\u003ehi\u003c/sup\u003eCD24\u003csup\u003elo\u003c/sup\u003e, CD44\u003csup\u003ehi\u003c/sup\u003eCD24\u003csup\u003ehi\u003c/sup\u003e, CD44\u003csup\u003elo\u003c/sup\u003eCD24\u003csup\u003ehi\u003c/sup\u003e, and CD44\u003csup\u003elo\u003c/sup\u003eCD24\u003csup\u003elo\u003c/sup\u003e immunophenotypes were quantified. Each bar represents the basal epithelial cells isolated from a single donor. \u003cstrong\u003e(C)\u003c/strong\u003e Same as B, using antibodies against P-cadherin and N-cadherin. Immunophenotypes were gated based on singly stained controls as P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e, and P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e. \u003cstrong\u003e(D)\u003c/strong\u003e P-cadherin, N-cadherin, Snail, Slug, Twist, and Zeb were analyzed by qPCR using cDNA produced from basal epithelial cells \u0026nbsp;(P5) directly sorted using FACS based on a P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e or\u0026nbsp; P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e immunophenotype. \u003cstrong\u003e(E)\u003c/strong\u003e P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e or\u0026nbsp; P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e\u0026nbsp; basal epithelial cells (P5) were expanded in MEGM and counted at each passage. \u003cstrong\u003e(F)\u003c/strong\u003e P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e, and P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e basal epithelial cells (P5) were sorted into 6-well plates, cultured for 14 days in MEGM, and re-analyzed by flow cytometry. The illustration on the right shows the proportion of cells that transitioned between states in three donors: P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e in blue, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e in purple, and P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e in red. \u003cstrong\u003e(G)\u003c/strong\u003e\u0026nbsp; P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e,\u003csup\u003e \u003c/sup\u003eand P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e basal epithelial cells (P5) were sorted into 6-well plates at concentrations of 100, 500, 1000, and 2000 cells and cultured in MEGM until variant colonies were observed. Estimates of variant frequency in the sorted populations and p-values were determined by Extreme Limiting Dilution Analysis [84] using the online resource:\u0026nbsp;\u003ca href=\"https://bioinf.wehi.edu.au/software/elda/\"\u003ehttps://bioinf.wehi.edu.au/software/elda/\u003c/a\u003e. \u003cstrong\u003e(H)\u003c/strong\u003e P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e, and P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e basal epithelial cells (P5) were sorted into ultra-low-adherence plates containing mammosphere medium [44] (DMEM/F12 supplemented with B27, EGF, FGF2, and heparin). The mammospheres were visualized after 14 days and counted using a phase-contrast microscope. \u003cstrong\u003e(I)\u003c/strong\u003e P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e and P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e cells were sorted from P8 (variant) basal epithelial cell cultures and re-analyzed by flow cytometry after 14 days. The illustration on the right shows the proportion of cells that transitioned between states in three donors: P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e-\u003c/sup\u003e in blue and P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e in purple. Note that very few P\u003csup\u003e-\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e cells were observed in variant cultures (see panel C).\u0026nbsp;\u003c/p\u003e","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/9472d2bf3989fe104f0478de.png"},{"id":65454163,"identity":"272e83b9-2932-4e90-b410-2201d41e5de6","added_by":"auto","created_at":"2024-09-27 15:50:04","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":564317,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003ePRC2-Mediated Epigenetic Regulation Governs Access to EMT States in Basal Epithelial Variants. \u003c/strong\u003e\u0026nbsp;\u003cstrong\u003e(A)\u003c/strong\u003e Sorted basal epithelial cells (P3) were transduced with lentiviral CRISPR constructs containing CAS9 and sgRNA targeting EZH2 or non-targeting controls (Cont.). These constructs included an IRES-GFP sequence, which facilitated the selection of cells incorporating the vector using FACS. EZH2 and control cells were expanded in MEGM; at each passage, cells were counted, and a growth curve of the cumulative population doubling over time was plotted. \u003cstrong\u003e(B) \u003c/strong\u003eSame as in A, starting with P8 (variant) basal epithelial cells. Cells were collected at P8+3, as indicated by the arrows. \u003cstrong\u003e(C)\u003c/strong\u003eWestern blot analysis of lysates from P8 EZH2 KO and control cells using antibodies against Zeb, N-cadherin, EZH2, H3K27me3, total Histone H3, and Actin ( loading control). \u003cstrong\u003e(D)\u003c/strong\u003e P-cadherin, N-cadherin, Slug, Snail, Twist, and Zeb were analyzed by qPCR using P8+3 EZH2 KO and control cells. \u003cstrong\u003e(E)\u003c/strong\u003e In the EMT continuum, cells exist in various states ranging from epithelial to mesenchymal phenotypes, including metastable intermediate transition states. Epigenetic processes involving PRC2 and likely other mechanisms are altered during adaptation to unfavorable cell culture conditions. These alterations influence the transitions between states and contribute to the stabilization of cells within intermediate EMT states.\u003c/p\u003e","description":"","filename":"Figure6.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/cff4a2428df54806a9be56c3.png"},{"id":65453035,"identity":"bf67b083-fd40-4e9c-aa38-210db8a20378","added_by":"auto","created_at":"2024-09-27 15:26:04","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":3933505,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eInhibition of the TGF-β pathway during cellular transformation reduces metaplastic phenotypes in basal epithelial-derived tumors\u003c/strong\u003e. \u003cstrong\u003e(A)\u003c/strong\u003e Activation of the TGF-β pathway and EMT collaborate to drive cellular dysfunction and growth arrest. Adapted variants that emerged from stasis cultures exhibit notable changes in gene expression associated with growth and differentiation and are stabilized in an intermediate EMT state. Inhibition of the TGF-β pathway eliminates these selective pressures during the expansion of basal epithelial cells. \u003cstrong\u003e(B)\u003c/strong\u003e Basal epithelial cell populations: (1) P3, cultured in MEGM + DMSO for three passages, (2) variant, cultured in MEGM + DMSO for eight passages, and (3) A83-01, cultured in MEGM + 500 nM A83-01 for three passages, were oncogenically transformed through sequential transduction with TERT and SV40, followed by selection (puromycin and hygromycin)\u0026nbsp; and lentiviral transduction of oncogenic KRAS\u003csup\u003eG12V\u003c/sup\u003e followed by a second round of selection (blasticidin and IRES-RFP). Western blot analysis was performed on lysates from these oncogenically transformed basal epithelial cells using antibodies against RAS and Actin (as a loading control). \u003cstrong\u003e(C)\u003c/strong\u003e At each passage, Oncogenically transformed basal epithelial cells were counted. A growth curve of the cumulative population doubling over time was plotted. \u003cstrong\u003e(D)\u003c/strong\u003e Transformed basal epithelial populations were suspended in soft agar in each well of a 6-well plate at a concentration of 20,000 cells/well and cultured in DMEM/F12 + 20% FBS. After four weeks, the resulting colonies were stained with crystal violet, imaged, and counted using the ImageJ software. \u003cstrong\u003e(E)\u003c/strong\u003e P-cadherin, N-cadherin, Slug, Snail, Twist, and Zeb were analyzed qPCR using cDNA produced from oncogenically transformed basal epithelial cell populations. \u003cstrong\u003e(F)\u003c/strong\u003e Variant and A83-01-treated oncogenically transformed basal epithelial cells were immunostained for P-cadherin and N-cadherin, then visualized by flow cytometry. \u003cstrong\u003e(G)\u003c/strong\u003e Variants and A83-01-treated oncogenically transformed basal epithelial cells (4 × 10\u003csup\u003e6\u003c/sup\u003e) were injected subcutaneously with 50% Matrigel into 8-to 12-week-old female NSG mice. Tumors formed over a period of 3-6 months. FFPE sections from five tumors were subjected to multiplexed immunohistochemical analysis for p63, cytokeratin 13 (CK13), Slug + Snail, and cytokeratin 5 + vimentin. Human-specific lamin A/C was used to identify the tumor cells.\u003c/p\u003e","description":"","filename":"Figure7.png","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/3417f9030473501c8e92cb3d.png"},{"id":72202733,"identity":"88064db2-45f0-4a3f-8b5c-c04423333fd9","added_by":"auto","created_at":"2024-12-23 16:15:48","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":15629116,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/9b388b5d-c37f-4fb7-bfb3-96a53a3554c5.pdf"},{"id":65453388,"identity":"49459030-15ea-4f85-836b-35b9c116fa98","added_by":"auto","created_at":"2024-09-27 15:34:04","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":6071399,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cbr\u003e\u003c/p\u003e","description":"","filename":"SupplementalFigs.docx","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/6385350f3fc027ac719f80b1.docx"},{"id":65453033,"identity":"6ae6125a-91c8-48f2-a985-9a068027455c","added_by":"auto","created_at":"2024-09-27 15:26:04","extension":"xlsx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":303138,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSupplemental Table 1: Quantitative Proteomic Analysis of Conditioned Media from HMECs Using Antibody Array. \u003c/strong\u003eQuantitative proteomic analysis was performed using the Quantibody array-based multiplex ELISA system to simultaneously quantify 1,000 protein targets, including cytokines, growth factors, proteases, soluble receptors, and other proteins. Conditioned media (CM) were collected from three different donors: passage 2 (P2) HMECs, stasis HMECs (S), and variant HMECs at two passages beyond stasis (S+2). The table presents the expression levels of the analyzed proteins in pg/mL. An unpaired t-test with Welch's correction was performed on each row of the dataset. Following this, a two-stage linear step-up procedure (Benjamini, Krieger, and Yekutieli) was applied to control the False Discovery Rate (FDR) and adjust for multiple comparisons.\u003c/p\u003e","description":"","filename":"SupplementaryTable1.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/b4ee44bbf910084b3be6f1c7.xlsx"},{"id":65453389,"identity":"6c43b678-8dee-4a00-8e86-0e65958c04a7","added_by":"auto","created_at":"2024-09-27 15:34:04","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":10007,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSupplemental Table 2: Characteristics of Donor Tissues Used in the Study. \u003c/strong\u003eThe table summarizes the tissue source, age, and ethnicity of the 12 donor tissues used in this study.\u003c/p\u003e","description":"","filename":"SupplementalTable2.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-4980285/v1/ffb133696c6ad64b146b89b5.xlsx"}],"financialInterests":"No competing interests reported.","formattedTitle":"An adaptive Epithelial-Mesenchymal Transition Program Enables Basal Epithelial Cells to Bypass Stress-Induced Stasis and Contributes to Metaplastic Breast Cancer Progenitor State","fulltext":[{"header":"Background","content":"\u003cp\u003eThe milk-producing lobules and ductal networks within a healthy human breast consist of an inner luminal layer of polarized columnar epithelial cells and an outer basal layer of contractile myoepithelial cells. Both luminal and basal epithelia contain hierarchies of progenitor cells [\u003cspan additionalcitationids=\"CR2 CR3 CR4 CR5\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e], whose transformation significantly contributes to the histological and molecular heterogeneity observed among breast cancers [\u003cspan additionalcitationids=\"CR8 CR9 CR10 CR11\" citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. These cancers are commonly classified into subtypes that can express markers of luminal or basal differentiation states [\u003cspan additionalcitationids=\"CR14\" citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]. However, most, including those with characteristically basal expression patterns, are believed to originate from progenitor cell populations in the luminal epithelium [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e, \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. Only rare metaplastic breast cancers (MBCs), representing 0.2-5 percent of all breast cancers [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e], are proposed to originate from basal progenitor populations [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e] primarily responsible for the maintenance of the myoepithelial layer [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eMBC is a histologically defined classification typified by differentiation into cell types that are not normally found in breast tissue, including squamous and mesenchymal (e.g., spindle cell, chondroid, and osseous metaplasia) phenotypes generating foci of bone, muscle and cartilage cells. These cancers predominantly fall within claudin-low or basal-like intrinsic molecular subtypes and are generally clinically classified as triple-negative breast cancers (TNBCs) [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Epithelial-mesenchymal transition (EMT) is closely associated with the pathogenesis of MBC and is believed to contribute to its comparatively poor prognosis, chemoresistance, and high metastatic potential. EMT is a developmental process in which epithelial cells acquire mesenchymal characteristics, including increased motility and invasiveness, which can be reactivated in response to injurious or pathological conditions, including carcinogenesis [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e, \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e]. Notably, EMT is not a binary switch from the epithelial to mesenchymal state. Instead, it produces a spectrum of intermediate cellular states that exhibit both epithelial and mesenchymal characteristics. These intermediate states are associated with an increase in stem cell properties, tumorigenicity, drug resistance, and metastasis [\u003cspan additionalcitationids=\"CR22 CR23 CR24 CR25 CR26 CR27\" citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThe human mammary epithelial cell (HMEC) culture system, developed more than 40 years ago [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e] by the Stampfer group, remains one of the most widely used platforms for studying epithelial cell biology and carcinogenesis. Establishing HMEC cultures requires mincing healthy human breast tissue, enzymatically digesting matrix components, and removing stromal cell types [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e, \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. In healthy tissues, cellular cooperation relies on complex communication networks involving a vast array of biochemical and biophysical interactions, including interactions between parenchymal and stromal cells, the extracellular matrix, and the systemic milieu [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. Constantly integrating information from these networks defines a cell's operational state and sets limits on cellular behavior [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e]. Indeed, after 5\u0026ndash;20 population doublings in culture, HMECs enter stasis, a 'stress-associated' barrier that is dependent on the engagement of the tumor suppressor retinoblastoma (Rb) and characterized by most cells adopting a senescence-like phenotype [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. Variants, exclusively derived from the basal compartment [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e], evade this fate and initiate a second phase of logarithmic growth as clonal outgrowths from stasis cultures. These variant basal epithelial cells (vHMECs) continue to experience telomere erosion, allowing for an additional 20\u0026ndash;70 population doublings before reaching a telomeric crisis state known as agonescence [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. Because of their extended lifespans, variants are more broadly available from commercial sources and are often marketed simply as HMECs. These cells have been immortalized and oncogenically transformed into cell lines that are widely used in laboratories worldwide to model breast epithelial biology and carcinogenesis. Thus, vHMECs have profoundly influenced many fields, including the study of stromal-epithelial interactions, aging, EMT, and stem cell biology.\u003c/p\u003e \u003cp\u003eVariant basal epithelial outgrowths that emerge from stasis and dominate post-selection HMEC cultures exhibit compromised growth control and acquire other characteristics associated with malignant phenotypes. Many of these alterations result from the activation of an epigenetic program, which is surprisingly consistent between donors [\u003cspan additionalcitationids=\"CR36 CR37 CR38 CR39 CR40\" citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. Characteristic epigenetic alterations include hypermethylation of specific genes and loci targeted by polycomb repressor complex 2 (PRC2), including the promoters of genes regulating growth control and differentiation, such as cyclin-dependent kinase inhibitor 2A (CDKN2A) and homeobox A9 [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e, \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. In addition to alterations in the chromatin landscape, cytogenetic analysis reveals an accumulation of chromosomal abnormalities in variant cultures, including translocations, deletions, telomeric associations, polyploidy, and aneuploidy [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. Collectively, modifications in the variant HMEC state are sufficient to allow these cells to tolerate oncogenic events, such as the introduction of oncogenic forms of RAS, without undergoing oncogene-induced senescence [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. Additional factors, such as the viral oncoprotein SV40 large T antigen and telomerase reverse transcriptase (TERT), are required to fully immortalize these cells, which is a prerequisite for the introduction of a transforming oncogene to achieve tumorigenicity in immunocompromised mice [\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e]. Notably, the malignant transformation of variant cultures characteristically results in the development of tumors that recapitulate the histopathological characteristics of the rare MBC subtype. Thus, many widely used breast cancer models derived from HMEC cultures, such as MCF-10AT, HMT-3522 T4-2, and HMLER, form tumors exhibiting metaplastic squamous and spindle cell differentiation [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e, \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e, \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn this study, we observed robust activation of the TGF-β pathway and EMT as HMEC cultures entered stasis, which was followed by the attenuation of these programs in emerging variant basal epithelial populations. In this system, activation of the TGF-β pathway and EMT was associated with cellular dysfunction, driving cells toward a senescent-like phenotype. Rare variants that spontaneously emerged from stasis cultures exhibited a restricted EMT trajectory, which prevented their complete transition to a mesenchymal state characterized by P-cadherin loss and irreversible growth arrest, instead confining them to dynamic transitions between epithelial and intermediate EMT states that maintained P-cadherin expression. The aberrant epigenetic program activated in variants not only compromised growth control and differentiation [\u003cspan additionalcitationids=\"CR36 CR37 CR38 CR39 CR40\" citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e] but also contributed to the restricted EMT program observed in variants. The variant phenotype was associated with increased susceptibility to oncogenic transformation compared to unselected, pre-stasis basal epithelial cells and produced tumors with metaplastic elements. Pharmacological inhibition of the TGF-β pathway in cultured basal epithelial cells limited EMT program activation, reduced the selective pressures underlying variant emergence, and enhanced tumorigenic potential in pre-stasis basal epithelial populations, which produced tumors with dramatically reduced metaplastic differentiation. Our findings suggest that in a perturbed microenvironment, the adaptability of the EMT program, favoring intermediate states over complete mesenchymal differentiation, is crucial for directing basal epithelial cells toward a progenitor state that can progress to MBC.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003eReagents were obtained from Thermo Fisher Scientific unless otherwise stated.\u003c/p\u003e \u003cp\u003e \u003cstrong\u003eTissue Dissociation\u003c/strong\u003e \u003cp\u003eDisease-free breast tissue from women undergoing reduction mammoplasty was provided by the Cooperative Human Tissue Network (CHTN) Western Division, Nashville, TN, and Kaiser Foundation Research Institute, Oakland, CA. Human breast tissue was minced and enzymatically dissociated in RPMI with 2.1 mM L-glutamine and 25 mM HEPES supplemented with 10% (v/v) FBS (Atlanta Biologicals), 10 U/mL penicillin, 10 \u0026micro;g/mL streptomycin, 2.5 \u0026micro;g/mL Amphotericin B, 50 \u0026micro;g/mL gentamycin, Collagenase Type-2 (Worthington), and 100 U/mL hyaluronidase (Sigma-Aldrich) at 37\u0026deg;C for 12\u0026ndash;16 hours. The digest was centrifuged at 300 g for 10 minutes and washed with RPMI-1640 supplemented with 10% FBS. Epithelial-enriched structures (referred to as tissue organoids) were recovered by filtration through a 150 \u0026micro;m nylon mesh filter and frozen for long-term storage in Dulbecco's Modified Eagle Medium (DMEM) containing 50% (v/v) FBS and 10% (v/v) DMSO. The age and race of the donors are listed in Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003eFluorescence-Activated Cell Sorting (FACS) and Flow Cytometry\u003c/span\u003e: Tissue organoids were dissociated into single cells using 5 U/mL Dispase (Stem Cell Technologies) containing 1 \u0026micro;g/mL DNAse I (Stem Cell Technologies) for 5\u0026ndash;10 minutes and then 0.05% trypsin with 0.14 g/L Ethylenediaminetetraacetic Acid (EDTA) for 5\u0026ndash;10 minutes. Cultured cells were collected from the tissue culture plates by incubating them for 10\u0026ndash;15 minutes with TrypLE. The cells were washed with Phosphate-Buffered Saline (PBS)\u0026thinsp;+\u0026thinsp;1% Bovine Serum Albumin (BSA), filtered through a 40-\u0026micro;m cell strainer (Falcon), resuspended at a concentration of 1 \u0026times; 10⁶ cells in 100 \u0026micro;L of PBS\u0026thinsp;+\u0026thinsp;1% BSA, and stained for 15 minutes at room temperature with the appropriate antibody combination. When necessary, the cells were washed with PBS\u0026thinsp;+\u0026thinsp;1% BSA and stained for 15 minutes at room temperature with 1 \u0026micro;L streptavidin conjugated with BV785 (BioLegend) and 0.5 \u0026micro;g/mL DAPI (Sigma). The cells were washed, resuspended at 1 \u0026times; 10⁶ cells in 100 \u0026micro;L PBS\u0026thinsp;+\u0026thinsp;1% BSA, filtered through a 40-\u0026micro;m cell strainer, and loaded onto a FACS Aria II Cell Sorter (BD) fitted with a 100 \u0026micro;m nozzle. To isolate basal epithelial cells from tissue organoids, we immunostained cells with 5 \u0026micro;L of anti-CD10 mouse monoclonal antibody (mAb) (clone: HI10a) conjugated with PE (Biolegend), 1.5 \u0026micro;L of anti-Epcam mouse mAb (clone VU1D9) conjugated with FITC (Stem Cell Technologies), 4 \u0026micro;L of anti-CD49f rat mAb (clone GoH3) conjugated with APC (Biolegend), 8 \u0026micro;L of anti-CD2 mouse mAb (clone RPA-2.10) conjugated with Biotin (Becton Dickinson [BD]), 8 \u0026micro;L of anti-CD3 mouse mAb (clone HIT3a) conjugated with Biotin (BD), 8 \u0026micro;L of anti-CD16 mouse mAb (clone 3G8) conjugated with Biotin (BD), 8 \u0026micro;L of anti-CD64 mouse mAb (clone 10.1) conjugated with Biotin (BD), 4 \u0026micro;L of anti-CD31 mouse mAb (clone MBC78.2) conjugated with Biotin (ThermoFisher Scientific), and 1 \u0026micro;L of anti-CD45 mouse mAb (clone HI30) conjugated with Biotin (Biolegend). To isolate EMT phenotypes from cultured basal epithelial cells, we used 15 \u0026micro;L of anti-P-cadherin mouse mAb (clone CSTEM29) conjugated with APC (Thermo Fisher Scientific) and 4 \u0026micro;L of anti-N-cadherin mouse mAb (clone 8C11) conjugated with PE (BioLegend). After gating out debris, doublets, dead cells (DAPI+), and, if applicable, unwanted cell types (BV785+), the remaining target populations expressing the desired immunoprofiles were collected into a 1.5 mL tube containing 0.5 mL cell culture medium. The HMECs were pulsed with 10 \u0026micro;M EdU for one hour before the proliferation index was determined. The Click-iT Plus EdU Alexa Fluor 488 Flow Cytometry Assay Kit and FxCycle PI/RNase Staining Solution (Thermo Fisher Scientific) were used according to the manufacturer\u0026rsquo;s instructions. The cells were analyzed using an LSRFortessa flow cytometer (BD Biosciences).\u003c/p\u003e \u003cp\u003e \u003cstrong\u003eCell Culture\u003c/strong\u003e \u003cp\u003eDissociated human breast tissue organoids were plated in T-75 flasks (Corning) in Mammary Epithelial Cell Growth Medium (MEGM, Lonza), which includes hydrocortisone, human epidermal growth factor, insulin, bovine pituitary extract (BPE), gentamicin, and amphotericin B. After the initial plating, differential trypsinization was performed to selectively detach and aspirate unwanted fibroblasts, followed by collection of epithelial cells for further passaging. Sorted populations were cultured in one well of a six-well plate, transferred to a p100 plate, and then to a T-75 flask. Subculturing was performed at 70\u0026ndash;80% confluence using TrypLE Express for cell dissociation. Cells were incubated at 37\u0026deg;C in a humidified atmosphere containing 5% CO₂. The medium was changed every two to three days. Cell viability and counts were assessed using trypan blue exclusion and a Countess automated cell counter (Thermo Fisher Scientific) before reseeding or experimental use. The number of population doubling (PD) per passage was determined using the equation PD\u0026thinsp;=\u0026thinsp;log [A/B]/log2, where A is the number of collected cells and B is the number of plated cells. Attachment efficiency was consistently greater than 95% during the exponential growth phase. Cell populations were frozen in cryopreservation medium composed of 90% FBS and 10% DMSO (Sigma-Aldrich), aliquoted into cryovials, and gradually cooled to -80\u0026deg;C before being transferred to liquid nitrogen for long-term storage. The TGF-β inhibitor A83-01 (Tocris) was dissolved in DMSO at 1 mM and used at a concentration of 0.5 \u0026micro;M; an equivalent amount of DMSO was added to control cultures.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cstrong\u003eMammosphere Culture\u003c/strong\u003e \u003cp\u003eSingle cells were directly sorted into 24-well ultra-low attachment plates (Corning) at a density of 10,000 cells per well in serum-free mammary epithelial basal medium (MEBM) (Lonza) supplemented with B27 (Thermo Fisher Scientific), 20 ng/mL Epidermal Growth Factor (EGF, Peprotech), 20 ng/mL Fibroblast Growth Factor 2 (FGF2, Peprotech), and 4 \u0026micro;g/mL heparin (Stem Cell Technologies). Mammospheres were counted and collected by centrifugation at 100 \u0026times; g after 10 days of culture, dissociated enzymatically for 5\u0026ndash;10 minutes in TrypLE Express, filtered through a 40-\u0026micro;m cell strainer, and cultured again in suspension [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e].\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003eLentiviral Transduction\u003c/span\u003e: Basal epithelial cell populations, 1 \u0026times; 10⁶ in a single well of a six-well plate, were incubated with concentrated lentiviral particles at a multiplicity of infection of three with 8 \u0026micro;g/mL polybrene overnight for 12\u0026ndash;16 hours. For transformation, the cells were infected sequentially with (1) TERT with a hygromycin selection marker (Genecopoeia), (2) SV40 small T and large T antigens with a puromycin selection marker (Genecopoeia), and (3) KRAS G12V with a blasticidin selection marker (Amsbio). Cells were selected using 200 \u0026micro;g/mL hygromycin (Sigma-Aldrich), 1 \u0026micro;g/mL puromycin (InvivoGen), and 10 \u0026micro;g/mL blasticidin (InvivoGen). Transformed cells were grown in DMEM/F12 (1:1) supplemented with 5% Newborn Calf Serum (Hyclone), 10 \u0026micro;g/mL insulin, 10 ng/mL EGF, and 1 \u0026micro;g/mL hydrocortisone. CRISPR-mediated knockout of core EMT transcription factors, Snail (SNAI1) and Slug (SNAI2), Zeb (ZEB1), and Twist (TWIST1) was carried out using transEDIT CRISPR gRNA target gene sets from Transomics. Three constructs were tested for each gene. To suppress the EMT program, we chose to continue with: SNAI1 TEVH-1091413-pCLIP-All-EFS-ZsGreen, SNAI2 TEVH-1091323-pCLIP-All-EFS-ZsGreen, and non-targeting control TELA1013-CRISPR-NT#1-pCLIP-All-EFS-ZsGreen. CRISPR-mediated knockout of EZH2 was carried out using the Edit-R Human hEF1a-EGFP All-in-one Lentiviral sgRNA construct (Horizon Discovery Biosciences): EZH2 VSGH12180-256507490 and non-targeting VSGC12063.\u003c/p\u003e \u003cp\u003e \u003cstrong\u003eSoft-Agar Assay\u003c/strong\u003e \u003cp\u003eA solution of 1% agar (Sigma-Aldrich) in PBS was prepared, and 1.5 mL was added to each well of a six-well plate and allowed to gel at room temperature for 30 minutes as the bottom layer. Cells were prepared by suspending 1 \u0026times; 10⁵ cells in 3 mL pre-warmed DMEM/F12 containing 20% (v/v) FBS, 10 U/mL penicillin, and 10 \u0026micro;g/mL streptomycin. The cell suspension was gently mixed with 2 mL of 1% agar in PBS to achieve a cell density of 20,000 cells/mL in 0.4% agar; 1 mL of this mixture was added to each well of the six-well plate and kept at 4\u0026deg;C for 10 minutes to allow quick gelling, followed by adding 0.5 mL culture medium on top of the agar gel. The cells were incubated at 37\u0026deg;C with 5% CO₂, and the medium was changed every four days for four weeks. Colonies were stained with 0.5 mL of 0.005% Crystal Violet (Sigma-Aldrich) in PBS containing 4% formaldehyde for one hour and then rinsed with distilled water before counting the visible colonies.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cstrong\u003eXenograft Study\u003c/strong\u003e \u003cp\u003eExponentially growing cultures of oncogenically transformed cells (4 \u0026times; 10⁶) were injected subcutaneously into 8\u0026ndash;12-week-old NSG female mice in 50% Matrigel. Tumors were allowed to grow for 3\u0026ndash;6 months. After the study, xenograft tumors were weighed, formalin-fixed, and paraffin-embedded at the UCSF Histology \u0026amp; Biomarker Core. The maximum tumor burden was not reached in any of the mice in this study.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cstrong\u003eSandwich ELISA\u003c/strong\u003e \u003cp\u003eAnalysis of conditioned media from cultured HMECs from three donors at P2 (pre-stasis), P4 (stasis), and P8 (post-stasis) was initially performed using the Raybiotech Quantibody Human Kiloplex array, which detects 1,000 human biomarkers, as a service by the manufacturer. We followed up on interesting candidates using sandwich ELISA kits from Raybiotech to detect TGF-β, Activin A, Serpin E1, Wnt4, R-Spondin2, DKK1, and Follistatin according to the manufacturer\u0026rsquo;s protocols. For these assays, conditioned media were generated by culturing 300,000 cells per well in a six-well plate in 1 mL of serum-free basal media for 24 h. Each analyte was measured in eight biological and three technical replicates. The concentration of each factor in the conditioned media was extrapolated from the corresponding standard curves.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003eImmunofluorescence\u003c/span\u003e: Cells were cultured on four-well glass chamber slides (Millipore). The cells were fixed with 4% paraformaldehyde in PBS for 15 minutes at room temperature and permeabilized for 10 minutes in PBS containing 0.25% Triton X-100. Non-specific binding was blocked with 10% donkey serum (Jackson Immunoresearch), 0.3M glycine, and 0.1% Tween-20 in PBS for one hour. Primary antibodies were applied overnight in 1% donkey serum and 0.1% Tween-20 in PBS. The following primary antibodies were used: p16 (E6H4, Roche) mouse mAb (pre-diluted), p21 Waf1/Cip1 (12D1, Cell Signaling Technologies) rabbit mAb (1:400), p63-α (D2K8X, Cell Signaling Technologies) XP rabbit mAb (1:200), α-tubulin (DM1A, Sigma-Aldrich) mouse mAb (1:4000), MUC1 (HMFG2, Abcam) mouse mAb (1:200), Lamin A\u0026thinsp;+\u0026thinsp;Lamin C (EPR4100, Abcam) rabbit mAb (1:500), cytokeratin 19 (EPR1579Y, Abcam), rabbit mAb (1:400), and cytokeratin 14 (LL002, Abcam) rabbit mAb (1:200). Alexa Fluor 555-conjugated phalloidin (1:400) was used to stain F-actin (Thermo Fisher Scientific). Secondary donkey anti-mouse 488 and anti-rabbit 555 antibodies (diluted 1:500) were added to 1% donkey serum and 0.1% Tween-20 in PBS for one hour. DAPI (0.2 \u0026micro;g/mL in PBS) was added for five minutes at room temperature. The slides were washed three times for three minutes between each step in PBS containing 0.1% Tween-20. The coverslips were mounted with Vectashield HardSet Mounting Medium (Vector Laboratories). Images were acquired using a Leica SP8 laser scanning confocal microscope (only Tubulin/F-actin and Lamin A/C) at a resolution of 1,024 \u0026times; 1,024, processed using LASX software (Leica Microsystems) or a BZ-X800 fluorescence microscope (Keyence), and processed using ImageJ (version 2.14.0/1.54f).\u003c/p\u003e \u003cp\u003e \u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003eQuantitative Polymerase Chain Reaction (qPCR)\u003c/span\u003e: Total RNA was isolated from cells and treated with DNase I using the RNeasy Mini Kit (Qiagen). cDNA synthesis was performed using the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific). Quantitative PCR was performed on a CFX-96 (Bio-Rad Laboratories) thermocycler using 2x SsoFast Master Mix (Bio-Rad Laboratories) and analyzed using the standard curve method. A standard curve was prepared using the cDNA produced from Human Reference RNA (Agilent Technologies). Pre-made TaqMan Gene Expression Assays (Thermo Fisher Scientific) were used: ZEB1: Hs00232783_m1, SNAI1: Hs00195591_m1, SNAI2: Hs00161904_m1, TWIST1: Hs01675818_s1, CDH3: Hs00999915_m1, VIM: Hs00958111_m1, CDH1: Hs01023895_m1, CDH2: Hs00983056_m1, CDKN1A: Hs00355782_m1, MMP-2: Hs01548727_m1, FN1: Hs01549976_m1, CHRDL2: Hs01060234_m1, GREM1: Hs01879841_s1, and WNT5a: Hs00998537_m1. Custom primer-probe sets were used for GUSB: Forward Primer: CTCATTTGGAATTTTGCCGATT, Reverse Primer: CCGAGGAAGATCCCCTTTTTA, and Probe: FAM-TGAACAGTCACCGACGAGAGTGCTGGTA-TAM and CDKN2A (specific for p16\u003csup\u003eINK4A\u003c/sup\u003e and not p15\u003csup\u003eARF\u003c/sup\u003e): Forward Primer: CCAACGCACCGAATAGTTACG, Reverse Primer: GAGTGGCGGAGCTGCT, and Probe: FAM-CCGATCCAGGTCATGATG-TAM produced by Integrated DNA Technologies. GUSB expression was used to normalize variance in the input cDNA.\u003c/p\u003e \u003cp\u003e \u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003eWestern Blotting\u003c/span\u003e: Cells were washed with PBS and lysed using radioimmunoprecipitation assay (RIPA) buffer containing the HALT protease and phosphatase inhibitor cocktail. For each sample, a cell extract (corresponding to 50 \u0026micro;g of protein as determined by a Micro BCA Protein Assay Kit) was prepared in NuPAGE LDS Sample Buffer (4X) with Reducing Agent (10X), heated at 95\u0026deg;C for 15 minutes, and electrophoresed in each lane of a NuPAGE 4\u0026ndash;12% Bis-Tris gel along with a BenchMark Pre-stained Protein Ladder (Thermo Fisher Scientific). Proteins were transferred onto a nitrocellulose membrane (Biorad) overnight at 4\u0026deg;C at 36 mV. The membrane was washed in TBS-T (0.05 M Tris-HCl pH 7.5, 0.15 M NaCl, 0.05% Tween-20), blocked for one hour at room temperature in blotto (5% nonfat dry milk in TBS-T), incubated with primary antibodies overnight at 4\u0026deg;C, washed in TBS-T, incubated with goat anti-rabbit horseradish peroxidase conjugate secondary antibody (Jackson Immunoresearch) 1:5000 in Blotto for one hour, washed in TBS-T, and developed with ECL Western Blotting Substrate. Chemiluminescent signals were detected by exposing the CL-XPosure Film to the membranes or using the KwikQuant Digital Western Blot Detection System. Primary antibodies used at a 1:1000 dilution from Cell Signaling Technologies were Slug (C19G7) rabbit mAb, ZEB1 (D80D3) rabbit mAb, Slug (C19G7) rabbit mAb, E-cadherin (24E10) rabbit mAb, N-cadherin (D4R1H) XP rabbit mAb, Vimentin (D21H3) XP rabbit mAb, RAS rabbit pAb, EZH2 (D2C9) XP rabbit mAb, Histone H3 (D1H2) XP rabbit mAb, and β-Actin (13E5) rabbit mAb; and from Active Motif Histone H3K27me3 (MABI 0323) mouse mAb. Full images of the western blots are shown in Supplemental Figs. S4 and S5.\u003c/p\u003e \u003cp\u003e \u003cspan type=\"ItalicUnderline\" class=\"ItalicUnderline\" name=\"Emphasis\"\u003eMultiplex Immunohistochemistry\u003c/span\u003e: Sections (5 \u0026micro;m thick) were cut from FFPE tissue blocks and placed on positively charged Superfrost Plus microscopy slides. Slides were baked at 60\u0026deg;C overnight, deparaffinized in xylene, and rehydrated in graded ethanol (100%, 100%, 95%, 85%, and 70%) in distilled H₂O. Endogenous peroxidases were quenched with 3% H₂O₂ (Sigma-Aldrich) diluted in phosphate-buffered saline (PBS). Heat-induced antigen retrieval was performed in citrate buffer pH 6.0 (Sigma-Aldrich) at 95\u0026deg;C for 15 minutes. Non-specific antibody binding was blocked using a Background Sniper (Biocare Medical). The tissue sections were incubated for one hour at room temperature (RT) with each primary antibody in 1% BSA and 30 minutes in pre-diluted MACH 2 conjugated anti-mouse or anti-rabbit secondary antibodies (Biocare Medical). The slides were washed in TNT buffer (0.1 M TRIS-HCL pH 7.5, 0.15 M NaCl, and 0.05% Tween-20) following blocking and three times for five minutes each after applying both primary and secondary antibodies. The signal was developed using a Tyramide Signal Amplification (TSA) solution (Akoya): FITC (two minutes), Cy3 (three minutes), or Cy5 (seven minutes). The antibody complex was removed by heating at 95\u0026deg;C in citrate buffer pH 6.0 for five minutes to allow for multiplex staining. Nuclei were counterstained with 3 \u0026micro;M DAPI in PBS for five minutes, washed in distilled H₂O, and mounted using Vectashield HardSet Mounting Medium (Vector Laboratories). The primary antibodies used from Cell Signaling Technologies were Slug (C19G7) rabbit mAb (1:1000), Snail (C15D3) rabbit mAb (1:1000), p63-α (D2K8X) XP rabbit mAb (1:2000), and vimentin (D21H3) rabbit mAb (1:5000); from AbCam Cytokeratin 5 (SP27) rabbit mAb (1:4000), Lamin A\u0026thinsp;+\u0026thinsp;Lamin C (EPR4100) rabbit mAb (1:10000), Cytokeratin 13 (EPR3671) rabbit mAb (1:2000), Ki67 (SP6) rabbit mAb (1:6000), and Cytokeratin 14 (SP53) rabbit mAb (1:3000).The slides were imaged using a BZ-X800 analyzer. Images were prepared for analysis using the ImageJ software (version 2.14.0/1.54f). Analysis was performed using the Qupath software (v0.5.0)[\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. Using the pre-trained models, nuclear segmentation was performed using StarDist[\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e], normalizePercentiles\u0026thinsp;=\u0026thinsp;1, 99, threshold\u0026thinsp;=\u0026thinsp;0.5, pixel size\u0026thinsp;=\u0026thinsp;0.5, and cell expansion\u0026thinsp;=\u0026thinsp;5.0. Individual classifiers were trained for each marker and combined to create a composite classifier for scoring the cells within the annotated regions.\u003c/p\u003e \u003cp\u003e \u003cstrong\u003eStatistical Analysis\u003c/strong\u003e \u003cp\u003eFor each graph, the mean and standard deviation are shown in red. Unless otherwise noted in the figure legend, each black circle indicates a single donor. Pairwise comparisons were performed using an unpaired Welch\u0026rsquo;s t-test (GraphPad Prism Software). Levels of significance used were *0.05 to 0.01, **0.01 to 0.001, ***0.001 to 0.0001, and ***\u0026lt;0.0001.\u003c/p\u003e \u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e \u003cem\u003eRobust Activation of the TGF-β Pathway and Induction of EMT in Stasis HMECs Precede Variant Emergence\u003c/em\u003e \u003c/p\u003e \u003cp\u003eWhen cultured in serum-free Mammary Epithelial Cell Growth Medium (MEGM), HMECs isolated from three tissue donors experienced growth cessation after 5\u0026ndash;10 population doublings. Stasis lasted between 20 and 40 days before the emergence of clonal variant populations, which re-initiated logarithmic growth (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA). Nearly complete proliferative arrest was observed at stasis (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB), marked by the upregulation of p16 and development of a senescence-like morphology (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC). These results are consistent with previously published studies that have shown the clonal emergence of rare variants with pre-malignant properties within a field of growth-arrested, senescent-like epithelial cells [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e, \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e, \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e, \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003ePrevious studies of this culture system have primarily focused on comparing the initial pre-selection HMEC population to post-selection variants. In contrast, we chose to characterize stasis cultures, specifically focusing on the factors that enforce growth arrest. Conditioned media (CM) from stasis cultures, diluted 1:1 with fresh MEGM, induced significant growth arrest in early passage HMECs compared to fresh MEGM or CM from pre-stasis HMECs at passage 2 (P2), mixed 1:1 with fresh MEGM (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD). Senescent and DNA double-strand breaks are associated with complex secretory programs, termed senescence-associated secretory phenotype (SASP) and stress-elicited extrinsic phenotype (SEEP), respectively, which can initiate or reinforce growth arrest in both autocrine and paracrine manners [\u003cspan additionalcitationids=\"CR50\" citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e]. To characterize alterations in the cellular secretory program during HMEC culture, we collected CM and extracted cellular RNA and protein from cultures of pre-stasis HMECs at P2, growth-arrested stasis HMECs (S), and variant HMECs at two passages beyond stasis (S\u0026thinsp;+\u0026thinsp;2). In stasis cultures, quantitative proteomic analysis using an antibody array (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e) indicated increased levels of transforming growth factor-beta (TGF-β) and Activin A, a member of the TGF-β superfamily. ELISA assays confirmed higher levels of TGF-β, Activin A, and Plasminogen Activator Inhibitor-1 (PAI-1), a biomarker of TGF-β pathway activation, in stasis cultures than in proliferating pre-stasis or variant HMEC cultures (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eE). Additionally, we detected an increase in matrix metalloproteinase 2 (MMP-2) and fibronectin (FN) transcript levels during stasis (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eF), which can be upregulated in response to the activation of the TGF-β pathway [\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e, \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. Several cytokines and chemokines typical of SASP (IL-6, IL-8) [\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e] remained unchanged in the stasis cultures (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003eA).\u003c/p\u003e \u003cp\u003eIn stasis HMECs, we observed significant changes in gene expression indicative of EMT activation, which is consistent with the fundamental role of the TGF-β pathway in this process [\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. At the transcript level, we observed decreased levels of the epithelial markers E-cadherin (CDH1, concentrated in luminal epithelial cells) and P-cadherin (CDH3, typically expressed by basal epithelial cells), and no change in N-cadherin (CDH2, typically expressed by mesenchymal cells) in stasis HMECs compared to pre-stasis HMECs (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eG). We also observed increased levels of core EMT transcription factors Slug (SNAI2), Snail (SNAI1), Twist (TWIST1), and Zeb (ZEB1). Transcript levels of E-cadherin, Snail, Twist, and Zeb decreased in variant cultures compared to those in stasis cultures and were not significantly different from those observed in pre-stasis HMECs (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eG). In contrast, Slug transcript levels were increased in stasis HMECs and maintained in the variants. The derivation of variants from the basal compartment [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e] likely accounts for the increased levels of P-cadherin observed in variant cultures compared with both pre-stasis and stasis cultures (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eG). Western blot analyses were consistent with the increased EMT activation, particularly during stasis. We observed decreased E-cadherin expression and increased Snail, Zeb, N-cadherin, and Vimentin expressions. Notably, protein expression was inconsistent with the changes in certain transcript levels (N-cadherin, Snail, and Slug), suggesting regulation by post-transcriptional processes commonly observed in EMT (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eH). In addition to TGF-β (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD), several paracrine and autocrine factors involved in EMT [\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e] were robustly increased in stasis HMECs compared with those in either of the two active growth phases (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003eB and S1C).\u003c/p\u003e \u003cp\u003eIn summary, our initial observations demonstrated transient activation of the TGF-β pathway and engagement of EMT programs in cultured HMECs as they entered stasis. Remarkably, the emergence of variants from these cultures was associated with the attenuation of the TGF-β and EMT pathways.\u003c/p\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eBasal Epithelial Cells Exhibit Enhanced EMT Activation Compared to Luminal Epithelial Cells\u003c/h2\u003e \u003cp\u003eThe human breast epithelium is composed of two distinct layers: an inner luminal layer of polarized columnar epithelial cells and an outer basal layer of contractile myoepithelial cells. Both layers comprise a heterogeneous mixture of mature cells and cells at various stages of differentiation, including stem and progenitor cells [\u003cspan additionalcitationids=\"CR2 CR3 CR4 CR5\" citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. For simplicity, we will refer to these complex populations as 'luminal' and 'basal' cells, acknowledging that these terms encompass a range of cellular states and phenotypes. To reduce the confounding variable of comparing a mixture of luminal and basal derivatives in pre-stasis cultures with exclusive basal derivatives in variant cultures [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e], we used fluorescence-activated cell sorting (FACS) to separate basal epithelial cells (EPCAM\u003csup\u003e\u0026minus;\u003c/sup\u003eCD49f\u003csup\u003e+\u003c/sup\u003eCD10\u003csup\u003e+\u003c/sup\u003e) from luminal epithelial cells (EPCAM\u003csup\u003e+\u003c/sup\u003e) (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA). When cultured in serum-free MEGM, the basal-enriched population from the three tissue donors entered stasis between 9 and 15 population doublings. These cultures consistently produced variant outgrowths 30 days after the cessation of growth at P5 (passage 5). In contrast, luminal cultures proliferated between 8 and 11 population doublings before growth arrest at P3, and did not produce variant colonies (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eB), consistent with a previous report [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. Basal epithelial cells contributed significantly greater levels of TGF-β but not Activin A compared to luminal epithelial cells when analyzed at stasis, P3 for isolated luminal cells, and P5 for isolated basal epithelial cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eAs early as P2, even before entering stasis, basal-enriched populations demonstrated a more dispersed and migratory phenotype and incorporated vimentin, a canonical marker of the mesenchymal phenotype, into their cytoskeleton, whereas luminal cells maintained tight colony morphologies and did not express significant levels of vimentin (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eD). At the transcriptional level, the basal epithelial population demonstrated a robust EMT response, as evidenced by significant increases in Twist, Zeb, and N-cadherin, and decreases in E-cadherin and P-cadherin between P2 and P5. Luminal epithelial cells demonstrated a significant increase in Slug and N-cadherin expression, and a decrease in E-cadherin expression (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eE). Thus, both luminal and basal epithelial cells demonstrated evidence of EMT; however, the magnitude of EMT-related changes in gene expression was markedly higher in the basal epithelial cells.\u003c/p\u003e \u003cp\u003eSimilar to our observations in bulk HMECs, basal epithelial cell cultures were characterized by high levels of TGF-β and Activin A at early passages (P3 and P5), followed by a significant reduction in expression of these cytokines at later passages (P8) following the emergence of variants (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eF). Western blot analysis revealed high levels of N-cadherin, Snail, and Zeb in early passage (P3) and stasis (P5) basal epithelial cells, indicating activation of EMT. However, these components of the EMT program were downregulated in the variant cultures (P8). Notably, among the EMT transcription factors analyzed, Slug levels were consistently high across all passages, including the variants (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eG). A role of Slug in maintaining basal epithelial gene expression has been observed previously [\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e]. In summary, only the basal epithelial population contains variants [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. This population recapitulated the transient activation of TGF-β and EMT phenotypes observed in bulk HMEC cultures, leading us to focus further on the basal epithelial cell population\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eInhibition of TGF-β Signaling and EMT Transcription Factors Prevents Stasis in Basal Epithelial Cells\u003c/h2\u003e \u003cp\u003eThe application of A83-01, a potent inhibitor of TGF-β type I receptor ALK5 kinase, type I activin/nodal receptor ALK4, and type I nodal receptor ALK7, was sufficient to prevent stasis in basal epithelial cells (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA). A83-01 inhibited EMT at the transcriptional level, as evidenced by the significantly increased expression of P-cadherin and the decreased expression of N-cadherin, Snail, Twist, and Zeb (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB). TGF-β inhibition was associated with a highly significant decrease in cell size (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eC), reduction in nuclear area, and elimination of micronuclei (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eD). A83-01 treatment also significantly reduced p21 and p16 levels (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eE). We observed that treatment with A83-01 maintained the protein expression of p63, an important transcriptional regulator of basal cell differentiation (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eF). Together, these results implicate the engagement of the TGF-β pathway in cellular dysfunction and the activation of tumor suppressor pathways during stasis.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eSimilar to the treatment of basal epithelial cells with A83-01, we found that CRISPR-mediated knockout (KO) of Snail or Slug eliminated the stasis barrier compared with non-targeting controls (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). Following selection, Snail and Slug KO basal epithelial cells generated tight colony morphologies similar to those produced by the variant basal epithelial cultures emerging from stasis, whereas non-targeting controls were more dispersed across the culture surface. In contrast, Zeb and Twist KO basal epithelial cells demonstrated a growth pattern similar to that of the non-targeting control, and were therefore not followed up in culture (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). Targeting Snail and Slug using CRISPR vectors reduced their transcript levels. Consistent with their role in the EMT network, Snail and Slug knockout significantly increased the transcript levels of P-cadherin and decreased those of N-cadherin, Twist, and Zeb compared to non-targeting controls (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). Basal epithelial cells with Snail or Slug KO demonstrated substantially lower transcription of p21 and p16 tumor suppressors than the controls (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC). We observed vimentin in the cytoskeleton of basal epithelial cells expressing p16, consistent with the idea that EMT is correlated with the senescence phenotype in this culture system (Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003eA). Notably, the reduced growth arrest observed upon Snail and Slug KO was associated with decreased levels of TGF-β and the TGF-β superfamily member activin A (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eD), which are part of the autocrine and paracrine loops that promote the EMT program[\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e] and perpetuate senescence[\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn contrast, overexpression (OE) of the EMT-related transcription factors Snail and Slug immediately suppressed the growth of basal epithelial cells, preventing their expansion (Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003eB). In this experiment, Snail and Slug were expressed at non-physiological levels, which were far higher than those observed in basal epithelial cells at stasis. Analysis of these cultures at the first passage showed a high degree of EMT, as evidenced by the downregulation of P-cadherin and upregulation of N-cadherin, Twist, and Zeb (Fig. \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003eC). Further examination of this model was not feasible because of the rapid and complete cessation of cell growth.\u003c/p\u003e \u003cp\u003eIn summary, TGF-β pathway inhibition or Snail or Slug knockdown prevented EMT activation and stasis in basal epithelial cells. Activation of the TGF-β pathway induces the expression of EMT-related transcription factors, such as Snail, Slug, Zeb, and TWIST, which promote the transition from an epithelial to a mesenchymal phenotype. In turn, mesenchymal cells resulting from EMT secrete additional TGF-β and Activin A creating a positive feedback loop that sustains and amplifies the EMT process. The activation of TGF-β and EMT in basal epithelial cells leads to growth suppression and cellular dysfunction, as is evident in conventional culture systems (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eE).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eEMT Generates a Spectrum of Basal Epithelial Cell States with Varying Proliferative and Self-Renewal Capacities\u003c/h2\u003e \u003cp\u003eLimiting Activation of the TGF-β pathway and EMT induction was sufficient to avoid stasis in the entire cultured basal epithelial cell population. However, rare basal epithelial variants were intrinsically capable of avoiding stasis. Variants can be visualized as epithelial colonies (tightly clustered and adhered to their neighbors) emerging among cells with spindle or senescent-like morphologies (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA). EMT is a dynamic process that enables epithelial cells to acquire mesenchymal properties, leading to increased plasticity and the ability to transition between multiple cellular states [\u003cspan additionalcitationids=\"CR22 CR23 CR24 CR25\" citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. This process creates a spectrum of intermediate phenotypes that allow cells to exhibit varying degrees of epithelial and mesenchymal characteristics.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eMany studies have used CD24 and CD44 to distinguish between epithelial and mesenchymal states [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. Indeed, almost all cells in stasis cultures expressed the CD44\u003csup\u003ehi\u003c/sup\u003eCD24\u003csup\u003elo\u003c/sup\u003e phenotype, consistent with the high levels of EMT at this time point. In contrast, pre-stasis HMECs and post-stasis variants mostly lacked CD44 (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eB). However, this method was unable to separate epithelial, intermediate, and more complete mesenchymal phenotypes. To better distinguish between these EMT states in basal epithelial cell cultures, we developed an alternative method that utilized antibodies against P-cadherin, an epithelial marker highly expressed by basal epithelial cells, and N-cadherin, a marker highly expressed in mesenchymal cell types. In early and stasis cultures, we observed complete EMT states expressing N-cadherin but lacking P-cadherin (P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e), an intermediate EMT state expressing both P- and N-cadherins (P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e), and an epithelial state expressing P-cadherin but lacking N-cadherin (P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e\u0026minus;\u003c/sup\u003e). In cultures dominated by variants, the complete EMT state, P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e, was mostly absent (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eC). We sorted basal epithelial cells in the P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e and P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e states from stasis cultures (P5). At the mRNA level, P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e cells exhibited significantly lower levels of P-cadherin and higher levels of N-cadherin, Snail, Slug, Twist, and Zeb than did P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e cells, confirming that P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e populations represent a higher degree of EMT (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eD). Notably, basal epithelial cells with a complete EMT phenotype were growth arrested and could not be passaged, although they remained adherent to the culture surface for several weeks (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eE).\u003c/p\u003e \u003cp\u003eEarly passage basal epithelial cells were sorted based on the immunophenotypes P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e\u0026minus;\u003c/sup\u003e, P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e, and P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e) corresponding to different states along the EMT spectrum, and cultured for 14 days. When the resulting populations were re-analyzed, both P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e\u0026minus;\u003c/sup\u003e (blue) and P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e (purple) cells were able to recapitulate the entire spectrum of EMT states, whereas P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e (red) cells were restricted to the \u0026lsquo;fully\u0026rsquo; mesenchymal state (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eF). When separated from stasis cultures before the emergence of variants, we demonstrated, using limiting dilution analysis, that basal epithelial cells in the P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e\u0026minus;\u003c/sup\u003e and P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e states could give rise to variant colonies, whereas cells in the P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e state could not (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eG). In addition to the ability to form variant colonies in monolayer culture, we observed that basal epithelial cells in the P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e\u0026minus;\u003c/sup\u003e and P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e states could form mammospheres under non-adherent conditions, which is suggestive of self-renewal and stem cell properties [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e]. In contrast, cells in the P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e state were unable to create mammospheres (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eH). The P\u003csup\u003e\u0026minus;\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e state was completely absent in the variant cultures, and cells could only move between the P\u003csup\u003e+\u003c/sup\u003eN \u003csup\u003e\u0026minus;\u003c/sup\u003e and P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e states, suggesting attenuation of the EMT program (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eI).\u003c/p\u003e \u003cp\u003eIn summary, isolated basal epithelial cells initially activated an EMT program that could progress to a fully mesenchymal state, ultimately leading to growth arrest. However, over several passages, we observed stabilization of the epithelial and intermediate EMT states, restricting access to the fully mesenchymal state. This restriction was associated with the emergence of variants that retained plasticity and could transition between epithelial and intermediate EMT states but not the fully mesenchymal state.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003ePRC2-Mediated Epigenetic Regulation Governs Access to EMT States in Basal Epithelial Variants\u003c/h2\u003e \u003cp\u003ePRC2 is known to play a critical role in epigenetic regulation of gene expression in variant basal epithelial cells. Indeed, by catalyzing the trimethylation of histone H3 on lysine 27 (H3K27me3), PRC2 silences the target genes involved in growth and differentiation. Over time in culture, silencing of many PRC2 targets can be reinforced by promoter hypermethylation [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. EZH2, the catalytic subunit of PRC2, is central to this function, and its knockout disrupts the ability of this complex to regulate gene expression. We found that the knockout of EZH2 in basal epithelial cells at P3 reliably prevented the formation of variant colonies, which was expected given that PRC2 plays a well-established role in silencing CDKN2A and removing other obstacles to proliferation (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eA). EZH2 knockout in the P8 variant basal epithelial cells also impaired cell proliferation (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eB).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eRecent studies have revealed that PRC2 plays a crucial role in controlling the EMT trajectory by repressing specific genes involved in EMT [\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. In EZH2 knockout basal epithelial variants, we confirmed a decrease in EZH2 expression and global H3K27me3 levels at the protein level, and observed increased levels of Zeb and N-cadherin (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eC). EZH2 knockout promoted EMT at the transcriptional level, as evidenced by decreased P-cadherin, but increased N-cadherin, Snail, Twist, and Zeb expression (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eD).\u003c/p\u003e \u003cp\u003eBy reducing the activation of the EMT program, PRC2 and potentially other epigenetic regulators maintain variant basal epithelial cells in an intermediate EMT state that favors continued proliferation and stemness, and disfavors the path to complete EMT states associated with growth arrest (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eE).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eTGF-β Inhibition Prior to Cellular Transformation Reduced Metaplastic Phenotypes in Basal Epithelial-Derived Tumors\u003c/h2\u003e \u003cp\u003eBasal epithelial progenitors have been proposed as the cells of origin for MBC [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e], a rare histological subtype associated with elevated levels of EMT compared with other breast cancer subtypes [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e] (Fig. S3A and S3B). In culture, basal epithelial cells strongly activated the TGF-β pathway and induced EMT. The selection of variants was associated with the stabilization of intermediate EMT states capable of continued proliferation, limiting access to a complete mesenchymal state subject to growth arrest. Pharmacological inhibition of the TGF-β pathway suppressed the activation of the EMT program and eliminated the selective pressure that promoted the outgrowth of the basal variants (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eA).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTo better understand the role of TGF-β and EMT-dependent selection of the variant phenotype in MBC carcinogenesis, we employed an established transformation protocol [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e] involving the sequential transduction of TERT, SV40 T antigen, and mutant K-Ras\u003csup\u003eG12V\u003c/sup\u003e. We compared the oncogenic transformation of P3 basal epithelial cell populations cultured in MEGM for three passages, poised to enter stasis, with variant basal epithelial cell populations cultured in MEGM for eight passages that had already overcome stasis. Additionally, we compared P3 basal epithelial cell populations cultured in MEGM\u0026thinsp;+\u0026thinsp;DMSO with those cultured in MEGM\u0026thinsp;+\u0026thinsp;A83-01 for three passages, which suppressed EMT and allowed them to avoid stasis and variant selection. Notably, both the oncogenically transformed variant and A83-01 basal epithelial populations expressed levels of the mutant K-Ras\u003csup\u003eG12V\u003c/sup\u003e equivalent to those of the P3 controls (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eB). Proliferation of oncogenically transformed P3 basal epithelial cells was suppressed, whereas the oncogenically transformed variant and A83-01 basal epithelial populations exhibited exponential growth kinetics (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eC). Additionally, oncogenically transformed P3 basal epithelial cells formed dramatically fewer colonies in soft agar, with most replicates showing no colony formation at all, indicating a reduced capacity for anchorage-independent growth compared to the variant and A83-01 basal epithelial populations, which showed no significant difference in their ability to form colonies in soft agar (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eD). A greater level of EMT activation was observed in oncogenically transformed P3 basal epithelial cells compared to both the variant and A83-01 populations, as evidenced by the lower levels of P-cadherin and higher levels of N-cadherin, Snail, Twist, and Zeb compared to the variant and A83-01 basal epithelial cell populations. In contrast, transformed variant populations exhibited intermediate levels of EMT, characterized by significantly higher levels of P-cadherin and lower levels of N-cadherin, Snail, Twist, and Zeb, compared to A83-01-treated basal epithelial populations (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eE). Flow cytometric analysis of P-cadherin and N-cadherin expression confirmed that oncogenically transformed variant populations predominantly exhibited intermediate EMT phenotypes (P\u003csup\u003e+\u003c/sup\u003eN\u003csup\u003e+\u003c/sup\u003e). In contrast, oncogenically transformed basal epithelial cells treated with A83-01 mostly expressed P-cadherin, with relatively little N-cadherin expression, thus mostly preserving their epithelial state (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eF).\u003c/p\u003e \u003cp\u003eOnly the oncogeneically transformed variant and A83-01 basal epithelial cell populations formed tumors when xenotransplanted subcutaneously into NSG mice (Fig. S3C). Tumors derived from both A83-01-treated and variant basal epithelial cells were highly proliferative based on the number of cells that stained positive for Ki67 (Fig. S3D). Both sets of tumors expressed p63, which is a typical phenotype of MBCs, and several even rarer histological subtypes likely arising from basal epithelial progenitors [\u003cspan additionalcitationids=\"CR58\" citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e]. Notably, A83-01 tumors exhibited a lower frequency of cells expressing the squamous differentiation marker cytokeratin 13, the EMT markers Snail and Slug, and the mesenchymal differentiation marker vimentin when compared to variant tumors (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eG).\u003c/p\u003e \u003cp\u003eIn summary, high levels of EMT in the malignantly transformed P3 basal epithelial group were associated with poor growth and tumorigenicity, whereas the variant basal epithelial group, which exhibited intermediate levels of EMT gene expression, demonstrated enhanced growth and tumorigenicity. Inhibition of the TGF-β pathway, which significantly reduced EMT gene expression, promoted growth and tumorigenicity in the P3 basal epithelial group but led to a loss of metaplastic squamous and metaplastic elements.\u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eIn conventional monolayer cultures [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e], HMECs are deprived of instructive biophysical and biochemical cues that maintain cellular homeostasis and cooperation within tissues. Instead, they form attachments with a non-compliant plastic substratum and proliferate in response to signals presented in an engineered medium. Consequently, HMECs enter stasis, the first, \u0026lsquo;stress-associated\u0026rsquo; barrier to indefinite proliferation, during which most cells adopt a senescence-like phenotype [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. To address this limitation, more recently developed culture systems optimized to support epithelial proliferation incorporate fibroblast feeders or hydrogels containing basement membrane components, mimicking key aspects of the \u003cem\u003ein vivo\u003c/em\u003e microenvironment lost in conventional monolayer cultures. These systems also employ chemical inhibitors to reduce the ability of the cells to engage in stress response pathways. Collectively, these modifications allow cultured breast epithelial cells to avoid both stasis and agonescence, the second barrier to indefinite proliferation resulting from telomere attrition, without the need for genetic manipulation. Notably, these culture systems support genetic and epigenetic stability during long-term propagation [\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e, \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn such conventional monolayer cultures, rare variants activate an aberrant yet remarkably consistent epigenetic program that enables them to adapt to unfavorable culture conditions and reinitiate proliferation after encountering the stasis barrier [\u003cspan additionalcitationids=\"CR36\" citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e]. Many epigenetic alterations favored in variant cultures, including hypermethylation of the CDKN2A promoter, have also been observed during human tumor progression [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e]. In intact tissues, etiological factors, such as aging, chronic inflammation, and environmental exposure, gradually disrupt homeostatic interactions. This progressive disruption contributes to a decline in the quality of the tissue microenvironments, which is hypothesized to favor the selection of clones (epigenetic progenitors) that are adaptive to altered microenvironments, potentially leading to the emergence of cancer over the course of many decades [\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e, \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e]. Arguably, conventional monolayer cultures represent an analogous process, albeit with a significantly accelerated time scale. Thus, variants within these cultures have proven useful in elucidating the mechanisms underlying the emergence of cancer progenitor states [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e, \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn this study, we report that robust activation of the TGF-β pathway and induction of the EMT program are drivers of stasis and, therefore, constitute an important selective pressure underlying the emergence of the variant phenotype. Inhibition of the TGFβ pathway allowed basal epithelial cells to bypass stasis and dramatically reduced cell size, chromosomal instability, and expression of both p21 and p16. TGF-β plays a direct role in enhancing cell size by activating the mammalian target of rapamycin (mTOR) [\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e]. Notably, increased cell size is not simply a consequence of cellular senescence, but also a contributing factor through mechanisms that include failure to scale macromolecule production with cell volume causing cytoplasmic dilution, altered nuclear organization, chromosomal missegregation, DNA damage, and mitochondrial dysfunction [\u003cspan additionalcitationids=\"CR66 CR67 CR68\" citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e]. TGF-β induces cellular senescence through SMADs, which can upregulate p15 and p21 and contribute to the establishment of SASP, which can reinforce senescence in an autocrine manner or induce senescence in neighboring cells in a paracrine manner [\u003cspan additionalcitationids=\"CR50\" citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e]. Knockout of the key EMT factors Snail and Slug also allowed basal epithelial cells to bypass stasis. The EMT program is associated with increased activation of the TGFβ pathway through autocrine and paracrine secretion [\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e] but can also have negative consequences through other mechanisms. Indeed, during EMT, if cells continue to proliferate, they exhibit mitotic abnormalities such as centrosome clustering, cytokinesis failure, and chromosome missegregation, which result in increased DNA damage, aneuploidy, binucleated cells, and micronuclei. These mitotic defects are associated with the suppression of nuclear envelope proteins like LaminB1, critical for mitotic regulation [\u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e70\u003c/span\u003e, \u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e71\u003c/span\u003e]. Thus, TGFβ pathway activation and EMT are key components of the stasis barrier, resulting in cellular arrest and dysfunction.\u003c/p\u003e \u003cp\u003eNotably, EMT is not a binary switch between epithelial and mesenchymal states but generates a spectrum of intermediates with both epithelial and mesenchymal traits. When basal epithelial cells underwent a complete transition to the mesenchymal state, which peaked at stasis, they were subject to proliferative arrest. In contrast, the re-initiation of proliferation within stasis cultures by variants was associated with the stabilization of intermediate states. These hybrid states are characterized by enhanced cellular plasticity and adaptability, which can lead to increased stem cell properties and tumorigenic potential [\u003cspan additionalcitationids=\"CR22 CR23 CR24 CR25\" citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. PRC2 is known to contribute to the epigenetic program underlying the emergence of the variant phenotype by silencing genes involved in growth control and differentiation [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e, \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e]. PRC2 also suppresses the expression of EMT-inducing transcription factors, such as ZEB1 and ZEB2 [\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. Indeed, in variant basal epithelial cells, knockout of EZH2, the catalytic subunit of PRC2, increased EMT gene expression, influencing the trajectory of the EMT program. Significantly, epigenetic alterations in the TGF-β pathway have been previously observed in variants [\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e, \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e72\u003c/span\u003e]. EZH2 knockout in basal epithelial cells and variants impaired PRC2 function and diminished basal epithelial proliferation, limiting our study. In a prior study, fully transformed HMLER cells, which were derived from variant cultures, were not subject to growth arrest in response to compromised PRC2 activity; instead, they transitioned into a quasi-mesenchymal state associated with high metastatic potential. In contrast, the loss of another epigenetic regulator, KMT2D, results in a highly mesenchymal state [\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. Further research is warranted to elucidate the interplay between epigenetic regulators and EMT during the emergence of the variant phenotype. Nevertheless, our findings suggest that the adaptation of EMT toward favorable intermediate EMT states is a consequence of the epigenetic program, which facilitates the emergence of variants by modulating gene expression involved in growth control, differentiation, and other hallmark programs characteristic of a cancer progenitor state.\u003c/p\u003e \u003cp\u003eConsistent with the idea that variants are representative of a cancer progenitor state, the variant state confers tolerance to oncogenic insults such as RAS activation, thereby bypassing oncogene-induced senescence [\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. In this study, we also observed the increased vulnerability of variants to oncogenic transformation when compared to pre-stasis basal epithelial cells, which had not undergone the same epigenetic reprogramming characteristic of the variant population. TGF-β inhibition, which bypassed stasis in basal epithelial cultures, improved the transformation efficiency of pre-stasis basal epithelial cells, suggesting that overcoming the cellular arrest and dysfunction resulting from activation of the TGFβ pathway and EMT is necessary for achieving malignant transformation in these cells. This result was consistent with a previous report that found that either a dominant-negative TGFβ type II receptor or a TGFβ pathway inhibitor was able to bypass RAS-induced growth arrest variants and HMECs with compromised p16/Rb and p53. This study reported no difference in oncogenic transformation efficiency between pre-stasis HMEC and variant HMEC populations [\u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e73\u003c/span\u003e]. However, they compared pre-stasis HMEC (specimen 48R, batch T) with post-selection HMEC (specimen 48R, batch S), which were developed in different media systems\u0026mdash;serum-free MCDB170 medium (a non-commercially available version of MEGM) versus serum-containing M87A supplemented with cholera toxin and oxytocin [\u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e74\u003c/span\u003e], which differed from our analysis of basal epithelial cells cultured in MEGM. Strikingly, the oncogenic transformation of basal epithelial cells resulted in metaplastic tumors with squamous and mesenchymal components in immunocompromised mice, which is consistent with previous reports [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e, \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. However, treatment with TGF-β pathway inhibitors generated tumors with a marked reduction in metaplastic elements, even though the TGF-β pathway inhibitor was not maintained in culture after selection for KRAS\u003csup\u003eG12V\u003c/sup\u003e or in mice.\u003c/p\u003e \u003cp\u003eBasal cells adapt to conventional HMEC culture systems by restricting the EMT program to a hybrid epithelial/mesenchymal state that maintains a high proliferation, stemness, and tumorigenicity. This cellular plasticity is required to achieve a basal progenitor state that gives rise to rare MBCs. Notably, the MBC-like tumor phenotypes observed in the variant tumors show similarities to squamous cell carcinomas originating from basal epithelial cells in stratified epithelia throughout the body, such as lung, head and neck, esophageal, and cutaneous squamous cell carcinomas [\u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e]. A subset of these squamous tumors are associated with mesenchymal differentiation, akin to that sometimes observed in MBCs [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e], often referred to as spindle cells or sarcomatous elements [\u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e76\u003c/span\u003e]. Both MBC and squamous cell carcinoma are associated with inactivation of CDKN2A [\u003cspan additionalcitationids=\"CR78 CR79 CR80\" citationid=\"CR77\" class=\"CitationRef\"\u003e77\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR81\" class=\"CitationRef\"\u003e81\u003c/span\u003e]. Chronic inflammation is often a contributory factor in the development of squamous cell carcinoma, often linked to tobacco use in lung, head, neck, and esophageal cancers, as well as to ultraviolet radiation in cutaneous squamous cell carcinomas [\u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e75\u003c/span\u003e]. MBC has been observed in the clinical context of chronic inflammation, such as breast abscess associated with a recurrent or long-standing breast cyst or implant capsule [\u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e82\u003c/span\u003e]. Evidence of intermediate EMT states has been identified in single-cell RNA-sequencing and immunohistochemical studies of squamous cell carcinomas, which are hypothesized to play a role in promoting aggressive tumor properties, including the formation of circulating tumor cell clusters that are highly metastatic and exhibit cancer stem cell characteristics [\u003cspan citationid=\"CR83\" class=\"CitationRef\"\u003e83\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eVariants and their derivatives have been widely used to model breast epithelial cell biology and carcinogenesis, facilitating many influential observations. However, it is essential to note that these cultures selectively favor basal epithelial cells, which are thought to be the origin of a small fraction of breast cancers, mainly MBCs [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. Despite this limitation, findings from HMEC culture systems are often extrapolated to breast cancer as a whole, although they may not accurately represent the biology of the more common breast cancer subtypes.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eThis study illuminated dynamic processes operative in the HMEC model system, which is widely used to study epithelial cell biology and carcinogenesis. Following disruption of their native microenvironment, activation of the TGF-β pathway and induction of EMT drove stasis in HMEC cultures, exerting selective pressure that favored the emergence of rare variant basal epithelial cells. The epigenetic stabilization of intermediate EMT states in variants contributed to their enhanced proliferative capacity and avoided complete mesenchymal differentiation, which was associated with irreversible growth arrest. The ability of basal epithelial cells to navigate the spectrum of EMT states and stabilize favorable intermediate states contributed to the emergence of a cancer progenitor state that gives rise to tumors with histopathological features consistent with rare MBCs and may overlap with tumor types initiated from basal epithelial cells, particularly squamous cell carcinomas, in other tissues.\u003c/p\u003e"},{"header":"Abbreviations","content":"\u003cp\u003eHMEC: Human mammary epithelial cell; CDH1: cadherin1 (E-cadherin); CDH2: cadherin2 (N-cadherin); CDH3: cadherin 3 (P-cadherin); CDKN2A: Cyclin-dependent kinase inhibitor 2A; CM: Conditioned medium; EMT: Epithelial-mesenchymal transition; EZH2: Enhancer of zeste homolog 2; H3K27me3: Trimethylation of lysine 27 on histone H3; MBC: Metaplastic breast cancer; MEGM: Mammary epithelial cell growth medium; PRC2: Polycomb repressive complex 2; SASP: Senescence-associated secretory phenotype; SEEP: Stress-elicited extrinsic phenotype; SNAI1: Snail family transcriptional repressor 1(Snail); SNAI2: Snail family transcriptional repressor 2 (Slug); TGF-\u0026beta;: Transforming growth factor-beta; TNBC: Triple-negative breast cancer; TWIST1: Twist family BHLH transcription factor 1 (Twist); ZEB1: Zinc finger E-box binding homeobox 1 (Zeb);\u0026nbsp;\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgments:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors would like to acknowledge Dr. Philippe Gascard for his scientific input and assistance in editing the manuscript, Dr. Deng Pan for his scientific input,\u0026nbsp;Chira Chen-Tanyolac\u0026nbsp;for processing donor material,\u0026nbsp;Xianhong Wang for assistance with experiments,\u0026nbsp;the staff within the Biological Imaging Development CoLab (BIDC) at UCSF Parnassus Heights for their training and support in using the Leica SP8 Confocal microscope, and the Parnassus Flow Cytometry Core (PFFC) for their training and support in using the FACS Aria II Cell Sorter and\u0026nbsp;LSRFortessa flow cytometer.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the California Breast Cancer Research Program (CBCRP) Translational Research Award 22OB-0032 to T.D.T., National Cancer Institute R35 Outstanding Investigator Award 1R35CA197694-01 to T.D.T, and Grand Challenge Cancer Research UK award A27145 to T.D.T. The UCSF Parnassus Flow Cytometry Core (RRID: \u003cem\u003eSCR_018206\u003c/em\u003e) was supported in part by Grant NIH P30 DK063720\u0026nbsp;and by the NIH S10 Instrumentation Grant S10 1S10OD021822-01.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll data generated for this study\u0026nbsp;are available in the paper and its supplementary information files.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor Contributions:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eJ.A.C. and T.D.T Designed research; J.A.C. Performed research; J.A.C. and T.D.T. Analyzed data; J.A.C. and T.D.T. Wrote the paper\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate:\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll specimens were obtained with patient consent under the Institutional Review Board-approved protocol #10-01532. The samples were identified using unlinked codes, following the Federal Health Insurance Portability and Accountability Act (HIPAA) guidelines. Animal experiments were conducted following protocols approved by the University of California, San Francisco Institutional Animal Care and Use Committee.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eVan Keymeulen A, Rocha AS, Ousset M, Beck B, Bouvencourt G, Rock J, Sharma N, Dekoninck S, Blanpain C: \u003cstrong\u003eDistinct stem cells contribute to mammary gland development and maintenance\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2011, \u003cstrong\u003e479\u003c/strong\u003e(7372):189-193.\u003c/li\u003e\n\u003cli\u003eDavis FM, Lloyd-Lewis B, Harris OB, Kozar S, Winton DJ, Muresan L, Watson CJ: \u003cstrong\u003eSingle-cell lineage tracing in the mammary gland reveals stochastic clonal dispersion of stem/progenitor cell progeny\u003c/strong\u003e. \u003cem\u003eNat Commun \u003c/em\u003e2016, \u003cstrong\u003e7\u003c/strong\u003e:13053.\u003c/li\u003e\n\u003cli\u003eStingl J, Eaves CJ, Zandieh I, Emerman JT: \u003cstrong\u003eCharacterization of bipotent mammary epithelial progenitor cells in normal adult human breast tissue\u003c/strong\u003e. \u003cem\u003eBreast Cancer Res Treat \u003c/em\u003e2001, \u003cstrong\u003e67\u003c/strong\u003e(2):93-109.\u003c/li\u003e\n\u003cli\u003eVilladsen R, Fridriksdottir AJ, Ronnov-Jessen L, Gudjonsson T, Rank F, LaBarge MA, Bissell MJ, Petersen OW: \u003cstrong\u003eEvidence for a stem cell hierarchy in the adult human breast\u003c/strong\u003e. \u003cem\u003eJ Cell Biol \u003c/em\u003e2007, \u003cstrong\u003e177\u003c/strong\u003e(1):87-101.\u003c/li\u003e\n\u003cli\u003eGudjonsson T, Villadsen R, Nielsen HL, Ronnov-Jessen L, Bissell MJ, Petersen OW: \u003cstrong\u003eIsolation, immortalization, and characterization of a human breast epithelial cell line with stem cell properties\u003c/strong\u003e. \u003cem\u003eGenes Dev \u003c/em\u003e2002, \u003cstrong\u003e16\u003c/strong\u003e(6):693-706.\u003c/li\u003e\n\u003cli\u003eEirew P, Stingl J, Raouf A, Turashvili G, Aparicio S, Emerman JT, Eaves CJ: \u003cstrong\u003eA method for quantifying normal human mammary epithelial stem cells with in vivo regenerative ability\u003c/strong\u003e. \u003cem\u003eNat Med \u003c/em\u003e2008, \u003cstrong\u003e14\u003c/strong\u003e(12):1384-1389.\u003c/li\u003e\n\u003cli\u003eLim E, Vaillant F, Wu D, Forrest NC, Pal B, Hart AH, Asselin-Labat ML, Gyorki DE, Ward T, Partanen A\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eAberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers\u003c/strong\u003e. \u003cem\u003eNat Med \u003c/em\u003e2009, \u003cstrong\u003e15\u003c/strong\u003e(8):907-913.\u003c/li\u003e\n\u003cli\u003eInce TA, Richardson AL, Bell GW, Saitoh M, Godar S, Karnoub AE, Iglehart JD, Weinberg RA: \u003cstrong\u003eTransformation of different human breast epithelial cell types leads to distinct tumor phenotypes\u003c/strong\u003e. \u003cem\u003eCancer Cell \u003c/em\u003e2007, \u003cstrong\u003e12\u003c/strong\u003e(2):160-170.\u003c/li\u003e\n\u003cli\u003eKeller PJ, Arendt LM, Skibinski A, Logvinenko T, Klebba I, Dong S, Smith AE, Prat A, Perou CM, Gilmore H\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eDefining the cellular precursors to human breast cancer\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2012, \u003cstrong\u003e109\u003c/strong\u003e(8):2772-2777.\u003c/li\u003e\n\u003cli\u003eVan Keymeulen A, Lee MY, Ousset M, Brohee S, Rorive S, Giraddi RR, Wuidart A, Bouvencourt G, Dubois C, Salmon I\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eReactivation of multipotency by oncogenic PIK3CA induces breast tumour heterogeneity\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2015, \u003cstrong\u003e525\u003c/strong\u003e(7567):119-123.\u003c/li\u003e\n\u003cli\u003eMolyneux G, Geyer FC, Magnay FA, McCarthy A, Kendrick H, Natrajan R, Mackay A, Grigoriadis A, Tutt A, Ashworth A\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eBRCA1 basal-like breast cancers originate from luminal epithelial progenitors and not from basal stem cells\u003c/strong\u003e. \u003cem\u003eCell Stem Cell \u003c/em\u003e2010, \u003cstrong\u003e7\u003c/strong\u003e(3):403-417.\u003c/li\u003e\n\u003cli\u003eCancer Genome Atlas N: \u003cstrong\u003eComprehensive molecular portraits of human breast tumours\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2012, \u003cstrong\u003e490\u003c/strong\u003e(7418):61-70.\u003c/li\u003e\n\u003cli\u003eSorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eGene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2001, \u003cstrong\u003e98\u003c/strong\u003e(19):10869-10874.\u003c/li\u003e\n\u003cli\u003ePrat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, He X, Perou CM: \u003cstrong\u003ePhenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer\u003c/strong\u003e. \u003cem\u003eBreast Cancer Res \u003c/em\u003e2010, \u003cstrong\u003e12\u003c/strong\u003e(5):R68.\u003c/li\u003e\n\u003cli\u003ePerou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eMolecular portraits of human breast tumours\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2000, \u003cstrong\u003e406\u003c/strong\u003e(6797):747-752.\u003c/li\u003e\n\u003cli\u003eSantagata S, Thakkar A, Ergonul A, Wang B, Woo T, Hu R, Harrell JC, McNamara G, Schwede M, Culhane AC\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eTaxonomy of breast cancer based on normal cell phenotype predicts outcome\u003c/strong\u003e. \u003cem\u003eJ Clin Invest \u003c/em\u003e2014, \u003cstrong\u003e124\u003c/strong\u003e(2):859-870.\u003c/li\u003e\n\u003cli\u003eReddy TP, Rosato RR, Li X, Moulder S, Piwnica-Worms H, Chang JC: \u003cstrong\u003eA comprehensive overview of metaplastic breast cancer: clinical features and molecular aberrations\u003c/strong\u003e. \u003cem\u003eBreast Cancer Res \u003c/em\u003e2020, \u003cstrong\u003e22\u003c/strong\u003e(1):121.\u003c/li\u003e\n\u003cli\u003eTaube JH, Herschkowitz JI, Komurov K, Zhou AY, Gupta S, Yang J, Hartwell K, Onder TT, Gupta PB, Evans KW\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eCore epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2010, \u003cstrong\u003e107\u003c/strong\u003e(35):15449-15454.\u003c/li\u003e\n\u003cli\u003eThiery JP, Acloque H, Huang RY, Nieto MA: \u003cstrong\u003eEpithelial-mesenchymal transitions in development and disease\u003c/strong\u003e. \u003cem\u003eCell \u003c/em\u003e2009, \u003cstrong\u003e139\u003c/strong\u003e(5):871-890.\u003c/li\u003e\n\u003cli\u003eKalluri R, Weinberg RA: \u003cstrong\u003eThe basics of epithelial-mesenchymal transition\u003c/strong\u003e. \u003cem\u003eJ Clin Invest \u003c/em\u003e2009, \u003cstrong\u003e119\u003c/strong\u003e(6):1420-1428.\u003c/li\u003e\n\u003cli\u003eKroger C, Afeyan A, Mraz J, Eaton EN, Reinhardt F, Khodor YL, Thiru P, Bierie B, Ye X, Burge CB\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eAcquisition of a hybrid E/M state is essential for tumorigenicity of basal breast cancer cells\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2019, \u003cstrong\u003e116\u003c/strong\u003e(15):7353-7362.\u003c/li\u003e\n\u003cli\u003eJolly MK, Boareto M, Huang B, Jia D, Lu M, Ben-Jacob E, Onuchic JN, Levine H: \u003cstrong\u003eImplications of the Hybrid Epithelial/Mesenchymal Phenotype in Metastasis\u003c/strong\u003e. \u003cem\u003eFront Oncol \u003c/em\u003e2015, \u003cstrong\u003e5\u003c/strong\u003e:155.\u003c/li\u003e\n\u003cli\u003eGrosse-Wilde A, Fouquier d\u0026apos;Herouel A, McIntosh E, Ertaylan G, Skupin A, Kuestner RE, del Sol A, Walters KA, Huang S: \u003cstrong\u003eStemness of the hybrid Epithelial/Mesenchymal State in Breast Cancer and Its Association with Poor Survival\u003c/strong\u003e. \u003cem\u003ePLoS One \u003c/em\u003e2015, \u003cstrong\u003e10\u003c/strong\u003e(5):e0126522.\u003c/li\u003e\n\u003cli\u003eMorel AP, Lievre M, Thomas C, Hinkal G, Ansieau S, Puisieux A: \u003cstrong\u003eGeneration of breast cancer stem cells through epithelial-mesenchymal transition\u003c/strong\u003e. \u003cem\u003ePLoS One \u003c/em\u003e2008, \u003cstrong\u003e3\u003c/strong\u003e(8):e2888.\u003c/li\u003e\n\u003cli\u003eMani SA, Guo W, Liao MJ, Eaton EN, Ayyanan A, Zhou AY, Brooks M, Reinhard F, Zhang CC, Shipitsin M\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eThe epithelial-mesenchymal transition generates cells with properties of stem cells\u003c/strong\u003e. \u003cem\u003eCell \u003c/em\u003e2008, \u003cstrong\u003e133\u003c/strong\u003e(4):704-715.\u003c/li\u003e\n\u003cli\u003eChaffer CL, Brueckmann I, Scheel C, Kaestli AJ, Wiggins PA, Rodrigues LO, Brooks M, Reinhardt F, Su Y, Polyak K\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eNormal and neoplastic nonstem cells can spontaneously convert to a stem-like state\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2011, \u003cstrong\u003e108\u003c/strong\u003e(19):7950-7955.\u003c/li\u003e\n\u003cli\u003eQin S, Jiang J, Lu Y, Nice EC, Huang C, Zhang J, He W: \u003cstrong\u003eEmerging role of tumor cell plasticity in modifying therapeutic response\u003c/strong\u003e. \u003cem\u003eSignal Transduct Target Ther \u003c/em\u003e2020, \u003cstrong\u003e5\u003c/strong\u003e(1):228.\u003c/li\u003e\n\u003cli\u003eCorgnac S, Damei I, Gros G, Caidi A, Terry S, Chouaib S, Deloger M, Mami-Chouaib F: \u003cstrong\u003eCancer stem-like cells evade CD8(+)CD103(+) tumor-resident memory T (T(RM)) lymphocytes by initiating an epithelial-to-mesenchymal transition program in a human lung tumor model\u003c/strong\u003e. \u003cem\u003eJ Immunother Cancer \u003c/em\u003e2022, \u003cstrong\u003e10\u003c/strong\u003e(4).\u003c/li\u003e\n\u003cli\u003eHammond SL, Ham RG, Stampfer MR: \u003cstrong\u003eSerum-free growth of human mammary epithelial cells: rapid clonal growth in defined medium and extended serial passage with pituitary extract\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e1984, \u003cstrong\u003e81\u003c/strong\u003e(17):5435-5439.\u003c/li\u003e\n\u003cli\u003eLabarge MA, Garbe JC, Stampfer MR: \u003cstrong\u003eProcessing of human reduction mammoplasty and mastectomy tissues for cell culture\u003c/strong\u003e. \u003cem\u003eJ Vis Exp \u003c/em\u003e2013(71).\u003c/li\u003e\n\u003cli\u003eBissell MJ, Hall HG, Parry G: \u003cstrong\u003eHow does the extracellular matrix direct gene expression?\u003c/strong\u003e \u003cem\u003eJ Theor Biol \u003c/em\u003e1982, \u003cstrong\u003e99\u003c/strong\u003e(1):31-68.\u003c/li\u003e\n\u003cli\u003eDischer DE, Mooney DJ, Zandstra PW: \u003cstrong\u003eGrowth factors, matrices, and forces combine and control stem cells\u003c/strong\u003e. \u003cem\u003eScience \u003c/em\u003e2009, \u003cstrong\u003e324\u003c/strong\u003e(5935):1673-1677.\u003c/li\u003e\n\u003cli\u003eBrenner AJ, Stampfer MR, Aldaz CM: \u003cstrong\u003eIncreased p16 expression with first senescence arrest in human mammary epithelial cells and extended growth capacity with p16 inactivation\u003c/strong\u003e. \u003cem\u003eOncogene \u003c/em\u003e1998, \u003cstrong\u003e17\u003c/strong\u003e(2):199-205.\u003c/li\u003e\n\u003cli\u003eRomanov SR, Kozakiewicz BK, Holst CR, Stampfer MR, Haupt LM, Tlsty TD: \u003cstrong\u003eNormal human mammary epithelial cells spontaneously escape senescence and acquire genomic changes\u003c/strong\u003e. \u003cem\u003eNature \u003c/em\u003e2001, \u003cstrong\u003e409\u003c/strong\u003e(6820):633-637.\u003c/li\u003e\n\u003cli\u003eLocke WJ, Zotenko E, Stirzaker C, Robinson MD, Hinshelwood RA, Stone A, Reddel RR, Huschtscha LI, Clark SJ: \u003cstrong\u003eCoordinated epigenetic remodelling of transcriptional networks occurs during early breast carcinogenesis\u003c/strong\u003e. \u003cem\u003eClin Epigenetics \u003c/em\u003e2015, \u003cstrong\u003e7\u003c/strong\u003e(1):52.\u003c/li\u003e\n\u003cli\u003eNovak P, Jensen TJ, Garbe JC, Stampfer MR, Futscher BW: \u003cstrong\u003eStepwise DNA methylation changes are linked to escape from defined proliferation barriers and mammary epithelial cell immortalization\u003c/strong\u003e. \u003cem\u003eCancer Res \u003c/em\u003e2009, \u003cstrong\u003e69\u003c/strong\u003e(12):5251-5258.\u003c/li\u003e\n\u003cli\u003eHuschtscha LI, Noble JR, Neumann AA, Moy EL, Barry P, Melki JR, Clark SJ, Reddel RR: \u003cstrong\u003eLoss of p16INK4 expression by methylation is associated with lifespan extension of human mammary epithelial cells\u003c/strong\u003e. \u003cem\u003eCancer Res \u003c/em\u003e1998, \u003cstrong\u003e58\u003c/strong\u003e(16):3508-3512.\u003c/li\u003e\n\u003cli\u003eHinshelwood RA, Huschtscha LI, Melki J, Stirzaker C, Abdipranoto A, Vissel B, Ravasi T, Wells CA, Hume DA, Reddel RR\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eConcordant epigenetic silencing of transforming growth factor-beta signaling pathway genes occurs early in breast carcinogenesis\u003c/strong\u003e. \u003cem\u003eCancer Res \u003c/em\u003e2007, \u003cstrong\u003e67\u003c/strong\u003e(24):11517-11527.\u003c/li\u003e\n\u003cli\u003eDumont N, Crawford YG, Sigaroudinia M, Nagrani SS, Wilson MB, Buehring GC, Turashvili G, Aparicio S, Gauthier ML, Fordyce CA\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eHuman mammary cancer progression model recapitulates methylation events associated with breast premalignancy\u003c/strong\u003e. \u003cem\u003eBreast Cancer Res \u003c/em\u003e2009, \u003cstrong\u003e11\u003c/strong\u003e(6):R87.\u003c/li\u003e\n\u003cli\u003eBerman H, Zhang J, Crawford YG, Gauthier ML, Fordyce CA, McDermott KM, Sigaroudinia M, Kozakiewicz K, Tlsty TD: \u003cstrong\u003eGenetic and epigenetic changes in mammary epithelial cells identify a subpopulation of cells involved in early carcinogenesis\u003c/strong\u003e. \u003cem\u003eCold Spring Harb Symp Quant Biol \u003c/em\u003e2005, \u003cstrong\u003e70\u003c/strong\u003e:317-327.\u003c/li\u003e\n\u003cli\u003eReynolds PA, Sigaroudinia M, Zardo G, Wilson MB, Benton GM, Miller CJ, Hong C, Fridlyand J, Costello JF, Tlsty TD: \u003cstrong\u003eTumor suppressor p16INK4A regulates polycomb-mediated DNA hypermethylation in human mammary epithelial cells\u003c/strong\u003e. \u003cem\u003eJ Biol Chem \u003c/em\u003e2006, \u003cstrong\u003e281\u003c/strong\u003e(34):24790-24802.\u003c/li\u003e\n\u003cli\u003eElenbaas B, Spirio L, Koerner F, Fleming MD, Zimonjic DB, Donaher JL, Popescu NC, Hahn WC, Weinberg RA: \u003cstrong\u003eHuman breast cancer cells generated by oncogenic transformation of primary mammary epithelial cells\u003c/strong\u003e. \u003cem\u003eGenes Dev \u003c/em\u003e2001, \u003cstrong\u003e15\u003c/strong\u003e(1):50-65.\u003c/li\u003e\n\u003cli\u003eRizki A, Weaver VM, Lee SY, Rozenberg GI, Chin K, Myers CA, Bascom JL, Mott JD, Semeiks JR, Grate LR\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eA human breast cell model of preinvasive to invasive transition\u003c/strong\u003e. \u003cem\u003eCancer Res \u003c/em\u003e2008, \u003cstrong\u003e68\u003c/strong\u003e(5):1378-1387.\u003c/li\u003e\n\u003cli\u003eDontu G, Abdallah WM, Foley JM, Jackson KW, Clarke MF, Kawamura MJ, Wicha MS: \u003cstrong\u003eIn vitro propagation and transcriptional profiling of human mammary stem/progenitor cells\u003c/strong\u003e. \u003cem\u003eGenes Dev \u003c/em\u003e2003, \u003cstrong\u003e17\u003c/strong\u003e(10):1253-1270.\u003c/li\u003e\n\u003cli\u003eBankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, McQuaid S, Gray RT, Murray LJ, Coleman HG\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eQuPath: Open source software for digital pathology image analysis\u003c/strong\u003e. \u003cem\u003eSci Rep \u003c/em\u003e2017, \u003cstrong\u003e7\u003c/strong\u003e(1):16878.\u003c/li\u003e\n\u003cli\u003eSchmidt U, Weigert M, Broaddus C, Myers G: \u003cstrong\u003eCell Detection with Star-Convex Polygons\u003c/strong\u003e. In\u003cem\u003e: 2018; Cham\u003c/em\u003e: Springer International Publishing; 2018: 265-273.\u003c/li\u003e\n\u003cli\u003eHolst CR, Nuovo GJ, Esteller M, Chew K, Baylin SB, Herman JG, Tlsty TD: \u003cstrong\u003eMethylation of p16(INK4a) promoters occurs in vivo in histologically normal human mammary epithelia\u003c/strong\u003e. \u003cem\u003eCancer Res \u003c/em\u003e2003, \u003cstrong\u003e63\u003c/strong\u003e(7):1596-1601.\u003c/li\u003e\n\u003cli\u003eCrawford YG, Gauthier ML, Joubel A, Mantei K, Kozakiewicz K, Afshari CA, Tlsty TD: \u003cstrong\u003eHistologically normal human mammary epithelia with silenced p16(INK4a) overexpress COX-2, promoting a premalignant program\u003c/strong\u003e. \u003cem\u003eCancer Cell \u003c/em\u003e2004, \u003cstrong\u003e5\u003c/strong\u003e(3):263-273.\u003c/li\u003e\n\u003cli\u003eAcosta JC, Banito A, Wuestefeld T, Georgilis A, Janich P, Morton JP, Athineos D, Kang TW, Lasitschka F, Andrulis M\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eA complex secretory program orchestrated by the inflammasome controls paracrine senescence\u003c/strong\u003e. \u003cem\u003eNat Cell Biol \u003c/em\u003e2013, \u003cstrong\u003e15\u003c/strong\u003e(8):978-990.\u003c/li\u003e\n\u003cli\u003eFordyce CA, Patten KT, Fessenden TB, DeFilippis R, Hwang ES, Zhao J, Tlsty TD: \u003cstrong\u003eCell-extrinsic consequences of epithelial stress: activation of protumorigenic tissue phenotypes\u003c/strong\u003e. \u003cem\u003eBreast Cancer Res \u003c/em\u003e2012, \u003cstrong\u003e14\u003c/strong\u003e(6):R155.\u003c/li\u003e\n\u003cli\u003eCoppe JP, Desprez PY, Krtolica A, Campisi J: \u003cstrong\u003eThe senescence-associated secretory phenotype: the dark side of tumor suppression\u003c/strong\u003e. \u003cem\u003eAnnu Rev Pathol \u003c/em\u003e2010, \u003cstrong\u003e5\u003c/strong\u003e:99-118.\u003c/li\u003e\n\u003cli\u003eMiettinen PJ, Ebner R, Lopez AR, Derynck R: \u003cstrong\u003eTGF-beta induced transdifferentiation of mammary epithelial cells to mesenchymal cells: involvement of type I receptors\u003c/strong\u003e. \u003cem\u003eJ Cell Biol \u003c/em\u003e1994, \u003cstrong\u003e127\u003c/strong\u003e(6 Pt 2):2021-2036.\u003c/li\u003e\n\u003cli\u003eKim ES, Kim MS, Moon A: \u003cstrong\u003eTGF-beta-induced upregulation of MMP-2 and MMP-9 depends on p38 MAPK, but not ERK signaling in MCF10A human breast epithelial cells\u003c/strong\u003e. \u003cem\u003eInt J Oncol \u003c/em\u003e2004, \u003cstrong\u003e25\u003c/strong\u003e(5):1375-1382.\u003c/li\u003e\n\u003cli\u003eScheel C, Eaton EN, Li SH, Chaffer CL, Reinhardt F, Kah KJ, Bell G, Guo W, Rubin J, Richardson AL\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eParacrine and autocrine signals induce and maintain mesenchymal and stem cell states in the breast\u003c/strong\u003e. \u003cem\u003eCell \u003c/em\u003e2011, \u003cstrong\u003e145\u003c/strong\u003e(6):926-940.\u003c/li\u003e\n\u003cli\u003ePhillips S, Prat A, Sedic M, Proia T, Wronski A, Mazumdar S, Skibinski A, Shirley SH, Perou CM, Gill G\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eCell-state transitions regulated by SLUG are critical for tissue regeneration and tumor initiation\u003c/strong\u003e. \u003cem\u003eStem Cell Reports \u003c/em\u003e2014, \u003cstrong\u003e2\u003c/strong\u003e(5):633-647.\u003c/li\u003e\n\u003cli\u003eZhang Y, Donaher JL, Das S, Li X, Reinhardt F, Krall JA, Lambert AW, Thiru P, Keys HR, Khan M\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eGenome-wide CRISPR screen identifies PRC2 and KMT2D-COMPASS as regulators of distinct EMT trajectories that contribute differentially to metastasis\u003c/strong\u003e. \u003cem\u003eNat Cell Biol \u003c/em\u003e2022, \u003cstrong\u003e24\u003c/strong\u003e(4):554-564.\u003c/li\u003e\n\u003cli\u003eTavassoli FA: \u003cstrong\u003eMyoepithelial lesions of the breast. Myoepitheliosis, adenomyoepithelioma, and myoepithelial carcinoma\u003c/strong\u003e. \u003cem\u003eAm J Surg Pathol \u003c/em\u003e1991, \u003cstrong\u003e15\u003c/strong\u003e(6):554-568.\u003c/li\u003e\n\u003cli\u003eTse GM, Tan PH, Chaiwun B, Putti TC, Lui PC, Tsang AK, Wong FC, Lo AW: \u003cstrong\u003ep63 is useful in the diagnosis of mammary metaplastic carcinomas\u003c/strong\u003e. \u003cem\u003ePathology \u003c/em\u003e2006, \u003cstrong\u003e38\u003c/strong\u003e(1):16-20.\u003c/li\u003e\n\u003cli\u003eKoker MM, Kleer CG: \u003cstrong\u003ep63 expression in breast cancer: a highly sensitive and specific marker of metaplastic carcinoma\u003c/strong\u003e. \u003cem\u003eAm J Surg Pathol \u003c/em\u003e2004, \u003cstrong\u003e28\u003c/strong\u003e(11):1506-1512.\u003c/li\u003e\n\u003cli\u003eRosenbluth JM, Schackmann RCJ, Gray GK, Selfors LM, Li CM, Boedicker M, Kuiken HJ, Richardson A, Brock J, Garber J\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eOrganoid cultures from normal and cancer-prone human breast tissues preserve complex epithelial lineages\u003c/strong\u003e. \u003cem\u003eNat Commun \u003c/em\u003e2020, \u003cstrong\u003e11\u003c/strong\u003e(1):1711.\u003c/li\u003e\n\u003cli\u003eLiu X, Ory V, Chapman S, Yuan H, Albanese C, Kallakury B, Timofeeva OA, Nealon C, Dakic A, Simic V\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eROCK inhibitor and feeder cells induce the conditional reprogramming of epithelial cells\u003c/strong\u003e. \u003cem\u003eAm J Pathol \u003c/em\u003e2012, \u003cstrong\u003e180\u003c/strong\u003e(2):599-607.\u003c/li\u003e\n\u003cli\u003eFeinberg AP, Koldobskiy MA, Gondor A: \u003cstrong\u003eEpigenetic modulators, modifiers and mediators in cancer aetiology and progression\u003c/strong\u003e. \u003cem\u003eNat Rev Genet \u003c/em\u003e2016, \u003cstrong\u003e17\u003c/strong\u003e(5):284-299.\u003c/li\u003e\n\u003cli\u003eFeinberg AP, Ohlsson R, Henikoff S: \u003cstrong\u003eThe epigenetic progenitor origin of human cancer\u003c/strong\u003e. \u003cem\u003eNat Rev Genet \u003c/em\u003e2006, \u003cstrong\u003e7\u003c/strong\u003e(1):21-33.\u003c/li\u003e\n\u003cli\u003eLamouille S, Derynck R: \u003cstrong\u003eCell size and invasion in TGF-beta-induced epithelial to mesenchymal transition is regulated by activation of the mTOR pathway\u003c/strong\u003e. \u003cem\u003eJ Cell Biol \u003c/em\u003e2007, \u003cstrong\u003e178\u003c/strong\u003e(3):437-451.\u003c/li\u003e\n\u003cli\u003eNeurohr GE, Terry RL, Lengefeld J, Bonney M, Brittingham GP, Moretto F, Miettinen TP, Vaites LP, Soares LM, Paulo JA\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eExcessive Cell Growth Causes Cytoplasm Dilution And Contributes to Senescence\u003c/strong\u003e. \u003cem\u003eCell \u003c/em\u003e2019, \u003cstrong\u003e176\u003c/strong\u003e(5):1083-1097 e1018.\u003c/li\u003e\n\u003cli\u003eManohar S, Estrada ME, Uliana F, Vuina K, Alvarez PM, de Bruin RAM, Neurohr GE: \u003cstrong\u003eGenome homeostasis defects drive enlarged cells into senescence\u003c/strong\u003e. \u003cem\u003eMol Cell \u003c/em\u003e2023, \u003cstrong\u003e83\u003c/strong\u003e(22):4032-4046 e4036.\u003c/li\u003e\n\u003cli\u003eLanz MC, Zatulovskiy E, Swaffer MP, Zhang L, Ilerten I, Zhang S, You DS, Marinov G, McAlpine P, Elias JE\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eIncreasing cell size remodels the proteome and promotes senescence\u003c/strong\u003e. \u003cem\u003eMol Cell \u003c/em\u003e2022, \u003cstrong\u003e82\u003c/strong\u003e(17):3255-3269 e3258.\u003c/li\u003e\n\u003cli\u003eDemidenko ZN, Blagosklonny MV: \u003cstrong\u003eGrowth stimulation leads to cellular senescence when the cell cycle is blocked\u003c/strong\u003e. \u003cem\u003eCell Cycle \u003c/em\u003e2008, \u003cstrong\u003e7\u003c/strong\u003e(21):3355-3361.\u003c/li\u003e\n\u003cli\u003eCrozier L, Foy R, Adib R, Kar A, Holt JA, Pareri AU, Valverde JM, Rivera R, Weston WA, Wilson R\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eCDK4/6 inhibitor-mediated cell overgrowth triggers osmotic and replication stress to promote senescence\u003c/strong\u003e. \u003cem\u003eMol Cell \u003c/em\u003e2023, \u003cstrong\u003e83\u003c/strong\u003e(22):4062-4077 e4065.\u003c/li\u003e\n\u003cli\u003eComaills V, Kabeche L, Morris R, Buisson R, Yu M, Madden MW, LiCausi JA, Boukhali M, Tajima K, Pan S\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eGenomic Instability Is Induced by Persistent Proliferation of Cells Undergoing Epithelial-to-Mesenchymal Transition\u003c/strong\u003e. \u003cem\u003eCell Rep \u003c/em\u003e2016, \u003cstrong\u003e17\u003c/strong\u003e(10):2632-2647.\u003c/li\u003e\n\u003cli\u003eKnouse KA, Lopez KE, Bachofner M, Amon A: \u003cstrong\u003eChromosome Segregation Fidelity in Epithelia Requires Tissue Architecture\u003c/strong\u003e. \u003cem\u003eCell \u003c/em\u003e2018, \u003cstrong\u003e175\u003c/strong\u003e(1):200-211 e213.\u003c/li\u003e\n\u003cli\u003eLopez-Diaz FJ, Gascard P, Balakrishnan SK, Zhao J, Del Rincon SV, Spruck C, Tlsty TD, Emerson BM: \u003cstrong\u003eCoordinate transcriptional and translational repression of p53 by TGF-beta1 impairs the stress response\u003c/strong\u003e. \u003cem\u003eMol Cell \u003c/em\u003e2013, \u003cstrong\u003e50\u003c/strong\u003e(4):552-564.\u003c/li\u003e\n\u003cli\u003eGarbe JC, Vrba L, Sputova K, Fuchs L, Novak P, Brothman AR, Jackson M, Chin K, LaBarge MA, Watts G\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eImmortalization of normal human mammary epithelial cells in two steps by direct targeting of senescence barriers does not require gross genomic alterations\u003c/strong\u003e. \u003cem\u003eCell Cycle \u003c/em\u003e2014, \u003cstrong\u003e13\u003c/strong\u003e(21):3423-3435.\u003c/li\u003e\n\u003cli\u003eCipriano R, Kan CE, Graham J, Danielpour D, Stampfer M, Jackson MW: \u003cstrong\u003eTGF-beta signaling engages an ATM-CHK2-p53-independent RAS-induced senescence and prevents malignant transformation in human mammary epithelial cells\u003c/strong\u003e. \u003cem\u003eProc Natl Acad Sci U S A \u003c/em\u003e2011, \u003cstrong\u003e108\u003c/strong\u003e(21):8668-8673.\u003c/li\u003e\n\u003cli\u003eJohnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR: \u003cstrong\u003eHead and neck squamous cell carcinoma\u003c/strong\u003e. \u003cem\u003eNat Rev Dis Primers \u003c/em\u003e2020, \u003cstrong\u003e6\u003c/strong\u003e(1):92.\u003c/li\u003e\n\u003cli\u003eFletcher CD, Beham A, Schmid C: \u003cstrong\u003eSpindle cell haemangioendothelioma: a clinicopathological and immunohistochemical study indicative of a non-neoplastic lesion\u003c/strong\u003e. \u003cem\u003eHistopathology \u003c/em\u003e1991, \u003cstrong\u003e18\u003c/strong\u003e(4):291-301.\u003c/li\u003e\n\u003cli\u003eYarbrough WG, Aprelikova O, Pei H, Olshan AF, Liu ET: \u003cstrong\u003eFamilial tumor syndrome associated with a germline nonfunctional p16INK4a allele\u003c/strong\u003e. \u003cem\u003eJ Natl Cancer Inst \u003c/em\u003e1996, \u003cstrong\u003e88\u003c/strong\u003e(20):1489-1491.\u003c/li\u003e\n\u003cli\u003eReed AL, Califano J, Cairns P, Westra WH, Jones RM, Koch W, Ahrendt S, Eby Y, Sewell D, Nawroz H\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eHigh frequency of p16 (CDKN2/MTS-1/INK4A) inactivation in head and neck squamous cell carcinoma\u003c/strong\u003e. \u003cem\u003eCancer Res \u003c/em\u003e1996, \u003cstrong\u003e56\u003c/strong\u003e(16):3630-3633.\u003c/li\u003e\n\u003cli\u003eGonzalez-Zulueta M, Shibata A, Ohneseit PF, Spruck CH, 3rd, Busch C, Shamaa M, El-Baz M, Nichols PW, Gonzalgo ML, Elbaz M\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eHigh frequency of chromosome 9p allelic loss and CDKN2 tumor suppressor gene alterations in squamous cell carcinoma of the bladder\u003c/strong\u003e. \u003cem\u003eJ Natl Cancer Inst \u003c/em\u003e1995, \u003cstrong\u003e87\u003c/strong\u003e(18):1383-1393.\u003c/li\u003e\n\u003cli\u003eFoulkes WD, Flanders TY, Pollock PM, Hayward NK: \u003cstrong\u003eThe CDKN2A (p16) gene and human cancer\u003c/strong\u003e. \u003cem\u003eMol Med \u003c/em\u003e1997, \u003cstrong\u003e3\u003c/strong\u003e(1):5-20.\u003c/li\u003e\n\u003cli\u003eBartels S, van Luttikhuizen JL, Christgen M, Magel L, Luft A, Hanzelmann S, Lehmann U, Schlegelberger B, Leo F, Steinemann D\u003cem\u003e et al\u003c/em\u003e: \u003cstrong\u003eCDKN2A loss and PIK3CA mutation in myoepithelial-like metaplastic breast cancer\u003c/strong\u003e. \u003cem\u003eJ Pathol \u003c/em\u003e2018, \u003cstrong\u003e245\u003c/strong\u003e(3):373-383.\u003c/li\u003e\n\u003cli\u003eBehranwala KA, Nasiri N, Abdullah N, Trott PA, Gui GP: \u003cstrong\u003eSquamous cell carcinoma of the breast: clinico-pathologic implications and outcome\u003c/strong\u003e. \u003cem\u003eEur J Surg Oncol \u003c/em\u003e2003, \u003cstrong\u003e29\u003c/strong\u003e(4):386-389.\u003c/li\u003e\n\u003cli\u003ePal A, Barrett TF, Paolini R, Parikh A, Puram SV: \u003cstrong\u003ePartial EMT in head and neck cancer biology: a spectrum instead of a switch\u003c/strong\u003e. \u003cem\u003eOncogene \u003c/em\u003e2021, \u003cstrong\u003e40\u003c/strong\u003e(32):5049-5065.\u003c/li\u003e\n\u003cli\u003eHu Y, Smyth GK: \u003cstrong\u003eELDA: extreme limiting dilution analysis for comparing depleted and enriched populations in stem cell and other assays\u003c/strong\u003e. \u003cem\u003eJ Immunol Methods \u003c/em\u003e2009, \u003cstrong\u003e347\u003c/strong\u003e(1-2):70-78.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"breast-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"brcr","sideBox":"Learn more about [Breast Cancer Research](http://breast-cancer-research.biomedcentral.com)","snPcode":"13058","submissionUrl":"https://submission.nature.com/new-submission/13058/3","title":"Breast Cancer Research","twitterHandle":"@BCRJournal","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Human mammary epithelial cells, Basal, Myoepithelial, TGF-β Pathway, Epithelial-mesenchymal transition, Epigenetic Regulation, Metaplastic Breast Cancer","lastPublishedDoi":"10.21203/rs.3.rs-4980285/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4980285/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e Human mammary epithelial cell (HMEC) cultures encounter a stress-associated barrier termed stasis, during which most cells adopt a senescence-like phenotype. From these cultures, rare variants emerge from the basal epithelial population, re-initiating growth. Variants exhibit pre-malignant properties, including an aberrant epigenetic program that enables continued proliferation and acquisition of genetic changes. Following oncogenic transformation, variants produce tumors that recapitulate the histopathological characteristics of metaplastic breast cancer (MBC), a rare subtype characterized by squamous and mesenchymal differentiation.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMethods:\u003c/strong\u003e Using the conventional serum-free HMEC culture system, we probed the capacity for phenotypic plasticity inherent to basal epithelial cell populations from human breast tissue as they navigated stasis and emerged as variant populations.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e We observed robust activation of a TGF-β-dependent epithelial-mesenchymal transition (EMT) program in basal epithelial cells during stasis, followed by subsequent attenuation of this program in emerging variants. Inhibiting the TGF-β pathway or depleting the EMT regulators Snail or Slug allowed basal epithelial cells to collectively bypass stasis, demonstrating that cellular dysfunction and arrest resulting from TGF-β and EMT activation are central to this \u003cem\u003ein vitro\u003c/em\u003e barrier. The spontaneous emergence of variants from stasis cultures was associated with a restricted EMT trajectory, which diverted cells away from a complete mesenchymal state characterized by irreversible growth arrest, and instead limited variants to epithelial and intermediate EMT states associated with greater proliferative capacity and stemness. Epigenetic mechanisms, which contributed to the dysregulated growth control characteristic of the variant phenotype, also contributed to the constrained EMT program in variants. By overcoming the cellular dysfunction and growth arrest resulting from TGF-β and EMT activation, variants exhibited increased oncogenic transformation efficiency compared to pre-stasis basal epithelial cells. Inhibiting the TGF-β pathway prior to stasis significantly reduced EMT in the basal epithelial population, alleviated selective pressure driving variant emergence, and enhanced oncogenic transformation efficiency, resulting in tumors with markedly diminished metaplastic differentiation.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConclusions:\u003c/strong\u003e This study reveals how adaptive EMT reprogramming governs basal epithelial cell fate decisions and contributes to the development of MBC progenitors by restricting access to terminal mesenchymal states that induce growth arrest and, instead, favoring intermediate states with enhanced tumorigenic potential.\u003c/p\u003e","manuscriptTitle":"An adaptive Epithelial-Mesenchymal Transition Program Enables Basal Epithelial Cells to Bypass Stress-Induced Stasis and Contributes to Metaplastic Breast Cancer Progenitor State","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-27 15:25:59","doi":"10.21203/rs.3.rs-4980285/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-10-22T22:10:33+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-10-22T20:00:56+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-10-02T17:42:46+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"5270817372473980230338103506374546026","date":"2024-09-24T07:17:52+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"149745722092046641409612186451411574356","date":"2024-09-23T16:55:23+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"26137803603677095963702289703791667740","date":"2024-09-13T07:58:55+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-09-05T14:26:24+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-08-29T13:00:53+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-08-29T07:23:57+00:00","index":"","fulltext":""},{"type":"submitted","content":"Breast Cancer Research","date":"2024-08-26T21:52:46+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"breast-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"brcr","sideBox":"Learn more about [Breast Cancer Research](http://breast-cancer-research.biomedcentral.com)","snPcode":"13058","submissionUrl":"https://submission.nature.com/new-submission/13058/3","title":"Breast Cancer Research","twitterHandle":"@BCRJournal","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"acc1d16d-0b80-48d0-89cc-2608fc0e90b4","owner":[],"postedDate":"September 27th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-12-23T16:09:48+00:00","versionOfRecord":{"articleIdentity":"rs-4980285","link":"https://doi.org/10.1186/s13058-024-01920-8","journal":{"identity":"breast-cancer-research","isVorOnly":false,"title":"Breast Cancer Research"},"publishedOn":"2024-12-18 15:58:35","publishedOnDateReadable":"December 18th, 2024"},"versionCreatedAt":"2024-09-27 15:25:59","video":"","vorDoi":"10.1186/s13058-024-01920-8","vorDoiUrl":"https://doi.org/10.1186/s13058-024-01920-8","workflowStages":[]},"version":"v1","identity":"rs-4980285","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4980285","identity":"rs-4980285","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}