Clinical diagnosis of skin cancers is based on several morphological criteria, among which is the presence of microstructures (e.g. dots and nests) sparsely distributed within the tumour lesion. In this study, we demonstrate that these patterns might originate from a phase separation process. In the absence of cellular proliferation, in fact, a binary mixture model, which is used to represent the mechanical behaviour of skin cancers, contains a cell–cell adhesion parameter that leads to a governing equation of the Cahn–Hilliard type. Taking into account a reaction–diffusion coupling between nutrient consumption and cellular proliferation, we show, with both analytical and numerical investigations, that two-phase models may undergo a spinodal decomposition even when considering mass exchanges between the phases. The cell–nutrient interaction defines a typical diffusive length in the problem, which is found to control the saturation of a growing separated domain, thus stabilizing the microstructural pattern. The distribution and evolution of such emerging cluster morphologies, as predicted by our model, are successfully compared to the clinical observation of microstructural patterns in tumour lesions.
Read full abstract