Spatial variation in boundary conditions can govern selection and location of eyespots in butterfly wings

Despite being the subject of widespread study, many aspects of the development of eyespot patterns in butterfly wings remain poorly understood. In this work, we examine, through numerical simulations, a mathematical model for eyespot focus point formation in which a reaction-diffusion system is assumed to play the role of the patterning mechanism. In the model, changes in the boundary conditions at the veins at the proximal boundary alone are capable of determining whether or not an eyespot focus forms in a given wing cell and the eventual position of focus points within the wing cell. Furthermore, an auxiliary surface reaction-diffusion system posed along the entire proximal boundary of the wing cells is proposed as the mechanism that generates the necessary changes in the proximal boundary profiles. In order to illustrate the robustness of the model, we perform simulations on a curved wing geometry that is somewhat closer to a biological realistic domain than the rectangular wing cells previously considered, and we also illustrate the ability of the model to reproduce experimental results on artificial selection of eyespots.