﻿Seasonal and microclimatic effects on leaf beetles (Coleoptera, Chrysomelidae) in a tropical forest fragment in northeastern Mexico

﻿Abstract Leaf beetles (Coleoptera: Chrysomelidae) constitute a family of abundant, diverse, and ecologically important herbivorous insects, due to their high specificity with host plants, a close association with vegetation and a great sensitivity to microclimatic variation (factors that are modified gradually during the rainy and dry seasons). Therefore, the effects of seasonality (rainy and dry seasons) and microclimate on the community attributes of chrysomelids were evaluated in a semideciduous tropical forest fragment of northeastern Mexico. Monthly sampling was conducted, between March 2016 and February 2017, with an entomological sweep net in 18 plots of 20 × 20 m, randomly distributed from 320 to 480 m a.s.l. Seven microclimatic variables were simultaneously recorded during each of the samplings, using a portable weather station. In total, 216 samples were collected at the end of the study, of which 2,103 specimens, six subfamilies, 46 genera, and 71 species were obtained. The subfamily Galerucinae had the highest number of specimens and species in the study area, followed by Cassidinae. Seasonality caused significant changes in the abundance and number of leaf beetle species: highest richness was recorded in the rainy season, with 60 species, while the highest diversity (lowest dominance and highest H’ index) was obtained in the dry season. Seasonal inventory completeness of leaf beetles approached (rainy season) or was higher (dry season) than 70%, while the faunistic similarity between seasons was 0.63%. The outlying mean index was significant in both seasons; of the seven microclimatic variables analyzed, only temperature, heat index, evapotranspiration and wind speed were significantly related to changes in abundance of Chrysomelidae. Association between microclimate and leaf beetles was higher in the dry season, with a difference in the value of importance of the abiotic variables. The results indicated that each species exhibited a different response pattern to the microclimate, depending on the season, which suggests that the species may exhibit modifications in their niche requirements according to abiotic conditions. However, the investigations must be replicated in other regions, in order to obtain a better characterization of the seasonal and microclimatic influence on the family Chrysomelidae.


Introduction
Accelerated loss of biological diversity, as well as the alterations in native ecosystems as a result of human activities, are among the most important environmental issues at a global level (Challenger and Dirzo 2009). These include land cover fragmentation, overexploitation of natural resources, pollution, and climate change (Hautier et al. 2015).
Abiotic modification produces direct effects on organisms, affecting physiology, behavior, and reproduction (Uribe-Botero 2015). Changes in precipitation and increased environmental temperature (Schaefer et al. 2008) are likely to cause alterations in abundance and even loss of species (Brook et al. 2008), as well as changes in their geographical distribution (Parmesan and Yohe 2003;Root et al. 2003). However, these responses are variable, based on the type of organism and its niche breadth (Vié et al. 2009). Therefore, changes in climatic abiotic variables are key factors in the composition and structure of biological communities, besides other ecological aspects (Pimm 2007), such as seasonal changes during wet and dry seasons (Wolda 1988;Rzedowski 2006).
An aspect of greatest influence on these communities is the microclimate (Cloudsley-Thompson 1962). This is the result of local spatial and seasonal variations in climate and has been shown to play an important role in the dynamics of metapopulations (Checa et al. 2014). Likewise, it is essential for the survival and development of the species, affecting larval diapause or growth, or indirectly modifying the availability of food resources (Currano et al. 2008;DeLucia et al. 2008). The microclimate is related to seasonal variations in the communities of phytophagous insects (Chen et al. 1999), but its specific influence has been scarcely studied.
Phytophagous insects are among the most important trophic groups that respond significantly to climatic changes. Their presence is key in natural or anthropic ecosystems, either playing a relevant role in nutrient cycling processes, or in the diet of other organisms (Iannacone and Alvariño 2006). Furthermore, their physiological processes are determined by the conditions of the environment (Régnière 2009).
Leaf beetles (Coleoptera: Chrysomelidae) constitute a model family to evaluate the seasonal effects of abiotic variation on herbivorous insect communities, since they occupy one of the first places in worldwide diversity (Santiago-Blay 1994). Most chrysomelid species exhibit phytophagous feeding habits and a close relationship with their host plants, as well as a great sensitivity to microclimatic variation . Also, they are considered to be a group with important potential for monitoring natural areas (Furth et al. 2003).
The present study was carried out in a semideciduous tropical forest (STF) fragment in the municipality of Victoria, Tamaulipas, in northeastern Mexico. The area is included in the biogeographic province of the Sierra Madre Oriental and is located within one of the 15 panbiogeographic nodes of Mexico (Morrone and Márquez 2008). Therefore, it constitutes a region with a high priority for conservation (CONABIO et al. 2007). Despite this, there are no studies in STF evaluating the effect of seasonality and microclimate on the family Chrysomelidae. It is important to recognize the factors that restrict the distribution of the species, and thus further delimit efficient conservation strategies of this important area. Based on the above, the objectives of this study were 1) to prepare a faunistic list of chrysomelid species, 2) to compare their richness, abundance, and diversity between seasons, 3) to define the abiotic variables seasonally related to the presence and abundance of the species, and 4) to delimit the breadth niche and categorize the leaf beetles as specialists or generalists, based on their variation related to the seasonal abiotic environment.

