The radial growth phase of malignant melanoma: multi-phase modelling, numerical simulations and linear stability analysis

CD133-Dependent Activation of Phosphoinositide 3-Kinase /AKT/Mammalian Target of Rapamycin Signaling in Melanoma Progression and Drug Resistance

Melanoma frequently harbors genetic alterations in key molecules leading to the aberrant activation of PI3K and its downstream pathways. Although the role of PI3K/AKT/mTOR in melanoma progression and drug resistance is well documented, targeting the PI3K/AKT/mTOR pathway showed less efficiency in clinical trials than might have been expected, since the suppression of the PI3K/mTOR signaling pathway-induced feedback loops is mostly associated with the activation of compensatory pathways such as MAPK/MEK/ERK. Consequently, the development of intrinsic and acquired resistance can occur. As a solid tumor, melanoma is notorious for its heterogeneity. This can be expressed in the form of genetically divergent subpopulations including a small fraction of cancer stem-like cells (CSCs) and non-cancer stem cells (non-CSCs) that make the most of the tumor mass. Like other CSCs, melanoma stem-like cells (MSCs) are characterized by their unique cell surface proteins/stemness markers and aberrant signaling pathways. In addition to its function as a robust marker for stemness properties, CD133 is crucial for the maintenance of stemness properties and drug resistance. Herein, the role of CD133-dependent activation of PI3K/mTOR in the regulation of melanoma progression, drug resistance, and recurrence is reviewed.

3D Bioprinting: An Enabling Technology to Understand Melanoma

Melanoma is a potentially fatal cancer with rising incidence over the last 50 years, associated with enhanced sun exposure and ultraviolet radiation. Its incidence is highest in people of European descent and the ageing population. There are multiple clinical and epidemiological variables affecting melanoma incidence and mortality, such as sex, ethnicity, UV exposure, anatomic site, and age. Although survival has improved in recent years due to advances in targeted and immunotherapies, new understanding of melanoma biology and disease progression is vital to improving clinical outcomes. Efforts to develop three-dimensional human skin equivalent models using biofabrication techniques, such as bioprinting, promise to deliver a better understanding of the complexity of melanoma and associated risk factors. These 3D skin models can be used as a platform for patient specific models and testing therapeutics.

3D Melanoma Cocultures as Improved Models for Nanoparticle-Mediated Delivery of RNA to Tumors

Cancer therapy is an emergent application for mRNA therapeutics. While in tumor immunotherapy, mRNA encoding for tumor-associated antigens is delivered to antigen-presenting cells in spleen and lymph nodes, other therapeutic options benefit from immediate delivery of mRNA nanomedicines directly to the tumor. However, tumor targeting of mRNA therapeutics is still a challenge, since, in addition to delivery of the cargo to the tumor, specifics of the targeted cell type as well as its interplay with the tumor microenvironment are crucial for successful intervention. This study investigated lipoplex nanoparticle-mediated mRNA delivery to spheroid cell culture models of melanoma. Insights into cell-type specific targeting, non-cell-autonomous effects, and penetration capacity in tumor and stroma cells of the mRNA lipoplex nanoparticles were obtained. It was shown that both coculture of different cell types as well as three-dimensional cell growth characteristics can modulate distribution and transfection efficiency of mRNA lipoplex formulations. The results demonstrate that three-dimensional coculture spheroids can provide a valuable surplus of information in comparison to adherent cells. Thus, they may represent in vitro models with enhanced predictivity for the in vivo activity of cancer nanotherapeutics.

Optical coherence elastography to evaluate depth-resolved elasticity of tissue

Skin-elasticity measurements can assist in the clinical diagnosis of skin diseases, which has important clinical significance. Accurately determining the depth-resolved elasticity of superficial biological tissue is an important research direction. This paper presents an optical coherence elastography technique that combines surface acoustic waves and shear waves to obtain the elasticity of multilayer tissue. First, the phase velocity of the high-frequency surface acoustic wave is calculated at the surface of the sample to obtain the Young's modulus of the top layer. Then, the shear wave velocities in the other layers are calculated to obtain their respective Young's moduli. In the bilayer phantom experiment, the maximum error in the elastic estimation of each layer was 2.2%. The results show that the proposed method can accurately evaluate the depth-resolved elasticity of layered tissue-mimicking phantoms, which can potentially expand the clinical applications of elastic wave elastography.

The Impact of Hyaluronan on Tumor Progression in Cutaneous Melanoma

