Monarch butterfly-breeding habitat restoration: how movement ecology research can inform best practices for site selection.

Population dynamics, persistence, and distribution are emergent properties of animal movement behavior and the spatial configuration of resources. Monarch butterflies are a vagile species with an open-population structure. Selecting locations for monarch butterfly- breeding habitat restoration that aligns with natural movement behavior will facilitate efficient habitat utilization across the landscape, increase realized fecundity, and ultimately support increases in the overwintering population size in Mexico. Obtaining and interpreting empirical movement and space-use data through field and laboratory studies are fundamental to this effort. To gain insights into population responses at larger, spatially explicit landscape scales, the results from empirical studies can be incorporated into simulation models. Together, empirical and simulation studies can inform options for creating functional connectivity of monarch butterfly-breeding habitats. Given currently available information, we synthesize studies for the eastern monarch butterfly to illustrate how an improved understanding of movement ecology can assist in planning conservation practices.

Monarch Butterfly Ecology, Behavior, and Vulnerabilities in North Central United States Agricultural Landscapes.

The North American monarch butterfly (Danaus plexippus) is a candidate species for listing under the Endangered Species Act. Multiple factors are associated with the decline in the eastern population, including the loss of breeding and foraging habitat and pesticide use. Establishing habitat in agricultural landscapes of the North Central region of the United States is critical to increasing reproduction during the summer. We integrated spatially explicit modeling with empirical movement ecology and pesticide toxicology studies to simulate population outcomes for different habitat establishment scenarios. Because of their mobility, we conclude that breeding monarchs in the North Central states should be resilient to pesticide use and habitat fragmentation. Consequently, we predict that adult monarch recruitment can be enhanced even if new habitat is established near pesticide-treated crop fields. Our research has improved the understanding of monarch population dynamics at the landscape scale by examining the interactions among monarch movement ecology, habitat fragmentation, and pesticide use.

Neonicotinoids can cause arrested pupal ecdysis in Lepidoptera.

Recently, we reported a novel mode of action in monarch butterfly (Danaus plexippus) larvae exposed to neonicotinoid insecticides: arrest in pupal ecdysis following successful larval ecdysis. In this paper, we explore arrested pupal ecdysis in greater detail and propose adverse outcome pathways to explain how neonicotinoids cause this effect. Using imidacloprid as a model compound, we determined that final-instar monarchs, corn earworms (Helicoverpa zea), and wax moths (Galleria mellonella) showed high susceptibility to arrested pupal ecdysis while painted ladies (Vanessa cardui) and red admirals (Vanessa atalanta) showed low susceptibility. Fall armyworms (Spodoptera frugiperda) and European corn borers (Ostrinia nubilalis) were recalcitrant. All larvae with arrested ecdysis developed pupal cuticle, but with incomplete shedding of larval cuticle and unexpanded pupal appendages; corn earworm larvae successfully developed into adults with unexpanded appendages. Delayed initiation of pupal ecdysis was also observed with treated larvae. Imidacloprid exposure was required at least 26 h prior to pupal ecdysis to disrupt the molt. These observations suggest neonicotinoids may disrupt the function of crustacean cardioactive peptide (CCAP) neurons, either by directly acting on their nicotinic acetylcholine receptors or by acting on receptors of inhibitory neurons that regulate CCAP activity.

Estimating Perceptual Range of Female Monarch Butterflies (Danaus plexippus) to Potted Vegetative Common Milkweed (Asclepias syriaca) and Blooming Nectar Resources.

