Introduction

Cancer has traditionally been understood to be a consequence of deregulated cell proliferation driven by genetic alterations that enable clonal expansion and tumor formation [1]. Interestingly, genomic analysis reveals that many tumor suppressor genes and oncogenes exhibit tissue specificity, meaning they are altered in some cancers but not others [2-4]. Current research suggests that a combination of intrinsic biological factors [4], such as differential epigenetic changes [3,5] and aneuploidy patterns [6-8], as well as extrinsic factors like the tumour microenvironment, including the presence of tissue-specific immune cells [9-10] influences tissue-specific cancer outcomes. However, a relatively underexplored yet crucial aspect of tissue-specific oncogenesis is the mechanical interaction between normal and cancerous cells. To this end, recent studies have highlighted a tumor suppressive behavior in epithelial tissues known as epithelial defense against cancer (EDAC) [11], where transformed cells are expelled from the tissue by surrounding normal cells [12-14]. However, EDAC can fail [14], allowing mutant cells to evade extrusion when the microenvironment becomes tumour-permissive, such as in conditions of increased stiffness [15] or inflammation due to a high-fat diet [16]. The success or failure of EDAC is closely linked to changes in cellular contractility [12-14] and cell-cell adhesions [17], suggesting that the underlying biophysical properties of the tissue play a decisive role in determining the outcome of cell competition. Since these properties can vary across different epithelial tissues, tissue-specific differences in cancer susceptibility may be explained by the varying ability of epithelia to manage newly transformed cells. However, despite these advances, the mechanical basis of tissue-specific oncogenesis remains missing. Specifically, how differences in cellular mechanics at the interface between oncogenic and wild-type populations influence tumor initiation and progression is poorly understood. Addressing this gap could not only enhance our understanding of cancer biology but also open new avenues for tissue-specific therapeutic interventions by identifying what makes certain tissues more permissive to tumor formation.

To this end, we investigated and compared the bio-physical factors influencing early tumorigenesis in human mammary and bronchial epithelial monolayers by sporadically transfecting a subset of healthy cells into oncogenic cells expressing a constitutively active HRasV12 protein. We then tracked the fate of these transfected cells and their interactions with the wild-type counterparts in real time via live imaging. Strikingly, HRasV12-mutant cells exhibited contrasting behaviors in the two tissues. In breast epithelium, isolated transformed cells were frequently extruded. However, when present in groups, extrusion of oncogenic cells was not seen, but instead, they remained spatially confined in a jammed state. In contrast, in lung epithelium, single as well as groups of HRasV12 expressing cells persisted and exhibited unrestricted expansion over time, regardless of cluster size. To explain these findings, we employed a bi-disperse vertex model of epithelia [18-21] and demonstrated that differences in interfacial tension between wild-type and transformed populations play a crucial role in determining whether mutant cells are eliminated, restrained, or are allowed to persist and expand. Taken together, these findings provide novel insights into the differential outcomes of competitive cellular interactions in the precancer stages and uncover a mechanical basis for tissue-specific oncogenesis.

Results

1. Singlets and groups of HRasV12 oncogenic mutants show tissue-specific fates in epithelia

To explore the differences in the premalignant stages across epithelial tissues, we sporadically transfected monolayers of the two non-transformed human epithelial cell lines, Mcf10a (mammary) and Beas2b (bronchial), with a plasmid construct that constitutively expresses HRasV12 along with the GFP reporter. We then monitored the cellular dynamics over 90 hours. Consistent with previous studies reporting extrusion of single HRasV12 mutants in MDCK monolayers [11-15, 22], we observed that isolated HRasV12 transfected cells (singlets) surrounded by wild-type neighbours were extruded from the mammary epithelium [Figure. 1b, Supplementary video 1]. This process correlated with a reduction in shape indices of oncogenic mutants as the density of wild-type cells around them increased [Figure. 1d (upper panel)]. Interestingly, this well-known mechanism of epithelial defense against cancer was impaired in Beas2b cells, where HRasV12 mutants evaded extrusion [Figure. 1c, e], and wild-type cell compaction did not restrict their growth [Figure. 1c, Supplementary video 1]. Instead, the mutants in the bronchial epithelium continued to proliferate and exhibited increased elongation with long protrusions [Figure. 1c (lower panel), Supplementary video 1]. Given that cancers frequently arise from fields of oncogenic mutants rather than isolated mutated cells [23-25], we also investigated whether these differential dynamics of mutant singlets would also extend to groups of oncogenic mutants, potentially enabling more aggressive HRasV12 growth in lung tissue compared to breast tissue. By adjusting cell seeding conditions, we generated larger groups of HRasV12 transfected cells in both Mcf10a and Beas2b monolayers and then tracked the dynamics of these groups over time. In the mammary epithelium, wild-type cells were unable to extrude groups of HRasV12-transfected cells but instead confined them into compact, circular clusters with smooth, rounded interfaces [Figure 1g, (upper panel), i (upper panel)] Supplementary video 2]. As a result, oncogenic clusters gradually demixed from the surrounding wild-type Mcf10a population and remained spatially constrained [Figure 1j]. In contrast, HRasV12 clusters in the bronchial epithelium gradually spread outward [Figure. 1h, j, Supplementary video 2], formed protrusive lamellipodia, grew on top of one another [Figure. 2b] and developed an irregular interface with surrounding wild-type cells [Figure 1h, i (lower panel)]. Together, these striking differences in the dynamics of HRasV12 oncogenic cells demonstrate that activation of the same oncogene can lead to distinct outcomes across tissues in terms of growth, spreading, and extrusion. This prompted us to further investigate the unique bio-physical mechanisms driving these variations.

