Round-trip across the Sahara: Afrotropical Painted Lady butterflies recolonize the Mediterranean in early spring

Biomass burning smoke pollution stimulates painted lady butterflies (Vanessa cardui L.) to increase flight speed

Smoke from biomass burning significantly degrades air quality due to high concentrations of particulate matter (PM2.5) and trace gases. While the ecological and health impacts of smoke pollution are well documented, its effects on insect migration remain poorly understood. In this study, we conducted two experiments to investigate the flight performance of Vanessa cardui butterflies under varying smoke conditions and identify the mechanisms influencing their behaviour. Butterflies were tethered to flight mills (TFMs) for 6 h, during which flight speed, distance, and duration were recorded across clean-air conditions and three levels of PM2.5 concentrations. Statistical analysis revealed that flight speed increases significantly as smoke concentration increases, although the increased range decreases. At a mean PM2.5 concentration of 120 μg m-3, flight speed increased by 52 % compared to clean-air conditions. To determine whether particulate matter was driving this response, individuals were exposed to smoke with and without particulates. In smoke with particulates retained, butterflies exhibited nearly double the flight speed compared to filtered smoke, indicating that particulates play a key role in altering flight behaviour. Scanning electron microscopy revealed significant deposition of smoke particulates on the antennae and abdomen, suggesting a sensory or physical response triggering accelerated flight. We interpret these findings as evidence that Vanessa cardui accelerates flight in smoky environments as an escape response. This study highlights the remarkable sensitivity of butterflies to smoke pollution and provides novel insights into the ecological consequences of biomass burning, particularly its potential impacts on insect behaviour and migration dynamics.

Trans-Atlantic Dispersal and Introgression Explain Holarctic Disjunct Distributions in Vanessa Butterflies

Species with disjunct distributions have long puzzled evolutionary biologists and biogeographers. Long-distance dispersal can play a pivotal role in generating intra-specific disjunct distributions, initiating early stages of allopatric speciation and leading to eventual interspecific disjunctions. Vanessa butterflies exhibit diverse movement behaviours, from low-dispersal species with restricted distributions to others that engage in annual extensive migratory cycles. The biogeographic history of Vanessa presents intriguing cases of both intra- and interspecific disjunctions. Vanessa atalanta is present in the Nearctic and Western Palearctic but is absent in Asia, while its sister species V. tameamea is endemic to Hawaii. Vanessa indica occurs only in Asia, and its sister species, V. vulcania, is endemic to Macaronesia. Here, we investigate this conundrum through population genomics and demographic analyses of Vanessa atalanta using ddRAD data from 70 samples across its entire distributional range, identifying two genetically differentiated populations separated by the Atlantic Ocean. Demographic simulations and phylogenetic analyses suggest that these originated via long-distance dispersal from the Nearctic to Europe around the Last Glacial Maximum. Hybridisation tests revealed introgression between the Palearctic population of V. atalanta and V. indica, indicating that their distributions overlapped during V. atalanta's colonisation of Europe. We hypothesise that V. atalanta caused a species displacement of V. indica from Europe to Asia, explaining their current allopatric distributions-a scenario that is supported by ecological niche modelling. Together, our results illustrate the role of long-distance dispersal and species interactions in shaping complex biogeographic patterns.

Trans-Saharan migratory patterns in Vanessa cardui and evidence for a southward leapfrog migration

Some insects, such as the painted lady butterfly Vanessa cardui, exhibit complex annual migratory cycles spanning multiple generations. Traversing extensive seas or deserts is often a required segment of these migratory journeys. We develop a bioavailable strontium isoscape for Europe and Africa and then use isotope geolocation combining hydrogen and strontium isotopes to estimate the natal origins of painted ladies captured north and south of the Sahara during spring and autumn, respectively. Our findings reveal moderate migratory connectivity across the Sahara characterized by a broad-front, parallel migration. We also report evidence of a leapfrog migration, wherein early autumn migrants from higher latitudes cover greater distances southward than their late autumn counterparts. This work represents a major advancement in understanding insect migratory patterns and connectivity, particularly across extensive barriers, which is essential for understanding population dynamics and predicting the impacts of global change on insect-mediated ecosystem services.

Strontium isotope and trajectory method elucidating overseas migration of Mythimna separata to Japan

