A model for selection of eyespots on butterfly wings

Sekimura, Toshio, Venkataraman, Chandrasekhar and Madzvamuse, Anotida (2015) A model for selection of eyespots on butterfly wings. PLoS ONE, 10 (11). pp. 1-23. ISSN 1932-6203

[img] PDF (Article) - Published Version
Available under License Creative Commons Attribution No Derivatives.

Download (4MB)


The development of eyespots on the wing surface of butterflies of the family Nympalidae is one of the most studied examples of biological pattern formation.However, little is known about the mechanism that determines the number and precise locations of eyespots on the wing. Eyespots develop around signaling centers, called foci, that are located equidistant from wing veins along the midline of a wing cell (an area bounded by veins). A fundamental question that remains unsolved is, why a certain wing cell develops an eyespot, while other wing cells do not. We illustrate that the key to understanding focus point selection may be in the venation system of the wing disc. Our main hypothesis is that changes in morphogen concentration along the proximal boundary veins of wing cells govern focus point selection. Based on previous studies, we focus on a spatially two-dimensional reaction-diffusion system model posed in the interior of each wing cell that describes the formation of focus points. Using finite element based numerical simulations, we demonstrate that variation in the proximal boundary condition is sufficient to robustly select whether an eyespot focus point forms in otherwise identical wing cells. We also illustrate that this behavior is robust to small perturbations in the parameters and geometry and moderate levels of noise. Hence, we suggest that an anterior-posterior pattern of morphogen concentration along the proximal vein may be the main determinant of the distribution of focus points on the wing surface. In order to complete our model, we propose a two stage reaction-diffusion system model, in which an one-dimensional surface reaction-diffusion system, posed on the proximal vein, generates the morphogen concentrations that act as non-homogeneous Dirichlet (i.e., fixed) boundary conditions for the two-dimensional reaction-diffusion model posed in the wing cells. The two-stage model appears capable of generating focus point distributions observed in nature.

We therefore conclude that changes in the proximal boundary conditions are sufficient to explain the empirically observed distribution of eyespot focus points on the entire wing surface. The model predicts, subject to experimental verification, that the source strength of the activator at the proximal boundary should be lower in wing cells in which focus points form than in those that lack focus points. The model suggests that the number and locations of eyespot foci on the wing disc could be largely controlled by two kinds of gradients along two different directions, that is, the first one is the gradient in spatially varying parameters such as the reaction rate along the anterior-posterior direction on the proximal boundary of the wing cells, and the second one is the gradient in source values of the activator along the veins in the proximal-distal direction of the wing cell.

Item Type: Article
Keywords: Eyespot pattern formation, bulk-surface reaction-diffusion systems, finite elements, diffusion-driven instability, butterfly wings
Schools and Departments: School of Mathematical and Physical Sciences > Mathematics
Subjects: Q Science
Q Science > QA Mathematics
Q Science > QA Mathematics > QA0297 Numerical analysis
Q Science > QA Mathematics > QA0299 Analysis. Including analytical methods connected with physical problems
Q Science > QH Natural history > QH0001 Natural history (General)
Q Science > QH Natural history > QH0001 Natural history (General) > QH0075 Nature conservation
Q Science > QH Natural history > QH0301 Biology > QH0359 Evolution
Q Science > QH Natural history > QH0301 Biology
Depositing User: Anotida Madzvamuse
Date Deposited: 09 Nov 2015 10:01
Last Modified: 03 Jul 2019 01:46
URI: http://sro.sussex.ac.uk/id/eprint/57772

View download statistics for this item

📧 Request an update
Project NameSussex Project NumberFunderFunder Ref
Mathematical Modelling and Analysis of Spatial Patterning on Evolving SurfacesG0872EPSRC-ENGINEERING & PHYSICAL SCIENCES RESEARCH COUNCILEP/J016780/1
InCeM: Research Training Network on Integrated Component Cycling in Epithelial Cell MotilityG1546European Union642866
Unravelling new mathematics for 3D cell migrationG1438LEVERHULME TRUSTRPG-2014-149