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

doi:10.21203/rs.3.rs-4980285/v1

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

https://doi.org/10.21203/rs.3.rs-4980285/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

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

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.

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.

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. 1A). Nearly complete proliferative arrest was observed at stasis (Fig. 1B), marked by the upregulation of p16 and development of a senescence-like morphology (Fig. 1C). 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. 1D). 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. 1E). Additionally, we detected an increase in matrix metalloproteinase 2 (MMP-2) and fibronectin (FN) transcript levels during stasis (Fig. 1F), 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. S1A).

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. 1G). 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. 1G). 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. 1G). 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. 1H). In addition to TGF-β (Fig. 1D), 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. S1B 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. 2A). 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. 2B), 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. 2C).

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. 2D). 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. 2E). 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. 2F). 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. 2G). 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. 3A). 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. 3B). TGF-β inhibition was associated with a highly significant decrease in cell size (Fig. 3C), reduction in nuclear area, and elimination of micronuclei (Fig. 3D). A83-01 treatment also significantly reduced p21 and p16 levels (Fig. 3E). We observed that treatment with A83-01 maintained the protein expression of p63, an important transcriptional regulator of basal cell differentiation (Fig. 3F). 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. 4A). 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. 4A). 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. 4B). Basal epithelial cells with Snail or Slug KO demonstrated substantially lower transcription of p21 and p16 tumor suppressors than the controls (Fig. 4C). 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. S2A). 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. 4D), 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. S2B). 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. S2C). 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. 4E).

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. 5A). 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^hiCD24^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. 5B). 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. 5C). 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. 5D). 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. 5E).

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. 5F). 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. 5G). 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. 5H). 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. 5I).

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. 6A). EZH2 knockout in the P8 variant basal epithelial cells also impaired cell proliferation (Fig. 6B).

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. 6C). EZH2 knockout promoted EMT at the transcriptional level, as evidenced by decreased P-cadherin, but increased N-cadherin, Snail, Twist, and Zeb expression (Fig. 6D).

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. 6E).

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. 7A).

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. 7B). 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. 7C). 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. 7D). 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. 7E). 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. 7F).

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. 7G).

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.

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.

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.

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);

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.

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.

No competing interests reported.

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.

Download PDF

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

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

Status:

Version 1

Abstract

Figures

Background

Methods

Results

Basal Epithelial Cells Exhibit Enhanced EMT Activation Compared to Luminal Epithelial Cells

Inhibition of TGF-β Signaling and EMT Transcription Factors Prevents Stasis in Basal Epithelial Cells

EMT Generates a Spectrum of Basal Epithelial Cell States with Varying Proliferative and Self-Renewal Capacities

PRC2-Mediated Epigenetic Regulation Governs Access to EMT States in Basal Epithelial Variants

TGF-β Inhibition Prior to Cellular Transformation Reduced Metaplastic Phenotypes in Basal Epithelial-Derived Tumors

Discussion

Conclusions

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1