A Model for Population Dynamics of the Mimetic Butterfly Papilio polytes in Sakishima Islands, Japan (II)

Based on recent progresses of both experiment and mathematical analysis, we present an extension of the model for population dynamics of the mimetic swallowtail butterfly Papilio polytes in Sakishima Islands, Japan (Sekimura et al. J Theor Biol 361:133–140, 2014). The model includes four major variables, that is, population densities of three kinds of butterflies (two female forms f. cyrus and f. polytes and the unpalatable butterfly Pachliopta aristolochiae) and their predator. In this extension, we introduce difference in the predation rate between two forms f. cyrus and f. polytes. We still assume that both the benefit of mimicry for the mimic f. polytes and the cost for the model are dependent on their relative frequencies, i.e., the mortality of the mimic by predation decreases with increase in frequency of the model, while the mortality of the model increases as the frequency of the mimic increases. Taking the density-dependent effect by carrying capacity into account, we set up an extended model system consisting of three ordinary differential equations (ODEs), analyze it mathematically, and provide computer simulations that confirm the analytical results.


Introduction
P. polytes is a mimetic swallowtail butterfly species widely distributed across India and Southeast Asia, including Southeast of China, the Philippines, Taiwan, and the Ryukyu Islands of Japan (Clarke and Sheppard 1972). P. polytes exhibits the female limited polymorphism, that is, the female is polymorphic, whereas the male is monomorphic and exhibits a white bar on the black hindwing. In the Ryukyu Islands located in the southwest of Japan, the female of P. polytes has two different forms, the mimetic form f. polytes and the nonmimetic form f. cyrus resembling the monomorphic male in appearance. The form f. polytes mimics the unpalatable butterfly Pachliopta aristolochiae as a mimetic model, which has a large white area in the center and a row of submarginal red spots on the black hindwing ( Fig. 12.1). Mimicry in the female of P. polytes is known to be Batesian mimicry. Yamauchi (1994) built a population dynamic model of Batesian mimicry, in which two populations of both model and mimic species were considered. The dynamic model has two components, growth at intrinsic growth rate and carrying capacity and reduction by predation. The probability of a predator catching prey on an encounter was assumed to depend on the frequency of the mimic. He applied the dynamic model to field records of butterflies in Ryukyu Islands, and his model has successfully explained some features, e.g., multiple dynamic equilibria between the model and the mimic in the field. However, his model did not account for realistic model-mimic systems such as polymorphism in P. polytes and intraspecific competition between the mimic f. polytes and the non-mimic f. cyrus. Sekimura et al. (2014) presented a mathematical model for population dynamics of P. polytes observed in the Sakishima Islands (i.e., the Miyako Islands and the Yaeyama Islands), which are the southernmost island group of the Ryukyu Islands, Japan ( Fig. 12.2). The model system consists of three ordinary differential equations (ODEs), in which variables are population densities of three butterflies, the unpalatable butterfly P. aristolochiae and two female forms of P. polytes. The model was constructed on the basis of field data in the islands and also experimental data in the laboratory. Using mathematical analysis and computer simulations of the system equations, they clarified the logical relationship hidden behind field data on the population dynamics of the mimetic butterfly P. polytes in Sakishima Islands. In particular, they discussed both temporal change in the relative abundance (RA) since 1975 in Miyako-jima Island and variation in the RA in Sakishima Islands.
Before going into details of our extended model and results, we summarize field data briefly in Sect. 12.2 and describe main features of mimicry of P. polytes in Sect. 12.3.  Uesugi (1992) observed temporal change of the relative abundance of the mimic f. polytes in all the females of P. polytes for 14 years from 1975 to 1989 after the establishment of the model P. aristolochiae in the land. Here, the relative abundance (RA) of the mimic denotes the population ratio of f. polytes to all the two forms (f. cyrus and f. polytes) and is defined as follows:

Field Records of
RA ¼ population of the mimic f: polytes population of the non-mimic f: cyrus þ population of the mimic f: polytes Â 100 % ð Þ The observational result revealed that the RA increases with the date like the sigmoidal curve, and in 1985 about 10 years after starting observation, the RA reached at a saturated value (or the value of equilibrium) of about 50% (Fig. 12.5a). Map of the Sakishima Islands, Japan. The Sakishima Islands are located at the southeast end of Japan and part of the Ryukyu Islands, which include both the Miyako Islands and the Yaeyama Islands. The Miyako Islands include Miyako-jima Is., Tarama-jima Is., etc., and the Yaeyama Islands include Ishigaki-jima Is., Hateruma-jima Is., Iriomote-jima Is., Taketomi-jima Is., Kohama-jima Is., etc.

