# Morphology of melanocytic lesions in situ

Thibaut Balois, Martine Ben Amar

► To cite this version:

Thibaut Balois, Martine Ben Amar. Morphology of melanocytic lesions in situ. Scientific Reports, 2014, 4, pp.3622. 10.1038/srep03622 . hal-01359311

HAL Id: hal-01359311

https://hal.sorbonne-universite.fr/hal-01359311

Submitted on 2 Sep 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Morphology of melanocytic lesions in situ

Thibaut Balois¹,² & Martine Ben Amar¹,²

SUBJECT AREAS:
BIOLOGICAL PHYSICS
PHYSICS

¹Laboratoire de Physique Statistique, Ecole Normale Supérieure, UPMC Univ Paris 06, Université Paris Diderot, CNRS, 24 rue Lhomond, 75005 Paris, France, ²Institut Universitaire de Cancérologie, Faculté de médecine, Université Pierre et Marie Curie-Paris 6, 91 boulevard de l'hôpital, 75013 Paris, France.

Received
3 September 2013

Accepted
25 November 2013

Published
10 January 2014

Correspondence and
requests for materials
should be addressed to
M.B.A. (benamar@lps.ens.fr)

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.

Melanoma is a rare form of skin cancer, less than 5% of the cases, but the most lethal and more than 75% of deaths¹. However, recent drug discoveries and targeted therapies indicate evidence for durable remissions but no guaranty can be given yet on drug resistance². It is why early detection and diagnosis remain of utmost importance³ and the evolution of melanocytic lesion morphology remains an important diagnostic tool.

In order to distinguish benign lesions, nevi, from melanomas, physicians have several visual features at their disposal, the most famous one being the ABCDE criteria: Asymmetry, Border, Color, Diameter, Evolution³. They are based on medical statistical studies coupled to the malignancy probability of a skin lesion. Nevertheless, the mechanism of their apparition and evolution remains poorly understood.

Since the work of Greenspan in 1976⁴, several approaches have been developed to model tumour growth⁵, continuum models at the tissue scale⁶, discrete models at the cell scale⁷,⁸, and hybrid continuum-discrete models⁹. But, in spite of their complexity, none of these models take into account the specificity of melanoma. Even in the avascular growth phase, most of them do not exhibit a necrotic core. Introducing only the minimal set of physical parameters, we develop a simple model that can explain how changes at the cell level result in melanoma morphology. We aim to justify the absence of a necrotic core despite the absence of angiogenesis. But first we need to recall briefly the bio-structure and the biology of the skin.

Skin is composed of three main layers, the hypodermis, the dermis and the epidermis. In this work, we focus on skin cancers located in the upper layer of the skin, the epidermis, where tumors can easily be seen and are well referenced. Beneath the epidermis, the dermis contains blood and lymphatic networks providing nutrient supply. It is separated from the epidermis by a basement membrane composed of fibers like collagen. This basement membrane plays the role of a barrier between the dermis and the epidermis, but also the support necessary for cell mitosis. The epidermis is mainly composed of keratinocytes (~90%). They proliferate on the basal layer, and progressively migrate to the skin surface creating a protective barrier called the stratum corneum (see Fig. 1). In order to protect the body from UV damages, melanocytes (~5-10% of epidermis cells), also located in the basal layer of the epidermis, strongly bind to several keratinocytes (~35¹⁰) delivering them melanin, an UV absorbing molecule. Melanocytes are bound to the basement membrane via collagen IV and laminin-mediated adhesion, and to keratinocytes via E-cadherin and Desmoglein (see the Supplementary). In order to divide, melanocytes need to decouple from the basement membrane and from the keratinocytes, and finally they retract their dendrites. After dividing, they migrate before re-binding to the matrix and to keratinocytes.

Through mutations, melanocytic cells change their behaviour resulting in different types of lesions. The most benign one is the nevus, also known as “regular mole”. On the contrary, the dysplastic nevus known as atypical mole display most of the visual characteristics of a melanoma (irregular border, asymmetry) but they are also benign, because of an extension arrest through a senescence process. When a lesion becomes malignant, it is called a melanoma. Usually melanomas undergo different phases. In the radial growth phase (RGP) melanomas grow
<figure><img src="image_3.png" /><figcaption>Figure 1 | Schematic representation of the two outmost layers of the skin, the dermis and the epidermis.</figcaption></figure>

horizontally and are restricted to the epidermis, then in the vertical growth phase (VGP) melanomas invade the dermis rapidly and drastically lowering the survival rate¹¹. (see Fig. 2).