Study area
The study area of semideciduous tropical forest (STF) is located in the Ejido Santa Ana, municipality of Victoria, in the center of the state of Tamaulipas,northeastern Mexico 23°52'4.27"N,99°13'51.37"W and 23° 47'23.06"N,99°18'10.22"W (DMS) (Fig. 1). It is included in the biogeographic province of the Sierra Madre Oriental (SMO), converging to the south with Peregrina Canyon, within the Natural Protected Area (NPA) "Altas Cumbres." Two climate groups characteristic of Tamaulipas were observed in the area: 1) Semiwarm, sub-humid, with summer rains, averaging temperatures between 16.4 °C and 29.2 °C, and 2) Semi-warm, semi-dry subtype, with average temperatures from 15.1 °C to 22.9 °C. The average annual precipitation is 577 mm, with May to October as the wettest months (rainy season) and November to April having the lowest precipitation (dry season) (Gobierno del Estado 2015). Regarding the semideciduous tropical forest, it is the second richest ecosystem in plant species of the state of Tamaulipas and is located between 350 and 500 m a.s.l, comprising areas adjacent to the margin of rivers and streams. Therefore, this habitat conserves higher environmental humidity for most of the year, protecting it from sudden climatic changes, such as sudden temperature differences (García-Morales et al. 2014).

Sampling
A total of 18 plots measuring 20 × 20 m (400 m 2 ) was randomly established over an approximate land area of 5 km 2 . Plots were distributed in areas of dense herbaceous and shrub vegetation, separated at least by 10 meters from the main road, in order to minimize anthropogenic influence. Each plot was measured and delimited with a 50 m tape, using trunks, trees, or branches as vertices; the center of the plot was georeferenced with a Garmin Etrex 30 GPS and then marked with a brightly colored ribbon to facilitate its location in the field.
Beetles were sampled with an entomological sweep net of 60 cm length and 40 cm rim diameter. In each plot (sample unit), 200 net beats were made, covering all the sampling area zigzagging on the understory vegetation. The contents of the net were placed inside a polyethylene bag with 70% alcohol and a collecting label data. All of the 18 plots were sampled from 10:00 to 17:00 hours, once a month, from March 2016 to February 2017.
Sample bags were processed in the Entomology Laboratory of the Facultad de Ingeniería y Ciencias, Universidad Autónoma de Tamaulipas. Each sample was placed in a tray with water, plant debris were then removed using entomological forceps, and the insect specimens were afterwards placed in small bottles with 70% alcohol. Later, the contents of each bottle were analyzed in a Petri dish, using a stereoscopic microscope to identify the specimens; chrysomelids were dried on absorbent paper and mounted in opaline triangles, following the methodology of Triplehorn and Johnson (2005). Taxonomic determination of subfamilies was carried out using the keys of Triplehorn and Johnson (2005), while genera and/or species were identified by consulting various authors (Wilcox 1972;Scherer 1983;White 1993;Flowers 1996;Riley et al. 2002;Staines 2002), as well as by comparison with previously identified specimens.
Microclimatic variables were recorded using a Kestrel 3500 portable meteorological station, with which the following variables were evaluated: maximum wind speed (m/s), average wind speed (m/s), temperature (°C), relative humidity (%), heat index (°C), dew point (°C) and evapotranspiration (°C). Abiotic data collection was carried out in each plot, simultaneously with the sampling of leaf beetles (once a month for each plot, during the period from March 2016 to February 2017).

