Next Article in Journal
Molecular Characterization of Heat-Induced HSP11.0 and Master-Regulator HSF from Cotesia chilonis and Their Consistent Response to Heat Stress
Previous Article in Journal
Molecular Characteristics of Fat Body Protein 1 in the Oriental Fruit Fly, Bactrocera dorsalis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evidence for Range Expansion and Origins of an Invasive Hornet Vespa bicolor (Hymenoptera, Vespidae) in Taiwan, with Notes on Its Natural Status

1
Taiwan Forestry Research Institute, Council of Agriculture, Executive Yuan, Taipei City 100051, Taiwan
2
Faculty of Life Sciences, Kyoto Sangyo University, Kyoto City 603-8555, Japan
3
Endemic Species Research Institute, Council of Agriculture, Executive Yuan, Nantou County 552203, Taiwan
4
Department of Plant Medicine, National Chiayi University, Chiayi City 600355, Taiwan
*
Authors to whom correspondence should be addressed.
Insects 2021, 12(4), 320; https://doi.org/10.3390/insects12040320
Submission received: 1 March 2021 / Revised: 23 March 2021 / Accepted: 30 March 2021 / Published: 2 April 2021

Abstract

:

Simple Summary

The invasive hornet Vespa bicolor Fabricius was first discovered in Taiwan in 2003 and was not confirmed to have been established until 2014. This study was conducted in order to (1) assess the current status of V. bicolor abundance, dispersal, seasonality, and possible impact on honeybee (Apis mellifera Linnaeus) in Taiwan; (2) and to trace the origins of Taiwan’s V. bicolor population. To assess V. bicolor abundance, we used visual surveys, sweep netting, and hornet traps in four known ranges in northern and central Taiwan from 2016 to 2020. Additionally, to understand V. bicolor dispersion, we mapped environmental data using ArcGIS, and to predict future V. bicolor range, we used ecological niche modeling. The results show that V. bicolor has stable populations in three areas in northern and central Taiwan, and mainly preys on Apis mellifera. Our analyses suggest samples from Southeastern China as having the closest relation in DNA sequences with Taiwan’s V. bicolor population. Due to the negative economic and ecological impacts of V. bicolor in Taiwan, our findings shed light on the value of monitoring and controlling its populations, rather than working exclusively towards elimination.

Abstract

The invasive alien species (IAS) Vespa bicolor is the first reported hornet that has established in Taiwan and is concerning as they prey on honeybee Apis mellifera, which leads to colony losses and public concerns. Thus, the aim of this study was to assess the current status of V. bicolor abundance, dispersal, and impact and to trace the origins of Taiwan’s V. bicolor population. Our studies took place in five areas in northern to central Taiwan. We used mtDNA in the phylogenetic analyses. Field survey and ecological niche modeling (ENM) were used to understand the origins and current range of the invasive species. Two main subgroups of V. bicolor in the phylogenetic tree were found, and a clade with short branch lengths in Southeastern China and Taiwan formed a subgroup, which shows that the Taiwan population may have invaded from a single event. Evidence shows that V. bicolor is not a severe pest to honeybees in the study area; however, using ENM, we predict the rapid dispersion of this species to the cooler and hilly mountain areas of Taiwan. The management of V. bicolor should also involve considering it a local pest to reduce loss by beekeepers and public fear in Taiwan. Our findings highlight how the government, beekeepers, and researchers alike should be aware of the implications of V. bicolor’s rapid range expansion in Taiwan, or in other countries.

1. Introduction

Invasive alien species (IAS) are non-native species that often cause adverse consequences in their new environment [1]. IAS should be monitored and controlled to prevent negative impact on agriculture and ecosystem function, as well as social and economic impacts on human beings [2]. The sheer diversity and niche breadth of insects allow them to play a significant role as IAS globally. Despite this, they are far understudied compared with plants and vertebrates. In addition, a recent review found that among the limited studies on IAS insects, nearly 18% were on the fire ant, Solenopsis invicta Buren, and up to 70% of the studies were conducted in North America [3]. These findings beg the need for studying IAS in Asia, especially in the North Asian region. Among insects, hornets are dominant and increasingly concerning players, as they can (1) be aggressively harmful to humans and other animals [4] and (2) spread easily due to nesting substrate availability during international trade [5,6].
Hornet has often become an example of an invasive species abroad. The yellow-legged hornet Vespa velutina Lepeletier is a widely distributed species indigenous to Southeast Asia, and it invaded Korea, Europe, and Japan [7,8,9,10,11]. Recently, there has been a lot of attention given to the giant Asian hornet Vespa mandarinia Smith, as it has invaded parts of North America [12]. This hornet species is the largest in Asia and has sparked public concern due to its painful sting and detriment to honeybee (Apis mellifera Linnaeus) populations. Hornet belongs to the family of social wasps in the genus Vespa (family: Vespidae), which comprises 22 species, all which were restricted to Asia until introduction and invasions occurred [13]. These species are important for the control of other insects in their native range [14], yet can be destructive to those outside it; thus it is critical to understand the biology and status of local species to avoid such perils.
In 2003, Taiwan first reported the invasive hornet species Vespa bicolor Fabricius [15], known as two-colored hornet or black shield hornet in Chinese [16,17]. V. bicolor is widely distributed in India, Bhutan, Cambodia, China, Laos, Myanmar, Nepal, Thailand, and Vietnam [18]. The existence of V. bicolor in Taiwan is concerning as it is known to prey on A. mellifera, which leads to colony loss and economic unfavorability [15]. Numerous hornets of the Vespa spp. have been reported to prey on honeybees in Asia [19,20], such as the hornets Vespa affinis (Linnaeus), Vespa basalis Smith, V. bicolor, Vespa ducalis Smith, V. mandarinia, and Vespa velutina in Taiwan [15]. These hornets are very adaptable and exhibit a wide range of prey capture strategies, from group hunting to sit and wait [4,15]. In addition, hornets, honeybees, and honey buzzards evolve complex interactions between prey and predator. For example, the migrant Oriental honey-buzzard Pernis ptilorhynchus orientalis Taczanowski preys on hornets (V. velutina), wasps (Polistes spp., Vespidae), and honeybees (A. mellifera) [21]; and a few populations are adapting to become summer residents in Taiwan [21,22]. Besides, the European honey buzzard Pernis apivorus (Linnaeus) has adapted to prey on hornet V. velutina [23].
Knowledge of IAS origin through the analysis of DNA and ecological landscape structures can help decision makers to draft risk assessment and management strategies [24,25]. The studies used mitochondrial and microsatellite DNA for V. velutina samples, and the origins of Japan, South Korea, and Europe populations were traced [10,11,26]. Furthermore, an intraspecies genetic study on V. velutina nigrithorax du Buysson that tested microsatellite variation in colonizing populations that were introduced via a very small number of propagules provided valuable information [26]. Moreover, ecological niche modeling (ENM) has been widely used for predicting the spread of IAS and for simulating spatial distribution. Using ENM and data from native and invaded ranges to assess the potential invasion risk of the invasive hornet V. velutina nigrithorax has predicted a noticeable remark on beekeeping activities that are threatened factors [27].
This study followed our 2014 report of V. bicolor establishment in Taiwan [15]. We used phylogenetic analyses, field survey, and ENM to understand the origins and current range of the invasive species. For the phylogenetic analyses, we sequenced mitochondrial DNA from V. bicolor samples to compare the genetic variations across their range. For the field survey, we conducted proactive monitoring of the V. bicolor activities to know the current temporal and spatial occurrences. For the ENM, the geographical locations were included in the model. These records contained long-term field surveys: (1) collection records from 2011 to 2013 of our study [15], (2) intermittent notifications of V. bicolor activity ranges from beekeepers or civilians from 2014 to 2020, (3) and survey data in this study from 2016 to 2020. We intend for our study to inform future control and management strategies.

