Broad-Scale Latitudinal Variation in Female Reproductive Success Contributes to the Maintenance of a Geographic Range Boundary in Bagworms (Lepidoptera: Psychidae)

Background Geographic range limits and the factors structuring them are of great interest to biologists, in part because of concerns about how global change may shift range boundaries. However, scientists lack strong mechanistic understanding of the factors that set geographic range limits in empirical systems, especially in animals. Methodology/Principal Findings Across dozens of populations spread over six degrees of latitude in the American Midwest, female mating success of the evergreen bagworm Thyridopteryx ephemeraeformis (Lepidoptera: Psychidae) declines from ∼100% to ∼0% near the edge of the species range. When coupled with additional latitudinal declines in fecundity and in egg and pupal survivorship, a spatial gradient of bagworm reproductive success emerges. This gradient is associated with a progressive decline in local abundance and an increased risk of local population extinction, up to a latitudinal threshold where extremely low female fitness meshes spatially with the species' geographic range boundary. Conclusions/Significance The reduction in fitness of female bagworms near the geographic range limit, which concords with the abundant centre hypothesis from biogeography, provides a concrete, empirical example of how an Allee effect (increased pre-reproductive mortality of females in sparsely populated areas) may interact with other demographic factors to induce a geographic range limit.


Introduction
Understanding how species' geographic distributions arise and are maintained constitutes one of the central goals of ecology. The 'abundant center' hypothesis from biogeography [1] predicts that local population density should decline as one moves from the core of a species' distribution toward the outer fringes, but many ecological mechanisms could give rise to such a pattern. Indeed across species, a broad array of causal factors are known to influence the positions and characteristics of geographic boundaries, but most studies, and especially those dealing with terrestrial animals, have evaluated only one factor at a time in isolation from other determinants [2]. Even among insects, where spatially replicated populations are often more tractable than in other animals, few studies have documented broad scale variation in reproductive success, nor how such variation may limit distribution range [3,4]. This is unfortunate because spatial variation in birth rate (reproduction) is probably the most critical determinant of geographic range boundaries [2].
Premature mortality of adult females (when females die before they lay their full complement of eggs) has long been known as a determinant of insect population dynamics [5,6,7], but the demographic consequences of reproductive failure at low popu-lation density (''demographic'' Allee effect) [8] have been documented for few insect species [9,10]. Theoretical models predict that a demographic Allee effect can contribute to maintenance of distributional range both with [11] and without [12] a strong environmental gradient. However, neither scenario has received strong empirical support.
Efforts to understand the structure, maintenance, and dynamics of animal range boundaries in a synthetic way are currently hamstrung by the lack of a model system for which multiple demographic parameters can be concurrently estimated in natural populations starting in the interior of a species range and moving out toward the edge of distribution [2]. Here, we answer the call for a model system for the study of geographic range limits. The bagcentred lifestyle of bagworms (Lepidoptera: Psychidae) makes them ideal animals for investigating geographic variation in demography. Multiple components of female fitness can be assessed using preand postmortem dissection of bags, including mortality during the pupal stage, timing of adult stage, mating success of adults, fecundity, overwintering survival of eggs, and reproductive output (Table 1). We provide herein a clear demonstration of how the interplay among a variety of demographic factors, including a striking spatial gradient in reproductive success, contributes to bagworm' geographic range boundary.

Study system
The bagworm Thyridopteryx ephemeraeformis (Haworth) is a univoltine, polyphagous moth widespread in the United States. Throughout its range, T. ephemeraeformis is broadly distributed as a pest in urban and agricultural landscapes on ornamental trees, predominantly juniper (Juniperus sp) and arborvitae (Thuja occidentalis).
Thyridopteryx ephemeraeformis possesses a suite of life history traits that make it an ideal candidate for a holistic approach to understanding what factors set and maintain range limits. Females are flightless as adults and reproduce within their bags, two traits that greatly facilitate studies of lifetime reproductive success and spatial population dynamics [13,14]. The bags are conspicuous on their host plant and infestations tend to occur in discrete patches on isolated plants or clusters of plants [14], thus facilitating sampling of local population even at low population density. Populations can be sampled along a broad latitudinal range (32-42uN) in the Midwest, but the species features a distinct geographic limit corresponding to northern Indiana that apparently has been stable for decades [15].
First instars construct a self-enclosing bag from host-plant material and enlarge this bag throughout their development. Upon completion of feeding, larvae tightly attach their bag to the host plant to pupate. Adults emerge in the fall. Males, which are typical winged moths, actively forage for sexually receptive females. The females are paedomorphic (neotenous), flightless, and do not leave their bag before the end of their life. Females attract mates during a 'calling stage' in which they disseminate setae impregnated with pheromone. Shortly after mating, the female oviposits a single clutch of eggs inside her pupal case and bag [16,17]; upon oviposition, the females drop to the ground and die. Females that fail to mate do not oviposit and eventually die within their bag, usually outside of the pupal case. The eggs laid by mated females overwinter inside the maternal bag, and neonates emerge in the spring.