Habitat loss in the summer breeding range contributes to eastern North American monarch (Lepidoptera: Nymphalidae) population decline. Habitat restoration efforts include increasing native prairie plants for adult forage and milkweed (Asclepias spp.) for oviposition and larval development. As the monarch is a vagile species, habitat establishment at a grain that matches the monarch perceptual range will facilitate efficient movement, decrease fitness costs of dispersal, and increase oviposition. We released 188 experimental monarch females 5, 25, 50, and 75 m downwind from potted milkweed and blooming forbs in 4-32 ha sod fields. Perceptual range was estimated from monarchs that flew towards and landed on the milkweed and forbs. Flight patterns of 49 non-experimental monarchs that landed on the resources were also observed. In our experimental, resource-devoid setting, wind-facilitated movement occurred most frequently. Monarchs performed direct displacement as evidenced by shallow turn angles and similarity of Euclidian and total distances traveled. We hypothesize similar monarch flight behavior when traveling over other resource-devoid areas, such as crop fields. Although the majority of experimental monarchs flew downwind, eight experimental and 49 non-experimental monarchs were observed flying upwind toward, and landing on, the potted resources from distances ranging from 3 to 125 m (mean = 30.98 m, median = 25 m, mode = 25 m). A conservative estimate of the perceptual range is 125 m, as longer distances cannot be precluded; however, the majority of observations were ≤50 m. Our findings suggest establishing habitat patches ~ 50 m apart would create functional connectivity across fragmented agricultural landscapes.

Evaluating toxicity of Varroa mite (Varroa destructor)-active dsRNA to monarch butterfly (Danaus plexippus) larvae.

Varroa mites (Varroa destructor) are parasitic mites that, combined with other factors, are contributing to high levels of honey bee (Apis mellifera) colony losses. A Varroa-active dsRNA was recently developed to control Varroa mites within honey bee brood cells. This dsRNA has 372 base pairs that are homologous to a sequence region within the Varroa mite calmodulin gene (cam). The Varroa-active dsRNA also shares a 21-base pair match with monarch butterfly (Danaus plexippus) calmodulin mRNA, raising the possibility of non-target effects if there is environmental exposure. We chronically exposed the entire monarch larval stage to common (Asclepias syriaca) and tropical (Asclepias curassavica) milkweed leaves treated with concentrations of Varroa-active dsRNA that are one- and ten-fold higher than those used to treat honey bee hives. This corresponded to concentrations of 0.025-0.041 and 0.211-0.282 mg/g leaf, respectively. Potassium arsenate and a previously designed monarch-active dsRNA with a 100% base pair match to the monarch v-ATPase A mRNA (leaf concentration was 0.020-0.034 mg/g) were used as positive controls. The Varroa mite and monarch-active dsRNA's did not cause significant differences in larval mortality, larval or pupal development, pupal weights, or adult eclosion rates when compared to negative controls. Irrespective of control or dsRNA treatment, larvae that consumed approximately 7500 to 10,500-mg milkweed leaf within 10 to 12 days had the highest pupal weights. The lack of mortality and sublethal effects following dietary exposure to dsRNA with 21-base pair and 100% base pair match to mRNAs that correspond to regulatory genes suggest monarch mRNA may be refractory to silencing by dsRNA or monarch dsRNase may degrade dsRNA to a concentration that is insufficient to silence mRNA signaling.

Oviposition Response of Monarch Butterfly (Lepidoptera: Nymphalidae) to Imidacloprid-Treated Milkweed.

Monarch butterfly (Danaus plexippus) populations have declined over the last two decades, attributable in part to declines in its larval host plant, milkweed (Asclepias spp.), across its breeding range. Conservation efforts in the United States call for restoration of 1.3 billion milkweed stems into the Midwestern landscape. Reaching this goal will require habitat establishment in marginal croplands, where there is a high potential for exposure to agrochemicals. Corn and soybean crops may be treated with neonicotinoid insecticides systemically or through foliar applications to provide protection against insect pests. Here, we investigate whether ovipositing monarchs discriminate against milkweed plants exposed to the neonicotinoid insecticide imidacloprid, either systemically or through foliar application. In our first experiment, we placed gravid females in enclosures containing a choice of two cut stems for oviposition: one in 15 ml of a 0.5 mg/ml aqueous solution of imidacloprid and one in 15 ml water. In a second experiment, females were given a choice of milkweed plants whose leaves were treated with 30 µl of a 0.825 mg/ml imidacloprid-surfactant solution or plants treated with surfactant alone. To evaluate oviposition preference, we counted and removed eggs from all plants daily for 3 d. We also collected video data on a subset of butterflies to evaluate landing behavior. Results indicate that neither systemic nor foliar treatment with imidacloprid influenced oviposition behavior in female monarchs. The implications of these findings for monarch conservation practices will be informed by the results of ongoing egg and larval toxicity studies.