The incidence of cutaneous melanoma is rapidly increasing worldwide. Cutaneous melanoma is an aggressive type of skin cancer, which originates from malignant transformation of pigment producing melanocytes. The main risk factor for melanoma is ultraviolet (UV) radiation, and thus it often arises from highly sun-exposed skin areas and is characterized by a high mutational burden. In addition to melanoma-associated mutations such as BRAF, NRAS, PTEN and cell cycle regulators, the expansion of melanoma is affected by the extracellular matrix surrounding the tumor together with immune cells. In the early phases of the disease, hyaluronan is the major matrix component in cutaneous melanoma microenvironment. It is a high-molecular weight polysaccharide involved in several physiological and pathological processes. Hyaluronan is involved in the inflammatory reactions associated with UV radiation but its role in melanomagenesis is still unclear. Although abundant hyaluronan surrounds epidermal and dermal cells in normal skin and benign nevi, its content is further elevated in dysplastic lesions and local tumors. At this stage hyaluronan matrix may act as a protective barrier against melanoma progression, or alternatively against immune cell attack. While in advanced melanoma, the content of hyaluronan decreases due to altered synthesis and degradation, and this correlates with poor prognosis. This review focuses on hyaluronan matrix in cutaneous melanoma and how the changes in hyaluronan metabolism affect the progression of melanoma.

Nuclear Transport Factor 2 (NTF2) suppresses WM983B metastatic melanoma by modifying cell migration, metastasis, and gene expression

While changes in nuclear structure and organization are frequently observed in cancer cells, relatively little is known about how nuclear architecture impacts cancer progression and pathology. To begin to address this question, we studied Nuclear Transport Factor 2 (NTF2) because its levels decrease during melanoma progression. We show that increasing NTF2 expression in WM983B metastatic melanoma cells reduces cell proliferation and motility while increasing apoptosis. We also demonstrate that increasing NTF2 expression in these cells significantly inhibits metastasis and prolongs survival of mice. NTF2 levels affect the expression and nuclear positioning of a number of genes associated with cell proliferation and migration, and increasing NTF2 expression leads to changes in nuclear size, nuclear lamin A levels, and chromatin organization. Thus, ectopic expression of NTF2 in WM983B metastatic melanoma abrogates phenotypes associated with advanced stage cancer both in vitro and in vivo, concomitantly altering nuclear and chromatin structure and generating a gene expression profile with characteristics of primary melanoma. We propose that NTF2 is a melanoma tumor suppressor and could be a novel therapeutic target to improve health outcomes of melanoma patients.

Oxidative Stress and Autophagy as Key Targets in Melanoma Cell Fate

Melanoma originates from the malignant transformation of melanocytes and is one of the most aggressive forms of cancer. The recent approval of several drugs has increased the chance of survival although a significant subset of patients with metastatic melanoma do not show a long-lasting response to these treatments. The complex cross-talk between oxidative stress and the catabolic process autophagy seems to play a central role in all aspects of melanoma pathophysiology, from initiation to progression and metastasis, including drug resistance. However, determining the fine role of autophagy in cancer death and in response to redox disruption is still a fundamental challenge in order to advance both basic and translational aspects of this field. In order to summarize the interactions among reactive oxygen and nitrogen species, autophagy machinery and proliferation/growth/death/apoptosis/survival, we provide here a narrative review of the preclinical evidence for drugs/treatments that modulate oxidative stress and autophagy in melanoma cells. The significance and the potential for pharmacological targeting (also through multiple and combination approaches) of these two different events, which can contribute independently or simultaneously to the fate of melanoma, may help to define new processes and their interconnections underlying skin cancer biology and unravel new reliable approaches.

Homo and heterobimetallic palladium and platinum complexes bearing μ-diphosphane bridges involved in biological studies

Given the increasing reports of well-defined bimetallic molecular complexes as potential anticancer agents in the last decades, along with the prevalence of platinum in anticancer therapy, we report here a detailed survey of bimetallic platinum and palladium complexes investigated as potential anticancer agents. Specifically, we will concentrate on the synthesis, characterisation and biological (anticancer) studies of a sub-class of these agents, namely homo and heterobimetallic complexes bearing a bridging phosphane ligand of the type: [L_nM₁(μ-R₂P(CH₂)_nPR₂)M₂L_m] (where M₁ is platinum or palladium, M₂ is any other transition metal, R = alkyl or aryl substituents, L_n or L_m are co-ligands, n = 1-6). We will review the in vitro and in vivo activities and any mechanistic anticancer studies of these complexes with a view of trying to delineate patterns in biological activity and structure-activity relationships (SAR). We do not include the review of bimetallic complexes in this class that have not undergone any anticancer testing, nor those that have been involved in other biological investigations unrelated to cancer studies.

