Abstract
We combined two models to investigate the theoretical dynamics of five exotic and native blowfly species in response to the Allee effect by using demographic parameters estimated from experimental populations. Most of the results suggest stabilization of dynamic behavior in response to the Allee effect. However, the results depended on the magnitude of the demographic parameters of each species, and also indicated chaotic fluctuations and limit cycles. The results are discussed in the context of larval aggregation, an important biological process for blowflies, which naturally incorporates the Allee effect.
Calliphoridae; ecological model; population dynamics
ECOLOGY, BEHAVIOR AND BIONOMICS
Allee effect in exotic and introduced blowflies
H Serra^{I}; Mis Costa^{II}; Wac Godoy^{III}
^{I}Univ Federal da Grande Dourados, Faculdade de Educação, Dourados, MS, Brasil
^{II}Coordenação de Sistemas e Controle, Lab Nacional de Computação Científica, MCT, Petrópolis, RJ, Brasil
^{III}Depto de Entomologia e Acarologia, Escola Superior de Agricultura "Luiz de Queiroz", USP, Piracicaba, SP, Brasil
^{Correspondence} Correspondence: Wesley A C Godoy Depto de Entomologia e Acarologia Escola Superior de Agricultura "Luiz de Queiroz" 13418900, Piracicaba, SP, Brasil wacgodoy@esalq.usp.br
ABSTRACT
We combined two models to investigate the theoretical dynamics of five exotic and native blowfly species in response to the Allee effect by using demographic parameters estimated from experimental populations. Most of the results suggest stabilization of dynamic behavior in response to the Allee effect. However, the results depended on the magnitude of the demographic parameters of each species, and also indicated chaotic fluctuations and limit cycles. The results are discussed in the context of larval aggregation, an important biological process for blowflies, which naturally incorporates the Allee effect.
Keywords: Calliphoridae, ecological model, population dynamics
Introduction
Larval aggregation is an important biological mechanism that takes place during the initial period of development in blowflies (Goodbrod & Goff 1990, Reis et al 1999). After the eggs hatch, firstinstar larvae tend to remain aggregated inside the food substrate, in order to increase the efficiency of the feeding process by secreting saliva and proteolytic enzymes (Goodbrod & Goff 1990, Reis et al 1999). These enzymes act together with bacteria present in the food substrate to decompose food (Hobson 1932, Mackerras 1933, Mackerras & Freney 1933, Goodbrod & Goff 1990).
This process is essential to create initial conditions for the larvae to exploit food, and has direct implications for the survival and fecundity of the species (Goodbrod & Goff 1990, Gião & Godoy 2007). It is a specific case of the Allee effect, because it requires a minimum number of larvae to decompose the food substrate and make the population viable (Wertheim et al 2002).
The Allee effect can be understood as a positive relationship between fitness and density of conspecifics, providing advantages to the individual (Stephens et al 1999). It reflects the mechanism of isolation, resulting in a benefit of aggregation, since for some densities the positivedensity dependence dominates, but for other densities it does not (Stephens et al 1999). For example, the Allee effect may operate in mating systems, not only of obligate cooperators but also of other species (Boukal & Berec 2002). In these systems, the spatial distribution generally determines the number of females that males can monopolize (Allee 1931, 1938, Boukal & Berec 2002). Other studies have discussed the diversity of mechanisms of the Allee effect (Amarasekare 1998, Dennis 1989, Courchamp et al 1999), including beneficial effects such as mitigation of predation, antipredator vigilance, social thermoregulation, reduction of inbreeding, and exploitation of resources (Stephens et al 1999, Liebhold & Tobin 2008).
The exotic blowfly species (Diptera: Calliphoridae) Chrysomya albiceps (Wiedemann 1819), C. megacephala (Fabricius 1794), and C. putoria (Wiedemann 1830) were introduced into Brazil about 30 years ago, and have displaced native species such as Cochliomyia macellaria (Fabricius 1775) and, to a lesser extent, Lucilia eximia (Wiedemann 1819) (Guimarães et al 1978, Silva et al 2003). This episode has motivated several studies on the population dynamics of introduced and native blowfly species, in an attempt to better understand their ecological patterns of fluctuation as well as the potential impact of biological invasion on the structure of the local fauna (Reis et al 1996, Serra et al 2007). The most recent studies include investigations of population dynamics by connecting laboratory experiments with mathematical models, focusing on density dependence, spatial structure, and stochasticity (Castanho et al 2006, Serra et al 2007).
Blowfly populations have been investigated with experiments and mathematical modelling since the classical studies on resource limitation carried out by Nicholson (1954, 1957) and revisited by Gurney et al (1980). The results obtained by Nicholson (1954, 1957) have been regularly analyzed because of the cyclic tendency found in the blowfly populations. The causes of cycles in insect populations have received considerable attention from ecologists and entomologists because they are often associated with densitydependent mechanisms (Hassell 1975, Muller 1986). A densitydependent mathematical model developed by Prout & McChesney (1985) to study population dynamics in experimental populations of Drosophila melanogaster has been successfully employed in studies of blowfly population dynamics (Godoy et al 2001, Godoy 2007, Serra et al 2007).
The model has raised important points in the context of blowfly biological invasions. The analysis indicates that the exotic species introduced into the Americas from the Old World exhibit dynamic behavior described by a twopoint limit cycle, characterized by a fluctuation between two points. The native species show a monotonic stable equilibrium, described by a constant population size (Godoy et al 2001, Silva et al 2003).
Previous investigations on the abundance of exotic and native blowfly species in the Americas reveals that if on the one hand, exotic species have successfully expanded their geographic area, on the other hand, native species have become much less abundant in all areas surveyed (Baumgartner & Greenberg 1984, Moura et al 1997, Centeno et al 2002, Carvalho et al 2004). These observations indicate that the demographic attributes of the several species are strongly influenced by interaction mechanisms.
In spite of the importance of these findings, an important point about the invasion has received no attention until now. The colonization process in blowflies differs according to areas and species, i.e., the four introduced species have reached different places, at different times, and with very different abundances (Guimarães et al 1978, Souza & Linhares 1997, Tomberlin et al 2001, Wolff et al 2001, Rosati & VanLaerhoven 2007). The most interesting example is the invasion of Central and North America by C. rufifacies and the introduction of C. albiceps into South America (Wells & Sperling 1999). The two species are ecologically similar, but have not overlapped in the same areas (Wells & Sperling 1999). Even among the species introduced into South America, their distribution and colonization differ among areas (Prado & Guimarães 1982, Carvalho et al 2000).
There is no specific study yet designed to investigate the causes of these differences, but there is evidence indicating that they may be associated with ecological factors governed by density dependence, such as demography, and interactions (Wells & Sperling 1999, Godoy et al 2001). Density dependence is a very influential mechanism in populations and acts on lifehistory parameters in insects (Muller 1986). The theoretical structure of the model employed in the previous studies includes negative densitydependence and delay effects (Prout & McChesney 1985). However, the Allee effect, recognized as a positive densitydependent effect, is not incorporated into this formulation yet.
Fowler & Ruxton (2002) exploited the Hassell model, including the Allee effect, and assumed that the contribution of an individual to the next generation is given by the Hassell function multiplied by an Allee effect. The Hassell function describes the effect of population size on individual reproductive output, by considering the reproductive potential in the absence of competition, the population size associated with the carrying capacity of the habitat, and the strength and form of competition (Hassell 1975). The authors demonstrated the changes that population dynamics could undergo when a factor representing an Allee effect is introduced.
The Allee effect has conventionally been investigated in an attempt to analyze low densities of endangered species (Tobin et al 2009). However, exotic species are generally at low density during an invasion process, the conditions in which the Allee effect takes place (Tobin et al 2009). The mathematical modeling scenario previously studied in blowflies indicates that exotic species exhibit high demographic values, tending to exhibit cyclic dynamics, and native species trend toward stability (Godoy et al 2001, Serra et al 2007).
In the present study, we used the Prout & McChesney model combined with the Fowler & Ruxton Alleeeffect component to investigate the theoretical dynamic behavior of five native and exotic blowfly species, L. eximia, C. macellaria, C. megacephala, C. albiceps, and C. putoria. We believe that the Allee effect may impact the distribution of insects living in ephemeral substrates such as carcasses, influencing their ecological patterns of dynamic behavior important to their persistence.
Our main concern about this point is, could the Allee effect significantly influence the blowfly population dynamics, stabilizing or not populations of both exotic and native species? Native blowfly species in Brazil have apparently shown strong resistance to biological invasion since the exotic species introduction, in spite of their decline in population numbers (Serra et al 2007). This population performance could be associated with the stability previously mentioned.
The primary objective of this study was to investigate the theoretical influence of the Allee effect on the quantitative and qualitative behavior of native and exotic blowfly species. Secondly, we were interested to know how fecundity and survival can change the stability of population in association with different magnitudes of the Allee effect.
Material and Methods
Mathematical models
A mathematical model developed by Prout & McChesney (1985), combined with the model of Fowler & Ruxton (2002), was used to show how the theoretical variation of the strength of the Allee effect can influence the dynamic behavior of L. eximia, C. macellaria, C. albiceps, C. megacephala, and C. putoria. The Prout & McChesney (1985) model is based on a finite difference equation that models population dynamics, taking into account the number of immatures, eggs, or larvae in succeeding generations, n_{t+1} and n_{t}. The model also incorporates two densitydependent processes, fecundity (F) and survival (S), as a function of the density of immatures, n_{t} , which were estimated from experimental populations (Table 1). The recursion is expressed by the nonlinear finite difference equation
where F* and S* are the intercepts in the exponential regression analysis of fecundity and survival as a function of larval density (Table 2). These parameters describe the theoretical values for maximum fecundity and survival respectively. The factor ½ indicates that only half of the population consists of adult females that contribute eggs to the next generation. The values of f and s are regression coefficients that estimate the slope of fecundity and survival on the density of immatures (Table 2).
The exponential function was used because it fits the blowfly data as well as or better than linear and hyperbolic functions (Godoy et al 1993, Von Zuben et al 1993). In addition, linear regression produces steeper slopes (in absolute magnitude) that in turn produce larger eigenvalues, which do not accurately describe the model dynamics at carrying capacity (Mueller 1986). Furthermore, the decrease in fecundity as a function of the density of immatures can be viewed biologically as a Poisson process that is described by an exponential function (Rodriguez 1989).
A component of fitness that increases with population size, termed the Allee effect (Fowler & Ruxton 2002), is expressed by
where A and γ are positive constants, and A is restricted to [0,1]. In this specific formulation, it was assumed that A is a strength parameter of the Allee effect and γ is the population subject to the Allee effect. Parameter a is a coefficient to scale the population size to the carrying capacity. In a blowfly context, A = 0 means no Allee effect, i.e., the larval density of blowflies is above a minimum necessary to make the population viable. With A = 1 the larval population is not able to become adults.
The new equation, including both the Allee effect and the demographic functions F and S, is given by
Using Matlab 7.0 (Hanselman & Littlefield 1997), simulations were run based on bifurcation theory (May & Oster 1976). Then, the results obtained were analyzed by bifurcation diagrams. A bifurcation diagram is a way of summarizing ranges of dynamic behaviors, since it describes locations and stability properties of periodic states (Vandermeer & Goldberg 2003). The variation of parameters may cause qualitative changes in populations, which can be viewed as different patterns of fluctuation or stability.
In bifurcation diagrams, the horizontal axis gives the parameter value (A in the current study) and the vertical axis represents the magnitudes of the steady state of the equation. Bifurcations represent new steady states (Vandermeer & Goldberg 2003). Two stable branches indicate the existence of a stable period 2 orbit, four stable branches a period 4 orbit, and so on. In the simulations we used the values obtained from Table 2, which were statistically estimated from data acquired in laboratory experiments (Table 1).
Experiments
Laboratory populations of L. eximia, C. macellaria, C. albiceps, C. megacephala, and C. putoria were founded from specimens collected from the Campus of the Universidade Estadual de Campinas, Campinas, São Paulo, Brazil. Adult flies were maintained at 25 ± 1°C in cages (30 x 30 x 48 cm) covered with nylon, and were fed water and sugar ad libitum.
Eggs were obtained by providing females with fresh beef liver. Larval densities (Table 1) were established using newly hatched larvae in 50 g of rearing medium prepared according to Leal et al (1982), with two replicates. Fecundity was estimated by randomly removing 30 females per replicate and counting the number of eggs per female. Survival was estimated as the number of adults emerging from each vial. The data presented in Table 1 were obtained from an experimental setting of increasing larval densities designed to estimate fecundity and survival as decreasing functions of the density.
With the results obtained from the experiments, exponential regressions were performed in order to analyze the decreasing functions of density, survival, and fecundity, and to obtain the parameters to be used in the Prout equation. The estimates of the parameters resulted in two theoretical values for fecundity (F) and survival (S). These values are the intercept of the regression. The regression coefficients (s) and (f) estimate the variation of survival and fecundity, respectively, as a function of density. The results of this analysis are shown in Table 2. It is important to emphasize that the F and S estimates were obtained in laboratory experiments designed to study intraspecific competition without considering the Allee effect. Then, the densities used in the experiments were suitable only to provide negative densitydependence. The Allee effect was taken into account only during simulations run with the combination between Prout & McChesney (1985) and Fowler & Ruxton (2002) models.
In order to investigate the magnitude of parameter values required to exhibit different dynamic behaviors, such as stable equilibrium, cycles, and unpredictable oscillations, arbitrary values of F and S were used in the simulations. The values were employed taking into account the maximum and minimum boundaries of survival and fecundity observed in natural populations (Ullyett 1950, Reigada & Godoy 2005, Gião & Godoy 2006, Riback & Godoy 2008). This procedure allowed us to determine the critical values for a change of equilibrium in each species. The constants a and γ were the same values employed by Fowler & Ruxton (2002), γ = 1 and a = 0.005.
The values used by Fowler & Ruxton (2002) guaranteed values that would produce chaotic dynamics in the equivalent Hassell model. Hence, a gradual introduction of these components resulted in stabilization of the dynamics, a typical result found with an influence from the Allee effect. An evaluation of the sensitivity in all parametrical space of the Fowler & Ruxton parameters ( γ and a ) was performed in order to understand how this can extrapolate into qualitative changes in the equilibrium dynamics of experimental populations of the five blowfly species, and consequently what are the consequences for the population viability.
The sensitivity analysis was conducted by employing the bifurcation theory following the procedures indicated by May & Oster (1976), with results indicating no qualitative change in all parametrical space of γ and a, and viability for all blowfly species. The results found guarantee that the Fowler & Ruxton (2002) parameter values applied to γ and a in this study do not influence qualitatively the dynamic behaviors observed in the current study and the quantitative changes observed do not drive the blowfly populations to extinction.
Results
Bifurcation diagrams obtained from experimentally estimated parameters
Lucilia eximia and C. macellaria exhibited a monotonic stable equilibrium within the entire parametric space of A, with L. eximia trending to extinction for values of A higher than 0.8 (Fig 1a, b). Chrysomya albiceps showed a twopoint limit cycle for initial values of A, and a stable equilibrium above 0.15 (Fig 1c). Chrysomya megacephala and C. putoria exhibited a twopoint limit cycle during the entire parametric space of A (Fig 1d, e).
Bifurcation diagrams obtained from varying parameters (native species)
The neighborhood of fecundity and survival was investigated by varying parameter values in the simulations, in order to determine what demographic values are necessary to find the exact point of bifurcation, taking into account the Allee effect. For high values of fecundity in L. eximia, twice the experimentally obtained value, the bifurcation diagram showed a twopoint limit cycle for the initial values of A and a stabilizing effect for values of A higher than 0.25 (Fig 2a). A very slight increase of fecundity in C. macellaria, from 18.33 to 18.65, was sufficient to show the bifurcation point, where the twopoint limit cycle changed to stable equilibrium (Fig 2b). A decrease in survival from the estimate value of 1 (100% survival) produced stable equilibrium, but an increase in the survival of C. macellaria from 0.78 to 0.8 was also sufficient to show the bifurcation, trending to stability for A values closer to 1 (Fig 2c).
Bifurcation diagrams obtained from varying parameters (exotic species)
Increasing the fecundity of C. albiceps from 27.11 to 48, it was possible to observe a bifurcation diagram containing two and fourpoint limit cycles (Fig 2d). The same result was obtained by changing survival from 0.56 to 1 (Fig 2e). The increase of fecundity in C. megacephala and C. putoria produced destabilizing effects for increasing values of A, with results indicating a succession of bifurcations leading to unpredictable fluctuations, alternating with periodicity windows (Fig 2f, h). Nevertheless, C. putoria showed periodic behaviors at the parametric extremes of A, differing from C. megacephala, which showed periodicity only for initial values of A.
Because the survival values of C. megacephala and C. putoria were very close to 1 (0.91 and 0.97 respectively), an increase in survival would produce no qualitative change in the dynamic behavior of the two species. Otherwise, a decrease in survival in C. megacephala and C. putoria produced a monotonic stable equilibrium (Fig 2g, i), leading both species to the same result obtained for C albiceps in the simulations with experimental data.
Discussion
In L. eximia and C. macellaria, only the arbitrary increase in fecundity resulted in the appearance of two behaviors, a twopoint limit cycle and a stable equilibrium. This result can be explained by considering the values of Λ for both species, which are more distant from 1 compared to the values found for C. albiceps (Reis et al 1996, Silva et al 2003). The Λ value has a critical change when it is equal to 1, i.e., when Λ is lower than 1 the result is a monotonic stable equilibrium, but if Λ is higher than 1 the result is a limit cycle (Prout & McChesney 1985). Lucilia eximia requires at least twice its estimated value to exhibit changes, since its experimentally obtained fecundity was very low. Nevertheless, natural populations of this species show fecundity values significantly higher than the numbers found in laboratory, justifying the arbitrary change in the simulations (Gião & Godoy 2006).
Cochliomyia macellaria had been previously recognized as a species capable of exhibiting only a monotonic stable equilibrium, but in our simulations it requires only a slight change in the fecundity (F) value to exhibit the twopoint limit cycle, i.e. from F = 18.33 to F = 18.65. This result was found because although Λ is less than 1, it is very close to 1, differing from C. albiceps, a species that exhibits Λ slightly higher than 1 (Reis et al 1996, Godoy et al 2001).
By using the experimentally obtained demographic values, only C. albiceps showed more than one type of dynamic behavior, i.e., a twopoint limit cycle and a stable equilibrium. This result could be understood as a response to its eigenvalue (Λ), which is very close to 1, but higher than 1 (Godoy et al 2001), making the species prone to changes in terms of equilibrium. This result suggests that the Allee effect was responsible for the change in Λ and consequently for the appearance of the two different equilibria. These findings make sense if the results are compared to the previous studies by Godoy et al (2001), where only one type of dynamic behavior was found.
Chrysomya albiceps differs from the other species in several respects. It is cannibal and predator during larval stage, and its predatory habit has strongly influenced the current Brazilian necrophagous fauna (Faria et al 1999). Its strategy to obtain food is probably different from the other blowflies, in response to larval habits, which could affect the larval aggregation process, making the species more prone to the Allee effect than other species (Faria et al 2007). The survival of C. albiceps was probably influenced by cannibalism. This habit is certainly the reason for the difference between its survival pattern and the other blowfly species. Cannibalistic behavior could explain, at least in part, the poor performance of C. albiceps for intraspecific densities.
The variation of fecundity in C. albiceps, i.e. the increase in fecundity, resulted in only two types of behavior, differing from C. megacephala and C. putoria, which showed a cascade of bifurcations. There is a very interesting point here. The changes in fecundity resulted in different patterns of bifurcation for C. albiceps and the other two Chrysomya species. In C. albiceps, cycles appeared first, subsequently trending to stability. For C. megacephala and C. putoria the result was the reverse, starting with periodic cycles, but changing to more complex cycles and chaos, including the occurrence of periodicity windows. Thus, for C. albiceps the increase of the Allee effect is a stabilizing strength, but for C. megacephala and C. putoria it is not.
These results probably reflect the differences among these species with respect to demographic parameter values. However, the differences can only be observed by analyzing the combination of the two mathematical models. In spite of their biological differences, C. albiceps seems closer to L. eximia and C. macellaria than the other Chrysomya species, at least in terms of susceptibility to changes from cyclic to stable behavior. Two levels of stability were observed by Godoy et al (2001), with implications for the equilibrium dynamics of these species. High fecundity and survival values led to eigenvalues higher than 1.3, exhibited by C. megacephala and C. putoria, and low fecundity and survival values led to eigenvalues lower than 1.3, exhibited by L. eximia, C. macellaria, and C. albiceps (Godoy et al 2001). The combination found between high and low demographic values, associated with different values of A (Allee effect), may result in different patterns of bifurcation, explaining the results observed for exotic and native blowflies in this report. For survival, all species underwent a stabilizing effect with increasing values of this parameter, showing a monotonic stable equilibrium.
Although the laboratory data do not wholly reproduce the parameter values in natural populations, the results found represent a good approach for the blowfly system. In addition, the results from previous studies designed to investigate the variation of demographic parameters in natural populations of blowflies show a spectrum of parameter values compatible with the numbers used in this study (Reigada & Godoy 2005, Gião & Godoy 2006, Riback & Godoy 2008).
Blowflies are insects that easily exhibit cyclic fluctuations (Nicholson 1954, 1957, Von Zuben et al 1993), a behavior characterized by demographic parameter values of high magnitude (EdelsteinKeshet 1988). Unfortunately there is no comparable longterm study focusing on realworld observations of population patterns. However, the spectrum of oscillations in natural populations of exotics is probably wider than native species because of the magnitude of their survival and fecundity in natural conditions (Ullyett 1950, Reigada & Godoy 2005, Gião & Godoy 2006, Riback & Godoy 2008).
In a previous study, Godoy et al (1996) showed that the succession of increases both in fecundity and survival results in destabilizing the population equilibrium. On the other hand, several studies have shown that the Allee effect influences populations with a stabilizing effect on their dynamics (Scheuring 1999, Fowler & Ruxton 2002). It was clear in our study that the Allee effect may invert tendencies toward both stabilizing and destabilizing populations, and this depends on the magnitude of and combinations among demographic parameters.
The inversions caused by the Allee effect address questions of fundamental importance in the context of biological invasions, since the ability of invader organisms to colonize new niches depends essentially on the magnitude of their demographic parameters and the strength of the Allee effect. On the other hand, for native species, the presence of other species that share the same resources may cause serious population declines, decreasing population sizes and increasing the possibility of local extinction (Reis et al 1996, Serra et al 2007).
In blowflies, the level of larval aggregation is directly associated with the females' oviposition behavior. Blowfly females usually lay eggs in batches of 100300, frequently among eggs of several different species (Smith 1986). The larval distribution in carcasses depends on the dispersal of the females, and can have implications for their population dynamics (De Jong 1976, Blackith & Blackith 1990, Turchin 1998). However, the initial larval density, the factor that is subject to the Allee effect, also depends on the existing assemblage of blowfly species, because the previous presence of some species, as for example the predatory larvae of C. albiceps, significantly influences the diversity and abundance of blowflies (Gião & Godoy 2007). Hence, for blowflies, the Allee effect depends on a set of factors important to the larval aggregation, such as: carrion size, attractiveness of substrates, ovarian development stage in females, and number of eggs laid.
The strategies of exploiting food resources differ significantly in blowflies (Goodbrood & Goff 1990, Godoy et al 1996, Faria et al 1999). Intra and interspecific interactions in blowflies have been investigated, and have suggested different larval aggregation patterns, probably associated with the type of food item and the stage of larval development (Reis et al 1999, Gião & Godoy 2007). We believe that the blowfly system deserves special attention, considering all these aspects, in an attempt to analyze which exogenous and endogenous forces might explain the inter and intraspecific variability in terms of population dynamics of exotic and native blowfly species.
We have shown that the differences in terms of population behavior found among blowfly species play an important role in determining dynamics, resulting in clear population patterns. Several studies have documented the causes of Allee effects and their association with biological invasions (Taylor & Hastings 2005). For invader insects such as Chrysomya species, the colonization and trophic interactions could be subject to different mechanisms of an Allee effect (Tobin et al 2009). In conclusion, our results show clearly that the Allee effect may stabilize or destabilize the dynamic behavior of blowflies, depending on the parametric space of demographic parameters.
Acknowledgments
This research was supported by a grant from the Fundação de Amparo à Pesquisa do Estado de São Paulo. Work by HS was supported by a postgraduate fellowship from the Conselho Nacional de Desenvolvimento Científico e Tecnológico. WACG was partially supported by a research fellowship from the Conselho Nacional de Desenvolvimento Científico e Tecnológico. The authors also thank Janet W. Reid for revising the English text. The manuscript was improved significantly by constructive criticisms from an anonymous reviewer.
Received 24 February 2010 and accepted 06 April 2010
Edited José R Trigo  UNICAMP
 Allee WC (1931) Animal aggregations: a study in general sociology. Chicago, University of Chicago Press, 431p.
 Allee WC (1938) The social life of animals. London, Heinemann, 293p.
 Amarasekare P (1998) Allee effects in metapopulation dynamics. Am Nat 152: 298302.
 Blackith RE, Blackith RM (1990) Insect infestations of small corpses. J Nat Hist 24: 699709.
 Baumgartner D L, Greenberg B (1984) The genus Chrysomya (Diptera: Calliphoridae) in the New World. J Med Entomol 21: 105113.
 Boukal DS Berec LK (2002) Singlespecies models of the Allee effect: extinction boundaries, sex ratios and mate encounters. J Theor Biol 218: 375394.
 Carvalho LML, Thyssen PJ, Goff ML, Linhares AX (2004) Observations on the succession patterns of necrophagous insects on a pig carcass in an urban area of Southeastern Brazil. Agg Inter J For Med Tox 5: 3339.
 Carvalho LML, Thyssen PJ, Linhares AX, Palhares FB (2000) A checklist of arthropods associated with carrion and human corpses in southeastern Brazil. Mem Inst Oswaldo Cruz 95: 135138.
 Castanho MJP, Magnago KF, Bassanezi RC, Godoy WAC (2006) Fuzzy subset approach in coupled population dynamics of blowflies. Biol Res 39: 341352.
 Centeno N, Maldonato M, Olivia A (2002) Seasonal patterns of arthropods occurring on sheltered and unsheltered pig carcasses in Buenos Aires province (Argentina). For Sci Int 126: 6370.
 Courchamp F, CluttonBrock T, Grenfell B (1999) Inverse density dependence and the Allee effect. TREE 14: 405410.
 De Jong GA (1976) A model of competition for food. I. Frequencydependent viabilities. Am Nat 110: 10131027.
 Dennis B (1989) Allee effects: population growth, critical density, and chance of extinction. Nat Res Mod 3: 481538.
 EdelsteinKeshet L (1988) Mathematical models in biology. New York, Random House, 586p.
 Faria LDB, Orsini L, Trinca LA, Godoy WAC (1999) Larval predation by Chrysomya albiceps on Cochliomyia macellaria, Chrysomya megacephala and Chrysomya putoria Entomol Exp Appl 90: 149155.
 Faria LDB, Reigada C, Trinca LA, Godoy WAC (2007) Foraging behaviour by an intraguild predator blowfly, Chrysomya albiceps (Diptera: Calliphoridae). J Ethol 25: 287294.
 Faria LDB, Trinca LA, Godoy WAC (2004) Cannibalistic behavior and functional response in Chrysomya albiceps (Diptera: Calliphoridae). J Insect Behav 17: 251261.
 Fowler MS, Ruxton GD (2002) Population dynamics consequences of Allee effects. J Theor Biol 215: 3946.
 Friedenberg NA, Powell JA, Ayres MP (2007) Synchrony's double edge: transient dynamics and the Allee effect in stagestructured populations. Ecol Lett 10: 564573.
 Gião JZ, Godoy WAC (2006) Seasonal population dynamics in Lucilia eximia (Wiedemann) (Diptera: Calliphoridae). Neotrop Entomol 35: 753756.
 Gião JZ, Godoy WAC (2007) Ovipositional behavior in predator and prey blowflies. J Insect Behav 20: 77.
 Godoy WAC (2007) Dynamics of blowfly populations. Func Ecosys Comm 1: 129139.
 Godoy WAC, Reis SF, Von Zuben CJ, Ribeiro OB (1993) Population dynamics of Chrysomya putoria (Wied.) (Dipt., Calliphoridae). J Appl Entomol 116: 163169.
 Godoy WAC, Von Zuben CJ, Reis SF, Von Zuben FJ (1996) Dynamics of experimental blowflies (Diptera: Calliphoridae): Mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations. Mem Inst Oswaldo Cruz 91: 641648.
 Godoy WAC, Von Zuben CJ, Reis SF, Von Zuben FJ (1997) The spatial dynamics of native and introduced blowflies (Diptera: Calliphoridae). J Appl Entomol 121: 305309.
 Godoy WAC, Von Zuben CJ, Reis SF, Von Zuben FJ (2001) Spatiotemporal dynamics and transition from asymptotic equilibrium to bounded oscillations in Chrysomya albiceps (Diptera: Calliphoridae). Mem Inst Oswaldo Cruz 96: 627634.
 Goodbood JR, Goff ML (1990) Effects of larval population density on rates of development and interactions between two species of Chrysomya (Diptera: Calliphoridae) in laboratory culture. J Med Entomol 27: 338343.
 Guimarães JH, Prado AP, Linhares AX (1978) Three newly introduced blowfly species in Southern Brazil (Diptera: Calliphoridae). Revta Bras Entomol 22: 5360.
 Gurney WSC, Blythe SP, Nisbet RM (1980) Nicholson's blowflies revisited. Nature 287: 1721.
 Hanselman D, Littlefield B (1997) The student edition of Matlab. Upper Saddle River, Prentice Hall, 429p.
 Hassell MP (1975) Densitydependence in singlespecies populations. J Anim Ecol 44: 283295.
 Hobson RP (1932) Studies on the nutrition on the blowfly larvae. III. The liquefaction of muscle. J Exp Biol 9: 359365.
 Leal TTS, Prado AP, Antunes AJ (1982) Rearing the larvae of the blowfly Chrysomya chloropyga (Wiedemann) (Diptera: Calliphoridae) on oligidic diets. Revta Bras Zool 1: 4144.
 Liebhold AM, Tobin PC (2008) Population ecology of insect invasions and their management. Ann Rev Entomol 55: 387408.
 Mackerras MJ (1933) Observations on the lifehistories, nutritional requirements and fecundity of blowflies. Bull Entomol Res 24: 353362.
 Mackerras MJ, Freney MR (1933) Observations on the nutrition of maggots of Australian blowflies. J Exp Biol 10: 237246.
 May RM, Oster GF (1976) Bifurcations and dynamic complexity in simple ecological models. Am Nat 110: 573599.
 Mueller LD (1986) Densitydependent rates of population growth: estimation in laboratory populations. Am Nat 128: 282293.
 Moura MO, Carvalho CJB, Monteiro ELA (1997) A preliminary analysis of insects of medicolegal importance in Curitiba, State of Paraná. Mem Inst Oswaldo Cruz 93: 269274.
 Nicholson AJ (1954) An outline of the dynamics of animal populations. Aust J Zool 2: 965.
 Nicholson AJ (1957) The selfadjustment of populations to change. Cold Spring Harbor Symp Quant Biol 22: 153173.
 Prado AP, Guimarães JH (1982) Estado atual de dispersão e distribuição do gênero Chrysomya RobineauDesvoidy na Região Neotropical (Diptera, Calliphoridae). Revta Bras Entomol 26: 225231.
 Prout T, McChesney F (1985) Competition among immatures affects their adult fertility: population dynamics. Am Nat 126: 521558.
 Reigada C, Godoy WAC (2005) Seasonal fecundity and body size in Chrysomya megacephala (Fabricius) (Diptera: Calliphoridae). Neotrop Entomol 34: 163168.
 Reis SF, Godoy WAC, Von Zuben CJ (1999) Larval aggregation and competition for food in experimental populations of Chrysomya putoria (Wied.) and Cochliomyia macellaria (F.) (Dipt., Calliphoridae). J Appl Entomol 123: 485489.
 Reis SF, Teixeira MA, Von Zuben FJ, Godoy WAC, Von Zuben CJ (1996) Theoretical dynamics of experimental populations of introduced and native blowflies (Diptera, Calliphoridae). J Med Entomol 33: 537544.
 Riback TIS, Godoy WAC (2008) Fecundity, body size and population dynamics of Chrysomya albiceps (Wiedemann, 1819) (Diptera: Calliphoridae). Braz J Biol 68: 131136.
 Rodriguez DJ (1989) A model of population dynamics for the fruit fly Drosophila melanogaster with density dependence in more than one life stage and delayed effects. J Anim Ecol 58: 349365.
 Rosa GS, Carvalho LR, Reis SF, Godoy WAC (2006) The dynamics of intraguild predation in Chrysomya albiceps Wied. (Diptera: Calliphoridae): interactions between instars and species under different abundances of food. Neotrop Entomol 35: 775780.
 Rosati J, VanLaerhoven SL (2007) New record of Chrysomya rufifacies (Macquart) in Canada: current distribution and the possibility of future range expansion. Can Entomol 139: 670677.
 Scheuring I (1999) Allee effect increases the dynamical stability of populations. J Theor Biol 199: 407414.
 Serra H, Silva ICR, Mancera PFA, Faria LDB, Von Zuben CJ, Von Zuben FJ, Reis SF, Godoy WAC (2007) Stochastic dynamics in exotic and native blowflies: an analysis combining laboratory experiments and a twopatch metapopulation model. Ecol Res 22: 686695.
 Silva ICR, Mancera PFA, Godoy WAC (2003) Population dynamics of Lucilia eximia (Diptera: Calliphoridae). J Appl Entomol 127: 26.
 Smith KGV (1986) A manual of forensic entomology. London, Oxford, University Printing House, 205p.
 Souza AM, Linhares AX (1997) Diptera and Coleoptera of potential forensic importance in Southeastern Brazil: relative abundance and seasonality. Med Vet Entomol 11: 812.
 Stephens PA, Sutherland WJ, Freckleton RP (1999) What is the Allee effect? Oikos 87: 185190.
 Taylor CM, Hastings A (2005) Allee effects in biological invasions. Ecol Lett 8: 895908.
 Tobin PC, Robinet C, Johnson DM, Whitmire SL, BjØrnstad ON, Liebhold AM (2009) The role of Allee effects in gypsy moth (Lymantria dispar (L.)) invasions. Popul Ecol 51: 373384.
 Tomberlin JK, Reeves W K, Sheppard DC (2001) First record of Chrysomya megacephala (Diptera: Calliphoridae) in Georgia, USA. Fla Entomol 84: 300301.
 Turchin P (1998) Quantitative analysis of movement: measuring and modelling population redistribution in animals and plants. Sunderland, Sinauer Associates, 396p.
 Ullyett GC (1950) Competition for food and allied phenomena in sheepblowfly populations. Phil Trans Royal Soc London Ser B Biol Sci 234: 77174.
 Vandermeer JH, Goldberg DH (2003) Population ecology: first principles. Princeton, Princeton University Press, 280p.
 Von Zuben CJ, Reis SF, Val JBR, Godoy WAC, Ribeiro OB (1993) Dynamics of a mathematical model of Chrysomya megacephala (Diptera: Calliphoridae). J Med Entomol 30: 443448.
 Wells JD, Sperling FAH (1999) Molecular phylogeny of Chrysomya albiceps and C rufifacies (Diptera: Calliphoridae). J Med Entomol 36: 222226.
 Wertheim B, Marchais J, Vet LEM, Dicke M (2002) Allee effect in larval resource exploitation in Drosophila: an interaction among density of adults, larvae, and microorganisms. Ecol Entomol 27: 608617.
 Wolff M, Uribe A, Ortiz A, Duque P (2001) A preliminary study of forensic entomology in Medellín, Colombia. For Sci Inter 120: 5359.
Correspondence:
Publication Dates

Publication in this collection
07 Nov 2011 
Date of issue
Oct 2011
History

Received
24 Feb 2010 
Accepted
06 Apr 2010