Tissue-specific outcomes of HRasV12 oncogenic mutants in epithelial monolayers

(a) Schematic representation of the experimental setup and outcome (b) Representative images showing extrusion of single oncogenic cells (singlets) in mammary epithelium (upper panel), but protrusive behavior of the same oncogenic singlets in bronchial epithelium (bottom panel) (c) Quantification of shape indices of HRasV12 singlets showing a reduction in shape indices, with increasing wild-type cell density in mammary epithelium (upper panel), but the opposite trend in bronchial epithelium (bottom panel), where HRasV12 cells continues to spread. Shaded regions indicate the results simple linear regression analysis with 95% confidence intervals (d) Extrusion rates of HRasV12 singlets are significantly higher in mammary epithelium, while (e) Percentage of HRasV12 cells forming protrusions is markedly higher in bronchial epithelium. Representative data are plotted from one of three independent experiments with the median shown as a bold dashed line and the first and third quartiles are shown as thin dashed lines. Statistical significance was calculated using Unpaired t-test with Welch’s correction. (f) Behaviour of HRasV12 oncogenic clusters in the two tissues-In mammary epithelium (upper panel), clusters become spatially confined with a smooth, circular interface with the wild-type population, while in bronchial epithelium (bottom panel), HRasV12 clusters expand, forming long protrusions over time (g) Immunofluorescence images of HRasV12 clusters in mammary and bronchial epithelia, highlighting differences in spreading patterns in the two epithelia (h) Quantification of cluster segregation: mammary epithelium exhibits a higher segregation index and reduced eccentricity compared to bronchial epithelium, indicating greater spatial confinement of oncogenic clusters. The segregation index (SI), defined as the average ratio of homotypic and all cell neighbors, quantifies the degree of demixing and is averaged over all cells inside an oncogenic cluster. Data are mean±sem and plotted from different clusters from one representative experiment. (Scale bars= 50 μm)

Distinct mechanical signatures of HRasV12 clusters in mammary and bronchial epithelia

(a) Representative images showing regions of interest (ROIs) selected for analysis in mammary (upper panel) and bronchial epithelia (lower panel) (b) Immunofluorescence images of HRasV12 clusters stained for F-actin in mammary (upper panel) and bronchial epithelia (lower panel) (c) Oncogenic cluster–wild-type interface showing distinct actin belt in mammary marked by white arrows (upper panel) and absence of it in bronchial epithelia (lower panel) (d) ROI-based F-actin stained regions of wild-type and oncogenic cells in mammary (upper panel) and bronchial epithelia (lower panel) revealing distinct shape differences in different regions-in mammary epithelia, wild type cells at the interface show elongated shapes (ROI2), and oncogenic cells inside the cluster show jammed shapes (ROI4), in comparison to wild type cells at a random location (ROI1), while in bronchial epithelia, oncogenic cells (ROI 3 and 4) show more elongated shapes compared to wild-type beas2b cells (e), (f) Quantifications of F-actin intensities and shape indices, in mammary epithelia (e), and bronchial epithelia (f), across the four ROIs confirming jammed oncogenic clusters, with unjamming at the interface, in mammary epithelia and unjammed clusters in bronchial epithelia (g) Schematic illustrating the Bayesian Force Inference pipeline used to estimate relative intracellular pressure and cell–cell edge tension (h) Heatmaps depicting relative cell pressures and cell–cell edge tensions in mammary epithelium, revealing localized mechanical heterogeneity (i) Heatmaps of relative cell pressures and edge tensions in bronchial epithelium, showing a distinct mechanical landscape compared to mammary epithelium (j) Quantification of relative intracellular pressures and cell–cell edge tensions across the four ROIs in mammary (upper panel) and bronchial epithelia (lower panel) highlighting tissue-specific differences. All data represented as mean±sem in box and whisker plots and are plotted from one out of three independent experiments with whiskers extending from the minimum to the maximum values, showing the full range of the dataset including outliers. Line represents the median. Statistical significance was calculated using Unpaired t-test with Welch’s correction. Scale bars= 20 μm

2. Biophysical analysis reveals demixing via jamming of oncogenic clusters in Mcf10A but unjamming and protrusive growth in Beas2b

To understand the mechanisms driving the differential fates of oncogenic fields in the two epithelial tissues, we mapped the bio-physical signatures of wild-type-mutant interactions by subdividing the images into four regions of interests for both monolayers-ROI1: wild type cells distant from oncogenic cells, ROI2: interfacial wild type cells in direct contact with oncogenic cells, ROI3: oncogenic cells at the cluster boundaries in direct contact with wild type cells, and ROI4: oncogenic cells in contact with other oncogenic cells [Figure. 2a]. We then compared the F-actin levels and shapes of both mutant and wild-type populations across these ROIs. Additionally, we employed Bayesian force inference, as previously established [26], to infer the differential tissue stresses associated with cluster formation from the observed cell geometries. Interestingly, the differences across these parameters were apparent not only between wild-type and oncogenic counterparts but also based on the spatial localization of the cells, with cells at the oncogenic-wild type interface exhibiting distinct behaviors, regardless of their origin.