2. Materials and Methods

2.1. Genetic Variations across V. bicolor Ranges

Since previous reports of V. bicolor hornets were predominantly from apiaries or from their own nests in northern and central Taiwan, the samples were collected for phylogenetic analyses there. The collections were carried out in Taoyuan City and the counties of Hsinchu and Miaoli from 2016 to 2020 discontinuously (Table 1). We collected a total 23 female individuals from different locations, 7 from the apiaries, 2 from flowers, and 14 from their own nests (Table 1). The samples were preserved in 95% EtOH and stored at 4 °C for subsequent DNA analysis.
The genomic DNA was extracted from a single midleg of the V. bicolor samples using an DNA Extraction Kit (Fisher Biotec, Wembley, Western Australia). The target DNA was amplified using the universal LCO1490/HCO2198 primers to get the COI region [28]. However, not all samples were successful in using the recommended protocol with LCO/HCO primers, as mentioned in the studies [29,30]. We also used the primers suggested by van Houdt, L1440d (5- TYT CAA CWA ATC ATA ARG ATA TTG G -3) and H2123d (5- TAW ACT TCW GGR TGW CCA AAR AAT CA -3), which can improve the amplification efficiency [30]. The primers used in the amplification were effective for our difficult samples, and we obtained the same fragment of the COI region [31].
Polymerase chain reaction (PCR) was performed with a thermocycler (Applied Biosystems 2720 Thermal Cycler, Foster, CA, USA). Each PCR reaction mixture contained 0.5 µL of Vas Taq (5 units/µL) (Bionovas Biotechnology Co., Ltd., Toronto, Canada), 5 µL of 10× Taq buffer (Bionovas Biotechnology Co., Ltd., Toronto, Canada), 4 µL of 2.5 mM dNTPs (Won-Won Biotechnology Co., Ltd., New Taipei City, Taiwan), 1 µL of 10 µM forward and reverse primers, 1 µL of template DNA, and deionized water (up to 50 µL). The PCR amplification consisted of an initial denaturation at 94 °C for 5 min, followed by 40 cycles of denaturation at 94 °C for 30 s, annealing at 50 °C for 30 s, extension at 72 °C for 30 s, and final extension at 72 °C for 10 min. PCR products were separated by electrophoresis in 1% agarose gel, stained with ethidium bromide, and visualized with a UV transilluminator (Major Science, CA, USA) The genes were sequenced with Sanger dideoxy sequencing by a commercial genome sequencing company (Genomics, New Taipei City, Taiwan).
The software MEGA X was used to compare [32], align, and trim the double-stranded base pair sequences of the same arrangement of COI sequences. The supplemental sequences from Huangshan in China, Hong Kong, Nepal, and Vietnam were provided by one of our authors (J.T.). We also included sequence data on V. bicolor from GenBank (accession nos. KT257112, KF933079, Vietnam; KJ735511, Beijing) [33] and used as outgroups [34] Vespa simillima (accession no. NC046020), V. velutina (accession no. NC035146), and V. mandarinia (accession no. KR059904). The model of maximum likelihood method and Tamura 3-parameter, 1000 bootstrap numbers were used to construct phylogenetic trees. DnaSP v6 software was used to detect V. bicolor COI haplotypes [35]. For this analysis, we also included data on COI sequences from GenBank (V. bicolor accession nos. KT257112, KF933079, Vietnam; KJ735511, Beijing). The haplotype network was calculated using Network v10 and plotted by hand [36].

2.2. Occurrences, Daily and Seasonal Activities of V. bicolor

The studies took place in five areas from northern to central Taiwan: the counties Hsinchu, Miaoli, and Nantou and the cities Taoyuan and Taichung (Figure 1). Investigation by visual surveys was carried out at 13 apiaries in the cities of Taoyuan and Taichung and the counties of Hsinchu, Miaoli, and Nantou from April to November discontinuously in 2019 and 2020 (Table 2). In addition, 7 apiaries were set up with hornet traps, 1 to 2600 cc PET bottles with two entrances and filled 50% with alcohol-based bait attractant (Table 2) [37]. Moreover, beekeepers and affiliated persons were trained to distinguish V. bicolor hunting behaviors and activities, and we recorded their findings and notifications. The apiary (ca. 30 hives) of V. bicolor occurred in Tongluo, Miaoli, where we used a visual survey monthly from May to November in 2019 and 2020. Each V. bicolor count lasted for 1 h by two researchers. In addition, the apiaries (each apiary with 30 hives) of V. bicolor were observed to occur in Tongluo, Miaoli, and Beipu, Hsinchu, using mark and recapture for 5 days (twice a day) during the peak occurrence of the year (i.e., from September to October 2020). Each survey was done by two researchers for 1 h, and the interval was 2 h between the first and second survey. The V. bicolor captured during the first survey was marked with a color acrylic pen marker, and then the V. bicolor captured during the second survey, and the number of captured was counted and recorded. The daily and seasonal activities of V. bicolor at the apiary were estimated.

2.3. Ecological Niche Modeling

