Skip to main content
  • Research article
  • Open access
  • Published:

Three-dimensional experiments and individual based simulations show that cell proliferation drives melanoma nest formation in human skin tissue

Abstract

Background

Melanoma can be diagnosed by identifying nests of cells on the skin surface. Understanding the processes that drive nest formation is important as these processes could be potential targets for new cancer drugs. Cell proliferation and cell migration are two potential mechanisms that could conceivably drive melanoma nest formation. However, it is unclear which one of these two putative mechanisms plays a dominant role in driving nest formation.

Results

We use a suite of three-dimensional (3D) experiments in human skin tissue and a parallel series of 3D individual-based simulations to explore whether cell migration or cell proliferation plays a dominant role in nest formation. In the experiments we measure nest formation in populations of irradiated (non-proliferative) and non-irradiated (proliferative) melanoma cells, cultured together with primary keratinocyte and fibroblast cells on a 3D experimental human skin model. Results show that nest size depends on initial cell number and is driven primarily by cell proliferation rather than cell migration.

Conclusions

Nest size depends on cell number, and is driven primarily by cell proliferation rather than cell migration. All experimental results are consistent with simulation data from a 3D individual based model (IBM) of cell migration and cell proliferation.

Background

Clusters of melanoma cells, called nests, are a prominent feature of melanoma progression [1, 2]. Identifying the presence and characteristics of melanoma nests in human skin is an important diagnostic tool [3, 4]. Furthermore, nest size is an important characteristic because larger melanoma nests are associated with more aggressive melanoma [3]. Recent 3D experimental work by Wessels et al. [5] suggests that melanoma nest formation in Matrigel is driven by cell migration. However, nest formation might be different in human skin, where melanoma cells are in contact with other cell types [1, 6]. We hypothesise that two different mechanisms could lead to nest formation: (i) cell proliferation, where clusters of melanoma cells are formed primarily through mitosis (Fig. 1a); and (ii) cell migration, where clusters of adhesive melanoma cells form primarily through melanoma cell migration (Fig. 1b). Cell migration occurs over a short time scale of hours, whereas cell proliferation takes place over a much longer time scale of days. Since our work is focused on the role of proliferation, we perform experiments over a period of four days so that we are able to observe and quantify the role of cell proliferation. This choice of experimental time scale means that our experimental observations do not resolve the details of cell migration, which would require a much finer time resolution in the experiments.

Fig. 1
figure 1

Mechanisms that drive melanoma nest formation. Schematics illustrating: (a) proliferation-driven nests; and (b) migration-driven nests. In both cases the schematic shows an initially-uniform distribution of cells that lead to the formation of a nest either by the action of proliferation (a) or migration (b)

We use a 3D human skin experimental model [7, 8] to discriminate between these two conceptual models by performing a suite of experiments in which we systematically vary the initial density of proliferative melanoma cells placed on 3D human skin. This initial series of experiments allow us to examine the role of initial cell number in driving nest formation. All experiments are repeated using non-proliferative, gamma-irradiated melanoma cells. We find that higher initial numbers of melanoma cells lead to larger nests, and that cell proliferation leads to dramatically-larger nests. All experimental outcomes are consistent with a series of 3D simulations from an IBM [9]. These results provide insight into the mechanisms driving nest formation, showing that the mechanisms in 3D human skin are different to monoculture experiments performed in Matrigel.

Results and discussion

Confirmation that irradiated melanoma cells do not proliferate and are capable of migrating in a two-dimensional barrier assay

Experiments involving populations of proliferative melanoma cells are performed using non-irradiated SK-MEL-28 cells [10]. Experiments where melanoma cell proliferation is suppressed are performed using irradiated, but otherwise identical SK-MEL-28 cells [11, 12]. The melanoma cells are gamma-irradiated to inhibit mitosis. It is possible that irradiation may have other impacts on cellular behaviour and could also influence DNA functioning [12, 13]. We perform a series of live assays to show that irradiation does not affect the adherence or morphology of melanoma cells. These live cell assays involve placing populations of irradiated and non-irradiated melanoma cells on a two-dimensional tissue culture plate and making observations of cell numbers of a period of 24 h [14]. Therefore, these assays provide quantitative information about whether the populations of melanoma cells are capable of proliferating. Results confirm that irradiated melanoma cells do not proliferate. Furthermore, these assays show that irradiation does not cause the cells to die and does not affect cell morphology [see Additional file 1].

Two-dimensional (2D) barrier assays confirm that irradiated melanoma cells survive and migrate. Populations of irradiated melanoma cells are monitored over four days to confirm that irradiation does not impede the ability of cells to migrate. We use circular barrier assays to compare the spatial expansion of irradiated and non-irradiated melanoma cell populations. The leading edge of the spreading populations is detected using ImageJ [15], which also provides an estimate of the area occupied by the spreading population of cells. Since the spreading populations of cells maintain an approximately circular shape, we convert the estimates of area into an equivalent diameter and we report data in terms of the diameter of the spreading population. Results are obtained in triplicate. Images in Fig. 2a-b show the increase in the diameter of the spreading cell populations for both irradiated and non-irradiated melanoma cells over four days. The upper row of images in Fig. 2a-b, show increased spatial expansion of the population of non-irradiated cells compared to the population of irradiated melanoma cells in the lower row. Since irradiated melanoma cells do not proliferate, we expect that the size of the expanding population of irradiated cells will be smaller than the size of the expanding population of non-irradiated cells [16]. However, the area occupied by the population of irradiated melanoma cells increases over the four-day period, and this increase is due to cell migration alone. To confirm these visual observations, we also quantify the spatial spreading of irradiated and non-irradiated melanoma cell populations.

Fig. 2
figure 2