Monarch Butterfly (Danaus plexippus) Life-Stage Risks from Foliar and Seed-Treatment Insecticides.

Conservation of North America's eastern monarch butterfly (Danaus plexippus) population would require establishment of milkweed (Asclepias spp.) and nectar plants in the agricultural landscapes of the north central United States. A variety of seed-treatment and foliar insecticides are used to manage early- and late-season pests in these landscapes. Thus, there is a need to assess risks of these insecticides to monarch butterfly life stages to inform habitat conservation practices. Chronic and acute dietary toxicity studies were undertaken with larvae and adults, and acute topical bioassays were conducted with eggs, pupae, and adults using 6 representative insecticides: beta-cyfluthrin (pyrethroid), chlorantraniliprole (anthranilic diamide), chlorpyrifos (organophosphate), imidacloprid, clothianidin, and thiamethoxam (neonicotinoids). Chronic dietary median lethal concentration values for monarch larvae ranged from 1.6 × 10-3 (chlorantraniliprole) to 5.3 (chlorpyrifos) µg/g milkweed leaf, with the neonicotinoids producing high rates of arrested pupal ecdysis. Chlorantraniliprole and beta-cyfluthrin were generally the most toxic insecticides to all life stages, and thiamethoxam and chlorpyrifos were generally the least toxic. The toxicity results were compared to insecticide exposure estimates derived from a spray drift model and/or milkweed residue data reported in the literature. Aerial applications of foliar insecticides are expected to cause high downwind mortality in larvae and eggs, with lower mortality predicted for adults and pupae. Neonicotinoid seed treatments are expected to cause little to no downslope mortality and/or sublethal effects in larvae and adults. Given the vagile behavior of nonmigratory monarchs, considering these results within a landscape-scale context suggests that adult recruitment will not be negatively impacted if new habitat is established in close proximity of maize and soybean fields in the agricultural landscapes of the north central United States. Environ Toxicol Chem 2021;40:1761-1777. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Conservation risks and benefits of establishing monarch butterfly (Danaus plexippus) breeding habitats close to maize and soybean fields in the north central United States: A landscape-scale analysis of the impact of foliar insecticide on nonmigratory monarch butterfly populations.

Establishing habitat in agricultural landscapes of the north central United States is critical to reversing the decline of North America's eastern monarch butterfly (Danaus plexippus) population. Insecticide use could create population sinks and threaten recovery. Discouraging habitat establishment within a 38-m zone around crop fields is a suggested risk mitigation measure. In Story County, Iowa, United States, this mitigation would discourage habitat establishment in 84% of roadsides and 38% of noncrop land. It is unclear if the conservation benefits from establishing habitat close to crop fields outweigh suppression of population growth owing to insecticide exposure. Consequently, monarch conservation plans require spatially and temporally explicit landscape-scale assessments. Using an agent-based model that incorporates female monarch movement and egg laying, the number and location of eggs laid in Story County were simulated for four habitat scenarios: current condition, maximum new establishment, moderate establishment, and moderate establishment only outside a 38-m no-plant zone around crop fields. A demographic model incorporated mortality from natural causes and insecticide exposure to simulate adult monarch production over 10 years. Assuming no insecticide exposure, simulated adult production increased 24.7% and 9.3%, respectively, with maximum and moderate habitat establishment and no planting restrictions. A 3.5% increase was simulated assuming moderate habitat establishment with a 38-m planting restriction. Impacts on adult production were simulated for six representative insecticides registered for soybean aphid (Aphis glycines) management. Depending on the frequency of insecticide applications over a 10-year period, simulated production increased 8.2%-9.3%, assuming moderate habitat establishment with no planting restrictions. Results suggest that the benefits of establishing habitat close to crop fields outweigh the adverse effects of insecticide spray drift; that is, metapopulation extirpation is not a concern for monarchs. These findings are only applicable to species that move at spatial scales greater than the scale of potential spray-drift impacts. Integr Environ Assess Manag 2021;17:989-1002. © 2021 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Prediction of mitochondrial genome-wide variation through sequencing of mitochondrion-enriched extracts.

