Niry T Dianzinga

Marie-Ludders Moutoussamy

Joëlle Sadeyen

Lala H R Ravaomanarivo

Enric Frago
email: enric.frago@cirad.fr

The interacting effect of habitat amount, habitat diversity and fragmentation on insect diversity along elevational gradients

Aim: Elevational gradients are a useful approach to evaluate how environmental factors affect animal diversity. Decades of studies on the elevation-diversity gradient have revealed that this gradient varies greatly with taxa and geographic regions. One potential explanation for this may be the dependence of the relationship on landscape features. We explore the impact of fragmentation, habitat diversity and habitat amount on insect diversity (alpha and beta) and abundance along elevational gradients. We hypothesize that insect diversity and abundance will relate negatively with elevation, but positively with these landscape features. We also hypothesize that landscape features will interact in a way that the positive effect of a given variable on insect diversity may be offset by the others.

Location: Reunion Island (Indian Ocean).

Taxon: The insect order thrips (Thysanoptera).

Methods: Insects were sampled along replicated elevational gradients, and at each sampling plot landscape features and abiotic variables were estimated within buffers surrounding the site. Insect alpha diversity was estimated using abundance-based rarefaction methods, whereas beta diversity was estimated calculating the "Local contributions to beta diversity" metric. The effect of elevation, rainfall, landscape features and their interactions was assessed on insect alpha and beta diversity and abundance during two consecutive seasons using linear mixed effects models.

Results:

We found that thrips alpha and beta diversity was negatively related with elevation, but the relationship varied between seasons and rainfall regimes. Among the different landscape features considered, we found that habitat amount had the strongest effect on diversity. The effect of habitat amount on diversity, however, was offset in areas of low habitat (or land cover) diversity.

Main conclusions:

Generalizing the factors that underlie the elevation diversity gradient has become a cornerstone in ecological theory because it can help to understand the impact of human activities on diversity. Here we show that taking landscape information into account may help to fulfil this objective because landscape effects co-vary with elevation with often intricate consequences for diversity.

| INTRODUC TI ON

It is nowadays widely accepted that human activities are triggering unprecedented biodiversity losses in most terrestrial ecosystems.

Human activities usually have a negative impact on biodiversity by directly degrading and transforming natural habitats, but also indirectly through changes in global climatic conditions [START_REF] Fe R E N C E S Barnosky | Has the Earth's sixth mass extinction already arrived?[END_REF][START_REF] Cardinale | Biodiversity loss and its impact on humanity[END_REF][START_REF] Hooper | A global synthesis reveals biodiversity loss as a major driver of ecosystem change[END_REF][START_REF] Urban | Accelerating extinction risk from climate change[END_REF].

By comparing diversity across space and assessing its changes in relationship with climate, relief and landscape configuration ecologists seek to understand and predict large-scale impacts of human activities on biodiversity. To achieve this, elevational gradients are useful because assemblages of species found at very different conditions can be studied within small geographic distances [START_REF] Hodkinson | Terrestrial insects along elevation gradients: Species and community responses to altitude[END_REF][START_REF] Sundqvist | Community and ecosystem responses to elevational gradients: Processes, mechanisms, and insights for global change[END_REF][START_REF] Tylianakis | Ecological networks across environmental gradients[END_REF].

Several environmental factors vary with elevation leading to lower plant species richness and structural complexity at higher elevations [START_REF] Kreft | Global patterns and determinants of vascular plant diversity[END_REF]. This elevation-diversity gradient is observed in insects because they are affected by environmental factors too, but also because they tightly depend on plants. This gradient, however, often varies with climatic regions and insect taxa [START_REF] Chamberlain | Alpine bird distributions along elevation gradients: The consistency of climate and habitat effects across geographic regions[END_REF][START_REF] Mccain | Global analysis of bird elevational diversity[END_REF][START_REF] Sundqvist | Community and ecosystem responses to elevational gradients: Processes, mechanisms, and insights for global change[END_REF]. Most studies have explored how species richness and evenness (i.e. alpha diversity) change with elevation, and there is an increasing interest in how species composition (i.e. beta diversity) changes along with these gradients. While alpha diversity expresses net diversity differences among species assemblages, through the lens of beta diversity, differences in the composition of insect assemblages can be assessed, and areas where unique species are found may be pinpointed [START_REF] Mori | β-Diversity, community assembly, and ecosystem functioning[END_REF][START_REF] Socolar | How should beta-diversity inform biodiversity conservation?[END_REF].

The elevational-diversity gradient was already noticed by Alexander von Humboldt in the 18th century, and since then it has been intensively studied with mixed and often contradictory results.