In mammary epithelium, oncogenic clusters were surrounded by a prominent actin belt [Figure. 2b (upper panel), c (upper panel)], coupled with drastically low F-actin levels inside the clusters [Figure. 2e (left panel)], suggesting that the de-mixing of oncogenic cells in breast tissue is likely driven by changes in cellular mechanics. Further, oncogenic clusters in mammary epithelium (ROI4) also exhibited low cell shape indices [Figure. 2d (upper panel), f (left panel)], high internal pressure [Figure. 2h (upper panel), j (upper panel)], and low line tension [Figure. 2h (upper panel), j (lower panel)]-characteristics indicative of jamming [27]. Notably, both oncogenic cells and wildtype Mcf10a cells at the cluster interface (ROI 2 and ROI3) showed F-actin enrichment, while the interfacial wild shape cells show elongated shapes [Figure. 2e (left panel), f (left panel)], indicating that interfacial cellular mechanics play a key role in the de-mixing process as the wild-type cells encircle the oncogenic mutants. This interfacial region enriched with F-actin also contained phosphomyosin [Supplementary Figure. 1a], indicating that a contractile actomyosin belt drives the compaction of oncogenic clusters in mammary epithelia. In contrast, bronchial epithelium exhibited no actin belt formation at the wild type-oncogenic cell interface, and no traits of jamming were observed in the oncogenic clusters [Figure. 2b (lower panel), c (lower panel)]. Instead, oncogenic cells continued to proliferate, displaying elongation measured via higher shape indices and high F-actin expression [Figure. 2d (lower panel), e (right panel), f (right panel)]-indicative of unjamming behaviour [27,28]. The relative internal cell pressure and cell-cell line tension inferred with Bayesian force inference showed slight or no modulation in the oncogenic clusters, given the general elongated shapes of wild-type cells in Beas2b cells [Figure. 2i, j]. Together, these findings highlight distinct interfacial mechanics between oncogenic cells and their wild-type neighbors in both Mcf10A and Beas2b monolayers, resulting in differing outcomes in mutant colony growth.

3. PIV analysis reveals kinetic arrest of oncogenic clusters in mcf10a coupled with tangential motion of interfacial cells

To further probe the cellular dynamics that lead to the contrasting fates of the oncogenic mutant colonies in the two tissues, we analyzed cellular velocities during the de-mixing of oncogenic cells, using Particle Image Velocimetry (PIV) on consecutive image pairs from mammary and bronchial monolayers. Interestingly, velocity maps in samples imaged for mutant cluster – wild type dynamics in Mcf10a [Supplementary video 3] revealed an unusual movement pattern at confluency: wild type Mcf10a cells exhibited tangential motion around oncogenic clusters [Figure. 3b], which coupled with a drastic reduction in the velocity of oncogenic cells [Figure. 3a], likely caused their compaction and jamming into clusters. Notably, this tangential movement was unique to wild-type Mcf10a cells surrounding oncogenic clusters and was absent in other regions across the monolayer [Figure. 3b] as well as under control conditions-i.e., in the absence of oncogenic mutants [Supplementary video 4], in both of which a typical velocity drop was observed in Mcf10a cells due to density-induced jamming [27-30] [Figure. 3d (upper panel), Supplementary 2]. These results aligned with the observed differences in cell shape indices, where oncogenic clusters underwent jamming, while wild-type cells at the interface exhibited unjamming behavior [29,30]. In contrast, in the bronchial epithelium, mutant clusters remained unjammed, even as surrounding wild-type cells became more crowded [Figure. 3c, d (lower panel), Supplementary video 5]. Additionally, no tangential cellular movement was detected around the oncogenic clusters in the bronchial epithelium. These velocity maps showed negligible deviations from wild type (control) beas2b velocity maps [Supplementary video 6]. Together, these findings suggest that mutant clusters carrying the same oncogene can exhibit drastically different behaviors depending on the epithelial tissue type. Specifically, differences in the dynamics of the interfacial cells drive the demixing and jamming of HRasV12 oncogenic cells in mammary epithelium while promoting protrusive, unjammed growth in bronchial epithelium.

Particle Image Velocimetry (PIV) reveals distinct cellular movement patterns in HRasV12 clusters in the two monolayers

(a), (b) Representative snapshots of HRasV12 clusters, and the corresponding PIV maps, in mammary (a) and bronchial epithelium (b), as the density of wild type cells increase showing wild type cells along the cluster-wild-type interface exhibit a higher velocity attributed to tangential motion in mammary epithelium (a), but undergo a more kinetically arrested state in bronchial epithelium (b) (c) Oncogenic cluster in mcf10a monolayer with concentric circles drawn from the centre of the cluster, radially outward till its boundary (upper panel) with the red (highlighted) circle representing the one along which tangential velocity was the highest and plot showing tangential motion (expressed as a function of the circles), revealing highest tangential motion (data point highlighted with red) along the interface of oncogenic cluster with the wild-type cells in mammary epithelia (lower panel). No significant tangential motion in random wild-type region (c) Representative snapshots of HRasV12 clusters in bronchial epithelium at the same two time points and the corresponding PIV velocity maps (d) Quantification of cellular velocities in local regions as cell densities reach confluency showing oncogenic cells clustered in mammary epithelium (top panel) jam and show a significant reduction in velocities. In bronchial epithelium (bottom panel) while wild type cells show reduced movement as cell density increased-transformed cells continued to show higher motion, signifying jamming of wild type cells have no effect on mutants. All data are represented as mean±sem plotted from different clusters from one representative experiment. Scale bars= 50 μm

4. Bi-disperse vertex model reveals differential interfacial tension between mutant and wild-type population governs tissue-specific response towards oncogenic mutants