Preparation, Drug Treatment, and Immunohistological Analysis of Tri-Culture Spheroid 3D Melanoma-Like Models

Most currently available three-dimensional melanoma models have either focused on simplicity or were optimized for physiological relevance. Accordingly, these paradigms have been either composed of malignant cells only or they were sophisticated human skin equivalents featuring multiple cell types and skin-like organization. Here, an intermediate spheroid-based assay system is presented, which uses tri-cultures of human CCD-1137Sk fibroblasts, HaCaT keratinocytes, and SK-MEL-28 melanoma cells. Being made of cell lines, these spheroids can be reliably reproduced without any special equipment using standard culture procedures, and they feature different aspects of skin and early stage melanoma. Therefore, this kind of model can be useful for lead-compound testing or addressing fundamental principles of early melanoma formation.

An integrative network-driven pipeline for systematic identification of lncRNA-associated regulatory network motifs in metastatic melanoma

BACKGROUND

Melanoma phenotype and the dynamics underlying its progression are determined by a complex interplay between different types of regulatory molecules. In particular, transcription factors (TFs), microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) interact in layers that coalesce into large molecular interaction networks. Our goal here is to study molecules associated with the cross-talk between various network layers, and their impact on tumor progression.

RESULTS

To elucidate their contribution to disease, we developed an integrative computational pipeline to construct and analyze a melanoma network focusing on lncRNAs, their miRNA and protein targets, miRNA target genes, and TFs regulating miRNAs. In the network, we identified three-node regulatory loops each composed of lncRNA, miRNA, and TF. To prioritize these motifs for their role in melanoma progression, we integrated patient-derived RNAseq dataset from TCGA (SKCM) melanoma cohort, using a weighted multi-objective function. We investigated the expression profile of the top-ranked motifs and used them to classify patients into metastatic and non-metastatic phenotypes.

CONCLUSIONS

The results of this study showed that network motif UCA1/AKT1/hsa-miR-125b-1 has the highest prediction accuracy (ACC = 0.88) for discriminating metastatic and non-metastatic melanoma phenotypes. The observation is also confirmed by the progression-free survival analysis where the patient group characterized by the metastatic-type expression profile of the motif suffers a significant reduction in survival. The finding suggests a prognostic value of network motifs for the classification and treatment of melanoma.

Computational models of melanoma

Genes, proteins, or cells influence each other and consequently create patterns, which can be increasingly better observed by experimental biology and medicine. Thereby, descriptive methods of statistics and bioinformatics sharpen and structure our perception. However, additionally considering the interconnectivity between biological elements promises a deeper and more coherent understanding of melanoma. For instance, integrative network-based tools and well-grounded inductive in silico research reveal disease mechanisms, stratify patients, and support treatment individualization. This review gives an overview of different modeling techniques beyond statistics, shows how different strategies align with the respective medical biology, and identifies possible areas of new computational melanoma research.

Uncertainty quantification on a spatial Markov-chain model for the progression of skin cancer

A spatial Markov-chain model is formulated for the progression of skin cancer. The model is based on the division of the computational domain into nodal points, that can be in a binary state: either in 'cancer state' or in 'non-cancer state'. The model assigns probabilities for the non-reversible transition from 'non-cancer' state to the 'cancer state' that depend on the states of the neighbouring nodes. The likelihood of transition further depends on the life burden intensity of the UV-rays that the skin is exposed to. The probabilistic nature of the process and the uncertainty in the input data is assessed by the use of Monte Carlo simulations. A good fit between experiments on mice and our model has been obtained.

Effect of Mechanical Compression on Invasion Process of Malignant Melanoma Using In Vitro Three-Dimensional Cell Culture Device

Malignant melanoma in the plantar surface of the foot is subjected to various mechanical stimuli generated by daily human activity such as walking. Some studies have reported that mechanical compression affects the development and progression of melanoma. However, little is known about how mechanical compression affects the behavior of malignant melanoma cells in a physiological condition due to the complexity of the invasion mechanisms. In this study, we developed an in vitro three-dimensional cell culture device using microporous membrane in order to evaluate the effects of mechanical compression on the invasion process of malignant melanoma. Our results suggest that the invasion of melanoma cells under the compressive stress for 8 h of culture was promoted with the elongation of F-actin filaments compared to control groups, whereas there was no significant difference between both groups at 32 h of culture, with increasing cell death associated with promoting melanin synthesis. The results of this study contribute to the elucidation of the invasion mechanisms of malignant melanoma caused by mechanical stimulation.