2.3.1. Species Occurrence and Environmental Data

Twenty-seven geographical locations in the counties and cities of Taoyuan, Hsinchu, Miaoli, and Taichung, where V. bicolor was most likely to occur, were included in the model (Figure 1). Environmental predictors refer to the critical factors that affect species distribution. To determine which environmental variables influence the distribution of V. bicolor, we included in our model 19 bioclimatic variables [38] and 1 biophysical variable (elevation) with a 30 s (ca. 1 km) spatial resolution, downloaded from a WorldClim dataset (www.worldclim.org) (accessed on 3 November 2020). The mean slope and aspect were calculated from the elevation layer using a surface analysis tool from the Spatial Analyst toolbox in ArcGIS 10.6 (Esri) (Environmental Systems Research Institute, Redlands, CA, USA). Then, we used the Geomorphometry and Gradient Metrics toolbox (version 2.0) [39] to generate compound topographic indices (representing humidity) and heat load indices (representing solar radiation intensity). To reduce the multicollinearity among the 24 variables (Table S1), the variance inflation factor (VIF) implemented in the “usdm” package in the R platform version 3.6.3 was used to exclude predictors with VIF values >10 [40].

2.3.2. Ensemble Modelling

We employed an ensemble modelling approach to predict V. bicolor distribution using R with the “sdm” package [40]. We selected five algorithms: generalized linear models (GLMs), multiple adaptive regression splines (MARS), random forest (RF), maximum entropy (MaxEnt), and support vector machine (SVM). As all these models require background data (pseudo-absence points), we generated a random sampling of 5000 locations across the study area. Models were evaluated using k-fold cross-validation with 10 folds and 10 replications for each algorithm. We projected each of the models using 70% of the training data and 30% for evaluation. The performance of the models was evaluated using threshold-independent indices of the area under the curve (AUC) of a receiver operating characteristic and the true skill statistic (TSS) [41,42]. The AUC ranged from 0.5 to 1, with a higher value denoting greater prediction accuracy. An AUC value ≥0.7 represents reliable prediction accuracy [43]. The relative importance of predictor variables was estimated using the function getVarImp of the “sdm” package. Moreover, to obtain the consensus predictions, we used the function “ensemble” based on weighted TSS values. Finally, the output map represents the occurrence probability for each pixel. We extracted the pixels with maximum occurrence probability, representing probabilities of more than 75% as suitable habitats [44].

3. Results

3.1. Genetic Variations across V. bicolor Ranges

The phylogenetic tree of the COI region near the 5′ end showed that all V. bicolor hornets from Beijing and Huangshan, Hong Kong, Nepal, Taiwan, and Vietnam were closely related to form a monophyletic clade (Figure 2A). Two main subgroups diverged: a clade with short branch lengths in Huangshan, Hong Kong, and Taiwan formed the first subgroup, while the Nepal samples and a clade in Beijing and Vietnam formed the second. All individuals collected in Miaoli, Hsinchu, and Taoyuan in Taiwan have identical sequences. The Beijing and Vietnam group had a relatively distant relationship with the Nepal group. The phylogenetic results show a V. bicolor population in Taiwan, suggesting it may come from the same source near Huangshan and Hong Kong in the Southeastern China area. The Taiwan sequences were submitted to GenBank (accession nos. MW455061–MW455084) (Table 1). Moreover, we found seven haplotypes in the COI sampling of V. bicolor (from Hap_1 to Hap_7) (Figure 2B and Figure S1). All Taiwan samples were detected in haplotype Hap_3, with three mutations linked to Hong Kong samples (Hap_5) and Huangshan samples (Hap_6) to form a branch. The second and third branches were Beijing samples (Hap_2) and Nepal samples (Hap_4), respectively. The Vietnam haplotype includes Hap_1 and Hap_7, which was linked to two median vectors. The haplogroups radiate from Hap_1 and Hap_7 to form a tripodlike pattern, indicating that the Vietnam samples represent ancestral haplotypes.

3.2. Current Status, Daily and Seasonal Activities of V. bicolor

The field survey recorded six hornet species, V. affinis, V. basalis, V. bicolor, V. ducalis, V. mandarinia, and V. velutina, in 13 apiaries in the city of Taoyuan and the counties of Hsinchu, Miaoli, and Nantou. The body size of V. bicolor (15–19 mm) is smallest, but it is a fast-flying hornet among these Vespa species, so it is not easy to spot and catch around the hives. According the beekeepers and affiliated persons’ notifications, a native Taiwanese wasp, Polistes rothneyi Cameron, was easy to confuse with V. bicolor in body color and size, which gave false information. There were 125 notifications of V. bicolor from 2016 to 2020, of which 82 (65.6%) notifications were invalid records. A total of 43 (34.4%) notifications were recorded as V. bicolor (Figure 3A), and most of them were individuals. Among them, 12 nests were located and removed (Figure 3B). The number and proportion of these notifications by city and county included 2 in Taoyuan (4.6%), 11 in Hsinchu (25.6%), and 30 in Miaoli (69.8%). We also confirmed that 6 of 13 apiaries in the city of Taoyuan and the counties of Hsinchu and Miaoli had V. bicolor records (Table 2). Seven apiaries failed to confirm V. bicolor occurrences in our monitoring, where there were 4 apiaries in the city of Taichung and the county of Nantou. Based on beekeeper witnesses, there was 3 apiaries in the counties of Hsinchu and Miaoli and the city of Taichung (Table 2). The years of these findings by city and county included 2020 (Taoyuan), 2016–2020 (Hsinchu), 2011–2020 (Miaoli), and 2019 (Taichung). Until 2020, the occurrences of V. bicolor confirmed by our research team were in the city of Taoyuan and the counties of Hsinchu and Miaoli.
There were 6 V. bicolor nests removed in Hsinchu County, 5 in Miaoli County, and 1 in Taoyuan City. There were 3 nests removed in Hsinchu and Miaoli from 2016 to 2019. However, in 2020, the number of removals reached 9 nests in Taoyuan, Hsinchu, and Miaoli (Figure 3A), showing a rapidly increasing trend. The months for the removal were mainly from July to November, but 2 nests were discovered in January 2019 (Figure 3B). Among the nest sites, 7 of 12 (58%) nests were built in the tree holes of Acacia confusa Merr., Cinnamomum camphora (L.) J. Presl, Cunninghamia sp., and Lagerstroemia subcostata Koehne. There were 3 of 12 (25%) nests that were built in man-made places, such as discarded tires, abandoned water towers, and wooden loudspeakers. There were 2 of 12 (17%) nests that were built in underground caves formed by the root spaces of plants (Figure 3C). Most of the nests in tree holes and underground caves cannot be taken out completely, and can only be removed with insecticides or by blocking the entrance of the nest. The adults in 9 nests were caught and counted as much as possible, which contained 29 to 1030 individuals.
Based on visual surveys, seasonal activities of V. bicolor were found from July to October at the Tongluo apiary in 2019 and 2020 (Figure 3D). The peak activities in the Tongluo and Beipu apiaries were from September to October in 2020, and there were nine active surveys during sunny, cloudy, and rainy days, with the exception of a single day due to rain conditions at Beipu on 22 October 2020 (Figure 3E). The number of V. bicolor hornets that attacked honeybees ranged from 1 to 28 per hour; in addition, it was found that the marked V. bicolor was captured in the second investigation at the Tongluo apiary (Figure 3E). The mean hourly numbers and standard deviations of V. bicolor hunters were about 9.8 ± 8.1 in Tongluo and 4.4 ± 2.9 in Beipu.