Reproductive success of females for the 2008 reproductive season
Bagworm bags were sampled in March and April 2009, before the hatching of neonates. Study sites were located by driving through Indiana, Kentucky, and Tennessee along a 155 km wide north-south corridor (Fig. 1) and inspecting junipers for the presence of bagworms. Approximately 50 bags were collected on infested junipers at each of 110 sites. The bags were dissected to determine the mating status of females and the weight of egg masses. Postmortem assessments of the bags were conducted to evaluate the mating status of females based on the presence or absence of eggs (mated and unmated females, respectively) inside the pupal cases of females that had emerged as functional adults [18]. Adult emergence was diagnosed using as criteria the anterior split of the female pupal case and the presence of pheromoneimpregnated setae in the lower portion of the bag.
At the time of sampling in early spring, all egg masses laid by mated females appeared healthy (whitish colour, smooth shape).
To determine whether the eggs had survived the winter, pupal cases with egg masses were individually marked by location, kept in SoloH cups on a laboratory bench, and monitored daily to determine hatching of early instars. Early instars hatched over a 12 week period from the time of collection to 1 June. After a period of 7 days without emergence, the remaining egg masses that did not yield live larvae were visually inspected. All unhatched eggs had shrunk and turned black, indicative of overwintering mortality [19]. For each site, the mating success of females and proportion of females that produced live progeny (i.e., those females whose eggs overwintered successfully to hatch) were estimated as described in table 1.  Reproductive success for the 2009 reproductive season Sampling was conducted in Indiana ( Fig. 2) throughout the emergence period of adults in late summer and early fall (26 sites for arborvitae and 24 sites for juniper). Sampling was initiated at the onset of pupation and terminated when all females had emerged (10 August to 18 November). Because the sites were distributed across a broad latitudinal range, they could not all be sampled on the same day. Each sampling interval lasted 3-4 days and gaps between sampling intervals lasted 5-10 days. Sampling at a given site ceased when all females had emerged. For each site and sampling interval, between 5 and 34 females (usually .10) were collected during the emergence period, depending on the relative availability of bagworms. Females were removed from their bag and classified as either in the pupal stage (subclassified as live or dead pupae) or emerged adults (subclassified as mated, unmated or calling females). The weight of egg masses laid by mated females was determined for different sites and sampling intervals. Female survival during the pupal stage, mating success, and reproductive success were evaluated at different sites using the equations listed in table 1.

Rate of extinction of local populations
Because a substantial number of bagworm larvae remain and develop on their natal host [14], it was assumed that host plants that harbored bagworms during the 2008 generation but had no live larvae in 2009 represent local extinction events. The study was conducted at 56 sites previously sampled in Indiana in March 2008 to determine the reproductive success of females. Each site consisted of a group of trees infested with bagworms that was at least 10 m away from other infested trees, with a distance between sites .2 km. Each site was sampled a second time in June 2009 to determine the presence (sustained infestation) or absence (local extinction) of live larvae on the juniper plants. For each site, an index of between-year reproductive output (RS) was tabulated taking into account the probability that females mated and the probability of overwintering survival of eggs. The probability of local extinction across sites was estimated for two classes of female reproductive output [RS = 0 (complete reproductive failure); RS .0 (some females successfully reproduced)] and four latitudinal classes [,39uN; 39-40uN; 40-41uN; .41uN].

Data analysis
Statistical analysis was conducted with the SAS statistical package (version 9.1, SAS Institute, Cary, NC). Partition of variance analysis was used to evaluate the variance associated with latitude and longitude (one-degree bands) and host plant (juniper or arborvitae). Linear and logistic regression was used to evaluate the effect of latitude on different parameters of fitness. Unless otherwise stated, all the relationships reported are highly significant (P,0.0001).