Variation in the Relative Abundance
(RA) in Sakishima Islands Uesugi (1992) also observed variation in the RA in Sakishima Islands. In order to investigate the relationship between two populations of the model P. aristolochiae and the mimic f. polytes, he recorded three populations of the model and two female forms of P. polytes in seven islands from island to island in the Sakishima Islands in 1982, 14 years after the establishment of the model in the Yaeyama Islands. The horizontal axis in Fig. 12.5b shows the advantage index (AI) of Batesian mimicry, which is defined to be the population ratio of the model to all the related butterflies as follows: AI ¼ population of the model population of the model þ populations of f: cyrus and f: polytes Â 100 % ð Þ The vertical axis in Figs. 12.5a and 12.5b is the relative abundance (RA), which denotes the population rate of the form f. polytes to all two forms (f. cyrus and f. polytes). Field data, which is shown by solid circles in Fig. 12.5b, clearly shows the positive correlation between the AI and the RA, which means that the higher the population ratio of the model butterfly in an island, the higher the ratio of the mimetic female to all females in the island.

Extended Mathematical Model for Population
Dynamics of P. polytes We first summarize fundamental facts on the mimicry, which will allow us to design the content of the mathematical model of three variables, that is, three kinds of populations of the model P. aristolochiae, the mimic f. polytes, and the non-mimic f. cyrus. The butterfly P. polytes has two female forms: one is the non-mimic f. cyrus that resembles the monomorphic male with a white bar on the black hindwing, while the other is the mimic f. polytes that resembles the unpalatable butterfly species P. aristolochiae ( Fig.12.1). The mimetic form f. polytes is considered to have an advantage over f. cyrus with respect to protection from predation, when it lives in sympatry with the model P. aristolochiae. In reality, Uesugi (1996) examined the idea positively by learning experiments. Unexperienced birds, brown-eared bulbuls, Hypsipetes amaurotis oryeri as predators, were first trained to take food from two feeders in a cavity and then offered P. aristolochiae in one of the feeders. After experiencing an uncomfortable encounter with this butterfly, the birds reduced the frequency of taking regular food from the feeder where the butterfly had been placed. On the other hand, Ohsaki (1995) paid attention to the rates of beak marks by predators on wings of both palatable and unpalatable butterflies. By analyzing the number of beak marks on butterfly wings caught in Borneo, he found that by comparing the species P. polytes and the model species P. aristolochiae, nonmimetic females were selectively attacked, while males, mimetic females, and models were attacked less. Thus, the mimic f. polytes is considered to gain benefit of reduced predation risk by living in sympatry with the model P. aristolochiae. However, we should note an experimental result on survival rates of the mimetic form f. polytes, the nonmimetic form f. cyrus, and male in Taketomi-jima Is., where the model P. aristolochiae and the mimetic female f. polytes are both absent (Uesugi 1997). Uesugi (1997) used a releasing and re-catching method, that is, he released at first 300 butterflies into the field and then caught again released butterflies day after day for a week in the field. The result was remarkable, that is, the survival rate of the mimetic f. polytes was statistically significantly lower than those of the nonmimetic f. cyrus and male. This means that the mimetic f. polytes has received an apparent disadvantage for survival in Taketomi-jima Is. in comparison with the nonmimetic f. cyrus and male.

Males Prefer the Non-mimic f. cyrus to the Mimic f. polytes?
The mimetic female f. polytes might be less distinguishable than the nonmimetic female f. cyrus by the male, because the form f. polytes resembles a different species P. aristolochiae, while the form f. cyrus resembles the male of the same species. This means that f. cyrus could get more chances to mate with the male than f. polytes. In reality, Uesugi (1997) counted mating times of young butterflies of both f. cyrus and f. polytes just after emergence. The result was striking to show that f. polytes had no approach from the male and no count, while f. cyrus could mate with the male. From the viewpoint of making offspring, the non-mimic f. cyrus has an advantage over the mimic f. polytes. Thus, viewing from the preference of the male to females in the species, the mimicry is not always beneficial.