The oriental armyworm, Mythimna separata, generally migrates from eastern to northeastern China in early summer, and some individuals are believed to migrate overseas to Japan depending on meteorological conditions. This potential migration was investigated with the immigrants' strontium radiogenic isotope ratio 87Sr/86Sr and backward flight trajectories from Japanese trapping sites. The results showed that the 87Sr/86Sr ratios of Chinese reared M. separata were significantly higher than those of reared insects of Japanese immigration areas. As some individuals trapped in western Japan had 87Sr/86Sr ratios higher than a statistical discriminating ratio, they likely originated in China. Trajectory analysis also indicated those individuals might have originated from the East Asian continent, such as the first-generation outbreak region in China and their migration waypoint regions. Our analysis, thus, suggests direct or multistep overseas migration of individual M. separata from the East Asian continent to Japan, giving insight into migration pathways and population dynamics.

A trans-oceanic flight of over 4,200 km by painted lady butterflies

The extent of aerial flows of insects circulating around the planet and their impact on ecosystems and biogeography remain enigmatic because of methodological challenges. Here we report a transatlantic crossing by Vanessa cardui butterflies spanning at least 4200 km, from West Africa to South America (French Guiana) and lasting between 5 and 8 days. Even more, we infer a likely natal origin for these individuals in Western Europe, and the journey Europe-Africa-South America could expand to 7000 km or more. This discovery was possible through an integrative approach, including coastal field surveys, wind trajectory modelling, genomics, pollen metabarcoding, ecological niche modelling, and multi-isotope geolocation of natal origins. The overall journey, which was energetically feasible only if assisted by winds, is among the longest documented for individual insects, and potentially the first verified transatlantic crossing. Our findings suggest that we may be underestimating transoceanic dispersal in insects and highlight the importance of aerial highways connecting continents by trade winds.

Pollen metabarcoding reveals the origin and multigenerational migratory pathway of an intercontinental-scale butterfly outbreak

Migratory insects may move in large numbers, even surpassing migratory vertebrates in biomass. Long-distance migratory insects complete annual cycles through multiple generations, with each generation's reproductive success linked to the resources available at different breeding grounds. Climatic anomalies in these grounds are presumed to trigger rapid population outbreaks. Here, we infer the origin and track the multigenerational path of a remarkable outbreak of painted lady (Vanessa cardui) butterflies that took place at an intercontinental scale in Europe, the Middle East, and Africa from March 2019 to November 2019. Using metabarcoding, we identified pollen transported by 264 butterflies captured in 10 countries over 7 months and modeled the distribution of the 398 plants detected. The analysis showed that swarms collected in Eastern Europe in early spring originated in Arabia and the Middle East, coinciding with a positive anomaly in vegetation growth in the region from November 2018 to April 2019. From there, the swarms advanced to Northern Europe during late spring, followed by an early reversal toward southwestern Europe in summer. The pollen-based evidence matched spatiotemporal abundance peaks revealed by citizen science, which also suggested an echo effect of the outbreak in West Africa during September-November. Our results show that population outbreaks in a part of species' migratory ranges may disseminate demographic effects across multiple generations in a wide geographic area. This study represents an unprecedented effort to track a continuous multigenerational insect migration on an intercontinental scale.

A hydrogen isoscape for tracing the migration of herbivorous lepidopterans across the Afro-Palearctic range

RATIONALE

Many insect species undertake multigenerational migrations in the Afro-tropical and Palearctic ranges, and understanding their migratory connectivity remains challenging due to their small size, short life span and large population sizes. Hydrogen isotopes (δ2 H) can be used to reconstruct the movement of dispersing or migrating insects, but applying δ2 H for provenance requires a robust isotope baseline map (i.e. isoscape) for the Afro-Palearctic.

METHODS

We analyzed the δ2 H in the wings (δ2 Hwing ) of 142 resident butterflies from 56 sites across the Afro-Palearctic. The δ2 Hwing values were compared to the predicted local growing-season precipitation δ2 H values (δ2 HGSP ) using a linear regression model to develop an insect wing δ2 H isoscape. We used multivariate linear mixed models and high-resolution and time-specific remote sensing climate and environmental data to explore the controls of the residual δ2 Hwing variability.

RESULTS