Stability and Roughness of Interfaces in Mechanically Regulated Tissues

Cell division and death can be regulated by the mechanical forces within a tissue. We study the consequences for the stability and roughness of a propagating interface by analyzing a model of mechanically regulated tissue growth in the regime of small driving forces. For an interface driven by homeostatic pressure imbalance or leader-cell motility, long and intermediate-wavelength instabilities arise, depending, respectively, on an effective viscosity of cell number change, and on substrate friction. A further mechanism depends on the strength of directed motility forces acting in the bulk. We analyze the fluctuations of a stable interface subjected to cell-level stochasticity, and find that mechanical feedback can help preserve reproducibility at the tissue scale. Our results elucidate mechanisms that could be important for orderly interface motion in developing tissues.

Nonlinear studies of tumor morphological stability using a two-fluid flow model

We consider the nonlinear dynamics of an avascular tumor at the tissue scale using a two-fluid flow Stokes model, where the viscosity of the tumor and host microenvironment may be different. The viscosities reflect the combined properties of cell and extracellular matrix mixtures. We perform a linear morphological stability analysis of the tumors, and we investigate the role of nonlinearity using boundary-integral simulations in two dimensions. The tumor is non-necrotic, although cell death may occur through apoptosis. We demonstrate that tumor evolution is regulated by a reduced set of nondimensional parameters that characterize apoptosis, cell-cell/cell-extracellular matrix adhesion, vascularization and the ratio of tumor and host viscosities. A novel reformulation of the equations enables the use of standard boundary integral techniques to solve the equations numerically. Nonlinear simulation results are consistent with linear predictions for nearly circular tumors. As perturbations develop and grow, the linear and nonlinear results deviate and linear theory tends to underpredict the growth of perturbations. Simulations reveal two basic types of tumor shapes, depending on the viscosities of the tumor and microenvironment. When the tumor is more viscous than its environment, the tumors tend to develop invasive fingers and a branched-like structure. As the relative ratio of the tumor and host viscosities decreases, the tumors tend to grow with a more compact shape and develop complex invaginations of healthy regions that may become encapsulated in the tumor interior. Although our model utilizes a simplified description of the tumor and host biomechanics, our results are consistent with experiments in a variety of tumor types that suggest that there is a positive correlation between tumor stiffness and tumor aggressiveness.

The influence of soluble fragments of extracellular matrix (ECM) on tumor growth and morphology

A major challenge in matrix-metalloproteinase (MMP) target validation and MMP-inhibitor-drug development for anti-cancer clinical trials is to better understand their complex roles (often competing with each other) in tumor progression. While there is extensive research on the growth-promoting effects of MMPs, the growth-inhibiting effects of MMPs has not been investigated thoroughly. So we develop a continuum model of tumor growth and invasion including chemotaxis and haptotaxis in order to examine the complex interaction between the tumor and its host microenvironment and to explore the inhibiting influence of the gradients of soluble fragments of extracellular matrix (ECM) density on tumor growth and morphology. Previously, it was shown both computationally (in one spatial dimension) and experimentally that the chemotactic pull due to soluble ECM gradients is anti-invasive, contrary to the traditional view of the role of chemotaxis in malignant invasion [1]. With two-dimensional numerical simulation and using a level set based tumor-host interface capturing method, we examine the effects of chemotaxis on the progression and morphology of a tumor growing in nutrient-rich and nutrient-poor microenvironments which was not investigated before. In particular we examine how the geometry of the growing tumor is affected when placed in different environments. We also investigate the effects of varying ECM degradation rate, the production rate of matrix degrading enzymes (MDE), and the conversion of ECM into soluble ECM. We find that chemotaxis due to ECM-fragment gradients strongly influences tumor growth and morphology, and that the instabilities caused by tumor cell proliferation and haptotactic movements can be prevented if chemotaxis is sufficiently strong. The influence of chemotaxis and the above factors on tumor growth and morphology are found to be more prominent in nutrient-poor environments than in nutrient-rich environments. So we extend our investigations of these antinvasive chemotactic influences by examining the effects of cell-cell and cell-ECM adhesion and low proliferation rate for tumors growing in low-nutrient environments. We find that as the extent of chemotaxis increases, the effects of adhesion on tumor growth and shape become negligible. Under conditions of low cell mitosis, chemotaxis may cause the tumor to shrink, as the extent of chemotaxis increases. Both stable and unstable tumor shrinkage are predicted by our model. Unexpectedly, in some cases chemotaxis may contribute toward developing instability where haptotaxis alone induces stable growth.