To explain the differential dynamics and degrees of demixing of oncogenic clusters observed in the two tissues, we drew inspiration from the differential interfacial tension hypothesis (DITH) [31]. This hypothesis suggests that cell sorting in tissues could occur due to the variations in interfacial tension, as seen in embryonic cell rearrangement, where cells regulate contractility, as well as adhesions to demix from their neighbors [31-33]. To test whether such a mechanism could be applied to the de-mixing of oncogenic clusters from wild-type cells-leading to tissue-specific outcomes, we employed the cell-based bi-disperse Vertex Model (VM). Our model allowed us to simulate an epithelial monolayer and capture the dynamic interactions between wild-type and mutant cell populations. We analyzed the mechanics underlying the segregation dynamics of mutant clusters by examining the interfacial line tension between the oncogenic and wild-type cells [Figure. 4a]. In our model, Λij = 0, for homotypic interfaces, i.e. wild-type-wild type or mutant-mutant interactions, and Λij = Λ ≠ 0, for heterotypic interfaces between oncogenic and wild type cells. When Λ > 0, our model ensures that the interface between oncogenic and healthy cells contracts to minimize energy, leading to a reduction in total interfacial length, and when Λ < 0, our model promotes a rough interface by favouring increased interface between healthy and oncogenic cells [Supplementary video 7]. The final interface shape is determined by the balance between interfacial tension and the intrinsic perimeter and area elasticity of individual cells [Figure. 4b].

Bi-disperse vertex model reveals interfacial tension-driven segregation of oncogenic clusters

(a) Schematic of the bi-disperse vertex model, where interfacial line tension (red) is applied at the boundary between mutant (green) and wild-type (grey) (b) The value of Λ determines whether the mutant cluster remains compact or spreads into the wild-type population, where the red line gives the outline of the morphology of the mutant cluster interface (c) Difference in shape index between interfacial and bulk cells as a function of Λ (interfacial tension), shown for different stress thresholds. Experimental average (red) overlaid for comparison (d) Eccentricity of mutant clusters (given by the red lines in (b)) over time for different values of Λ (with Π0 = 0.3), indicating changes in cluster morphology (e) Demixing parameter at t = 300 as a function of Λ forΠ0 = 0.3, quantifying the extent of cluster segregation (f) Time evolution of the demixing parameter for different values of Λ (with Π0 = 0.3), demonstrating interfacial tension-driven segregation (g) Hydrostatic pressure maps of mutant and wild-type cells at t = 300, with mutant cells labeled by green directors and wild-type cells by black directors (h) Velocity field maps of epithelial monolayers for different values of Λ, revealing flow patterns around oncogenic clusters (i) Tangential velocity loop integral for circles of varying radius centered at the mutant cluster (Region 1) and centered away from the cluster (Region 2), showing tissue-specific velocity variations (j) Extrusion probability of a single mutant cell as a function of Λ, demonstrating how interfacial tension influences mutant cell fate.

To explain the observed differences in interfacial cellular dynamics, we introduced an active contribution to interfacial tension, which was coupled to the shape of cells at the interface, called shape-tension coupling[34]. This active interfacial tension is extensile, meaning it elongates cells along the interface depending on their orientation relative to the interface. It originates from the anisotropic distributions of stress-fibres in when cells are under external stress. Hence, we assume that it becomes significant only when the hydrostatic pressure on the interfacial cells exceeds a threshold value Π0 [Figure. 4b, g].

Remarkably, our model robustly captured the experimentally observed behaviors. As a first check, we compared the experimentally observed differences in the shape indices between interfacial and non-interfacial wild-type and oncogenic cells in both mammary and lung epithelium with our model’s predictions. For Γ = 0.2, Π0 = {0, 0.1, 0.3}, Λ = {0 − 1.6}, our model accurately captured the experimental observations [Figure. 4c]. For stress threshold values lower than Π0 = 0.3, we find that the differences in shape indices are non-zero even for Λ = 0, which is an unphysical situation. This is unphysical because we expect there to be no difference in the shape indices because in principle, there is only one cell type for Λ = 0. Thus, we have taken Π0 = 0.3 for further analyses. To make the comparison between the model and the experiment more robust, we compared other metrics as well.

To quantify the extent of mixing between wild-type and mutant cells, we used the demixing parameter (DP), defined as: , where Ns represents the number of same-type neighbours, and Nt is the total number of neighbours of a given cell[35]. Our simulations showed that when Λ < 0 the demixing parameter rapidly decreases over time, indicating outward protrusion of the mutant cluster. In contrast, for Λ > 0 the demixing parameter remained close to 1, reflecting a compact, jammed cluster [Figure. 4f]. Additionally, we observed that when Λ > 0, cells within the cluster experience higher internal pressure compared to when Λ < 0 [Figure 4g], consistent with experimental findings [Figure. 2h,i]. Furthermore, flow field analysis from our simulations [Figure 4h,i] revealed that for Λ > 0, cell movement aligns tangentially along the cluster boundaryconsistent with the Mcf10A behaviourwhereas for Λ < 0 cell motion was uniform throughout the monolayer, resembling the Beas2b phenotype [Figure. 3a, b]. These comparions further establish that different values of Λ correspond to different cell types. The same oncogene thus might affect different cell lines differently. We also estimated the probability of single mutant cell extrusion as a function of Λ [Figure. 4j], consistent with the findings from our experiments which show that isolated HRasV12 transfected cells are frequently extruded in mammary tissue [Figure. 1b]. Together, our findings demonstrate how interfacial tension-driven mechanics could play a crucial role in determining the competitive dynamics between wild-type and mutant populations, thus providing compelling evidence that cellular and tissue mechanics are key determinants of cell competition.

5. Disrupting actomyosin belt using blebbistatin prevents demixing and compaction of mutant clusters in mammary epithelia