Although mitochondrial DNA (mtDNA) haplotype variation is often applied for estimating population dynamics and phylogenetic relationships, economical and generalized methods for entire mtDNA genome enrichment prior to high-throughput sequencing are not readily available. This study demonstrates the utility of differential centrifugation to enrich for mitochondrion within cell extracts prior to DNA extraction, short-read sequencing, and assembly using exemplars from eight maternal lineages of the insect species, Ostrinia nubilalis. Compared to controls, enriched extracts showed a significant mean increase of 48.2- and 86.1-fold in mtDNA based on quantitative PCR, and proportion of subsequent short sequence reads that aligned to the O. nubilalis reference mitochondrial genome, respectively. Compared to the reference genome, our de novo assembled O. nubilalis mitochondrial genomes contained 82 intraspecific substitution and insertion/deletion mutations, and provided evidence for correction of mis-annotated 28 C-terminal residues within the NADH dehydrogenase subunit 4. Comparison to a more recent O. nubilalis mtDNA assembly from unenriched short-read data analogously showed 77 variant sites. Twenty-eight variant positions, and a triplet ATT codon (Ile) insertion within ATP synthase subunit 8, were unique within our assemblies. This study provides a generalizable pipeline for whole mitochondrial genome sequence acquisition adaptable to applications across a range of taxa.

Optimization of QuEChERS Method for Simultaneous Determination of Neonicotinoid Residues in Pollinator Forage.

Consistent with the large-scale use of pesticide seed treatments in U.S. field crop production, there has been an increased use of neonicotinoid-treated corn and soybean seed over the past decade. Neonicotinoids can move downwind to adjacent off-field pollinator habitats in dust from planting and/or move downslope to habitats in surface water. The extent of potential neonicotinoid exposure to pollinators from neonicotinoid movement into these adjacent pollinator habitats is unclear. Pollen and leaf tissue extractions were completed using a quick, easy, cheap, effective, rugged, and safe (QuEChERS) extraction procedure. Samples were subjected to a clean-up step using dispersive solid-phase extraction (dSPE) techniques prior to analysis. The compounds in the extracts were separated on a reversed-phase column with gradient elution and confirmed with tandem mass spectrometry. The extraction method showed acceptable recoveries of analytes ranging from 78.4 to 93.6% and 89.4 to 101% for leaf tissue and pollen, respectively. The method's detection limits ranged from 0.04 to 0.3 ng/g in milkweed leaf tissue and 0.04 to 1.0 ng/g in pollen. The method is currently being employed in ongoing studies surveying pollen from a diversity of forbs and milkweed leaves obtained from habitat patches established within fields with a history of using neonicotinoid-treated seeds.

Employing Very High Frequency (VHF) Radio Telemetry to Recreate Monarch Butterfly Flight Paths.

The overwintering population of eastern North American monarch butterflies (Danaus plexippus) has declined significantly. Loss of milkweed (Asclepias sp.), the monarch's obligate host plant in the Midwest United States, is considered to be a major cause of the decline. Restoring breeding habitat is an actionable step towards population recovery. Monarch butterflies are highly vagile; therefore, the spatial arrangement of milkweed in the landscape influences movement patterns, habitat utilization, and reproductive output. Empirical studies of female movement patterns within and between habitat patches in representative agricultural landscapes support recommendations for habitat restoration. To track monarch movement at distances beyond human visual range, we employed very high frequency radio telemetry with handheld antennae to collect movement bearings on a biologically relevant time scale. Attachment of 220-300 mg transmitters did not significantly affect behavior and flight capability. Thirteen radio-tagged monarchs were released in a restored prairie, and locations were estimated every minute for up to 39 min by simultaneous triangulation from four operators. Monarchs that left the prairie were tracked and relocated at distances up to 250 m. Assuming straight flights between locations, the majority of steps within the prairie were below 50 m. Steps associated with exiting the prairie exceeded 50 m with high directionality. Because butterflies do not fly in straight lines between stationary points, we also illustrate how occurrence models can use location data obtained through radio telemetry to estimate movement within a prairie and over multiple land cover types.