A nonlinear poroelastic theory of solid tumors with glycosaminoglycan swelling

Mechanics plays a crucial role in the growth, development, and therapeutics of tumors. In this paper, a nonlinear poroelastic theory is established to describe the mechanical behaviors of solid tumors. The free-swollen state of a tumor is chosen as the reference state, which enables us to avoid pursuing a dry and stress-free state that is hard to achieve for living tissues. Our results reveal that the compression resistance of a tumor is primarily attributed to glycosaminoglycan (GAG) swelling, and the compactness of cell aggregates is found to affect tumor consolidation. Over-expressed GAGs and dense cell aggregates can stiffen the tumor, a remodeling mechanism that makes the tumor with higher elastic modulus than its surrounding host tissues. Glycosaminoglycan chains also influence the transient mechanical response of the tumor by modulating the tissue permeability. The theoretical results show good agreement with relevant experimental observations. This study may not only deepen our understanding of tumorigenesis but also provide cues for developing novel anticancer strategies.

Incompressible Limit of a Mechanical Model for Tissue Growth with Non-Overlapping Constraint

A mathematical model for tissue growth is considered. This model describes the dynamics of the density of cells due to pressure forces and proliferation. It is known that such cell population model converges at the incompressible limit towards a Hele-Shaw type free boundary problem. The novelty of this work is to impose a non-overlapping constraint. This constraint is important to be satisfied in many applications. One way to guarantee this non-overlapping constraint is to choose a singular pressure law. The aim of this paper is to prove that, although the pressure law has a singularity, the incompressible limit leads to the same Hele-Shaw free boundary problem.

Tissue-scale, personalized modeling and simulation of prostate cancer growth

Recently, mathematical modeling and simulation of diseases and their treatments have enabled the prediction of clinical outcomes and the design of optimal therapies on a personalized (i.e., patient-specific) basis. This new trend in medical research has been termed "predictive medicine." Prostate cancer (PCa) is a major health problem and an ideal candidate to explore tissue-scale, personalized modeling of cancer growth for two main reasons: First, it is a small organ, and, second, tumor growth can be estimated by measuring serum prostate-specific antigen (PSA, a PCa biomarker in blood), which may enable in vivo validation. In this paper, we present a simple continuous model that reproduces the growth patterns of PCa. We use the phase-field method to account for the transformation of healthy cells to cancer cells and use diffusion-reaction equations to compute nutrient consumption and PSA production. To accurately and efficiently compute tumor growth, our simulations leverage isogeometric analysis (IGA). Our model is shown to reproduce a known shape instability from a spheroidal pattern to fingered growth. Results of our computations indicate that such shift is a tumor response to escape starvation, hypoxia, and, eventually, necrosis. Thus, branching enables the tumor to minimize the distance from inner cells to external nutrients, contributing to cancer survival and further development. We have also used our model to perform tissue-scale, personalized simulation of a PCa patient, based on prostatic anatomy extracted from computed tomography images. This simulation shows tumor progression similar to that seen in clinical practice.

Adhesion-induced instabilities and pattern formation in thin films of elastomers and gels

A hydrostatically stressed soft elastic film circumvents the imposed constraint by undergoing a morphological instability, the wavelength of which is dictated by the minimization of the surface and the elastic strain energies of the film. While for a single film, the wavelength is entirely dependent on its thickness, a co-operative energy minimization dictates that the wavelength depends on both the elastic moduli and thicknesses of two contacting films. The wavelength can also depend on the material properties of a film if its surface tension has a pronounced effect in comparison to its elasticity. When such a confined film is subjected to a continually increasing normal displacement, the morphological patterns evolve into cracks, which, in turn, govern the adhesive fracture behavior of the interface. While, in general, the thickness provides the relevant length scale underlying the well-known Griffith-Kendall criterion of debonding of a rigid disc from a confined film, it is modified non-trivially by the elasto-capillary number for an ultra-soft film. Depending upon the degree of confinement and the spatial distribution of external stress, various analogs of the canonical instability patterns in liquid systems can also be reproduced with thin confined elastic films.

Metastatic melanoma cells evade immune detection by silencing STAT1