3.3. Potential Range for V. bicolor

The 24 initial variables were reduced to 8 after applying the standard of VIF values >10. The obtained values of TSS and AUC for five modelling algorithms indicated good predictive performance (TSS, >0.71 and AUC, >0.83), which can be considered an accurate validation. The results showed that the selected variables described the distribution of V. bicolor well. Among the eight environmental variables (Figure S1), temperature seasonality (BIO4, 48.8% contribution) and temperature of driest quarter (BIO9, 12.8% contribution) were the two most important predictors determining the distribution (Figure S2).
The total amount of pixels with maximum occurrence probability (suitable habitats) covered 440 km2. This area was mainly distributed to cover 232 (52.0%) and 168 (38.0%) km2 in the counties of Miaoli and Hsinchu (Figure 4). In addition, these areas were distributed in hilly or shallow mountains below 1000 m above sea level (Figure 4).

4. Discussion

Biological invasions often cause improper affection to the natural ecological environment, such as the decline of biodiversity or deterioration of ecosystem functions [45,46,47], threat to the native flora [48], and threat to agricultural and urban environments [49,50,51], and the levels involved can also harm the human and social systems [52]. Additionally, they posed a risk to aviation safety in Australia recently [53]. Due to multiple factors, hornet has become an invasive alien species (IAS), not only causing harm to the ecological environment, but also bringing a lot of challenges to human society, economy, and agricultural production safety. Since hornets have no nesting preference, they often establish colonies in anthropogenic areas, posing a potential threat to humans and other animals. Every year in Taiwan, the number of treated patients who experience Vespidae stings has been greater than those who experience Formicidae stings, causing life hazards [54]. To prevent this spread of the wasp species, government officials often adopt by incentivizing elimination programs, or hiring hornet hunters to collect nests for hornet liquor or larva dishes [55]. However, hornets are omnivores that can use various animal body or botanical products as their food, such as honeybees, rotten fruits, and plant nectars [14]. They hunt both beneficial insect pollinators and pest species, such as those in the orders Diptera and Hemiptera [14]. They have also been found to be a biological control for Polistes wasps [56]. We considered these elimination actions to be contentious due to the potential for negative ecological impact or harm to those proximal to the removal effort.
Although there were few samples out of Taiwan in this study, according to the phylogenetic analysis, the V. bicolor collected from different regions of Taiwan have high similarity. Only the samples from Southeastern China had high similarity with those from Taiwan. The database of this species in GenBank was not complete, so we can infer from the current results that the invasion of V. bicolor in Taiwan may have resulted from a single event. The earliest record of this species in Taiwan can be traced back to 2003 in Taichung City according to personal observations by J. T. Chao [15]. In previous studies concerning the Taiwanese hornet fauna, V. bicolor was never mentioned [57,58]; thus we reasonably confirmed our previous findings [15].
It was found that, until 2020, the population of V. bicolor was quite stable in the counties of Hsinchu and Miaoli. Miaoli is a county where V. bicolor has steadily occurred in all discovered regions, and monitoring has been done since 2011. Since 2016, Hsinchu County has noticed V. bicolor activities in apiaries, and nests have been occasionally found. In 2019 and 2020, there was evidence of nest activities and dispersal to Taoyuan City, and a suspected case was found in Taichung. The result shows that the trend of occurrence of V. bicolor is developing towards the northern or central areas of Taiwan. In Hong Kong, V. bicolor is widely distributed, and the activity period and colony cycle are extremely long and very tolerant of cold weather [17]. Moreover, nests were built in all kinds of environments, including in underground crevices and in and around buildings, all sorts of assorted objects, small trees, and shrubs [17]. These situations are similar to our findings.
During the peak activity in a suitable weather, the number of V. bicolor preying on honeybees was up to 28 at the Tongluo apiary per hour. The average number ranged from 4.4 to 9.8 per hour at the Tongluo and Beipu apiaries. An individual V. bicolor hornet was reported to prey on honeybees twice a day in the peak activity season. Previous studies found that the activity range of V. velutina depended on resources, but could expand up to 1 km from the nest [59]. Here we find that food resources and altitude are the two most important environmental factors determining the location of V. bicolor nests. The predation behavior of V. bicolor is similar to that of V. velutina [15]. The number of V. bicolor individuals was between 45 and 98 per day (for a 10 h activity calculation). Apiaries were reported to have up to 350 V. velutina individuals per day [60], indicating that V. bicolor has not been an important pest to honeybees until now. Therefore, the relevant data could provide some reference. This study provides a preliminary understanding of the situation of V. bicolor preying on honeybees.
With respect to ENM, the obtained AUC and TSS values were all larger than 0.7, which indicates that the predicted results were reliable [45,61]. We assessed the contributions of various environmental variables in the prediction of the V. bicolor distribution, and results indicate that temperature seasonality and climate were the key factors of species spatial distributions. Previous studies also found that the effects of temperature on hornets were significantly greater than those of other environmental factors [62]. The ecological niche of V. bicolor can greatly aid our understanding of changes to their geographical ranges. According to the result of the simulation, V. bicolor is distributed in the hilly or shallow mountain areas in Hsinchu and Miaoli counties. These areas, potentially due to having a suitable climate, provide nesting spaces and food resources for V. bicolor.
Our study found that there is not much genetic variation in the known existing populations. In fact, in the known positions of the existing records, the shortest distance between two points is not far away (not more than 15 km). Based on the estimation of the longest foraging distance of V. velutina, it is 11 km [63]. We suppose that there is no significant genetic variation in the existing V. bicolor populations. The reason is either there is a single source of invasion or there is no geographical isolation or both. In the future, more information can be learned from a large population of nuclear microsatellite DNA analysis [26]. Relevant information will be carried over to future experiments.