Physiological Life Span of Two Forms f. cyrus and f. polytes
Ohsaki (2005) measured physiological life spans of flying three types (the male, the non-mimic f. cyrus, and the mimic f. polytes) in Itami City Museum of Insects, Japan, and found that the order of the life span length is as follows: f. cyrus > the male > f. polytes. The problem is why the life span of the mimic f. polytes is the shortest. To explain the reason, Ohsaki proposed a hypothesis that the mimic f. polytes pays an additional cost for producing red colored pigments (carotenoid) by activating biochemical reaction networks in cells. The hypothesis coincides with the result that the bigger the red colored area on the wing that the butterfly has, the shorter the physiological life span that the butterfly has.
On the other hand, by using P. polytes caught in Okinawa Island, Kinjyo (2000) showed that the physiological life span of f. cyrus and f. polytes is dependent on the feeding condition. For example, there is no statistically significant difference in the life span when both f. cyrus and f. polytes are fed in a good condition, under which honey is given every day (i.e., the mean life span, 23.5 days for f. cyrus and 22 days for f. polytes, respectively). When butterflies are fed in other conditions, under which plain water or no honey is given, mean life spans of the non-mimic f. cyrus become 8.8 days (with plain water) and 5.3 days (with no food), and those of the mimic f. polytes become 7.8 days (with plain water) and 4.6 days (with no food), respectively. This result shows that in nature, when environmental condition or food conditions (e.g., long rainy weather, dry weather) become worse, the life span of both f. cyrus and f. polytes becomes shorter rapidly. In any case, the result shows that the life span of f. cyrus is somewhat longer than that of f. polytes.
For the last several years, Sekimura measured the life span of both f. cyrus and f. polytes in the laboratory and found that in a regular condition at 25 C room temperature, there was no statistically significant difference in the life span of both f. cyrus and f. polytes under different feeding conditions (Fig. 12.3).