Generalizing the factors underlying this gradient across taxa and geographic regions has become a cornerstone in ecological theory, particularly because it can help to assess and to forecast the impact of human activities on diversity [START_REF] Hodkinson | Terrestrial insects along elevation gradients: Species and community responses to altitude[END_REF][START_REF] Sundqvist | Community and ecosystem responses to elevational gradients: Processes, mechanisms, and insights for global change[END_REF][START_REF] Tylianakis | Ecological networks across environmental gradients[END_REF]. The elevational-diversity gradient varies greatly with taxonomic groups and climatic regions because the impact of elevation on diversity depends on complex interactions among environmental variables and habitat and landscape characteristics [START_REF] Chamberlain | Alpine bird distributions along elevation gradients: The consistency of climate and habitat effects across geographic regions[END_REF][START_REF] Mccain | Global analysis of bird elevational diversity[END_REF][START_REF] Sundqvist | Community and ecosystem responses to elevational gradients: Processes, mechanisms, and insights for global change[END_REF]. This variation also depends on human activities, which are usually greater at lower elevations where human populations concentrate [START_REF] Luck | A review of the relationships between human population density and biodiversity[END_REF]. Landscape features have known impacts on insect diversity and are thus important to understand the elevation-diversity gradient. The most important landscape features with known impacts on animal diversity are habitat diversity, habitat amount and fragmentation. As revealed in a meta-analysis by [START_REF] Stein | Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales[END_REF], diverse or heterogeneous habitats often contain a larger number of species because they provide a larger range of abiotic conditions, available niches and refuges. These conditions ultimately allow species coexistence, reduce extinctions and promote speciation. Habitat diversity has been strongly reduced by human activities particularly due to the homogenizing effect of extensive cropping, forestry and livestock production [START_REF] Kremen | Landscapes that work for biodiversity and people[END_REF]. Despite the results found in the meta-analysis, evidence of negative or neutral effects of habitat diversity on species diversity also exist [START_REF] Bertrand | Spatial and temporal heterogeneity of the crop mosaic influences carabid beetles in agricultural landscapes[END_REF][START_REF] Marini | Waterenergy, land-cover and heterogeneity drivers of the distribution of plant species richness in a mountain region of the European Alps[END_REF]. The effect of habitat amount on diversity, on the other hand, is clearly positive. When the amount of habitat that species can exploit is reduced, species struggle to persist because resources are of less quality and less abundant, and because critical resources may be lacking [START_REF] Samways | Insect conservation: A handbook of approaches and methods[END_REF]. Habitat loss is one of the landscape features most strongly influenced by human activities, and is considered as one of the main drivers of current biodiversity loss [START_REF] Cardinale | Biodiversity loss and its impact on humanity[END_REF]Mantykapringle, Martin, & Rhodes, 2012). Fragmentation too affects species diversity, although whether the effect is positive or negative remains a controversial topic. Human activities are important drivers of fragmentation. Fragmented landscapes are often composed of small patches of natural habitat surrounded by urban, agricultural or timber-producing lands. Most authors consider that fragmentation has a negative effect on species diversity [START_REF] Fletcher | Is habitat fragmentation good for biodiversity?[END_REF] because it increases the division and isolation of natural habitats thus exposing them to disturbing human land uses [START_REF] Haddad | Habitat fragmentation and its lasting impact on Earth's ecosystems[END_REF]. [START_REF] Fahrig | Ecological responses to habitat fragmentation per se[END_REF][START_REF] Fahrig | Habitat fragmentation: A long and tangled tale[END_REF], however, suggests that since habitat fragmentation often comes together with habitat loss, when the amount of habitat available for organisms to exploit is taken into account (an estimate known as fragmentation per se [sensu [START_REF] Haila | Methodology for studying the effect of habitat fragmentation on land birds[END_REF]]), the negative effect of fragmentation is lost. This idea was proposed as the habitat amount hypothesis [START_REF] Fahrig | Rethinking patch size and isolation effects: The habitat amount hypothesis[END_REF] [but see criticism by [START_REF] Fletcher | Is habitat fragmentation good for biodiversity?[END_REF] and [START_REF] Hanski | Correspondence-Habitat fragmentation and species richness[END_REF]], where she suggests that habitat amount and fragmentation can be measured independently so that their relative contribution to diversity can be estimated. A meta-analysis that encompasses 35 studies from different locations around the globe has recently provided support for this hypothesis [START_REF] Watling | Support for the habitat amount hypothesis from a global synthesis of species density studies[END_REF].

Human impacts on natural habitats vary in intensity with elevation, but their effect on different landscape features is often intricate. For instance, while most human populations concentrate at lower elevations [START_REF] Luck | A review of the relationships between human population density and biodiversity[END_REF], food, timber and livestock production usually concentrates at particular elevational ranges, with concomitant effects on habitat diversity and fragmentation. Even if diversity usually decreases with elevation, speciation and endemism may be

K E Y W O R D S
alpha diversity, beta diversity, elevational gradient, environmental gradient, multi-scale analysis, barcoding, seasonality, insect, thrips, Thysanoptera larger at high elevations where habitats are more isolated and, in temperate areas, more fragmented due to glacial advances and retreats [START_REF] Steinbauer | Topography-driven isolation, speciation and a global increase of endemism with elevation[END_REF][START_REF] Vetaas | Distribution of vascular plant species richness and endemic richness along the Himalayan elevation gradient in Nepal[END_REF]. In addition, landscape features are usually interconnected. The positive effect of fragmentation on species diversity, for instance, may be linked to habitat diversity because fragmented landscapes are usually more heterogeneous and structurally more complex [START_REF] Fahrig | Ecological responses to habitat fragmentation per se[END_REF]. A recent theoretical model showed that habitat amount modulated fragmentation effects on diversity so that fragmentation effects were positive when the amount of habitat was large, but negative when habitat amount was reduced [START_REF] Rybicki | Habitat fragmentation and species diversity in competitive communities[END_REF].

In another example, [START_REF] Corcos | Predator and parasitoid insects along elevational gradients: Role of temperature and habitat diversity[END_REF] found that the diversity of four different insect guilds was affected both by elevation and habitat diversity, but the effect of one variable never modulated the effect of the other. Despite all these recent reports, few studies have so far measured the relative contribution of different landscape attributes to diversity [START_REF] Fahrig | Habitat fragmentation: A long and tangled tale[END_REF][START_REF] Hanski | Correspondence-Habitat fragmentation and species richness[END_REF], and how they interact with elevation. This knowledge is important to evaluate whether diversity changes along elevation are mostly driven by abiotic conditions or the local landscape. Exploring interacting effects among landscape attributes is particularly important to better understand whether the impact of habitat loss may be exacerbated or buffered by changes in habitat diversity or fragmentation. This knowledge may ultimately help better transferring biogeography theory into conservation and management practice.

In this study we explore the relative contribution of elevation, environmental variables and landscape features to Thysanoptera (i.e. thrips) diversity. Thrips are slender insects, of approximately 1 to 3 mm long with about 6.000 extant species described worldwide [START_REF] Stork | How many species of insects and other terrestrial arthropods are there on earth?[END_REF]. Thrips larvae can be herbivorous, predatory or fungivorous, and many species disperse to flowers as adults to feed on pollen or to mate [START_REF] Mound | Thyssanoptera: Diversity and interactions[END_REF]. This study was performed in Reunion, a small volcanic island located in the south-western Indian Ocean. Local-landscape characteristics were obtained by estimating habitat diversity, habitat amount and fragmentation using highly-resolved vector layers and establishing buffers around sampling points. Habitat diversity was assessed as land cover diversity [START_REF] Stein | Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales[END_REF], fragmentation as the length of the perimeter of the different polygons found within the buffer, and habitat amount as the proportion of potentially suitable habitat within the buffer [START_REF] Watling | Support for the habitat amount hypothesis from a global synthesis of species density studies[END_REF]. Since insects were sampled from flowers and during two consecutive seasons, the role of the plant, season and rainfall were also taken into consideration. The following hypotheses are tested (Figure 1): (a) We hypothesize that alpha diversity (i.e. species richness and evenness), insect abundance and beta diversity will decrease with elevation. (b) We hypothesize that habitat amount and habitat diversity will relate positively with insect diversity and abundance. Since the proportion of suitable habitat is taken into account in our study, our estimate of fragmentation can be considered as fragmentation per se, and we thus hypothesize that this feature too will relate positively with insect diversity and abundance. (c) We also explore interactions among the different landscape features, and in particular whether the relationship between habitat amount and variables related to insect diversity and abundance is modulated by either habitat diversity or fragmentation. We hypothesize that the positive effect of habitat amount on insect diversity may be offset in low diversity habitats Fragmentation, on the other hand, may have an either positive or negative effect on the relationship between habitat amount and species diversity. The interaction between habitat diversity and fragmentation may provide a better understanding F I G U R E 1 Conceptual diagram of the experimental design and variables tested in this study. Arrows point from predictor to response variables, whereas double-headed arrows indicate interactions between landscape variables. Positive and negative symbols represent our a priori hypotheses on the effect of elevation,landscape features, and their interactions on insect diversity and abundance. Larger positive symbols represent expected stronger effects on diversity and abundance, whereas question marks are shown when no directional hypothesis is proposed. A map of Reunion island with 500m contour lines and the location of transects is also shown. of whether fragmentation effects operate via changes in habitat diversity.

