Decision Letter: Windborne Migration Amplifies Insect-Mediated Pollination Services

Youngsung Joo,Jason Chapman,Gao Hu

DOI: https://doi.org/10.7554/elife.76230.sa1

2022-01-01

Abstract:Article Figures and data Abstract Editor's evaluation Introduction Results Discussion Materials and methods Data availability References Decision letter Author response Article and author information Metrics Abstract Worldwide, hoverflies (Syrphidae: Diptera) provide crucial ecosystem services such as pollination and biological pest control. Although many hoverfly species exhibit migratory behavior, the spatiotemporal facets of these movement dynamics, and their ecosystem services implications are poorly understood. In this study, we use long-term (16-year) trapping records, trajectory analysis, and intrinsic (i.e., isotope, genetic, pollen) markers to describe migration patterns of the hoverfly Episyrphus balteatus in northern China. Our work reveals how E. balteatus migrate northward during spring–summer and exhibits return (long-range) migration during autumn. The extensive genetic mixing and high genetic diversity of E. balteatus populations underscore its adaptive capacity to environmental disturbances, for example, climate change. Pollen markers and molecular gut analysis further illuminate how E. balteatus visits min. 1012 flowering plant species (39 orders) over space and time. By thus delineating E. balteatus transregional movements and pollination networks, we advance our understanding of its migration ecology and facilitate the design of targeted strategies to conserve and enhance its ecosystem services. Editor's evaluation Hoverflies are a group of insects which provide crucial ecosystem services such as pollination and crop protection. Their migratory behavior in western countries is well characterized, but in eastern Asia, the annual summer monsoon provides a 'highway' of favorable winds for the airborne transport of migratory organisms, and the migration of hoverflies in this large region has not been well studied. This study addresses hoverfly migration in East Asia and its consequences using a variety of suitable methods and will be of great interest to insect migration biologists and pollination ecologists. https://doi.org/10.7554/eLife.76230.sa0 Decision letter Reviews on Sciety eLife's review process Introduction Migration plays a key role in the evolution and life history of many organisms, with insects being the most abundant, speciose, and economically important group of terrestrial migrants (Chapman et al., 2011). Across the globe, billions of insects annually undertake long-range movements. By transporting energy, nutrients, and other organisms between distant regions, insect migrants provide a multitude of ecosystem services and disservices (Chapman et al., 2012; Hu et al., 2016; Satterfield et al., 2020). Despite the important socioecological consequences of insect migration, research has primarily centered on a handful of large-bodied charismatic species or agricultural pests (e.g., monarch butterflies, locusts; Chapman et al., 2015). For most other taxa, there is a critical dearth of information. Hoverflies (Diptera: Syrphidae) are a speciose family of beneficial insects – deemed to be the second most important pollinators after bees (Branquart and Hemptinne, 2000; Rader et al., 2020). The larval stages of many hoverflies are effective predators of homopteran feeders, providing natural biological control across geographies and farming contexts (Tenhumberg, 1995; Tenhumberg and Poehling, 1995). Evidence to date suggests that hoverfly species are abundant diurnal migrants that deliver ecosystem services in both natural and man-made habitats (e.g., David, 1951; Aubert and Tiefenau, 1981; Wotton et al., 2019). Moreover, given that (migratory) hoverflies exhibit comparatively stable population numbers and transport pollen over long distances (Wotton et al., 2019), these species potentially can sustain pollination and pest control services in the face of a global insect decline (Sánchez‐Bayo and Wyckhuys, 2020; Powney et al., 2019). Yet, though (long-range) dispersal is a central determinant of their survival, hoverfly migration has only been intermittently studied since the 1950s. In order to effectively conserve these organisms and to raise their contribution to (agro-)ecosystem functioning, a more in-depth understanding needs to be gained of hoverfly migration. In recent years, several new technologies have helped to uncover insects’ seasonal migration patterns and population genetic structure. Stable isotope analysis, molecular genetics, tethered flight mill assays, insect radar, and aerial trapping have all yielded insights into the migration behavior of hoverflies (Raymond et al., 2014; Dällenbach et al., 2018; Wotton et al., 2019; Gao et al., 2020). Attempts have equally been made to capture the geographical extent and ecological impacts of hoverfly migration (Wotton et al., 2019). Most of these studies however originate from (a small area within) Europe, while virtually no information is available from other parts of the world. Also, as the nutritional ecology of most species waits to be deciphered, little is known about hoverfly–plant associations and how those are modulated by (long-range) migration dynamics. In China, approx. 580 hoverfly species have been described. These include the marmalade hoverfly Episyrphus balteatus (DeGeer) (Li et al., 2009), a common flower visitor in urban and agricultural settings across the Palearctic realm. Locally, (insect) migration primarily takes place within the East Asia monsoon climatic zone (Drake and Farrow, 1988). Owing to its geographical range, complex topography, and diverse agroecological conditions, this climatic zone constitutes an exceptional setting to study broad-scale migration dynamics of hoverfly species, for example, as compared to other parts of the globe (Wotton et al., 2019; Menz et al., 2019; Finch and Cook, 2020). In northern China, numerous insect species annually undertake (two-way) migration across the Bohai Strait, influenced by the East Asian monsoon cycle (Feng et al., 2003). In this study, we employed a suite of novel methodologies to characterize the migration dynamics of E. balteatus in China. More specifically, we conducted long-term (16-year) searchlight trapping on Beihuang (BH), a small isolated island with sparse vegetation in the center of the Bohai Strait, to clarify whether E. balteatus engages in long-range migration and to describe the ensuing migration patterns. Second, we deployed backward trajectory analysis and stable isotope analysis to infer the E. balteatus migration routes and source areas. Third, we employed a population genetics approach to compare the genetic makeup and demographic history of migrant and field-collected individuals throughout China. Fourth, we described E. balteatus host–plant associations by identifying the pollen grains attached to hoverfly bodies. Lastly, we paired molecular gut content analysis with high-throughput sequencing (HTS) to investigate the spatiotemporal distribution of its (flower) host plants across a broad geographic range. As such, our work characterized E. balteatus migration behavior and captured its broader ecological relevance, for example, in terms of flower visitation networks, pollination, or natural biological control services. Results Migration dynamics Despite the local availability of weedy host plants, intensive field surveys did not detect a presence of E. balteatus larvae on Beihuang (BH) island. Yet, night-time trapping consistently yielded E. balteatus adults from late April to October (Figure 1A) throughout the 16-year sampling period. Hoverfly abundance (or trap capture rate) showed important annual variation (Marginal R2 = 0.163, χ2 = 75.4, df = 15, p < 0.001) (Figure 1—figure supplement 1), with peak population sizes of 5068 and 2381 individuals in 2009 and 2010, respectively. Overall, E. balteatus trap capture rate exhibited a bimodal pattern with peak abundance from May to June and August to September, thus comprising two distinct migration stages. On an annual basis, E. balteatus migration covered a period of 151.2 ± 17.9 days (Supplementary file 1e). Figure 1 with 1 supplement see all Download asset Open asset Seasonal migration patterns of Episyrphus balteatus in Eastern China. (A) Annual migration dynamics, expressed as nightly searchlight-trap catches, of E. balteatus on Beihuang Island (BH; Bohai Gulf, China) from April to October 2003–2018. (B) Endpoints of backward trajectories of BH-caught E. balteatus individuals during mass migration events over 2003–2018 for a 12-hr flight duration. Darker colors indicate a higher density of endpoints at a particular location. The left panel (i.e., orange colors) reveals the possible source areas of late spring and summer immigrants, while the right panel (i.e., blue colors) indicates those of autumn immigrants. Figure 1—source data 1 Nightly catches of E. balteatus in the searchlight trap on BH from April to October 2003–2018. https://cdn.elifesciences.org/articles/76230/elife-76230-fig1-data1-v2.xlsx Download elife-76230-fig1-data1-v2.xlsx Next, the possible origin of E. balteatus migrants on BH was identified through backward trajectory analysis using the HYSPLIT model. For ‘mass migration’ events over 2003–2018 (Supplementary file 1a), spring–summer migrants primarily originated in southern areas, while autumn migrants arrived from northern areas. During spring–summer, a total of 9069 valid endpoints endpoints were identified and more than 90% endpoints where located south of BH, that is, in Shandong, Jiangsu, Henan, and Anhui provinces; during autumn, of 4945 valid endpoints, 85% endpoints were distributed in areas north of BH such as Liaoning (79.78%) likely acted as key source areas (Figure 1B; Supplementary file 1f). To corroborate the above patterns, wings of BH migrants and wild-caught individuals from sites across China were subject to hydrogen isotope analysis. On BH, greater variability in δD values was recorded for autumn versus spring–summer migrants over the entire study period (Wilcoxon rank-sum test; W = 11.93, df = 254.13, p < 0.0001) (Figure 2B). By comparing the above δD values with its established precipitation gradient, source areas of E. balteatus migrants were identified. Considering how δD values in field-collected populations divert from latitudinal gradients (Figure 2C), seasonal differences became apparent upon a geographical grouping of E. balteatus populations. Overall, spring-caught adults showed higher δD values than those captured during autumn (Wilcoxon rank-sum test; NE subgroup: W = 3232.5, p < 0.001; NP: W = 13804, p < 0.001; YzP: W = 3440, p < 0.001) (Figure 2D). Isotope patterns reflected how hoverfly adults engage in bidirectional migration over hundreds of kilometers. Figure 2 Download asset Open asset Stable hydrogen isotope analysis of E. balteatus. (A) Sampling locations. (B) Seasonal incidence of mean δD values in the wings of migratory E. balteatus adults, as recorded for different migration stages over 2014–2018 on Beihuang Island (BH). Double asterisks (***) indicate a statistically significant difference (p < 0.01, Wilcoxon test). (C, D) Seasonal incidence of mean δD values in the wings of wild-caught individuals from each sampling location or geographic region during April to October 2017–2018, respectively. Locations are ranked as per their longitudinal position. PE, Puer; NC, Nanchang; WH, Wuhan; TZ, Taizhou; YZ, Yangzhou; YL,Yangling; XX, Xinxiang; XJ, Xiajin; CZ,Cangzhou; YT, Yantai; LF, Langfang; CD, Changdao; SY, Shenyang; CC, Changchun; HB, Harbin. Population genetics To gain genetic evidence of its regional migration, we described E. balteatus genetic diversity and population structure using one mitochondrial DNA gene (Cytb) and two nuclear DNA genes (i.e., 18s rRNA, 28s rRNA). Upon analysis of 670 field-collected specimens and 133 light-trapped individuals (representing a respective 16 and 2 populations), high haplotype diversity and low nucleotide diversity was recorded (Table 1). Based upon Cytb sequences, 83 haplotypes were identified among 530 individuals, with haplotype diversity (Hd) ranging from 0.0800 (XJ) to 0.857 (SH) (total = 0.351) and nucleotide diversity (π) from 0.000110 to 0.00536 (total = 0.00178), respectively. Conversely, the concatenated nuclear gene possesses improbably high haplotype diversity and low nucleotide diversity. Up to 145 haplotypes were detected among 260 individuals with Hd ranging from 0.810 to 1 (total = 0.961) and π ranging from 0.00392 to 0.01384 (total = 0.000008). Table 1 Genetic diversity indices of 18 E. balteatus populations based on Cytb and 18S-28S rRNA gene. For each population (and sampling location), the following metrics are reported: N, sample size; S, number of segregating sites; H, number of haplotypes; Hd, haplotype diversity; K, average number of differences; Pi, nucleotide diversity. For the site names see Figure 2. Site codeCytB18S-28S rRNANSHHdKPiPiJCNSHHdKPiPiJCLF233770.462453.920950.005360.00543122580.848485.515150.004520.00454XX42650.336820.452960.000620.000621631100.8254.791670.003920.00394HB23120.086960.086960.000120.000121828100.810465.607840.004590.00461WH341050.22460.588240.00080.00081123290.909096.909090.005660.00569CC162870.6253.50.004780.004821122100.981825.636360.004620.00464NC211250.352381.657140.002260.002281545120.942869.619050.007880.00793SY451660.250510.753540.001030.001042149180.980959.26190.007590.00764CF272180.458692.028490.002770.0027916391619.008330.007380.00742XJ25120.080.080.000110.000111831130.954256.718950.00550.00553YT292250.317731.950740.002660.002681634120.958337.80.006390.00643YP18830.215690.888890.001210.001228388112.892860.010560.01064YL30460.310340.331030.000450.0004514351418.560440.007010.00705TZ28230.140210.142860.00020.0002N/AN/AN/AN/AN/AN/AN/APE18740.470591.797390.002460.00247119811122.963640.018810.01922YZ31450.245160.258060.000350.00035107290.9777815.622220.012790.01309SH7850.857142.285710.003120.00313N/AN/AN/AN/AN/AN/AN/ACDI5243160.524132.301660.003140.0031625103210.9766716.90.013840.01413CDII6148170.482511.700550.002320.002343765280.981989.151650.00750.00754 Though phylogenetic analyses showed four and five distinct clades among the Cytb and nuclear haplotypes, haplotype and geographical origins were not linked (Figure 3—figure supplement 1). Median-joining also did not reveal geographical clustering. Instead, a star-like pattern was displayed with the most common, ancient haplotypes in the center (Figure 3—figure supplement 2). Most haplotypes were unique to individuals and populations, while only 9 (out of 83 Cytb haplotypes) and 30 (out of 145 nuclear haplotypes) were shared. In each population, shared haplotypes occurred at 63–100% frequencies for Cytb and 33–88% for nuclear genes (Figure 3A, B). Only 1.8% and 11.1% of genetic variation for the respective Cytb and nuclear gene could be attributed to variation among populations (Supplementary file 1g). Low pairwise FST values between different localities equally reflected low levels of genetic differentiation (Figure 3—figure supplement 3). When describing effective population sizes and migration rates between (geographically grouped) populations, high levels of interpopulation gene flow were recorded (Figure 3C, D). Migration between E. balteatus populations was asymmetrical and exhibited a general migration trend toward the Yangtze basin, with north > center, south > center, and west > east as dominant directions. Migration thus enables genetic mixing among E. balteatus populations from remote origins. Given the negative Tajima’s D and Fu’s FS (Cytb: Tajima’s D = −2.72034, Fu’s FS = −27.4859; concatenated 18S-28S rRNA: Tajima’s D = −2.33710, Fu’s FS = −34.456) and unimodal distribution for both markers (Figure 3—figure supplement 4), E. balteatus populations possibly experienced recent population expansion. Figure 3 with 4 supplements see all Download asset Open asset Population genetic structure of E. balteatus populations sampled in 18 locations across China. Analyses are either based on one mitochondrial Cytb gene (upper panels, i.e., A, C) or on combined 18s and 28s rRNA nuclear genes (lower panels, i.e., B, D). (A, B), Spatial distribution of E. balteatus haplotypes. At each given location, a pie chart shows the proportional abundance of haplotypes. (C, D) Migration dynamics of different populations as estimated through MIGRATE-N. Individual arrows represent the prevailing migration direction and arrow thickness is proportional to the number of migrants. SW, Southwestern Region; YzP, Yangtze Plain; NP, Northern Region; and NE, Northeastern Region, as defined in Figure 2. Palynological analysis Morphological characteristics and barcode markers were used to identify pollen species attached to migrating E. balteatus adults, and to infer the associated movement and flower visitation patterns. Among 1014 BH adult migrants collected during 2015–2018, 32% had pollen grains adhering to the body surface. Using a combination of DNA sequences, pollen morphology and plant distribution records, a total of 46 pollen species representing at least 42 plant genera and 26 families were identified (Figure 4; Supplementary file 1h). Out of these, 10 were identified to genus level and the remained to species level. Few adults carried pollen from multiple plant species. Pollen-bearing plants mainly pertained to Asteraceae (12), Moraceae (4), and Celastraceae (3), and were primarily herbaceous as compared to woody plants (χ2 = 112.26, df = 1, p < 0.001). It should be noted that most plant hosts were identified from one single hoverfly, with only few adults carrying pollen from different plant species. Figure 4 Download asset Open asset Scanning electron microscopy (SEM) microphotographs of pollen grains attached to E. balteatus migrants on BH during 2014–2018. 1. Ailanthus altissima; 2. Cotinus coggygria; 3. Forsythia suspensa; 4. Prunus avium; 5. Brassica L.; 6. Morus alba; 7. Citrus sinensis; 8. Descurainia sophia; 9. Euonymus L.; 10. Taraxacum L.; 11. Sedum japonicum; 12. Populus cathayana; 13. Celastrus orbiculatus; 14. Daucus carota; 15. Chenopodium L.; 16. Castanea mollissima; 17. Amorpha fruticosa; 18. Diospyros lotus; 19. Ziziphus jujuba; 20. Cirsium setosum; 21. Neoshirakia japonica; 22. Flueggea L.; 23. Maclura pomifera; 24. Rumex L.; 25. Euonymus L.; 26. Schisandra chinensis; 27. Eleusine indica; 28. Actinidia kolomikta; 29. Cannabis sativa; 30. Humulus scandens; 31. Helianthus annuus; 32. Persicaria orientalis; 33. Adenophora trachelioides; 34. Gypsophila paniculata; 35. Artemisia L.; 36. Rubia cordifolia; 37. Rubia cordifolia; 38. Artemisia L.; 39. Artemisia L.; 40. Artemisia L.; 41. Allium tuberosum; 42. Tripolium vulgare; 43. Ambrosia trifida; 44. Sorghum bicolor; 45. Aster tataricus; 46. Chrysanthemum zawadskii. The scale bar is shown on the bottom of each photograph. Furthermore, temporal patterns were recorded in the pollen adherence ratio, with the respective highest and lowest levels recorded in October (50%) and April (17%). Also, the identity of pollen grains that adhered to E. balteatus bodies equally exhibited seasonal variation. Specifically, 29 pollen species were identified from spring–summer migrants, with Taraxacum mongolicum, Ailanthus altissima, Amorpha fruticose, and Chenopodium giganteum the most common plant hosts (i.e., accounting for a 72% total carrier rate). Conversely, 17 pollen species were identified from autumn migrants, with Artemisia L., Chrysanthemum zawadskii and Ambrosia trifida accounting for a 82.2% total carrier rate (Supplementary file 1c). Lastly, specific plant taxa were associated with certain migration stages. Plants endemic to certain ecological zones helped to pinpoint migration origins (Jones and Jones, 2001). During spring–summer, the presence of Citrus L., Sedum japonicum, or Euonymus myrianthus endemic to central and southern China hinted at migration origins in south-central China. Conversely, during autumn, the presence of C. zawadskii mirrored potential migration origins in northeastern China (Figure 5). These results further reflected northward E. balteatus migration flows during spring with return movements in autumn. Figure 5 Download asset Open asset Geographical distribution of pollen-bearing plant species that help delineate E. balteatus migration patterns. Maps show the (district-level) geographical distribution of individual plant species. Molecular gut content analysis Under laboratory conditions, a polymerase chain reaction PCR-based assay was developed to identify the plant species that are consumed by E. balteatus. Using DNA extracts from excised syrphid abdomens and ITS2-targeted plant primers, we successfully amplified the expected band with specific size of all extracts. BLAST analysis confirmed that sequencing products were indeed the three target plant species (M. sativa, H. scandens, and H. annuus) that were consumed by E. balteatus adults. Additionally, DNA of all three ingested plants could be detected (at high frequency) in E. balteatus guts for up to 9-day postfeeding, the maximum digestion time in our experiments. Moreover, qPCR methods further yielded standard curves to quantify the degradation of the plant DNA through digestion, and a negative exponential equation to compute the time to reach a 50% detection probability (Figure 6). Detection times ranged between 5.14 and 5.79 hr for the three plant species. Once more than 8.9 copy numbers/μl of ingested plant DNA were present in E. balteatus guts, amplification product could also be visualized by ethidium bromide staining. Figure 6 with 1 supplement see all Download asset Open asset Detection of three different plant DNA in the guts of E. balteatus adults at different times after ingestion by qPCR analysis. Error bars at each point on the curves represent the standard error of replicates. Spatiotemporal dietary shifts By pairing gut content analysis with HTS (Illumina MiSeq), (floral) diet profiles of 180 light-trapped E. balteatus migrants on BH and 436 field-collected individuals at 19 different locations in China were elucidated. MiSeq paired-end sequencing of the ITS2 gene yielded 3,348,773 raw reads from 616 samples. After quality filtering and reference-based chimera removal, a total of 2,952,448 sequences remained, ranging from 100 to 28,582 sequences per sample. These valid reads were clustered into Operational Taxonomic Units (OTUs) for further analysis (Supplementary file 2). According to the taxonomic assignment, almost all sequences were identified at least to the genus level. Overall, 1012 plant species belonging to 39 orders, 91 families, and 429 genera were identified, with Asteraceae (84), Poaceae (37), Apiaceae (21), and Fabaceae (21) (count at genus level) the dominant families accounting for 38.0% of all plant taxa (Supplementary file 3). Further analysis of plant characteristics of identified plant taxa revealed important similarities with the above palynological analysis, in which more Angiosperm and Dicotyledon plants were recorded than Gymnosperms (χ2 = 393.76, df = 1, p < 0.0001) or Monocotyledons (χ2 = 253.04, df = 1, p < 0.0001). Woody plant hosts were more common than herbaceous ones (χ2 = 61.94, df = 1, p < 0.0001) (Figure 7). These results indicated the hoverflies might have visited herbaceous plants more often than woody plants. Figure 7 with 8 supplements see all Download asset Open asset Sankey diagram summarizing the host plants of E. balteatus, as identified via DNA metabarcoding (i.e., for BH migrant individuals) or DNA-based gut content analysis (i.e., for field-collected individuals from 19 sites). Taxonomic structure is shown for 429 plant genera, with individual numbers referring to the number of counts of particular taxa. Figure 7—source data 1 List of the host plants of E. balteatus, as identified via DNA metabarcoding (i.e., for BH migrant individuals) or DNA-based gut content analysis (i.e., for field-collected individuals from 19 sites). https://cdn.elifesciences.org/articles/76230/elife-76230-fig7-data1-v2.xlsx Download elife-76230-fig7-data1-v2.xlsx The diversity and abundance of host plant communities varied across geographical sites and sampling times (Figure 7—figure supplements 1 and 2). In particular, significantly more plant taxa were recorded in the BH migratory population (36 orders, 76 families, and 320 genera) as compared to field-collected populations. Moreover, Venn diagram displayed that only a small number of plant taxa were unique to particular (geographical) populations at both the phylum and genus levels, most of the these were widespread among two or more regions, and there are several taxa that included six families Apiaceae, Poaceae, Asteraceae, Brassicaceae, Fabaceae, Ranunculaceae, and one genera Brassica observed in all of the populations across different geographical sites and sampling times (Figure 7—figure supplement 3a, b). The results of hierarchical cluster analysis revealed that the overall diet composition of this species exhibited weak temporal differences upon a spatiotemporal grouping of E. balteatus populations at the genus level (Figure 7—figure supplement 4). These results strongly indicated that geographical location thus only affected pollen transport networks to small extent. Importantly, the host plant composition of the BH migratory population almost equals the sum of 19 field-collected populations (Figure 7—figure supplement 3c, d), which suggested the migratory behavior of some insects has important role in shaping flower visitation network structures. Next, upon analyzing host plants derived from light-trapped E. balteatus migrants on BH (36 orders, 76 families, and 320 genera), we further assessed how sampling time did affect the primary (plant) foraging resource. The results showed that the overall plant community composition exhibited important within-year variability, as reflected by the host-associated community separation visualized in the principal coordinates analysis (PCA) profile and the significance test based on analysis of similarities (ANOSIM; p < 0.001; Figure 7—figure supplement 5a). Moreover, Venn diagrams indeed displayed a certain number of unique plant species in all E. balteatus groups, with a clear shift in plant taxa over time (Figure 7—figure supplement 5b). Based on the average relative abundance, Asteraceae was the most important plant families during the whole study period (Figure 7—figure supplement 6a). At the genus level, during spring (April to May), the genera Pinus and Capsella are the important food resources; during summer (June to August), Flueggea and Maclura became more common pollen hosts; during autumn (September–October), Ajania and Artemisia was the most represented species (Figure 7—figure supplement 6b), the importance of these main plant taxa during specific seasons was also confirmed by random forest machine learning (Figure 7—figure supplement 7). Also, in addition to the above five plants that are informative of migration identified by pollen analysis, more geographically confined species such as Acer pictum, Pittosporum truncatum, Pseudotsuga sinensis, and Xylosma congesta, have been identified by applying this method (Figure 7—figure supplement 8). Discussion Several European studies have unveiled how hoverflies engage in long-distance migration, yet little is known about such phenomenon in other parts of the world. In this study, we relied upon different analytical approaches to delineate (long-range) dispersal patterns of the widely distributed E. balteatus in Asia. Specifically, combining spatially explicit data from searchlight trapping, trajectory analysis, molecular palynology, and metabarcoding, we showed how E. balteatus annually engages in bidirectional migration over hundreds of kilometers. Our study further revealed high genetic diversity and extensive genetic mixing at a continental scale, hinting at the superior adaptability and plasticity of this species. Pollen marker and next-generation sequencing NGS-based gut content approaches further illuminated how E. balteatus exploits a spatiotemporally diverse community of pollen-bearing plants. In view of the unrelenting global environmental degradation and the precipitous decline of insect populations, our work helps to guide interventions to conserve this beneficial insect species and to safeguard its vital ecosystem services. BH island has proven to be a highly suitable site to monitor insect migration in eastern Asia, because of its unique geographical location and local habitat condition. BH island is positioned within a major migration corridor at considerable distance from mainland China, and contains no arable land or extensive vegetation. From 2003 onwards, BH has lent itself to decipher the migration behavior of more than 100 insect species (e.g., Guo et al., 2020), including agricultural pests such as Helicoverpa armigera and natural enemies such as the dragonfly Pantala flavescens (Feng et al., 2009; Cao et al., 2018). In the present study, drawing upon (BH-specific) light-trapping records, we demonstrate how E. balteatus regularly migrates across the Bohai Gulf from late April to October each year. To our best knowledge, this constitutes the first evidence for transregional migration of hoverflies in eastern Asia. Meanwhile, its population dynamics are in accordance to those of other hoverfly species in Europe and elsewhere (e.g., Gao et al., 2020). As equally noted for other migratory insects (e.g., Guo et al., 2020; Stefanescu et al., 2013), important inter- and intraannual variation occurred in E. balteatus trap catches. Such variation can possibly be ascribed t