Since our model showed that differences in interfacial tension regulate mutant cluster dynamics, with a higher interfacial tension driving cluster compaction [Supplementary video 7], we reasoned that this tension is mediated by the supracellular actomyosin belt [36] that drives the compaction and segregation of oncogenic clusters in mammary epithelium. Therefore, we investigated whether disrupting the cytoskeletal organization and impairing actomyosin belt formation would affect the demixing and compaction of oncogenic clusters. To test this, we treated Mcf10a monolayers with 5 µM blebbistatin, a myosin II inhibitor, to disrupt belt formation.

As expected, blebbistatin treatment prevented mutant cluster compaction [Figure. 5a], leading to a loss of their rounded morphology and restoring actin levels and shape indices of mutant cells to wild-type levels [Figure. 5b]. This shift resulted in mutant clusters resembling those observed in bronchial epithelia, suggesting that cytoskeletal organization plays a key role in tissue-specific differences in mutant cluster dynamics. This was also consistent with the fact that interfacial tension between cells arises from the interplay between the cortical actomyosin network and cell-cell adhesion [36].

Actomyosin belt disruption prevents mutant cluster compaction in mammary epithelia

(a) Representative actin staining image upon blebbistatin treatment, done after oncogenic induction in mammary epithelial monolayer, showing an absence of an actin belt (b) Quantification of F-actin intensities (left) and shape indices (right), showing difference between oncogenic clusters and wild-type cells in these parameters, is getting small – although not completely vanished because of low concentration of blebbistatin used for the experiment (c) Representative actin and E-cadherin staining images of the wild type monolayers showing shorter stress fibers, higher E-cadherin intensities in Mcf10a wild-type monolayer, and the opposite in Beas2b (d) Plots comparing F-actin stress fiber length (left panel) and shape indices (right panel) in wild-type Mcf10a and Beas2b monolayers. Stress fiber lengths were calculated manually in ImageJ (e) Plots comparing F-actin and E-cadherin levels in the two wild-type tissues. Representative data are plotted from one of three independent experiments with the median shown as a bold line. Statistical significance was calculated using Unpaired t-test with Welch’s correction. Scale bars= 50 μm.

Additionally, we looked into the inherent differences in cytoskeletal mechanics of the host tissues and observed distinct differences in actin organization between the two tissues. Wild-type bronchial epithelial cells exhibited higher F-actin intensities [Figure. 5c (upper panel), e] and had longer stress fibers [Figure. 5c (middle panel), d (left panel)] compared to wild-type mammary epithelial cells. These cytoskeletal differences were also coupled with significant differences in cell shapes between confluent Mcf10a and Beas2b monolayers with wild-type beas2b cells exhibited significantly higher shape indices, indicating elongated, loosely packed morphologies and a lower degree of jamming [Figure 5d (right panel). In contrast, mammary epithelial cells had lower shape indices, reflecting a more compact and tightly packed arrangement [Figure. 5d (right panel)]. These differences were also in agreement with the PIV maps of the two wild-type tissues [Supplementary figure 2, Supplementary videos 4 and 6] where the mammary monolayers underwent a sharper drop in velocities and therefore kinetic arrest, compared to the bronchial epithelia. Bronchial epithelial cells also exhibited weaker cell-cell adhesions, as indicated by reduced E-cadherin staining intensity at cell junctions compared to mammary epithelial cells [Figure. 5c (bottom panel), e]. These weaker adhesions likely contributed to increased cellular motility and a diminished ability to constrain mutant clusters, in contrast to the more stable, compacted clusters observed in mammary epithelium. Together, these findings demonstrate that variations in mechanical properties-including cell shape, actin organization, and adhesion strength-govern the differential responses of mammary and bronchial epithelia to oncogenic mutants and may in turn, influence tumor growth and progression.

Discussion

Our findings reveal a striking tissue-specific response to HRasV12 oncogene activation in epithelial monolayers, highlighting key differences in the growth, segregation, and mechanical behavior of mutant clusters in mammary and bronchial epithelia. These insights shed light on the differential susceptibility of epithelial tissues to oncogenic transformation and uncover a mechanical basis for tissue-specific oncogenesis. By demonstrating how interfacial tension at the mutant-wild-type boundary governs whether oncogenic cells are restrained or overgrow, our work shifts the focus from purely molecular pathways to the biophysical constraints that regulate tumor initiation.

While previous studies have emphasized the contextual nature of tumorigenesis [2-4] and highlighted the critical role of interactions between healthy and oncogenic cells during epithelial cancer initiation [14, 22, 38], our work integrates and extends these perspectives by demonstrating how interfacial tension at the mutant-wild-type boundary determines whether transformed cells remain constrained or undergo unchecked expansion.

HRasV12 mutant clusters exhibited contrasting behaviors in mammary and bronchial epithelia. In the mammary epithelium, mutant clusters underwent demixing, and got restrained into compact, circular aggregates with smooth interfaces whereas in bronchial epithelia, they dispersed, exhibited protrusive growth, and formed jagged interfaces with surrounding wild-type cells. Bio-physical characterizations and PIV analysis revealed tissue specificity in cytoskeletal changes, and, subsequently cell shape changes, movement patterns, and internal pressures in both wild type and oncogenic populations, coupled with distinct mechanical behaviours of interfacial cells at the boundary of the two populations.

To probe the mechanistic basis of these differences, we applied the Differential Interfacial Tension Hypothesis (DITH) [31,32]and simulated mutant-wild-type interactions using a bi-disperse Vertex model. Our simulations demonstrated that interfacial tension is a key determinant of cluster morphology and compaction. When interfacial tension was positive, mutant clusters compacted into circular aggregates, recapitulating Mcf10A behavior. In contrast, negative interfacial tension (Λ<0) resulted in irregular, protrusive clusters that failed to compact, mirroring the Beas2b phenotype. Agreeing with our model, we showed that disrupting the interfacial actomyosin belt in the mammary epithelia disrupted the restraint on oncogenic clusters, preventing their demixing and compaction and shifting mutant behavior to a more unjammed state. These findings underscore the role of interfacial cellular mechanics in shaping oncogenic cell fate, with mammary epithelium imposing physical constraints that limit mutant expansion, while bronchial epithelium lacks such mechanical barriers, permitting unchecked growth.

By integrating mechanical constraints into models of cancer initiation, our study provides a framework for understanding how different epithelial tissues resist or permit tumorigenesis. However, while interfacial tension emerges as a key regulator of oncogenic segregation, the molecular pathways driving these mechanical differences remain to be elucidated. Future studies should investigate how targeting specific signaling networks or cytoskeletal components influence interfacial tension and whether modulating these factors could provide therapeutic benefits.

Additionally, extending this framework to other oncogenes and tissue types could further elucidate the mechanical principles governing tumor initiation and progression.

Materials and methods

Cell culture

Non-transformed human mammary epithelial cell line Mcf10a were maintained in complete medium composed of phenol-free Dulbecco’s modified Eagle’s medium (DMEM–F12, Gibco) supplemented with 5% charcoal-stripped horse serum (Gibco), 10 U ml−1 penicillin and 10 μg ml−1 streptomycin (Pen-Strep, Invitrogen), epidermal growth factor (20 ng/ml; PeproTech), hydrocortisone (0.5 mg/ml; Sigma-Aldrich), cholera toxin (100 ng/ml; Sigma-Aldrich), and insulin (10 μg/ml; Sigma-Aldrich) at 37°C in a humidified incubator with 5% CO2.

Non-transformed human bronchial epithelium epithelial cell line Beas2b were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco) supplemented with GlutaMax (Gibco) and 5% fetal bovine serum (FBS, Gibco) along with 10 U ml−1 penicillin and 10 μg ml−1 streptomycin (Pen-Strep, Invitrogen). Cells were maintained at 37 °C and 5% CO2 unless mentioned otherwise.