| MATERIAL S AND ME THODS

| Study system

In this study, adult insects belonging to the order Thysanoptera (i.e. thrips) were sampled from flowers along replicated elevational gradients that run across different habitat types. Thrips are small insects with about 6.000 extant species [START_REF] Stork | How many species of insects and other terrestrial arthropods are there on earth?[END_REF], and were sampled from flowers because even if larvae can be herbivorous, predatory or fungivorous, many species feed on pollen or disperse to flowers for mating [START_REF] Mound | Thyssanoptera: Diversity and interactions[END_REF]. This study was performed in Reunion, which is located at 55°39′ E, 21°00′ S in the south-west of the Indian Ocean between Madagascar and Mauritius with 2,512 km 2 of surface. The island has a maximum diameter of 70 km.

Its rugged terrain reaches 3,070 m.a.s.l. and 2,631 m.a.s.l. in the Piton des Neiges and Piton de la Fournaise, respectively. The island has a tropical climate with two seasons: summer and winter. The former occurs from November to April and it is rainy and hot, while the latter occurs from June to September and it is less hot and dryer.

Climatic conditions are strongly dependent on the topography and trade winds (i.e. regular ocean winds blowing from East to West).

Mean annual rainfall is higher in the eastern part (windward coast) with 1,500 mm to > 8,000 mm, whereas in the western part (leeward coast) annual rainfall varies between <500 mm to <1,500 mm.

Mean annual temperatures are of 24°C near the coastline, and they can descend to <12°C at 2,000 m.a.s.l. This island is considered a global hot-spot of diversity (IUCN, 2008) and is typified by a high level of endemism, with 46% of its 1.712 species of vascular plants endemic and with eight endemic genera. Among arthropods, it is estimated that 40% of beetles and 25% of spiders are also endemic (IUCN, 2008). Urbanization and agricultural activities are mostly concentrated below 1,000 m, while the largest proportion of natural habitats (c. 30%) is mostly concentrated above 1,200 m in areas that are less inhabited and often protected under the umbrella of the National Park [START_REF] Strasberg | An assessment of habitat diversity and transformation on La Réunion Island (Mascarene Islands, Indian Ocean) as a basis for identifying broad-scale conservation priorities[END_REF]. As in most tropical islands [START_REF] Harter | Impacts of global climate change on the floras of oceanic islands -Projections, implications and current knowledge[END_REF][START_REF] Russell | Island biodiversity in the anthropocene[END_REF], the main threats to this biodiversity are global changes, urbanization, agriculture and alien plant invasions, the latter three being particularly important at low elevations [START_REF] Strasberg | An assessment of habitat diversity and transformation on La Réunion Island (Mascarene Islands, Indian Ocean) as a basis for identifying broad-scale conservation priorities[END_REF].

| Sampling design and insect collection

Thrips were collected along seven replicated elevational transects (Figure 1). Transects were selected to capture landscape diversity and the different vegetation zones, and to ensure that landscape features varied with elevation but not collinearly. Reunion is a small island, and to avoid higher elevation points to converge, sampling took place below a mid elevation point of 1,600 m. Three transects In this study the unit of analysis was the assemblage of insects obtained from a given plant. Site was not the unit of analysis because the number of plants (and hence sampling effort) varied between sites. Plant flowers were thus microhabitats from where adults were obtained because even if adults of most species disperse to flowers to feed on pollen, or to mate, their larvae feed on other plant species, fungi, or on other arthropods [START_REF] Mound | Thyssanoptera: Diversity and interactions[END_REF]. Taking the habitat amount hypothesis as an example [START_REF] Fahrig | Rethinking patch size and isolation effects: The habitat amount hypothesis[END_REF], diversity patterns should be studied by exploring "the number of species in a plot of fixed size" and how they vary with "total habitat area in the 'local landscape' surrounding the plot". In our study, flowers represent the plots of fixed size, whereas the information extracted from buffers (see below) represents the local landscape.

| Morphological and molecular insect identification

A short description of the methods used to identify insects is provided here, but more details can be found at Appendix S2. Insects were first identified morphologically by N.T.D and a subset of specimens barcoded to ensure identifications and to detect potential cryptic species. Based on their abundance and the potential presence of cryptic species [as in Frankliniella schultzei [START_REF] Tyagi | DNA barcoding studies on thrips in India: Cryptic species and species complexes[END_REF]], a total of 223 specimens representing four known endemic species [Thrips bourbonensis, T. candidus, T. quilicii and T. reunionensis n. sp [START_REF] Goldarazena | A new species of the Genus Thrips (Thysanoptera, Thripidae) from the Afro-tropical Region[END_REF]] and six species with a global distribution (T. parvispinus, T. florum, Megalurothrips sjostedti, F. schultzei, Hercinothrips pattersoni and H. gowdeyi) were barcoded.