Results
Partition of variance analysis indicated that the variance component associated with latitude was consistently larger than that associated with longitude or host plant for all parameters, usually by a factor .5 (Table 2). This variance structure justified our use of latitude as a key variable across which we quantified the fitness parameters of female bagworms.
For the 2008 generation of bagworms, logistic regression revealed a significant latitudinal decline in female mating success and egg overwintering survival, with steep declines at latitudes corresponding to central-northern Indiana (above 39uN for mating success and above 41.5uN for overwintering survival). Female fecundity declined linearly with latitude ( Fig. 1).
For the 2009 generation, survival during the pupal stage, female mating success and egg biomass all declined linearly with latitude (Fig. 2). Pupal mortality was primarily associated with Hymenoptera and Diptera generalist parasitoids.
The variation in female reproductive output (mating success * egg biomass; Table 1) in 2008 and 2009 was evaluated using only latitudes above 38.4uN so that the data could be compared for different years. Analysis of covariance revealed a highly significant effect of latitude on reproductive success (F = 204.98, df = 1,107, P,0.0001), but no significance of year either alone (F = 0.58, Female fitness declined linearly with latitude in both years and was particularly low above 41uN (Fig. 3).
The proportion of infested trees declined non-linearly with latitude, exhibiting a steep decline above 41u N; no infested trees (out of 109 sampled) were observed above 42u N ( Fig. 3; Table 3). The abundance of potential host plants was relatively constant between 39-41uN and increased above 41uN, thus the northern range limit of bagworms cannot be attributed to a lack of potential host plants. Furthermore, the latitudinal decline in abundance is not due to interspecific competition because very few defoliators other than bagworms were observed on arborvitae or juniper in the study area.
Local extinction of populations between the 2008 and 2009 bagworm generations was observed at 9 of 56 sites (16.1%). The probability of extinction was 100% (N = 4) at sites where females experience complete reproductive failure (RS = 0); all these sites occurred above 41uN (Table 4). At sites where some females reproduced (RS .0), no extinction event was observed below 38uN, and the probability of local extinction was roughly constant further north (13.6-15.4%) ( Table 4). The high rate of extinction events above 41uN (6 of 17 sites, or 32.6%) was due to the high proportion of sites where females experienced complete reproductive failures (4 of 17 sites, or 23.5%) ( Table 4). Logistic regression revealed a significant increase in the rate of extinction (y) as a function of latitude (x) [y = e (41.8720.996x) /(1+e (41.8720.996x) ), x = 3.98, P = 0.046].

Discussion
Demographic Allee effects, defined as positive impacts of density on the total fitness of individuals (e.g., high rate of mortality and low mating success in sparse populations), have been hypothesized to strongly influence population dynamics and to help constrain geographic range boundaries [8,11,12]. Unfortunately, empirical data that could be used to understand how the interplay between spatial gradients and population demography influences the establishment and maintenance of geographic range limits are  (table 3). For comparison, the effect of latitude (x) on female reproductive output (y, in mg eggs; see Table 1  rare in terrestrial animals [8,9]. This lack of data stems, in part, from the difficulty of adequately sampling low population density toward the edge of a species' distributional range and also from the lack of a model animal system for which multiple demographic parameters can be concurrently estimated [2]. We report here extremely low mating success of female bagworms in undisturbed, natural populations toward the edge of the distribution range, including the occurrence of total mating failure (0% mated female) at some sites. Of particular interest are the coincident latitudinal decline in bagworm abundance and female mating success above 41uN, and the apparent robustness of the latitudinal trends in 2008 and 2009. Because restricted mobility of females constrains their mating ability [10], species with flightless females may be particularly susceptible to low female mating success at low population density (mate encounter Allee effect) which may in turn influence the species' distributional range, as reported in the gypsy moth, Lymantria dispar [20,21,22]. Low vagility may further influence the interface between climate change and geographic range limits [3,23], particularly when species are unable to keep pace with changing landscapes through dispersal [24]. The abundance of potential host plants per se does not set the bagworm's range limit, as indicated by the increasing abundance of junipers and arborvitae with latitude (Table 3) and the absence of latitudinal variation in foliar nutrient content of the two main host plants of bagworms, junipers and arborvitae [25]. Because of the limited dispersal ability of bagworms and the naturally fragmented distribution of their host plants in urban and rural landscapes, we suggest that local reproductive success of females helps drive regional persistence for the species in concert with other demographic factors. Indeed, several components of female fitness declined toward the edge of the bagworms' distributional range, including survival during the pupal stage, mating success, fecundity, and overwintering survival of progeny, resulting in an overall reduction in reproductive success of females at northern locations, and in extreme cases, total reproductive failure. The increasing probability of extinction of local populations with declining female reproductive output (Table 4) indicates that patchy bagworm populations toward the range limit are temporally unstable. Such a demographic structure would be consistent with an ''invasion pinning'' scenario in which an Allee effect limits spatial spread [12].
The 'abundant-centre' hypothesis proposes that in a species range, a larger percentage of individuals of a population will be present in the centre of the range, where conditions are more favorable. Likewise, the hypothesis proposes a reduced density near the edge of a species range due to the interplay between numerous biotic and abiotic aspects of the habitat that worsens or becomes more intense as a species range boundary is approached [26,27,28]. In some cases, researchers question the validity of the abundant-centre hypothesis on the ground that apparent empirical support stems from reduced sampling near the edge of species range [27,28]. However, in other cases -such as the present studythe sampling designs for assessing population density across space are robust, and evidence for an abundant centre (and low density edges) is strong [29,30,31]. Predictive models suggest that the interplay between dispersal and demography can result in species that are 2 to 30 times denser in the centre of the range than at the edges [32]. If a species is going to be subject to Allee type dynamics, the likelihood for or intensity of such effects would be much greater in the vicinity of the range boundaries where densities are lower.
Several mechanisms may be simultaneously at work to restrict the reproductive output of female bagworms toward the edge of the distribution range, including an increased abundance of generalist pupal parasitoids (the most common natural enemies of bagworm pupae), low summer temperatures and short growing seasons (which together restrict the body size and fecundity of females; 33), low population density (which constrains female mating success), and low winter temperatures (which elevate egg mortality). Teasing apart the relative importance of these factors for the maintenance of the bagworms' geographic range limit will require a detailed population model that is parameterized from field data.