Long-distance dispersal of the gypsy moth (Lepidoptera: Lymantriidae) facilitated its initial invasion of Wisconsin

Size-dependent flight capacity and propensity in a range-expanding invasive insect

For capital-breeding insects, all resources available for adult metabolic needs are accumulated during larval feeding. Therefore, body size at adult eclosion represents the total energetic capacity of the individual. For female capital breeders, body size is strongly correlated with lifetime fecundity, while in males, body size, which correlates with fitness, is less understood. In capital-breeding species with wingless, flightless, or dispersal-limited females, flight potential for male Lepidoptera has important implications for mate-finding and may be correlated with body size. At low population densities, failure to mate has been identified as an important Allee effect and can drive the success or failure of invasive species at range edges and in species of conservation concern. Th capital-breeding European subspecies of Lymantria dispar (L.), was introduced to North America in 1869 and now ranges across much of eastern North America. In L. dispar, females are flightless and mate-finding is entirely performed by males. We quantified male L. dispar flight capacity and propensity relative to morphological and physiological characteristics using fixed-arm flight mills. A range of male body sizes was produced by varying the protein content of standard artificial diets while holding other dietary components constant. Wing length, a proxy for body size, relative thorax mass, and forewing aspect were all important predictors of total flight distance and maximum speed. These results have important implications for mate-finding and invasion dynamics in L. dispar and may apply broadly to other capital-breeding insects.

Universal Mitochondrial Multi-Locus Sequence Analysis (mtMLSA) to Characterise Populations of Unanticipated Plant Pest Biosecurity Detections

Simple Summary

Agricultural and environmental sustainability requires effective biosecurity responses that prevent the establishment or spread of exotic insect pests. Understanding where new detections may have come from or if recurrent detections are connected contributes to this. Suitable population genetic markers use relatively rapidly evolving gene regions which render the PCR method species-specific at best. Because resource limitations mean these are pre-emptively developed for the highest risk species, populations of other exotic pests are unable to be characterised at the time. Here we have developed a generic method that is useful across species within the same taxonomic Order, including where there is little or no prior knowledge of their gene sequences. Markers are formed by concomitant sequencing of four gene regions. Sequence concatenation was shown to retrieve higher resolution signatures than standard DNA barcoding. The method is encouragingly universal, as illustrated across species in ten fly and 11 moth superfamilies. Although as-yet untested in a biosecurity situation, this relatively low-tech, off-the-shelf method makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for impromptu species-specific research and development.

Abstract

Biosecurity responses to post-border exotic pest detections are more effective with knowledge of where the species may have originated from or if recurrent detections are connected. Population genetic markers for this are typically species-specific and not available in advance for any but the highest risk species, leaving other less anticipated species difficult to assess at the time. Here, new degenerate PCR primer sets are designed for within the Lepidoptera and Diptera for the 3′ COI, ND3, ND6, and 3′ plus 5′ 16S gene regions. These are shown to be universal at the ordinal level amongst species of 14 and 15 families across 10 and 11 dipteran and lepidopteran superfamilies, respectively. Sequencing the ND3 amplicons as an example of all the loci confirmed detection of population-level variation. This supported finding multiple population haplotypes from the publicly available sequences. Concatenation of the sequences also confirmed that higher population resolution is achieved than for the individual genes. Although as-yet untested in a biosecurity situation, this method is a relatively simple, off-the-shelf means to characterise populations. This makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for unprepared species-specific research and development.

Evaluating the success of treatments that slow spread of an invasive insect pest

BACKGROUND

Treatments for the suppression and eradication of insect populations undergo substantial testing to ascertain their efficacy and safety, but the generally limited spatial and temporal scope of such studies limit knowledge of how contextual factors encountered in operational contexts shape the relative success of pest management treatments. These contextual factors potentially include ecological characteristics of the treated area, or the timing of treatments relative to pest phenology and weather events. We used an extensive database on over 1000 treatments of nascent populations of Lymantria dispar (L.) (gypsy moth) to examine how place-based and time-varying conditions shape the success of management treatments.

RESULTS

We found treatment success to vary across states and years, and to be highest in small treatment blocks that are isolated from other populations. In addition, treatment success tended to be lower in treatment blocks with open forest canopies, possibly owing to challenges of effectively distributing treatments in these areas.

CONCLUSIONS

Our findings emphasize the importance of monitoring for early detection of nascent gypsy moth colonies in order to successfully slow the spread of the invasion. Additionally, operations research should address best practices for effectively treating with patchy and open forest canopies. © 2021 Society of Chemical Industry.