Here, we focus on growth restricted to the epidermis, nevi and melanoma in situ or RGP melanoma. In this phase, melanomas are treated by excision with a very good survival rate (more than 99% of cases¹¹), the metastasis spreading being unlikely. At this point, growth is avascular, and even though angiosis does not occur, the peculiarity of melanoma is the absence of necrotic core during the horizontal development: this is contrary to other solid tumours and to usual modelling of avascular tumour growth where a necrotic core is commonly found. Nonetheless, a necrotic core can appear during the vertical growth of the lesion. For instance, nodular melanomas (~15% of cutaneous melanoma) grow vertically from the beginning and present a necrotic core¹².

Focusing on these main differences and cell micro-environment, we now construct a model in order to explain the different features of melanoma, and its different stages of progression as seen in Fig. 2.

## Results

### Modeling the morphogenesis of melanocytic lesion at the tissue level.
We present here a continuous model of melanoma growth as a two phase mixture, a cancerous phase, $\phi_c$, containing the proliferating melanocytic cells and a second healthy phase containing the interstitial fluid and the keratinocytes, $\phi_h = 1 - \phi_c$. Both phases contain nutrients ensuring cell development. Hereafter, we explicit

<figure><img src="image_4.png" /><figcaption>Figure 2 | Pigmented lesions given their malignancies. The nevi are senescent and benign. The RGP melanomas are immortal but restricted to the epidermis, and the VGP are invasive. This drawing was inspired by¹³.</figcaption></figure>

the equations governing each population, beginning by the nutrient equations.

Nutrients are small molecules which diffuse inside the tissues and their concentration $n$ satisfies a diffusion equation with consumption. Many kinds of nutrients intervene to maintain the homeostasis and growth of the tumour such as glucose³⁶, oxygen³⁴,³⁵ and growth factors. They satisfy the same equation with different biochemical factors. Only oxygen, the slowest one, that is the one with the longest consumption time will control the growth dynamic. Due to skin hypoxia (small concentration of oxygen in the skin), consumption is described by a linearized Michaelis-Menten law, and we get:

$$0 = D_n \Delta n - \delta_n \phi_c n - \kappa (1 - \phi_c) n + S_n (n_s - n) \quad (1)$$

where $\Delta$ is the 2D Laplacian (see the Supplementary), $n$ and $\phi_c$ are now the concentration averaged along the epidermis depth, $n_s$ is a typical nutrient concentration necessary to maintain the homeostatic state of the skin. $S_n(n_s - n)$ is the orthogonal flux of nutrients which is generated by the disequilibrium from the homeostatic state, coming from the dermis (see Fig. 3). This flux regulates the oxygen consumption insuring the nutrient provision and so the natural development of cells on the basement membrane in healthy tissue, including keratinocytes. It is mainly responsible for the absence of a necrotic core inside a growing melanoma in situ. Let us consider now the cancerous cell equation.

Assuming the same density for the two phases, $\phi_c$ and $\phi_h$, the mass balance equation can be written for each phase:

$$\begin{cases} \partial_t \phi_c + \nabla \cdot (\phi_c \mathbf{v}_c) = \Gamma(\phi_c, n) \\ \partial_t \phi_h + \nabla \cdot (\phi_h \mathbf{v}_h) = \Gamma(\phi_h, n) \end{cases} \quad (2)$$

where $\Gamma$ is the mass exchange rate between each phase, and $\Gamma(\phi_c, n) = -\Gamma(\phi_h, n)$. $\mathbf{v}_c$ and $\mathbf{v}_h$ are the velocities of the cancerous and healthy phase. The closure condition imposes the proliferation rate of each phase to be the opposite of each other. We choose the simplest form for the mass exchange between the two phases, $\Gamma(\phi_c, n)$. Cancerous cell death rate is taken constant, while the growth rate depends linearly on nutrient and cancer cell concentrations. Finally it reads:

$$\Gamma(\phi_c, n) = \gamma_c \phi_c (n / n_s - \delta_c) \quad (3)$$

In order to ascertain the velocity $\mathbf{v}$, mechanical ingredients are needed. A variational principle (see the Supplementary) in conjunction with the elimination of $\phi_h = 1 - \phi_c$¹⁴,¹⁵ leads to a Darcy law for the velocity of the cancerous phase $\mathbf{v}_c$.

$$\mathbf{v}_c = -K(1 - \phi_c)^2 \nabla \Sigma(\phi_c) = -K(1 - \phi_c)^2 \nabla(f(\phi_c)) - e^2 \Delta \phi_c \quad (4)$$