Two-dimensional spatial expansion of irradiated and non-irradiated melanoma cell monocultures. a Experimental images show barrier assays initialised with approximately 10,000 melanoma cells. The upper row of images shows non-irradiated (proliferative) melanoma cells, and the lower row shows irradiated (non-proliferative) melanoma cells. The images show the spreading of cell populations at zero, two and four days, respectively. The scale bar is 2 mm in each image. b Experimental images from (a) analysed by ImageJ. Results show the position of the leading edge of the spreading population (red) superimposed on images of the spreading populations. The upper row of images corresponds to non-irradiated melanoma cells, and the lower row of images show irradiated melanoma cells. The images show the spreading of cell populations at zero, two and four days, respectively. The scale bar in each image is 3 mm. c Data shows the average diameter of the spreading populations as a function of time (n = 3). All data generated using non-irradiated melanoma cells is in blue, and data generated using irradiated melanoma cells is in red. Plots in (c) also show the variability. The error bars correspond to the sample standard deviation (n = 3)

Data in Fig. 2c shows the increase in diameter of both irradiated and non-irradiated melanoma cell populations over four days. At all times considered, the average diameter of the irradiated cell population is less than the average diameter of the non-irradiated cell population. This is expected because the irradiated melanoma cells do not proliferate, and it is known that proliferative populations of cells expand and invade the surrounding empty space faster than non-proliferative populations of cells [9, 16]. Most importantly, the experiments initialised with irradiated melanoma cells show an increase in the diameter of the spreading population, confirming that irradiation does not prevent migration. All experiments are performed in triplicate and the averaged results are presented. We now use both, irradiated and non-irradiated melanoma cells in 3D experiments to identify the mechanism that drives melanoma nest formation.

Identifying the dominant mechanism driving melanoma nest formation

Nests of melanoma cells are well-characterised histological features of melanoma progression. Early identification of these nests is critical for successful melanoma treatment. However, in addition to examining the presence of melanoma nests, it is important to identify the biological mechanisms that lead to nest formation as this information might be relevant to the development of new drugs. To examine these pathways we use a 3D experimental skin model.

Irradiated and non-irradiated melanoma cells are cultured with primary keratinocytes and primary fibroblasts in the 3D experimental skin model for four days. From this point we refer to keratinocyte and fibroblast cells as skin cells. All cells are initially placed onto the 3D experimental skin model as a monolayer, as uniformly as possible. MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) assays highlight the metabolic activity of all cells, and show the spatial extent and spatial structure of cells on the top surface of the 3D experimental skin model. Images in Fig. 3a-b show prominent dark purple clusters on the surface of some 3D experimental skin models. Control studies, where 3D experiments are constructed without melanoma cells, show a complete absence of nests [see Additional file 1] suggesting that the dark purple clusters in Fig. 3a-b are melanoma nests. We make the natural assumption that higher densities of metabolically active cells are associated with darker purple colouration.

Fig. 3
figure 3

Proliferation drives melanoma nest formation. a MTT assays show all metabolically active cells (light purple) on the surface of the 3D experimental skin model initialised with different numbers of proliferating melanoma cells, as indicated. b Equivalent results with irradiated melanoma cells. Melanoma nests are in dark purple (arrows). Scale bars are 1 mm. c-d Box plots showing nest area as a function of initial number of melanoma cells. Outliers are indicated by red crosses. Inset in (d) shows details in the range 0–0.04 mm2

Images in Fig. 3a show that larger nests are associated with higher initial numbers of melanoma cells. To quantify this we measure the area of individual nests using ImageJ [15], and data in Fig. 3c confirms our visual observation. Interestingly, larger initial numbers of melanoma cells lead to a smaller number of larger nests [see Additional file 2]. This is consistent with smaller sized nests coalescing into a smaller number of larger nests over time. These results suggest smaller nests might coalescence into larger nests over time. To confirm this conjecture we would need to analyse our experiments using time-lapse imaging. Since our results show that cell number plays a critical role, we now examine the role of proliferation by suppressing mitosis.

We examine the role of cell proliferation by constructing 3D experimental skin models with irradiated melanoma cells. Images in Fig. 3b show that this leads to the formation of dramatically smaller nests. To quantify our results, the area of individual nests is measured using ImageJ [15] [see Additional file 2]. Data in Fig. 3d shows a similar trend to data in Fig. 3c as the nest area increases with initial cell number. However, comparing results in Fig. 3c-d shows that proliferation plays a dominant role in nest formation. For example, experiments initialised with 8500 proliferative melanoma cells leads to a median nest area of 0.15 mm2, whereas the median nest area is just 0.027 mm2 when proliferation is suppressed. These measurements of nest area do not provide direct estimates of the number of cells present in each nest. However, it is reasonable to assume that larger nests contain more cells than smaller nests.

Our results are different to previous 3D studies that show melanoma nests are formed by cell migration [5]. We anticipate that the difference in our outcome could be due to: (i) differences between the melanoma cell lines used; (ii) the interaction of melanoma cells with the surrounding skin cells in the 3D experiments; or, (iii) differences in the material used to construct the 3D model described in [5] and the 3D model used in this study. Since our experiments are performed in 3D materials derived from human skin, and our experiments involve culturing melanoma cells together with primary human skin cells, we feel that our results are more realistic than examining nest formation in monoculture experiments in Matrigel. We now perform immunohistochemistry to confirm that irradiated melanoma cells survive in the 3D experimental human skin model over a period of four days.

Irradiated melanoma cells survive in a 3D experimental skin model

Here, we perform a series of experiments using a specific melanoma marker to provide additional evidence that nests observed on the 3D experimental human skin models are clusters of melanoma cells, and that irradiated melanoma cells survive in a 3D environment over four days. The 3D experimental skin models are constructed using both irradiated and non-irradiated melanoma cells. Vertical cross-sections through the 3D experimental skin models initialised with melanoma cells are stained using S100, which is a reliable melanoma cell marker [17]. Both irradiated and non-irradiated melanoma cells are found in the 3D experimental skin model after four days. Images in Fig. 4a-f show positive S100 staining of melanoma cells. In particular, Fig. 4b–f show positive S100 staining of irradiated melanoma cells after four days. This immunostaining confirms that irradiation does not alter the antigen properties of melanoma cells for this marker, and the irradiated melanoma cells survive in a 3D experimental skin model for four days. Our experimental results use skin cells and skin dermis from one donor. Additional results using cells and dermis from two other donors show little variability between them.

Fig. 4
figure 4

