Behavioral consequences of exotic host plant adoption and the differing roles of male harassment on female movement in two checkerspot butterflies

Disease from leaves to landscapes: viral hotspots are determined by spatial arrangement and phytochemistry of host plants in specialist caterpillars

Although infectious diseases play a critical role in population regulation, our knowledge of complex drivers of disease for insects is limited. We conducted a field study on Baltimore checkerspot caterpillars (Euphydryas phaeton), chemical specialists on plants containing iridoid glycosides (IGs), to investigate the roles of host plant, phytochemistry, ontogeny and spatial associations in determining viral prevalence. We analysed individuals for viral presence and loads, quantified leaf IG concentrations from their native and novel host plants, and sequestered IGs in caterpillars. We found proximate caterpillar groups had greater similarity in infection prevalence, with areas of high prevalence indicating viral hotspots. Underlying variation in host plant chemistry corresponded to differences in viral prevalence. Furthermore, we used structural equation modeling to examine causal drivers of infection prevalence and loads. Advanced ontogeny was associated with increased viral prevalence and loads, as well as decreased sequestration of IGs. Infection loads were lower on the novel host plant, but prevalence was slightly higher, partially explained by decreased sequestration of IGs. Altogether, our findings reveal that spatial proximity, ontogeny, host plant species and secondary phytochemistry can all contribute to structuring infection risk, and thus offer insight into causal drivers of disease prevalence in complex plant-insect systems.

Plant Pathogen Invasion Modifies the Eco-Evolutionary Host Plant Interactions of an Endangered Checkerspot Butterfly

Simple Summary

Taylor’s checkerspot butterfly is a critically endangered species of northwestern North America that has become dependent on an exotic food plant, English plantain, which was acquired over a century ago. In the mid-2000s, a non-native plant pathogen from Europe, invaded Taylor’s checkerspot populations causing English plantain leaves to die in the winter when Taylor’s checkerspot larvae are feeding. We characterized butterfly and larval food plant (native and non-native) timing in Oregon and Washington populations and discovered that the invasive plant disease is active for ~60 days when larvae are feeding in January, February and March. Only one native larval foodplant, the annual Collinsia parviflora, can provide food for caterpillars during the time the plantain disease is common. However, this plant is rare in Taylor’s checkerspot habitat and may only be suitably timed to Washington checkerspot populations. Other native perennial larval food plants (Castillejalevisecta and likely C. hispida) do not appear suitably timed to provide resources throughout the entire Taylor’s checkerspot lifecycle in the low-elevation English plantain dependent populations. Understanding and accounting for the plant population disease dynamics is essential for the long-term conservation of Taylor’s checkerspot butterfly.

Abstract

New plant pathogen invasions typified by cryptic disease symptoms or those appearing sporadically in time and patchily in space, might go largely unnoticed and not taken seriously by ecologists. We present evidence that the recent invasion of Pyrenopeziza plantaginis (Dermateaceae) into the Pacific Northwest USA, which causes foliar necrosis in the fall and winter on Plantago lanceolata (plantain), the primary (non-native) foodplant for six of the eight extant Taylor’s checkerspot butterfly populations (Euphydryas editha taylori, endangered species), has altered eco-evolutionary foodplant interactions to a degree that threatens butterfly populations with extinction. Patterns of butterfly, larval food plant, and P. plantaginis disease development suggested the ancestral relationship was a two-foodplant system, with perennial Castilleja spp. supporting oviposition and pre-diapause larvae, and the annual Collinsia parviflora supporting post-diapause larvae. Plantain, in the absence of P. plantaginis disease, provided larval food resources throughout all butterfly life stages and may explain plantain’s initial adoption by Taylor’s checkerspot. However, in the presence of severe P. plantaginis disease, plantain-dependent butterfly populations experience a six-week period in the winter where post-diapause larvae lack essential plantain resources. Only C. parviflora, which is rare and competitively inferior under present habitat conditions, can fulfill the post-diapause larval feeding requirements in the presence of severe P. plantaginis disease. However, a germination timing experiment suggested C. parviflora to be suitably timed for only Washington Taylor’s checkerspot populations. The recent invasion by P. plantaginis appears to have rendered the ancestrally adaptive acquisition of plantain by Taylor’s checkerspot an unreliable, maladaptive foodplant interaction.

Preference, performance, and chemical defense in an endangered butterfly using novel and ancestral host plants

Adoption of novel host plants by herbivorous insects can require new adaptations and may entail loss of adaptation to ancestral hosts. We examined relationships between an endangered subspecies of the butterfly Euphydryas editha (Taylor’s checkerspot) and three host plant species. Two of the hosts (Castilleja hispida, Castilleja levisecta) were used ancestrally while the other, Plantago lanceolata, is exotic and was adopted more recently. We measured oviposition preference, neonate preference, larval growth, and secondary chemical uptake on all three hosts. Adult females readily laid eggs on all hosts but favored Plantago and tended to avoid C. levisecta. Oviposition preference changed over time. Neonates had no preference among host species, but consistently chose bracts over leaves within both Castilleja species. Larvae developed successfully on all species and grew to similar size on all of them unless they ate only Castilleja leaves (rather than bracts) which limited their growth. Diet strongly influenced secondary chemical uptake by larvae. Larvae that ate Plantago or C. hispida leaves contained the highest concentrations of iridoid glycosides, and iridoid glycoside composition varied with host species and tissue type. Despite having largely switched to a novel exotic host and generally performing better on it, this population has retained breadth in preference and ability to use other hosts.

Mechanism matters: the cause of fluctuations in boom-bust populations governs optimal habitat restoration strategy