where $\Sigma(\phi_c)$ is the derivative of the free-energy of cancerous cell-cell interactions with respect to $\phi_c$. As the cell-cell interaction is weak at low concentration, becoming repulsive at high concentration and attractive in between, $f(\phi_c)$ in Eq. (4) has a simple form reminding the Lennard-Jones potential (see Fig. 4). A surface term $e^2 \Delta \phi_c$ is

<figure><img src="image_5.png" /><figcaption>Figure 3 | Schematic representation of the epidermis.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>Figure 4 | Σ, derivative of the free-energy as a function of concentration, $f(\phi_c)$ representing the attraction between melanocytic cells. $\phi_e$ is the concentration at equilibrium with the surrounding tissue.</figcaption></figure>

added to $\Sigma(\phi_c)$ penalizing large concentration gradients.

Gathering the former equations gives us a master equation for the cancerous phase $\phi_c$:

$$\partial_t \phi_c + \nabla \cdot (-K \phi_c (1 - \phi_c)^2 \nabla (f(\phi_c) - \epsilon^2 \Delta \phi_c) = \gamma \phi_c n - \delta \phi_c \quad (5)$$

This model takes into account most of the cellular pathways affected by melanoma: apoptosis in the term $\delta$, proliferation in $\gamma$, and cell-cell interactions in $\Sigma(\phi)$.

The governing equations for the nutrient concentration and cancerous phase, (1) and (5) can be rewritten in terms of dimensionless quantities: $\hat{x}=x/l_n$, $\hat{t}=t/\gamma n_s$, $\hat{\delta}=\delta/(\gamma n_s)$, $\hat{n}=n/n_s$, $\hat{f}=f/\chi$, $\hat{\epsilon}^2=\epsilon^2/(\chi l_n^2)$, $D=\chi/(M_n^2\gamma n_s)$, $\beta = S_n/\delta n$, $\hat{\kappa}=\kappa/\delta n$ where $l_n = \sqrt{D_n/\delta n}$ is the nutrient penetration length and $\chi$ is the typical value of the excess pressure generated by the cell (typical value of $\Sigma$). Notice that uncertainties in the precise knowledge of biomechanical factors do not change our treatment and main results. As an example, different values of the oxygen diffusion coefficient can be found in the literature$^{34,35}$ depending on the location inside the skin and on the orientation (see Table 1). These uncertainties in the oxygen diffusion coefficient (a factor two is commonly found) enter only in the pattern length-scale (by a factor 1.4) without changing the results drastically. Dropping the hats and the index for the cancerous phase for simplification, the dimensionless equations are restricted to:

$$\Delta n - \kappa(1 - \phi)n - \phi n + \beta(1 - n) = 0 \quad (6)$$

$$\frac{\partial \phi}{\partial t} = D \nabla \cdot (\phi(1-\phi)^2 \nabla (f(\phi) - \epsilon^2 \Delta \phi)) + \phi(n-\delta) \quad (7)$$

**Analytical predictions.** Although simplified in comparison with avascular tumour growth modeling$^5$, these two equations remain difficult to solve explicitly and numerics is required. It is why we begin by an analysis of possible explicit cases.

One can grasp the main results of Eq. (6) and (7) by focusing on the simplest solutions compatible with these equations. Then, we study the dynamical stability of such solutions. We consider first the time evolution of the tumour concentration.

This set of Eq.(6,7) accepts obvious time and space independent solutions both for $\phi$ and $n$ given by

$$\begin{cases} \phi_{c.s.h.} = \frac{\beta(1-\delta) - \kappa\delta}{\delta(1-\kappa)} \\ n_{c.s.h.} = \delta \end{cases} \quad \text{for the cancerous phase,} \quad (8)$$

$$\begin{cases} \phi_0 = 0 \\ n_0 = \frac{\beta}{\beta + \kappa} \end{cases} \quad \text{for the healthy phase.} \quad (9)$$

This set exists if and only if $\beta$ does not vanish. Let us remind that $\beta$ measures the strength of the vertical flux of nutrients, being a specificity of the melanoma geometry in the epidermis. It does not exist for spheroid as example. Let us consider now the possibility of small micro-structures and averaged quantities on a scale larger than their size but smaller than the tumour size. In this hypothesis, neglecting the mass creation at the interface and all border effects, integrating equations (6) and (7) on a surface immersed in the lesion, and taking nutrients at equilibrium give us:

$$\hat{\partial}_t \langle \phi \rangle_x + \delta \langle \phi \rangle_x = \frac{\beta(1 - \langle n \rangle_x) - \kappa \langle n \rangle_x}{1 - \kappa}$$

This equation can be solved by an exponential law, whose decreasing time constant is $\tau = 1/\delta$. So, inside a wide enough tumour, after a typical time about 5 $\tau$, the average concentrations $\langle \phi \rangle_x$ and $\langle n \rangle_x$ should be given by (8) and should be constant over time. As expected the nutrient diffusion from the dermis insures the existence of such homogeneous state for the cancer cell population which is in

<table><caption>Table 1 | Estimation of the model parameters from experimental data on healthy and diseased skin</caption><tr><td>φ<sub>e</sub> cell volume fraction at mechanical equilibrium</td><td>0.6–0.9</td><td>28</td></tr><tr><td>M interphase friction</td><td>963–11571 mm<sup>-2</sup> Pa day</td><td>29, 30</td></tr><tr><td>χ interstitial fluid pressure in healthy skin</td><td>133 Pa</td><td>28</td></tr><tr><td>χ interstitial fluid pressure in skin carcinoma</td><td>1330 Pa</td><td>28</td></tr><tr><td>(ε/√χ) melanoma cell size</td><td>6–20 μm</td><td>31</td></tr><tr><td>γ<sub>c</sub> melanoma cell proliferation rate</td><td>0.2 day<sup>-1</sup></td><td>32</td></tr><tr><td>δ<sub>c</sub> threshold for cell death rate due to anoxia</td><td>0.1–0.33</td><td>33</td></tr><tr><td>D<sub>n</sub> Oxygen lateral diffusion coefficient in the stratum corneum</td><td>39.7 mm<sup>2</sup> day<sup>-1</sup></td><td>34</td></tr><tr><td>D<sub>n</sub> Oxygen perpendicular diffusion coefficient</td><td>18.5–26.6 mm<sup>2</sup> day<sup>-1</sup></td><td>35</td></tr><tr><td>D<sub>n</sub> Sugar diffusion coefficient</td><td>4.4–6.4 mm<sup>2</sup> day<sup>-1</sup></td><td>36</td></tr><tr><td>δ<sub>n</sub> Oxygen consumption rate of the skin</td><td>1190–3310 day<sup>-1</sup></td><td>35</td></tr><tr><td>S<sub>n</sub></td><td>10<sup>4</sup> day<sup>-1</sup></td><td>14</td></tr><tr><td>n<sub>s</sub> Oxygen partial pressure in the skin</td><td>3320–10400 Pa</td><td>35</td></tr><tr><td>l<sub>n</sub> nutrient penetration length</td><td>0.04–0.18 mm</td><td></td></tr><tr><td>D</td><td>1.7–1725</td><td></td></tr><tr><td>ε</td><td>0.02–0.5</td><td></td></tr><tr><td>δ</td><td>0.1–0.33</td><td></td></tr><tr><td>β</td><td>0.1–8.4</td><td></td></tr><tr><td>κ</td><td>0.1–10</td><td></td></tr></table>
<figure><img src="image_7.png" /></figure>

contradiction with the existence of a stable necrotic core ($\Gamma = 0$ and $\phi_c \neq 0$ inside the tumor). However, our average does not consider the border, the tumour grows through the border, in analogy with epithelial cell colony experiments where mitosis is localised at the border¹⁶.

We show that the large scale concentration inside a lesion is given by (8), but the structure at small scale could be more complicated with emergence of microstructures.

**Space distribution analysis.** An exact solution of the concentration inside the tumour cannot be found analytically, a numerical simulation is needed. Nevertheless, it is possible to predict mathematically the existence of a phase separation. An homogeneous distribution inside the tumour with $\langle\phi\rangle_x = \phi_{c,s,h}$ exists if this solution is stable when perturbations occur. It is why we perform a linear stability analysis in the vicinity of the stationary homogeneous solution. An infinitesimal perturbation of wavevector $k$ can be written as:

$$ \phi = \phi_{c,s,h} + \tilde{\phi} \exp(\lambda t) \cos(kx) \quad (10) $$

$$ n = n_{c,s,h} + \tilde{n} \exp(\lambda t) \cos(kx) \quad (11) $$

The stability of the uniform concentration inside the tumour (8) is insured by a negative growth rate $\lambda < 0$ for all values of the spatial wavevector $k$. Therefore, substituting these expressions in equations (6) and (7), we obtain the growth rate of the perturbation:

$$ \lambda = -\phi_{c,s,h} \left[ D(1 - \phi_{c,s,h})^2 (f'(\phi_{c,s,h})k^2 + \epsilon^2 k^4) + \frac{\delta(1-\kappa)}{k^2 + \phi_{c,s,h} + \beta + \kappa(1-\phi_{c,s,h})} \right] \quad (12) $$

The diagram of $\lambda(k)$ displayed in Fig. 5 shows that the nutrient contribution in the last term of Eq. (12) stabilizes the long wavelength perturbations with small $k$. This insures that the average concentration inside the tumour $\langle\phi\rangle_x$ is given by (8)¹⁷. Because of this stabilization, an homogeneous distribution of the concentration inside the tumour may be stable at all wavelengths (dashed curve of Fig. 5) and microstructures may not automatically appear. The domain size of microstructures is determined by the nutrient penetration length, showing an effect of the micro-environment of the lesion.

## Discussion

Our simulations (see Fig. 6 and Materials and Methods) exhibit several behaviours that are also seen in vivo. Microstructural patterns appear in the second and third simulations. Their origin is a phase

<figure><img src="image_8.png" /><figcaption>**Figure 5 | Growth rate $\lambda$ of a spatial instability inside the tumour versus the wavenumber $k$ of a spatial perturbation.** Long-wavelengths, small $k$ are stabilized due to nutrients. The homogeneous distribution of the mixture is linearly unstable for $f'(\phi_0) = 0.1$ (continuous line), while it becomes stable for $f'(\phi_0) = -0.03$ (dashed line). Other material parameters are set at $\epsilon = 0.08$, $D = 2.0$, $\beta = 0.2$, $\phi_0 = 0.3$, $\kappa = 0$. The parameters taken for the continuous line lead to a phase separation then predicting micro-structures, whereas for the dashed line the concentration should remain uniform.</figcaption></figure>

separation process. It corresponds to the dots observed in dermoscopy and to the nests seen in histopathology (see Fig. 7b and Fig. 7c). Since the concentration inside microstructures is close to the mechanical equilibrium concentration $\phi_e$, the simulations show that with this separation process cells insure their mechanical equilibrium at small scales. In the same time, they insure their equilibrium at large scales, because the long-scale concentration is given by the interaction with nutrients and the difference between proliferation and death as predicted previously (see Analytical predictions).

Moreover, this phase separation leads to a growth asymmetry, whereas in the first picture the radial symmetry is preserved. Therefore, this model explains for the clinical differences between nevi and melanomas. In addition and especially a necrotic core is not found. But, the asymmetry and irregularity displayed in the simulation can be attributed to either a melanoma or a dysplastic nevus that actually show the same aspect at the tissue level. While the micro-environment can be responsible for the different characteristics between a nevus, a dysplastic nevus or a melanoma, the difference between dysplastic nevus and melanoma is the existence of growth arrest for the nevus which stops its final size.

For a range of biological parameters, this model leads to $\phi_{c,s,h} > 1$ which seems to be inconsistent. The choice of $\Sigma$ with a divergence as $\phi \to 1$, should prevent the concentration to exceed one. Nevertheless, because of the geometry chosen, restriction to a 2D model for the epidermis and a continuous flux of incoming nutrient, the concentration may diverge in the simulation. The growth of a circular tumour of size $R$ is proportional to the size of the border:

$$ \Phi_{growth} \approx 2\Pi R \nabla \Sigma \leq \frac{2\Pi R f'_{\max}}{\epsilon} $$

whereas the development inside the tumour is proportional to the size of the tumour itself so to its surface:

$$ \Phi_{reproduction} \approx \Pi R^2 \Gamma. $$

As $\Phi_{growth} < \Phi_{reproduction}$, above a critical radius $R_c$, the geometry of the avascular growth can lead to a divergence of the cell concentration. This seems like a failure of the model, with no biological consistency. But, the system can either change its shape in order to increase the size border compared to the volume creating a shape instability¹⁸,¹⁹, or grow in depth. A type of melanoma called *nodular melanoma* grows vertically from the first stages of its development. It invades the dermis quickly and keeps a round and symmetric border, and usually exhibits a exophytic nodule (small bump)²⁰ (see Fig. 7d). These findings are shown in the simulation of Fig. 6, where the cell concentration is really high and diverges. Taking into account the cellular repulsion that could break the basement membrane, it is clear that the model cannot be restricted to the epidermis and that a 3D model should be developed to investigate the vertical growth phase.

Melanogenesis cannot be restricted to a small number of mutations. Moreover, the nature and number of these mutations differ from a patient to another. Nevertheless, these mutations affect mainly 4 pathways: proliferative, anti-apoptotic, senescence, and cell interactions pathways. For instance, among the most common mutations (80–90% of cases), the p16 mutation leads to evade senescence¹³,²¹, an irreversible arrest of proliferation. The most common trigger for senescence is the telomere shortening, limiting the number of population doubling (pd) from 2 to 40 pd for adult melanocytes²¹. This number can differ for dysplastic nevi if the senescence pathways are affected, how they are affected, or if their surrounding medium changes²². All cellular activities stop after the lesion reaches a certain size when senescence is triggered. Another switch is the shift from E-cadherin (a protein which links together the melanocytes to keratinocytes) to N-cadherin that insures cell-cell adhesion between melanoma cells with the production of melanoma cells adhesion molecules called Mel-CAM²³,²⁴. This shift results in a radical change
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 6** | Several behaviours shown in the simulations on a 512 × 512 grid. The stepsize *dx* is 0.12, the steptime is 10⁻⁴. The same expression for Σ is chosen for all pictures, *p* = 2 and *φ*<sub>e</sub> = 0.6. Different values of *φ*<sub>c,s,h</sub> are taken 0.65 for (a), 0.33 for (b), 0.25 for (c) and 1.3 for (d).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 7** | a) Picture of a nevus presenting regular border. b) Picture of a dysplastic nevus presenting irregular border and microstructures.
c) Histopathological cut of a melanoma presenting melanoma nests.
d) Picture of a nodular melanoma. Figures a), b) and d) were kindly given by Dr. Pascale Guitera. Figure c) was taken from the work of Viros²⁵.</figcaption></figure>