Transcriptional activation of major histocompatibility complex (MHC) I and II molecules by the cytokine, interferon γ (IFN-γ), is a key step in cell-mediated immunity against pathogens and tumors. Recent evidence suggests that suppression of MHC I and II expression on multiple tumor types plays important roles in tumor immunoevasion. One such tumor is malignant melanoma, a leading cause of skin cancer-related deaths. Despite growing awareness of MHC expression defects, the molecular mechanisms by which melanoma cells suppress MHC and escape from immune-mediated elimination remain unknown. Here, we analyze the dysregulation of the Janus kinase (JAK)/STAT pathway and its role in the suppression of MHC II in melanoma cell lines at the radial growth phase (RGP), the vertical growth phase (VGP) and the metastatic phase (MET). While RGP and VGP cells both express MHC II, MET cells lack not only MHC II, but also the critical transcription factors, interferon response factor (IRF) 1 and its upstream activator, signal transducer and activator of transcription 1 (STAT1). Suppression of STAT1 in vitro was also observed in patient tumor samples, suggesting STAT1 silencing as a global mechanism of MHC II suppression and immunoevasion.

Patterns in melanocytic lesions: impact of the geometry on growth and transport inside the epidermis

In glabrous skin, nevi and melanomas exhibit pigmented stripes during clinical dermoscopic examination. They find their origin in the basal layer geometry which periodically exhibits ridges, alternatively large (limiting ridges) and thin (intermediate ridges). However, nevus and melanoma lesions differ by the localization of the pigmented stripes along furrows or ridges of the epidermis surface. Here, we propose a biomechanical model of avascular tumour growth which takes into account this specific geometry in the epidermis where both kinds of lesions first appear. Simulations show a periodic distribution of tumour cells inside the lesion, with a global contour stretched out along the ridges. In order to be as close as possible to clinical observations, we also consider the melanin transport by the keratinocytes. Our simulations show that reasonable assumptions on melanocytic cell repartition in the ridges favour the limiting ridges of the basal compared with the intermediate ones in agreement with nevus observations but not really with melanomas. It raises the question of cell aggregation and repartition of melanocytic cells in acral melanomas and requires further biological studies of these cells in situ.

Wrinkling pattern evolution of cylindrical biological tissues with differential growth

Three-dimensional surface wrinkling of soft cylindrical tissues induced by differential growth is explored. Differential volumetric growth can cause their morphological stability, leading to the formation of hexagonal and labyrinth wrinkles. During postbuckling, multiple bifurcations and morphological transitions may occur as a consequence of continuous growth in the surface layer. The physical mechanisms underpinning the morphological evolution are examined from the viewpoint of energy. Surface curvature is found to play a regulatory role in the pattern evolution. This study may not only help understand the morphogenesis of soft biological tissues, but also inspire novel routes for creating desired surface patterns of soft materials.

A tumor growth model with deformable ECM

Existing tumor growth models based on fluid analogy for the cells do not generally include the extracellular matrix (ECM), or if present, take it as rigid. The three-fluid model originally proposed by the authors and comprising tumor cells (TC), host cells (HC), interstitial fluid (IF) and an ECM, considered up to now only a rigid ECM in the applications. This limitation is here relaxed and the deformability of the ECM is investigated in detail. The ECM is modeled as a porous solid matrix with Green-elastic and elasto-visco-plastic material behavior within a large strain approach. Jauman and Truesdell objective stress measures are adopted together with the deformation rate tensor. Numerical results are first compared with those of a reference experiment of a multicellular tumor spheroid (MTS) growing in vitro, then three different tumor cases are studied: growth of an MTS in a decellularized ECM, growth of a spheroid in the presence of host cells and growth of a melanoma. The influence of the stiffness of the ECM is evidenced and comparison with the case of a rigid ECM is made. The processes in a deformable ECM are more rapid than in a rigid ECM and the obtained growth pattern differs. The reasons for this are due to the changes in porosity induced by the tumor growth. These changes are inhibited in a rigid ECM. This enhanced computational model emphasizes the importance of properly characterizing the biomechanical behavior of the malignant mass in all its components to correctly predict its temporal and spatial pattern evolution.

Laser induced surface acoustic wave combined with phase sensitive optical coherence tomography for superficial tissue characterization: a solution for practical application

Mechanical properties are important parameters that can be used to assess the physiologic conditions of biologic tissue. Measurements and mapping of tissue mechanical properties can aid in the diagnosis, characterisation and treatment of diseases. As a non-invasive, non-destructive and non-contact method, laser induced surface acoustic waves (SAWs) have potential to accurately characterise tissue elastic properties. However, challenge still exists when the laser is directly applied to the tissue because of potential heat generation due to laser energy deposition. This paper focuses on the thermal effect of the laser induced SAW on the tissue target and provides an alternate solution to facilitate its application in clinic environment. The solution proposed is to apply a thin agar membrane as surface shield to protect the tissue. Transient thermal analysis is developed and verified by experiments to study the effects of the high energy Nd:YAG laser pulse on the surface shield. The approach is then verified by measuring the mechanical property of skin in a Thiel mouse model. The results demonstrate a useful step toward the practical application of laser induced SAW method for measuring real elasticity of normal and diseased tissues in dermatology and other surface epithelia.