A strong linear relationship was found between δ2 Hwing and δ2 HGSP values (r2  = 0.53). The resulting isoscape showed strong patterns across the Palearctic but limited variation and high uncertainty for the Afro-tropics. Positive residuals of this relationship were correlated with dry conditions for the month preceding sampling whereas negative residuals were correlated with more wet days for the month preceding sampling. High intra-site δ2 Hwing variance was associated with lower relative humidity for the month preceding sampling and higher elevation.

CONCLUSION

The δ2 Hwing isoscape is applicable for tracing herbivorous lepidopteran insects that migrate across the Afro-Palearctic range but has limited geolocation potential in the Afro-tropics. The spatial analysis of uncertainty using high-resolution climatic data demonstrated that many African regions with highly variable evaporation rates and relative humidity have δ2 Hwing values that are less related to δ2 HGSP values. Increasing geolocation precision will require new modeling approaches using more time-specific environmental data and/or independent geolocation tools.

The frequency of wing damage in a migrating butterfly

The ability to fly is crucial for migratory insects. Consequently, the accumulation of damage on the wings over time can affect survival, especially for species that travel long distances. We examined the frequency of irreversible wing damage in the migratory butterfly Vanessa cardui to explore the effect of wing structure on wing damage frequency, as well as the mechanisms that might mitigate wing damage. An exceptionally high migration rate driven by high precipitation levels in their larval habitats in the winter of 2018-2019 provided us with an excellent opportunity to collect data on the frequency of naturally occurring wing damage associated with long-distance flights. Digital images of 135 individuals of V. cardui were collected and analyzed in Germany. The results show that the hindwings experienced a greater frequency of damage than the forewings. Moreover, forewings experienced more severe damage on the lateral margin, whereas hindwings experienced more damage on the trailing margin. The frequency of wing margin damage was higher in the painted lady butterfly than in the migrating monarch butterfly and in the butterfly Pontia occidentalis following artificially induced wing collisions. The results of this study could be used in future comparative studies of patterns of wing damage in butterflies and other insects. Additional studies are needed to clarify whether the strategies for coping with wing damage differ between migratory and nonmigratory species.

Haemolymph viscosity in hawkmoths and its implications for hovering flight

Viscosity determines the resistance of haemolymph flow through the insect body. For flying insects, viscosity is a major physiological parameter limiting flight performance by controlling the flow rate of fuel to the flight muscles, circulating nutrients and rapidly removing metabolic waste products. The more viscous the haemolymph, the greater the metabolic energy needed to pump it through confined spaces. By employing magnetic rotational spectroscopy with nickel nanorods, we showed that viscosity of haemolymph in resting hawkmoths (Sphingidae) depends on wing size non-monotonically. Viscosity increases for small hawkmoths with high wingbeat frequencies, reaches a maximum for middle-sized hawkmoths with moderate wingbeat frequencies, and decreases in large hawkmoths with slower wingbeat frequencies but greater lift. Accordingly, hawkmoths with small and large wings have viscosities approaching that of water, whereas hawkmoths with mid-sized wings have more than twofold greater viscosity. The metabolic demands of flight correlate with significant changes in circulatory strategies via modulation of haemolymph viscosity. Thus, the evolution of hovering flight would require fine-tuned viscosity adjustments to balance the need for the haemolymph to carry more fuel to the flight muscles while decreasing the viscous dissipation associated with its circulation.

Environmental stress during larval development induces DNA methylation shifts in the migratory painted lady butterfly (Vanessa cardui)

Seasonal environmental fluctuations provide formidable challenges for living organisms, especially small ectotherms such as butterflies. A common strategy to cope with harsh environments is to enter diapause, but some species avoid unsuitable conditions by migrating. Despite a growing understanding of migration in the life cycles of some butterfly species, it remains unknown how individuals register and store environmental cues to determine whether and where to migrate. Here, we explored how competition and host plant availability during larval development affect patterns of DNA methylation in the migratory painted lady (Vanessa cardui) butterfly. We identify a set of potentially functional methylome shifts associated with differences in the environment, indicating that DNA methylation is involved in the response to different conditions during larval development. By analysing the transcriptome for the same samples used for methylation profiling, we also uncovered a non-monotonic relationship between gene body methylation and gene expression. Our results provide a starting point for understanding the interplay between DNA methylation and gene expression in butterflies in general and how differences in environmental conditions during development can trigger unique epigenetic marks that might be important for behavioural decisions in the adult stage.