Genetic diversity and relatedness in aquaculture and marina populations of the invasive tunicate Didemnum vexillum in the British Isles

Introductions of invasive, non-native species in the marine environment are increasing as human activity within coastal areas rises. Genetic datasets are useful tools to identify source populations, track routes of invasions, and illuminate the role of genetic variation in the establishment and subsequent spread of novel introductions. Here, a microsatellite dataset is used to estimate the genetic diversity and population structure of 7 introduced Didemnum vexillum populations in Britain and Ireland, 4 of which are associated with aquaculture and 3 with marinas. Genetic differentiation observed between these populations indicates human-mediated transport as the main mechanism underlying the population structure of D. vexillum in Britain and Ireland. In addition to elucidating patterns of population structure we found that aquaculture sites showed significantly higher genetic diversity (measured as allelic richness) in comparison to the marina sites. We discuss these findings in relation to the history of each invasion, the complex life history of D. vexillum, and available evidence of the relative invasiveness of these populations. Our results show numerous interesting patterns which highlight further research avenues to elucidate the complex factors underlying the global spread of this successful invader.

Potential Differences and Methods of Determining Gypsy Moth Female Flight Capabilities: Implications for the Establishment and Spread in Novel Habitats

The introduction of the Asian gypsy moth into novel environments continues with frequent interceptions in North America. There is a concern that these subspecies will pose a greater threat to the forests and urban environments of North America than the established gypsy moths (Lymantria dispardispar L.), due to their greater capacity for female flight. Asian gypsy moth populations vary in many key traits, including female flight capabilities. The potential impacts of female flight, in combination with the other key traits, on the ecology and spread of this insect are first discussed in this communication. This also provides the first review of most of the current literature on the variations in flight capability and flight distance of gypsy moth populations, as well as variation in other traits of concern and the potential methods of identification, with special attention paid to the Asian subspecies Lymantria dispar japonica Motschulsky and Lymantria dispar asiatica Vinkovskij. There are currently good tools for identifying the general origin of introduced gypsy moth populations, but these do not provide enough information to effectively manage introductions. Gypsy moth key traits differ among populations, even within each subspecies of the gypsy moth, so introduction of gypsy moths from other world areas into locations where the gypsy moth is already present could result in unwanted changes in gypsy moth biology. It also appears that the introduction of flight-capable females could enhance a population’s dispersal capability and require modifications to management protocols used for flightless females. Therefore, rapid tools to assess key traits in introduced populations are needed to adequately plan for, or deal with, new introductions into novel habitats.

Population spatial synchrony enhanced by periodicity and low detuning with environmental forcing

Explaining why fluctuations in abundances of spatially disjunct populations often are correlated through time is a major goal of population ecologists. We address two hypotheses receiving little to no testing in wild populations: (i) that population cycling facilitates synchronization given weak coupling among populations, and (ii) that the ability of periodic external forces to synchronize oscillating populations is a function of the mismatch in timescales (detuning) between the force and the population. Here, we apply new analytical methods to field survey data on gypsy moth outbreaks. We report that at timescales associated with gypsy moth outbreaks, spatial synchrony increased with population periodicity via phase locking. The extent to which synchrony in temperature and precipitation influenced population synchrony was associated with the degree of mismatch in dominant timescales of oscillation. Our study provides new empirical methods and rare empirical evidence that population cycling and low detuning can promote population spatial synchrony.

Iterative Models for Early Detection of Invasive Species across Spread Pathways

Species distribution models can be used to direct early detection of invasive species, if they include proxies for invasion pathways. Due to the dynamic nature of invasion, these models violate assumptions of stationarity across space and time. To compensate for issues of stationarity, we iteratively update regionalized species distribution models annually for European gypsy moth (Lymantria dispar dispar) to target early detection surveys for the USDA APHIS gypsy moth program. We defined regions based on the distances from the invasion spread front where shifts in variable importance occurred and included models for the non-quarantine portion of the state of Maine, a short-range region, an intermediate region, and a long-range region. We considered variables that represented potential gypsy moth movement pathways within each region, including transportation networks, recreational activities, urban characteristics, and household movement data originating from gypsy moth infested areas (U.S. Postal Service address forwarding data). We updated the models annually, linked the models to an early detection survey design, and validated the models for the following year using predicted risk at new positive detection locations. Human-assisted pathways data, such as address forwarding, became increasingly important predictors of gypsy moth detection in the intermediate-range geographic model as more predictor data accumulated over time (relative importance = 5.9%, 17.36%, and 35.76% for 2015, 2016, and 2018, respectively). Receiver operating curves showed increasing performance for iterative annual models (area under the curve (AUC) = 0.63, 0.76, and 0.84 for 2014, 2015, and 2016 models, respectively), and boxplots of predicted risk each year showed increasing accuracy and precision of following year positive detection locations. The inclusion of human-assisted pathway predictors combined with the strategy of iterative modeling brings significant advantages to targeting early detection of invasive species. We present the first published example of iterative species distribution modeling for invasive species in an operational context.