in cell interaction, favoring melanoma cell-melanoma cell contacts to melanoma cell-keratinocyte.

An important question remains unsolved by our model: the growth arrest which occurs for dysplastic nevi for example. This growth can be limited by senescence. In the model, senescence is responsible of the maximum tumour size corresponding to a thresh- old in population doublings (pd), after which there is no more exchange of mass between phases, Γ(φ<sub>c</sub>, *n*) = 0. The number of pd is not rigorously determined, varying from 2 up to 40<sup>13</sup>. Moreover, it surely varies with the micro-environment and it differs in a tumour<sup>2</sup> from the one estimated *in vitro*. In the work of Zapperi<sup>8</sup> a pd number is estimated to be 39 which is probably over-estimated. So consider- ing a maximum pd of 32, the number of cells at the limit of sen- escence is then of order 10<sup>9</sup>. The cell size of a dysplastic melanocytic nevus has been measured with precision in the work of Rhodes<sup>26,27</sup>: the height is *h*<sub>c</sub> ~ 6 µm while the average diameter is about 11 µm (a value bigger than for ordinary melanocyte cells) and we arrive to a volume per cell of *V* ~ 570(µm<sup>3</sup>). Taking an epithelium of order 1 mm thick we arrive to a maximum size for radius of order 1.6 cm if the nevus is entirely located above the basal membrane and if only one primitive melanoma cell is present in the lesion permanently. Otherwise the maximal nevus size varies as the square root of the number of primitive melanoma cells. In any case, sen- escence is consistent with the limiting size (*R*<sub>m</sub> ~ 2 cm) of a dys- plastic nevus if we exclude giant congenital nevi.
<figure><img src="image_12.png" /></figure>