The Afrotropical breeding grounds of the Palearctic-African migratory painted lady butterflies (Vanessa cardui)

Migratory insects are key players in ecosystem functioning and services, but their spatiotemporal distributions are typically poorly known. Ecological niche modeling (ENM) may be used to predict species seasonal distributions, but the resulting hypotheses should eventually be validated by field data. The painted lady butterfly (Vanessa cardui) performs multigenerational migrations between Europe and Africa and has become a model species for insect movement ecology. While the annual migration cycle of this species is well understood for Europe and northernmost Africa, it is still unknown where most individuals spend the winter. Through ENM, we previously predicted suitable breeding grounds in the subhumid regions near the tropics between November and February. In this work, we assess the suitability of these predictions through i) extensive field surveys and ii) two-year monitoring in six countries: a large-scale monitoring scheme to study butterfly migration in Africa. We document new breeding locations, year-round phenological information, and hostplant use. Field observations were nearly always predicted with high probability by the previous ENM, and monitoring demonstrated the influence of the precipitation seasonality regime on migratory phenology. Using the updated dataset, we built a refined ENM for the Palearctic-African range of V. cardui. We confirm the relevance of the Afrotropical region and document the missing natural history pieces of the longest migratory cycle described in butterflies.

Temporal stability of δ2H in insect tissues: Implications for isotope-based geographic assignments

Hydrogen isotope geolocation of insects is based on the assumption that the chitin in the wings of adult migratory insects preserves the hydrogen isotope composition (δ2H) of the larval stages without influence of adult diet. Here, we test this assumption by conducting laboratory feeding experiments for monarch butterflies (Danaus plexippus) including: (1) a starvation treatment where adults were not fed and (2) an enriched treatment where adults were fed a diet isotopically enriched in deuterium (~ +78‰) compared to the larval diet. The δ2H values of adult wings were measured at different time steps along the 24-day experiment. We also investigated intra-wing differences in δ2H values caused by wing pigmentation, absence of wing scales, and presence of major wing veins. We conclude that, although the magnitude of the changes in δ2H values are small (~6‰), wing δ2H values vary based on adult diet and insect age, particularly early after eclosion (i.e., 1–4 days). We found that wing shade, wing pigmentation, and the presence of wing scales do not alter wing δ2H values. However, wing samples containing veins had systematically higher δ2H values (~9‰), suggesting that adult diet influences the hemolymph that circulates in the wing veins. We hypothesise that there is a stronger influence of adult diet on the isotope signal of wings during early adult life relative to later life because of increased metabolic and physiologic activity in young insect wings. We argue that the influence of the isotopic contribution of adult diet is generally small and is likely minimal if the wings are carefully sampled to avoid veins. However, we also demonstrated that wings are not inert tissues, and that adult feeding contributes to some of the intra-population δ2H variance. We conclude that δ2H geolocation using insect wings remains valid, but that adult feeding, butterfly age and wing vein sampling generate an inherent uncertainty limiting the precision of geolocation.

Migratory behaviour is positively associated with genetic diversity in butterflies

Migration is typically associated with risk and uncertainty at the population level, but little is known about its cost-benefit trade-offs at the species level. Migratory insects in particular often exhibit strong demographic fluctuations due to local bottlenecks and outbreaks. Here, we use genomic data to investigate levels of heterozygosity and long-term population size dynamics in migratory insects, as an alternative to classical local and short-term approaches such as regional field monitoring. We analyse whole-genome sequences from 97 Lepidoptera species and show that individuals of migratory species have significantly higher levels of genome-wide heterozygosity, a proxy for effective population size, than do nonmigratory species. Also, we contribute whole-genome data for one of the most emblematic insect migratory species, the painted lady butterfly (Vanessa cardui), sampled across its worldwide distributional range. This species exhibits one of the highest levels of genomic heterozygosity described in Lepidoptera (2.95 ± 0.15%). Coalescent modelling (PSMC) shows historical demographic stability in V. cardui, and high effective population size estimates of 2-20 million individuals 10,000 years ago. The study reveals that the high risks associated with migration and local environmental fluctuations do not seem to decrease overall genetic diversity and demographic stability in migratory Lepidoptera. We propose a "compensatory" demographic model for migratory r-strategist organisms in which local bottlenecks are counterbalanced by reproductive success elsewhere within their typically large distributional ranges. Our findings highlight that the boundaries of populations are substantially different for sedentary and migratory insects, and that, in the latter, local and even regional field monitoring results may not reflect whole population dynamics. Genomic diversity patterns may elucidate key aspects of an insect's migratory nature and population dynamics at large spatiotemporal scales.