5. Conclusions

In the prediction of the dispersal model and monitoring survey, it was found that Hsinchu and Miaoli had stable V. bicolor populations due to the climate condition in these counties. According to the prediction results, the suitable habitat for V. bicolor is about 440 km2. Additionally, many local populations spread to the northern area of Taoyuan City; thus we assess that V. bicolor adapts to cooler weather in subtropical areas. In addition, the geographical distribution range was from hilly to shallow mountain areas of Taoyuan, Hsinchu, and Miaoli, potentially due to increased spaces and resources for nesting. Our results show that the V. bicolor population in Taiwan might have derived from the same invasion source. The predation of V. bicolor was mainly on honeybees, and the range of forage activities was highly connected to the presence of apiaries. ENM techniques allowed the spread of the invasive species after its recent arrival to be predicted with a good level of accuracy; hence, we recommend focusing on control efforts in these areas. Although the continuous removal of the IAS is often necessary, the management of V. bicolor should also involve considering it a local pest to reduce loss by beekeepers and public fear. Our findings highlight that we should be aware of the V. bicolor’s rapid range expansion in Taiwan. The value of monitoring and controlling its populations brings to light further control plan improvement.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12040320/s1: Figure S1: Sequence variation of seven haplotypes of V. bicolor in the COI region, Figure S2: Variables’ contribution for modeling, Table S1: List of environmental variables used in the model development.

Author Contributions

Designed the research, I-H.S., M.-L.L., and S.-S.L.; conducted field samplings and experiments, I-H.S., S.-S.L., Y.-J.L., and W.-C.Y.; provided samples other than Taiwan, J.T.; analyzed data, I-H.S., J.-Y.H., M.-L.L., Y.-J.L., and J.T.; wrote the draft manuscript, I-H.S., M.-L.L., S.-S.L., and Y.-J.L.; wrote the final version of the manuscript, I-H.S., J.T., and S.-S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Council of Agriculture, Executive Yuan, Taiwan, grant numbers 108 AS-1.2.1-S-aF and 109 AS-1.1.5-S-aM.

Data Availability Statement

All data presented in this study are available in the article.

Acknowledgments