Limited Genetic Structure of Gypsy Moth Populations Reflecting a Recent History in Europe

The gypsy moth, Lymantria dispar, a prominent polyphagous species native to Eurasia, causes severe impacts in deciduous forests during irregular periodical outbreaks. This study aimed to describe the genetic structure and diversity among European gypsy moth populations. Analysis of about 500 individuals using a partial region of the mitochondrial COI gene, L. dispar was characterized by low genetic diversity, limited population structure, and strong evidence that all extant haplogroups arose via a single Holocene population expansion event. Overall 60 haplotypes connected to a single parsimony network were detected and genetic diversity was highest for the coastal populations Croatia, Italy, and France, while lowest in continental populations. Phylogenetic reconstruction resulted in three groups that were geographically located in Central Europe, Dinaric Alps, and the Balkan Peninsula. In addition to recent events, the genetic structure reflects strong gene flow and the ability of gypsy moth to feed on about 400 deciduous and conifer species. Distinct genetic groups were detected in populations from Georgia. This remote population exhibited haplotypes intermediate to the European L. dispar dispar, Asian L. dispar asiatica, and L. dispar japonica clusters, highlighting this area as a possible hybridization zone of this species for future studies applying genomic approaches.

Genetic signals of artificial and natural dispersal linked to colonization of South America by non-native Chinook salmon (Oncorhynchus tshawytscha)

Genetics data have provided unprecedented insights into evolutionary aspects of colonization by non-native populations. Yet, our understanding of how artificial (human-mediated) and natural dispersal pathways of non-native individuals influence genetic metrics, evolution of genetic structure, and admixture remains elusive. We capitalize on the widespread colonization of Chinook salmon Oncorhynchus tshawytscha in South America, mediated by both dispersal pathways, to address these issues using data from a panel of polymorphic SNPs. First, genetic diversity and the number of effective breeders (Nb) were higher among artificial than natural populations. Contemporary gene flow was common between adjacent artificial and natural and adjacent natural populations, but uncommon between geographically distant populations. Second, genetic structure revealed four distinct clusters throughout the Chinook salmon distributional range with varying levels of genetic connectivity. Isolation by distance resulted from weak differentiation between adjacent artificial and natural and between natural populations, with strong differentiation between distant Pacific Ocean and Atlantic Ocean populations, which experienced strong genetic drift. Third, genetic mixture analyses revealed the presence of at least six donor geographic regions from North America, some of which likely hybridized as a result of multiple introductions. Relative propagule pressure or the proportion of Chinook salmon propagules introduced from various geographic regions according to government records significantly influenced genetic mixtures for two of three artificial populations. Our findings support a model of colonization in which high-diversity artificial populations established first; some of these populations exhibited significant admixture resulting from propagule pressure. Low-diversity natural populations were likely subsequently founded from a reduced number of individuals.

Landscape-Level Patterns of Elevated FS1 Asian Allele Frequencies in Populations of Gypsy Moth (Lepidoptera: Erebidae) at a Northern U.S. Boundary

From a regulatory perspective, Asian gypsy moth is a species complex consisting of three species of Lymantria and two subspecies of Lymantria dispar (L.), differing from the European subspecies, L. dispar dispar (L.), by having consistently flight-capable females. As such, the invasion potential in North America is thought to exceed that of European gypsy moth. USDA-APHIS therefore has a monitoring program to detect Asian gypsy moth at high-risk introduction pathways. Molecular markers are used to improve the diagnosis of Asian gypsy moth. One such marker, which targets the FS1 locus, detects an allele, FS1-A, prevalent in Asian populations but occurring at low frequencies (3-6%) throughout the European gypsy moth's range in North America. However, some locales, such as Minnesota, exhibit elevated FS1-A frequencies. We studied the distribution of the FS1-A allele in northern Minnesota, 2013-2014, assessing spatial patterns in the distribution of the FS1-A allele using Moran's I and using spatial regression techniques to examine if the FS1-A allele was associated with putative movement pathways. We also used time series analysis to discern if temporal patterns in FS1-A or possible introduction events occurred. Our results indicated that FS1-A occurred randomly in space and time. We found no evidence that elevated FS1-A frequencies were associated with movement pathways or possible immigration events into this region over the two years. Elevated frequencies of the FS1-A allele within this region could be due to genetic drift and allelic surfing along the expanding population front, or to selection of physiological or behavioral traits.