Community Science Reveals High Diversity of Nectaring Plants Visited by Painted Lady Butterflies (Lepidoptera: Nymphalidae) in California Sage Scrub

California's sage scrub habitats support a diversity of nectar and host plants for migrating and resident populations of painted lady butterflies (Vanessa cardui) throughout all seasons. North America experiences spring V. cardui migrations involving butterflies totaling in the millions in some years. These irruptive years are thought to be driven by winter weather patterns at breeding grounds near the US-Mexico border and due to their irregularity, it is difficult to study floral resource use along the migration route. Here we used the community science platform iNaturalist to quantify patterns in V. cardui nectar resource use in sage scrub over time and space during irruptive and nonirruptive years. We identified over 329 different nectaring plant species of varying functional types (72% native to California) visited by adult V. cardui, 195 of which had not been previously identified as known nectar plants for V. cardui. Vanessa cardui butterflies were observed in similar locations regardless of whether an irruptive migration occurred, indicating the presence of either sparse migrants or resident populations across California. Moreover, irruptive years were positively correlated with warmer and wetter local conditions at observation locations. Our results provide new insights into patterns of floral resource use by North American V. cardui by harnessing the power of community science data and while highlighting the factors associated with its North American migration.

Linkage mapping and genome annotation give novel insights into gene family expansions and regional recombination rate variation in the painted lady (Vanessa cardui) butterfly

Characterization of gene family expansions and crossing over is crucial for understanding how organisms adapt to the environment. Here, we develop a high-density linkage map and detailed genome annotation of the painted lady butterfly (Vanessa cardui) - a non-diapausing, highly polyphagous species famous for its long-distance migratory behavior and almost cosmopolitan distribution. Our results reveal a complex interplay between regional recombination rate variation, gene duplications and transposable element activity shaping the genome structure of the painted lady. We identify several lineage specific gene family expansions. Their functions are mainly associated with protein and fat metabolism, detoxification, and defense against infection - critical processes for the painted lady's unique life-history. Furthermore, the detailed recombination maps allow us to characterize the regional recombination landscape, data that reveal a strong effect of chromosome size on the recombination rate, a limited impact of GC-biased gene conversion and a positive association between recombination and short interspersed elements.

The genome sequence of the painted lady, Vanessa cardui Linnaeus 1758