Data analysis
Statistical differences in abundance and number of species between seasons were calculated with a non-parametric Mann-Whitney test and a diversity permutation test, respectively. Both analyses were conducted using PAST 3.17 software (Hammer et al. 2017).
Seasonal estimated richness was determined using Chao 1, Chao 2, Jackknife 1 and ACE non-parametric estimators. These indices are recommended for the minimum estimate of richness and useful as a complementary measure in biodiversity analyzes (Gotelli and Colwell 2011). Chao 1 considers the abundance of rare species (singletons and doubletons). Chao 2 is robust for presence-absence data. Jackknife 1 is a conservative index based on incidence data of those species found only in a single sample, while ACE is an index that considers the abundance of species represented by 1-10 individuals (Magurran 2004). The estimators were calculated by means of 100 randomizations without replacement in the software EstimateS 9.1.0 (Colwell 2013), based on the abundance of the recorded species. In addition, the Clench model was used to calculate the estimated species richness, following the methods proposed by Jiménez-Valverde and Hortal (2003). This procedure was performed in STATISTICA 8.0 (StatSoft, Inc. 2007).
Alpha diversity was estimated using Shannon's entropy index (H') and Simpson's dominance index (D). Both values were transformed to the effective number of species (true diversity), through the Hill numbers of order (q) 1 and 2, respectively (Jost 2006). To measure beta diversity, the Bray-Curtis similarity index was used, which relates the abundance of the shared species with the total abundance in two samples. Therefore, it constitutes a robust measure for the analysis of biotic similarity between communities (Magurran 2004). All diversity analyses were carried out with PAST software.
Association between leaf beetle species and the environmental abiotic variables, as well as the measure of niche breadth, were calculated with the Outlying Mean Index (OMI). This index identifies the niche of the species, or marginality, according to the average distance between the abiotic resources used by each species (centroid) with respect to the total resources available (microclimate) in the area. It gives a more even weight to all sampling units, including those with a low number of species or individuals (Dolédec et al. 2000). First, the OMI assesses the contribution of the abiotic variables to the niche separation of the species by computing a Principal Component Analysis, and higher correlation values (loadings) are interpreted at each of the most important axes. Then, a total Inertia (InerO) value is obtained, which is a measure proportional to the average marginality of the species and represents a quantification of the influence of environmental variables on the separation of the species niche. Lastly, the analysis decomposes the inertia associated with the distribution of a species (InerO) into three main parameters: Marginality, Tolerance (T1), and Residual Tolerance (T2) (Dolédec et al. 2000).
Marginality represents the deviation of the environmental conditions used by a species with respect to the average environment for the entire study area. Species with high OMI values have marginal niches (occur in atypical habitats, and are influenced by a specific subset of environment variables), while those with low values have non-marginal niches (common species occurring in typical habitats, without a specific response to environment variables). Tolerance (T1) measures the dispersion of the assessment units that contain a species along an environmental gradient (the range of habitat of the species), and it is analogous to the concept of niche breadth: high tolerance values represent greater niche breadth, and the species are distributed in habitats with widely variable conditions (generalist); contrarily, low tolerance values indicate a smaller niche width where a species is distributed in habitats with a limited range of conditions (specialists). Finally, T2 is defined as the variance in the species niche that is not considered by the marginality axes, and it is useful for determining the reliability of a set of environmental conditions for the definition of the niche of each species (Dolédec et al. 2000).
Statistical significance of the OMI was determined with a Monte Carlo test, in which the observed marginalities are compared with 10,000 random permutations, in order to reject the null hypothesis that species are equally distributed in relation to (not influenced by) environmental variables (Dolédec et al. 2000). All OMI analyses were carried out in ADE-4 software (Thioulouse et al. 1997), and they were calculated separately for the rainy and dry seasons. Data input consisted of a matrix with the abundances of each of the species in each month/season and a matrix with the values of the seven environmental variables registered in each of the sampling plots. Ordination graphics of centroids and loadings were generated in the same software and later exported to CorelDRAW X3 to be edited. Environmental ranges of species were calculated for each of the significant variables using the Kriging interpolation technique, which is a geostatistical method that quantifies spatial autocorrelation for the prediction and generation of continuous surfaces (Murillo et al. 2012). Procedures were carried out in ArcGis 10.2.2 (ESRI 2014).