Morphology of melanocytic lesions in situ

Melanoma is a solid tumour with its own specificity from the biological and morphological viewpoint. On one hand, numerous mutations are already known affecting different pathways. They usually concern proliferation rate, apoptosis, cell senescence and cell behaviour. On the other hand, several visual criteria at the tissue level are used by physicians in order to diagnose skin lesions. Nevertheless, the mechanisms between the changes from the mutations at the cell level to the morphology exhibited at the tissue level are still not fully understood. Using physical tools, we develop a simple model. We demonstrate analytically that it contains the necessary ingredients to understand several specificities of melanoma such as the presence of microstructures inside a skin lesion or the absence of a necrotic core. We also explain the importance of senescence for growth arrest in benign skin lesions. Thanks to numerical simulations, we successfully compare this model to biological data.

Multiple types of data are required to identify the mechanisms influencing the spatial expansion of melanoma cell colonies

BACKGROUND

The expansion of cell colonies is driven by a delicate balance of several mechanisms including cell motility, cell-to-cell adhesion and cell proliferation. New approaches that can be used to independently identify and quantify the role of each mechanism will help us understand how each mechanism contributes to the expansion process. Standard mathematical modelling approaches to describe such cell colony expansion typically neglect cell-to-cell adhesion, despite the fact that cell-to-cell adhesion is thought to play an important role.

RESULTS

We use a combined experimental and mathematical modelling approach to determine the cell diffusivity, D, cell-to-cell adhesion strength, q, and cell proliferation rate, λ, in an expanding colony of MM127 melanoma cells. Using a circular barrier assay, we extract several types of experimental data and use a mathematical model to independently estimate D, q and λ. In our first set of experiments, we suppress cell proliferation and analyse three different types of data to estimate D and q. We find that standard types of data, such as the area enclosed by the leading edge of the expanding colony and more detailed cell density profiles throughout the expanding colony, does not provide sufficient information to uniquely identify D and q. We find that additional data relating to the degree of cell-to-cell clustering is required to provide independent estimates of q, and in turn D. In our second set of experiments, where proliferation is not suppressed, we use data describing temporal changes in cell density to determine the cell proliferation rate. In summary, we find that our experiments are best described using the range D=161-243μm2 hour-1, q=0.3-0.5 (low to moderate strength) and λ=0.0305-0.0398 hour-1, and with these parameters we can accurately predict the temporal variations in the spatial extent and cell density profile throughout the expanding melanoma cell colony.

CONCLUSIONS

Our systematic approach to identify the cell diffusivity, cell-to-cell adhesion strength and cell proliferation rate highlights the importance of integrating multiple types of data to accurately quantify the factors influencing the spatial expansion of melanoma cell colonies.

Senescent fibroblasts in melanoma initiation and progression: an integrated theoretical, experimental, and clinical approach

We present an integrated study to understand the key role of senescent fibroblasts in driving melanoma progression. Based on the hybrid cellular automata paradigm, we developed an in silico model of normal skin. The model focuses on key cellular and microenvironmental variables that regulate interactions among keratinocytes, melanocytes, and fibroblasts, key components of the skin. The model recapitulates normal skin structure and is robust enough to withstand physical as well as biochemical perturbations. Furthermore, the model predicted the important role of the skin microenvironment in melanoma initiation and progression. Our in vitro experiments showed that dermal fibroblasts, which are an important source of growth factors in the skin, adopt a secretory phenotype that facilitates cancer cell growth and invasion when they become senescent. Our coculture experiments showed that the senescent fibroblasts promoted the growth of nontumorigenic melanoma cells and enhanced the invasion of advanced melanoma cells. Motivated by these experimental results, we incorporated senescent fibroblasts into our model and showed that senescent fibroblasts transform the skin microenvironment and subsequently change the skin architecture by enhancing the growth and invasion of normal melanocytes. The interaction between senescent fibroblasts and the early-stage melanoma cells leads to melanoma initiation and progression. Of microenvironmental factors that senescent fibroblasts produce, proteases are shown to be one of the key contributing factors that promoted melanoma development from our simulations. Although not a direct validation, we also observed increased proteolytic activity in stromal fields adjacent to melanoma lesions in human histology. This leads us to the conclusion that senescent fibroblasts may create a prooncogenic skin microenvironment that cooperates with mutant melanocytes to drive melanoma initiation and progression and should therefore be considered as a potential future therapeutic target. Interestingly, our simulations to test the effects of a stroma-targeting therapy that negates the influence of proteolytic activity showed that the treatment could be effective in delaying melanoma initiation and progression.