We aim to construct a simple and unique model based on physical equations that could exhibit the different features specific to melanomas and nevi. Taking into account the changes at the cellular level due to mutations and the importance of the micro-environment, this model manages to display some of the observations usually used to distinguish nevi from melanomas. It also explains the absence of a necrotic core in melanoma, even in the avascular solid tumour phase. Melanoma differs from other solid tumors by many aspects during its evolution and metastasis but the absence of necrotic core is really its specificity. In this work this is explained by the growth confined in the epidermis which modifies the usual spheroid geometry and the nutrient distribution as well.

## Methods

**Derivation of the main equations.** Eq.(6) is a modified diffusion equation in the static limit with a consumption term proportional to the cancerous cell concentration $\phi$ (according to the Michaelis-Merten law) and a source flux of nutrients. It results from an average of the 3-dimensional diffusion equation in the vertical direction, which is justified by the tiny thickness of the epidermis. It transforms the 3-dimensional diffusion equation into a simplified effective 2D diffusion equation (more details can be found in the Supplementary information). The velocity field of the cancerous cells Eq.(4) is derived from a variational principle of energy minimization including the cell-cell adhesion and the friction between phases. Elimination of the liquid phase is made possible due to conservation relations such as $\phi_c + \phi_i = 1$ (see the Supplementary) and the theoretical treatment requires only two coupled partial differential equations depending on (x,y,t): Eq.(1) and Eq.(5), instead of the 3 expected. However remaining enough complicated, inspection of these equations is done first analytically and reveals the existence of a simple homogeneous solution which allows a linear perturbation analysis as commonly done in nonlinear physics. Depending on the range of parameters entering the model and according to the dispersion relation given the growth rate of a perturbation Eq.(12) with a wavenumber $k$ (plotted in Fig. 6), one can expect or not a phase segregation inside the growing tumour.

