Recent tumor growth models are often based on the multiphase mixture framework. Using bifurcation theory techniques, we show that such models can give contour instabilities. Restricting to a simplified but realistic version of such models, with an elastic cell-to-cell interaction and a growth rate dependent on diffusing nutrients, we prove that the tumor cell concentration at the border acts as a control parameter inducing a bifurcation with loss of the circular symmetry. We show that the finite wavelength at threshold has the size of the proliferating peritumoral zone. We apply our predictions to melanoma growth since contour instabilities are crucial for early diagnosis. Given the generality of the equations, other relevant applications can be envisaged for solving problems of tissue growth and remodeling. Introduction.-Soft tissue growth models in ideal geometries have shown shape instabilities with a special focus on morphogenesis of living systems [1][2][3]. Anisotropic or inhomogeneous growth and constraints due to boundaries are responsible for these shape instabilities. Most models explore the mechanical properties with perhaps an excessive simplification of the growth process itself. On the contrary, after several decades since the pioneering work of Greenspan [4], tumor growth models have explored all the facets of possible elementary biological processes known to date [5] with different modeling approaches at different scales: continuum models at the tissue scale, discrete models at the cells scale, and eventually hybrid continuum-discrete models [6]. The recent review by [6] reports more than 550 citations but few of them explore the possibility of shape instabilities except in simplified models [4,7]. However, models based on nonlinear coupled partial-differential equations (PDEs) require numerical investigations [8,9] with lack of understanding of the parameters at the origin of shape instabilities.Nevertheless, these instabilities are often correlated with the beginning of aggressive neoplasia: they are of utmost importance for melanoma where symmetry and regularity are key criteria in clinical methods [10]. Our aim is then to identify the main biophysical interactions at the origin of these contour instabilities. It has already been shown, both in vitro [11,12] and theoretically [13], that mechanical interactions between cells and their microenvironment are crucial for tumor evolution and can determine some macroscopic tumor properties [11].We discuss here how they can control tumor shape evolution as well, dictating the condition for the development of contour instabilities on a circular growing tumor, possibly giving rise to the patterns observed in vivo. We consider a continuous multiphase model describing the tumor as a two component mixture [6]. Although simplificative, such a model already contains the main