Assessing Field-Scale Risks of Foliar Insecticide Applications to Monarch Butterfly (Danaus plexippus) Larvae.

Establishment and maintenance of milkweed plants (Asclepias spp.) in agricultural landscapes of the north central United States are needed to reverse the decline of North America's eastern monarch butterfly (Danaus plexippus) population. Because of a lack of toxicity data, it is unclear how insecticide use may reduce monarch productivity when milkweed habitat is placed near maize and soybean fields. To assess the potential effects of foliar insecticides, acute cuticular and dietary toxicity of 5 representative active ingredients were determined: beta-cyfluthrin (pyrethroid), chlorantraniliprole (anthranilic diamide), chlorpyrifos (organophosphate), and imidacloprid and thiamethoxam (neonicotinoids). Cuticular median lethal dose values for first instars ranged from 9.2 × 10-3 to 79 µg/g larvae for beta-cyfluthrin and chlorpyrifos, respectively. Dietary median lethal concentration values for second instars ranged from 8.3 × 10-3 to 8.4 µg/g milkweed leaf for chlorantraniliprole and chlorpyrifos, respectively. To estimate larval mortality rates downwind from treated fields, modeled insecticide exposures to larvae and milkweed leaves were compared to dose-response curves obtained from bioassays with first-, second-, third-, and fifth-instar larvae. For aerial applications to manage soybean aphids, mortality rates at 60 m downwind were highest for beta-cyfluthrin and chlorantraniliprole following cuticular and dietary exposure, respectively, and lowest for thiamethoxam. To estimate landscape-scale risks, field-scale mortality rates must be considered in the context of spatial and temporal patterns of insecticide use. Environ Toxicol Chem 2020;39:923-941. © 2020 SETAC.

Application of the tissue residue approach in ecological risk assessment.

The objective of this work is to present a critical review of the application of the tissue residue approach (TRA) in ecological risk and/or impact assessment (ERA) of chemical stressors and environmental criteria development. A secondary goal is to develop a framework for integrating the TRA into ecological assessments along with traditional, exposure concentration-based assessment approaches. Although widely recognized for its toxicological appeal, the utility of the TRA in specific applications will depend on numerous factors, such as chemical properties, exposure characteristics, assessment type, availability of tissue residue-response data, and ability to quantify chemical exposure. Therefore, the decision to use the TRA should include an evaluation of the relative strengths, limitations, and uncertainties among exposure and residue-based methods for characterizing toxicological effects. Furthermore, rather than supplanting exposure concentration-based toxicity assessments, the TRA can be highly effective for evaluating and reducing uncertainty when used in a complementary manner (e.g., when evaluating multiple lines of evidence in field studies). To address limitations with the available tissue residue-response data, approaches for extrapolating residue-based toxicity data across species, tissues, and exposure durations are discussed. Some of these approaches rely on predicted residue-response relationships or toxicological models that have an implicit residue-response basis (e.g., biotic ligand model). Because risk to an organism is a function of both its exposure potential and inherent sensitivity (i.e., on a residue basis), bioaccumulation models will be required not only for translating tissue residue criteria into corresponding water and sediment criteria, but also for defining the most vulnerable species in an assemblage (i.e., highly exposed and highly sensitive species). Application of the TRA in ecological assessments and criteria development are summarized for bioaccumulative organic chemicals, TBT, and in situ bioassays using bivalve molluscs.

Overview of data and conceptual approaches for derivation of quantitative structure-activity relationships for ecotoxicological effects of organic chemicals.