Seasonal variation
Seasonal differences in abundance of the leaf beetle community were statistically significant (Mann-Whitney U = 4039; p ≤ 0.0001). The highest number of specimens was recorded during the rainy season (1,242 specimens, involving 41 genera), followed by the dry season (861, involving 30 genera). According to the permutation test, significant differences were also found in the number of species and diversity. Highest species richness was recorded in the rainy season. In contrast, the lowest dominance and highest diversity were obtained in the dry season (Table 1).
Estimated species richness according to the non-parametrical estimators in the rainy season ranged between 85 and 100 species; therefore, the observed richness represents between 59.66 and 69.96% of completeness. For the dry season, the estimated richness varied from 48 to 56 species, indicating a completeness from 70.49 to 82.85% (Table 2). Inventory reliability with Clench's model was higher during the dry season, with a completeness of 81% and a lower slope value, compared with the rainy season ( Table 2).
The best represented subfamily during the rainy season was Galerucinae (943 specimens, 32 species), followed by Cassidinae (260, 11 species). This same pattern was reflected in the dry season: Galerucinae with 685 specimens (21 species), followed by Cassidinae with 150 specimens (7 species). The remainder of the subfamilies had lower abundances and number of species for both seasons (Table 3).

Response of Chrysomelidae to seasonal microclimatic variation
The OMI analysis for the rainy season indicated a significant deviation between the abiotic conditions used by the leaf beetles and the average total microclimatic conditions (Monte Carlo test, p = 0.047). Of the 60 species registered in this season, only six showed a significant association. Centralaphthona diversa and M. bumeliae obtained low marginality values, which represents a wider niche breadth, and they were thus considered to be generalist species (Table 4); abundance of these species was equally distributed in almost all samples (Fig. 2). The rest of the species presented high marginality and lower tolerance values, which indicates a smaller niche breadth, and they were therefore categorized as specialists. Labidomera suturella Guérin-Méneville, 1838 was the species with the highest marginality and the lowest tolerance, followed by Walterianella sp. 1, Zenocolaspis inconstans (Lefèvre, 1878) and Alagoasa trifasciata (Fabricius, 1801) ( Table 4). The aforementioned species had lower abundance, 1-15 specimens, in a minor number of samples (Fig. 2).
In the case of the dry season, marginality was significant (Monte Carlo test, p = 0.017) for only seven of the 40 registered species. Two were considered as generalists, with low marginality values; of these, B. pumila presented the highest tolerance, while  C. diversa showed the lowest marginality (Table 4). Abundance of both species was uniformly distributed in almost all samples (Fig. 3). The other five chrysomelids had high marginality and low tolerance values (specialists): the highest marginality and lowest tolerance occurred in A. trifasciata, and it was consequently the species most specialized to microclimatic conditions during the dry season in the semideciduous tropical forest. In descending order, Syphrea sp. 1, Chaetocnema sp. 1, Acallepitrix sp. 7, and Epitrix sp. 1 (Table 4) were species recorded in few samples, with abundances between four and 18 specimens (Fig. 3). Heat index, evapotranspiration and temperature were the microclimatic variables most related with the abundance of leaf beetle species during the rainy season and were represented in Axis 1 of the OMI analysis (Eigenvalue = 4.9077, inertia = 55.74%). In Axis 2 (Eigenvalue = 2.6344, inertia = 29.92%) the most important variable was the average wind speed (Table 5). For the dry season, evapotranspiration, temperature, and heat index in Axis 1 (Eigenvalue = 7.9982, inertia = 75.67%) were the microclimatic variables most associated with the changes in abundance of leaf beetles. Maximum wind speed had the highest correlation in Axis 2 (Eigenvalue = 1.7084, inertia = 0.1616%) ( Table 5).
The association of the species with the environmental variables was determined based on the positions of the centroids and their closeness with respect to Axes 1 and 2. Those species that were located very close to the origin of both axes were considered to be related to average microclimatic values. For the rainy season, A. trifasciata and Z. inconstans, were related with low values of average wind speed ( (18.10-20.57 °C), but higher wind speed (1.06-2.12 m/s). Lastly, C. diversa and M. bumeliae did not follow a specific pattern in relation to the significant variables in any axis since they were at the origin of the niche dispersion (Fig. 4).
During the dry season, the average distribution of Syphrea sp. 1, Acallepitrix sp. 7, and Epitrix sp. 1 was correlated with areas of lower evapotranspiration (13-16.82 °C), temperature (16.30-19.69 °C) and heat index (16.60-22.20 °C) in Axis 1. Similarly on Axis 2, these species predominated under conditions of low to average maximum wind speed (1.42-2.84 m/s). Chaetocnema sp. 1 occurred in conditions of low evapotranspiration (13-14.91 ° C) and low temperature (21.39-23.09 °C), as well as low heat index (19.40-22.20 °C), but this species was associated with high values of maximum wind speed (1.42-2.13 m/s). Alagoasa trifasciata was the species with the lowest tolerance value; so, its centroid was positioned in areas with high evapotranspiration values (22.56-24.47 °C), high temperature (24.79-26.49 °C), high heat index (30.60-33.40 °C), and low maximum wind speed (0-0.71 m/s). Finally, the centroid of the distribution of B. pumila and C. diversa was significantly associated with average microclimatic conditions, since their distribution included areas with high and low values for the heat index, as well as for the other variables (Fig. 5).

Faunistic inventory and chrysomelid biodiversity
Prior to this study, 2,660 species of Chrysomelidae had been recorded from Mexico (Niño-Maldonado and Sánchez-Reyes 2017) and 257 from the state of Tamaulipas ). Accordingly, results in the STF of the study area represent 2.7% of the leaf beetle biodiversity reported for the country, and 27.6% for the state. Our study revealed Diachus chlorizans (Suffrian, 1852) as a new country record for Mexico, and Diabrotica biannularis Harold, 1875 as a new state record for Tamaulipas. These records were previously published in preliminary works from the study area (  research is above other STF fragments in northeastern Mexico, and it represents 68.2% of the most biodiverse site. Regarding true diversity, the numbers of equally dominant (1/D) and typical (e H ) species in this study were lower than those observed in Peregrina Canyon ), although they were higher than those observed in STF fragments from Cañón del Novillo (Martínez-Sánchez 2016) or Cerro El Diente (Sánchez-Reyes et al. 2015b).
Galerucine dominance as observed in our study has also been reported in other studies in northeastern Mexico (Niño-Maldonado et al. 2005;Furth 2009;Furth 2013;Sánchez-Reyes et al. 2013;Sánchez-Reyes et al. 2014;Bouzan et al. 2015;Flinte et al. 2017;Lucio-García et al. 2019, 2020, and this may be due to the subfamily's high number of species (Riley et al. 2002), with specimens found in all ecosystems during most parts of the year (Furth 2013;Sánchez-Reyes et al. 2013). In contrast, the subfamily composition of this study is quite different from that in tropical forests. In the Chamela region, on the Pacific side of Mexico, 49 species of Cassidinae were listed (Noguera 1988).
As a whole, the aforementioned results highlight the great importance of the study area, since it was possible to find a large percentage of species within a smaller expanse when compared to larger space-temporal gradients or natural protected areas. This can be attributed to the geographic location of the studied STF within a region with a high conservation priority (CONABIO et al. 2007). The area, although adjacent to the Altas Cumbres Natural Protected Area, constitutes a mosaic with fragments of different durations since last disturbance, and this may favor the presence of a complex community of species . Furthermore, the STF is one of the ecosystems with the highest biodiversity of plants (Rzedowski 2006;García-Morales et al. 2014), and it is one of the most important in terms of chrysomelid species richness in Mexico (Noguera 1988;Burgos-Solorio and Anaya-Rosales 2004;Niño-Maldonado et al. 2005;Lucio-García et al. 2019). The combination of environmental factors in the STF results in a great diversification of plants, providing a wide range of food resources, which could lead to the high number of leaf beetles in this plant community in Tamaulipas and other states of Mexico.

Seasonal variation
On a temporal scale, the chrysomelid community followed a seasonal pattern, where the rainy season was the most favorable for the presence of this group in the study area. Increase in abundance and species richness during this season has also been found in numerous studies worldwide, including studies in Tamaulipas and other parts of Mexico (Petitpierre et al. 2000;Esker et al. 2002;Burgos-Solorio and Anaya-Rosales 2004;Koji and Nakamura 2006;Furth 2009;Martínez-Sánchez et al. 2009;Furth 2013;Sánchez-Reyes et al. 2015b;Sánchez-Reyes et al. 2016a;Sandoval-Becerra et al. 2016;Şen and Gök 2016;Miwa and Meinke 2017;Lucio-García et al. 2019). Results of the richness estimators support these patterns, because the percentage of completeness during rains is lower when compared to the dry season. Thus, in certain areas, the highest activity of chrysomelids is restricted to the rainy season, while inactivity increases during drought conditions (Noguera 1988;Furth 2013). This is due to the association of chrysomelids with the quality and availability of their host plants (Řehounek 2002;Şen and Gök 2016), which are some of the most important elements in their diet (Ávila and Postali-Parra 2003), as well as with the abundance of young foliage (Basset and Samuelson 1996), variables that are increased during the period of highest rainfall. In addition, there is more vegetation cover producing shade, creating microenvironments that could be more favorable to maintaining a high population density (Hill and Hill 2001), with the climatic conditions of humidity necessary for the adult beetles to emerge and fly (Yanes-Gómez and Morón 2010).
However, in other geographic regions, such as the subtropical areas of Brazil, the highest abundance has occurred in the dry season, specifically within the subfamilies Galerucinae, Cassidinae and Chrysomelinae (Linzmeier and Ribeiro-Costa 2008;Flinte et al. 2011;Bouzan et al. 2015;Flinte et al. 2017). In addition, in some areas of northeastern Mexico, greater numbers of species and specimens have also been recorded during the dry season ). These discrepancies can be attributed to the climatic and biogeographic differences between plant communities. For example, in cloud forests, dry periods are shorter and less intense, causing a favorable increase in specimens of some Coleoptera families (Pedraza et al. 2010). Although soil moisture and precipitation are reduced in these areas, the cloudiness in the form of mist reduces evaporation, providing water during periods of low rain; in consequence, marked deficiency of humidity in these forests is rare. In other tropical forests near the study area, the dry season is not as severe, for example in the Peregrina Canyon, where a higher abundance of adult chrysomelids often occurs concentrated in refuges during this season, while the larval stages are more abundant during the rains ). On the contrary, differences in geographic position, latitude and elevation influence the contrast that exists between the dry and rainy seasons in other fragments of the same type of vegetation in northeastern Mexico. In the study area, there are well-defined periods of high temperature and precipitation, in addition to a non-continuous flow of water currents during the year, which lead to a more severe dry season. Similar and more extreme cases exist in tropical dry forests from the south or Pacific coast of Mexico, where the plants lose their leaves completely during the dry season, resulting in a notable absence of chrysomelids (Noguera 1988). Likewise, these climatic variations and their effects on the phenology of the host plants are probably the main drivers of the temporal dynamics in these beetles (Flinte et al. 2017).
Unlike other investigations where the greatest diversity also occurs in the wet season (Sánchez-Reyes et al. 2016a), in this work, the low abundance and species richness resulted in a high diversity in the dry season, by decreasing the dominance and increasing the effective number of species (Magurran 2004). Therefore, the dry season is of great importance for the chrysomelid community in the STF of the study area, since the prevailing conditions increase the evenness of the chrysomelid community. Species may exploit food resources in a more efficient way during this season, achieving a balance in their populations and reducing the dominance of most species, thus suggesting an adaptation of Chrysomelidae to acute drought conditions. This could be noted also when observing the high percentage of faunistic similarity between seasons, which indicates that most of the leaf beetles are the same in both periods. Therefore, it is possible that their resource acquirement changes and consequently their abundances are modified during the seasonal variations. Moreover, 31 species were registered exclusively for the rainy season, while only 11 for the dry season. Together, these results highlight the relevance of areas where there is a marked temporal or seasonal heterogeneity, since it can generate unique species compositions.