Irradiated and non-irradiated melanoma cells survive in 3D experimental skin models. S100 identifies melanoma cells (brown), and the arrows indicate positive staining. a, c and e Cross-sections through 3D experimental skin models initialised with 1250, 5000 and 8500 non-irradiated melanoma cells, as indicated. b, d, and f Cross-sections through 3D experimental skin models initialised with 1250, 5000 and 8500 irradiated melanoma cells, respectively. Scale bar in each image is 100 Î¼m

Variability between skin samples

We now examine whether there is any important variability in our results between skin samples from different donors. To examine this we perform additional experiments using dermis and primary skin cells from three different donors, which we denote as donor A, donor B and donor C. We show MTT assays on the 3D experimental skin models initialised with non-irradiated and irradiated melanoma cells in Fig. 5. The upper row of images in Fig. 5a-c show 3D experimental skin models initialised with 1250, 5000 and 8500 non-irradiated melanoma cells, respectively. In each case, we see that larger nests are associated with higher initial number of melanoma cells. A similar trend is observed for the images in the lower row of images in Fig. 5a-c where the experiments are initialised with an equivalent number of irradiated melanoma cells. However, regardless of whether we consider results from donor A, donor B or donor C, we always see that nest formation is dramatically reduced when we consider irradiated, non-proliferative melanoma cells.

Fig. 5
figure 5

Donor variability in 3D experimental skin models with melanoma cells. Experimental images show metabolically active cells (light purple) on the 3D experimental skin model after four days. The skin models are constructed using primary skin cells and dermis from three different donor samples denoted a; b; and c. The scale bars are 1 mm. The melanoma nests are shown in dark purple. In each set of subfigures, (a-c), the images in the upper row show experiments initialised with 1250, 5000 and 8500 non-irradiated melanoma cells, respectively. In each set of subfigures, (a-c), the images in the lower row show experiments initialised with 1250, 5000 and 8500 irradiated melanoma cells, respectively

Visual inspection of the images in Fig. 5 suggests that the size, shape and number of individual nests does vary slightly between the three donors. However, the influence of the initial cell number and the action of cell proliferation on nest formation remains consistent between the skin samples from the three different donors. That is, larger initial numbers of cells produces larger nests, and the action of cell proliferation leads to dramatically larger nests. To provide additional evidence we also measure the area of individual nests on skin samples from all donors using ImageJ [15]. Data provided [see Additional file 2] confirm that the relationship between initial cell number and the action of cell proliferation holds for all three donor samples.

The nests on the 3D experimental skin model initialised with 1250 irradiated melanoma cells are very small. Most experimental replicates of this particular experiment do not lead to any visually observable nests, as shown in the lower row of images in Fig. 5a. Therefore, data for nest area in these experiments is omitted [see Additional file 2]. We now use an IBM to verify our experimental outcomes.

Modelling melanoma nest formation using an individual based model

To corroborate our experimental findings, we use a random walk-based IBM to simulate the key features of the experiments. The IBM describes the spatial distribution of simulated cells on a 3D square lattice [18]. We use a 3D lattice of cross section 3 mm × 3 mm, and depth 2 mm, to represent the central region of each experimental 3D skin model (Fig. 6a). The lattice spacing is 20 Î¼m. Simulated cells are called agents. We consider non-adhesive skin agents (green, Fig. 6b) and adhesive melanoma agents (blue, Fig. 6b). Note that the domain shown in Fig. 6b is a small subregion within the overall domain so that we visualise just the upper portion of the lattice where the majority of agents are located.

Fig. 6
figure 6

IBM simulations corroborate experiments. a Experimental image showing all metabolically active cells (light purple) on a 3D experimental skin model initialised with 5000 proliferating melanoma cells. The magnified 3 mm × 3 mm region shows melanoma nests (dark purple). b Sub-region of the 3D simulated skin model with simulated skin agents (green) and simulated melanoma agents (blue). The dimension of the upper surface is 0.8 mm × 0.8 mm, and the depth is 0.4 mm. c Upper surface of the simulated skin model. d Cross-section through the simulated skin model. e-f Experimental and simulated nests initiated with varying numbers of melanoma cells, as indicated, and an equivalent density of simulated melanoma agents, respectively. All IBM results correspond to a typical realisation of the stochastic model. Results in (e) correspond to non-irradiated (proliferative) melanoma cells/agents. Results in (f) correspond to irradiated (non-proliferative) melanoma cells/agents. Images in (e-f) have dimensions 3 mm × 3 mm, and the depth is 2 mm. IBM parameters are τ = 0.01 h; ∆ = 20 Î¼m; P p (m) = 0.0004; P m (m) = 0.075; P p (s) = 0.00025; P m (s) = 0.3; and q = 0.7. These parameters correspond to cell proliferation rates of λ(m) = 0.04 /h and λ(s) = 0.025 /h, and cell diffusivities of D(m) = 500 Î¼m2/hour and D(s) = 2000 Î¼m2/hour for melanoma and skin cells, respectively [19]. Simulations with suppressed melanoma proliferation use P p (m) = 0

It is well-known that it can be difficult to quantitatively calibrate stochastic IBMs to match complicated experimental data precisely [19, 20]. Therefore we simply use parameters in the IBM that are adapted from previous work [9, 21]. These previous studies report estimates of the proliferation rate of SK-MEL-28 melanoma cells, the proliferation rate of primary human fibroblast cells, the cell diffusivity of SK-MEL-28 melanoma cells and the cell diffusivity of primary human fibroblast cells [21]. We make a reasonable assumption that the proliferation rate and cell diffusivity of keratinocyte cells are the same as the proliferation rate and the cell diffusivity of the fibroblast cells, respectively. Our estimate of the strength of cell-to-cell adhesion is also adapted from a previous study where this parameter was determined using a series of two-dimensional barrier assays with a metastatic melanoma cell line [9]. This approach of using previously-reported parameter estimates allows us to focus on understanding the roles of the key underlying biological features, such as the role of cell migration and cell proliferation, without being distracted by the secondary task of obtaining precise parameter estimates. We achieve this by using previously determined parameter estimates and simply comparing simulation results where melanoma cell proliferation is present, with simulation results where melanoma cell proliferation is suppressed.