Foliage maturity of Quercus ilex affects the larval development of a Croatian coastal population of Lymantria dispar (Lepidoptera: Erebidae)

Gypsy moth (Lymantria dispar) is one of the most important forest pests in the world. Numerous previous studies focused only on different host tree species, but small number of them on foliage age. Since recent genetic analyses showed that there are significant differences between Croatian continental and coastal population we investigated coastal population since there was no previous research. For this research juvenile and mature foliage of Holm oak (Quercus ilex) was used. Larval development was investigated in two laboratory experiments. One experiment involved rearing trials consisting of 50 individual larvae per treatment while a parallel experiment investigated gregarious feeding conditions using 120 larvae in a rearing treatment. Larval mortality was lower and development time shorter for individuals reared on juvenile foliage. High pupation success in both the individual and group rearing experiment for larvae reared on juvenile foliage was also observed. To conclude, our results showed high mortality, poor larval development and low pupation success in larvae reared on mature foliage. This research is significant because in the aspect of ongoing climate changes there is a possibility that gypsy moth will move to the north and shift its distribution by expanding into new climatic area.

Simulation Modeling to Interpret the Captures of Moths in Pheromone-Baited Traps Used for Surveillance of Invasive Species: the Gypsy Moth as a Model Case

When pheromone traps are used for detection of an invasive pest and then delimitation of its distribution, an unresolved issue is the interpretation of failure to capture any target insects. Is a population present but not detected, a so-called false negative? Using the gypsy moth (Lymantria dispar) as an exemplar, we modeled the probability of males being captured in traps deployed at densities typical for surveillance (1 per 2.6 km² or 1 per mi²) and delimitation (up to 49 per 2.6 km²). The simulations used a dynamic wind model generating a turbulent plume structure and varying wind direction, and a behavior model based on the documented maneuvers of gypsy moths during plume acquisition and along-plume navigation. Several strategies of plume acquisition using Correlated Random Walks were compared to ensure that the generated dispersions over three days were not either overly clumped or ranged many km. Virtual moths were released into virtual space with patterns mimicking prior releases of gypsy moth males in Massachusetts at varying distance from a baited trap. In general, capture rates of virtual and real moths at varying trap densities were similar. One application of this approach was to estimate through bootstrapping the probabilities of not detecting populations having densities ranging from 1 to 100 moths per 2.6 km² and using traps that varied from 25 to 100 % in their efficiencies of capture. Low-level populations (e.g., 20-30 per 2.6 km²) often were not detected with one trap per 2.6 km², especially when traps had low efficiencies.

Multiple dispersal vectors drive range expansion in an invasive marine species

The establishment and subsequent spread of invasive species is widely recognized as one of the most threatening processes contributing to global biodiversity loss. This is especially true for marine and estuarine ecosystems, which have experienced significant increases in the number of invasive species with the increase in global maritime trade. Understanding the rate and mechanisms of range expansion is therefore of significant interest to ecologists and conservation managers alike. Using a combination of population genetic surveys, environmental DNA (eDNA) plankton sampling and hydrodynamic modelling, we examined the patterns of introduction of the predatory Northern Pacific seastar (Asterias amurensis) and pathways of secondary spread within southeast Australia. Genetic surveys across the invasive range reveal some genetic divergence between the two main invasive regions and no evidence of ongoing gene flow, a pattern that is consistent with the establishment of the second invasive region via a human-mediated translocation event. In contrast, hydrodynamic modelling combined with eDNA plankton sampling demonstrated that the establishment of range expansion populations within a region is consistent with natural larval dispersal and recruitment. Our results suggest that both anthropogenic and natural dispersal vectors have played an important role in the range expansion of this species in Australia. The multiple modes of spread combined with high levels of fecundity and a long larval duration in A. amurensis suggests it is likely to continue its range expansion and significantly impact Australian marine ecosystems.

Interpretation of gypsy moth frontal advance using meteorology in a conditional algorithm

