Research Article Journal of Orthoptera Research 2024, 33(2): 217-227 Analysis of potential niche shifts in alien pairs of mantis species (Insecta, Mantodea) with comments on the current taxonomic and ecological knowledge MattiA DE Vivo!2 1 Biodiversity Program, Taiwan International Graduate Program, Academia Sinica and National Taiwan Normal University, No. 28, Lane 70, Section 2, Yanjiuyuan Road, Nangang District, Taipei 115, Taipei, Taiwan. 2 Department of Biogeography, Trier University, Universitatsring 15, Trier D-54296, Trier, Germany. Corresponding author: Mattia De Vivo (mattiadevivopatalano@gmail.com) Academic editor: Matan Shelomi | Received 12 August 2023 | Accepted 5 February 2024 | Published 11 June 2024 https://zoobank.org/89E1F112-C575-4D75-9FD3-8081EDA7B8FD Citation: De Vivo M (2024) Analysis of potential niche shifts in alien pairs of mantis species (Insecta, Mantodea) with comments on the current taxonomic and ecological knowledge. Journal of Orthoptera Research 33(2): 217-227. https://doi.org/10.3897/jor.33.111057 Abstract Due to the pet and goods trade, several animals are now present in regions outside of their traditional native ranges. A peculiar situation has arisen in mantises, insects that are becoming more popular as pets: two genera (Hierodula and Tenodera) have begun to spread around the world, with two Hierodula species overlapping in Europe and two Tenodera species doing the same in North America. Such an event can lead to possible competition with both local taxa and alien congeneric sister species; the latter may reduce the impact of one of the invaders. Additionally, the situation allows the comparisons of niche shifts in displaced mantises, allowing us to understand whether such animals respect general patterns shown in terrestrial ectothermic invasive species. To do this, I adapted scripts from previous publications for analyzing niche overlap (Schoener’s D), niche expansion (E), and unfilling (U) through the centroid shift, overlap, unfilling, and expansion (COUE) scheme using presence records from GBIF and iNaturalist Research-Grade observations and bioclimatic variables available in BIOCLIM, selected according to variance inflation factor (VIF) values. I also evaluated the overlap between the sister species in the non-native range with D. Overall, there was relatively high niche expansion and unfilling patterns shared among the taxa, although species tended to have low abiotic overlap between native and alien ranges, and a relatively high niche overlap was present among congeneric species in the shared non-native area. However, such analyses may be biased due to chosen variables, taxonomic uncertainty, and lack of information on mantises’ ecology; particularly, the situation regarding H. tenuidentata/transcaucasica should be monitored and clarified, given the higher potential invasion risk of these species compared to other mantises and the uncertainties regarding which populations have reached Europe. Additionally, the biology of alien mantises should be studied in more detail in both native and non-native environments given the current critical lack of information. Keywords COUE, generalist predator, Hierodula, overlap, Tenodera Introduction The increase in global trade has allowed several species to be present outside their native ranges, leading to the presence of more than 37,000 alien species globally (IPBES 2023) and to an exponential increase of invasive taxa (Meurisse et al. 2019, Mormul et al. 2022). The establishment of an alien species starts with overcoming the dispersal barrier, which is aided by humans through accidental transportation or species traded for commercial reasons (Seebens et al. 2015, Meurisse et al. 2019, Sinclair et al. 2020, Mormul et al. 2022, IPBES 2023, Pili et al. 2023). For example, several animals are imported through the pet market, which is known to have a strong influence on their invasiveness (Lockwood et al. 2019, Maceda-Veiga et al. 2019, Gippet and Bertelsmeier 2021). Insects are usually introduced in non-native environments by accidental means (e.g., Meurisse et al. 2019). Among them, man- tises (Mantodea) were probably introduced in different parts of the world by both accidental transportation (Battiston et al. 2018, 2020, Moulin 2020, Connors et al. 2022) and pet mismanage- ment (Battiston et al. 2022, Connors et al. 2022), although the introduction pathway of each species may be different (Battiston et al. 2022). Non-native mantises may compete with native Man- todea taxa (e.g., Fea et al. 2013) and may also have secondary eco- system effects, such as feeding on pollinators (Anderson 2019), herbivores, and intra-guild species (Snyder and Evans 2006, Crow- der and Snyder 2010). At the same time, alien (i.e., non-native, e.g., Almena et al. 2023) mantises sharing non-native environments may also compete with each other (e.g., Snyder and Hurd 1995). Competition between introduced species has been shown to po- tentially reduce the impact of an invasive species (e.g., Van Riel et al. 2009, Russell et al. 2014). Copyright Mattia De Vivo. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 218 Although it is known that several mantises are now present in regions outside their known native ranges (e.g., Fea et al. 2013, Battiston et al. 2017, 2018, 2020, 2022, Schwarz and Ermann 2018, Anderson 2019, Shcherbakov and Govorov 2020, 2021, Connors et al. 2022, Martinovi¢ et al. 2022), to the best of my knowledge, it has never been analyzed if their alien niches are similar to their native ones or if these niches have shifted. Generally, such analyses are based on the centroid shift, overlap, unfilling, and expansion (COUE) scheme, which was first presented by Broennimann et al. (2012) and later modified by Petitpierre et al. (2012) (Liu et al. 2020). This scheme is based on first reducing the environmental space in two dimensions through principal component analysis (PCA) and then splitting said space into three parts: stability (S), which is the space that the species occupies in both the native and introduced range; unfilling (U), which is the space occupied in the native range and also present but not occupied in the alien one; and expansion (E), which is the area occupied only in the alien range and also present but not occupied in the native range (Broen- nimann et al. 2012, Liu et al. 2020). To avoid issues caused by dif- ferent sampling efforts and strategies, a kernel density function is used to smooth the occurrence densities in both native and intro- duced environments (Broennimann et al. 2012, Liu et al. 2020). The COUE scheme is important for the quantification of niche space and for understanding the niche shift of invasive species. E provides the niche shift magnitude, given that U shows a potential condition that species may colonize in the alien range after some time (Liu et al. 2020). Furthermore, the COUE scheme allows for the testing of niche conservatism as a binary pattern and checks if both native and alien niches are significantly similar to each other (niche equivalency test; Bates et al. 2020, Liu et al. 2020) and if such resemblance is caused by chance by randomly allocating the occurrences in one range (niche similarity test) (Liu et al. 2020). COUE analyses can be helpful for understanding the potential ecological patterns of alien species. For example, it is known that ectothermic terrestrial species considered as invasive may tend to conserve their climatic niche compared to those regarded as alien (introduced outside their native range) but non-invasive (Bates et al. 2020, Liu et al. 2020). This pattern is probably explained by the biotic, abiotic, movement (BAM) model introduced by Soberon and Peterson (2005); according to this model, species distribu- tion is constrained by factors linked to biotic (B) and abiotic (A) interactions, coupled with the need of the species to reach certain places (M) (Soberon and Peterson 2005). Therefore, species might need more overlapping of such factors in alien environments to become invasive (Bates et al. 2020, Mormul et al. 2022). There- fore, species that are introduced to climates similar to their native ones may be able to establish populations in such environments (Liu et al. 2020). At the same time, methods for analyzing niche shifts may also be biased by several factors. For example, the cho- sen metric thresholds and variables have been shown to hugely influence results and interpretations (Bates and Bertelsmeier 2021; Lo Parrino et al. 2023). It is also important to note that incom- plete ecological niches may be suitable for one population of a non-native species to colonize while excluding a second, closely related population (Liu et al. 2020). In fact, it is recognized that different populations or subspecies of the same species may have different ecological niches (e.g., Weaver et al. 2022) or plasticity to colonize different environments (Bates and Bertelsmeier 2021). A curious situation arose in recent years with at least two spe- cies of giant Asian mantises (genus Hierodula Burmeister, 1838) colonizing environments in Europe (Battiston et al. 2018, 2020, Moulin 2020, Martinovi¢ et al. 2022, Sevgili and Yilmaz 2022, M. DE VIVO Moulin and Rouard 2023) and two species of Tenodera Burmeister, 1838 colonizing North America, particularly Northeastern United States (Snyder and Hurd 1995, Snyder and Evans 2006, Anderson 2019). Such a situation may allow us not only to compare the pos- sibility of niche shifts in such insects but also to evaluate potential differences between congeneric taxa sharing a non-native range, which may lead to potential competition between the invaders (e.g., Snyder and Hurd 1995). It could be possible that the native conditions for species spreading more easily, such as one of the Hierodula species in Europe (Sevgili and Yilmaz 2022), are similar to those in its non-native one, while species with more difficulty in dispersal could exhibit less overlap if the generic pattern observed for ectothermic species are respected (Bates et al. 2020, Liu et al. 2020). In addition, different climatic requirements may lead the species to occupy different environments, as happened with other sister pairs in a different taxonomic group (Summers et al. 2023). Iran niche shift analyses on the two clades, Hierodula and Teno- dera, to evaluate possible patterns in alien Mantodea sister species. First, I considered the native and alien ranges of a single species and calculated overlap and niche shifts, and then I compared the alien ranges shared by sister taxa. Based on previous observations (Anderson 2019, Battiston et al. 2020, Sevgili and Yilmaz 2022), I expected that the most widespread alien species of each genus would have high overlap between their own alien and native niche, as happens in other ectothermic terrestrial animals (Bates et al. 2020, Liu et al. 2020), and that can lead to a competitive advantage, as has happened in other taxa. Such tasks were com- plicated and influenced by taxonomic (Battiston et al. 2018, 2021, Shcherbakov and Battiston 2020, Liu et al. 2021, Shim et al. 2021, van der Heyden and Schwarz 2021) and ecological uncertainties (Anderson 2019). Furthermore, the analyses relied on citizen sci- ence data, which are very helpful in tracking invasions (e.g., Bat- tiston et al. 2020, Moulin 2020, Connors et al. 2022, Martinovic et al. 2022, Sevgili and Yilmaz 2022, Moulin and Rouard 2023) and generating observational records (iNaturalist Blog 2023, GBIEorg 2023a, 2023b) but may have limitations in terms of identifying insects at the species level (Gardiner et al. 2012, Ratnieks et al. 2016). I briefly review these issues here given that they will influ- ence possible future studies. Materials and methods Considered species and caveats.—I conducted analyses on the fol- lowing taxa: the Transcaucasian Giant Mantis Hierodula transcau- casica Brunner de Wattenwyl, 1878 and/or the Giant Asian Mantis H. tenuidentata Saussure, 1869 (see sections below); the Indochina mantis H. patellifera Serville, 1839; the narrow-winged mantis Ten- odera angustipennis Saussure, 1869; and the Chinese mantis Teno- dera sinensis (Saussure, 1871). Species of the genus Hierodula are naturally distributed in the Oriental, Palaearctic, and Australian biogeographic regions but mostly in Asia (Leong 2009, Liu et al. 2021). H. tenuidentata, in particular, is thought to be distributed from Indonesia to Caucasus, with alien populations from Northern Ukraine to Spain (Pushkar and Kavurka 2016, Shcherbakov and Battiston 2020, Sevgili and Yilmaz 2022, Mirzaee et al. 2023, Pintilioaie et al. 2023). However, the presence of this species in Indonesia is doubtful (Shcherbakov and Battiston 2020). Battiston et al. (2018) officially synonymized H. transcaucasica with H. tenuidentata, given the high morphologi- cal similarities between the two. H. transcaucasica is generally de- fined as present in the Caucasus region from Iran to Georgia (Kol- negari 2023). While synonymy has been used in some works (e.g., JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) M. DE VIVO Shcherbakov and Battiston 2020, Sevgili and Yilmaz 2022, Mirzaee et al. 2023, Pintilioaie et al. 2023), it does not seem to be univer- sally accepted (e.g., van der Heyden and Schwarz 2021, Kharabadze et al. 2022, Kolnegari 2023, Moulin and Rouard 2023). In some instances, authors have stated that they were not able to disprove synonymy or split in such taxon (Liu et al. 2021, van der Heyden and Schwarz 2021). Given such uncertainties, | compared the nich- es between the two potential native ranges of said species/popula- tions and found that their ecological overlap is very low (D = 0.091, niche equivalency test’s p < 0.0001; Suppl. material 1: figs S1-S3). Therefore, I split the datasets and checked for potential niche shifts between the invading European populations and H. transcaucasica sensu stricto and H. tenuidentata sensu stricto or sensu Battiston et al. (2018). To perform this split, I checked the presence records and di- vided H. transcaucasica and H. tenuidentata sensu stricto according to the lack of records in some Central Asian areas, which was already noted by Pintilioaie et al. (2023) (Fig. 1). The split was also used for the overlap analysis described above. Following Shcherbakov and Battiston (2020), I considered the Balkan population as alien and the Crimean one as native, although their status (i.e., whether they are a product of a natural expansion or artificial introduc- tion) remains uncertain (Pushkar and Kavurka 2016, Battiston et al. 2018, Shcherbakov and Battiston 2020, Pintilioaie et al. 2023). It is known that H. patellifera actually represents a species complex distributed in several South, South-East, and East Asia regions (Liu et al. 2021, Chih-Ting Hsu pers. comm.). It also pre- sents high intra-specific morphological variation, even in terms of the white coaxal marks that are usually used for species identifica- tion (Leong 2009, Battiston et al. 2020, Oshima 2021, Shim et al. 2021). Therefore, the limits of the taxon are unclear. After discuss- ing with a taxonomic expert (Chih-Ting Hsu), I decided to run the intra-species analyses with two datasets: one in which I removed potentially problematic localities (India, Nepal, Philippines) and another in which I kept the data from these countries (Suppl. ma- terial 1: fig. S4). It must be noted that the niche similarity between the India, Nepal, and Philippines datasets and the one with the rest of the native area was low (D = 0. 294, p = 0; Suppl. material 1: figs S5-S7). It should also be noted that I considered points from areas outside of Europe (e.g., Seychelles and Hawaii) for the intra- species comparisons in this taxon, given that such observations seem to represent H. patellifera sensu stricto (Nguyen and Maxwell 50 - 40- °® ® t + 30- & 20 - 10 - 0 20 40 60 long 219 2008, Seychelles Nation 2009, Chih-Ting Hsu pers. comm.). In addition, given that the population from Christmas Island is un- likely to represent H. patellifera sensu stricto (Connors 2023, Chih- Ting Hsu pers. comm.), those data were not included in my analy- ses (Suppl. material 1: fig. $4). Furthermore, given that the status of H. patellifera in Korea is unclear (Taewoo Kim, pers. comm.), I considered the Korean population as native, although this species was not recorded in the region before 1999 (Shin et al. 2023). The species of the genus Tenodera are naturally distributed in the eastern hemisphere in various biogeographical realms from Africa to Australasia (Jensen et al. 2009; Anderson 2019). T. sin- ensis was introduced in North America around 1896 and T. an- gustipennis around 1926 in the same region (Snyder and Hurd 1995, Anderson 2019). Both taxa are naturally distributed in East Asia, with T: sinensis being more spread in both the native and non-native areas (Anderson 2019, GBIEorg 2023b) and known to be more biologically competitive in North America (Hurd 1988, Snyder and Hurd 1995, Anderson 2019). It is unclear if these spe- cies are capable of establishing populations outside of the North- eastern area of the United States, although some individuals can overcome such geographical barrier (Anderson 2019, iNaturalist 2023). Therefore, I analyzed the data first with all the available points and again omitting the points according to the distribu- tion maps made by Anderson (2019), which considered the conti- nental United States and Canada (Suppl. material 1: figs S8-S11). Taxonomic issues are also present within this genus, given that mostly citizen-science observations (e.g., iNaturalist records) were used and that such mantises can be morphologically similar to each other with few characters available for distinguishing them Jensen et al. 2009, Anderson 2019). Given the considerations in the paragraphs above, identifica- tion of such animals based on citizen science might be biased. However, since 2020, iNaturalist has provided at least 64% of the insect data present in GBIF (iNaturalist Blog 2023), which is one of the most used archives for biological data (Heberling et al. 2021). Conversely, if all the Mantodea data with coordinates and uncertainty less than or equal to 10 km (limits set given the specifi- cations used for “spThin”; see sections below) are considered, iN- aturalist Research-Grade observations are responsible for approxi- mately 81,647 observational data out of the 111,728 available in GBIF (GBIEorg 2023a), accounting for roughly 73.08% of such ee Se ; Supposed population & S 3 ® H. tenuidentata sensu stricto - OP pants af e Ali e 8e ien range eT) H. transcaucasica wi" e @ od 80 Fig. 1. Presence records used in the analyses for H. tenuidentata sensu Battiston et al. (2018). The maps for the other species are available in Suppl. material 1: figs $4, S8-S11. JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 220 data in the archive. Critically, such a percentage increases to ap- proximately 98.1% if only the evaluated taxa are taken into con- sideration (GBIEorg 2023b). Therefore, citizen-science records might currently be regarded as the best method for obtaining as many observation points for mantises as possible. Presence records and environmental variables.—Presence records from the focal species were downloaded from Research-Grade iNatural- ist (iNaturalist 2023) observations, which means that at least two- thirds of the users who tried to identify the animal agreed on the species ID (https://www.inaturalist.org/pages/help#quality), and from GBIF (GBIF.org 2023b). Such records were then filtered using the R package “spThin” (version 0.2.0: Aiello-Lammens et al. 2015) to reduce potential sampling bias and also remove duplicated co- ordinates, with each point set at least 10 km apart, one maximum possible output file, and 100 repetitions for getting an optimal number of presence points. The records were then converted to spa- tial points with the “sp” package (version 1.4-7; Bivand et al. 2013). Potential erroneous points, points representing potential single es- caped pets, or records of single instances of species in regions with no literature available (e.g., a single point of H. patellifera in Cali- fornia: https://www.inaturalist.org/observations/93123522) were removed from the datasets. In the case of doubtful records (e.g., the ones for Oceanian Hierodula), the possible species ID was checked in collaboration with a mantis taxonomist (Chih-Ting Hsu). For the analyses regarding the statuses of H. tenuidentata, I also added Chinese observation points from Liu et al. (2021). For building the climatic niche of the species, I downloaded the 19 bioclimatic variables from WorldClim2 (Fick and Hijmans 2017) at 30s (~1 km’) resolution. Then, the value of each climate variable for each spatial point was extracted with the extract() function of the “raster” package (version 3.5-15; Hijmans 2022). I then removed all the records with at least one missing environ- mental variable with the na.omit() function of the “stats” pack- age (version 4.1.1; R Core Team 2022). After that, I merged all the values of the variables together with the rbind() function of the “base” package (version 4.1.1; R Core Team 2022) and sub- setted them by removing collinear ones; collinearity was evalu- ated by variance inflation factor (VIF) matrices, with 5 as a thresh- old, using the vifstep() function implemented in the “usdm” R package (version 2.1-6; Naimi et al. 2014), and I chose to use all the metrics with VIF < 5 (Table 2). It must be noted that the cho- sen variables can influence the results of the analyses (Bates and Bertelsmeier 2021, Lo Parrino et al. 2023). However, the major- ity of the variables used here were also used for niche modeling of mantises in past publications (e.g., Steger et al. 2020, De Vivo and Huang 2022, Pintilioaie et al. 2023, Shin et al. 2023), and they should also be relevant for the studied taxa. For example, it is known that temperature can influence the biology of H. patellifera (Shin et al. 2023) and Tenodera species (Snyder and Hurd 1995, Hurd et al. 2020) and therefore including a variable such as BIO5 (max temperature of warmest month) is biologically meaningful. Potential spatial autocorrelation was evaluated by Mantel correlo- grams through the ecospat.mantel.correlogram() function of the ecospat” package (version 3.2.2; Di Cola et al. 2017). Niche shift and overlap analyses.—To identify potential niche shifts and congeneric competition, while also checking on potential pat- terns shared by those taxa, I adapted and used scripts from Bates et al. (2020) and ran the analyses in R (version 4.1.3; R Core Team 2022). The variables were reduced through the use of a principal component analysis (PCA) with functions from the “ade4” pack- M. DE VIVO age (version 1.7-19; Dray and Dufour 2007). Then, I ran a be- tween-class analysis with 10 axes from said PCA using a priori classes (i.e., native vs non-native for all the taxa or H. tenuidentata sensu stricto vs H. transcaucasica sensu stricto) to identify which axis separated the two categorical ranges the most for each of the con- sidered groups (i.e., the same species in two different ranges or two species sharing the alien range). These axes were then con- verted into densities of occurrences using the “ecospat” package (Di Cola et al. 2017). To understand which factors were the most impactful on the difference in native and alien distributions of the mantises, I plotted the contributions of the variables for the first two PCs with the fviz_pca_var() and the fviz_contrib() functions of the “factoextra” package (version 1.0.7: Kassambara and Mundt 2020). Particularly, the last function also detects which variables are significantly important in the considered PCs with a reference line; every variable above said line should be regarded as impor- tant (Kassambara and Mundt 2020). Overlap of occurrence densities in the environmental space was evaluated by Schoener’s D (D; Schoener 1968), which ranges from 0 (no overlap) to 1 (totally overlapping). Niche expansion (E) followed the definition of Bates et al. (2020), which regards it as the percentage of abiotic features in the non-native range that are not present in the colonized one. I considered E as significant if its value was > 10%. Unfilling (U) followed the definition of Petitpierre et al. (2012), and I considered it significant if it was > 10%, as done with E in previous studies (Petitpierre et al. 2012, Lo Parrino et al. 2023). To evaluate if similarities/differences were caused by random effects, I ran niche equivalency tests following Bates et al. (2020). I calculated a potential correlation between the number of native or alien points and D overlap for each intra-specific analysis using Kendall's rank correlation (t; Kendall 1938), following Bates et al. (2020), using the “stats” package (version 4.1.1; R Core Team 2022). Results If analyzed as H. transcaucasica sensu stricto, the alien Hiero- dula taxon seemed to have moderately high overlap with its na- tive range, low levels of expansion, and no unfilling (D = 0.5, E = 3.06%, and U = 0%; Table 1, Fig. 2A, Suppl. material 1: fig. S12A). When considered as H. tenuidentata sensu stricto, there was low overlap between native and alien ranges, high expansion, and very high unfilling (D = 0.119, E = 21.55%, and U = 84.37%; Table 1, Suppl. material 1: figs S13, $14). When considered as H. tenui- dentata sensu Battiston et al. (2018), the native-alien niche overlap was lower than the one considered for H. transcaucasica sensu stricto, and there was no expansion, while unfilling was high (D = 0.476, E=0%, and U = 47.77%; Table 1, Suppl. material 1: figs $15, S16). For H. patellifera, overlap was low and both expansion and unfilling were high when using both the dataset with no dubious (i.e., the previously mentioned uncertain data from Nepal, India, and Phil- ippines) points (D = 0.18, E = 24.05%, and U = 20.97%; Table 1, Fig. 2B, Suppl. material 1: fig. $12B) and the one with all the points of specimens regarded as the same species (D = 0.179, E = 18.39%. and U = 23.72%; Table 1, Suppl. material 1: figs $17, $18). T. angustipennis showed almost no overlap between native and alien ranges and extremely significant expansion and unfill- ing with both the full dataset (D = 0.01, E = 96.56%, and U = 97.35%; Table 1, Fig. 3A, Suppl. material 1: fig. $19A) and with the one using the Anderson (2019) distribution map (D = 0.009, E = 98.51%, and U = 98.85%; Table 1, Suppl. material 1: figs S20, $21). T. sinensis also had very low overlap with its native range and sig- JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) M. DE VIVO 221 Table 1. Number of points and metric values per each evaluated taxon. D is rounded up to the third decimal number. The p and 1 values reported in the native and alien points’ columns come from the Kendall test, while the p values in the last columns come from the niche equivalency tests. Taxon Alien points (p = 1,T=0) H. transcaucasica sensu stricto 347 H. tenuidentata sensu stricto 347 H. tenuidentata sensu Battiston et al. (2018) 347 H. patellifera (no dubious points) 89 H. patellifera (all) 89 T. angustipennis (full dataset) 292 T. angustipennis (trimmed) 262 T. sinensis (full dataset) 3315 T. sinensis (trimmed) 3065 A H. transcaucasica density of occurrence H. patellifera (no dubious points) 0.4 0.6 0.8 1.0 density of occurrence 0.2 0.0 -6 -4 -2 0 2 4 Axis 1 Fig. 2. Intra-species niche dynamics analyses for the Hierodula spe- cies. A. Results for H. transcaucasica sensu stricto; B. Results for H. patellifera without dubious points. Blue: niche stability. Red: niche expansion. Green: unfilling. U and E values are in Table 1. The analyses for H. tenuidentata sensu stricto and H. tenuidentata sensu Battiston et al. (2018) are reported in Suppl. material 1: figs $14, S16. The analyses for the whole H. patellifera dataset are available in Suppl. material 1: fig. S18. Native points (p = 0.3991, D E U p t= 0.229) 173 O75 3.06% 0% <0.0001 180 0.119 21.55% 84.37% <0.0001 353 0.476 0% 47.77% <0.0001 835 0.18 24.05% 20.97% <0.0001 852 0.179 18.39% 23.72% <0.0001 161 0.01 96.56% 97.35% <0.0001 161 0.009 98.51% 98.85% <0.0001 462 0.022 36.62% 66.59% <0.0001 462 0.02 54.04% 95.97% <0.0001 Table 2. Chosen bioclimatic variables and their VIF values. The descriptions were taken from the WorldClim website. Variable Description VIF BIO2 Mean Diurnal Range (Mean of monthly 2.439 (max temp - min temp)) BIO3 Isothermality (BIO2/BIO7) (x100) 2.54 BIO5 Max Temperature of Warmest Month Ql Be BIO8 Mean Temperature of Wettest Quarter 2.606 BIO9 Mean Temperature of Driest Quarter 3.886 BIO13 Precipitation of Wettest Month 3.567 BIO14 Precipitation of Driest Month 2.609 BIO18 Precipitation of Warmest Quarter 3.415 BIO19 Precipitation of Coldest Quarter 3.014 nificant expansion and unfilling in the alien one, both when the whole possible range was evaluated (D = 0.022, E = 36.62%, and U = 66.59%; Table 1, Fig. 3B, Suppl. material 1: fig. S19B) and when only the trimmed distribution was considered (D = 0.02, E = 54.04%, and U = 95.97%; Table 1, Suppl. mate- rial 1: figs S22, $23). All the intra-species analyses were significant according to niche equivalency tests (p < 0.0001; Table 1). Therefore, their D values were not generated by random chance. Furthermore, none of the t values were significantly different from random, both for native points-D correlation (t = 0.229, p = 0.3991) and for alien points-D correlation (t = 0, p = 1; Table 1). These results highlight the lack of effect the number of occurrence points had on the re- sults. In addition, the Mantel tests proved that there was no spatial autocorrelation at 10 km (Suppl. material 1: figs S24-S27). In terms of the importance of the variables in the intra-spe- cies analyses, BIO19 (Precipitation of Coldest Quarter) was the most important variable if we considered one of the Hierodula species as H. transcaucasica (Table 3, Suppl. material 1: figs $28, S29), followed by BIO18 (precipitation of warmest quarter), BIO 13 (precipitation of wettest month), bio14 (precipitation of driest month), BIO9 (mean temperature of driest quarter), BIO3 (isothermality), and BIO5 (Table 3, Suppl. material 1: fig. $28). When the taxon was regarded as H. tenuidentata sensu stricto, BIO9 was the most important variable (Table 3, Suppl. material 1: figs S30, S31), followed by BIO8 (mean temperature of wet- test quarter), BIO5, BIO14, and BIO3 (Table 3, Suppl. material 1: fig. S30). BIO5 was the most important for the difference be- tween native and alien niches when one of the Hierodula species was considered as H. tenuidentata sensu Battiston et al. (2018) (Table 3, Suppl. material 1: figs S32, S33), followed by BIO18, JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 222 Table 3. Important variables in the intra-species analyses. M. DE VIVO Taxon Most important variable Other important variables H. transcaucasica sensu stricto BIOI9 BIO18, BIO13, BIO14, BIOS, BIO3, BIO5 H. tenuidentata sensu stricto BIOS BIO8, BIO5, BIO14, BIO3 H. tenuidentata sensu Battiston et al. (2018) BIOS BIO 18, BIO3, BIO13, BIO14 H. patellifera (no dubious points) BIO9 BIOS, BIO14, BIO8, BIO19 H. patellifera (all) BIO9 BIOS, BIO14, BIO8, BIO19 T. angustipennis (full dataset) BIO18 BIO13, BIO19, BIO2, BIO9 T. angustipennis (trimmed) BIO2 BIO19, BIO14, BIO13, BIO18, BIO9, BIO3 T. sinensis (full dataset) BIO13 BIO18, BIO19, BIO9, BIO2 T. sinensis (trimmed) BIO13 BIO18, BIO19, BIOS, BIO14, BIO2, BIO3 A T. angustipennis (all) 0.8 1.0 0.6 density of occurrence 0.4 0.2 = \ 0.0 Axis 1 T. sinensis (all) 1.0 0.8 0.6 density of occurrence 0.4 0.2 0.0 Axis 1 Fig. 3. Intra-species niche dynamics analyses for the Tenodera spe- cies with the full datasets. A. Results for T: angustipennis; B. Results for T. sinensis. Blue: niche stability. Red: niche expansion. Green: unfilling. U and E values are in Table 1. The analyses with the trimmed points are available in Suppl. material 1: figs $21, $23. BIO3, BIO13, and BIO14 (Table 3, Suppl. material 1: fig. $32). For H. patellifera, the dataset trimming did not lead to a change in inferences given that BIO9 was the most important variable; the other significant ones were, in order of importance, BIO5, BIO14, BIO8, and BIO19 for both datasets (Table 3, Suppl. material 1: figs S34-S37). For T. angustipennis, if the full dataset was considered, BIO18 was the most important variable (Table 3, Suppl. material 1: figs S38, S39), while the other important ones were BIO13, BIO19, BIO2 (mean diurnal range), and BIO9 (Table 3, Suppl. material 1: fig. $38). If only the Anderson (2019) distribution was con- sidered, BIO2 was the most important variable (Table 3, Suppl. material 1: figs S40, S41), followed by BIO19, BIO14, BIO13, BIO18, BIO9, and BIO3 (Table 3, Suppl. material 1: fig. S40). Concerning T: sinensis, BIO13 was the most important variable in both the full and trimmed datasets (Table 3, Suppl. material 1: figs S42-S45). In the analyses with the full dataset, BIO18 followed suit, with a seemingly negligible difference between it and the most important variable, while the other significant ones were BIO19, BIO9, and BIO2 (Table 3, Suppl. material 1: fig. S42). In the trimmed analysis, the other important variables were BIO18, BIO19, BIO9, BIO14, and BIO3 (Table 3, Suppl. material 1: fig. S44). In the inter-taxa analyses, the two Hierodula taxa in Europe had high overlap (D = 0.57), and the analysis was not significantly influenced by random chance (p = 0.008; Fig. 4A). The Tenodera species also had high overlap according to both whole (D = 0.665, p < 0.0001; Fig. 4B) and trimmed datasets (D = 0.524, p < 0.0001; Suppl. material 1: fig. S44). BIO13 was the most important vari- able for the comparison between the two Hierodula followed by BIO18, BIOS, BIO19, and BIO 14 (Table 4, Suppl. material 1: figs S47, S48). In the comparison between the two Tenodera species, BIO 13 was the most important variable for both analyses (Table 4, Suppl. material 1: figs S49-S52). When the full dataset was used, BIOS, BIO2, BIO9, and BIO3 were the other important variables, while BIO8, BIO3, BIO14, BIO2, BIO9, and BIO 13 were the other important variables in the trimmed dataset (Table 4, Suppl. mate- rial 1: figs S49, S51). Table 4. Important variables in the inter-species analyses. Comparison Mostimportant Other important variables variable Hierodula in Europe BIO13 BIO 18, BIOY, BIO19, BIO14 Tenodera in North BIO19 BIO5, BIO2, BIOY, BIO3 America (full dataset) Tenodera in North BIO19 BIO8, BIO3, BIO14, BIO2, America (trimmed) BIOS, BIO13 JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) M. DE VIVO Hierodula in Europe T : fo Te 0.008 ap) o a ra ae 5 O° QO S co) ro) anil ] mint | [ T T T -4 -2 0 2 4 6 Between-Class PCA 1 Tenodera in North America (whole dataset) co ro) 0.665 \<0.000 Ww oS Tt s o ra 7 ao N ro) S io) oO Between-Class PCA 1 Fig. 4. Inter-species analyses for the niche overlap in displaced congeneric species. A. Results for Hierodula species in shared Euro- pean areas. Yellow: the Hierodula species with taxonomic contro- versial status (H. tenuidentata/H. transcaucasica). Green: H. patel- lifera. Grey: shared niche. B. Results for Tenodera species in North America with the whole dataset. Green: T: angustipennis. Yellow: T. sinensis. Grey: shared niche. The analyses with the trimmed points for the Tenodera species are available in Suppl. material 1: fig. S46. Discussion Intra-species analyses and the impact of taxonomic impediments.—It seems that there is both little overlap between the native and alien ranges in the considered mantises and high levels of unfilling and expansion if one of the Hierodula taxa is regarded as H. tenuidentata sensu stricto. Potentially, a low D overlap might be caused by recent introduction (Liu et al. 2020), which would explain the alien-na- tive overlap and high unfilling for H. tenuidentata sensu stricto and H. patellifera (Battiston et al. 2018, 2020). However, given the time period in which T. angustipennis and T. sinensis were introduced in continental North America (Snyder and Hurd 1995, Ander- son 2019), it seems unlikely that time is a factor for the low D in Tenodera. In addition, the Tenodera species showed high levels of unfilling and different degrees of expansion, with T. angustipennis showing an E larger than that of T. sinensis, although the metric 223 was significant for both taxa. This may seem counterintuitive given that T. sinensis is generally thought to have a higher impact in the United States than T. angustipennis (Snyder and Hurd 1995, Hurd et al. 2004, Anderson 2019). However, it must be noted that potential impacts may not be related to expansion (Petitpierre et al. 2012, Bates et al. 2020), and species with smaller native ranges may ex- perience greater niche shifts (Bates et al. 2020). At the same time, it seems that T: sinensis is constrained by abiotic and movement vari- ables in certain areas given the doubts about its ability to establish populations outside of the northeast of the United States (Ander- son 2019) and the reports of very high rates of nymph mortality in dry conditions (Hurd et al. 2004). The last point also explains the importance of precipitation variables such as BIO13 for this species in the intra-species analyses. Thus, a low D overlap and high unfill- ing may imply that Tenodera species are constrained by abiotic and movement factors in the United States, and the high expansion may be explained by a relatively small native range (Bates et al. 2020) or by the chosen variables (see Lo Parrino et al. 2023). However, it is also recognized that T. sinensis shifted its maturation period over the years in the United States, showing a degree of biological plas- ticity potentially triggered by temperature changes, which may also lead to a range expansion in this species (Hurd et al. 2020). In terms of T. angustipennis, the two datasets showed different variables as the most important (BIO18, a precipitation-related variable, for the full dataset and BIO2, a temperature-related vari- able, for the trimmed dataset). This might mean that this species is constrained by more factors than the congeneric species. However, two precipitation variables always ranked among the three most important for both Tenodera, highlighting the importance of rain to this genus in North America. Given the potential increase in rain- fall in some eastern areas of the United States in the future (Harp and Horton 2023), it is possible that both species may increase their non-native range. It must be noted that T: sinensis also feeds on T. angustipennis (Hurd 1988, Snyder and Hurd 1995), and it seems to be more prolific in nymph production (Anderson 2019). These biological factors, together with movement ones, may re- duce the impact and the dispersal of T. angustipennis in the United States and show how the competition between two alien species shapes the distribution of sister species in non-native ranges. It is possible that the two Hierodula taxa occupy niches that are not occupied by other European mantises, given their arboreal ecol- ogy (Battiston et al. 2020). At the same time, these niches might be only a small part of their native ones, and therefore they may survive in Europe only by expanding them (if one of the taxa is considered as H. tenuidentata sensu stricto). However, such results may be influenced by taxonomic status. For example, considering the most taxonomi- cally problematic Hierodula taxon as H. transcaucasica sensu stricto, H. tenuidentata sensu stricto, or H. tenuidentata sensu Battiston et al. (2018) clearly changed the value of the examined metrics and the inference about this taxon. Specifically, when one of the Hierodula taxon was considered as H. tenuidentata sensu Battiston et al. (2018) or H. transcaucasica sensu stricto, its alien niche tended to be relatively similar to its native one, hinting at higher potential invasiveness if mantises respect the general pattern in terrestrial ectotherms (Bates et al. 2020, Liu et al. 2020). However, the higher D and the smaller E in H. transcaucasica sensu stricto is surprising given the smaller range compared to H. tenuidentata sensu Battiston et al. (2018). This, to- gether with the lack of U, might imply a potential natural expansion of said species/population that started in the past, as suggested by Pintilioaie et al. (2023), which was also considered as a hypothesis for the Crimean and Balkan populations (Pushkar and Kavurka 2016, Battiston et al. 2018, Shcherbakov and Battiston 2020). Critically, the JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 224 taxon (regardless of the controversial taxonomy) has been shown to be more widespread in Europe than H. patellifera (Sevgili and Yilmaz 2022), potentially due to its better ability to circumvent movement barriers (Battiston et al. 2020, Moulin 2020, William Di Pietro pers. comm.). Therefore, the higher D, lower E, and lack of U may high- light a higher BAM overlap than that of the sister species. However, it is unclear if the population that is traditionally regarded as H. tran- scaucasica is the source one for the European alien range given that such a hypothesis can only be tested by molecular means (Bates and Bertelsmeier 2021) and the fact that the lower D in H. patellifera or H. tenuidentata sensu stricto may also be explained by recent introduction (Battiston et al. 2018, 2020, Liu et al. 2020). In addition, it is unclear whether the observations of this Hierodula complex belong to one or more lineages (i.e., populations from both South Asia and Caucasus). In terms of the abiotic factors, the importance of the variables in the controversial Hierodula taxon seemed to vary according to the tax- onomic assignment. Generally, the Hierodula species had temperature variables as the most important, which makes sense particularly for H. patellifera (Shin et al. 2023). However, if we consider the Hierodula taxon to be H. transcaucasica, BIO19 was the most relevant variable for explaining the difference between native and non-native ranges. Previous modeling attempts have shown that both temperature and precipitation are important for H. tenuidentata sensu Battiston et al. (2018), and it is thought that the taxon could spread further in Eu- rope due to climate change (Pintilioaie et al. 2023). Particularly, it is recognized that precipitation in cold months has been increasing due to increased temperature in several European areas (Hyn¢ica and Huth 2019), and therefore both factors (increased temperature and increased precipitation in coldest months) might have been impor- tant in the dispersal in the mantis, regardless of the source population. Inter-species analyses.—There is relatively high overlap between sis- ter taxa occupying the same alien range. This may not be surprising given that such species potentially share areas in their native ranges as well (Snyder and Hurd 1995, Anderson 2019, Liu et al. 2021). As stated above, the natural spreading ability of H. patellifera in Europe, where it seems to be most likely to be in places where it is transported through anthropogenic means (Battiston et al. 2020, Moulin 2020, William Di Pietro pers. comm.), appears to be less than the spreading of the other Hierodula taxon, which seems to be able to become more widespread in this continent (Battiston et al. 2018, Shcherbakov and Battiston 2020, van der Heyden and Schwarz 2021, Sevgili and Yilmaz 2022, Moulin and Rouard 2023, Pintilio- aie et al. 2023), although H. patellifera was probably able to reach Croatia from Italy by natural means (Martinovicé et al. 2022). Pre- cipitation variables mostly explained the difference between the two non-native species in Europe, particularly BIO13. Therefore, precipi- tation might be a factor explaining the difference between the two. The ecology of these two species/taxa is mostly known from labo- ratory observations in areas in which they do not co-exist (Leong 2009, Kharabadze et al. 2022, Mirzaee et al. 2023), and therefore it is challenging to make assumptions about potential competition and differences in abiotic requirements between the taxa, particular- ly in non-native environments where different populations of two species from non-shared environments may meet for the first time. However, given that it is also observed in Tenodera species (Anderson 2019), it would not be surprising if BAM factors shape the distribu- tion of Hierodula species in Europe. Critically, given what is observed by the distribution of H. patellifera, it seems that the other Hierodula taxa might have a competitive movement advantage that may result in a reduction in the impact of the former species, as observed in other taxa (e.g., Van Riel et al. 2009, Russell et al. 2014). M. DE VIVO Concerning the Tenodera species, the high overlap is consist- ent with the observations regarding their distribution (Anderson 2019) and competition (Hurd 1988, Snyder and Hurd 1995). Therefore, a high overlap in the area they share is expected, al- though T. sinensis is more widespread for biological reasons (An- derson 2019), highlighting the effect of competition on the dis- tribution of alien sister species. In both the datasets used, BIO19 was the variable that explained more of the difference between the two, which may be a limiting factor for T. angustipennis given its reduced geographical range in both native and non-native envi- ronments compared to T: sinensis coupled with the biological and movement constraints discussed in the previous section. Limitations of the study.—As I stated in the previous sections, taxonom- ic statuses and misidentification could radically change the results of this study. In the future, an integrative taxonomic approach in which both morphological and molecular characters are used (e.g., Shim et al. 2021) is needed to understand which species and which popula- tion is present as an alien taxon, which is very important for tracking potential changes between the original native niche and the non-na- tive one (Liu et al. 2020, Bates and Bertelsmeier 2021). Another issue is the lack of information on mantis ecology. Specifically, although several mantises have been called “invasive” in some studies (e.g., Snyder and Evans 2006, Crowder and Snyder 2010, Fea et al. 2013, Battiston et al. 2018, Schwarz and Ehrmann 2018, Moulin 2020, Shcherbakov and Govorov 2020, Moulin and Rouard 2023, Pintilioaie et al. 2023), no Mantodea species is cur- rently present in the Global Invasive Species Database (GSID) made by the Invasive Species Specialist Group (ISSG) of IUCN (GISD 2023). IUCN defines species as invasive if they are “nega- tively impacting native biodiversity, ecosystem services or human economy and well-being” (IUCN 2023). However, the impact of alien mantises on ecosystems is not very well understood (Fea et al. 2013, Battiston et al. 2020, 2022, Connors et al. 2022), and few studies have evaluated the impact of displaced Mantodea species on native taxa in the same order or niche (see Snyder and Evans 2006 and Fea et al. 2013). This causes serious problems with the possible listing of mantises as alien or invasive alien; specifically, a lack of baseline may cause concerns in analyzing potential patterns shared by non-native mantises and differentiating potential inva- sive ones (Liu et al. 2020). That being said, several alien mantises (e.g., all the species evaluated in this study) would be regarded as invasive if looser definitions are used (e.g., Gippet and Bertelsmei- er (2021) regard all the species “with at least one established popu- lation outside of the native range” as invasive, even if no potential impact is known), but finding a unified definition is still challeng- ing, and only recently has a framework been proposed for find- ing a shared terminology among ecologists (Almena et al. 2023). Critically, by using shared terminology, we should consider all the mantises analyzed here as “invasive non-native” species, which are “established non-native species” spreading and establishing popu- lations further from their introduction points (Almena et al. 2023). In addition, there are biases toward some countries for both research efforts (Battiston et al. 2022) and presence records. For example, there were almost nine times more observation points for T. sinensis in the United States than in Asia (Table 1). Critically, some areas of the world need to be studied due to potential inva- sions of mantises or to better understand the ecology of invad- ing mantises. For example, I found several points on iNaturalist of a mantis identified as H. patellifera sensu stricto in the Seychelles (identification by Chih-Ting Hsu), but the only alternative source I could find about its presence in the nation was an online news JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) M. DE VIVO source (Seychelles Nation 2009). Therefore, it is not known wheth- er such species had any impact on the Seychellois ecosystems. An- other potential example is Sphodromantis viridis (Forskal 1775), which is becoming widespread in Europe (Battiston et al. 2017, Martinovic et al. 2022). Despite its increasing presence in this area, we still know very little about its ecology (Battiston et al. 2017), and in both iNaturalist and GBIF, there are a very limited number of records from its native range (GBIF 2023a, iNaturalist 2023). Conclusions Overall, the analyzed mantises did not exhibit high overlap between their native and alien niches, while showing high expan- sion and unfilling, which shows that mantises can establish popu- lations in environments dissimilar to their native ones while hav- ing limits in their spread at the same time. However, the situation for H. tenuidentata/transcaucasica must be clarified given the po- tential change of results and inferences caused by the taxonomic issues in such species complex. Given these considerations, I argue that future taxonomic and ecological studies are critically needed to better understand mantis ecology and taxonomy, especially in less studied areas of the world and in natural environments where alien species are present, which would also help to understand whether the invasion was started by different or specific popula- tions of the same species with different requirements. However, given my results, coupled with previous modeling attempts and inferences by other research teams, it is highly likely that the un- clear Hierodula taxon could be the population/species referred to as H. transcaucasica, given its niche overlap that seems to resemble the general pattern of invasive ectothermic terrestrial species, al- though only genetic studies could prove these assumptions. Data availability statement The scripts and the elaborated presence points with the values per each variable are present in Zenodo (https://zenodo.org/re- cords/10101384, DOI: 10.5281/zenodo. 10101384). Acknowledgements I thank Chih-Ting Hsu (National Chung Hsing University, Tai- chung, Taiwan) for discussion and suggestions on the taxonomy of the evaluated species, William Di Pietro (World Biodiversity Asso- ciation Onlus, Verona, Italy) for discussion regarding the presence of H. patellifera in Italy, Howon Rhee (University of Trier, Trier, Germany) and Taewoo Kim (National Institute of Biological Resources, Incheon, South Korea) for discussion regarding H. patellifera in Korea, Trevor Padgett (Taiwan International Graduate Program on Biodiversity and Tunghai University, Taichung, Taiwan) for revising both the grammar and the content of the manuscript, Jen-Pan Huang (Academia Sinica, Taipei, Taiwan) for suggestions on the manuscript, and the Orthop- terists’ Society for the fee waiver. For this study, I was supported by the Taiwan International Graduate Program (TIGP) through the TIGP Research Performance Fellowship 2022 and by internal research sup- port from Academia Sinica and the University of Trier. References Aiello-Lammens ME, Boria RA, Radosavljevic A, Vilela B, Anderson RP (2015) spThin: An R package for spatial thinning of species occur- rence records for use in ecological niche models. Ecography 38: 541-545. https://doi.org/10.1111/ecog.01132 225 Almena IS, Balzani P, Carneiro L, Cuthbert RN, Macedo R, Tarkan AS, Ahmed DA, Bang A, Bacela-Spychalska K, Bailey SA, Baudry T, Balles- teros L, Bortolus A, Briski E, Britton JR, Buric M, Camacho-Cervantes M, Cano-Barbacil C, Copilas-Ciocianu D, Coughlan N, Courtois P, Csabai Z, Dalu T, Santis VD, Dickey JWE, Dimarco R, Falk-Andersson J, Fernandez R, Florencio M, Franco ACS, Garcia-Berthou E, Gian- netto D, Glavendekic M, Grabowski M, Heringer G, Herrera I, Wei H, Kamelamela KL, Kirichenko NI, Kouba A, Kourantidou M, Kurtul I, Laufer G, Liptak B, Liu C, Lopez-Lopez E, Lozano V, Mammola S, Mar- chini A, Meshkova V, Meyerson L, Milardi M, Musolin DL, Nunez M, Oficialdegui FJ, Patoka J, Pattision Z, Petrusek A, Pincheira-Donoso D, Piria M, Probert A, Rasmussen JJ, Renault D, Ribeiro E Rilov G, Robinson TB, Sanchez A, Schwindt E, South J, Stoett P, Verreycken H, Vilizzi L, Wang Y-J, Watari Y, Wehi PM, Weiperth A, Wiberg-Larsen P, Yapici S, Yogurtcuoglu B, Zenni R, Galil BS, Dick JTA, Russell J, Ric- ciardi A, Simberloff D, Bradshaw CJA, Haubrock PJ (2023) Taming the terminological tempest in invasion science. https://ecoevorxiv. org/repository/view/5904/ [Accessed October 23, 2023] Anderson K (2019) Praying mantises of the United States and Canada. 2"¢ edn. Independently published, 297 pp. Bates OK, Ollier S, Bertelsmeier C (2020) Smaller climatic niche shifts in invasive than non-invasive alien ant species. Nature Communications 11: 5213. https://doi.org/10.1038/s41467-020-19031-1 Bates OK, Bertelsmeier C (2021) Climatic niche shifts in introduced spe- cies. Current Biology 31: R1252-R1266. https://doi.org/10.1016/j. cub.2021.08.035 Battiston R, Andria S, Ruzzante G (2017) The silent spreading of a giant mantis: a critical update on the distribution of Sphodromantis viridis (Forskal, 1775) in the Mediterranean islands (Mantodea: Mantidae). Onychium 13: 25-30. https://doi.org/10.5281/zenodo.546318 Battiston R, Leandri F, Di Pietro W, Andria S (2018) The giant Asian mantis, Hierodula tenuidentata, spreads in Italy: a new invasive alien species for the European fauna? Biodiversity Journal 9: 399-404. https://doi. org/10.31396/Biodiv.Jour.2018.9.4.399.404 Battiston R, Amerini R, Di Pietro W, Guariento LA, Bolognin L, Moretto, E (2020) A new alien mantis in Italy: is the Indochina mantis Hierodula patellifera chasing the train for Europe? Biodiversity Data Journal 8: e50779. https://doi.org/10.3897/BDJ.8.e50779 Battiston R, Di Pietro W, Anderson K (2022) The pet mantis market: a first overview on the praying mantis international trade (Insecta, Mantodea). Journal of Orthoptera Research 31: 63-68. https://doi. org/10.3897/jor.31.71458 Bivand R, Pebesma EJ, G6mez-Rubio V (2013) Applied spatial data analy- sis with R. Second edition. Springer, New York, 405 pp. https://doi. org/10.1007/978-1-4614-7618-4 Broennimann O, Fitzpatrick MC, Pearman PB, Petitpierre B, Pellissier L, Yoccoz NG, Thuiller W, Fortin M-J, Randin C, Zimmermann NE, Graham CH, Guisan A (2012) Measuring ecological niche over- lap from occurrence and spatial environmental data. Global Ecol- ogy and Biogeography 21: 481-497. https://doi.org/10.1111/j.1466- 8238.2011.00698.x Brunner de Wattenwyl K (1878) Naturwissenschaftliche Betrage zur Kennt- nis der Kaukasuslander auf Grund seiner Sammelbeute. Naturwissen- schatliche Beitrage zum Kaukasus 5: 87-90. Burmeister HC (1838) Handbuch der Entomologie. Fangschrecken, Man- todea. Handbuch der Entomologie 2 (V-VIII): 517-552. Connors MG (2023) Studies in Australian Mantodea: New synonyms, notes on distributions, and an updated checklist. Zootaxa 5239: 477-499. https://doi.org/10.11646/zootaxa.5239.4.2 Connors MG, Chen H, Li H, Edmonds A, Smith KA, Gell C, Clitheroe K, Miller IM, Walker KL, Nunn JS, Nguyen L, Quinane LN, Andreoli CM, Galea JA, Quan B, Sandiford K, Wallis B, Anderson ML, Canziani EV, Craven J, Hakim RRC, Lowther R, Maneylaws C, Menz BA, New- man J, Perkins HD, Smith AR, Webber VH, Wishart D (2022) Citi- zen scientists track a charismatic carnivore: Mapping the spread and impact of the South African Mantis (Miomantidae, Miomantis caffra) in Australia. Journal of Orthoptera Research 31: 69-82. https://doi. org/10.3897/jor.31.79332 JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) 226 Crowder DW, Snyder WE (2010) Eating their way to the top? Mechanisms underlying the success of invasive insect generalist predators. Biological Invasions 12: 2857-2876. https://doi.org/10.1007/s10530-010-9733-8 De Vivo M, Huang J-P (2022) Modeling the geographical distributions of Chordodes formosanus and its mantis hosts in Taiwan, with consid- erations for their niche overlaps. Ecology and Evolution 12: e9546. https://doi.org/10.1002/ece3.9546 Di Cola V, Broennimann O, Petitpierre B, Breiner FI, D‘Amen M, Randin C, Engler R, Pottier J, Pio D, Dubuis A, Pellissier L, Mateo RG, Hordijk W, Salamin N, Guisan A (2017) ecospat: an R package to support spa- tial analyses and modeling of species niches and distributions. Ecog- raphy 40: 774-787. https://doi.org/10.1111/ecog.02671 Dray S, Dufour A-B (2007) The ade4 Package: Implementing the Duality Diagram for Ecologists. Journal of Statistical Software 22. https://doi. org/10.18637/jss.v022.i04 Fea MP, Stanley MC, Holwell GI (2013) Fatal attraction: Sexually canni- balistic invaders attract native mantids. Biology Letters 9: 20130746. https://doi.org/10.1098/rsbl.2013.0746 Fick SE, Hijmans RJ (2017) WorldClim 2: new 1-km spatial resolution cli- mate surfaces for global land areas. International Journal of Climatol- ogy 37: 4302-4315. https://doi.org/10.1002/joc.5086 Forskal P (1775) Descriptiones Animalium Avium, Amphibiorum, Pisci- um, Insectorum, Vermium; quae in Itinere Orientall observati Petrus Forskal. Prof. Haun. Post morten Acutoris editt Carsten Nieburhr, 164 pp. https://doi.org/10.5962/bh1.title.2154 Gardiner MM, Allee LL, Brown PM, Losey JE, Roy HE, Smyth RR (2012) Lessons from lady beetles: accuracy of monitoring data from US and UK citizen-science programs. Frontiers in Ecology and the Environ- ment 10: 471-476. https://doi.org/10.1890/110185 GBIForg (2023a) Occurrence Download. https://doi.org/10.15468/ dl.y59wdq [Accessed 30 October 2023] GBIForg (2023b) Occurrence Download. https://doi.org/10.15468/ dl.cx7em7 [Accessed 18 September 2023] Gippet JMW, Bertelsmeier C (2021) Invasiveness is linked to greater com- mercial success in the global pet trade. Proceedings of the National Academy of Sciences 118(14): e2016337118. https://doi.org/10.1073/ pnas.2016337118 GISD (2023) Global Invasive Species Database. http://www.iucngisd.org/ gisd/ [Accessed March 30, 2023] Harp RD, Horton DE (2023) Observed Changes in Interannual Precipita- tion Variability in the United States. Geophysical Research Letters 50: €2023GL104533. https://doi.org/10.1029/2023GL104533 Heberling JM, Miller JI, Noesgaard D, Weingart SB, Schigel D (2021) Data integration enables global biodiversity synthesis. Proceedings of the National Academy of Sciences 118: e2018093118. https://doi. org/10.1073/pnas.2018093118 Hijmans RJ (2022) raster: Geographic Data Analysis and Modeling. R package version 3.5-15. https://CRAN.R-project.org/package=raster Hurd LE (1988) Consequences of divergent egg phenology to predation and coexistence in two sympatric, congeneric mantids (Orthoptera: Man- tidae). Oecologia 76: 549-552. https://doi.org/10.1007/BF00397868 Hurd LE, Mallis RE, Bulka KC, Jones AM (2004) Life History, Environment, and Deme Extinction in the Chinese Mantid Tenodera aridifolia sinen- sis (Mantodea: Mantidae). Environmental Entomology 33: 182-187. https://doi.org/10.1603/0046-225X-33.2.182 Hurd LE, Cheng KX, Abcug J, Calhoun LV, Geno ME, Merhige RR, Rosenthal IH (2020) Climate Change, Succession, and Reproductive Success of a Praying Mantid. Annals of the Entomological Society of America 113: 202-206. https://doi.org/10.1093/aesa/saz070 Hyn¢éica M, Huth R (2019) Long-term changes in precipitation phase in Europe in cold half year. Atmospheric Research 227: 79-88. https://doi.org/10.1016/j.atmosres.2019.04.032 iNaturalist (2023) iNaturalist. https://www.inaturalist.org/ [Accessed 18 September 2023] iNaturalist Blog (2023) Thank you for helping generate most GBIF records for most species since 2020. https://www.inaturalist.org/blog/76606- thank-you-for-helping-generate-most-gbif-records-for-most-species- since-2020 [Accessed March 24, 2023] M. DE VIVO IPBES (2023) IPBES Invasive Alien Species Assessment: Summary for Policymakers. IPBES secretariat, Bonn, Germany, 56 pp. https://doi. org/10.5281/zenodo. 10004329 IUCN (2023) Invasive Alien Species. https://www.iucn.org/our-work/top- ic/invasive-alien-species [Accessed March 30, 2023] Jensen D, Svenson GJ, Song H, Whiting MF (2009) Phylogeny and evo- lution of male genitalia within the praying mantis genus Tenodera (Mantodea:Mantidae). Invertebrate Systematics 23: 409. https://doi. org/10.1071/1S09004 Kassambara A, Mundt F (2020) factoextra: Extract and Visualize the Results of Multivariate Data Analyses. R package version 1.0.7. https://CRAN.R-project.org/package=factoextra. Kendall MG (1938) A new measure of rank correlation. Biometrika 30: 81-93. https://doi.org/10.1093/biomet/30.1-2.81 Kharabadze N, Chkhaidze N, Abramishvili T, Goktiirk T, Lobjanidze M, Burjanadze M (2022) First report of Hierodula transcaucasica (Brun- ner von Wattenwyl, 1878) predation on the Halyomorpha halys (Stal, 1855) in Georgia. International Journal of Tropical Insect Science 42: 3283-3292. https://doi.org/10.1007/s42690-022-00826-2 Kolnegari M (2023) Mantodea of Iran: A review-based study. Journal of Or- thoptera Research 32: 177-188. https://doi.org/10.3897/jor.32.97388 Leong T (2009) Oviposition and hatching in the praying mantis, Hierodula patellifera (Serville) in Singapore (Mantodea: Mantidae: Paramanti- nae). Nature in Singapore 2: 55-61. Liu C, Wolter C, Xian W, Jeschke JM (2020) Most invasive species largely conserve their climatic niche. Proceedings of the National Academy of Sciences 117: 23643-23651. https://doi.org/10.1073/ pnas.2004289117 Liu Q-P, Liu Z-J, Wang G-L, Yin Z-X (2021) Taxonomic revision of the pray- ing mantis subfamily Hierodulinae of China (Mantodea: Mantidae). Zootaxa 4951: 401-433. https://doi.org/10.11646/zootaxa.4951.3.1 Lo Parrino E, Falaschi M, Manenti R, Ficetola GF (2023) All that changes is not shift: methodological choices influence niche shift detection in freshwater invasive species. Ecography 2023: e06432. https://doi. org/10.1111/ecog.06432 Lockwood JL, Welbourne DJ, Romagosa CM, Cassey P, Mandrak NE, Strecker A, Leung B, Stringham OC, Udell B, Episcopio-Sturgeon DJ, Tlusty MF, Sinclair J, Springborn MR, Pienaar EF Rhyne AL, Keller R (2019) When pets become pests: the role of the exotic pet trade in producing invasive vertebrate animals. Frontiers in Ecology and the Environment 17: 323-330. https://doi.org/10.1002/fee.2059 Maceda-Veiga A, Escribano-Alacid J, Martinez-Silvestre A, Verdaguer I, Mac Nally R (2019) What's next? The release of exotic pets continues virtu- ally unabated 7 years after enforcement of new legislation for manag- ing invasive species. Biological Invasions 21: 2933-2947. https://doi. org/10.1007/s10530-019-02023-8 Martinovié M, Cato S, Lengar M, Skejo J (2022) First records of three exotic giant mantid species on the Croatian coast. Journal of Orthoptera Re- search 31: 55-61. https://doi.org/10.3897/jor.31.76075 Meurisse N, Rassati D, Hurley BP, Brockerhoff EG, Haack RA (2019) Com- mon pathways by which non-native forest insects move internation- ally and domestically. Journal of Pest Science 92: 13-27. https://doi. org/10.1007/s10340-018-0990-0 Mirzaee Z, Sadeghi S, Battiston R (2022) Biology and Life Cycle of the Pray- ing Mantid Hierodula tenuidentata Saussure, 1869 (Insecta: Mantodea). Iranian Journal of Science and Technology, Transactions A: Science 46: 1163-1169. https://doi.org/10.1007/s40995-022-01325-2 Mormul RP, Vieira DS, Bailly D, Fidanza K, da Silva VFB, da Graca WJ, Pontara V, Bueno ML, Thomaz SM, Mendes RS (2022) Invasive alien species records are exponentially rising across the Earth. Biological In- vasions 24: 3249-3261. https://doi.org/10.1007/s10530-022-02843-1 Moulin N (2020) When Citizen Science highlights alien invasive species in France, the case of Indochina mantis, Hierodula patellifera (In- secta, Mantodea). Biodiversity Data Journal 8: e46989. https://doi. org/10.3897/BDJ.8.e46989 Moulin N, Rouard J (2023) Hierodula transcaucasica continues its invasion of Western Europe (Mantodea, Mantidae). Bulletin de la Société ento- mologique de France 128: 103-107. https://doi.org/10.32475/bsef_2265 JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2) M. DE VIVO Naimi B, Hamm NAS, Groen TA, Skidmore AK, Toxopeus AG (2014) Where is positional uncertainty a problem for species distribution modelling? Ecography 37: 191-203. https://doi.org/10.1111/j.1600- 0587.2013.00205.x Nguyen DT, Maxwell MR (2008) Stalking of Stationary Prey by a Pray- ing Mantid (Hierodula patellifera Serville) (Mantodea: Mantidae). Entomological News 119: 425-427. https://doi.org/10.3157/0013- 872X-119.4.425 Oshima K (2021) A new subspecies of the mantis Hierodula patellifera (Mantodea: Mantidae) from the Daito Islands, the Ryukyus, Japan. Journal of Orthoptera Research 30: 95-98. https://doi.org/10.3897/ jor.30.62022 Petitpierre B, Kueffer C, Broennimann O, Randin C, Daehler C, Guisan A (2012) Climatic Niche Shifts Are Rare Among Terrestrial Plant Invaders. Science 335: 1344-1348. https://doi.org/10.1126/sci- encen2 15933 Pili AN, Tingley R, van Winkel D, Maria L, Chapple DG (2023) The es- calating global problem of accidental human-mediated transport of alien species: A case study using alien herpetofauna interceptions in New Zealand. Biological Conservation 278: 109860. https://doi. org/10.1016/j.biocon.2022.109860 Pintilioaie A-M, Sfica L, Baltag ES (2023) Climatic Niche of an Invasive Mantid Species in Europe: Predicted New Areas for Species Expan- sion. Sustainability 15: 10295. https://doi.org/10.3390/su151310295 Pushkar T, Kavurka V (2016) New data about the distribution of Hierodula transcaucasica in Ukraine. Ukrainska Entomofaunistyka 7: 77-78. R Core Team (2022) R: A language and environment for statistical comput- ing. R Foundation for Statistical Computing. https://www.R-project. org/ Ratnieks FLW, Schrell F Sheppard RC, Brown E, Bristow OE, Garbuzov M (2016) Data reliability in citizen science: learning curve and the ef- fects of training method, volunteer background and experience on identification accuracy of insects visiting ivy flowers. Methods in Ecology and Evolution 7: 1226-1235. https://doi.org/10.1111/2041- 210X.12581 Russell JC, Sataruddin NS, Heard AD (2014) Over-invasion by function- ally equivalent invasive species. Ecology 95: 2268-2276. https://doi. org/10.1890/13-1672.1 Saussure H de (1869) Essai d‘un Systeme des Mantides. Mittheilungen der Schweizer Entomologischen Gesellschaft 3: 49-73. Saussure H de (1871) Memoires de la Societe de Physique et d'Histoire naturelle de Geneve 21(2): 214 pp. Schoener TW (1968) The Anolis Lizards of Bimini: Resource Parti- tioning in a Complex Fauna. Ecology 49: 704-726. https://doi. org/10.2307/1935534 Schwarz CJ, Ehrmann R (2018) Invasive Mantodea species in Europe. Ar- ticulata 33: 73-90. https://www.biodiversidadvirtual.org/taxofoto/ sites/default/files/schwarz_ehrmann_2018_invasive_mantodea_spe- cies_in_europe.pdf Seebens H, Ess] F Dawson W, Fuentes N, Moser D, Pergl J, PySek P, van Kleunen M, Weber E, Winter M, Blasius B (2015) Global trade will accelerate plant invasions in emerging economies under cli- mate change. Global Change Biology 21: 4128-4140. https://doi. org/10.1111/gcb.13021 Serville JGA (1839) Histoire naturelle des insectes. Orthopteres. Librairie encylopédique de Roret, Paris, [i-xviii,] 776 pp. [pl. 1-14] Sevgili H, Yilmaz K (2022) Contributions of citizen scientists to monitor- ing alien species: the case study on Giant Asian Mantes, Hierodula tenuidentata and H. patellifera (Mantodea: Mantidae). Zoology in the Middle East 68: 350-358. https://doi.org/10.1080/09397140.2022.2 145802 Seychelles Nation (2009) The praying predator — a new arrival raises ques- tions. _https://www.nation.sc/archive/225382/island-conservation- the-praying-predator-a-new-arrival-raises-questions [Accessed March 30, 2023] Shcherbakov E, Battiston R (2020) Hierodula tenuidentata. The IUCN Red List of Threatened Species 2020: e.T118892125A118892175. https://doi. org/10.2305/IUCN.UK.2020-2.RLTS.T118892125A118892175.en 227 Shcherbakov E, Govorov V (2020) Statilia maculata (Thunberg, 1784) - the first invasive praying mantis (Mantodea, Mantidae) in the fauna of Russia. Annales de la Société entomologique de France (N.S.) 56: 189-202. https://doi.org/10.1080/00379271.2020.178 5941 Shcherbakov E, Govorov V (2021) Riders on the storm? A short note on the biology of Severinia turcomaniae (Saussure, 1872) (Mantodea: Toxo- deridae). Annales de la Société entomologique de France (N.S.) 57: 372-378. https://doi.org/10.1080/00379271.2020.1785941 Shim J, Park H, Kim S, Ju H-J, Song J-H (2021) Species delimitation of the praying mantis Hierodula patellifera (Audinet-Serville) based on mor- phological and molecular characters (Mantodea: Mantidae). Zootaxa 4951: 147-158. https://doi.org/10.11646/zootaxa.4951.1.7 Shin S, Kang D, Lee J, Seock Do M, Gu Kang H, Suh J-H, Kyung Oh H, Woo Kim T (2023) Climate warming induces the activity period pro- longation and distribution range expansion of the Asian mantis Hi- erodula patellifera in South Korea. Journal of Asia-Pacific Entomology: 102162. https://doi.org/10.1016/j.aspen.2023.102162 Sinclair JS, Lockwood JL, Hasnain S, Cassey P, Arnott SE (2020) A frame- work for predicting which non-native individuals and species will en- ter, survive, and exit human-mediated transport. Biological Invasions 22: 217-231. https://doi.org/10.1007/s10530-019-02086-7 Snyder WE, Evans EW (2006) Ecological Effects of Invasive Arthropod Generalist Predators. Annual Review of Ecology, Evolution, and Systematics 37: 95-122. https://doi.org/10.1146/annurev.ecol- sys.37.091305.110107 Snyder WE, Hurd LE (1995) Egg-hatch phenology and intraguild preda- tion between two mantid species. Oecologia 104: 496-500. https:// doi.org/10.1007/BF00341347 Soberon J, Peterson AT (2005) Interpretation of Models of Fundamental Ecological Niches and Species’ Distributional Areas. Biodiversity In- formatics 2: 1-10. https://doi.org/10.17161/bi.v2i0.4 Steger J, Schneider A, Brandl R, Hotes S (2020) Effects of projected cli- mate change on the distribution of Mantis religiosa suggest expan- sion followed by contraction. Web Ecology 20: 107-115. https://doi. org/10.5194/we-20-107-2020 Summers J, Lukas D, Logan CJ, Chen N (2023) The role of climate change and niche shifts in divergent range dynamics of a sister-species pair. Peer Community Journal 3: e23. https://doi.org/10.24072/pcjour- nal.248 van der Heyden T, Schwarz CJ (2021). New data on the presence of Hi- erodula transcaucasica Brunner von Wattenwyl, 1878 on Crete/Greece (Mantodea: Mantidae). BV news Publicaciones Cientificas 10: 6-12. https://www.biodiversidadvirtual.org/taxofoto/sites/default/files/ new_data_on_the_presence_of_hierodula_transcaucasica_brunner_ von_wattenwyl_1878_on_crete_greece.pdf Van Riel MC, Van Der Velde G, Bij De Vaate A (2009) Interference compe- tition between alien invasive gammaridean species. Biological Inva- sions 11: 2119. https://doi.org/10.1007/s10530-009-9486-4 Weaver S, McGaugh SE, Kono TJY, Macip-Rios R, Gluesenkamp AG (2022) Assessing genomic and ecological differentiation among subspecies of the rough-footed mud turtle, Kinosternon hirtipes. Journal of Hered- ity 113: 538-551. https://doi.org/10.1093/jhered/esac036 Supplementary material 1 Author: Mattia De Vivo Data type: pdf Explanation note: 52 supplementary images for the study. Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/ odbl/1.0/). The Open Database License (ODDbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://doi.org/10.3897/jor.33.111057.suppl1 JOURNAL OF ORTHOPTERA RESEARCH 2024, 33(2)