Mathematical Model of Three ODEs for Population Dynamics of P. polytes with Intraspecific Competition
In order to analyze field data summarized in Sect. 12.2, we present here a mathematical model for population dynamics of P. polytes. The model system consists of three ordinary differential equations (ODEs), in which variables are population densities of three butterflies, the unpalatable butterfly P. aristolochiae and two female forms of P. polytes. The male of P. polytes is not included directly in the system but included indirectly through parameter values such as intraspecific competition coefficients.
We denote four population densities of the model P. aristolochiae, the mimic f. polytes, the non-mimic f. cyrus, and predator as n 1 , n 2 , n 3 , and p, respectively. Changes of the densities in time become the following three ODEs: Phylogenetic life span of two female forms f. cyrus and f. polytes under the feeding condition. The horizontal axis indicates the feeding condition. "A" is the condition under which Calpis (or Calpico) water (i.e., Japanese milk-based soft drink) is given twice a day; "B" is the condition under which the Calpis water with doubling dilution is given twice a day; "C" is the condition under which only water is given twice a day; "D" is the condition under which nothing is given after hatching. The left side (orange-colored) column of each condition corresponds to the average life span of f. cyrus, and the right side (blue colored column) is that of f. polytes. The number upon each column indicates the average life span. For example, the number 16.50 upon the left column of the condition "A" is the average life span (days) of 10 (f. cyrus) individuals. Numbers of individuals used in the experiment are as follows: numbers of (f. cyrus) individuals, 10 for the feeding condition "A," 7 for "B," 14 for "C," and 14 for "D," respectively, and numbers of (f. polytes) individuals: 13 for the feeding condition "A," 7 for "B," 10 for "C," and 10 for "D," respectively. According to our statistical analysis of the data, there is no statistically significant difference in the physiological life span between two forms f. cyrus and f. polytes for all feeding conditions A-D. In contrast, the data analysis has made it clear that there exists statistically significant effect of the feeding condition (A-D) on difference in the life span where the population density of predator p is given as a fixed parameter value. In the system Eqs. (12.1), (12.2), and (12.3), we evaluated the density effect or saturation effect by using carrying capacities K 1 , K 2 , K 3 (¼K 2 ) for P. aristolochiae, f. polytes, and f. cyrus, respectively. We assume that K 2 , K 3 have the same value since both f. polytes and f. cyrus are females in the same species. Growth rates of three butterflies are denoted by r 1 , r 2 , r 3 , respectively. In Eqs. (12.2) and (12.3), we introduced an intraspecific competition between f. polytes and f. cyrus through competition coefficients (a 23 , a 32 ). The intraspecific competition effect includes competitions for resources such as nectar and indirectly the male as noted in Sect. 12.3.1.2. The term (n 2 /(n 1 + n 2 )) multiplied by p in Eqs. (12.1) and (12.2) represents the effect of mimicry, that is, the negative density effect implying that the increase in the density of f. polytes n 2 causes the decrease in both densities of P. aristolochiae and f. polytes, n 1 and n 2 . Parameters α , β 2 , β 3 represent difference in the predation rate among P. aristolochiae, P. polytes f. polytes, and f. cyrus, respectively, and it would be reasonable to assume the inequality α < β 2 , β 3 , because P. aristolochiae is an unpalatable butterfly species. In Sect. 12.3.1.1, we noted the result on beak marks by predators showing that nonmimetic females were selectively attacked, while males, mimetic females, and model butterflies were attacked less. This fact is evaluated mathematically in the second term of Eqs. (12.2) and (12.3) by multiplying p by (n 2 /(n 1 + n 2 )) (<1) for the mimic f. polytes, while by 1 for the non-mimic f. cyrus.
We put one more assumption on the survival of f. cyrus, r 3 > βp, i.e., the growth rate of f. cyrus is larger than the predation rate. Note that the population n 3 tends to zero when r 3 βp, since dn 3 /dt 0.
[Result C1-1] The change of relative abundance (RA) of the mimic in the female of P. polytes with respect to the carrying capacity K 1 of the model is positive, that is, The change of the advantage index (AI) with respect to the carrying capacity K 1 of the model is positive: Simulation results on temporal changes of all related quantities (a) The blue line denotes the population density n 3 of the non-mimic f. cyrus, the green line is the density n 2 of the mimic f. polytes, the purple line is the RA, the orange line is the AI, and the ocher line denotes the density n 1 of the model P. aristolochiae, respectively. Parameter values used in numerical simulations are as follows, and the stability condition for the positive equilibrium is satisfied with these parameter values: r 1 ¼ 0:5, r 2 ¼ 1:0, r 3 ¼ 2:0; K 1 ¼ 50, K 2 ¼ 50; a 23 ¼ 23 0:5, a 32 ¼ 0:35; p ¼ 0:5; α ¼ 0:6, β 2 ¼ β 3 ¼ 1:0; n 10 ¼ 0:01, n 20 ¼ 7:2, n 30 ¼ 34 initial values ð Þ (b) Numerical simulation results of the positive dependence of the RA on the AI. Parameter values used in the numerical simulation are all the same as in Fig. 12.4a From inequalities (12.4) and (12.5), we see that the ratio of change of the relative abundance (RA) to change of the advantage index (AI) with respect to the carrying capacity K 1 of the model is also positive: The inequality (12.6) provides the analytical evidence for the field record on the positive dependence of the RA on the AI in Sakishima Islands in Fig.12.5b (Uesugi 1992) noted in Sect. 12.2.2. Parameter values used in numerical simulations of Figs. 12.4, 12.5a and 12.5b are chosen as follows, so as to satisfy the existence and stability condition for the equilibrium of all the three population densities E 123 ¼ ( n 1 , n 2 , n 3 ), r 1 ¼ 0.5,   Fig. 12.4a, and the stability condition for the positive equilibrium is satisfied with these parameter values: r 1 ¼ 0:5, r 2 ¼ 1:0, r 3 ¼ 2:0; K 1 ¼ 50, K 2 ¼ 50; a 23 ¼ 23 0:5, a 32 ¼ 0:35; p ¼ 0:5; α ¼ 0:6, β 2 ¼ β 3 ¼ 1:0; n 10 ¼ 0:01, n 20 ¼ 7:2, n 30 ¼ 34 initial values ð Þ 12.4.1.2 Case 2: r 2 ¼ r 3 and β 2 > β 3 This is the case corresponding to experimental results noted in the last paragraphs of Sects. 12.3.1.3 and 12.3.1.1, that is, (a) there is no statistically significant difference in the life span of both f. cyrus and f. polytes (i.e., r 2 ¼ r 3 ), and (b) the survival rate of mimetic f. polytes is statistically significantly lower than that of nonmimetic f. cyrus (i.e., β 2 > β 3 ) (this result was obtained by an experiment in Taketomi-jima Is. (Uesugi 1997), where the model P. aristolochiae and the mimetic female f. polytes are both absent) (Uesugi 1991(Uesugi , 1992. [Result C2-1] First consider the system Eqs. (12.1), (12.2), and (12.3) under the condition n 10 ¼ 0 , n 20 > 0 and n 30 > 0. The setting of initial parameter values corresponds to the situation of Taketomi-jima Is. The uniqueness of the solution of the system implies that n 1 (t) ¼ 0 , n 2 (t) > 0 and n 3 (t) > 0 for any t > 0. Then the solution n 2 (t) and n 3 (t) satisfies The horizontal axis shows the advantage index (AI), which is defined as the number ratio of the model to all the related butterflies (the model and two female forms of P. polytes). The vertical axis is the relative abundance (RA). Solid circle ( ) represents the RA corresponding to the AI on each Island of 7 Islands (i.e., Taketomi-jima Is., Tarama-jima Is., Miyako-jima Is., Ishigaki-jima Is., Hateruma-jima Is., Obama-jima Is., and Iriomote-jima Is.) in Sakishima Islands (From Uesugi 1992). The solid line (i.e., Fig. 12.4b) represents numerical simulation results by use of the Eqs. (12.1), (12.2), and (12.3), which clearly show the positive dependence of the RA on the AI with a convexity. The inequality (6), which is an analytical result of the system equations, provides the theoretical basis on the positive dependence. Parameter values used in the numerical simulation are all the same as in Fig. 12.5a. dn 2 dt ¼ n 2 r 2 À β 2 p À r 2 n 2 þ a 23 n 3 K 2 ð12:7Þ dn 3 dt ¼ n 3 r 3 À β 3 p À r 3 n 3 þ a 32 n 2 K 3 , ð12:8Þ which is a traditional Lotka-Volterra competition model. Suppose that r i > β i p for i ¼ 2 , 3 (otherwise species i always tends to zero). It is easy to show that: 1. n 2 (t) ! 0 for any n 20 > 0 and n 30 > 0, when a 23 > 1 À β 2 p=r 2 1 À β 3 p=r 3 and a 32 < 1 À β 3 p=r 3 1 À β 2 p=r 2 ð12:9Þ 2. n 2 (t) ! 0 for some n 20 > 0 and n 30 > 0, when a 23 > 1 À β 2 p=r 2 1 À β 3 p=r 3 and a 32 > 1 À β 3 p=r 3 1 À β 2 p=r 2 : ð12:10Þ Mathematically we can prove that the equilibrium point E 3 ¼ 0, K 3 r 3 ðr 3 À β 3 pÞ is globally stable for system (12.7) and (12.8) under the condition (12.9). Also E 3 and E 2 ¼ K 2 r 2 r 2 À β 3 p ð Þ ; 0 are locally stable under (12.10). Note that under the parameter values for case 1, a 23 < 1 À β 2 p=r 2 1 À β 3 p=r 3 and a 32 > 1 À β 3 p=r 3 1 À β 2 p=r 2 ð12:11Þ are satisfied, and f. polytes and f. cyrus can coexist at the positive equilibrium point when P. aristolochiae is absent.
[Result C2-2] Now let us choose the parameter values for case 2 as r 1 ¼ 0.5 , , n 10 ¼ 0.01 , n 20 ¼ 1 , and n 30 ¼ 34. This case corresponds to an expectation of the population dynamics of butterflies after introducing the model P. aristolochiae into Taketomi-jima Is. The different choices between this and case 1 are r 2 ¼ r 3 ¼ 1.0, β 2 > β 3 , and the smaller initial value for f. polytes, since the above parameters satisfy condition (12.9) and only f. cyrus remains when P. aristolochiae is absent. Figure 12.6a shows that the introduction of the model yields the stable coexistence of P. aristolochiae, f. polytes, and f. cyrus. Further Fig. 12.6b shows the positive dependence of the RA and the AI.
[Result C2-3] The above results imply that the qualitative properties obtained for case 1 still hold true for case 2.
12.4.1.3 Case 3: r 2 ¼ r 3 and β 2 ¼ β 3 (¼β ) Now consider case 3. We adopt the parameter values a 23 ¼ 1.1 , a 32 ¼ 0.1 , β ¼ 1 and the same values for the remaining as case 2. Similar case 2, the parameters satisfy condition (12.9). Figure 12.7 shows the similar property as Fig. 12.6. Finally, we note mathematical results in a compact way in Sect. 12.4 as follows, Let us define the right-hand sides of the inequality in (12.9), (12.10), and (12.11) as a c 23 ¼ 1 À β 2 p=r 2 1 À β 3 p=r 3 and a c 32 ¼ 1 À β 3 p=r 3 1 À β 2 p=r 2 : ð12:12Þ It is easy to check that a c 23 ¼ a c 32 ¼ 1 for case 3 and a c 23 < 1, a c 32 > 1 for case 1 and 2.Since the condition a 23 > a c 23 implies that f. polytes goes extinct under the competition with f. cyrus (n 2 (t) ! 0) when P. aristolochiae is absent, both cases 1 and 2 enlarge the possibility for the extinction of f. polytes when P. aristolochiae is absent. The results obtained in Sect. 12.4 show that f. polytes can coexist with f. cyrus under the invasion of P. aristolochiae even f. polytes has disadvantageous property as case 1 (relatively small intrinsic growth rate) or case 2 (relatively small survival rate under that predation). Note that the same result can be obtained the case where β 2 r 2 > β 3 r 3 .

Summary and Discussions
Based on new experimental results, we presented an extension of the model for population dynamics of the mimetic butterfly P. polytes in Sakishima Islands, Japan (Sekimura et al. 2014). We introduced here difference in the predation rate between two female forms f. polytes and f. cyrus by parameters β 2 , β 3 , respectively. The new model system (12.1), (12.2), and (12.3) still includes three major effects: (a) self-density effect by carrying capacity; (b) mimetic effect, that is, the probability that a predator attacks prey on an encounter is proportional to the relative frequency of the mimic among the model and the mimic; and (c) intraspecific competition between two forms f. polytes and f. cyrus for resources such as nectar and, indirectly, the male. As to growth rate and predation rate, we took account two cases (r 2 < r 3 ) and (r 2 ¼ r 3 ) into consideration and assumed the inequality α < β 2 , β 3 among predation rates of the model P. aristolochiae, f. polytes, and f. cyrus, respectively.
Using mathematical analysis and computer simulations of the system equations, we extended the possibility of the logical relationship on the population dynamics Fig. 12.7 Simulation results on temporal changes of all related quantities (a) The blue line denotes the population density n 3 of the non-mimic f. cyrus, the green line is the density n 2 of the mimic f. polytes, the purple line is the RA, the orange line is the AI, and the ocher line denotes the density n 1 of the model P. aristolochiae, respectively. Parameter values used in numerical simulations are as follows, and the stability condition for the positive equilibrium is satisfied with these parameter values: r 1 ¼ 0:5, r 2 ¼ 1:0, r 3 ¼ 1:0; K 1 ¼ 50, K 2 ¼ 50; a 23 ¼ 23 1:1, a 32 ¼ 0:1; p ¼ 0:5; α ¼ 0:6, β 2 ¼ β 3 ¼ 1:0; n 10 ¼ 0:01, n 20 ¼ 1, n 30 ¼ 34 initial values ð Þ (b) Numerical simulation results of the positive dependence of the RA on the AI. Parameter values used in the numerical simulation are all the same as in Fig. 12.7a of the mimetic butterfly P. polytes in Sakishima Islands. In particular, we discussed both temporal change in the relative abundance (RA) and variation in the RA in Sakishima Islands. We estimated conditions for existence of equilibrium solutions of butterfly populations by making a comparison between experimental results and mathematical analyses of the model equations. Our results show that one of key factors to understand field data is the carrying capacity K 1 of the model in each island. The positive dependence of the RA on the AI originates from the result that changes of both the relative abundance (RA) and the advantage index (AI) with respect to the carrying capacity K 1 are positive.
The results in Sect. 12.4 have shown that both cases with respect to production rate: (r 2 < r 3 ) and (r 2 ¼ r 3 ) could be possible to reproduce experimental data on the population dynamics of the mimetic butterfly P. polytes in Sakishima Islands. The first case: (r 2 < r 3 ) means that the mimicry of P. polytes requires a kind of genetic change in production rate of the mimetic form f. polytes to reproduce the data. On the other hand, the second case (r 2 ¼ r 3 ) means that in order to reproduce the data, changes in ecological factors such as intraspecific competition coefficients (a 23 , a 32 ) and predation rate β 2 , β 3 are required for the mimicry of P. polytes without any genetic change. We think that it is not enough to determine at the moment which case is the real case that occurred in Sakishima Islands, because experiments on butterflies in both the field and the laboratory are somewhat subtle. We hope that much reliable experiments will be done to understand the reality in the future.
Finally, we hope that the mathematical analysis and computer simulations in the paper provide the theoretical basis on the female limited polymorphism of P. polytes in Sakishima Islands. Uesugi K (1996)  Open Access This chapter is licensed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this chapter are included in the chapter's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the chapter's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.