Response of Chrysomelidae to seasonal microclimatic changes
In this research, the niches of chrysomelid species were examined by means of the Outlying Mean Index. This showed that the variations in the abundance of leaf beetles were significantly related to the microclimatic changes in each season. Factors that influence the distribution of phytophagous insects are a combination of geographic and environmental elements (Wąsowska 2004;Andrew and Hughes 2005;Lassau et al. 2005;Baselga and Jiménez-Valverde 2007). It has also been shown that leaf beetles present different degrees of association with the microclimatic conditions of the habitats where they develop (Sánchez-Reyes et al. 2016b;Sandoval-Becerra et al. 2017), and this is demonstrated in our study. However, the variation explained by the analysis and the correlation values of the variables were higher in the dry season, suggesting a stronger association between the microclimate and the chrysomelid community with respect to the rainy season. This can be attributed to more heterogeneous environment values during low precipitation months. For example, in tropical forests it has been observed that lower microclimatic variability occurs through the rainy season (Checa et al. 2014;Sánchez-Reyes et al. 2019), which could be due to a higher homogeneity in the vegetation structure. Therefore, chrysomelid populations are more variable in relation to the seasonal microclimate prevailing during the dry season, so that the effects, particularly of precipitation, determine strong positive or negative responses in these insects (Pinheiro et al. 2002); this pattern also occurs in other phytophagous groups, such as Curculionidae or Cicadidae (Novotny et al. 1999;Silva et al. 2017).
Significant microclimatic variables were very similar between seasons (environmental temperature, heat index, evapotranspiration and, to a lesser extent, wind speed), although there were differences in the order of importance and in their contribution to the variations in abundance of leaf beetles. In the rainy season, the most important variable to characterize the niche of the species was the heat index, which is considered to be a combination of humidity and temperature in the same value and represents the thermal sensation (Lee and Brenner 2015). In physiological terms, phytophagous insects must accumulate a certain amount of heat to be able to hatch and accelerate their development rate, thereby increasing the number of generations (Marco 2001;Mejía 2005). However, in the dry season, the variable of greatest importance was evapotranspiration. Such variation can be attributed to the environmental humidity stress to which the host plants are exposed after rainfall, modifying the moisture content of leaves and stems and thereby affecting feeding patterns of chrysomelids. This suggests that direct effects on trophic networks may occur during drought periods, which influence the development of phytophagous insects, particularly due to desiccation (Martínez et al. 2010;Anderson et al. 2016).
A similar set of microclimatic variables has been associated with Chrysomelidae in other works, specifically temperature, heat index, maximum wind speed and evapotranspiration (Stewart et al. 1996;Flinte and Valverde de Macedo 2004;Isard et al. 2004;Baselga and Jiménez-Valverde 2007;Linzmeier and Ribeiro-Costa 2013;Aneni et al. 2014;Sánchez-Reyes et al. 2016b;Oliveira et al. 2017;Sandoval-Becerra et al. 2017). There are other studies where the most important abiotic variables were solar radiation, precipitation, relative humidity, photoperiod and condensation point (Flinte and Valverde de Macedo 2004;Isard et al. 2004;Linzmeier and Ribeiro-Costa 2013). It should be mentioned that differences compared to the present study arise due to various factors, including the type of study, geographic location, ecosystems evaluated and method for measuring microclimatic variables, as well as the specific response of taxa to the variables.
Regarding the individual response of leaf beetles to the variables, it was observed that only 11 of the 71 species registered a significant variation between their niche and the average microclimatic conditions in the STF. It has been observed that the number of chrysomelid species that present a significant relationship with abiotic parameters is variable, although previous studies have focused on the effect of disturbance (Sandoval-Becerra et al. 2017) and elevation gradients (Sánchez-Reyes et al. 2016b).
The study of chrysomelid species associated with abiotic variables has been useful in recognizing part of their biology and ecology, specifically their reproductive cycle or their potential for biological control (Stewart et al. 1996;Flinte and Valverde de Macedo 2004;Isard et al. 2004;Oliveira et al. 2017). Also, such study has been applied to know their niches (Sánchez-Reyes et al. 2015a) and to identify indicator species of conservation status (Ohsawa and Nagaike 2006) or disturbance (Sandoval-Becerra et al. 2017). Other studies have focused on the influence of climate in the distribution of species (Sánchez-Reyes et al. 2016b;Wang et al. 2017). The present work, on the other hand, is one of the first to address the influence of microclimatic variation on Chrysomelidae from a seasonal perspective.
Specifically, in the dry season, seven significant species were recorded, while in the rainy season there were only six. In both seasons, C. diversa was categorized as a generalist species, since it presented a low marginality and a high tolerance, which indicates a wide distribution in the study area associated with average microclimatic values. This response is similar to that observed in the same and other species within the genus, but in different areas (Sánchez-Reyes et al. 2016b;Sandoval-Becerra et al. 2017). The second common species in both seasons was A. trifasciata although it was categorized as a specialist due to high marginality and low tolerance values. A similar response pattern was previously recorded in the Sierra de San Carlos for this species (Sánchez-Reyes 2014). Association between variables and A. trifasciata was higher in the dry season, since the marginality parameters were higher, while the tolerance values were lower; that is, the distribution of A. trifasciata appears to be more restricted during the dry conditions. The rest of the species demonstrated seasonal differences. For example, during the rainy season, L. suturella presented the highest marginality value and had a low tolerance to the microclimatic environment; similar responses were observed in Walterianella sp. 1 and Z. inconstans. In the dry season, Syphrea sp. 1 was the species with the lowest tolerance, followed by Chaetocnema sp. 1, Acallepitrix sp. 7, Epitrix sp. 1 (Galerucinae), and B. pumila (Cassidinae). However, some of these species also occurred throughout the year, despite being significantly associated with only one season. The above observations provide evidence that leaf beetles have seasonal modifications in their niche requirements. Influence of the microclimate may be more important in the rainy season, while in the dry season (or vice versa) the variables that determine niches are different, or they may have a non-significant contribution to the distribution of the species (Basset et al. 1992;Martínez et al. 2010;García-Atencia et al. 2015). These seasonal changes may be associated with the synchronization of the reproductive cycles of the phytophagous insects, particularly depending on the precipitation and temperature provided by the forest structure, which is not constant throughout the year and tends to be increasingly variable (Basset et al. 1992;Soler et al. 2002;García-Atencia et al. 2015).
The broad microclimatic tolerance of C. diversa and the abiotic specialization of A. trifasciata represent a first approach to the analysis of the generalized environmental response of chrysomelids, even though both have been documented in other studies. In this way, it is probable that the behavior of the species is similar and constant in other geographical areas, which would allow the use of such taxa in environmental monitoring. New studies on chrysomelid niches would allow us to elucidate these effects. It is also important to recognize that phytophagous insects and specialist taxa with a small niche breadth could be negatively influenced by the possible effects of climate change (Williams et al. 2007;Dormann et al. 2008;Hill et al. 2011), which will impact the structure and functioning of the communities (Hegland et al. 2009;Stuble et al. 2013;Luna-Castellanos et al. 2017). Effects extend to plant-insect interactions (mutualism, predation, competition, etc.), either due to phenological changes (synchronization in the interaction) or distribution of species (Hódar et al. 2004;Luna-Castellanos et al. 2017), with some species even being susceptible to local extinction (Tscharntke et al. 2002;Petermann et al. 2010). Furthermore, the present results and similar evidence suggest that climate variability can lead to significant biodiversity losses (Parmesan et al. 1999;Hill et al. 2002;Konvicka et al. 2003;Wilson et al. 2007). However, despite having knowledge about possible consequences, little information is available on the effects that the changing microclimate can have on biodiversity, its populations, biological communities, and the ecosystems that harbor them.