We present a genome assembly from an individual female Vanessa cardui (the painted lady; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 425 megabases in span. The majority of the assembly is scaffolded into 32 chromosomal pseudomolecules, with the W and Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 12,821 protein coding genes.

Unraveling the World’s Longest Non-stop Migration: The Indian Ocean Crossing of the Globe Skimmer Dragonfly

Insect migration redistributes enormous quantities of biomass, nutrients and species globally. A subset of insect migrants perform extreme long-distance journeys, requiring specialized morphological, physiological and behavioral adaptations. The migratory globe skimmer dragonfly (Pantala flavescens) is hypothesized to migrate from India across the Indian Ocean to East Africa in the autumn, with a subsequent generation thought to return to India from East Africa the following spring. Using an energetic flight model and wind trajectory analysis, we evaluate the dynamics of this proposed transoceanic migration, which is considered to be the longest regular non-stop migratory flight when accounting for body size. The energetic flight model suggests that a mixed strategy of gliding and active flapping would allow a globe skimmer to stay airborne for up to 230–286 h, assuming that the metabolic rate of gliding flight is close to that of resting. If engaged in continuous active flapping flight only, the flight time is severely reduced to ∼4 h. Relying only on self-powered flight (combining active flapping and gliding), a globe skimmer could cross the Indian Ocean, but the migration would have to occur where the ocean crossing is shortest, at an exceptionally fast gliding speed and with little headwind. Consequently, we deem this scenario unlikely and suggest that wind assistance is essential for the crossing. The wind trajectory analysis reveals intra- and inter-seasonal differences in availability of favorable tailwinds, with only 15.2% of simulated migration trajectories successfully reaching land in autumn but 40.9% in spring, taking on average 127 and 55 h respectively. Thus, there is a pronounced requirement on dragonflies to be able to select favorable winds, especially in autumn. In conclusion, a multi-generational, migratory circuit of the Indian Ocean by the globe skimmer is shown to be achievable, provided that advanced adaptations in physiological endurance, behavior and wind selection ability are present. Given that migration over the Indian Ocean would be heavily dependent on the assistance of favorable winds, occurring during a relatively narrow time window, the proposed flyway is potentially susceptible to disruption, if wind system patterns were to be affected by climatic change.

An insect isoscape of UK and Ireland

RATIONALE

The study of insect migration is problematic due to the small size of insects. Stable isotope analysis can be used to elucidate movement, either by geographic assignment of location of a species, or by simply distinguishing migrant from resident populations. There are few isoscapes of any kind in the UK/Ireland available for interrogation. Thus, I have measured stable isotope ratios (of H, C, N and S) of 299 individuals of the non-migratory Brimstone moth (Opisthograptis luteolata) collected from 93 locations around the UK and Ireland by citizen scientists.

METHODS

After removing lipids, stable isotope ratios were measured by continuous flow isotope ratio mass spectrometry, using either a conventional elemental analyser (C, N and S) or a high-temperature, thermal conversion elemental analyser in reductive mode.

RESULTS

Maps (isoscapes) were constructed that illustrate the stable isotope spatial distribution of this insect. These are the first isoscapes of H, C, N and S of biological samples covering both UK and Ireland.

CONCLUSIONS

The insect isoscape patterns can be explained from what we know of moth diet, climate and geology. Sulfur isotopes may be of particular use for distinguishing individuals from areas of unique geology. Isoscape patterns may (with care) predict isotope compositions of other, herbivorous, non-aquatic, chitinous taxa. Such isoscapes, when extended beyond the UK and Ireland, would provide a useful tool to elucidate insect migration.

Environmental drivers of annual population fluctuations in a trans-Saharan insect migrant

The painted lady butterfly is an annual migrant to northern regions, but the size of the immigration varies by more than 100-fold in successive years. Unlike the monarch, the painted lady breeds year round, and it has long been suspected that plant-growing conditions in winter-breeding locations drive this high annual variability. However, the regions where caterpillars develop over winter remained unclear. Here, we show for the European summer population that winter plant greenness in the savanna of sub-Saharan Africa is the key driver of the size of the spring immigration. Our results show that painted ladies regularly cross the Sahara Desert and elucidate the climatic drivers of the annual population dynamics.

Many latitudinal insect migrants including agricultural pests, disease vectors, and beneficial species show huge fluctuations in the year-to-year abundance of spring immigrants reaching temperate zones. It is widely believed that this variation is driven by climatic conditions in the winter-breeding regions, but evidence is lacking. We identified the environmental drivers of the annual population dynamics of a cosmopolitan migrant butterfly (the painted lady Vanessa cardui) using a combination of long-term monitoring and climate and atmospheric data within the western part of its Afro-Palearctic migratory range. Our population models show that a combination of high winter NDVI (normalized difference vegetation index) in the Savanna/Sahel of sub-Saharan Africa, high spring NDVI in the Maghreb of North Africa, and frequent favorably directed tailwinds during migration periods are the three most important drivers of the size of the immigration to western Europe, while our atmospheric trajectory simulations demonstrate regular opportunities for wind-borne trans-Saharan movements. The effects of sub-Saharan vegetative productivity and wind conditions confirm that painted lady populations on either side of the Sahara are linked by regular mass migrations, making this the longest annual insect migration circuit so far known. Our results provide a quantification of the environmental drivers of large annual population fluctuations of an insect migrant and hold much promise for predicting invasions of migrant insect pests, disease vectors, and beneficial species.

Seasonal spatial dynamics of butterfly migration

Understanding the seasonal movements of migratory species underpins ecological studies. Several hundred butterfly species show migratory behaviour, yet the spatial pattern of these migrations is poorly understood. We developed climatic niche models for 405 migratory butterfly species globally to estimate patterns of seasonal movement and the distribution of seasonal habitat suitability. We found strong seasonal variation in habitat suitability for most migratory butterflies with >75% of pixels within their distributions showing seasonal switching in predicted occupancy for 85% of species. The greatest rate of seasonal switching occurred in the tropics. Several species showed extreme range fluctuations between seasons, exceeding 10-fold for 53 species (13%) and more than 100-fold for nine species (2%), suggesting that such species may be at elevated extinction risk. Our results can be used to search for the ecological processes that underpin migration in insects, as well as to design conservation interventions for declining migratory insects.

Long-term temporal patterns in flight activities of a migrant diurnal butterfly

Recent studies demonstrated that the Painted Lady (Vanessa cardui), a cosmopolitan diurnal butterfly performs long-range migration between subtropical Africa and north-western Europe, covered by individuals belonging to up to six generations. Here we analyze temporal patterns of complete annual migratory activity of the Painted Lady in Hungary, located in its Central European migratory route, almost completely unstudied before. To do so, we used field occurrence data collected between 2000 and 2019 and estimated temporal patterns in migratory activity by fitting kernel density functions on the daily mean number of individuals and observation frequency. The temporal distributions of kernel density estimates were analyzed as a function of time and key climatic predictors of the study area. We found that (i) the timing of spring arrivals has been advancing; (ii) the relative intensity of the first and last migratory peaks of the Painted Lady significantly increased during the past decades; and (iii) intensity of the last migratory peak is related to the mean temperature of the previous month, inferring that the migration is shifting to earlier dates and their volume of the migration has substantially intensified, evoking mutually nonexclusive, competing hypotheses. Our study indicates the strengthening migration activities of a southerly distributed, long-distance migrant diurnal butterfly, most probably linked to the northward shift of wintering areas induced by warming trends of the southern parts of Europe. However, the complexity of the likely processes leading to changing migratory strategies calls up for further research in both breeding and wintering areas.

Migration in butterflies: a global overview

Insect populations including butterflies are declining worldwide, and they are becoming an urgent conservation priority in many regions. Understanding which butterfly species migrate is critical to planning for their conservation, because management actions for migrants need to be coordinated across time and space. Yet, while migration appears to be widespread among butterflies, its prevalence, as well as its taxonomic and geographic distribution are poorly understood. The study of insect migration is hampered by their small size and the difficulty of tracking individuals over long distances. Here we review the literature on migration in butterflies, one of the best-known insect groups. We find that nearly 600 butterfly species show evidence of migratory movements. Indeed, the rate of 'discovery' of migratory movements in butterflies suggests that many more species might in fact be migratory. Butterfly migration occurs across all families, in tropical as well as temperate taxa; Nymphalidae has more migratory species than any other family (275 species), and Pieridae has the highest proportion of migrants (13%; 133 species). Some 13 lines of evidence have been used to ascribe migration status in the literature, but only a single line of evidence is available for 92% of the migratory species identified, with four or more lines of evidence available for only 10 species - all from the Pieridae and Nymphalidae. Migratory butterflies occur worldwide, although the geographic distribution of migration in butterflies is poorly resolved, with most data so far coming from Europe, USA, and Australia. Migration is much more widespread in butterflies than previously realised - extending far beyond the well-known examples of the monarch Danaus plexippus and the painted lady Vanessa cardui - and actions to conserve butterflies and insects in general must account for the spatial dependencies introduced by migratory movements.

Analysis of monitoring data where butterflies fly year-round

Butterfly Monitoring Schemes (BMSs) engage the public in conservation and provide data sets that cover broad geographical areas over long timescales. Most existing BMSs are in temperate climates; however, the Israeli Butterfly Monitoring Scheme (BMS-IL), established in 2009, is a notable exception as it encompasses a large climatic gradient from Euro-Siberian through Mediterranean to hyper-arid regions. Israel's climate poses challenges in analyzing data from year-round butterfly activity, as in other tropical or arid countries. The Regional Generalized Additive Model (Regional GAM) is a butterfly phenology and abundance model based on repeat visits throughout species' flight season. We tested the applicability of Regional GAM for species with complex flight seasonality (e.g., multivoltine) by comparing estimated abundance and seasonal indices for the full data set and rarefied subsets. We assessed the reliability of modeled flight seasons and compared abundance estimates per site resulting from biologically plausible and unreliable seasonal models. The reliability of Regional GAM rises with the number of observations, and the model tends to produce more biologically plausible models for species with simple phenologies (e.g., univoltine with a single peak in activity). Abundance estimates based on unreliable models produce values with inter-quartile ranges of 90%-153% compared with biologically plausible models, while peak time changes with an interquartile range of 0-22.5 d when comparing all rarefied models with the full data set. Regional GAM should be applied with great caution for rare species and those with a complex flight season, and the date of year start needs to be carefully chosen for species that are active year-round. We identified the key sources of error and propose an operational workflow to address them. With few adaptations, Regional GAM can support new BMSs in analyzing data where butterflies are active year-round, including tropical climates. We propose guidelines for analyzing BMS data for species or regions with long activity periods and complex phenologies.

Spatio-temporal ecological niche modelling of multigenerational insect migrations

Modelling ecological niches of migratory animals requires incorporating a temporal dimension, in addition to space. Here, we introduce an approach to model multigenerational migratory insects using time-partitioned environmental variables (by months and years) and time- and behaviour-partitioned records (breeding records to model reproductive habitat). We apply this methodology to modelling the Palearctic-African migratory cycle of the Painted Lady butterfly (Vanessa cardui), based on data encompassing 36 years (646 breeding sites from 30 countries). Each breeding record is linked to a particular time (month and year), and the associated values of the bioclimatic variables are used for an ensemble modelling strategy, to finally obtain monthly projections. The results show obligated movements, mostly latitudinal, for the species' successive generations across the overall range, and only scattered locations show high probabilities of reproduction year-round. The southernmost reproductive areas estimated for the Palearctic-African migratory pool reach equatorial latitudes from December to February. We thus propose a potential distribution for the winter 'missing generations' that would expand the V. cardui migration cycle to encompass about 15 000 km in latitude, from northernmost Europe to equatorial Africa. In summer, Europe represents the major temporary resource for V. cardui, while January and February show the lowest overall suitability values, and they are potentially the most vulnerable period for the species to suffer yearly bottlenecks. In summary, we demonstrate the potential of the proposed niche modelling strategy to investigate migratory movements of insects.

Tracking dragons: stable isotopes reveal the annual cycle of a long-distance migratory insect

Insect migration is globally ubiquitous and can involve continental-scale movements and complex life histories. Apart from select species of migratory moths and butterflies, little is known about the structure of the annual cycle for migratory insects. Using stable-hydrogen isotope analysis of 852 wing samples from eight countries spanning 140 years, combined with 21 years of citizen science data, we determined the full annual cycle of a large migratory dragonfly, the common green darner ( Anax junius). We demonstrate that darners undertake complex long-distance annual migrations governed largely by temperature that involve at least three generations. In spring, the first generation makes a long-distance northbound movement (further than 650 km) from southern to northern range limits, lays eggs and dies. A second generation emerges and returns south (further than 680 km), where they lay eggs and die. Finally, a third resident generation emerges, reproducing locally and giving rise to the cohort that migrates north the following spring. Since migration timing and nymph development are highly dependent on temperature, continued climate change could lead to fundamental changes in the biology for this and similar migratory insects.

Pollen metabarcoding as a tool for tracking long-distance insect migrations

Insects account for a large portion of Earth's biodiversity and are key players for ecosystems, notably as pollinators. While insect migration is suspected to represent a natural phenomenon of major importance, remarkably little is known about it, except for a few flagship species. The reason for this situation is mainly due to technical limitations in the study of insect movement. Here, we propose using metabarcoding of pollen carried by insects as a method for tracking their migrations. We developed a flexible and simple protocol allowing efficient multiplexing and not requiring DNA extraction, one of the most time-consuming part of metabarcoding protocols, and apply this method to the study of the long-distance migration of the butterfly Vanessa cardui, an emerging model for insect migration. We collected 47 butterfly samples along the Mediterranean coast of Spain in spring and performed metabarcoding of pollen collected from their bodies to test for potential arrivals from the African continent. In total, we detected 157 plant species from 23 orders, most of which (82.8%) were insect-pollinated. Taxa present in Africa-Arabia represented 73.2% of our data set, and 19.1% were endemic to this region, strongly supporting the hypothesis that migratory butterflies colonize southern Europe from Africa in spring. Moreover, our data suggest that a northwards trans-Saharan migration in spring is plausible for early arrivals (February) into Europe, as shown by the presence of Saharan floristic elements. Our results demonstrate the possibility of regular insect-mediated transcontinental pollination, with potential implications for ecosystem functioning, agriculture and plant phylogeography.