Many populations exhibit boom-bust dynamics in which abundance fluctuates dramatically over time. Past research has focused on identifying whether the cause of fluctuations is primarily exogenous, e.g., environmental stochasticity coupled with weak density dependence, or endogenous, e.g., over-compensatory density dependence. Far fewer studies have addressed whether the mechanism responsible for boom-bust dynamics matters with respect to at-risk species management. Here, we ask whether the best strategy for restoring habitat across a landscape differs under exogenously vs. endogenously driven boom-bust dynamics. We used spatially explicit individual-based models to assess how butterfly populations governed by the two mechanisms would respond to habitat restoration strategies that varied in the level of resource patchiness, from a single large patch to multiple patches spaced at different distances. Our models showed that the restoration strategy that minimized extinction risk and boom-bust dynamics would be markedly different depending on the governing mechanism. Exogenously governed populations fared best in a single large habitat patch, whereas for endogenously driven populations, boom-bust dynamics were dampened and extinction risk declined when the total restored area was split into multiple patches with low to moderate inter-patch spacing. Adding environmental stochasticity to the endogenous model did not alter this result. Habitat fragmentation lowered extinction risk in the endogenously driven populations by reducing their growth rate, precluding both "boom" phases and, more importantly, "bust" phases. Our findings suggest that (1) successful restoration will depend on understanding the causes of fluctuations in at-risk populations, (2) the level and pattern of spatiotemporal environmental heterogeneity will also affect the ideal management approach, and (3) counterintuitively, for at-risk species with endogenously governed boom-bust dynamics, lowering the intrinsic population growth rate may decrease extinction risk.

Minimum area requirements for an at-risk butterfly based on movement and demography

Determining the minimum area required to sustain populations has a long history in theoretical and conservation biology. Correlative approaches are often used to estimate minimum area requirements (MARs) based on relationships between area and the population size required for persistence or between species' traits and distribution patterns across landscapes. Mechanistic approaches to estimating MAR facilitate prediction across space and time but are few. We used a mechanistic MAR model to determine the critical minimum patch size (CMP) for the Baltimore checkerspot butterfly (Euphydryas phaeton), a locally abundant species in decline along its southern range, and sister to several federally listed species. Our CMP is based on principles of diffusion, where individuals in smaller patches encounter edges and leave with higher probability than those in larger patches, potentially before reproducing. We estimated a CMP for the Baltimore checkerspot of 0.7-1.5 ha, in accordance with trait-based MAR estimates. The diffusion rate on which we based this CMP was broadly similar when estimated at the landscape scale (comparing flight path vs. capture-mark-recapture data), and the estimated population growth rate was consistent with observed site trends. Our mechanistic approach to estimating MAR is appropriate for species whose movement follows a correlated random walk and may be useful where landscape-scale distributions are difficult to assess, but demographic and movement data are obtainable from a single site or the literature. Just as simple estimates of lambda are often used to assess population viability, the principles of diffusion and CMP could provide a starting place for estimating MAR for conservation.

Low relative error in consumer-grade GPS units make them ideal for measuring small-scale animal movement patterns

Consumer-grade GPS units are a staple of modern field ecology, but the relatively large error radii reported by manufacturers (up to 10 m) ostensibly precludes their utility in measuring fine-scale movement of small animals such as insects. Here we demonstrate that for data collected at fine spatio-temporal scales, these devices can produce exceptionally accurate data on step-length and movement patterns of small animals. With an understanding of the properties of GPS error and how it arises, it is possible, using a simple field protocol, to use consumer grade GPS units to collect step-length data for the movement of small animals that introduces a median error as small as 11 cm. These small error rates were measured in controlled observations of real butterfly movement. Similar conclusions were reached using a ground-truth test track prepared with a field tape and compass and subsequently measured 20 times using the same methodology as the butterfly tracking. Median error in the ground-truth track was slightly higher than the field data, mostly between 20 and 30 cm, but even for the smallest ground-truth step (70 cm), this is still a signal-to-noise ratio of 3:1, and for steps of 3 m or more, the ratio is greater than 10:1. Such small errors relative to the movements being measured make these inexpensive units useful for measuring insect and other small animal movements on small to intermediate scales with budgets orders of magnitude lower than survey-grade units used in past studies. As an additional advantage, these units are simpler to operate, and insect or other small animal trackways can be collected more quickly than either survey-grade units or more traditional ruler/gird approaches.

Apparent power-law distributions in animal movements can arise from intraspecific interactions

Lévy flights have gained prominence for analysis of animal movement. In a Lévy flight, step-lengths are drawn from a heavy-tailed distribution such as a power law (PL), and a large number of empirical demonstrations have been published. Others, however, have suggested that animal movement is ill fit by PL distributions or contend a state-switching process better explains apparent Lévy flight movement patterns. We used a mix of direct behavioural observations and GPS tracking to understand step-length patterns in females of two related butterflies. We initially found movement in one species (Euphydryas editha taylori) was best fit by a bounded PL, evidence of a Lévy flight, while the other (Euphydryas phaeton) was best fit by an exponential distribution. Subsequent analyses introduced additional candidate models and used behavioural observations to sort steps based on intraspecific interactions (interactions were rare in E. phaeton but common in E. e. taylori). These analyses showed a mixed-exponential is favoured over the bounded PL for E. e. taylori and that when step-lengths were sorted into states based on the influence of harassing conspecific males, both states were best fit by simple exponential distributions. The direct behavioural observations allowed us to infer the underlying behavioural mechanism is a state-switching process driven by intraspecific interactions rather than a Lévy flight.