Conclusions
The study of seasonal and microclimatic changes on species and communities is a topic of great importance in conservation ecology. Community attributes of the family Chrysomelidae and the beetles' response to microclimatic variation were evaluated for the first time from a seasonal perspective, in a semideciduous tropical forest fragment of northeastern Mexico. Overall, the observed results were similar to those from other faunistic studies of leaf beetles, although the number of species ranked third within tropical forest areas of the state of Tamaulipas. Seasonality induced significant changes in the parameters of abundance, diversity and faunistic composition in the chrysomelid community. The highest number of specimens and species were recorded in the rainy season, while the lowest dominance and highest diversity occurred in the driest period.
In this study, it was shown that Chrysomelidae were significantly associated with the microclimatic variation among seasons. However, the strength of this association and the number of significant species were different for each season. Changes in the abundance of the leaf beetles were influenced by the heat index, temperature, evapotranspiration, and average wind speed, reflected by specific conditions required for each species. Microclimatic and seasonal assessment could be useful for the evaluation of climate change, since niche analysis enables detection of specialized or vulnerable species, which are associated with a delimited set of environmental conditions. This characterization of the microclimate niche of Chrysomelidae from a seasonal perspective was conducted here for the first time in northeastern Mexico. However, additional studies are warranted to determine if the observed patterns are different when evaluating other abiotic factors or when evaluating other plant communities.