Evaluating elastic properties of heterogeneous soft tissue by surface acoustic waves detected by phase-sensitive optical coherence tomography

The combined use of surface acoustic wave (SAW) and phase-sensitive optical coherence tomography (PhS-OCT) is useful to evaluate the elasticity of layered biological tissues, such as normal skin. However, the pathological tissue is often originated locally, leading to the alternation of mechanical properties along both axial and lateral directions. We present a feasibility study on whether the SAW technique is sensitive to detect the alternation of mechanical property along the lateral direction within tissue, which is important for clinical utility of this technique to localize diseased tissue. Experiments are carried out on purposely designed tissue phantoms and ex vivo chicken breast samples, simulating the localized change of elasticity. A PhS-OCT system is employed not only to provide the ultra-high sensitive measurement of the generated surface waves on the tissue surface, but also to provide the real time imaging of the tissue to assist the elasticity evaluation of the heterogeneous tissue. The experimental results demonstrate that with PhS-OCT used as a pressure sensor, the SAW is highly sensitive to the elasticity change of the specimen in both vertical and lateral directions with a sensing depth of ∼5 mm with our current system setup, thus promising its useful clinical applications where the quantitative elasticity of localized skin diseases is needed to aid in diagnosis and treatment.

Thigmotropism of malignant melanoma cells

During malignant melanoma (MM) progression including incipient metastasis, neoplastic cells follow some specific migration paths inside the skin. In particular, they progress along the dermoepidermal basement membrane, the hair follicles, the sweat gland apparatus, nerves, and the near perivascular space. These features evoke the thigmotropism phenomenon defined as a contact-sensing growth of cells. This process is likely connected to modulation in cell tensegrity (control of the cell shape). These specifically located paucicellular aggregates of MM cells do not appear to be involved in the tumorigenic growth phase, but rather they participate in the so-called "accretive" growth model. These MM cell collections are often part of the primary neoplasm, but they may, however, correspond to MM micrometastases and predict further local overt metastasis spread.

Morphological changes in early melanoma development: influence of nutrients, growth inhibitors and cell-adhesion mechanisms

Current diagnostic methods for skin cancers are based on some morphological characteristics of the pigmented skin lesions, including the geometry of their contour. The aim of this article is to model the early growth of melanoma accounting for the biomechanical characteristics of the tumor micro-environment, and evaluating their influence on the tumor morphology and its evolution. The spatial distribution of tumor cells and diffusing molecules are explicitly described in a three-dimensional multiphase model, which incorporates general cell-to-cell mechanical interactions, a dependence of cell proliferation on contact inhibition, as well as a local diffusion of nutrients and inhibiting molecules. A two-dimensional model is derived in a lubrication limit accounting for the thin geometry of the epidermis. First, the dynamical and spatial properties of planar and circular tumor fronts are studied, with both numerical and analytical techniques. A WKB method is then developed in order to analyze the solution of the governing partial differential equations and to derive the threshold conditions for a contour instability of the growing tumor. A control parameter and a critical wavelength are identified, showing that high cell proliferation, high cell adhesion, large tumor radius and slow tumor growth correlate with the occurrence of a contour instability. Finally, comparing the theoretical results with a large amount of clinical data we show that our predictions describe accurately both the morphology of melanoma observed in vivo and its variations with the tumor growth rate. This study represents a fundamental step to understand more complex microstructural patterns observed during skin tumor growth. Its results have important implications for the improvement of the diagnostic methods for melanoma, possibly driving progress towards a personalized screening.

Contour instabilities in early tumor growth models

Recent tumor growth models are often based on the multiphase mixture framework. Using bifurcation theory techniques, we show that such models can give contour instabilities. Restricting to a simplified but realistic version of such models, with an elastic cell-to-cell interaction and a growth rate dependent on diffusing nutrients, we prove that the tumor cell concentration at the border acts as a control parameter inducing a bifurcation with loss of the circular symmetry. We show that the finite wavelength at threshold has the size of the proliferating peritumoral zone. We apply our predictions to melanoma growth since contour instabilities are crucial for early diagnosis. Given the generality of the equations, other relevant applications can be envisaged for solving problems of tissue growth and remodeling.