Climatic niche comparison between closely related trans-Palearctic species of the genus Orthocephalus (Insecta: Heteroptera: Miridae: Orthotylinae)

Previously climatic niche modelling had been studied for only a few trans-Palearctic species. It is unclear whether and to what extent those niches are different, and which climatic variables influence such a wide distribution. Here, environmental niche modelling is performed based on the Worldclim variables using Maxent for eight species of the genus Orthocephalus (Insecta: Heteroptera: Miridae: Orthotylinae). This group belongs to one of the largest insect families and it is distributed across Palearctic. Orthocephalus bivittatus, O. brevis, O. saltator and O. vittipennis are distributed across Europe and Asia; O. coriaceus, O. fulvipes, O. funestus, O. proserpinae have more limited distribution. Niche comparison using ENMTools was also undertaken to compare the niches of these species, and to test whether the niches of closely related species with trans-Palearctic distributions are more similar to each other, than to other congeners. It has been found that climatic niche models of all trans-Palearctic species under study are similar but are not identical to each other. This has been supported by niche geographic projections, climatic variables contributing to the models and variable ranges. Climatic niche models of all the trans-Palearctic Orthocephalus species are also very similar to two species having more restricted distribution (O. coriaceus, O. funestus). Results of this study suggest that trans-Palearctic distributions can have different geographic ranges and be shaped by different climatic factors.