The gypsy moth, Lymantria dispar, is a non-native species that continues to invade areas in North America. It spreads generally through stratified dispersal where local growth and diffusive spread are coupled with long-distance jumps ahead of the leading edge. Long-distance jumps due to anthropogenic movement of life stages is a well-documented spread mechanism. Another mechanism is the atmospheric transport of early instars and adult males, believed to occur over short distances. However, empirical gypsy moth population data continue to support the possibility of alternative methods of long-range dispersal. Such dispersal events seemed to have occurred in the mid- to late-1990s with spread across Lake Michigan to Wisconsin. Such dispersal would be against the prevailing wind flow for the area and would have crossed a significant physical barrier (Lake Michigan). The climatology of the region shows that vigorous cyclones can result in strong easterly winds in the area at the time when early instars are present. It is hypothesized that these storms would enable individuals to be blown across the Lake and explain the appearance of new population centers observed at several locations on the western shore of Lake Michigan nearly simultaneously. A synoptic climatology model coupled with population dynamics data from the area was parameterized to show an association between transport events and population spread from 1996 to 2007. This work highlights the importance of atmospheric transport events relative to the invasion dynamics of the gypsy moth, and serves as a model for understanding this mechanism of spread in other related biological invasions.

Dispersal behavior of neonate European corn borer (Lepidoptera: Crambidae) on Bt corn

European corn borer, Ostrinia nubilalis (Hübner) (Lepidoptera: Crambidae), has historically been a significant economically important insect pest of corn (Zea mays L.) in the United States and Canada. The development in the 1990s of genetically modified corn expressing genes derived from Bacillus thuringiensis (Bt) that encodes insecticidal crystalline (Cry) proteins has proven to be effective in controlling this insect as well as other corn pests. The purpose of this study was to assess the movement and dispersal behavior of neonate European corn borer on Bt corn. We examined differences in neonate European corn borer dispersal behavior for the first 4 h after eclosion in the field among a stacked pyramid (Cry1F X Cry1Ab X Cry34/35Ab1) Bt corn, a Cry1F Bt corn, and a non-Bt sweet corn; and in the laboratory among a Bt corn hybrid containing Cry1F, a hybrid containing Cry1Ab, a pyramid combining these two hybrids (Cry1F X Cry1Ab), and a non-Bt near isoline corn. In field experiments, we found that dispersal was significantly higher on Bt corn compared with sweet corn. In laboratory experiments, dispersal was significantly higher on Cry1Ab Bt corn and Cry1F X Cry1Ab Bt corn than on non-Bt near isoline corn. Results indicated that neonate dispersal may be significantly greater in Bt cornfields compared with non-Bt cornfields. The findings on dispersal behavior in this study will be useful in evaluating the efficacy of a blended seed refuge system for managing European corn borer resistance in Bt corn.

Dispersal strategies, secondary range expansion and invasion genetics of the nonindigenous round goby, Neogobius melanostomus, in Great Lakes tributaries

Dispersal strategies are important mechanisms underlying the spatial distribution and colonizing ability of all mobile species. In the current study, we use highly polymorphic microsatellite markers to evaluate local dispersal and colonization dynamics of the round goby (Neogobius melanostomus), an aquatic invader expanding its range from lake to river environments in its introduced North American range. Genetic structure, genotype assignment and genetic diversity were compared among 1262 round gobies from 20 river and four lake sites in three Great Lakes tributaries. Our results indicate that a combination of short-distance diffusion and long-distance dispersal, collectively referred to as 'stratified dispersal', is facilitating river colonization. Colonization proceeded upstream yearly (approximately 500 m/year; 2005-2009) in one of two temporal replicates while genetic structure was temporally stable. Contiguous dispersal from the lake was observed in all three rivers with a substantial portion of river fish (7.3%) identified as migrants. Genotype assignment indicated a separate introduction occurred upstream of the invasion front in one river. Genetic diversity was similar and relatively high among lake and recently colonized river populations, indicating that founder effects are mitigated through a dual-dispersal strategy. The remarkable success of round goby as an aquatic invader stresses the need for better diffusion models of secondary range expansion for presumably sessile invasive species.

Critical patch size generated by Allee effect in gypsy moth, Lymantria dispar (L.)

Allee effects are important dynamical mechanisms in small-density populations in which per capita population growth rate increases with density. When positive density dependence is sufficiently severe (a ‘strong’ Allee effect), a critical density arises below which populations do not persist. For spatially distributed populations subject to dispersal, theory predicts that the occupied area also exhibits a critical threshold for population persistence, but this result has not been confirmed in nature. We tested this prediction in patterns of population persistence across the invasion front of the European gypsy moth (Lymantria dispar) in the United States in data collected between 1996 and 2008. Our analysis consistently provided evidence for effects of both population area and density on persistence, as predicted by the general theory, and confirmed here using a mechanistic model developed for the gypsy moth system. We believe this study to be the first empirical documentation of critical patch size induced by an Allee effect.