We are grateful to Pei-Shou Hsu and Tzu-HsienWu at the Miaoli District Agricultural Research and Extension Station, Council of Agriculture, Executive Yuan, for their guidance on beekeeper information. We thank the beekeepers of Chin-Chong Hou and Yi-Chen Chen for providing experiment platforms. We are grateful to Chia-Hui Chao and Pin-Chen Yu at National Chiayi University, who assisted in the field collection. We are grateful to Ryo Inoue at Kyoto Sangyo University, Japan, for assisting with the test analysis. We are grateful to Erin Krichilsky for the first reading of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. EU. European Parliament and of the Council on the Prevention and Management of the Introduction and Spread of Invasive Alien Species. Regulation (EU) No. 1143. 2014. Available online: https://eur-lex.europa.eu/eli/reg/2014/1143/oj (accessed on 14 January 2021).
  2. SSC (Species Survival Commission). IUCN Guidelines of the Prevention of Biodiversity Loss Caused by Alien Invasive Species. Gland Switzerland. 2000. Available online: https://portals.iucn.org/library/efiles/documents/Rep-2000-052.pdf (accessed on 14 January 2021).
  3. Kenis, M.; Auger-Rozenberg, M.A.; Roques, A.; Timms, L.; Péré, C.; Cock, M.J.W.; Settele, J.; Augustin, S.; Lopez-Vaamonde, C. Ecological effects of invasive alien insects. In Ecological Impacts of Non-Native Invertebrates and Fungi on Terrestrial Ecosystems; Langor, D.W., Sweeney, J., Eds.; Springer: Dordrecht, The Netherlands, 2008; Volume 11, pp. 21–45. [Google Scholar] [CrossRef]
  4. Matsuura, M. Social Wasps of Japan in Color; Hokkaido University Press: Sapporo, Japan, 1995; p. 353. [Google Scholar]
  5. Haxaire, J.; Tamisier, J.P.; Bouguet, J.M. Vespa velutina Lepeletier, 1836, une redoutable nouveauté pour la faune de France (Hym., Vespidae). Bull. Soc. Entomol. Fr. 2006, 111, 194. [Google Scholar]
  6. Ueno, T. Establishment of the invasive hornet Vespa velutina (Hymenoptera: Vespidae) in Japan. Int. J. Chem. Environ. Biol. Sci. 2014, 2, 3. [Google Scholar]
  7. Kim, J.K.; Choi, M.B.; Moon, T.Y. Occurrence of Vespa velutina Lepeletier from Korea, and a revised key for Korean Vespa species (Hymenoptera: Vespidae). Entomol. Res. 2006, 36, 112–115. [Google Scholar] [CrossRef]
  8. Budge, G.E.; Hodgetts, J.; Jones, E.P.; Ostojá-Starzewski, J.C.; Hall, J.; Tomkies, V.; Semmence, N.; Brown, M.; Wakefield, M.; Stainton, K. The invasion, provenance and diversity of Vespa velutina Lepeletier (Hymenoptera: Vespidae) in Great Britain. PLoS ONE 2017, 12, e0185172. [Google Scholar] [CrossRef] [Green Version]
  9. Carvalho, J.; Hipólito, D.; Santarém, F.; Martins, R.; Gomes, A.; Carmo, P.; Rodrigues, R.; Grosso-Silva, J.; Fonseca, C. Patterns of Vespa velutina invasion in Portugal using crowdsourced data. Insect Conserv. Diver. 2020, 13, 501–507. [Google Scholar] [CrossRef]
  10. Takeuchi, T.; Takahashi, R.; Kiyoshi, T.; Nakamura, M.; Minoshima, Y.N.; Takahashi, J. The origin and genetic diversity of the yellow-legged hornet, Vespa velutina introduced in Japan. Insectes Soc. 2017, 64, 313–320. [Google Scholar] [CrossRef]
  11. Arca, M.; Mougel, F.; Guillemaud, T.; Dupas, S.; Rome, Q.; Perrard, A.; Muller, F.; Fossoud, A.; Capdevielle-Dulac, C.; Torres-Leguizamon, M.; et al. Reconstructing the invasion and the demographic history of the yellow-legged hornet, Vespa velutina, in Europe. Biol. Invasions 2015, 17, 2357–2371. [Google Scholar] [CrossRef]
  12. Wilson, T.M.; Takahashi, J.; Spichiger, S.E.; Kim, I.; van Westendorp, P. First reports of Vespa mandarinia (Hymenoptera: Vespidae) in North America represent two separate maternal lineages in Washington State, United States, and British Columbia, Canada. Ann. Entomol. Soc. Am. 2020, 20, 1–5. [Google Scholar] [CrossRef]
  13. Perrard, A.; Pickett, A.M.; Villemant, C.; Kojima, J.; CarpenterJ, J.M. Phylogeny of hornets: A total evidence approach (Hymenoptera, Vespidae, Vespinae, Vespa). J. Hymenopt. Res. 2013, 32, 1–15. [Google Scholar] [CrossRef]
  14. Barthélémy, C. A Provisional Identification Guide to the Social Vespids of Hong Kong (Hymenoptera: Vespidae). Available online: http://insectahk.com (accessed on 14 January 2021).
  15. Sung, I.H.; Lu, S.S.; Chao, J.T.; Yeh, W.C.; Lee, W.J. Establishment of Vespa bicolor in Taiwan (Hymenoptera, Vespidae). J. Insect Sci. 2014, 14, 231. [Google Scholar] [CrossRef] [Green Version]
  16. Chen, Y.; Tong, X. Biological habit of Vespa bicolor. J. Jishou Univ. Nat. Sci. 2004, 25, 80–84. [Google Scholar]
  17. Lee, J.X.Q. Potentially Dangerous Bees and Wasps of Hong Kong; Hong Kong Entomology Society: Hong Kong, China, 2009; pp. 50–53. [Google Scholar]
  18. Carpenter, J.M.; Kojima, J. Checklist of the species in the subfamily Vespinae (Insecta: Hymenoptera: Vespidae). Nat. Hist. Bull. Ibaraki Univ. 1997, 1, 51–92. [Google Scholar]
  19. Choi, M.B.; Kwon, O. Occurrence of Hymenoptera (Wasps and bees) and their foraging in the Southwestern part of Jirisan National Park, South Korea. J. Ecol. Environ. 2015, 38, 367–374. [Google Scholar] [CrossRef] [Green Version]
  20. Chantawannakul, P.; Guzman, L.D.; Li, J.; Williams, G.R. Parasites, pathogens, and pests of honeybees in Asia. Apidologie 2016, 47, 301–324. [Google Scholar] [CrossRef]
  21. Lee, V.J. The Oriental Honey-Buzzard of Ninety-Nine Peaks. Forestry Bureau, Council of Agriculture, Executive Yuan, Taiwan. Available online: https://youtu.be/KUeN6U2OnXI (accessed on 14 January 2021).
  22. Huang, K.Y.; Lin, Y.S.; Severinghaus, L.L. Nest provisioning of the Oriental Honey-buzzard (Pernis ptilorhyncus) in northern Taiwan. J. Raptor Res. 2004, 38, 367–371. [Google Scholar]
  23. Macià, F.X.; Menchetti, M.; Corbella, C.; Grajera, J.; Vila, R. Exploitation of the invasive Asian Hornet Vespa velutina by the European Honey Buzzard Pernis apivorus. Bird Study 2019, 66, 425–429. [Google Scholar] [CrossRef]
  24. With, K.A. The landscape ecology of invasive spread. Conserv. Biol. 2002, 16, 1192–1203. [Google Scholar] [CrossRef] [Green Version]
  25. Kamenova, S.; Bartley, T.; Bohan, D.A.; Boutain, J.R.; Colautti, R.I.; Domaizon, I.; Fontaine, C.; Lemainque, A.; Viol, I.L.; Mollot, G.; et al. Invasions toolkit: Current methods for tracking the spread and impact of invasive species. Adv. Ecol. Res. 2017, 56, 85–182. [Google Scholar]
  26. Choi, M.B.; Lee, S.A.; Suk, H.Y.; Lee, J.W. Microsatellite variation in colonizing populations of yellow-legged Asian hornet, Vespa velutina nigrithorax, in South Korea. Entomol. Res. 2013, 43, 208–214. [Google Scholar] [CrossRef]
  27. Villemant, C.; Barbet-Massin, M.; Perrard, A.; Muller, F.; Gargominy, O.; Jiguet, F.; Rome, Q. Predicting the invasion risk by the alien bee-hawking Yellow-legged hornet Vespa velutina nigrithorax across Europe and other continents with niche models. Biol. Conserv. 2011, 144, 2142–2150. [Google Scholar] [CrossRef]
  28. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar]
  29. Huang, W.; Xie, X.; Huo, L.; Liang, X.; Wang, X.; Chen, X. An integrative DNA barcoding framework of ladybird beetles (Coleoptera: Coccinellidae). Sci. Rep. 2020, 10, 10063. [Google Scholar] [CrossRef] [PubMed]
  30. van Houdt, J.K.; Breman, F.C.; Virgilio, M.; DeMeyer, M. Recovering full DNA barcodes from natural history collections of Tephritid fruitflies (Tephritidae, Diptera) using mini barcodes. Mol. Ecol. Resour. 2010, 10, 459–465. [Google Scholar] [CrossRef] [PubMed]
  31. Gibson, J.F.; Kelso, S.; Jackson, M.D.; Kits, J.H.; Miranda, G.F.G.; Skevington, J.H. Diptera-specific polymerase chain reaction amplification primers of use in molecular phylogenetic research. Ann. Entomol. Soc. Am. 2011, 104, 976–997. [Google Scholar] [CrossRef] [Green Version]
  32. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  33. Wei, S.J.; Niu, F.F.; Tan, J.L. The mitochondrial genome of the Vespa bicolor Fabricius (Hymenoptera: Vespidae: Vespinae). Mitochondrial DNA A DNA Mapp. Seq. Anal. 2016, 27, 875–876. [Google Scholar] [CrossRef]
  34. Clark, K.; Karsch-Mizrachi, I.; Lipman, D.J.; Ostell, J.; Sayers, E.W. GenBank. Nucleic Acids Res. 2016, 44, 67–72. [Google Scholar] [CrossRef] [Green Version]
  35. Rozas, J.; Ferrer-Mata, A.; Sánchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sánchez-Gracia, A. DnaSP v6: DNA sequence polymorphism analysis of large datasets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  36. Bandelt, H.J.; Forster, P.; Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Ecol. 1999, 16, 37–48. [Google Scholar] [CrossRef]
  37. An, J. Dances with Hornets; Independent Writer Press: Taipei, Taiwan, 2015; pp. 155–163. [Google Scholar]
  38. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  39. Evans, J.S.; Oakleaf, J.; Cushman, S.A.; Theobald, D. An ArcGIS Toolbox for Surface Gradient and Geomorphometric Modeling, Version 2.0-0. 2014. Available online: http://evansmurphy.wix.com/evansspatial (accessed on 14 January 2021).
  40. Naimi, B.; Araújo, M.B. Sdm: A reproducible and extensible R platform for species distribution modelling. Ecography 2016, 39, 368–375. [Google Scholar] [CrossRef] [Green Version]
  41. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  42. Hao, T.; Elith, J.; Guillera-Arroita, G.; Lahoz-Monfort, J.J. A review of evidence about use and performance of species distribution modelling ensembles like BIOMOD. Divers. Distrib. 2019, 25, 839–852. [Google Scholar] [CrossRef]
  43. Pearce, J.; Ferrier, S. Evaluating the predictive performance of habitat models developed using logistic regression. Ecol. Model. 2000, 133, 225–245. [Google Scholar] [CrossRef] [Green Version]
  44. Giannini, T.C.; Acosta, A.L.; Garófalo, C.A.; Saraiva, A.M.; Alves-dos-Santos, I.; Imperatriz-Fonseca, V.L. Pollination services at risk: Bee habitats will decrease owing to climate change in Brazil. Ecol. Model. 2012, 244, 127–131. [Google Scholar] [CrossRef]
  45. Pyšek, P.; Richardson, D.M. Invasive species, environmental change and management, and health. Annu. Rev. Environ. Resour. 2010, 35, 25–55. [Google Scholar] [CrossRef] [Green Version]
  46. Orlova-Bienkowskaja, M.J.; Bieńkowski, A.O. Alien Coccinellidae (Ladybirds) in Sochi National Park and its vicinity, Russia. Nat. Conserv. Res. 2017, 2, 96–101. [Google Scholar] [CrossRef] [Green Version]
  47. Ruchin, A.B.; Egorov, L.V.; Lobachev, E.A.; Lukiyanov, S.V.; Sazhnev, A.S.; Semishin, G.B. Expansion of Harmonia axyridis (Pallas, 1773) (Coleoptera: Coccinellidae) to European part of Russia in 2018–2020. Baltic J. Coleopterol. 2020, 20, 51–60. [Google Scholar]
  48. Dubovik, D.V.; Skuratovich, A.N.; Miller, D.; Spiridovich, E.V.; Gorbunov, Y.N.; Vinogradova, Y.K. The invasiveness of Solidago canadensis in the Sanctuary «Prilepsky» (Belarus). Nat. Conserv. Res. 2019, 4, 48–56. [Google Scholar] [CrossRef]
  49. Pimentel, D.; Zuniga, R.; Morrison, D. Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecol. Econ. 2005, 52, 273–288. [Google Scholar] [CrossRef]
  50. Milosavljević, I.; El-Shafie, H.A.; Faleiro, J.R.; Hoddle, C.D.; Lewis, M.; Hoddle, M.S. Palmageddon: The wasting of ornamental palms by invasive palm weevils, Rhynchophorus spp. J. Pest Sci. 2019, 92, 143–156. [Google Scholar] [CrossRef]
  51. Gaertner, M.; Wilson, J.R.U.; Cadotte, M.W.; MacIvor, J.S.; Zenni, R.D.; Richardson, D.M. Non-native species in urban environments: Patterns, processes, impacts and challenges. Biol. Invasions 2017, 19, 3461–3469. [Google Scholar] [CrossRef] [Green Version]
  52. Villamagna, A.M.; Murphy, B.R. Ecological and socio-economic impacts of invasive water hyacinth (Eichhornia crassipes): A review. Freshw. Biol. 2010, 55, 282–298. [Google Scholar] [CrossRef]
  53. House, A.P.N.; Ring, J.G.; Shaw, P.P. Inventive nesting behaviour in the keyhole wasp Pachodynerus nasidens Latreille (Hymenoptera: Vespidae) in Australia, and the risk to aviation safety. PLoS ONE 2020, 15, e0242063. [Google Scholar] [CrossRef]
  54. Lin, Y.Y.; Chiu, C.C.; Chang, H.S.; Kao, Y.H.; Hsiao, P.J.; Chuu, C.P. Comparison of clinical manifestations, treatments, and outcomes between Vespidae sting and Formicidae sting patients in the emergency department in Taiwan. Int. J. Environ. Res. Public Health 2020, 17, 6162. [Google Scholar] [CrossRef]
  55. Matsuura, M. Edible Wasps—Insect Food Culture Visitation; Hokkaido University Press: Sapporo, Japan, 2002; p. 356. [Google Scholar]
  56. Hozumi, S.; Yamane, S. Incubation ability of the functional envelope in paper wasp nests (Hymenopteta, Vespidae, Polistes): I. field measurements of nest temperature using paper models. J. Ethol. 2001, 19, 39–46. [Google Scholar] [CrossRef]
  57. Starr, C.K. The social wasps (Hymenoptera: Vespidae) of Taiwan. Bull. Nat. Mus. Nat. Sci. 1992, 3, 93–138. [Google Scholar]
  58. Kojima, J.I.; Saito, F.; Nguyen, L.T.P. On the species-group taxa of Taiwanese social wasps (Hymenoptera: Vespidae) described and/or treated by J. Sonan. Zootaxa 2011, 2920, 42–64. [Google Scholar] [CrossRef] [Green Version]
  59. Poidatz, J.; Monceau, K.; Bonnard, O.; Thiéry, D. Activity rhythm and action range of workers of the invasive hornet predator of honeybees Vespa velutina, measured by radio frequency identification tags. Ecol. Evol. 2018, 8, 7588–7598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Monceau, K.; Bonnard, O.; Moreau, J.; Thiery, D. Spatial distribution of Vespa velutina individuals hunting at domestic honeybee hives: Heterogeneity at a local scale. Insect Sci. 2014, 21, 765–774. [Google Scholar] [CrossRef]
  61. Lasram, F.B.R.; Guilhaumon, F.; Albouy, C.; Somot, S.; Thuiller, W.; Mouillot, D. The Mediterranean Sea as a ‘cul-de-sac’ for endemic fishes facing climate change. Glob. Change Biol. 2010, 16, 3233–3245. [Google Scholar] [CrossRef]
  62. Bessa, A.S.; Carvalho, J.; Gomes, A.; Santarém, F. Climate and land-use drivers of invasion: Predicting the expansion of Vespa velutina nigrithorax into the Iberian Peninsula. Insect Conserv. Diver. 2016, 9, 27–37. [Google Scholar] [CrossRef]
  63. Lioy, S.; Manino, A.; Porporato, M.; Laurino, D.; Romano, A.; Capello, M.; Bertolino, S. Establishing surveillance areas for tackling the invasion of Vespa velutina in outbreaks and over the border of its expanding range. NeoBiota 2019, 46, 51–69. [Google Scholar] [CrossRef]