**Simulations.** Numerical simulations of Eqs.(6,7) were performed in a two-dimensional $N \times N$ lattice, with $\phi = 0$ and $n = n_0$ at the lattice border. The initial conditions are taken as parabolic form for the concentration distribution and an initial white noise is added inside the parabola. The equations are discretized using centered second-order accurate differencing schemes in space, explicit first-order accurate in time. We use a FCTS (forward time centered space) scheme to represent Eq.(7). In order to express Eq.(6), we introduce a time dependence with a different time scale to insure that growth is slower than diffusion:

$$ \frac{n_{i,j}^{k+1} - n_{i,j}^k}{\Delta t} = \frac{n_{i+1,j}^k + n_{i-1,j}^k + n_{i,j+1}^k + n_{i,j-1}^k - 4n_{i,j}^k}{\Delta x} - \kappa(1 - \phi_{i,j}^k)n_{i,j}^k - n_{i,j}^k\phi_{i,j}^k + \beta(1 - n_{i,j}^k) $$

where $k$ is the time index, $i$ and $j$ are the grid position, and $\Delta t = 5\Delta t$. The analytical form of $f$ in $\Sigma$, (Eq.(4) in the simulations) is taken as:

$$ f(\phi) = \frac{\phi^p(\phi - \phi_c)}{1 - \phi} $$