INTRODUCTION
Environmental niches are important characteristics of species. Studying them can help to identify the environmental factors responsible for maximizing the species' fitness, and lead to a better understanding of how environment is connected to speciation and how closely related species are different in their ecological preferences. Studying climatic niches of widespread species also can help to reveal the climatic variables, connected with the species ability to adapt to different conditions. The Palearctic spans thousands of kilometers across different biomes and climatic zones. Species occupying large areas of Europe and Asia in this zoogeographic region are called ''trans-Palearctic'' and are often treated as having the same type of distribution (e.g., Abe, Melika & Stone, 2007;Hubenov, 2008;Potikha, 2015). However, the diversity of such distributions and corresponding climatic niches has never been quantitatively studied.
This current project describes the climatic niches and reveals their differences for closely related species of the Palearctic genus Orthocephalus Fieber, 1858 (Insecta: Heteroptera: Miridae: Orthotylinae). These species inhabit meadows and dry open areas, utilizing numerous species of Asteraceae, where at least the widely distributed taxa are likely to be polyphagous. Orthocephalus has been revised (Namyatova & Konstantinov, 2009), and a morphology-based phylogeny supporting its monophyly has been provided. Currently Orthocephalus includes 23 species, with O. bivittatus Fieber, 1864, O. brevis (Panzer, 1798, O. saltator (Hahn, 1835), and O. vittipennis (Herrich-Shaeffer, 1835) widely distributed in Europe and Asia (Namyatova & Konstantinov, 2009). This allows us to test whether the closely related species with wide distribution in the Palearctic occupy the same climatic niche, or whether their niches are significantly different.
Most of the Orthocephalus species have a low number (<15) of records, except for O. coriaceus (Fabricius, 1777), O. fulvipes Reuter, 1904, O. funestus Jakovlev, 1881, and O. proserpinae (Mulsant & Rey, 1852. The records of these species have also been used to build the climatic niches to find the factors influencing their distribution and compare them to those of the widely distributed species. Those niches will also allow testing whether the climatic niches of the species with trans-Palearctic distribution are more similar with each other than with those of the species having limited distribution. Comparisons of climatic niches of the Orthocephalus species will be the first step in answering the question, whether the niche conservatism or niche divergence or both can relate to speciation in this genus. The aims of the current work are to: (1) compare the niches of widely distributed species with each other and with those of species with limited distribution to determine the main climatic variables responsible for the trans-Palearctic distribution and limited distribution; (2) to test whether the trans-Palearctic species are significantly similar and whether they are more similar to each other or to the species with a limited distribution; and (3) to draw conclusions on the presence of the niche conservatism or niche divergence processes in speciation within Orthocephalus, based on the phylogeny provided in Namyatova & Konstantinov (2009).

Specimens and localities sources
Eight species (and number of unique records) have been analyzed in this work Orthocephalus bivittatus (171) Orthocephalus bivittatus is not recoded from the northern Europe and eastern Siberia, however, is common in Central Asia (Kerzhner & Josifov, 1999;Namyatova & Konstantinov, 2009). Orthocephalus coriaceus is mostly known from the middle and northern Europe with few specimens recorded from Kyrgyzstan. Orthocephalus funestus is known only from Northeast Asia. Orthocephalus proserpinae and O. fulvipes inhabit the Mediterranean region, O. fulvipes is additionally known from Arabian Peninsula and Iran (Kerzhner & Josifov, 1999;Namyatova & Konstantinov, 2009). Numerous collecting expeditions have been taken in Central Asia, Mongolia, Siberia, Russian Far East, European part of Russia by the Russian and Soviet entomologists, and their collections are mainly preserved at the Zoological Institute of the Russian Academy of Sciences (Konstantinov & Namyatova, 2019). Although many records from European countries were excluded (see below), this region is still well represented in the current analysis. Therefore, it is likely, that the known distribution for Orthocephalus species probably reflects the real distribution.

Maps
Layers in 5-arc minute (∼10 km) resolution representing different bioclimatic variables have been downloaded from Worldclim, Version 1.4 (https://www.worldclim.org/ version1). In general, the finer resolution leads to more accurate predictions, as the data are averaged within the grid cell. For this study coarser resolution was chosen for two reasons. First, it is a trade-off between the high resolution data across large geographic space and computational efficiency, as there is only a little difference between the models built using different scales especially in broadly distributed species (e.g., Araújo et al., 2005;Seo et al., 2009). Seo et al. (2009) also has shown that with a spatial grid size below 16 ×16 km, there is a good agreement among model area estimated for species of all range sizes. Second, the coordinates for many localities are approximate, and high resolution might also lead to the erroneous interpretations (Graham et al. 2005;Hanberry, 2013). Layers have been trimmed for Palearctic (20 • N-90 • N, −30 • W-180 • E) and converted to ASCII format in DIVA-GIS (https://www.diva-gis.org/). Those layers have been uploaded to QGIS 3.10 and converted to vector and used to create ''samples with data'' files (swd files).
Full-size DOI: 10.7717/peerj.10517/ fig-1 data and is considered to produce robust results with sparse, irregularly sampled data and minor location errors, which is applicable to museum data (Elith et al., 2006;Elith et al., 2011;Pearson et al., 2007;Kramer-Schadt et al., 2013). The models have been built using swd files and bioclim layers in ASCII format. For the datasets with >50 localities, bootstrap replicated run type with 25% of localities assigned for the random test percentage were applied. Overall, ten replicates were conducted. For the datasets with <50 localities (O. coriaceus, O. fulvipes and O. proserpinae), crossvalidation with the replicate number  corresponding to the locality number was used (Pearson et al., 2007;Shcheglovitova & Anderson, 2013).
The data used in this study are biased towards the easily accessed area, as most of the specimens were collected along major roads and railroads in the area currently corresponding to Russia, as well as Caucasus and Central Asia countries. This can exacerbate over-representation of some regions, which can lead to an inaccurate model. Sampling bias can be addressed by reducing the number of occurrence records in oversampled regions using spatial filtering (Kramer-Schadt et al., 2013); however, it can lead to a situation where the number of occurrences is too few to create a reliable model. Additionally, Maxent automatically discards the redundant records, appearing in the single grid cell (Fourcade et al. 2014;Chiarenza et al., 2019), which is ∼10 km in this study. This is a quite large area, which might accommodate different climate regimes. Assigning the larger grid cell was needed for the filtering, and leading to further loss of information on species climatic preferences. Alternatively, it is possible to manipulate the background data by choosing background data with the same bias as occurrences (Phillips et al., 2009;Elith et al., 2011;Kramer-Schadt et al., 2013). As soon as there are <20 occurrences for two of the analyzed species, the bias file approach was chosen. A previous study showed that using biased background data have increased the performance of the model and should be applicable for cases with small numbers of occurrence points (Kramer-Schadt et al., 2013). A bias file was created as a two-dimensional kernel density estimate, based on the coordinates of the occurrence points, using the kde2d function from the MASS package (Ripley et al., 2020) in R. This approach was applied in previous works (e.g., Filazzola, Sotomayor & Lortie, 2018;Mudereri et al., 2020). Bias files were converted to the raster ASCII format and have been implemented into the biasfile option in Maxent. Ten thousand background points, which is the default Maxent setting, was randomly selected from the area denoted in the bias file. The ''cloglog'' output was chosen for the visualization and further analysis.

Variable selection
Climatic niche modelling with two sets of variables was conducted. To avoid the model overfitting, it is possible to exclude the highly correlated variables and/or tune the model parameters (Merow, Smith & Silander Jr, 2013). It has been shown that Maxent can perform well with the correlating variables with tuned parameters (Merow, Smith & Silander Jr, 2013;Morales, Fernández & Baca-González, 2017;De Marco Junior & Nobrega, 2018); therefore, the first model type includes all bioclimatic variables (CF model). However, to test whether the model with all variables can be overfitted the highly correlated variables for each species were excluded (CR model) (see below for the details).

Parameter adjustment
There are two modifiable parameters in Maxent, which are feature classes and regularization multiplier. They should be adjusted for each particular case to avoid overfitting and/or over-complexity (Morales, Fernández & Baca-González, 2017). Feature classes correspond to the mathematical transformation of the variables and regularization multiplier (beta multiplier) limits the complexity of the model and generates a less localized prediction; i.e., smooths the model (Phillips & Dudík, 2008;Elith et al., 2011;Merow, Smith & Silander Jr, 2013). Best features and regularization multiplier set for the CF models has been selected using ENMeval package in R (Muscarella et al., 2014a;Muscarella, Kass & Galante, 2014b), using Maxent. The models have been tested against the regularization multipliers ranged from 0.5 to 6 and the default feature classes and their combinations, i.e., L, LQ, LQH, H, LQHP, LQHPT (L = linear, Q = quadratic, H = hinge, P = parameter, T = threshold). The method ''block'' was chosen, because it accounts for spatial autocorrelation (Muscarella et al., 2014a). This analysis can result in different parameter sets in different runs, so the analysis was done five times for each species. In case if the analyses resulted in different parameter sets, all of them were kept to run the environmental niche modelling in Maxent. All the parameter sets used for the modelling are provided in the Table S1.
For the CR model, the MaxentVariableSelection package in R (Jueterbock et al., 2016;Jueterbock, 2018) was used. It chooses the best set of variables, which has the lowest AICc value, based on the regularization multiplier and features. Comparisons were performed for the same parameters, as in the case of the ENMeval. A separate run was conducted to test the regularization multipliers for each feature class or combinations of classes. For the background data ten thousand background points were extracted from the bias file raster using R. For each feature (or combination of the features), its best regularization multiplier and variable set was kept for further analysis in Maxent, and they are provided in the Table  S1.
Maxent provides the list of the percent contribution (PC) and permutation importance (PI) for each variable in the model. The variables with PC and/or PI values higher than 10% are provided in Table 1 for CF and CR models for each species.

Model evaluation
For model evaluation, training and test AUC values are provided, which is valid for model comparison over the same study area (Bohl, Kass & Anderson, 2019) (see Table S1). The differences between training and test AUC values and omission error rates have been also compared. It has been shown that the model with high differences between AUC values and omission error rate > 0.1 is likely to be overfitted (Bohl, Kass & Anderson, 2019). In the case of each model type (CF and CR), the model with the relatively high AUC values, low differences between training and test AUC and low omission error rates was chosen for the visualization and niche comparisons.

Environmental niche projection area and climatic variable ranges
The obtained environmental niche models were thresholded using the ''Maximum training sensitivity plus specificity Cloglog threshold'', as the thresholds maximizing sensitivity and specificity perform well on presence only datasets (Liu, Newell & White, 2016). The total area of the thresholded niche projection was obtained using QGis 3.10. The thresholded maps were used as masks to trim the bioclim layers to obtain the climatic variables ranges for each model. These areas and variable ranges were used to compare the models. The correlation of climatic variables was estimated for each species separately using Pearson's correlations (PCor), as it is suitable for continuous variables. In this work PCor ≥ 0.9 is considered strong, as it is usually used to discriminate strongly correlated variables (e.g., Jezkova, Olah-Hemmings & Riddle, 2011;Dellicour et al., 2017). I also considered 0.7 ≤ PCor > 0.9 as significant. Tables with Pearson's correlations for the each species are provided in the Data S2.

Niche overlap
Testing for niche overlap was performed in ENMTools (Warren, Glor & Turelli, 2010). First, the niche overlap was conducted to get the Schoener's D (D) and Hellinger distance I (I) metrics, these measure similarities between species habitat suitability models. The Table 1 Variables, contributing to the models CF (first column for each species) and CR model (second column for each species). The variables used for modelling are marked with ''X''. PC and PI denote the variables having PC and PI higher than 10%. Total area of the suitable conditions projected area for each model is provided in the last raw.  values of both metrics ranged from 0 (the niches do not overlap) to 1 (the niches are identical). Comparisons were run between all species within each model type, CF and CR, separately. Second, the ''Identity test'' for each pair of species was conducted. This test randomizes the occurrences for two species, creating the pseudopopulations, and compares the environmental niches for those datasets, creating permuted D and I values. If the D and I values for the actual data are significantly lower than those of the randomized (permuted) data the niches are interpreted to be different (Warren, Glor & Turelli, 2010). The identity test can be performed only for the same set of environmental layers for both compared species; therefore, the identity test was used for the models with all environmental variables. Third, the background test was also performed. It measures the difference between the similarity of two species on one side and the similarity between species and background of another species on the other side. The test should be undertaken for two sides, as it can yield different results for the reversed comparison. If the D and I metrics for actual niche overlap and obtained with background test are similar, this means that the similarity of niches between two species is the same as expected from random data. If the D and I metrics of actual datasets are higher or lower than those from background test, this means that the niches are more similar or more different than expected from random data respectively (Warren, Glor & Turelli, 2010). The background for each species equals to its bias file.

Maps visualization
All maps have been prepared in QGis 10. The background for the Figs. 1-2 is the layer freely accessible at https://github.com/nvkelso/natural-earth-quick-start/tree/master/ 50m_raster/NE1_50M_SR_W, and it is not copyrighted. The maps for Figs. 3-6 have been created by the uploading the averaged maps resulted from the Maxent analysis to QGis 10. The country borders layer is freely accessible at https://github.com/petewarden/ openheatmap/tree/master/mapfileprocess/test_data/TM_WORLD_BORDERS-0.3, and it is not copyrighted.

Model evaluation
All the Maxent models have high discriminative power for the training datasets with high AUC. It is higher than 0.9 in all cases except for O. vittipennis, where AUC ranges vary from 0.87 to 0.89. The models are also able to predict the testing points with very similar AUC values, as in training datasets. The training AUC is higher than the test AUC, and the differences between the models chosen for the comparison vary from 0.002 ( The model with all variables is supposed to be overfitted because of the correlated variables. Therefore, this model is expected to predict smaller areas of suitable conditions Those areas of CF models are larger than or subequal to CR for all the species. In the case of variables ranges no model type is noticeably more restrictive than the other (Figs. 7-11).
In rare cases, the CF models show significantly more restricted ranges, rather than CR models (e.g., bio12, bio13 for O. bivittatus and O. proserpinae, bio18 for O. proserpinae).
The variables ranges for each model are shown on the Figs. 7-11. The climatic variables with high PC and PI for each model type and for each species are provided in the Table 1. Both models for the same species have different sets of climatic variables explaining their distribution, and each type of model (CF or CR) has different sets of climatic variables explaining the distribution in comparison between the species.
In most of the models both temperature related and precipitation related variables, are important for the species distribution, except for the CF models for O. saltator and O. vittipennis, having only temperature related variables significantly contributing. Most of the variables appear as important for at least one model, except for bio5 (max temperature of warmest month), bio8 (mean temperature of wettest quarter), bio12 (annual precipitation) and bio13 (precipitation of wettest month). In some cases, the same variable significantly contributes to the both type of models within the same species, i.e., bio1 (  Comparison of the variables for the species with similar environmental niches  annual range (bio7) and mean temperature of coldest quarter (bio11) highly contribute to at least one of the models in O. bivittatus, O. saltator and O. vittipennis, and they also significantly correlate with each other in all those species, as well as with bio1, bio6 and bio9 for most of the species. In the models of two species, widely occurring in Central Asia (O. bivittatus and O. vittipennis), mean diurnal range (bio2), bio9, mean temperature of coldest quarter (bio11) and precipitation of warmest quarter (bio18) significantly contribute to at least one of the models. In O. bivittatus and O. brevis precipitation seasonality (bio15) is important for at least one of the models. In O. brevis precipitation of driest month (bio14) significantly contributes to both models. Max temperature of warmest month (bio5), mean temperature of wettest quarter (bio8), annual precipitation (bio12), precipitation of wettest month (bio13), precipitation of driest quarter (bio17) only slightly contribute or do not contribute to the climatic models of the trans-Palearctic species. Suitable conditions for O. bivittatus are shifted to drier places than in other species, whereas suitable conditions for O. brevis are predicted for the places with higher precipitation than in other species. In contrast to other species, suitable conditions for O. vittipennis are predicted in areas with very low temperatures over the winter and very strong seasonality. The models for O. saltator are similar to O. brevis in precipitation levels and temperature changes around the year; however, suitable conditions of the former are also predicted for the areas with warmer temperatures over the summer, than in O. brevis.
The models of O. coriaceus are more similar than random with all the models of widespread species (Table 2). It is different from all of them in the lower margins for isothermality (bio3) limited with higher values (Fig. 7), and suitable conditions are predicted for the places with low temperatures over summer (bio5, bio10) (Figs. 8 and 9  (Figs. 10 and 11). The models of O. coriaceus differ from O. saltator in the upper margin of the precipitation of the driest month and quarter (bio14, bio17) range limited with higher values (Figs. 10 and 11), and the upper margin of seasonality (bio4) limited with lower values (Fig. 7).
The models for the northeastern O. funestus are more similar than random with O. brevis, O. saltator and O. vittipennis, and they are most similar to O. brevis and O. vittipennis. The models of O. funestus are different from those of abovementioned three widespread species in isothermality (bio3) range very narrow and shifted towards lower values (Fig. 7). In contrast to all other species, suitable conditions for O. funestus and O. vittipennis are modelled for the places with very low temperatures of coldest month, coldest and driest quarters (bio6, bio9, bio11) and strong seasonality (bio4) . The upper margins of the diurnal range (bio2), isothermality (bio3), temperature seasonality (bio4), temperature annual range (bio7) are limited with the lower values (Figs. 7-8), and the upper margins are limited with higher values for many precipitation variables (bio12-14, 16, 17-19) in O. funestus in comparison with O. vittipennis (Figs. 9-11).
Both O. fulvipes and O. proserpinae inhabit southern areas of European Palearctic (Figs. 1D, 2B), but the variables contributing to their models are different (Table 1). For both species in the CF model the mean temperature of driest quarter (bio9) is important. Precipitation over wettest or coldest quarters (bio16 and bio19) significantly contribute to the models of both species. They also have either precipitation of driest month (bio14) or precipitation of warmest quarter (bio18) with high PC. They have similar ranges for isothermality (bio3) and temperature seasonality (bio4), as well as higher temperatures of driest and coldest periods (boi6, bio9, bio11), and lower precipitations over the driest and warmest periods (bio14, bio17, bio18) (Figs. 7-11). The models of O. fulvipes differ from those of O. proserpinae in many temperature variables limited with higher values (bio5, bio6-bio11) (Figs. 8-9), as well as mean diurnal range and temperature seasonality and precipitation of driest month (bio2, bio4, bio14) (Figs. 7 and 10) .

Niche overlap, identity test and background test
The results for the niche overlap, identity test and background test are shown in the  model and statistical metric. The results also can differ for the pair of species, depending which species is used for the background. Generally, the CR models show more overlap with each other than the CF models, and background test is more often shows that the species are more similar to each other for the CR models. If two widespread species are compared, the background test results in very similar metrics values for both directions. If two species with very different areas of suitable conditions sizes are compared, the results depend on which of them is used for the background. In case when the species with larger distribution area is used as a background, the resulted metrics are lower, and therefore, the analysis shows that two species are more similar to each other than expected more often, rather than in the reversed comparison.

Climatic niches of Orthocephalus species comparison
The modelled environmental niches for the widespread species cover noticeably different areas. Areas in which conditions are suitable for O. bivittatus correspond with the other species the least. They are mostly projected on Central Asia and south of European part of Russia. (Figs. 3A, 4A). Orthocephalus brevis and O. saltator are more similar, having  (Figs. 5D, 6D), whereas the Mediterranean area is not suitable for this species. Those differences in the areas of suitable conditions are also supported by the differences in the variables with highest contribution to the climatic models ( Table  1) and comparisons of climatic variable ranges (see 'Results' for the details). Overall, all analyzed Orthocephalus species have different set of variables, most important for their climatic models and none of the pairs have identical climatic niches, which suggests that the climatic niche is species specific even for the closely related taxa. Finally, the differences in the environmental niches are also supported by the identity and background tests. Although in all widespread species niche overlap is high (I > 0.8, D > 0.6) and background test shows that the similarity between them is higher than that of random data, the identity test does not support the hypothesis that they are identical (Table 2). This supports the previous research, showing that the environmental niches in closely related species might be similar, but not identical (e.g., Wellenreuther, Larson & Svensson, 2012;López-Alvarez et al., 2015;Dellicour et al., 2017).
The climatic niche models of trans-Palearctic species are all more similar to each other than expected. The niches of O. coriaceus and O. funestus, distributed in Europe and Northeast Asia respectively, also highly overlap with some trans-Palearctic species. The climatic niche models of O. fulvipes and O. proserpinae, are not very similar with those of trans-Palearctic species. Schmitt (2007) made an overview of the European types of distribution, and delimited three main types. Species with the center of dispersal in Mediterranean regions are ''Mediterranean'', species having extra-Mediterranean center of dispersal belong to ''Continental'' type, and species with recent alpine or arctic distribution patterns are called ''Alpine'' or ''Arctic''. Testing the center of the distribution for the Orthocephalus species is beyond the scope of the current paper: a phylogeography analysis combined with modelling of distribution during the glacial period is necessary for such a study. However, it is very likely that the trans-Palearctic Orthocephalus species belong at least to two Schmitt's groupings (2007). Among the analyzed widespread species, suitable conditions in the Mediterranean region are generally not predicted for O. bivittatus, O. brevis and O. vittipennis and they probably belong to the ''Continental'' group. However, O. saltator has highly suitable climatic conditions in the Mediterranean region and might have its center of dispersal there (Figs. 5C, 6C).
Maps with the records and the areas of suitable conditions suggests that ''trans-Palearctic'' distribution is a term of convenience and might include many types of distributions connected with the different climatic conditions. Study of more species from different taxonomic groups, including closely related and unrelated species, are necessary to make conclusions on which types of environmental niches are suitable for the trans-Palearctic species and which climatic conditions are connected with such distributions.

Climatic variables important for the Orthocephalus distribution
Annual mean temperature (bio1), temperature annual range (bio7) and mean temperature of coldest quarter (bio11) are important for the models of at least three widely distributed Orthocephalus species. The results also correspond with the results of environment niche modelling of the trans-Palearctic beetle species Pterostichus oblongopunctatus (Fabricius, 1787), which showed that bio1 and bio11 had high contribution to the model of this species (Avtaeva et al. 2019). Modelling of the trans-Palearctic tick Ixodes ricinus (Linnaeus, 1758) showed that min temperature of coldest month (bio6) is among the variables with the highest PC in the climatic model of this species (Porretta et al., 2013) In all trans-Palearctic species bio6 and bio11 highly correlate with bio1 (PCor > 0.9), and this might mean that the distribution of the trans-Palearctic arthropods is limited at least partly by the winter temperatures. Overall, the set of the variables having high contribution to the climate models is unique for each Orthocephalus species, and analysis of more trans-Palearctic species is needed to make conclusions on how type of distribution correlates with the climatic variables.
Orthocephalus brevis, O. coriaceus, P. obongopunctatus and I. ricinus are mostly distributed in the Western Palearctic, and their climatic models have precipitation of driest month (bio14) or precipitation of driest quarter (bio17) with high PC. This reflects the fact that the most suitable conditions for all those species cover middle and northern Europe, but not Mediterranean zones, and therefore prefer the places with relatively high precipitation around the year. This is also supported by the variable ranges of bio14 and bio17 variables for O. brevis and O. coriaceus, which show that these species inhabit places with high precipitation values over the driest period. Similarly to four mentioned species, Orthocephalus saltator is mostly distributed in Western Palearctic, but the areas with the most suitable conditions cover Mediterranean zones. In case of this species, bio14 and bio17 do not contribute much to the climatic model (<10% of PC and PI).
Although O. coriaceus and O. funestus have very different distribution, isothermality (bio3) has high PC and PI for both models in both species. However, according to the temperature ranges the former tends to inhabit places with high isothermality, whereas the latter prefers the areas with low values of this variable. The study with more European and East Asian species and species from Northeast Asia is needed to confirm that isothermality is connected with such distributions.
The precipitation related variables over the different seasons (bio14, bio15, bio18, bio19) are important for models for two species mostly distributed in the Mediterranean regions, and mean temperature of driest quarter (bio9) has high PC and PI in the CF model. This is similar to the Mediterranean species Tomicus destruens (Wollaston, 1865) (Sánchez-García, Galián & Gallego, 2015). It was shown that the variables bio19 and bio9 are also important for the climatic models of different clades and haplotypes of this species. Study of climatic preferences of more insect species from the Mediterranean region is needed to draw conclusions on how the climate is connected with the species distribution in this region.

Comparison of environmental niches between closely related Orthocephalus species in the phylogenetic context
Although, a morphology-based phylogeny for Orthocephalus has been published (Namyatova & Konstantinov, 2009), there are many unresolved clades, and, therefore, cannot be used to adequately analyze potential phylogenetic signal in climatic related tolerances and environmental niches. However, a few conclusions still can be reached. Based on these relationships, it can hypothesized that there are at least three processes in this genus related to climatic niches which might be at play. First, there might be phylogenetic conservatism, at least for some climatic variables. For example, O. funestus and O. vittipennis are very similar morphologically. They mostly differ only in hemelytron coloration in males. According to the current analysis, their niches are also more similar to each other than to random data, and both those species can tolerate strong seasonality and very low winter temperatures. Another example is O. fulvipes, inhabiting southwestern Palearctic, which morphologically is very similar to O. saltator, and the latter is the only widespread species in which the Mediterranean region is suitable.
Second, climatic niche convergence is also observed. This study found that climatic niches are very similar in O. brevis and O. saltator, as well as O. brevis and O. vittipennis, and those pairs are not closely related.
Third, the analysis shows the possibility of the distinct niche divergence in the distantly related taxa. According to the background test, the niches are undoubtedly more different from each other than from random data only for the O. funestus and O. proserpinae pair. According to Namyatova & Konstantinov (2009), those two species are not closely related. Those significant niche differences might be explained by differences in habitats, occupied by those two species. Orthocephalus funestus inhabits places in Northeast Asia with wide ranges of temperatures and high precipitation, whereas O. proserpinae lives in the Mediterranean region with narrow temperature ranges and low precipitation (Figs. 5B, 5D, 6B, 6D). However, O. fulvipes, which is more closely related to O. funestus, also prefers dry conditions including Mediterranean and desert climates (Figs. 3D, 4D). The background test for this species pair provides inconsistent results, suggesting that the niches between O. funestus and O. fulvipes are more similar with each other, than those of O. funestus and O. prosepinae. This leads to the hypothesis that the closely related species in Orthocephalus cannot diverge very quickly.
However, the robust molecular-based phylogeny and niche models for other Orthocephalus species are needed to test all those hypotheses on the niche evolution in this genus.
Based on this phylogeny, it is unclear whether the ability for the wide distribution has phylogenetic signal. On one hand, it is very likely that two closely related species can similarly adapt to the climatic conditions (e.g., Losos, 2008;Wiens et al., 2010). On the other hand, even though both sister species can potentially tolerate wide range of climatic conditions, one species might have significantly limited realized niche and distribution because of the strong competition with its sister species.

CONCLUSIONS
The study on the climatic niche modelling for eight insect species with trans-Palearctic distribution from the genus Orthocephalus has been performed. The niches of widely distributed trans-Palearctic species (O. bivittatus, O. brevis, O. saltator, O. vittipennis) are very similar to each other, but not identical. The differences are confirmed by the ''Identity test'' and ''Background test'' in ENMTools, as well as the comparison of the climatic variables contributing to the modes and variable ranges for the areas, covered by preferable conditions of different species. The niches of the trans-Palearctic species are also similar to two species having more limited distribution (O. coriaceus, O. funestus). Overall, the similarity of the niches of widely distributed species is higher than in cases when the niches of widely distributed and locally distributed species or only locally distributed species are compared. The annual mean temperature significantly contributes (bio1) to the climatic models of all trans-Palearctic species. Other temperature related variables, i.e., min temperature of coldest month (bio6), temperature annual range (bio7), mean temperature of driest quarter (bio9), mean temperature of coldest quarter (bio11) are likely to be important for the climatic niches of the trans-Palearctic species. For the trans-Palearctic species, widely distributed in Central Asia, mean diurnal range (bio2), min temperature of coldest quarter (bio11) and precipitation of warmest quarter (bio18) are also important for at least one of the models. For the Western Palearctic species with the most suitable conditions corresponding to the areas outside the Mediterranean regions (O. brevis and O. coriaceus), precipitation of driest month (bio14) is important. Isothermality (bio3) has high PC and PI for European O. coriaceus and East Asian O. funestus. For the species, distributed in the Mediterranean region (O. fulvipes and O. proserpinae), precipitation related variables over different seasons (bio14, bio16, bio18, bio19), significantly contribute to at least one of the model each. The discussion of the results based on the phylogeny suggests that within Orthocephalus there might be different processes connected to the climate niche differentiation, such as niche conservatism, niche convergence and niche divergence. More studies of climatic niches of the species distributed in the Palearctic are needed to better understand the types of possible climatic niches of widespread species, the main climatic variables shaping the distribution of the Palearctic taxa and how the climatic niches are related to phylogenetic history.

ADDITIONAL INFORMATION AND DECLARATIONS Funding
This work was supported by the Russian Science Foundation grant (no 19-74-00077). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Grant Disclosures
The following grant information was disclosed by the author: Russian Science Foundation grant: 19-74-00077.