We initialise the IBM simulations to precisely mimic the way that cells are placed onto the upper surface of the 3D skin in the experiments. To initialise the simulations we randomly place a particular number of skin and melanoma agents onto the surface of the 3D lattice. The initial number of agents in each subpopulation is chosen to match to the initial cell density in the experiments. Figures 6b-d show smaller sub-regions of the 3D simulated skin to visualise the distribution of agents on the 3D lattice as clearly as possible. Results in Fig. 6b-c show that the IBM predicts the formation of clusters of adhesive melanoma agents on the surface of the 3D lattice. Results in Fig. 6d shows how the IBM predicts the downward movement of both skin and melanoma agents. Fig. 6d shows that skin agents move deeper into the 3D lattice than the melanoma agents, while nests of melanoma agents tend to remain on the surface. Overall, the spatial arrangement of skin and melanoma agents in the IBM (Fig. 6b-d) is similar to the spatial arrangement of cells in the 3D experiments (Figs. 3 and 4) [6].

To explore the role of initial melanoma cell number in nest formation, IBM results in Fig. 6e show that nests form on the surface of the 3D lattice, and that the trends in simulated nest area are qualitatively similar to those in the corresponding experiments. Therefore, the simulation outcomes in Fig. 6e confirm that initial melanoma cell number is an important factor in driving nest formation. We also explore the role of cell proliferation by repeating the simulations in Fig. 6e without any melanoma agent proliferation. Simulation results in Fig. 6f are comparable to the corresponding experimental results, as we observe similar trends in nest size and morphology. In conclusion, similar to the experiments, our 3D simulation results indicate that melanoma nest formation is driven by initial melanoma cell number, and that the presence of melanoma proliferation leads to dramatically-larger nests.

In addition to qualitatively visualising the trends in Fig. 6, we also use the IBM results to quantitatively examine trends in simulated nest size. Boxplots in Fig. 7 show data quantifying the size of nests predicted using the IBM under four different conditions. We measure the area of individual nests in the IBM using the Image Region Analyzer in MATLAB [22]. For model realisations where nests are not clearly defined we adjust the image manually by increasing the separation between neighbouring nests so that the Image Region Analyzer accurately measures nests separately. We exclude extremely small nests that are composed of less than four agents.

Fig. 7
figure 7

Mathematical model confirms trends in nest formation data and predicts new results. Box plots of nest size, assuming the area of each melanoma agent is 20 × 20 = 400 μm2, for a typical realisation of the IBM for (a) proliferative melanoma agents; (b) non-proliferative melanoma agents; (c) proliferative melanoma agents where the proliferation rate is reduced by 50%; and, (d) proliferative melanoma agents where adhesion strength is reduced by 50%. In each case, outliers are indicated by red crosses. Results in (a-b) confirm the trends in Figs. 3c-d and 6e-f. Note the difference in the vertical scale in (c-d) compared to (a-b)

Results in Fig. 7a confirm that larger initial number of melanoma agents leads to larger simulated nests. Results in Fig. 7b show that suppressing proliferation in the IBM leads to dramatically smaller nests. These results in Fig. 7a-b correspond to the experimental results in Fig. 3c-d. Both the boxplots in Figs. 3c-d and 7a-b report nest area in the same units, therefore this is a direct comparison of the experimental observation and the prediction of the computational model.

In addition to using the IBM to replicate the experimental results, it is also straightforward to adjust the parameters in the IBM to make some simple predictions that have not been experimentally validated. Additional results in Fig. 7c show the distribution of simulated nest size when the proliferation rate of melanoma cells is reduced by half. Noting the difference in the vertical scale in Fig. 7a and b, we see that reducing the proliferation rate of melanoma cells by half leads to a reduction in simulated nest size by a factor of ten. Similarly, additional results in Fig. 7d show the distribution of simulated nest size when the strength of cell-to-cell adhesion for the melanoma cells is reduced by half. Again, noting the difference in vertical scale in Fig. 7a and d shows that reducing the strength of melanoma adhesion by half reduces the size of simulated nests be a factor of ten.

Conclusion

Our combined experimental and simulation findings demonstrate that cell proliferation plays the dominant role in melanoma nest formation. While it is well-accepted that proliferation is important in the latter stages of tumour growth [23] and in the spatial spreading of cell populations [24], our work shows that proliferation is vitally important at the very earliest stages of melanoma progression. As far as we are aware, our work if the first to use a 3D experimental human skin model incorporating irradiated and non-irradiated melanoma cells and shows that cell proliferation is the dominant mechanism that drives melanoma nest formation.

Our results, pointing to the importance of cell proliferation, are interesting for a number of reasons: (i) previous monoculture experiments report that melanoma nests are formed by cell migration in Matrigel [5]. One potential explanation for this difference is that the Matrigel experiments are very different to our experiments since we study nest formation on 3D human tissues where melanoma cells are in contact with skin cells; (ii) some previous mathematical models of cluster/nest formation focus on cell migration only, e.g. [25], whereas we find that cell proliferation plays the most important role; and (iii) our findings about the importance of cell proliferation in melanoma progression are consistent with the fact that many promising melanoma drugs aim to suppress proliferation [26,27,28].

Our suite of experimental results can be extended in many ways. For example, one limitation of our work is that we group the keratinocytes and fibroblast cells together, and refer to these cells as skin cells. It would be interesting to repeat our work and use specific markers to differentiate between these two populations of skin cells [29]. Another interesting extension of our experimental work could be to examine nest formation in 3D experiment using a mixture of irradiated and non-irradiated melanoma cells. This condition could mimic a partial reduction in proliferation, whereas our results correspond to a total inhibition of melanoma cell proliferation. Additionally in these experiments, cell proliferation can be blocked using drug treatments, such as mitomycin-C or some other commercially available proliferation inhibitor. Another relevant extension could be to perform a series of 3D skin experiments where the melanoma cells are treated so that they are non-migratory but maintain their ability to proliferate. Finally, it could also be interesting to repeat the 3D skin experiments as described here, and to image the formation of nests on a much shorter timescale that is comparable to the timescale of cell migration. If we had access to such time course data, it might then be possible to compare this kind of transient data from the experiment with transient information from a mathematical model [30].