Cell seeding

105 Cells were seeded on fibronectin-coated (10 μg ml-1) glass bottom dishes (35 mm Cellvis) and transfected once they reached 80% confluency (for singlets) and 60% confluency to get transfected groups.

Transfection

Cells were transiently transfected with the HrasV12-GFP expressing plasmid DNA (Addgene #18780) using Lipofectamine 3000 (Invitrogen), following the manufacturer’s protocol. Post-transfection, confluent monolayers were either fixed and immuno-stained or imaged directly.

Blebbistatin addition

Blebbistatin was diluted to 5 μM final solution and added 48h post-transfection. Samples were washed after 24 hours and fixed right away.

Immunofluorescence

Cells were fixed with 4% paraformaldehyde (ThermoFisher) diluted in 1X Phosphate buffered saline (PBS, Sigma-Aldrich) for 15 minutes at room temperature. After washing away the fixative with 1X PBS, cells were permeabilized with 0.3% TritonX-100 (Sigma-Aldrich) in 1X PBS. Nonspecific antibody binding was blocked by incubating the samples in blocking solution (2% Bovine Serum Albumin (BSA Sigma-Aldrich) and 2% FBS in PBS) at room temperature for 45 minutes. Further, cells were incubated with primary antibody prepared in blocking/staining solution overnight at 4°C. Post primary antibodies incubation, cells were washed thrice with 1X PBS for 5 minutes per wash and incubated with secondary antibodies conjugated to Alexa fluor 555, Alexa fluor 594 conjugated phalloidin and 4′,6-diamidino-2-phenylindole (DAPI) to mark the cell nucleus, prepared in blocking/staining solution, for 2 hours at room temperature. Finally, the samples were washed four times with 1X PBS for 5 minutes per wash before proceeding to microscopy.

Microscopy

Fluorescence images were acquired using 20X objective and 63x oil-immersion objective of a Zeiss Axio Observer 7 Inverted Microscope with a scientific-grade sCMOS camera by Iberoptics. For live imaging, samples were set up on a stage-top incubator and maintained at 37°C and 5% CO2 throughout imaging.

Image Analysis

Image analysis for this study was performed using Fiji [39] except bayesian force inference and particle image velocimetry (PIV) analysis, which was performed in MATLAB (MathWorks). Cellpose [40] was used to segment cells to get cell count, shape indices and actin intensitites for each cell while Tissue Analyzer [41] was used to segment cells for force inference.

Quantification of cluster eccentricity and segregation index

The term ‘eccentricity’ expresses a size reduction in the contact with neighboring cells. It was quantified as the ratio of the length of the minor axis to the major axis of an ellipse, with a value closer to 0 indicating a more circular shape. Eccentricity quantifications were performed on fixed and live tissues using outlines that were manually drawn and measured using Fiji. The segregation index SI, defined as the average ratio of homotypic and all cell neighbors as done previously [42], quantifies the demixing degree. To quantify the SI, cells of each type were manually counted, and finally, to generate the plots, the SI was averaged over all GFP-labelled cells in one frame.

Bayesian Force Inference

Bayesian force inference was done as described previously [26] to determine the relative cellular pressures and cell-cell edge or line tension.

Particle Image Velocimetry (PIV)

Particle Image Velocimetry (PIV) was done using the PIV Lab package on MATLAB [43]. FFT window deformation (direct Fourier transform correlation with multiple passes and deforming windows) algorithm was used. Velocity maps were obtained by a 3-pass interrogation window of 128 × 128 pixels that halved after each pass. Tangential velocities were obtained by drawing concentric circles radially outward from the center to a point twice the radius of the mutant cluster and integrated over the entire circle using a custom python code. The output is the tangential velocity as a function of the circle.

Statistical analysis

Statistical analyses were carried out in GraphPad Prism 10. Statistical significance was calculated by Unpaired t-test with Welch’s correction. p-values greater than 0.05 were considered to be statistically not significant. Each plot has data pooled from multiple frames of one out of three independent experiments. Violin plots have median shown as a bold dashed line and the first and third quartiles shown as thin dashed lines. Box and whisker plots have whiskers extending from the minimum to the maximum values, showing the full range of the dataset including outliers. Scatter-bar plots were displayed as mean ± s.e.m. No statistical methods were used to set the sample size. Quantification was done using data from at least three independent biological replicates. For analysis involving live-imaging experiments, data were collected from three independent experiments. All the experiments with representative images were repeated at least three times.

Vertex Model

A bi-disperse model is used to demonstrate the role of interfacial tension between two types of cells, specifically the HRasV12 mutant and wild-type cells.

The first two sums are included in the standard vertex model, whereas the last term is introduced in our model to account for interfacial interactions between cells. Specifically Λij = 0 for interfaces between same cell types, and Λij = Λ between different cell types. Using a line tension term like Λlij in the vertex model energy, where the lij is the length of an interface edge, we can obtain both segregation and spreading dynamics just from energy minimization. This term has been described as heterotypic line tension (HLT)[42]. Due to energy minimization, heterotypic interface length is minimized for Λ > 0 and maximized for Λ < 0, leading to segregation of the two cell types in the former case and complete mixing of them in the latter case.

The mutant cluster in the experiments with Mcf10a cells showed an actin belt forming along the cluster interface, with the cells along the interface being more elongated than normal. This actin ring is important for the compaction of mutant clusters in the presence of wild-type cells in the epithelial monolayer. An isotropic interfacial tension, such as the Λlij term in our model cannot capture such anisotropic behaviour. Recently, it has been proposed that such anisotropic behaviour can arise from the anisotropic distribution of the cytoskeletal filaments inside the cell[41]. The resulting active interfacial tension, called the shape-tension coupling, can generate elongated shapes as observed in the experiments. The elongated cells in the experiments indeed have anisotropic distributions of the stress fibres. Hence, we included the effect of shape-tension in our model. The total interfacial tension is Λij + γij(t; Π0), where Λij is a time-independent, passive, isotropic tension and γij is a time-dependent active tension that depends on the shape of the cells sharing the interface and requires a threshold stress Π0 to be nonzero. The details of the shape-tension coupling is described in the interface. We find that a nonzero activation stress Π0 is necessary to reproduce experimentally observed results.

Supplementary figures

Myosin belt and per cell myosin intensities

(a) Images of HRasV12 cluster in mammary epithelia stained for phosphomyosin (left panel) and phosphomyosin belt around the cluster marked with white arrowheads (right panel) (b) Plot showing lower phosphomyosin levels inside oncogenic clusters in mcf10a. Representative data are plotted from one of three independent experiments with the median shown as a bold line. Statistical significance was calculated using Unpaired t-test with Welch’s correction. Scale bars= 50 μm.

PIV maps of wild-type mammary and bronchial epithelial monolayers with increasing cell density

(a) Representative time-lapse images of wild-type (control) mammary epithelium at two different time points, with corresponding PIV velocity maps showing a reduction in overall velocity as cell density increases (b) Representative time-lapse images of wild-type (control) bronchial epithelium at two different time points, with corresponding PIV velocity maps demonstrating a similar reduction in velocity with increasing cell density (c) Quantification of velocity reduction in wild-type mammary and bronchial epithelia as a function of cell density, highlighting tissue-specific differences in collective cell dynamics. All data are represented as mean±sem plotted from different clusters from one representative experiment. Scale bars= 50 μm.

Illustration of shape-tension coupling at the mutant (green)–wild-type (gray) interface (red), where angles θ1 and θ2 represent the orientation of boundary cells relative to their shared edge

Acknowledgements

We thank Sriram R. Ramaswamy, and Tamal Das for critical discussions and suggestions. M.V. is a partner group leader of the Max Planck Society (MPG), Germany, which has supported part of this work. Open access funding provided by Max Planck Society. This work is also supported by the Infosys foundation, Anusandhan National Research Foundation-previously called the Science and Engineering Research Board (project number: SERB SRG/2022/000534), and Indo German Science and Technology Centre (IGSTC WISER scheme). SS acknowledges funding from IISc, Axis Bank Center for Mathematics and Computing, and a startup grant from SERB-DST (SRG/2022/000163). We also acknowledge intramural funds at IISc Bangalore for providing support towards equipment and facilities and for salaries/fellowships of the authors.

Additional information

Author contributions

M.V conceived the project. M.V. and A.D. designed experiments. A.D performed all experiments except the staining on Beas2b cells, which were performed by A.A. Theoretical model was contributed by P.D and S.S. Analysis and interpretation of experimental data was done by A.D, T.C, T.T, A.A, S.M, A.R and M.V. A.D and M.V developed and wrote the manuscript with help from P.D and S.S. All authors read, discussed and commented on the manuscript.

Funding

Max Planck Society (Max Planck Partner Group)

Infosys (India) (Young investigator award)

Anusandhan National Research Foundation, India (SRG/2022/000534)

Indo German Science and Technology Centre (WISER Fellowship 2023)

Axis Bank Center for Mathematics and Computing

Anusandhan National Research Foundation India (SRG/2022/000163)

Additional files

Supplementary video 1: An HRasV12 singlet is extruded from Mcf10a monolayer but persists in Beas2b monolayer, multiply and form long protrusion. Scale bar= 20 μm.

Supplementary video 2: A group of HRasV12 transfected cells gradually gets confined into a compact cluster in Mcf10a monolayer and forms a smooth interface with the wild-type population, while a group of HRasV12 transfected cells gradually grows in size in Beas2b monolayer, forming a jagged interface with the wild-type population along with protrusions. Scale bar= 50 μm.

Supplementary video 3: PIV video during formation of a representative HRasV12 cluster in Mcf10a monolayer. Scale bar= 50 μm.

Supplementary video 4: PIV video showing the dynamics of a representative wild-type Mcf10a (control) monolayer without any oncogenic mutants. Scale bar= 50 μm.

Supplementary video 5: PIV video showing the dynamics of a representative HRasV12 cluster in Beas2b monolayer. Scale bar= 50 μm.

Supplementary video 6: PIV video showing the dynamics of a representative wild-type Beas2b (control) monolayer without any oncogenic mutants. Scale bar= 50 μm.

Supplementary video 7: Video showing output of the bi-disperse vertex model with cluster dynamics varying with interfacial tension (Λ).

Appendix

Vertex Model Simulation

To understand the cause and the dynamics of interaction between mutant and wild-type population in different epithelial monolayers during cancer initiation, we used a cell-based model of epithelial tissues called the vertex model. A monolayer in the vertex model has total energy given by:

In the above energy function, the constant parameters Kα is the area modulus which gives the area elasticity of each cell and Γα is the perimeter modulus which give the elasticity due to the actomyosin cortex. The gradient of this energy function gives the forces on each vertex of the cell, , where is the position of the ith vertex. The dynamics of each vertex position is given by the overdamped equation of motion:

where is any active force which might drive the system out of equilibrium. is a force which propels individual cells by F , where is the polarity of a given cell given by , where this angle θα performs rotational diffusion:

where ηα(t) is a Gaussian white noise with zero mean and correlation ⟦ηα(t)ηα (t′)⟧ = δ(tt′) δαα.

To capture the dynamics of cell segregation between two different types of cells in a monolayer of normal and oncogenic cells, we can use a bi-disperse vertex model with an interfacial line tension between the two types of cells - called mutant and wild type can be used.

The director of each cell corresponds to the polarization direction of the cell, which can be obtained as the larger eigenvector of its shape tensor, , where the sum is over N(α) vertices of the cell α, and is the position of the cell’s geometric center. The energy function for the vertex model used is given below:

Here, M denotes mutant cells and WT denotes wild-type cells. The index α denotes cells, and the indices i, j denote vertices and < ij > denotes an edge between vertices i and j. The parameters are the area moduli, are the perimeter moduli of the mutant and wild type respectively, and the Λij denotes the interfacial line tension. This interfacial line tension is only along junctions connecting two vertices and depends on the cell type on either side of the junction. Since the system we are considering is bi-disperse, the interfacial line tension can have three values, depending on the identities of cells on either side of a given edge < ij >:

Shape-tension coupling

In addition to this passive interfacial tension, we can also have an active contribution coming from the shape tension coupling. The shape-tension coupling is mathematically the same as the passive line tension force. The total force is given by:

γij(t) depends on the relative orientation of the cells director and the edge along which the tension acts,.(fig) Specifically, it has a relaxational dynamics of the form:

Here couples the tension along an edge to the shape of the cells on either side of the edge [Supplementary figure 3].

For ζ > 0, this tension is extensile and leads to the elongation of cells along the boundary. This happens because for the edges which have and aligned with the edge vector (θ or θ′ = 0, π, 2π), the shape-tension coupling contributes a negative line tension, which favours this alignment. When the edge vectors are not aligned with the edge vector, the contribution is positive, and that edge starts to shrink. Eventually, the interface consists of cells whose directors are always aligned with the interface. Energetically, it counteracts interfaces with Λ > 0, leading to the formation of elongated cells and tangential flows along the interface. Importantly, we find that cell elongation correlates with the formation of actin belt, leading us to postulate that the shape-tension coupling is activated above a threshold pressure, Π0, where the cell pressure is given by , When Λ < 0, such as in the Beas2b cells, we do not observe strong stress-fibre expression. Hence, we believe that it is not activated in those cells.

Simulation Setup

The vertex model energy can be non-dimensionalised by setting a length scale set by . In our simulation we have taken

and

Also,

and

The non-dimensionalization has been done, which defines the normalized parameters . (In further discussions, the normalized parameters will be represented by the symbols Γ and Λ.)

The monolayer modelled by the vertex model can be described by a single number, the shape index, defined by .Here too, in principle we could have different shape indices for the mutant and wild type cells, but we take them to be the same. Unless otherwise mentioned, we take ΛM/M = ΛWT/WT = 0 and ΛWT/M = Λ.

To simulate the elongation of cells along the interface, ζ > 0 is used for an active extensile contribution to the interfacial tension.

A box containing 20 × 20 cells has been used with mutant cells placed within a region (non-circular) from the center of the box. The vertices are arranged in a hexagonal lattice with random displacements added to the position of each vertex. Periodic boundary conditions were used to simulate a confluent monolayer.

Parameters

To prevent overlaps, a node switch operation is also added to the simulation, which changes the topology of the network [21].