where $p>1$, to satisfy the shape of Fig. 4. The results are not strongly dependent on the form taken for $f$. Here, we take $p=2$ and $\phi_c = 0.6$, the other parameters being chosen in the range of skin biological data given in Table 1. The program was implemented on a graphic processor unit, a GTX-580 Nvidia graphic card, using the CUDA parallel programming language, in order to reduce the computational time.

**Biological data.** An estimation of these parameters has been gathered in Table 1. We also gave an estimation of the dimensionless parameters. These quantities are used in our simulations.

1. Stern, R. S. Prevalence of a history of skin cancer in 2007: results of an incidence based model. *Arch. Dermatol* **146**, 279–282 (2010).

2. Pimiento, J. M. *et al*. Melanoma genotypes and phenotypes get personal. *Nat. Commun.* **13**, 858–867 (2013) doi:10.1038/labinvest2013.84.

3. Weinstock, M. A. Cutaneous melanoma: public health approach to early detection. *Dermatol. Ther.* **19**, 26–31 (2006).

4. Greenspan, H. P. On the growth and stability of cell cultures and solid tumors. *J. Theor. Biol.* **56**, 229–242 (1976).

5. Lowengrub, J. S. *et al*. Nonlinear modelling of cancer: Bridging the gap between cells and tumors. *Nonlinearity* **23**, 1–91 (2010).

6. Preziosi, L. & Vitale, G. A Multiphase model of tumor and tissue growth including cell adhesion and plastic reorganization. *Math.Mod. Meth. App.S.* **21**, 1901–1932 (2011).

7. Hoehme, S. & Drasdo, D. A cell-based simulation software for multi-cellular systems. *Bioinformatics* **26**, 2641–2642 (2010).

8. Zapperi, S., La Porta, C. A. & Sethna, J. P. Senescent cells in growing tumors: population dynamics and cancer stem cells. *PLoS Computational Biology* **57**, 1–13 (2012).

9. Bearer, E. L. *et al*. Multiparameter computational modeling of tumor invasion. *Cancer Res* **69**, 4493–4501 (2009).

10. Fitzpatrick, T. B. & Breathnach, A. S. Das epidermale Melanin-Einheit-System. *Dermatol. Wochenschr* **147**, 481–489 (1963).

11. Balch, C. M. *et al*. Final version of 2009 ajcc melanoma staging and classification. *J. Clin. Oncol.* **27**, 6199–6206 (2009).

12. Bachmann, I. M. *et al*. Tumor necrosis is associated with increased alphavbeta3 integrin expression and poor prognosis in nodular cutaneous melanomas. *BMC Cancer* **8**, 362 (2008).

13. Bennett, D. C. Human melanocyte senescence and melanoma susceptibility genes. *Oncogene* **22**, 3063–3069 (2003).

14. Chatelain, C., Ciarletta, P. & Ben Amar, M. Morphological changes in early melanoma development: Influence of nutrients, growth inhibitors and cell-adhesion mechanisms. *J. of Theoretical Biology* **290**, 46–59 (2011).

15. Doi, M. & Onuki, A. Dynamic coupling between stress and composition in polymer solutions and blends. *J. Physique II* **2**, 1631–1656 (1992).

16. Puliafito, A. *et al*. Collective and single cell behavior in epithelial contact inhibition. *PNAS* **109**, 739–744 (2012).