Our suite of modelling results can also be extended in many ways. In this work we choose to work with a relatively simple mathematical model that represents just the key processes of interest, namely a population of motile, proliferative and adhesive melanoma cells, and a population of motile and proliferative skin cells. This model is parameterised using previously-determined parameter estimates [9, 21]. While our model is useful in that it can both recapitulate our experimental results as well as generating new predictions that could be verified or challenged in future experimental studies, it would also possible to repeat all our simulation results using a more complicated mathematical model. For example, other modelling approaches such as continuous-space lattice-free models [31] or discrete models with force potentials between agents [32] could also be used in this context. While it is always tempting to use a more complicated mathematical model that incorporates additional biological detail, this approach is limited in that using more complex models requires additional parameters. We avoid this situation by always working with the simplest possible mathematical model that describes just the key features of interest.

Methods

Keratinocyte isolation and culture

Queensland University of Technology (QUT) human research ethics obtained written approval for the skin samples to be used in this study (approval number: QUT HREC #1300000063; UnitingCare Health 2003/46). Skin samples are collected from patients undergoing elective plastic surgery. Human keratinocyte cells are isolated from skin and cultured in full Green’s medium containing DMEM with Ham’s F12 (Thermo Scientific, Australia) in a 3:1 v/v ratio, 10% FCS, 2 mM L-glutamine, 50 U/ml of penicillin, 50 μg/ml of streptomycin, 180 mM adenine (Sigma Aldrich, Australia), 1 μg/ml insulin, 0.1 μg/ml cholera toxin (Sigma Aldrich), 0.01% non-essential amino acid solution (Thermo Scientific), 5 μg/ml transferrin (Sigma Aldrich), 0.2 μM triiodothyronine (Sigma Aldrich), 0.4 μg/ml hydrocortisone (Sigma Aldrich) and 10 ng/ml human recombinant EGF (Thermo Sceientific) following protocols described previously [17, 33]. Primary keratinocyte cells are cultured at 37 Â°C, in 5% CO2 and 95% air.

Fibroblast isolation and culture

Human fibroblast cells are isolated following protocols in Haridas et al. [17]. Primary fibroblast cells are cultured at 37 Â°C, in 5% CO2 and 95% air.

Melanoma cell culture

The human melanoma cell line SK-MEL-28 is cultured as described in Haridas et al. [17]. SK-MEL-28 melanoma cells are kindly donated by Professor Brian Gabrielli (Mater Research Institute-University of Queensland). Cells are cultured at 37 Â°C, in 5% CO2 and 95% air.

A batch of SK-MEL-28 melanoma cells is irradiated to prevent cell proliferation. Approximately 1 × 107 melanoma cells are gamma-irradiated using a Gammacell 40 research irradiator (Australia) at approximately 0.8 Gy/min for one hour resulting in a cumulative dose of 50 Gy. We refer to these non-proliferative cells as irradiated melanoma cells, and the proliferative cells as non-irradiated melanoma cells.

Identification of SK-MEL-28 cells is validated using short tandem repeat profiling (Cell Bank, Australia. January 2015).

Barrier assay

We perform circular barrier assays to observe and measure the spreading of populations of irradiated and non-irradiated melanoma cells. The protocol from Simpson et al. [16] is followed. Briefly, sterile stainless steel silicon barriers (Aix Scientific, Germany) are carefully placed in a 24-well tissue culture plate with 0.5 ml growth medium. The tissue culture plate containing cells is incubated for one hour at 37 Â°C, in 5% CO2 and 95% air. Viable cell suspensions of 20,000 cells/100 Î¼l of irradiated and non-irradiated melanoma cells are carefully introduced into the barriers to ensure an even distribution of cells. The tissue culture plates containing cell suspensions are incubated for a further two hours to allow cells to attach to the plate. The barriers are removed and the cell layers are washed with serum-free medium (culture medium without foetal calf serum) and replaced with fresh growth medium. Plates are then incubated at 37 Â°C, in 5% CO2 and 95% air for zero, two and four days. We replace the growth medium after two days to replenish the nutrients. Each assay is performed in triplicate.

Crystal violet staining

We use the staining technique described by Simpson et al. [16] to analyse the barrier assays. In brief, cell monolayers are washed with phosphate buffered saline (PBS; Thermo Scientific, Australia) and fixed using 10% neutral buffered saline (United Biosciences, Australia) for 20 min at room temperature. The fixed cells are stained using 0.01% v/v crystal violet (Sigma Aldrich, Australia) in PBS for 20 min at room temperature. Excess crystal violet stain is removed using PBS, and the plates are air-dried. Images of irradiated and non-irradiated cell populations are acquired using a Nikon SMZ 800 stereo microscope fitted with a Nikon digital camera.

Establishing 3D experimental skin model with melanoma cells

We establish 3D experimental skin models using the skin collected from donors undergoing elective plastic surgery. The protocol for establishing the 3D skin equivalent model with melanoma cells is adapted from previous work [7]. In brief, sterile stainless steel rings (Aix Scientifics) with a radius of 3 mm are placed on the papillary side of the de-epidermised dermis in a 24-well tissue culture plate (Nunc®, Australia). We refer to the de-epidermised dermis as dermis. Single cell suspensions of primary keratinocyte cells (20000), primary fibroblast cells (10000) and non-irradiated melanoma cells (1250; 5000; 8500), are seeded onto the dermis in full Green’s medium as uniformly as possible, and incubated at 37 Â°C, in 5% CO2 and 95% air for two days. We refer to the primary keratinocyte and fibroblast cells as skin cells. Subsequently, the stainless steel rings are removed and the dermis containing cells is submerged in full Green’s medium for a further two days. After this four-day pre-culture period, the spatial distribution of cells in the 3D experimental skin model is analysed. We also perform a series of equivalent experiments using irradiated melanoma cells.

All experiments are performed in triplicate. Furthermore, all experiments are repeated using primary skin cells and dermis from three separate donors to account for variability between different donors.

MTT assay

An MTT (Thermo Scientific) assay is performed to check the metabolic activity of cells on the 3D experimental skin models. These assays are imaged with a stereo microscope (Nikon SMZ 800) fitted with a Nikon digital camera. We follow the protocol from Haridas et al. [7].

Immunohistochemistry on 3D experimental skin models with melanoma cells

We use immunohistochemistry to identify melanoma cells in the 3D experimental skin models. 10% neutral buffered formalin (United Biosciences, Australia) is used to fix the 3D experimental skin models. The tissue is divided through the centre of the MTT positive region using a sterile blade. The two smaller pieces of tissue are processed and embedded in paraffin. These samples are sectioned into 5 Î¼m thick sections using a microtome. These sections are de-paraffinised, rehydrated and then subjected to heat-mediated antigen retrieval treatment using sodium citrate buffer (pH 6.0) in a decloaking chamber (Biocare Medical, USA) at 95 Â°C for 5 min. Skin sections are washed in PBS followed by immunostaining using the MACH 4â„¢ Universal HRP polymer kit (Biocare Medical). The primary antibody S100 (Dako, Australia) is diluted in DaVinci Green diluent (Biocare Medical) at 1:3000, and these sections are incubated with the primary antibody for one hour at room temperature. Positive immunoreactivity is visualized using 3,3-diaminobenzidine (DAB; Biocare Medical) and then counterstained with using Gill’s haematoxylin (HD Scientific, Australia). The sections are dehydrated, and mounted on coverslips using Pertex® mounting medium (Medite, Germany). All stained sections are imaged using an Olympus BX41 microscope fitted with an Olympus digital camera (Micropublisher, 3.3RTV, QImaging; Olympus, Q-Imaging, Tokyo, Japan).

IBM simulation methods

We use a 3D lattice-based IBM, with adhesion between some agents, to describe the 3D experiments. In the IBM, cells are treated as equally sized spheres, and referred to as agents. These agents are restricted to reside on a 3D square lattice, with no more than one agent per site. The lattice spacing, Δ, represents the approximate size of each simulated agent, or the minimum spacing between agents. Here, we set Δ = 20 Î¼m to match previous measurements [21]. We use a 3D lattice of dimension, 3 mm × 3 mm, and depth 2 mm, to represent the central region of each experimental skin model. This means that the number of lattice sites is 150 × 150 × 100. We choose the depth of the domain so that agents in the simulation never touch the bottom of the domain during the simulations. The parameters in the simulation model are adapted from previous studies [21]. Since we use the 3D lattice to represent the central region of the tissue, where cells are initialised uniformly across the surface, we apply periodic boundary conditions along all vertical boundaries. Since cells cannot leave the skin through the upper or lower surfaces, we apply no flux conditions on the upper and lower horizontal boundaries of the 3D lattice. We choose the depth of the 3D lattice to be large enough so that the agents never touch the bottom boundary of the lattice on the time scale of the simulations we consider.

To initialise simulations, we randomly place a particular number of simulated skin agents, N(s)(0), and a particular number of simulated melanoma agents, N(m)(0), onto the surface of the lattice. When the IBM is initialised we take care to ensure that no more than one agent occupies each lattice site. We always choose the initial number of agents in each subpopulation to match the equivalent initial density of cells in the experimental skin model. In the experiments, the initial populations of cells are uniformly placed inside a disc of radius 3 mm, whereas in the IBM the initial populations of agents are uniformly placed inside a square subregion of side length 3 mm. We set the initial number of skin agents to be N(s)(0) = 9549 to match the initial experimental population of 30,000 skin cells distributed in a disc of radius 3 mm. We vary the initial number of simulated melanoma agents to be N(m)(0) = 398, 1592 or 2706, to match the initial density of melanoma cells. This initial experimental density corresponds to 1250, 5000 and 8500 melanoma cells distributed in a disc of radius 3 mm. To match the experiments, the IBM simulations are run for four days.

At any time, t, there are N(t) = N(m)(t) + N(s)(t) agents on the lattice. In each discrete time step, of duration τ, we use a random sequential update method [34] to simulate motility and proliferation events. The algorithm involves executing two sequential steps:

  1. 1.

    N(t) agents are selected one at a time, with replacement and given the opportunity to move to a nearest neighbour lattice site with probability P m (s) ∈[0,1] and P m (m) ∈[0,1]. Here we can specify different motility probabilities for the skin cells and the melanoma cells, and this is important because previous work has shown that fibroblast cells are more motile than melanoma cells [21]. If the chosen agent is a melanoma agent, we incorporate adhesion into the model by examining the occupancy of the 26 nearest lattice sites in the 3D Moore neighbourhood. We count the number of those sites occupied by melanoma agents, a [18]. Potentially motile melanoma agents then attempt to move with a modified probability, P m * = (1 − q)a, which accounts for adhesion between neighbouring melanoma agents. The parameter q controls the strength of melanoma-melanoma agent adhesion, with q = 0 corresponding to no adhesion, and increasing q leading to increased adhesion [18]. Setting q = 1 corresponds to maximal adhesion, and this would prevent any motility of melanoma agents that are in contact with other melanoma agents. We do not include any adhesion between skin agents as fibroblast cells are known to be mesenchymal and act as individuals rather than being strongly affected by adhesion [16]. If a movement event is successful, the agent attempts to move to a nearest neighbour lattice site from the six sites in the 3D von Neumann neighbourhood. To simulate crowding effects, potential motility events that would place an agent on an occupied site are aborted.

  2. 2.

    N(t) agents are selected one at a time, with replacement and given the opportunity to proliferate with probability P p (s) ∈[0,1] and P p (m) ∈[0,1]. Again, this framework allows us to specify different proliferation probabilities for the skin cells and the melanoma cells [18, 21]. If a proliferation event is successful, a daughter agent is placed at a randomly chosen nearest site from the six sites in the 3D von Neumann neighbourhood. To simulate crowding effects, we abort the proliferation event if all six nearest neighbouring sites are occupied. In all cases where a proliferation event is successful, a proliferative melanoma agent will produce a daughter melanoma agent, and a proliferative skin agent will produce a daughter skin agent.

The parameters in the IBM are Δ, τ, P m (s), P m (m), P p (s) P p (m) and q. These IBM parameters are related to the cell proliferation rates (λ(s) = P p (s)/ τ, λ(m) = P p (m)/ τ) and cell diffusivities (D(s) = P m (s) Δ 2/ (6τ), D(m) = P m (m) Δ 2/ (6τ)).

References

  1. Beaumont KA, Mohana-Kumaran N, Haass NK. Modeling melanoma in vitro and in vivo. Healthcare. 2014;2:27–46. https://doi.org/10.3390/healthcare2010027.

    Article  Google Scholar 

  2. Meier F, Nesbit M, Hsu M, Martin B, Belle PV, Elder DE, Schaumburg-Lever G, Garbe C, Walz TM, et al. Human melanoma progression in skin reconstructs: biological significance of bFGF. Am J Pathol. 2000;156:193–200. https://doi.org/10.1016/S0002-9440(10)64719-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Balu M, Kelly KM, Zachary CB, Harris RM, Krasieva TB, Konig K, Durkin AJ, Tromberg BJ. Distinguishing between benign and malignant melanocytic nevi by in vivo multiphoton microscopy. Cancer Res. 2014;74:2688–97. https://doi.org/10.1158/0008-5472.CAN-13-2582.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Urso C, Rongioletti F, Innocenzi D, Batolo D, Chimenti S, Fanti PL, Filotico R, Gianotti R, Lentini M, Tomasini C, et al. Histological features used in the diagnosis of melanoma are frequently found in benign melanocytic naevi. J Clin Pathol. 2005;58:409–12. https://doi.org/10.1136/jcp.2004.020933.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Wessels D, Lusche DF, Voss E, Kuhl S, Buchele EC, Klemme MR, Russell KB, Ambrose J, Sol BA, Bossler A, et al. Melanoma cells undergo aggressive coalescence in a 3D Matrigel model that is repressed by anti-CD44. PLoS One. 2017;12:e0173400. https://doi.org/10.1371/journal.pone.0173400.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Eves P, Layton C, Hedley S, Dawson RA, Wagner M, Morandini R, Ghanem G, Mac NS. Characterization of an in vitro model of human melanoma invasion based on reconstructed human skin. Brit J Dermatol. 2000;142:210–22. https://doi.org/10.1046/j.1365-2133.2000.03287.x.

    Article  CAS  Google Scholar 

  7. Haridas P, McGovern JA, McElwain DLS, Simpson MJ. Quantitative comparison of the spreading and invasion of radial growth phase and metastatic melanoma cells in a three-dimensional human skin equivalent model. PeerJ. 2017;5:e3754. https://doi.org/10.7717/peerj.3754.

    Article  PubMed  PubMed Central  Google Scholar 

  8. MacNeil S, Eves P, Richardson B, Molife R, Lorigan P, Wagner M, Layton C, Morandini R, Ghanem G. Oestrogenic steroids and melanoma cell interaction with adjacent skin cells influence invasion of melanoma cells in vitro. Pigm Cell Melanoma Res. 2000;13:68–72. https://doi.org/10.1034/j.1600-0749.13.s8.13.x.

    Article  Google Scholar 

  9. Treloar KK, Simpson MJ, Haridas P, Manton KJ, Leavesley DI, McElwain DLS, Baker RE. Multiple types of data are required to identify the mechanisms influencing the spatial expansion of melanoma cell colonies. BMC Syst Biol. 2013;7:137. https://doi.org/10.1186/1752-0509-7-137.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Carey TE, Takahashi T, Resnick LA, Oettgen HF, Old LJ. Cell surface antigens of human malignant melanoma: mixed hemadsorption assays for humoral immunity to cultured autologous melanoma cells. P Natl Acad Sci USA. 1976;73:3278–82. https://doi.org/10.1073/pnas.73.9.3278.

    Article  CAS  Google Scholar 

  11. Deacon DH, Hogan KT, Swanson EM, Chianese-Bullock KA, Denlinger CE, Czarkowski AR, Schrecengost RS, Patterson JW, Teague MW, Slingluff CL Jr. The use of gamma-irradiation and ultraviolet-irradiation in the preparation of human melanoma cells for use in autologous whole-cell vaccines. BMC Cancer. 2008;8:360. https://doi.org/10.1186/1471-2407-8-360.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Todorovic D, Petrovic I, Todorovic M, Cuttone G, Ristic-Fira A. Early effects of gamma rays and protons on human melanoma cell viability and morphology. J Microsc-Oxford. 2008;232:517–21. https://doi.org/10.1111/j.1365-2818.2008.02151.x.

    Article  CAS  Google Scholar 

  13. Lahtz C, Bates SE, Jiang Y, Li AX, Wu X, Hahn MA, Pfeifer GP. Gamma irradiation does not induce detectable changes in DNA methylation directly following exposure of human cells. PLoS One. 2012;7:e44858. https://doi.org/10.1371/journal.pone.0044858.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Warne DJ, Baker RE, Simpson MJ. Optimal quantification of contact inhibition in cell populations. Biophys J. 2017;113:1920–4. https://doi.org/10.1016/j.bpj.2017.09.016.

    Article  CAS  PubMed  Google Scholar 

  15. Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2017;9:671–5. https://doi.org/10.1038/nmeth.2089.

    Article  Google Scholar 

  16. Simpson MJ, Treloar KK, Binder BJ, Haridas P, Manton KJ, Leavesley DI, McElwain DLS, Baker RE. Quantifying the roles of cell motility and cell proliferation in a circular barrier assay. J R Soc Interface. 2013;10:20130007. https://doi.org/10.1098/rsif.2013.0007.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Haridas P, McGovern JA, Kashyap AS, McElwain DLS, Simpson MJ. Standard melanoma-associated markers do not identify the MM127 metastatic melanoma cell line. Sci Rep. 2016;6:24569. https://doi.org/10.1038/srep24569.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Simpson MJ, Towne C, McElwain DLS, Upton Z. Migration of breast cancer cells: understanding the roles of volume exclusion and cell-to-cell adhesion. Phys Rev E. 2010;82:041901. https://doi.org/10.1103/PhysRevE.82.041901.

    Article  Google Scholar 

  19. Read MN, Alden K, Rose LM, Timmis J. Automated multi-objective calibration of biological agent-based simulations. J R Soc Interface. 2016;13:20160543. https://doi.org/10.1098/rsif.2016.0543.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Browning AP, McCue SW, Binny RN, Plank MJ, Shah ET, Simpson MJ. Inferring parameters for a lattice-free model of cell migration and proliferation using experimental data. J Theor Biol. 2018;437:251–60. https://doi.org/10.1016/j.jtbi.2017.10.032.

    Article  PubMed  Google Scholar 

  21. Haridas P, Penington CJ, McGovern JA, McElwain DLS, Simpson MJ. Quantifying rates of cell migration and cell proliferation in co-culture barrier assay reveals how skin and melanoma cells interact during melanoma spreading and invasion. J Theor Biol. 2017;423:13–25. https://doi.org/10.1016/j.jtbi.2017.04.017.

    Article  PubMed  Google Scholar 

  22. Image Acquisition Toolbox Documentation R2017b. In: Mathworks. [https://au.mathworks.com/help/images/calculate-region-properties-using-image-region-analyzer.html]. Accessed 17 Jan 2018.

  23. Gerlee P. The model muddle: in search of tumor growth laws. Cancer Res. 2013;73:2407–11. https://doi.org/10.1158/0008-5472.CAN-12-4355.

    Article  CAS  PubMed  Google Scholar 

  24. Vo BN, Drovandi CC, Pettitt AN, Simpson MJ. Quantifying uncertainty in parameter estimates for stochastic models of collective cell spreading using approximate Bayesian computation. Math Biosci. 2015;263:133–42. https://doi.org/10.1016/j.mbs.2015.02.010.

    Article  PubMed  Google Scholar 

  25. Green JEF, Waters SL, Whiteley JP, Edelstein-Keshet L, Shakesheff KM, Byrne HM. Non-local models for the formation of hepatocyte-stellate cell aggregates. J Theor Biol. 2010;267:106–20. https://doi.org/10.1016/j.jtbi.2010.08.013.

    Article  CAS  PubMed  Google Scholar 

  26. Chan KS, Koh CG, Li HY. Mitosis-targeted anti-cancer therapies: where they stand. Cell Death Dis. 2012;3:e411. https://doi.org/10.1038/cddis.2012.148.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Ramaraj P. In vitro inhibition of human melanoma (BLM) cell growth by progesterone receptor antagonist RU-486 (Mifprestone). J Caner Ther. 2016;7:1045–58. https://doi.org/10.4236/jct.2016.713101.

    Article  Google Scholar 

  28. Lai X, Friedman A. Combination therapy for melanoma with BRAF/MEK inhibitor and immune checkpoint inhibitor: a mathematical model. BMC Syst Biol. 2017;11:70. https://doi.org/10.1186/s12918-017-0446-9.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Krtolica A, Solorzano CO, Lockett S, Campisi J. Quantification of epithelial cells in coculture with fibroblasts by fluorescence image analysis. Cytometery. 2002;49:73–82. https://doi.org/10.1002/cyto.10149.

    Article  Google Scholar 

  30. Simpson MJ, Landman KA, Hughes BD, Fernando AE. A model for mesoscale patterns in motile populations. Physica A. 2010;389:1412–24. https://doi.org/10.1016/j.physa.2009.12.010.

    Article  CAS  Google Scholar 

  31. Bruna M, Chapman SJ. Diffusion of finite-size particles in confined geometries. Bull Math Biol. 2014;76:947–82. https://doi.org/10.1007/s11538-013-9847-0.

    Article  PubMed  Google Scholar 

  32. Matsiaka OM, Penington CJ, Baker RE, Simpson MJ. Continuum approximations for lattice-free multi-species models of collective cell migration. J Theor Biol. 2017;422:1–11. https://doi.org/10.1016/j.jtbi.2017.04.009.

    Article  CAS  PubMed  Google Scholar 

  33. Dawson RA, Upton Z, Malda J, Harkin DG. Preparation of cultured skin for transplantation using insulin-like growth factor I in conjunction with insulin-like growth factor binding protein 5, epidermal growth factor, and vitronectin. Transplantation. 2006;81:1668–76. https://doi.org/10.1097/01.tp.0000226060.51572.89.

    Article  CAS  PubMed  Google Scholar 

  34. Chowdhury D, Schadschneider A, Nishinari K, et al. Phys Life Rev. 2005;2:318–52. https://doi.org/10.1016/j.plrev.2005.09.001.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Brian Gabrielli for the SK-MEL-28 cell line. We thank four anonymous referees for their helpful suggestions.

Funding

We are supported by the Australian Research Council (DP170100474).

Availability of data and materials

All data generated and/or analysed supporting the outcome of this manuscript are included within this manuscript, in Additional files 1 and 2.

Author information

Authors and Affiliations

Authors

Contributions

PH and MJS conceived the study. PH, APB, JAM, DLSM and MJS designed the experiments. PH performed the experiments. APB performed the mathematical simulations. PH and APB analysed the data. PH and MJS wrote the manuscript. PH, APB, JAM, DLSM and MJS read, edited and approved the final manuscript.

Corresponding author

Correspondence to Matthew J. Simpson.

Ethics declarations

Ethics approval and consent to participate

Ethics approval and informed written consent from patients has been granted to conduct this study. Ethics approval number: QUT HREC #1300000063; UnitingCare Health 2003/46.

Consent for publication

Written consent from patients to publish data generated from this project has been obtained. All patient data generated in this study are unidentifiable, and there are no details about individuals reported within this manuscript.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Additional experimental results. (DOCX 2144 kb)

Additional file 2:

Experimental data (XLSX 22 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Haridas, P., Browning, A.P., McGovern, J.A. et al. Three-dimensional experiments and individual based simulations show that cell proliferation drives melanoma nest formation in human skin tissue. BMC Syst Biol 12, 34 (2018). https://doi.org/10.1186/s12918-018-0559-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12918-018-0559-9

Keywords