The use of quantitative structure-activity relationships (QSARs) in assessing potential toxic effects of organic chemicals on aquatic organisms continues to evolve as computational efficiency and toxicological understanding advance. With the ever-increasing production of new chemicals, and the need to optimize resources to assess thousands of existing chemicals in commerce, regulatory agencies have turned to QSARs as essential tools to help prioritize tiered risk assessments when empirical data are not available to evaluate toxicological effects. Progress in designing scientifically credible QSARs is intimately associated with the development of empirically derived databases of well-defined and quantified toxicity endpoints, which are based on a strategic evaluation of diverse sets of chemical structures, modes of toxic action, and species. This review provides a brief overview of four databases created for the purpose of developing QSARs for estimating toxicity of chemicals to aquatic organisms. The evolution of QSARs based initially on general chemical classification schemes, to models founded on modes of toxic action that range from nonspecific partitioning into hydrophobic cellular membranes to receptor-mediated mechanisms is summarized. Finally, an overview of expert systems that integrate chemical-specific mode of action classification and associated QSAR selection for estimating potential toxicological effects of organic chemicals is presented.

An overview of the use of quantitative structure-activity relationships for ranking and prioritizing large chemical inventories for environmental risk assessments.

Ecological risk assessments for chemical stressors are used to establish linkages between likely exposure concentrations and adverse effects to ecological receptors. At times, it is useful to conduct screening risk assessments to assist in prioritizing or ranking chemicals on the basis of potential hazard and exposure assessment parameters. Ranking of large chemical inventories can provide evidence for focusing research and/or cleanup efforts on specific chemicals of concern. Because of financial and time constraints, data gaps exist, and the risk assessor is left with decisions on which models to use to estimate the parameter of concern. In this review, several methods are presented for using quantitative structure-activity relationships (QSARs) in conducting hazard screening or screening-level risk assessments. The ranking methods described include those related to current regulatory issues associated with chemical inventories from Canada, Europe, and the United States and an example of a screening-level risk assessment conducted on chemicals associated with a watershed in the midwest region of the United States.

Influence of suspended solids on acute toxicity of carbofuran to Daphnia magna: I. Interactive effects.

This study explored the effects on Daphnia magna from exposure to the pesticide carbofuran in combination with stress from suspended solids exposure. Our objective was to assess whether suspended solids affects the toxicodynamic response of D. magna to carbofuran. A series of laboratory experiments was performed where animals were exposed to carbofuran concentrations ranging from 0 to 160 microg/l in combination with suspended solids concentrations ranging from 0 to 10000 mg/l. In the absence of suspended solids, effects of carbofuran were dose dependent and resulted in an EC(50) of 92 microg/l. Exposure to suspended solids, up to extreme levels that may be encountered in the environment and in the absence of carbofuran, showed no measurable toxicity. When D. magna were exposed to a constant carbofuran concentration, the numbers of affected organisms increased with increasing suspended solids concentrations. At a suspended solids concentration of 1000 mg/l, the EC(50) for carbofuran was reduced by half to 45 microg/l. The relationship between the toxicity of carbofuran (microg/l) and the concentration of suspended solids (mg/l) can be described with the following equation: carbofuran EC(50)=72 exp(-0.00014 [suspended solids]). An analysis of the data indicates that this relationship is consistent with a potentiated toxicity mechanism rather than an additive model.

Influence of suspended solids on acute toxicity of carbofuran to Daphnia magna: II. An evaluation of potential interactive mechanisms.

It has been demonstrated that simultaneous exposure of Daphnia magna to suspended solids and a carbamate pesticide potentiates the toxic response to the pesticide. The toxicodynamics between these stressors were investigated to determine possible mechanisms of interaction. Three experimental series were conducted with D. magna to determine: the effect of food availability on carbofuran toxicity; the effect of food availability on jointly administered carbofuran and suspended solids; and changes in the magnitude of effects which can occur with suspended solids of different composition. These experiments demonstrated that both carbofuran toxicity and the joint toxicity of carbofuran and suspended solids to D. magna can be modulated by food availability. While it is clear that food dilution could contribute to energy stress, it appears likely that additional interactive processes also contribute to the observed synergism between carbofuran and suspended solids. Additionally, decomposed peat was shown to be less of a stressor to these pelagic invertebrates than inorganic subsoil.