17. Chatelain, C., Balois, T., Ciarletta, P. & Ben Amar, M. Emergence of microstructural patterns in skin cancer: a phase separation analysis in a binary mixture. *New Journal of Physics* **13**, 115013 (2011).

18. Ciarletta, P., Foret, L. & Ben Amar, M. The radial growth phase of malignant melanoma. *J. Roy. Soc. Interface* **8**, 345–368 (2011).

19. Ben Amar, M., Chatelain, C. & Ciarletta, P. Contour instabilities in early tumor growth models. *Phys. Rev. Lett.* **106**, 148101 (2011).

20. Nordlund, J. J. Clinical appearance of cutaneous melanoma. *The Yale Journal of biology and medicine* **48**, 403–407 (1975).

21. Bennett, D. C. How to make a melanoma: what do we know of the primary clonal events? *Pigment Cell Melanoma Res.* **21**, 27–38 (2008).

22. Hsu, M. Y., Meier, F. & Herlyn, M. Melanoma development and progression: a conspiracy between tumor and host. *Differentiation* **70**, 522–536 (2002).

23. Haas, N. K. & Herlyn, M. Normal human melanocyte homeostasis as a paradigm for understanding melanoma. *J. Investg. Derm. Symp. Proc.* **10**, 153–163 (2005).

24. Haas, N. K. *et al*. Adhesion, migration and communication in melanoma and melanocytes. *Pigment Cell Res.* **18**, 150–159 (2005).

25. Viroš, A. *et al*. Improving Melanoma Classification by integrating Genetic and Morphologic Features. *PLoS Med.* **5**, 941–952 (2008).

26. Rhodes, A. R. *et al*. Increased intraepidermal melanocyte frequency and size in dysplastic melanocytic nevi and cutaneous melanoma. A comparative quantitative study of dysplastic melanocytic nevi, superficial spreading melanoma, nevocellular nevi, and solar lentigines. *J Invest Dermatol.* **80**, 452–459 (1983).

27. Rhodes, A. R. Melanocytic precursors of cutaneous melanoma. Estimated risks and guidelines for management. *Med Clin North Am.* **70**, 3–37 (1986).

28. Jain, R. K. Transport of molecules in the tumor interstitium: a review. *Cancer Res.* **47**, 3039–3051 (1987).

29. van Kemenade, P. M., Huyghe, J. M. & Douven, L. F. A. Triphasic FE modeling of skin water barrier. *Transp. Porous Media* **50**, 93–109 (2003).

30. Swabb, E. A., Wei, J. & Gullino, P. M. Diffusion and convection in normal and neoplastic tissues. *Cancer Res.* **34**, 2814–2822 (1974).

31. Clark, W. H., From, L., Bernardino, E. A. & Mihm, M. C. The histogenesis and biological behavior of primary human malignant melanoma of the skin. *Cancer Res.* **29**, 705–726 (1969).

32. Creasey, A. A. *et al*. Biological properties of human melanoma cells in culture. *In Vitro* **5**, 342–350 (1979).

33. Bedogni, B. & Powell, M. B. Hypoxia, melanocytes and melanoma - survival and tumor development in the permissive microenvironment of the skin. *Pigment Cell Melanoma Res.* **22**, 166–174 (2009).

34. Johnson, M. E., Blankschtein, D. & Langer, R. Evaluation of solute permeation through the stratum corneum: lateral bilayer diffusion as the primary transport mechanism. *J. Pharm. Sc.* **86**, 1162–1172 (1997).

35. Stucker, M. *et al*. The cutaneous uptake of atmospheric oxygen contributes significantly to the oxygen supply of human dermis and epidermis. *J. Physiol.* **538**, 985–994 (2002).

36. Cornelissen, L. H. *et al*. Diffusion measurements in epidermal tissues with fluorescent recovery after photobleaching. *Skin Res. Technol.* **14**, 462–467 (2008).

## Acknowledgments

This work was supported in part by AAP Physique Cancer 2012.

## Author contributions

T.B. and M.B.A. have contributed equally to this work. They wrote the main manuscript text, prepared the figures. All authors reviewed the manuscript.
<figure><img src="image_13.png" /></figure>

**Additional information**

**Supplementary information** accompanies this paper at http://www.nature.com/scientificreports

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Balois, T., Amar, M.B. & Amar, M.B. Morphology of melanocytic lesions *in situ*. *Sci. Rep.* 4, 3622; DOI:10.1038/srep03622 (2014).

<figure><img src="image_14.png" /></figure>

This work is licensed under a Creative Commons Attribution-
NonCommercial-NoDerivs 3.0 Unported license. To view a copy of this license,
visit http://creativecommons.org/licenses/by-nc-nd/3.0