Figure 1. Occurrence records of V. bicolor in Taiwan.
Figure 1. Occurrence records of V. bicolor in Taiwan.
Insects 12 00320 g001
Figure 2. Phylogenetic tree and haplotype network of V. bicolor from China, Hong Kong, Nepal, Taiwan, and Vietnam. (A) Phylogenetic tree computed on COI sequences. Each number on the branch represents node stats in maximum likelihood bootstrap values. (B) Haplotype network. Each bar on the branch represents a single nucleotide substitution. Circle size is equivalent to sample size. Median vectors are represented by black circles.
Figure 2. Phylogenetic tree and haplotype network of V. bicolor from China, Hong Kong, Nepal, Taiwan, and Vietnam. (A) Phylogenetic tree computed on COI sequences. Each number on the branch represents node stats in maximum likelihood bootstrap values. (B) Haplotype network. Each bar on the branch represents a single nucleotide substitution. Circle size is equivalent to sample size. Median vectors are represented by black circles.
Insects 12 00320 g002
Figure 3. Field survey of V. bicolor from Taiwan: (A) numbers of notifications and eliminated nests of V. bicolor; (B) months during which nests were eliminated; (C) nest site, number of nests, and proportion in V. bicolor; (D) months during which V. bicolor activities were observed; (E) daily activities of V. bicolor in the Tongluo and Beipu apiaries.
Figure 3. Field survey of V. bicolor from Taiwan: (A) numbers of notifications and eliminated nests of V. bicolor; (B) months during which nests were eliminated; (C) nest site, number of nests, and proportion in V. bicolor; (D) months during which V. bicolor activities were observed; (E) daily activities of V. bicolor in the Tongluo and Beipu apiaries.
Insects 12 00320 g003
Figure 4. Prediction of potential distribution area for V. bicolor in Taoyuan, Hsinchu, Miaoli, and Taichung areas.
Figure 4. Prediction of potential distribution area for V. bicolor in Taoyuan, Hsinchu, Miaoli, and Taichung areas.
Insects 12 00320 g004
Table 1. The sample of Vespa bicolor from Taiwan, with the collection site and accession number in GenBank.
Table 1. The sample of Vespa bicolor from Taiwan, with the collection site and accession number in GenBank.
Species
City/County
LocalityYear, DateCollection SiteGenBank Accession Number
Attacking Bees in ApiariesVisiting FlowersOwn Nests
Vespa bicolor
TaoyuanDaxi17 July 2020 MW455072
Fuxing29 August 2020 MW455074
HsinchuBeipu15 November 2019 MW455066
Beipu17 July 2020 MW455071
Guanxi22 July 2019 MW455065
Guanxi10 August 2020 MW455073
Hengshan22 October 2020 MW455078
Hengshan13 January 2016 MW455082
Hengshan13 January 2016 MW455083
Qionglin4 November 2018 MW455081
Wufeng9 September 2020 MW455076
Xinpu30 October 2020 MW455079
Zhudong3 January 2020 MW455068
MiaoliSanwan16 November 2017 MW455084
Sanyi3 January 2020 MW455067
Tongluo16 August 2019 MW455061
Tongluo2 August 2019 MW455062
Tongluo18 September 2019 MW455063
Tongluo16 July 2020 MW455069
Tongluo16 July 2020 MW455070
Tongluo19 August 2020 MW455075
Tongluo30 October 2020 MW455080
Touwu30 October 2020 MW455077
confirmed in this study.
Table 2. The field survey and occurrence records of V. bicolor in 13 apiaries.
Table 2. The field survey and occurrence records of V. bicolor in 13 apiaries.
City/CountyLocalitySurvey YearHornet TrapOccurrence
Taoyuan CityDaxi2020No
Fuxing2020No
Hsinchu CountyBeipu2019–2020Yes
Guanxi2019Yes
Hengshan2020No
Xinpu2020Yes
Miaoli CountyDahu2020No
Tongluo2019–2020Yes
Tongxiao2019Yes
Taichung CityHouli2020No
South District2020Yes
Wufeng2020Yes
Nantou CountyGuoxing2020No
○ confirmed in this study; — not confirmed in this study; △ beekeeper records.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, S.-S.; Takahashi, J.; Yeh, W.-C.; Lu, M.-L.; Huang, J.-Y.; Lin, Y.-J.; Sung, I.-H. Evidence for Range Expansion and Origins of an Invasive Hornet Vespa bicolor (Hymenoptera, Vespidae) in Taiwan, with Notes on Its Natural Status. Insects 2021, 12, 320. https://doi.org/10.3390/insects12040320

AMA Style

Lu S-S, Takahashi J, Yeh W-C, Lu M-L, Huang J-Y, Lin Y-J, Sung I-H. Evidence for Range Expansion and Origins of an Invasive Hornet Vespa bicolor (Hymenoptera, Vespidae) in Taiwan, with Notes on Its Natural Status. Insects. 2021; 12(4):320. https://doi.org/10.3390/insects12040320

Chicago/Turabian Style

Lu, Sheng-Shan, Junichi Takahashi, Wen-Chi Yeh, Ming-Lun Lu, Jing-Yi Huang, Yi-Jing Lin, and I-Hsin Sung. 2021. "Evidence for Range Expansion and Origins of an Invasive Hornet Vespa bicolor (Hymenoptera, Vespidae) in Taiwan, with Notes on Its Natural Status" Insects 12, no. 4: 320. https://doi.org/10.3390/insects12040320

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop