Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 3087 (2023)
1012
1
1
Metrics details
Cell mechanical properties have been proposed as label free markers for diagnostic purposes in diseases such as cancer. Cancer cells show altered mechanical phenotypes compared to their healthy counterparts. Atomic Force Microscopy (AFM) is a widely utilized tool to study cell mechanics. These measurements often need skilful users, physical modelling of mechanical properties and expertise in data interpretation. Together with the need to perform many measurements for statistical significance and to probe wide enough areas in tissue structures, the application of machine learning and artificial neural network techniques to automatically classify AFM datasets has received interest recently. We propose the use of self-organizing maps (SOMs) as unsupervised artificial neural network applied to mechanical measurements performed via AFM on epithelial breast cancer cells treated with different substances that affect estrogen receptor signalling. We show changes in mechanical properties due to treatments, as estrogen softened the cells, while resveratrol led to an increase in cell stiffness and viscosity. These data were then used as input for SOMs. Our approach was able to distinguish between estrogen treated, control and resveratrol treated cells in an unsupervised manner. In addition, the maps enabled investigation of the relationship of the input variables.
Mechanical forces are important for cell migration, interaction of cells with their environment, tissue morphogenesis and in various forms of disease1,2,3,4,5,6. In cancer, due to deposition of aligned extracellular matrix proteins solid, stiff tumours can be recognized on the macroscale by palpation and similar techniques7. Interestingly, on the single cell scale, cancer cells have been found to exhibit a softer, more fluid like phenotype8. This change is conserved over many types of cancer tissue origin, ranging from breast, brain, prostate, kidney to lung cancer. Besides other biochemical alterations in cancer cell metabolism, they are thought to have a distinctly reorganized cytoskeletal network which contributes to their mechanical phenotype9. This includes altered actin cytoskeleton dynamics such as large lamellar protrusions and actin-rich micro-spikes, or altered microtubule dynamics leading to amoeboid like increased cell motility. In addition, the softening of cells appears to correlate with the level of aggressiveness, as well as changes in the adhesive and motile properties of these cells10,11,12. Therefore, the mechanical phenotype of cancer cells is thought to be a promising, label-free biomarker13,14,15,16. Different studies have provided evidence for the softening of cancer cells and tissue, as well as enabled the differentiation of cancer progression in tissue sections11,17,18.
Breast cancer is the most commonly diagnosed cancer and the first cause of death from cancer in women worldwide19. A common therapeutic approach in hormone-dependent breast cancer, the most common type of breast cancer, is to target estrogen receptor (ER) signalling with drugs such as tamoxifen or aromatase inhibitors. ER signalling is widely recognized in playing active roles in cytoskeleton, motility and adhesion protein expression with multiple downstream targets20,21. A softening of breast cancer cells with increased aggressiveness has been shown in recent works18,22,23,24. The alterations in the mechanical phenotypes of breast cancer cells treated with drugs that either inhibit or activate ER were studied in vitro25,26.
Mechanical properties of cells and tissue can be measured either by active or passive methods27,28. The first include methods such as parallel plate rheometry29, optical and magnetic tweezers30,31, micropipette aspiration32, optical stretcher33, magnetic twisting cytometry34 and AFM35, while the latter include particle tracking, and spectroscopic methods27. AFM is a widely applied technique, enabling imaging with nanometric precision while measuring mechanical properties at a wide set of forces, strains, and frequencies. In AFM measurements, the probe (e.g., a micrometric particle or a nanometric tip) is brought into contact with the sample and the bending of the cantilever is measured. The bending in relation to the position and stiffness of the cantilever can be converted to determine the deformation of a sample under a given force (force-distance-curves). In addition, force and deformation can be monitored over time (force-time-curves). The mechanical properties of the sample can then be determined by fitting various mechanical models of different complexity to the data36. There is the notion that AFM based single cell and tissue mechanics could be applied in clinical settings to support medical professionals in classification of samples such as cancer progression through evaluation of biopsies18,37.
Recently, machine learning (ML) and neural network algorithms have been applied to surface probe microscopy data38, including image segmentation39, automatic data processing40,41,42,43, cancer cell classification and progression grade evaluation44,45,46,47,48. These approaches were performed in a supervised way, with manual classification of a data set to train the ML algorithms that were applied to classify data. In this work, we apply self-organizing maps (SOMs, also called Kohonen maps) as an artificial neural network approach to analyse mechanical measurements of breast cancer cells49,50,51,52. SOMs are unsupervised artificial neural networks that enable 2D visualization of multidimensional data space while still conserving the topology of the dataset. They have been widely applied in data science51, speech and document processing53, and recently also in life sciences and chemometrics54. One result of SOMs are U-matrix plots that show the similarity/dissimilarity between neurons, enabling post hoc cluster analysis51.
In this work we have used epithelial breast cancer (MCF-7) cells as model for breast cancer cell mechanics and performed stress relaxation measurements in the nuclear regions of cells via AFM using sharp pyramidal tips. We then applied a viscoelastic model to derive mechanical parameters that were used as input layer for a SOM, which was followed by cluster analysis. We show the principal ability of SOMs to differentiate between control, estrogen and resveratrol treated cells, while the SOMs placed control and tamoxifen treated cells in close neighbourhood.
Table 1 shows the values derived from the 5-element Maxwell model fitting for the stress relaxation segments of the cells measured via AFM after the three different treatments (plus the carrier DMSO, as control). Control cells showed moduli in the range of a few hundred Pa and relaxation times of 0.2 and 3.1 s respectively. Treatment with estrogen led to significant softening and less viscous mechanical phenotype, comparable to other published results (reduction of moduli of around 50%)26,55. The treatment with tamoxifen did not lead to significant changes in the mechanical properties of the cells, according to this analysis (increase of 10 to 15 % in moduli and viscosities). Finally, treatment with resveratrol led to a stark increase of moduli of up to 6 times, slightly longer relaxation times and threefold increased viscosities. These data overall agree well with published literature25. Note that some of the parameters showed strong correlations, as moduli reflect overall the elastic properties of the samples and are inverse proportional to indentations, while in the used model, viscosities are proportional to relaxation times and moduli.
Self-organizing maps were trained with the scaled cell viscoelastic properties (centred and divided by the standard deviation) obtained after the different treatments56. The performance of both iterative and batch algorithms was compared using a wide range of map parameters (iteration steps, size and shape of map, learning rates, and distance functions). These calculations were repeated multiple times to account for randomness in the initialisation and training. The maps shown from here on are the results of a batch trained map with 11 times 18 hexagonal nodes, a toroidal topography, a gaussian neighbourhood function, learning rates of (0.05, 0.01), Euclidean distance functions and 1000 iteration steps (map with smallest quantization and topographic error). Figure 1 shows the results of the SOM training.
Results of training the self-organizing map using a batch algorithm. (A) Unified distance matrix plot. The colour scale corresponds to Euclidean distance between nodes. The post hoc clustering was performed using kmeans clustering with 4 centroids. (B) Fan-diagram showing the spatial distribution of the nine different variables on the 2D map. (delta) is indentation, ({E}_{infty }) the equilibrium modulus, ({E}_{1}) and ({E}_{2}) the moduli of the springs in Maxwell arms, ({E}_{inst}) the instantaneous modulus, ({tau }_{1}) and ({tau }_{2}) the relaxation times of the dashpots and the viscosities are ({eta }_{1}) and ({eta }_{2}). (C) Counts plot of the number of observations per node. In grey nodes, zero observations are placed. (D) Mapping of the treatment as input factors. Cluster 1 includes 227 measurements, Cluster 2 69, Cluster 3 195 and Cluster 4 545.
In Fig. 1A the unified distance matrix of the SOM can be seen. The colour scale indicates the Euclidean distance between nodes, blue tones indicate small distances, while red tones indicate large ones. A distinct region with high dissimilarity can be observed at the bottom right of the map (visible as cluster). Note that the map is toroidal, thus reducing edge effects. The dark black lines indicate the post hoc clustering of the nodes. Figure 1B shows the distributions of the single codebook vectors and the determined clusters in greyscales. This analysis can be used to analyse the distribution of the single input variables. Even for the low number of input variables, this analysis can be quite overwhelming, as only regions with significant differences can be seen easily. As an example, the nodes that showed high dissimilarity in Fig. 1A on the right bottom side of the map showed large values of all moduli (large green and cyan fans), while only possessing small indentations of around 1 μm (small dark blue bars) and small relaxation times (red and orange bars) in the fan diagram in Fig. 1B. Measurements with high indentation values (2.5–3 μm) appeared to be placed in vicinity on the top left-hand side of the map. Figure 1C shows the number of measurements that were placed in each node. Finally, Fig. 1D uses the treatment factor (as these are known from the input data) and visualizes where measurements of respectively treated cells were placed on the map. The cluster on the bottom right side of the map is made up solely of resveratrol treated cells, while a large amount of estrogen treated cell measurements are located at the top left side of the map. Note that the number of clusters was defined as 4 to reflect the different treatments.
Tables 2 and 3 show a further analysis of the four clusters. Cluster 1 was predominantly (80%) made of cells from the control group (104) and tamoxifen treated cells (77). Cluster 2 only included resveratrol treated cells. Cluster 3 was a mixture of all cell types, mostly control and tamoxifen treated ones (50% control, 7% E2, 16% RESV and 27% TAM). Finally, cluster 4 was made up mostly of control and estrogen treated cells (39% control, 37% E2, 6% RESV and 18% TAM). With respect to mechanics, cells placed in cluster 1 and 3 showed intermediate stiffnesses and indentations. Cluster 3 represented higher relaxation times than all other clusters. Resveratrol treated cells, located in Cluster 2, presented high stiffness, low indentations, high viscosities, and intermediate relaxation times. The characteristics of cluster 4 were very soft cells, with large indentations, low moduli, intermediate relaxation times and low viscosities. The cluster analysis shows that resveratrol treated cells can be readily distinguished from the other groups, while estrogen treated cells could mostly be found in a cluster together with control and tamoxifen treated cells. Control cells and tamoxifen treated cells co-localised on the maps and the SOMs were not able to distinguish between them. This was expected, as the mechanical properties were very similar and the only used input for the neural network. Additional input parameters would be needed to distinguish between these treatments.
One strength of self-organizing maps is that they allow exploration of interconnectivity of different input variables. Figure 2 shows the component maps of the 9 input variables after training. In addition, the cluster analysis was superimposed on the maps. Note that the moduli and viscosities were plotted in logarithmic scale to allow for better comparison.
Component maps of the nine input parameters after training the SOMs. The colour coding indicates the magnitude of the nodes. The moduli and viscosities are logarithmically scaled to allow for comparison. The black lines show the borders of the four calculated clusters. Grey nodes have no measurements placed in them.
A correlation between all moduli and the viscosities was observed. In addition, the six parameters ({E}_{infty }), ({E}_{1}), ({E}_{2}), ({E}_{inst}), ({eta }_{1}) and ({eta }_{2}) were negatively correlated with the indentation (red colour tones indicating high values while the indentation shows blue colour tones indicating low values). This fits together with the mechanical model, as cells that are softer should show lower moduli, higher fluidity (lower viscosity) and larger deformations. Two distinct regions of deformation can be seen in the component maps: On the left top, a region with high indentation and low stiffness as well as viscosity was found, which corresponded mostly to estrogen treated cells (comparison of Figs. 2 and 1D). On the other hand, an island of low indentation, large moduli and large viscosities was detected on the bottom right part of the maps, which mostly displayed measurements made on resveratrol treated cells. Regarding the relaxation times, a large patch of small relaxation times was identified in the vertical middle region of the maps. Interestingly, these values appeared to correspond to intermediate deformation, moduli and viscosity values. A small patch of long relaxation times was noticed at the bottom left of the maps, which represented intermediate stiffness and deformation values.
In this work we apply self-organizing maps to viscoelastic data of breast cancer cells derived from stress relaxation atomic force microscopy measurements after treatment with different ER interacting drugs. We employ an unsupervised approach and use the output layer of this artificial neural network for cluster analysis. We show that resveratrol treatment leads to stiffening and higher viscosities, estrogen treatment to a softening and fluidization, while tamoxifen treatment does not appear to significantly affect these mechanical properties. We then provide evidence that this type of analysis can potentially be used to classify cells according to treatment in an unsupervised manner.
Cancer cells show altered mechanical phenotypes compared to non-malignant counterparts16. This results among other factors from changes in the mechanics of the environment, tumour hypoxia and aberrant gene expression, altering signalling pathways that affect cell metabolism, motility, and adhesion. Targeting the mechanical adaptation programme of cancer as part of therapy has been proposed, although its wide implications in cell function may lead to undesired side effects. Potential targets include disruption of actin filament assembly and dynamics, inhibition of myosin or targeting Rho/ROCK pathways, among others. Different drugs used as therapy such as paclitaxel, cisplatin, doxorubicin, and 5-fluorouracil, have been investigated for their effects on mechanical properties of cells from various origins. These treatments mostly led to cell stiffening, which contributed to altered cytoskeletal dynamics resembling a reversal of the epithelial to mesenchymal transition57,58,59. Such changes are believed to lead to reduced cell migration, partly inhibiting cancer growth. We provide further evidence that resveratrol leads to stiffening of cells, with novel data regarding the viscoelastic properties, while estrogen leads to a softening and fluidization of the cells. Interestingly, tamoxifen does not appear to significantly alter the mechanics of MCF-7 cells on solid substrates. A more thorough analysis regarding cytoskeletal arrangements and protein expression patterns resulting from such treatments is needed.
In the present approach we have used viscoelastic data derived from stress relaxation measurements using a five-element Maxwell model fitting. We have omitted the use of any further characteristics of force spectroscopy curves on cells (such as hysteresis, work of adhesion, tethers, tension). In addition, we have decided to ignore the measurement data itself, but rather evaluate derived mechanical properties. A critical parameter to consider is the computational time the analysis takes. While the computation of the self-organizing map approach chosen is relatively straightforward, the data processing steps to calculate the viscoelastic properties from raw force-distance-curves is the time intensive step of this analysis (the used framework takes around 5 s for curve processing, fitting, and property calculation per curve). Using simpler models (such as only elastic properties, indentations at a given force, hysteresis, adhesion properties) one can probably reduce the computational costs. In addition, we show that most of the variables used as inputs are correlated. Therefore, a priori reduction of data dimensionality using principal component analysis followed by statistically relevant principal components as input for the SOMs, will also speed up the analysis.
Most of the neural network or machine learning approaches that have been applied to AFM force spectroscopy data on biological materials are based on supervised methods. Recently, three different approaches have been provided in the literature: Using machine learning to classify the quality of curves41, training a machine learning algorithm with curve shapes44, or using curve characteristic properties as inputs for supervised training43,45,60. Such approaches have shown promising results in correctly classifying cancer cells from healthy cells, as well as cancerous tissue. Compared to these approaches, SOMs arguably show strengths in data exploration, enabling the simplification of a multidimensional data space to 2D representations. In addition, SOMs can also be performed supervised. With respect to the ability of this SOM approach to identify the treatment of breast cancer cells, we report mixed results. The addition of further independent variables, including more measurements, and morphological characteristics of the cells will probably help in training a SOM algorithm to successfully classify the different types of cells unsupervised.
Figure 3 shows the overall methodology performed in this study. Briefly, in step 1, cell mechanical properties were quantified using stress relaxation measurements via AFM. A scheme of such analysis can be seen in Fig. 3B. A 5-element Maxwell model was fitted to the stress relaxation curves, resulting in 9 fitting parameters (Fig. 3C). These mechanical parameters were used to train SOMs (Fig. 3D). Finally, these SOMs were used to visualize the result of the artificial neural network. To test the ability of the trained aNNs to distinguish between cells with different mechanical properties, MCF-7 breast cancer cells were treated with different substances used in breast cancer therapy that are known to influence cell mechanical properties.
Methodology of the present study. (A) Cell viscoelastic properties are measured using stress relaxation experiments with AFM. (B) and (C). The resulting stress relaxation segment (force decay over time while keeping the deformation constant) is fitted using a five-element Maxwell model. This results in 9 input variables. (D) A self-organizing map is trained unsupervised. (E). The complex interconnectivity of the data can be represented in 2D maps, allowing for data exploration.
MCF-7 cells were obtained from the American Type Culture Collection (ATCC) and grown as previously described61. Cells were seeded at a concentration of 200,000 cells/mL on plasma-treated glass slides in DMEM medium supplemented with 10% fetal bovine serum (FBS) and 1% penicillin/streptomycin at 37 °C in 5% CO2. Cells were either treated for 48 hours with DMSO (0.05 %), 100 nM estrogen (E2), 50 µM resveratrol (RESV) or 5 µM tamoxifen (TAM). For mechanical measurements, a JPK Nanowizard III with a CellHesion extension was used. Pyramidal cantilevers (DNP-S, B, Bruker), with nominal stiffness of 0.12 N/m, a resonance frequency of 23 kHz in air, an opening angle of 22° and a nominal tip radius of 10 nm were used. Calibration by thermal noise making use of the equipartition theorem was performed for each cantilever62. Measurements were done with a constant approach and retract rate of 5 µm/s, a maximum load of 1 nN, curve lengths of 50 µm, a constant deformation 10 s pause segment at an initial load of 1 nN and a sampling rate of 1024 Hz. Measurements were done in L15 medium at 37 °C and performed in the central region of the cell above the nucleus.
Data was extracted using the JPKSPM software (JPK, Bruker), and all further steps were performed in R. Data pre-processing was done making use of the R afmToolkit, a package for AFM force-distance and force-time analysis developed by our group63,64,65. Briefly, contact and detachment points were calculated, the curves were corrected for their baselines and the sample deformation ((delta )) was determined. In the present analysis, only the stress relaxation segments were considered. Those were fitted with a five-element Maxwell model as
Where (Eleft( t right)) is the relaxation modulus, (E_{infty }) the equilibrium modulus, (E_{1}) and (E_{2}) the moduli of the springs in Maxwell arms and (tau_{1}) and (tau_{2}) the relaxation times of the dashpots. The viscosity (eta_{i}) of the dashpots is defined as
Self-organizing maps were first introduced by Kohonen and have since received more and more attention in wider fields of statistics, chemometrics and life sciences. Their major use is found in data exploration and visualization, as they enable the breakdown of complex data matrices and the interaction of components via 2D similarity representations. In addition, the representation keeps the topological information of the input data. Here we describe the principles of the classical iterative and the batch algorithms. Any reader is referred to the extensive reviews on the theoretical backgrounds, application and usage of SOMs in statistical analysis, life sciences and industry51,52,66. Such applications are wide ranged, including analysis of spectroscopy data67,68, automated feature extraction from microscopy images69, segmentation and grading of tumours or application to omics data70,71.
SOMs are artificial neural networks that use unsupervised training procedures (competition, cooperation, and adaptation). In principle, a (mostly) 2D network of neurons (nodes, codebooks) is set up that is trained to represent the distance structure of the input data as closely as possible (see Fig. 3D). In the case of SOMs, nodes are defined by codebook vectors (weights) that encompass all variables of the input data. Prior to training, the network is initialized either randomly or linearly. One needs to a priori define network properties such as topology (hexagonal or rectangular), architecture of the grid (dimensions, shape, repetitiveness, dynamically growing SOMs), number of neurons (depends on input data, approach, and analysis route). Depending on the type of distance measure applied, input data needs to be normalized to account for different scaling.
In the original iterative approach, the following training steps are performed:
Randomly select an input vector.
Calculate the distance between the input and all neurons in the network (using Euclidean distance after normalization works for most cases).
Determine the best matching unit (BMU) with the lowest distance to the input vector (competition).
Calculate the neighbouring neurons that will also be updated according to the chosen neighbourhood function (gaussian, bubble) (cooperation).
Update the BMU and the neighbours applying learning rates, by moving them closer to the input vector (adaptation).
Repeat steps 1 to 6 until a given number of iterations is reached or the system converges to stability.
Both the neighbourhood distance and the learning rate are decreasing monotonically with each iteration. The learning rate is adjusted depending on the proximity of the neuron to the BMU. The learning phase is split in two main segments, first with a high learning rate (ordering phase) and then with a lower one (convergence phase). In batch SOMs, a similar approach is performed but rather all neurons are compared with all input data.
As input data, we used the 9 properties ((delta , { E}_{infty }, { E}_{1}, { E}_{2}, {E}_{inst}, {tau }_{1}, {tau }_{2}, {eta }_{1}, {eta }_{2})) derived from the fitting of the stress relaxation segments. Outliers were removed and data was normalized by mean value and standard deviation. For SOM calculations, the R package Kohonen was used56. Rectangular, toroidal maps with hexagonal nodes were used. The x-y-ratio of the maps was determined from the relationship of the eigenvalues of the first two principal components of the data. A Gaussian neighbourhood function and Euclidean distance functions were utilized. The learning rate was set to 0.05 in the beginning and decreased to 0.01. The neighbourhood function was defined as exponential decay starting with a value that covered 2/3 of all unit-to-unit distances. For the classical, iterative approaches, 50,000 iterations were performed, the calculations were repeated 50 times and the map with the lowest quantization error was chosen. Similarly, batch SOMs were performed with 1000 iterations and repeated 50 times. The results shown in this manuscript correspond to the batch SOM with the lowest quantization error.
Data is available per request from the corresponding author.
Vogel, V. & Sheetz, M. Local force and geometry sensing regulate cell functions. Nat. Rev. Mol. Cell Biol. 7, 265–275 (2006).
Article CAS PubMed Google Scholar
Vining, K. H. & Mooney, D. J. Mechanical forces direct stem cell behaviour in development and regeneration. Nat. Rev. Mol. Cell Biol. 18, 728–742 (2017).
Article CAS PubMed PubMed Central Google Scholar
Roca-Cusachs, P., Conte, V. & Trepat, X. Quantifying forces in cell biology. Nat. Cell Biol. 19, 742–751 (2017).
Article CAS PubMed Google Scholar
Petridou, N. I., Spiró, Z. & Heisenberg, C.-P. Multiscale force sensing in development. Nat. Cell Biol. 19, 581–588 (2017).
Article CAS PubMed Google Scholar
Huse, M. Mechanical forces in the immune system. Nat. Rev. Immunol. 17, 679–690 (2017).
Article CAS PubMed PubMed Central Google Scholar
Romani, P., Valcarcel-Jimenez, L., Frezza, C. & Dupont, S. Crosstalk between mechanotransduction and metabolism. Nat. Rev. Mol. Cell Biol. 22, 22–38 (2021).
Article CAS PubMed Google Scholar
Winkler, J., Abisoye-Ogunniyan, A., Metcalf, K. J. & Werb, Z. Concepts of extracellular matrix remodelling in tumour progression and metastasis. Nat. Commun. 11, 5120 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Wirtz, D., Konstantopoulos, K. & Searson, P. C. The physics of cancer: the role of physical interactions and mechanical forces in metastasis. Nat. Rev. Cancer 11, 512–522 (2011).
Article CAS PubMed PubMed Central Google Scholar
Fife, C. M., McCarroll, J. A. & Kavallaris, M. Movers and shakers: cell cytoskeleton in cancer metastasis. Br. J. Pharmacol. 171, 5507–5523 (2014).
Article CAS PubMed PubMed Central Google Scholar
Lekka, M. et al. Cancer cell recognition – mechanical phenotype. Micron 43, 1259–1266 (2012).
Article PubMed Google Scholar
Lin, H.-H. et al. Mechanical phenotype of cancer cells: cell softening and loss of stiffness sensing. Oncotarget 6, 20946–20958 (2015).
Article PubMed PubMed Central Google Scholar
Brill-Karniely, Y. et al. Triangular correlation (TrC) between cancer aggressiveness, cell uptake capability, and cell deformability. Sci. Adv. 6, eaax2861 (2022).
Article ADS Google Scholar
Di Carlo, D. A mechanical biomarker of cell state in medicine. J. Lab. Autom. 17, 32–42 (2012).
Article PubMed Google Scholar
Northcott, J. M., Dean, I. S., Mouw, J. K. & Weaver, V. M. Feeling stress: the mechanics of cancer progression and aggression. Front. Cell Dev. Biol. 6, 17 (2018).
Article PubMed PubMed Central Google Scholar
Kozminsky, M. & Sohn, L. L. The promise of single-cell mechanophenotyping for clinical applications. Biomicrofluidics 14, 31301 (2020).
Article CAS Google Scholar
Nguyen, L. T. S., Jacob, M. A. C., Parajón, E. & Robinson, D. N. Cancer as a biophysical disease: targeting the mechanical-adaptability program. Biophys. J. 121, 3573–3585 (2022).
Article ADS CAS PubMed Google Scholar
Lekka, M. Discrimination between normal and cancerous cells using AFM. Bionanoscience 6, 65–80 (2016).
Article PubMed PubMed Central Google Scholar
Plodinec, M. et al. The nanomechanical signature of breast cancer. Nat. Nanotechnol. 7, 757–765 (2012).
Article ADS CAS PubMed Google Scholar
Sung, H. et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249 (2021).
Article PubMed Google Scholar
Bouris, P. et al. Estrogen receptor alpha mediates epithelial to mesenchymal transition, expression of specific matrix effectors and functional properties of breast cancer cells. Matrix Biol. 43, 42–60 (2015).
Article CAS PubMed Google Scholar
Padilla-Rodriguez, M. et al. The actin cytoskeletal architecture of estrogen receptor positive breast cancer cells suppresses invasion. Nat. Commun. 9, 2980 (2018).
Article ADS PubMed PubMed Central Google Scholar
Smolyakov, G. et al. Elasticity, adhesion, and tether extrusion on breast cancer cells provide a signature of their invasive potential. ACS Appl. Mater. Interfaces 8, 27426–27431 (2016).
Article CAS PubMed Google Scholar
Aermes, C., Hayn, A., Fischer, T. & Mierke, C. T. Cell mechanical properties of human breast carcinoma cells depend on temperature. Sci. Rep. 11, 10771 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Li, Q. S., Lee, G. Y. H., Ong, C. N. & Lim, C. T. AFM indentation study of breast cancer cells. Biochem. Biophys. Res. Commun. 374, 609–613 (2008).
Article CAS PubMed Google Scholar
Iturri, J. et al. Resveratrol-induced temporal variation in the mechanical properties of MCF-7 breast cancer cells investigated by atomic force microscopy. Int. J. Mol. Sci. 20, 3275 (2019).
Article CAS PubMed PubMed Central Google Scholar
Bischoff, P. et al. Estrogens determine adherens junction organization and E-cadherin clustering in breast cancer cells via amphiregulin. iScience 23, 101683 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Moeendarbary, E. & Harris, A. R. Cell mechanics: principles, practices, and prospects. Wiley Interdis. Rev. Syst. Biol. Med. 6, 371–388 (2014).
Article Google Scholar
Wu, P.-H. et al. A comparison of methods to assess cell mechanical properties. Nat. Methods 15, 491–498 (2018).
Article CAS PubMed PubMed Central Google Scholar
Bufi, N., Durand-Smet, P. & Asnacios, A. Chapter 11 – Single-cell mechanics: The parallel plates technique. In Methods in Cell Biology (ed. Paluch, E. K.) 187–209 (Academic Press, 2015).
Google Scholar
Arbore, C., Perego, L., Sergides, M. & Capitanio, M. Probing force in living cells with optical tweezers: from single-molecule mechanics to cell mechanotransduction. Biophys. Rev. 11, 765–782 (2019).
Article CAS PubMed PubMed Central Google Scholar
Sarkar, R. & Rybenkov, V. V. A guide to magnetic tweezers and their applications. Front. Phys. 4, 48 (2016).
Article Google Scholar
Hochmuth, R. M. Micropipette aspiration of living cells. J. Biomech. 33, 15–22 (2000).
Article CAS PubMed Google Scholar
Guck, J. et al. The optical stretcher: a novel laser tool to micromanipulate cells. Biophys J 81, 767–784 (2001).
Article CAS PubMed PubMed Central Google Scholar
Puig-De-Morales, M. et al. Measurement of cell microrheology by magnetic twisting cytometry with frequency domain demodulation. J. Appl. Physiol. 91, 1152–1159 (2001).
Article CAS PubMed Google Scholar
Haase, K. & Pelling, A. E. Investigating cell mechanics with atomic force microscopy. J. R. Soc. Interface R. Soc. 12, 20140970 (2015).
Article Google Scholar
Efremov, Y. M., Okajima, T. & Raman, A. Measuring viscoelasticity of soft biological samples using atomic force microscopy. Soft Matter 16, 64–81 (2019).
Article ADS PubMed Google Scholar
Lekka, M. et al. Elasticity of normal and cancerous human bladder cells studied by scanning force microscopy. Eur. Biophys. J. 28, 312–316 (1999).
Article CAS PubMed Google Scholar
Azuri, I., Rosenhek-goldian, I., Regev-rudzki, N., Fantner, G. & Cohen, S. R. The role of convolutional neural networks in scanning probe microscopy : a review. Beilstein J. Nanotechnol. 12, 878–901 (2021).
Article CAS PubMed PubMed Central Google Scholar
Sotres, J., Boyd, H. & Gonzalez-Martinez, J. F. Enabling autonomous scanning probe microscopy imaging of single molecules with deep learning. Nanoscale 13, 9193–9203 (2021).
Article CAS PubMed Google Scholar
Braunsmann, C. & Schäffer, T. E. Note: artificial neural networks for the automated analysis of force map data in atomic force microscopy. Rev. Sci. Instrum. 85, 56104 (2014).
Article ADS Google Scholar
Müller, P. et al. nanite: using machine learning to assess the quality of atomic force microscopy-enabled nano-indentation data. BMC Bioinformat. 20, 465 (2019).
Article Google Scholar
Ilieva, N. I., Galvanetto, N., Allegra, M., Brucale, M. & Laio, A. Automatic classification of single-molecule force spectroscopy traces from heterogeneous samples. Bioinformatics 36, 5014–5020 (2020).
Article CAS PubMed Google Scholar
Sotres, J., Boyd, H. & Gonzalez-Martinez, J. F. Locating critical events in AFM force measurements by means of one-dimensional convolutional neural networks. Sci. Rep. 12, 12995 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Minelli, E. et al. A fully-automated neural network analysis of AFM force-distance curves for cancer tissue diagnosis. Appl. Phys. Lett. 111, 143701 (2017).
Article ADS Google Scholar
Sokolov, I. et al. Noninvasive diagnostic imaging using machine-learning analysis of nanoresolution images of cell surfaces: detection of bladder cancer. Proc. Natl. Acad. Sci. 115, 12920–12925 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Ciasca, G. et al. Efficient spatial sampling for AFM-based cancer diagnostics: a comparison between neural networks and conventional data analysis. Condens. Matter 4, 58 (2019).
Article CAS Google Scholar
Prasad, S. et al. Atomic force microscopy detects the difference in cancer cells of different neoplastic aggressiveness via machine learning. Adv. Nanobiomed. Res. 1, 2000116 (2021).
Article CAS Google Scholar
Tian, Y., Lin, W., Qu, K., Wang, Z. & Zhu, X. Insights into cell classification based on combination of multiple cellular mechanical phenotypes by using machine learning algorithm. J. Mech. Behav. Biomed. Mater. 128, 105097 (2022).
Article PubMed Google Scholar
Kohonen, T. Self-organized formation of topologically correct feature maps. Biol. Cybern 43, 59–69 (1982).
Article MathSciNet MATH Google Scholar
Kohonen, T., Oja, E., Simula, O., Visa, A. & Kangas, J. Engineering applications of the self-organizing map. Proc. IEEE 84, 1358–1384 (1996).
Article Google Scholar
Brereton, R. G. Self organising maps for visualising and modelling. Chem. Cent. J. https://doi.org/10.1186/1752-153X-6-S2-S1 (2012).
Article PubMed PubMed Central Google Scholar
Kohonen, T. Essentials of the self-organizing map. Neural Netw. 37, 52–65 (2013).
Article PubMed Google Scholar
Kohonen, T. et al. Self organization of a massive document collection. IEEE Trans. Neural Netw. 11, 574–585 (2000).
Article CAS PubMed Google Scholar
Marini, F., Bucci, R., Magrì, A. L. & Magrì, A. D. Artificial neural networks in chemometrics: history, examples and perspectives. Microchem. J. 88, 178–185 (2008).
Article CAS Google Scholar
Zbiral, B., Weber, A., Iturri, J., Vivanco, M. D. M. & Toca-Herrera, J. L. Estrogen modulates epithelial breast cancer cell mechanics and cell-to-cell contacts. Materials 14, 2897 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Wehrens, R. & Buydens, L. M. C. Self- and super-organizing maps in R: the Kohonen package. J. Stat. Softw. 21, 1–19 (2007).
Article Google Scholar
Raudenska, M. et al. Cisplatin enhances cell stiffness and decreases invasiveness rate in prostate cancer cells by actin accumulation. Sci. Rep. 9, 1660 (2019).
Article ADS PubMed PubMed Central Google Scholar
Rodríguez-Nieto, M. et al. Viscoelastic properties of doxorubicin-treated HT-29 cancer cells by atomic force microscopy: the fractional Zener model as an optimal viscoelastic model for cells. Biomech. Model. Mechanobiol. 19, 801–813 (2020).
Article PubMed Google Scholar
Kubiak, A., Zieliński, T., Pabijan, J. & Lekka, M. Nanomechanics in monitoring the effectiveness of drugs targeting the cancer cell cytoskeleton. Int. J. Mol. Sci. 21, 8786 (2020).
Article CAS PubMed PubMed Central Google Scholar
Darling, E. M., Topel, M., Zauscher, S., Vail, T. P. & Guilak, F. Viscoelastic properties of human mesenchymally-derived stem cells and primary osteoblasts, chondrocytes, and adipocytes. J. Biomech. 41, 454–464 (2008).
Article PubMed Google Scholar
Iriondo, O. et al. Distinct breast cancer stem/progenitor cell populations require either HIF1α or loss of PHD3 to expand under hypoxic conditions. Oncotarget 6, 31721–31739 (2015).
Article PubMed PubMed Central Google Scholar
Butt, H.-J. & Jaschke, M. Calculation of thermal noise in atomic force microscopy. Nanotechnology 6, 1–7 (1995).
Article ADS Google Scholar
Benítez, R., Moreno-Flores, S., Bolós, V. J. & Toca-Herrera, J. L. A new automatic contact point detection algorithm for AFM force curves. Microsc. Res. Tech 76, 870–876 (2013).
Article PubMed Google Scholar
Benítez, R., Bolós, V. J. & Toca-Herrera, J. L. afmToolkit: an R package for automated AFM force-distance curves analysis. R J. 9, 291–308 (2017).
Article Google Scholar
Weber, A., Benitez, R. & Toca-Herrera, J. L. Measuring (biological) materials mechanics with atomic force microscopy. 4. Determination of viscoelastic cell properties from stress relaxation experiments. Microsc. Res. Tech. 85, 3584–3295 (2022).
Article Google Scholar
Lloyd, G. R., Brereton, R. G. & Duncan, J. C. Self Organising Maps for distinguishing polymer groups using thermal response curves obtained by dynamic mechanical analysis. Analyst 133, 1046–1059 (2008).
Article ADS CAS PubMed Google Scholar
Beckonert, O., Monnerjahn, J., Bonk, U. & Leibfritz, D. Visualizing metabolic changes in breast-cancer tissue using 1H-NMR spectroscopy and self-organizing maps. NMR Biomed. 16, 1–11 (2003).
Article CAS PubMed Google Scholar
Riese, F. M., Keller, S. & Hinz, S. Supervised and semi-supervised self-organizing maps for regression and classification focusing on hyperspectral data. Remote Sens. 12, 1–7 (2020).
Google Scholar
Kang, M.-S., Kim, H.-R. & Kim, M.-H. Cell classification in 3D phase-contrast microscopy images via self-organizing maps. In Advances in Visual Computing: 10th International Symposium, ISVC 2014, Las Vegas, NV, USA, December 8–10, 2014, Proceedings, Part II (eds Bebis, G. et al.) 652–661 (Springer International Publishing, 2014). https://doi.org/10.1007/978-3-319-14364-4_63.
Chapter Google Scholar
Vijayakumar, C., Damayanti, G., Pant, R. & Sreedhar, C. M. Segmentation and grading of brain tumors on apparent diffusion coefficient images using self-organizing maps. Comput. Med. Imag. Graph 31, 473–484 (2007).
Article CAS Google Scholar
Binder, H. et al. Integrated multi-omics maps of lower-grade gliomas. Cancers (Basel) 14, 2797 (2022).
Article CAS PubMed Google Scholar
Download references
The authors want to acknowledge Amsatou Andorfer-Sarr for maintenance of cell culture facilities and Barbara Zbiral for critical discussion of the data and manuscript. This research was partially funded by Austrian Science Fund (FWF), Grant Number P35777-B and by the Spanish Ministry of Science and Innovation MCIN/AEI/10.13039/501100011033 (PID2020-118464RB-I00, CEX2021-001136-S) and Elkartek funding (KK-2020/00008) by the Basque Government (MdMV).
Institute of Biophysics, Department of Bionanosciences, University of Natural Resources and Life Sciences (BOKU), Vienna, Austria
Andreas Weber & José L. Toca-Herrera
CIC bioGUNE, Basque Research and Technology Alliance, BRTA, Technological Park of Bizkaia, Derio, Spain
Maria dM. Vivanco
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
A.W.: Conceptualization, Methodology, Software, Formal Analysis, Investigation, Writing – Original Draft, Writing – review & editing, visualization. M.d.M.: Writing – review & editing, supervision, project administration, funding acquisition. J.L.T.H.: Conceptualization, Writing – Original Draft, Writing – review & editing, visualization, supervision, project administration, funding acquisition.
Correspondence to José L. Toca-Herrera.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Weber, A., Vivanco, M.d. & Toca-Herrera, J.L. Application of self-organizing maps to AFM-based viscoelastic characterization of breast cancer cell mechanics. Sci Rep 13, 3087 (2023). https://doi.org/10.1038/s41598-023-30156-3
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-30156-3
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Scientific Reports (2023)
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing: Cancer newsletter — what matters in cancer research, free to your inbox weekly.