196 sequences were obtained and new ones imported from NCBI-GenBank (https://www.ncbi.nlm.nih.gov/genba nk/) to build a phylogenetic tree using the maximum likelihood method. All specimens, sequenced voucher specimens and their DNA are currently stored at the entomological collections of CIRAD-UMR PVBMT (CIRAD, Saint-Pierre, La Réunion) and at Centre for Biology and Management of Populations (CBGP-INRA, Montpellier/Montferrier-sur-Lez, France).

| Estimation of landscape features and environmental variables

Habitat diversity, habitat amount and fragmentation were obtained from five nested circular buffers surrounding sampling sites using a highly-resolved vector layer of Reunion Island land cover map (http:// aware.cirad.fr/layer s/geono de%3Acla ssif_gabir_2016_2017).

Vector-layer information was extracted using QGIS Desktop version 2.18.2 and then exported to R. The vector layer used contained nine habitat types: sugar cane, vegetable crops, forest, forest plantation, grassland, bare rock, savannah, orchards and urban habitats.

Habitat diversity was estimated with the Shannon index as the proportional contribution of the nine habitat types within the buffer [i.e. land cover diversity, [START_REF] Stein | Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales[END_REF]] using the function diversity from the package vegan]. Fragmentation was estimated as the length of the perimeter of the different polygons found within the buffer divided by the total surface of the buffer. Habitat amount was estimated as the proportion of the total surface of the buffer that was considered as habitat, a method recently proposed by [START_REF] Watling | Support for the habitat amount hypothesis from a global synthesis of species density studies[END_REF]. Even if the information provided by land cover data can be complemented with additional site data [START_REF] Betts | A species-centered approach for uncovering generalities in organism responses to habitat loss and fragmentation[END_REF] like plant diversity, land cover maps are ideal to capture the habitats potentially used by functionally-diverse animal groups. This is particularly true in this study given the different ecosystems studied, the abrupt landscape of the island, and the variable feeding habits of the insects studied (as explained above). Landscape features considered as non-habitats were urban areas, bare rock or the ocean (this latter case was common in sites located near the coast). Urban zones were not considered as thrips habitats because although some thrips species can be found in ornamental trees, the biomass of vegetation in these areas is usually low. If habitat amount is estimated considering urban zones as habitat, some minor differences are found, but the main results remain similar (results not shown). Environmental variables were obtained for each sampling site and season using the raster layers of average monthly data from 2007 to 2017 in Reunion Island (Météo-France, https://donne espub liques.meteo france.fr/).

The data obtained was mean, maximum and minimum annual temperature, and accumulated annual rainfall, which was extracted from raster layers using the function raster from the package raster version 2.7-15.

| Statistical analyses and multi-scale analysis

We performed all statistical analysis with R software version 3.4.2 using several packages (see R package details in Appendix S3). Alpha diversity was estimated using abundance-based rarefaction methods to estimate Hill numbers (or effective number of species) and to take into account potential undetected species of thrips [START_REF] Chao | Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies[END_REF][START_REF] Chao | Coverage-based rarefaction and extrapolation: Standardizing samples by completeness rather than size[END_REF][START_REF] Colwell | Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages[END_REF][START_REF] Jost | Entropy and diversity[END_REF]. Species richness and the Simpson index (i.e. evenness) were estimated as Hill numbers of order q = 0 and q = 2, respectively, using the iNEXT function from the package iNEXT. The package adespatial was used to assess the ecological uniqueness of each insect assemblage based on their contributions to global beta diversity, by calculating the 'local contributions to beta diversity' (LCBD) index [START_REF] Legendre | Beta diversity as the variance of community data: Dissimilarity coefficients and partitioning[END_REF]. To test the effect of elevation, landscape features, and environmental variables on diversity and abundance estimates, linear mixed effects models were built with a Gaussian error distribution using the function lmer from the package lme4. As explained above, in our study the unit of replication was the thrips obtained in a given plant. We assumed a linear relationship among variables as non-linear relationships were not observed (Appendix S4). An independent model was built for each of the following response variables: species richness (i.e. number of species), the Simpson index (or species evenness), beta diversity (the LCBD index), and insect total abundance. Given the large variation in environmental conditions during winter and summer, and that different sites were established during these two seasons, independent models were also built for the two different seasons studied. In each model plant identity, elevation, precipitation, fragmentation, habitat diversity, habitat amount, and pairwise interactions between landscape features were included as explanatory variables (Figure 1). As environmental variables, precipitation, minimum, mean and maximum temperatures were obtained, but only precipitation was retained in the models as the other variables were highly correlated with elevation (Pearson correlation R ≥ 0.95, p < 0.001, Appendix S5). Our study does not particularly focus on plant effects on thrips diversity and abundance, but this variable was included in all models as co-variable to account for this source of variability. Plants that were sampled less than five times during the whole study were considered as a single category (see details on sampled plants in Appendix S6). To account for the spatial autocorrelation between study sites, and for the non-independence of plants sampled in the same site, site nested within transect was included as a random effect. To further correct for any potential spatial autocorrelation a grid of cells of 0.09 × 0.09 decimal degrees was defined over the map. These cells were given a random identity, which was then included as a second random factor to ensure that the potential resemblance between nearby sites is accounted for.

To improve model fit and interpretation, all continuous explanatory variables were standardized to mean and standard deviation equal to zero and one, respectively [START_REF] Harrison | A brief introduction to mixed effects modelling and multi-model inference in ecology[END_REF]. In the models for insect abundance, richness and Simpson diversity, the response variable was log transformed (expressed as log[response + 0.5]) to satisfy model assumptions. In addition, the lmerControl function with the optmizer nloptwrap from the package nloptr was used to improve model performance [START_REF] Bates | Fitting linear mixed-effects models using lme4[END_REF].

Model fit was assessed by visual inspection of the residuals. In the models for insect richness and evenness, model residuals revealed slight evidence of non-homogeneity, which was considered as having negligible effects on our model estimates as recently suggested by [START_REF] Schielzeth | Robustness of linear mixed-effects models to violations of distributional assumptions[END_REF]. The absence of data points with high influence was checked by calculating Cook's distances with the function CookD from the predictmeans package. Our models did not have any significant spatial autocorrelation among data points as revealed by the function Moran-I from the ape package. p-values were obtained with a Chi-square test, using the function ANOVA from the package car. Significant results were plotted using the function visreg from the visreg package, which allows plotting partial residuals plots from mixed effects models. The non-linear multidimensional scaling (NMDS) was used to visualize how thrips diversity is influenced by elevation, rainfall and landscape features, and also to visualize how specific insect species associate with specific plants or with particular environmental or landscape features. This analysis was performed by transforming data with the Hellinger transformation [START_REF] Legendre | Distance-based redundancy analysis: Testing multispecies responses in multifactorial ecological experiments[END_REF][START_REF] Legendre | Ecologically meaningful transformations for ordination of species data[END_REF]) and using the capscale function from the package vegan version 2.5-1 .

Since we had little information on dispersal capabilities of the studied insects, and thus of the scale at which they may be affected by landscape features (i.e. the scale effect), a multi-scale analysis was also performed. As suggested by [START_REF] Fahrig | Rethinking patch size and isolation effects: The habitat amount hypothesis[END_REF], [START_REF] Jackson | Are ecologists conducting research at the optimal scale?[END_REF] and [START_REF] Watling | Support for the habitat amount hypothesis from a global synthesis of species density studies[END_REF] this analysis was done by estimating the effect of landscape features using five nested circular buffers with five different radii from the centre of each sampling site: 100, 300, 600, 1,000 and 3,000 m. These distances to the centre of the sampling site may represent short and long-distance dispersal of thrips. It is expected that the variance explained by the different landscape features in statistical models is highest for the scale that best fits the home ranges or life dispersal of the species studied. In this analysis, the variance explained by the models was obtained with the function r.squaredGLMM from the package MuMIn, and significance of a given landscape feature in the different models was tested by comparing simplified and complex models (see more details in the Appendix S12). For the models using the same response and predictor variables p-values were corrected for multiple testing using the false discovery rate method. This method is less stringent than other methods like Bonferroni, and is appropriate when a large number of comparisons are performed [START_REF] Benjamini | Controlling the false discovery rate-A practical and powerful approach to multiple testing[END_REF].

| RE SULTS

| Morphological and molecular identification of thrips

In this study, 4,280 specimens were collected and 4,261 identified to species level from flowers of 106 different plants: 86 exotic and 20 indigenous belonging to 44 families (Appendix S6). Insects were classified into 40 species, and we were unable to identify to species level 64 individuals (1.5%), which belonged to 19 different morphospecies. Seventeen of these morphospecies were identified to genus level (Appendix S7). The three known endemic species were well-represented with 1964 Thrips bourbonensis,150 T. quilicii and 114 T. candidus. A total of 329 specimens belonged to a species newly discovered in this study. This species has been named T. reunionensis n. sp [START_REF] Goldarazena | A new species of the Genus Thrips (Thysanoptera, Thripidae) from the Afro-tropical Region[END_REF], and is very similar morphologically and at the molecular level to T. bourbonensis. Some widespread (and probably exotic) species were also very abundant, including Frankliniella schultzei with 518 individuals, Haplothrips gowdeyi with 381 individuals, Megalurothrips sjostedti with 199 individuals, and Hercinothrips pattersoni with 198 individuals. Association between plants and particular thrips species can be visualized in the non-linear multidimensional scaling analysis (NMDS, Appendix S8). For example, the newly discovered T. reunionensis [START_REF] Goldarazena | A new species of the Genus Thrips (Thysanoptera, Thripidae) from the Afro-tropical Region[END_REF] was associated with the endemic plant Dombeya ficulnea, the endemic T. bourbonensis was mostly found at higher elevations, and the exotic and widespread F. schultzei was mostly found on the exotic plants Lantana camara, Solanum mauritianum and Ipomea indica. Maximum-likelihood trees (Appendix S9) mostly confirmed our morphological identifications. The endemic species T. bourbonensis, T. candidus, T. quilicii and T. reunionensis formed different clades, which were supported by 88, 100, 100 and 100 Ml bootstrap, respectively. However, two different evolutionary lineages were found for F. schultzei. Based on previous studies these two lineages were considered as different species and classified based on their colour [START_REF] Tyagi | DNA barcoding studies on thrips in India: Cryptic species and species complexes[END_REF]. These lineages are known as F. schultzei Ia1, which is close to F. schultzei from India, and F. schultzei IIa1, which is close to F. cf. schultzei/sulphurea from Australia.

| Landscape changes along elevation and multiscale analysis

Transects run from 0 to 1,600 m.a.s.l., and along this gradient important landscape changes were reported (Appendix S10 and S11).

The mean proportion of natural forest increased steadily with elevation from mean values of 25% at low elevations to more than 75% above 1,200 m.a.s.l. The proportion of urban areas followed the opposite pattern as low elevation sites were surrounded by landscapes containing between 40% and 50% of urban landscape, but this proportion decreased gradually to less than 10% at higher elevations. Although the mean proportion of sugar cane and of meadows never represented more than 25% of the area surrounding study sites, these two habitat types were important between 500 to 1,200 m.a.s.l. Savannas are common in the west of the island, and represented an average of around 25% at 200 m.a.s.l. The mean proportion of area occupied by orchards was more important at low elevations, but they never represented more than 10% of the area.

Other habitat types including forest plantations, vegetable crops or bare rock never occupied more than 5% of the area surrounding study sites (Appendix S10).

The different landscape features measured also changed with elevation without any remarkable difference when estimated using 300 and 1,000 m buffers (Appendix S11). Fragmentation changed with elevation but following a hump pattern that peaked TA B L E 1 Effect of elevation, plant of collection, accumulated rainfall, habitat diversity, fragmentation, habitat amount and pairwise interactions between habitat features on thrips alpha diversity (species richness and evenness), beta diversity (local contribution to beta diversity, LCBD) and abundance when extracting landscape data from buffers of 300m and 1000m of radius. A different linear mixed effects model was built for each buffer, and diversity and abundance estimates in both summer and winter. Significant P-values are presented in boldface type. Significant effects are graphically represented in Figure 4 as indicated near significant P-values. The total amount of variance explained by each model is also shown. Results using buffers of 100m, 600m and3000m of radius are shown in Appendix S13. .442

Richness

.503

Bold values are those that are lower than 0.05.

TA B L E 1 (Continued)

between 500 and 1,000 m.a.s.l., although fragmentation values were larger at elevations below 500 than above 1,000 m.a.s.l. The amount of habitat increased steeply from 0 to 400 m.a.s.l and then a plateau was reached, whereas habitat diversity decreased gradually with elevation. These two opposing effects probably reflect the dominance of natural forests at higher elevations, leading to a larger proportion of habitat available, but to a lower habitat diversity (or evenness).

In the multi-scale analysis, the deviance explained by the different landscape features varied greatly among landscape features and diversity metrics, and none of the buffer sizes used explained model variation greater than the others (Appendix S12). Based on this result, it was not possible to identify the buffer at which spatial effects were strongest and models were therefore built using landscape data from two different buffer sizes that may explain short and long-scale effects, 300 and 1,000 m, respectively (Table 1, but see Appendix S13 for analyses using the other three buffer sizes).

| Plant ef fec t s on insec t diver sit y an d abundance

The plant from which insects were collected was the variable with the strongest effect on insect evenness, beta diversity and abundance (Table 1, Figure 2). The plants H. ambavilla and L. camara had particularly large values of thrips evenness, D. ficulnea, H. ambavilla and Ipomoea sp. had the most different thrips assemblages (i.e. largest values of LCBD), whereas D. ficulnea, Ipomoea sp. and S. mauritianum were the plants from where the largest insect abundances were obtained. Species richness was only marginally affected by the plant in summer when using landscape data in the buffer of 1,000 m (Table 1). L. camara was the plant that contained the richest insect assemblage.

| Ef fec t of elevation an d r ainf all on insec t diver sit y an d abun dance

Of the different environmental variables considered, only rainfall was retained in the models because temperature was strongly correlated (R < 95%) with elevation (Appendix S5). In agreement with our hypotheses, elevation had a negative effect on insect alpha and beta diversity (Figure 3; Table 1). When models were built with data from 300 m buffers, the effect was observed in winter on species evenness (Table 1; Figure 4a). This effect was also observed on richness and evenness in winter when models were built with data from 1,000 m buffers (Table 1; Figure 4g andh). Elevation correlated negatively with beta diversity (LCBD) in summer in models using data from both 300 and 1,000 m buffers (Table 1, Figure 4b andi).

These results suggest that in Reunion Island thrips are richer, and the different species are more evenly represented and unique in lowlands. Rainfall had a negative effect on thrips species richness and evenness in the models built with landscape data from both 300 and 1,000 m buffers. This effect, however, was only observed in summer (Table 1; Figure 4c andj).

| Ef fec t of fr agment ation , habit at diver sit y an d habit at amount on insec t diver sit y and abundance

In winter, when including landscape data from both 300 or 1,000 m buffers, beta diversity and habitat amount related negatively (Table 1, Figure 4d andl). In the models built using data from 1,000 m buffers, insect abundance related positively with fragmentation in winter (Table 1, Figure 4k). In summer, and when using landscape data from 300 m buffers, for both species richness and evenness significant interactions between habitat amount and habitat diversity were found (Table 1). In both cases, interaction plots revealed that habitat amount had a positive effect on insect richness and evenness, but this effect was lost in areas of low habitat diversity (Figure 4e andf). This effect was not observed in winter.

| D ISCUSS I ON

In agreement with our hypotheses, insect richness and evenness decreased with elevation, although the effect was only observed in winter. In summer, accumulated monthly rainfall was the abiotic variable that dominated over thrips alpha diversity with a negative effect too. These contrasting seasonal effects may reflect the environmental conditions that dominate these two seasons. In the Indian Ocean, summer is rainy, warm and often with violent cyclones that may wash away adult thrips, an impact already observed by [START_REF] Boissot | Temporal analysis of western flower thrips (Thysanoptera: Thripidae) population dynamics on reunion island[END_REF] in this same area. In winter, however, low temperatures are likely to be the most important factor reducing alpha diversity at high elevations. As in many other insects in the tropics [START_REF] Kishimoto-Yamada | How much have we learned about seasonality in tropical insect abundance since Wolda[END_REF], abiotic variables are therefore important factors determining diversity patterns in our study. Insect beta diversity also decreased with elevation in summer, a result that we hypothesized. This effect is unlikely to be caused by the presence or absence of few influential species, or sampling effects because the effect was found in areas with particularly large insect richness. It is surprising, however, that in a hotspot of diversity like Reunion, lowland habitats that are dominated by invasive plant species (like L. camara) and urbanization [START_REF] Strasberg | An assessment of habitat diversity and transformation on La Réunion Island (Mascarene Islands, Indian Ocean) as a basis for identifying broad-scale conservation priorities[END_REF] host the most diverse and unique thrips assemblages. Many studies have already reported similar correlations between human populations and diversity [START_REF] Luck | A review of the relationships between human population density and biodiversity[END_REF] probably because human populations concentrate in areas with larger annual energy budgets that lead to more benign conditions for organisms to thrive [START_REF] Gaston | Biodiversity and extinction: Species and people[END_REF].

Thrips fauna in the Indian Ocean is poorly studied, and it is therefore difficult to know whether these lowland thrips assemblages are dominated by exotic species that reached the island together with exotic plants, or by native insects. The former is a likely situation because invasive insects often reach new territories together with their host plants [START_REF] Liebhold | Live plant imports: The major pathway for forest insect and pathogen invasions of the US[END_REF]. Knowing where native thrips are located would be particularly useful to identify hotspots of endemic diversity within the island [START_REF] Socolar | How should beta-diversity inform biodiversity conservation?[END_REF], and to unveil the role that invasive plants played in determining insect diversity patterns.

Thrips diversity was influenced by landscape features, when extracting landscape information from 300 and 1,000 m buffers. This suggests that these minute insects are able to disperse large distances as has already been reported for some species [START_REF] Mound | Thyssanoptera: Diversity and interactions[END_REF]. In summer, both species richness and evenness were significantly related to habitat amount, but this effect depended on habitat diversity. The relationship between habitat amount and insect diversity was positive in highly diverse landscapes, but this effect was offset in areas of low diversity. In agreement with our hypotheses, the negative effect of habitat loss on insect diversity may be thus exacerbated in low diversity habitats.

In a recent study, [START_REF] Corcos | Predator and parasitoid insects along elevational gradients: Role of temperature and habitat diversity[END_REF] explored whether the effect of temperature on the diversity of a range of different insect guilds was modulated by habitat diversity along elevational gradients. These authors found contrasting effects on the different insect guilds studied, but the interaction between landscape and elevation was never significant. As done by these authors, our measure of habitat diversity was based on the proportion of different habitat types within the buffer, also known as land cover diversity.

This measure mainly expresses topographic heterogeneity, and a meta-analysis by [START_REF] Stein | Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales[END_REF] revealed its strong effect on the richness of terrestrial plants and animals. In this meta-analysis, however, plant diversity and vegetation complexity also had strong impacts on diversity. Studies with cavity-nesting bees and wasps [START_REF] Staab | Tree phylogenetic diversity promotes host-parasitoid interactions[END_REF], and with insect predators [START_REF] Vehviläinen | Effects of stand tree species composition and diversity on abundance of predatory arthropods[END_REF], for instance, revealed that diversity at the plant layer can even cascade up to higher trophic levels, and that the effect may depend on plant composition and on phylogenetic diversity [reviewed in [START_REF] Moreira | Plant diversity effects on insect herbivores and their natural enemies: Current thinking, recent findings, and future directions[END_REF]]. It would be very interesting to obtain more detailed information from our study sites, and to explore how other habitat diversity measures affect the results obtained.

Relative to lowlands, highlands in Reunion had the lowest habitat diversity values because these areas are uniformly dominated by natural forests. These habitats, however, are likely to be far more diverse in terms of plant richness so the results observed here may be reversed. Going deeper into the diversity of our study sites, for example exploring plant diversity, may also help to understand the unexpected negative effect that habitat amount had on beta diversity. Exploring the effect of plant diversity may require expertise in plant taxonomy and intense fieldwork, and using highly-resolved landscape layers is likely to be the most feasible technique in areas where plant taxonomy is not fully resolved, or fieldwork is challenging. The method used here can therefore be useful to compare our results with other animal taxa because few studies have integrated a suite of landscape variables to identify and isolate the factors that drive alpha and beta diversity changes along environmental gradients [START_REF] Mori | β-Diversity, community assembly, and ecosystem functioning[END_REF].

There is a long and unsolved debate over whether fragmentation has a positive or negative effect on diversity [START_REF] Fahrig | Habitat fragmentation: A long and tangled tale[END_REF][START_REF] Fletcher | Is habitat fragmentation good for biodiversity?[END_REF], and on whether negative effects are mostly due to a reduction in the amount of habitat available. In our study, fragmentation had a positive effect on thrips total abundance in winter, but overall this variable had little influence on diversity estimates. Relative to habitat amount or habitat diversity, fragmentation is therefore a minor driver of thrips diversity in Reunion Island. In addition, if the models built are simplified and habitat amount removed (results not shown), fragmentation becomes a significant explanatory variable for diversity in some cases, which suggests that habitat amount may sometimes be the reason behind significant effects of fragmentation on diversity. Altogether, these results provide additional support to the habitat amount hypothesis [START_REF] Fahrig | Rethinking patch size and isolation effects: The habitat amount hypothesis[END_REF][START_REF] Watling | Support for the habitat amount hypothesis from a global synthesis of species density studies[END_REF], although more experimental work, and not just observational studies is needed to confirm this tendency over a range of taxa and ecosystem types.

F I G U R

F I G U R E 3

Scaterplots showing the relationship between elevation and landscape features with the different metrics of thrips diversity and abundance. Alpha diversity is expressed as species richness and evenness (Simpson index), and beta diversity as the local contribution to beta diversity metric (LCBD). Landscape features are estimated using 1000 m buffers. Black and grey dots represent values estimated in summer and winter, respectively.

Our diversity estimates were based on adult thrips collected from flowers, but their larvae are known to predate on other insects, or to feed on a variety of plants and fungi, and then to disperse as adults to feed on pollen or to mate [START_REF] Mound | Thyssanoptera: Diversity and interactions[END_REF]. Thrips are easy to sample and have a large functional and taxonomic diversity, which makes this insect order a good model to study diversity patterns. Even if some thrips species are known to be pollinators, this group is not usually considered as a pollinator guild, and most species are considered as passive dispersers with weak interactions with flowers as adults [START_REF] Mound | Thyssanoptera: Diversity and interactions[END_REF]. The multivariate analysis revealed that plants had a strong effect on thrips composition, thereby suggesting that this group can also be useful to study insect-plant interaction networks.

For example, the newly found species, T. reunionensis [START_REF] Goldarazena | A new species of the Genus Thrips (Thysanoptera, Thripidae) from the Afro-tropical Region[END_REF] was associated with the endemic plant Dombeya ficulnea, but only in winter when this plant was in full bloom.

One reason why the elevational diversity gradient in herbivores is not fully understood is that it depends on several interacting factors. Several environmental variables are hypothesized to be the driving force, but also their indirect effect via primary productivity and competitive interactions [START_REF] Kreft | Global patterns and determinants of vascular plant diversity[END_REF][START_REF] Vetaas | Principal factors controlling biodiversity along an elevation gradient: Water, energy and their interaction[END_REF]. Many hypotheses like the mid domain effect or the water-energy dynamics have been proposed. As we show here, finding general patterns may be challenging if habitat disturbance by human populations is not taken into account. Reunion is a small island in which the landscape changes dramatically with elevation because human settlements concentrate in lowlands, a large national park covers high elevation areas, and agricultural lands are found in between with strong differences between the east and the west of the island [START_REF] Strasberg | An assessment of habitat diversity and transformation on La Réunion Island (Mascarene Islands, Indian Ocean) as a basis for identifying broad-scale conservation priorities[END_REF]. Even if some landscape variables like habitat amount were collinear with elevation, elevation impacted the different habitat features measured differently. This situation offered a unique opportunity to explore how, in addition to elevation, landscape features affect diversity patterns.

Island diversity, especially in highly populated and remote islands like Reunion, is particularly threatened by habitat loss and invasive species. Given the disproportionate diversity that tropical islands host, understanding how diversity varies across space and time in these areas can help designing future conservation efforts [START_REF] Russell | Island biodiversity in the anthropocene[END_REF]. Many thrips species are pests, and studying the diversity of this taxa can also help to understand their spatio-temporal dynamics and allow the design of agricultural landscapes in which pest outbreaks are less likely.

ACK N OWLED G EM ENTS

We are grateful to Samuel Nibouche for insightful discussions F I G U R E 4 Plots showing partial residuals from the models on the effect of elevation, rainfall and landscape variables on thrips diversity and abundance when extracting landscape information from 300m (A-F) and 1000m (G-L) buffers. Only significant effects are plotted (see Table 1). The fitted lines are estimated from the linear mixed effects models, the points represent the partial residuals from the models, based on the visreg package in R, and the rugplot is shown along the abscissa axis. To fulfil model assumptions all analyses are performed after standardising all predictor continuous variables. A. The effect of elevation on thrips evenness in winter; B. The effect of elevation on thrips beta diversity (LCBD) in summer; C. The effect of accumulated rainfall on thrips species richness in summer; D. The effect of habitat amount on thrips beta diversity (LCBD) in winter. The effect of habitat amount and its interactions with habitat diversity on thrips species richness (E) and evenness (F) in summer. The interaction represents high (dark grey) and low (light grey) values of habitat diversity; The effect of elevation on thrips species richness (G) and evenness (H) in winter; I. The effect of elevation on thrips beta diversity (LCBD) in summer; J. The effect of accumulated rainfall on thrips species richness in summer; K. The effect of fragmentation on thrips abundance in winter; L.

The effect of habitat amount on thrips beta diversity (LCBD) in winter. For both the models built using landscape data from 300 and 1000m, the significant effect of accumulated rainfall on thrips evenness in summer has been omitted as the relationship is very similar to that with species richness. Raw data ranges can be seen in Figure 3.

  reached up to 1,400 m, and four to 1,600 m. Sampling sites were established every 200 m of elevation, these elevational ranges being approximate as sites were selected based on visual detection of abundant plants at the flowering stage. At each site, insects were sampled from all flowering plants in a circular area of 30 m of diameter using the beating sheet technique with a rectangular 40 × 30 cm white plastic tray and a mouth aspirator. To standardize sampling effort per plant, each was sampled for 10 min. Usually a single site per season, elevational range and transect was selected, but sometimes up to three sites per elevational range were established when few plants in bloom were found. Thrips were sampled in 2017 in summer (January to February, 64 sites) and winter (May to July, 56 sites) (Appendix S1) on sunny, dry and non-windy days from 9 a.m. to 3 p.m. Since most plants only bloom in one of the two seasons, sampling sites were not the same in winter and summer and the two seasons were therefore analysed independently. Sampling was carried out by N.T.D., M.-L.M. and E.F. To avoid any sampling bias, at each site sampling of the different plant species was randomly assigned to one of the three researchers. At each site, latitudinal and longitudinal coordinates were recorded, and elevation was obtained with the R function elevation from the package rgbif version 1.0.2.

  along the development of this project, to Gérard Lebreton and Eric Rivière (Pôle Forêt, Saint-Pierre) for help with plant identification, to Bruno Michel and Arturo Goldarazena for help with thrips identification, and to Christine Meynard, the many members of N.T.D's PhD committee, and two anonymous reviewers for insightful comments on earlier versions of the manuscript. This work was funded by the Regional Council of Reunion, the Departmental Council of the Region Reunion, the European Union (EAFRD and ERDF programs) and by the Centre de Coopération internationale en Recherche agronomique pour le Développement (CIRAD). N.T.D. was funded by Regional Council of Reunion via a PhD scholarship. EF was funded by the Regional Council of Reunion, the Departmental Council of the Region Reunion, the European Union (EAFRD), CIRAD and the Agence Nationale de la Recherche (ANR) via the PRIMA S2 2018 project INTOMED.

Evenness LCBD Abundance Chisq df p-val Chisq df p-val Chisq df p-val

																			4K							(Continues)
		p-val	p-val	<.001	.494	<.001 .569	.543 .319	.778 .695	.363 .872	.108 .745	.659 .764	.620 .694	.070 .184	.673	.222 <.001	.886	<.001 .433	.413 .215	.818 .042	.910 .438	.445 .719	.686 .281	.607 .283	.330	.056	.485
	Abundance	Chisq df	Chisq df	29.58 6	0.47 1	28.02 6 0.32 1	0.37 1 0.99 1	0.08 1 0.15 1	0.83 1 0.03 1	2.58 1 0.11 1	0.19 1 0.09 1	0.25 1 0.15 1	3.28 1	0.18 1	76.92 7	0.02 1	80.40 7 0.61 1	0.67 1 1.54 1	0.05 1 4.13 1	0.01 1 0.60 1	0.58 1 0.13 1	0.16 1 1.16 1	0.27 1 1.15 1	0.95 1	3.64 1
					4I		4B													4L		4D
		p-val		<.001	.001	<.001 .158	<.001 .428	.255 .883	.837 .435	.560 .469	.767 .471	.821 .346	.866 .415	.972	.405 <.001	.738	<.001 .772	.131 .302	.648 .093	.200 .007	.064 .919	.038 .794	.830 .733	.595	.452	.451
		df		6	1	6 1	1 1	1 1	1 1	1 1	1 1	1 1	1	1	7	1	7 1	1 1	1 1	1 1	1 1	1 1	1 1	1	1
	LCBD	Chisq		41.43	11.73	41.59 1.99	12.31 0.63	1.30 0.02	0.04 0.61	0.34 0.52	0.09 0.52	0.05 0.89	0.03	<0.01	86.65	0.11	94.23 0.08	2.28 1.06	0.21 2.81	1.64 7.18	3.43 0.01	4.31 0.07	0.05 0.12	0.28	0.56
													4F			4H		4A
		p-val		.063	.716	.089 .020	.684 .857	.026 .984	.515 .585	.401 .361	.243 .128	.597 .216	.013 .262	.998	.264 .002	.011	.002 .281	.005 .373	.296 .246	.546 .942	.832 .673	.724 .866	.567 .717	.960 .286	0.386	.281
		df		6	1	6 1	1 1	1 1	1 1	1 1	1 1	1 1	1	1	7	1	7 1	1 1	1 1	1 1	1 1	1 1	1 1	1	1
	Evenness	Chisq		11.96	0.13	10.98 5.40	0.17 0.03	4.95 <0.01	0.42 0.30	0.71 0.83	1.36 2.32	0.28 1.53	6.20	0.01	22.78	6.43	22.60 1.16	7.79 0.79	1.09 1.34	0.36 0.01	0.04 0.18	0.13 0.03	0.33 0.13	<0.01	0.75
						4J		4C					4E			4G
		p-val		.051	.920	.078 .004	.714 .464	.015 .795	.161 .566	.249 .499	.852 .050	.532 .152	.017 .321	.843	.331 .519	.044	.520 .296	.077 .891	.460 .307	.523 .279	.721 .878	.670 .766	.662 .381	.379 .097	.327	.094
		df		6	1	6 1	1 1	1 1	1 1	1 1	1 1	1 1	1	1	7	1	7 1	1 1	1 1	1 1	1 1	1 1	1 1	1	1
	Richness	Chisq		12.56	0.01	11.36 8.30	0.13 0.54	5.97 0.07	1.97 0.33	1.33 0.46	0.03 3.84	0.39 2.05	5.65	0.04	6.18	4.06	6.17 1.09	3.13 0.02	0.55 1.04	0.41 1.17	0.13 0.02	0.18 0.09	0.19 0.77	0.78	0.96
		1,000 m BUFFER	Summer	300 m BUFFER Host plant	Summer Elevation	Host plant Precipitation	Elevation Habitat diversity	Precipitation Fragmentation	Habitat diversity Habitat amount	Fragmentation Hab. div × Fragm.	Habitat amount Hab. div × Hab. am.	Hab. div × Fragm. Fragm. × Hab. am.	Hab. div × Hab. am. Total variance explained	Fragm. × Hab. am. Winter	Total variance explained Host plant	Winter Elevation	Host plant Precipitation	Elevation Habitat diversity	Precipitation Fragmentation	Habitat diversity Habitat amount	Fragmentation Hab. div × Fragm.	Habitat amount Hab. div × Hab. am.	Hab. div × Fragm. Fragm. × Hab. am.	Hab. div × Hab. am. Total variance explained	Fragm. × Hab. am.	Total variance explained

Funding information (EAFRD and ERDF programs);

DATA AVA I L A B I L I T Y S TAT E M E N T

The data sets generated during and/or analysed during the current study are available at Dryad doi:10.5061/dryad.18931zcvk. R code

S U PP O RTI N G I N FO R M ATI O N

Additional supporting information may be found online in the Supporting Information section.

How to cite this article