Acessibilidade / Reportar erro

High turnover of Chrysomelidae (Coleoptera) species in semideciduous forest remnants in an agricultural landscape

Abstract

Differences in species composition between sites (β diversity) may be the result of spatial species replacement (turnover) or nestedness (subgroups of species from a more diverse site). In fragmented landscapes, the environmental factors that lead to these differences may be spatially structured. Herein, our objective is to determine if the β diversity of Chrysomelidae (Coleoptera) is due to turnover or nestedness and whether the observed pattern is due to loss of forest cover or spatial processes in forest remnants immersed in a matrix dominated by intense agricultural practice. We used an incidence matrix of 99 species sampled from 16 forest remnants and found that the difference in species composition among the fragments is mostly determined by turnover and that this variation is not explained by forest cover or spatial variables. In regions where high habitat loss has generated landscapes containing small and islated forest fragments, structural features, related both to habitat (area, isolation, shape, etc.) and landscape (land use, landscape heterogeneity, etc.) could predict diversity patterns.

Key words
beetles; beta diversity; fragmentation; herbivory; insects

INTRODUCTION

The land-use conversion is one of the major threats to biodiversity, owing to its association with habitat loss and fragmentation. Currently, about 70% of the tropical forest remnants are isolated from original forests and subject to deleterious effects of fragmentation (Haddad et al. 2015HADDAD NM ET AL. 2015. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci Adv 1: 1-9.). In this anthropogenic process, the remaining native vegetation cover is generally fragmented into many small patches (Fahrig et al. 2019FAHRIG L ET AL. 2019. Is habitat fragmentation bad for biodiversity? Biol Conserv 230: 179-186.) and, as a result, many animal and plant populations are isolated in small, disconnected fragments surrounded by a matrix that is usually composed of agriculture and pasture, unsuitable for the survival of many native species (Fernández-Chacón et al. 2014FERNÁNDEZ-CHACÓN A, STEFANESCU C, GENOVART M, NICHOLS JD, HINES JE, PÁRAMO F, TURCO M & ORO D. 2014. Determinants of extinction-colonization dynamics in Mediterranean butterflies: The role of landscape, climate and local habitat features. J Anim Ecol 83: 276-285.). Isolation and size of fragments can lead to changes in the composition and structure of animal and plant communities over time, which may generate local extinctions and consequently affect essential ecosystem functions, such as decomposition (Didham et al. 1996DIDHAM RK, GHAZOUL J, STORK NE & DAVIS AJ. 1996. Insects in fragmented forests: a funcional approach. Tree 11: 225-260.). Small fragments are often considered of low conservation priority because they do not support viable populations and are very vulnerable to local extinctions. However, studies show that small fragments can increase beta diversity by creating irregularities in spatial distribution of species and different local extinction dynamics between fragments (Haddad et al. 2015HADDAD NM ET AL. 2015. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci Adv 1: 1-9.). In addition, a set of small habitat patches collectively may harbor more species than certain large patches with similar total area (Fahrig et al. 2019FAHRIG L ET AL. 2019. Is habitat fragmentation bad for biodiversity? Biol Conserv 230: 179-186.), resulting in high β diversity in fragmented landscapes (Tscharntke et al. 2012TSCHARNTKE T ET AL. 2012. Landscape moderation of biodiversity patterns and processes - eight hypotheses. Biol Rev 87: 661-685.).

Differences in species composition between assemblies, referred to as beta diversity (β) by Whittaker (1960)WHITTAKER RH. 1960. Vegetation of the Sisiyou Mountains, Orgeon and California. Ecol Monogr 30: 279-338, can be generated by two distinct phenomena: species spatial replacement (turnover) or nestedness (Harrison et al. 1992HARRISON S, ROSS SJ AND LAWTON JH. 1992. Beta Diversity on Geographic Gradients in Britain. J Anim Ecol 61: 151-158., Williams 1996WILLIAMS PH. 1996. Mapping variations in the strength and breadth of biogeographic transition zones using species turnover. Proc R Soc B Biol Sci 263: 579-588., Lennon et al. 2001LENNON JJ, KOLEFF P, GREENWOOD JJD & GASTON KJ. 2001. The geographical structure of British bird distributions : diversity , spatial turnover and scale. J Anim Ecol 70: 966-979.). Turnover is described as spatial substitutions of some species by other species along an ecological or temporal gradient, which implies simultaneous gain and loss of species due to environmental filters, competition and historical events (Leprieur et al. 2011LEPRIEUR F, TEDESCO PA, HUGUENY B, BEAUCHARD O, DÜRR HH, BROSSE S & OBERDORFF T. 2011. Partitioning global patterns of freshwater fish beta diversity reveals contrasting signatures of past climate changes. Ecol Lett 14: 325-334.). On the other hand, nestedness describes a condition in which sites with fewer species are subgroups of a site with a greater number of species (Baselga 2010BASELGA A. 2010. Partitioning the turnover and nestedness components of beta diversity. Glob Ecol Biogeogr 19: 134-143.). Nestedness, can reflect the diversity of niches present in the region or other ecological processes such as physical barriers (Legendre 2014LEGENDRE P. 2014. Interpreting the replacement and richness difference components of beta diversity. Glob Ecol Biogeogr 23: 1324-1334.). In this way, partitioning the β-diversity among turnover and nestedness is important for the development of biodiversity conservation strategies, as well as for accurate understanding of processes driving the observed patterns (Felinks et al. 2011FELINKS B, PARDINI R, DIXO M, FOLLNER K, METZGER JP & HENLE K. 2011. Effects of species turnover on reserve site selection in a fragmented landscape. Biodivers Conserv 20: 1057-1072.).

In the present study, our objective is to determine if turnover, nestedness or randomness would drive β-diversity of leaf beetle species (Chrysomelidae: Coleoptera) of forest remnants immersed in an agricultural matrix. We also aim to understand the contribution of spatial structures and forest cover on the variation of total beta diversity and its components turnover and nestedness.

Chrysomelidae is a family of beetles with wide geographical distribution and enormous diversity, constituted by small beetles that represent a large part of the herbivorous insect fauna in diverse biomes (Andrew & Hughes 2004ANDREW NR & HUGHES L. 2004. Species diversity and structure of phytophagous beetle assemblages along a latitudinal gradient: predicting the potential impacts of climate change. Ecol Entomol 29: 527-542.). More than 36,000 species are currently recognized in the family (Bouchard et al. 2009BOUCHARD P, GREBENNIKOV VV, SMITH ABT & DOUGLAS H. 2009. Diversity of Coleoptera. In: Foottit R & Adler P (Eds), Insect Biodiversity: Science and Society. Blackwell Publishing, p. 265-301.), with estimates of more than 60,000 species (Jolivet 2015JOLIVET P. 2015. Together with 30 years of Symposia on Chrysomelidae! Memories and personal reflections on what we know more about leaf beetles. Zookeys 547: 35-61.). In Brazil, more than 6,000 species are registered from more than 550 genera (Sekerka et al. 2017SEKERKA L, LINZMEIER AM, MOURA LA, RIBEIRO-COSTA CS, AGRAIN F, CHAMORRO ML, MANFIO D, MORSE GE & REGALIN R. 2017. Chrysomelidae. In: Catálogo Taxonômico da Fauna do Brasil. http://fauna.jbrj.gov.br/fauna/faunadobrasil/115540.
http://fauna.jbrj.gov.br/fauna/faunadobr...
), however, in many states there is still a lack of knowledge about the diversity and ecological patterns of these insects, as in Mato Grosso do Sul. These insects are sensitive to environmental disturbances, respond to differences in habitat structural complexity (Linzmeier & Ribeiro-Costa 2009LINZMEIER AM & RIBEIRO-COSTA CS. 2009. Spatio-temporal dynamics of Alticini (Coleoptera, Chrysomelidae) in a fragment of Araucaria Forest in the state of Parana, Brazil. Rev Bras Entomol 53(2): 294-299., Sandoval-Becerra et al. 2018SANDOVAL-BECERRA FM, SÁNCHEZ-REYES UJ, CLARK SM, VENEGAS-BARRERA CS, HORTA-VEGA JV & NIÑO-MALDONADO S. 2018. Influence of Habitat Heterogeneity on Structure and Composition of a Chrysomelidae (Coleoptera) Assemblage in a Temperate Forest in Northeast Mexico. Southwest Entomol 43(1): 115-130., Teles et al. 2019TELES TS, RIBEIRO DB, RAIZER J & LINZMEIER AM. 2019. Richness of Chrysomelidae (Coleoptera) depends on the area and habitat structure in semideciduous forest remnants. Iheringia Ser Zool 109: 1-8.) and have potential as bioindicators (Pimenta & De Marco 2015PIMENTA M & DE MARCO PJ. 2015. Leaf Beetle (Chrysomelidae: Coleoptera) Assemblages in a Mosaic of Natural and Altered Areas in the Brazilian Cerrado. Neotrop Entomol 44: 242-255., Sánchez-Reyes et al. 2019SÁNCHEZ-REYES UJ, NIÑO-MALDONADO S, CLARK SM, BARRIENTOS-LOZANO L & ALMAGUER-SIERRA P. 2019. Successional and seasonal changes of leaf beetles and their indicator value in a fragmented low thorn forest of northeastern Mexico (Coleoptera, Chrysomelidae). Zookeys 825: 71-103.).

MATERIALS AND METHODS

The study was carried out in the Dourados municipality, located in the central-southern portion of the State of Mato Grosso do Sul, Brazil (22 ° 13’15 “S, 54 ° 48’21” W, 430 m altitude) (Fig. 1). The region is an ecotone with natural vegetation composed mostly of Alluvial Seasonal Forest and Cerrado, with influences of Atlantic Forest and Meridional Forest (Rizzini 1997RIZZINI CT. 1997. Tratado de Fitogeografia do Brasil: aspectos ecológicos, sociológicos e florísticos, 2nd ed., Rio de Janeiro: Âmbito Cultural Edições Ltda.), resulting in very diverse vegetation. However, the process of territorial occupation of the region was characterized by a lack of planning, causing extensive destruction of natural resources, leading to the replacement of native vegetation by agricultural crops and pastures for cattle raising (Martins 2001MARTINS SV. 2001. Recuperação de matas ciliares. Aprenda Fácil: Viçosa, 219 p.). The climate of the region, according to Köpen classification, is Cwa type, humid mesothermic, with wet summers and dry winters (Fietz & Fisch 2008FIETZ CR & FISCH GF. 2008. O Clima da Região de Dourados, MS, 2nd ed., Embrapa Agropecuária Oeste, Documentos, 92, 32 p.) and a predominance of very clay-like Latosol (Amaral et al. 2000AMARAL JAM DO, MOTCHI EP, OLIVEIRA H DE, CARVALHO FILHO A, NAIME UJ & DOS SANTOS RD. 2000. Levantamento semidetalhado dos solos do campo experimental de Dourados, da Embrapa Agropecuária Oeste, município de Dourados, MS. Embrapa Agropecuária Oeste, Dourados.).

Figure 1
Collection sites (16 forest remnants) of Chrysomelidae in Dourados, Mato Grosso do Sul.

Between August 2012 and March 2013, we sampled 16 forest remnants with varying sizes (0.61-308ha) (Table I using Malaise traps (Townes 1972TOWNES H. 1972. A light-weight malaise trap. Entomol News 83: 239-247.). The sampling effort between the fragments was not the same, as it depended on the availability and access to privately owned lands/properties. In total, we collected 60 samples, with each sample equivalent to 14 consecutive days of trap exposure. We separated the leaf beetles in morphotypes and identified through comparisons with specimens deposited at the “Coleção Entomológica Pe. Jesus Santiago Moure”, Universidade Federal do Paraná. We consider morphotypes those individuals that were not possible to identify at species level, but hereafter we name them species. The specimens vouchers are deposited in the “Museu da Biodiversidade da Faculdade de Ciências Biológicas e Ambientais”, Universidade Federal da Grande Dourados.

Table I
Forest fragments of Chrysomelidae (Coleoptera) sampling in Dourados, Mato Grosso do Sul State, Brazil, with Malaise traps. S = species richness; N = abundance; Buffer (%) = percentage of forest cover in buffers of 250m and 1500m radius from the centroid of each fragment.

All analyzes were performed in the R program (R Core Team 2019R CORE TEAM. 2019. R: A language and environment for statistical computing. https://www.r-project.org/.
https://www.r-project.org/...
) with different statistical packages (see details below). In order to determine if the differences in species composition found in fragments were due to turnover (spatial replacement) or nestedness, we used a presence and absence matrix containing the 99 leaf beetle species. With this matrix, we calculated the total diversity, as well as its turnover and nestedness components, using the “beta.multi” function of the betapart R package (Baselga et al. 2018BASELGA A, ORME D, VILLEGER S, DE BORTOLI J & LEPRIEUR F. 2018. betapart: Partitioning Beta Diversity into Turnover and Nestdness Components. R package version: 1.5.1. https://CRAN.R-project.org/package=betapart.
https://CRAN.R-project.org/package=betap...
). We used the Sørensen dissimilarity index (βSOR), which was partitioned into the turnover components, calculated by the Simpson dissimilarity index (βSIM), and nestedness, expressed by the difference of the Sørensen dissimilarity index minus the dissimilarity index of Simpson (βSNE) (Baselga 2010BASELGA A. 2010. Partitioning the turnover and nestedness components of beta diversity. Glob Ecol Biogeogr 19: 134-143.). We tested the null hypothesis that leaf beetle diversity is randomly spatially distributed between the fragments by comparing the observed data against the results obtained by null models with the “oecosimu” function of the vegan package (Oksanen et al. 2018OKSANEN J ET AL. 2018. vegan: Community Ecology Package. R package version: 2.5-1. https: //cran.r-project.org/package=vegan.
https: //cran.r-project.org/package=vega...
) using the model construction method “r1” with 1000 permutations.

We used the percentage of forest cover in 250 and 1500 m radius buffers for each of the fragments as environmental variables. These variables were obtained through Landsat 8 images via the USGS EarthExplorer database (http://earthexplorer.usgs.gov/) in the ArcGis Desktop 10.5 program. The spatial variables were generated using the distance-based Moran Eigenvector Maps (dbMEM) method through the “dbmem” function of the adespatial package (Dray et al. 2019DRAY S ET AL. 2019. adespatial: Multivariate Multiscale Spatial Analysis. R package version: 0.3-4. https://CRAN.R-project.org/package=adespatial.
https://CRAN.R-project.org/package=adesp...
). This analysis creates spatial variables (axes) with the geographic coordinates of the centroids of the fragments, that represent spatial structures at different scales. We selected the first four axes that presented positive spatial correlation, so that the first two axes generated by dbMEM represent spatial variables referring to a broad scale, while the latter two represent a fine scale.

In order to determine the contribution of spatial variables and forest cover (environmental component) in the variation of beta diversity and its components, we first calculated three dissimilarities matrices that represent the variation of the total beta diversity, the turnover component and the nestedness component. We calculated these matrices using the “beta.pair” function of the betapart package (Baselga et al. 2018BASELGA A, ORME D, VILLEGER S, DE BORTOLI J & LEPRIEUR F. 2018. betapart: Partitioning Beta Diversity into Turnover and Nestdness Components. R package version: 1.5.1. https://CRAN.R-project.org/package=betapart.
https://CRAN.R-project.org/package=betap...
). We performed a distance-based redundancy analysis (db-RDA) with variation partitioning for each response matrix (Dray et al. 2012DRAY S ET AL. 2012. Community ecology in the age of multivariate spatial analysis. Ecol Monogr 82: 257-275., Legendre et al. 2012LEGENDRE P, BORCARD D & ROBERTS DW. 2012. Variation partitioning involving orthogonal spatial eigenfunction submodels. Ecology 93: 1234-1240.) to understand the importance of spatial and environmental variables in the variation of these matrices. This analysis divides the variation in the response matrix into four components: [a] pure environmental component; [b] component shared between environment and space; [c] pure spatial component; and [d] unexplained variation.

RESULTS

We found 450 leaf beetles belonging to 99 species (Table II. The most abundant and speciose subfamily was Galerucinae with 273 individuals and 43 species. Fifty-one species (51.5% of total) presented only one individual and 10 species had 10 or more individual. The most abundant species (≥ 10 individual) were Galerucinae sp.1 (53 individuals), followed by Wanderbiltiana sejuncta (Harold, 1880) (36 individuals), Wanderbiltiana sp.1 (35 individuals), Galerucinae sp.4 (32 individuals), Systena sp.1 (32 individuals) and Costalimaita ferruginea (Fabricius, 1801) (24 individuals) (Table II). Only four species were distributed in more than five fragments. Wanderbiltiana sp.1 and W. sejuncta were the most frequent species, found in nine of 16 fragments, and Maecolaspis sp.1 and Systena sp.1 found in five fragments.

Table II
Frequency and abundance (number of leaf beetles) of most frequent (> 5 fragments) and abundant (≥ 10 leaf-beetles) Chrysomelidae (morpho)species* sampled with Malaise traps in 16 remnants of semi-deciduous forest (sites) in Dourados, Mato Grosso do Sul, Brazil.

The beta diversity found in fragments was high (βSOR = 0.94) and due almost exclusively to turnover (βSIM = 0.90), compared to nestedness (βSNE = 0.04). This pattern differs significantly from what would be expected by randomness, indicating that leaf beetle species are not randomly distributed among the remaining forest fragments (Table III.

Table III
Values of beta diversity (β SOR) with turnover (β SIM) and nestedness (β NES), calculated according to Baselga (2010), of Chrysomelidae (Coleoptera) from forest remnants of Dourados, Mato Grosso do Sul State, Brazil. Obs = observed value; SES = standardized effect size; Est = average of the estimated values based on 1000 simulations of null models (“r1” method); 2.5%, 50% and 97.5% = percentiles of the distribution of the estimated values; P = probability based on the comparisons between the observed values and the mean of the estimated values.

The variation partitioning showed that the variation of general beta diversity (βSOR) is explained almost exclusively by forest cover, however this explanation is low (2%) and insignificant (df = 2, F = 1.1872, P = 0.158; Fig. 2a). In the same way, forest cover explained little (3.8%) and was insignificant (df = 2, F = 1.2315, P = 0.214) for the variation in spatial replacement (turnover component - βSIM; Fig. 2b). Regarding the variation of the nestedness component (βSNE), forest cover and spatial component have no influence. On the other hand, unlike βSOR and βSIM, the shared fraction between forest cover and spatial component accounts for about 21% of the variation in nestedness (Fig. 2c). The variation explained by the spatial variables in βSOR, βSIM and βSNE, as well as the shared fraction between forest cover and spatial variables in βSIM presented negative values ​​and are interpreted as zero (Legendre 2008LEGENDRE P. 2008. Studying beta diversity: ecological variation partitioning by multiple regression and canonical analysis. J Plant Ecol 1: 3-8.), as they indicate that the explanatory variables explain less variation than the normal random variables.

Figure 2
Venn diagrams showing the proportions of variation in general beta diversity (βSOR; a), turnover (βSIM; b) and nestedness (βNES; c) in relation to forest cover (VegCov) and spatial variables (Space). Res = variation not explained by factors; “-” = variation explained ≤ 0.

DISCUSSION

Our results show that in forest remnants immersed in an agricultural matrix, the beta diversity of Chrysomelidae is high and differs significantly from what would be expected by randomness. This fact indicates that species are not randomly distributed among the fragments. We also observed that the difference in species composition between fragments is almost exclusively due to the turnover process, with great spatial substitution of species between the areas and a low contribution of the nestedness to this difference. However, although we see that species are not randomly distributed among fragments, this variation is neither caused by vegetation cover (environmental variables) nor by spatial variables.

The pattern of spatial species replacement (turnover) has been reported for several insect groups (e.g. Tenebrionidae, Fattorini & Baselga 2012FATTORINI S & BASELGA A. 2012. Species richness and turnover patterns in european tenebrionid beetles. Insect Conserv Divers 5: 331-345., Carabidae, Gañán et al. 2008GAÑÁN I, BASELGA A & NOVOA F. 2008. Diversity patterns in Iberian Calathus (Coleoptera, Carabidae: Harpalinae): species turnover shows a story overlooked by species richness. Environ Entomol 37: 1488-1497., Zygoptera, Brasil et al. 2018BRASIL LS, OLIVEIRA-JÚNIOR JM, CALVÃO LB, CARVALHO FG, MONTEIRO-JÚNIOR CS, DIAS-SILVA K & JUEN L. 2018. Spatial, biogeographic and environmental predictors of diversity of Amazonian Zygoptera. Insect Conserv Diver 11(2): 174-184.) including leaf beetles (Freijeiro & Baselga 2016FREIJEIRO A & BASELGA A. 2016. Spatial and environmental correlates of species richness and turnover patterns in European cryptocephaline and chrysomeline beetles. Zookeys 2016: 81-99.) and different factors have been attributed as drivers of such patterns. For example, environmental factors, such as geographic constraints, indicate that they have played an important role in the compositional gradient of tenebrionids throughout the European continent (Fattorini & Baselga 2012FATTORINI S & BASELGA A. 2012. Species richness and turnover patterns in european tenebrionid beetles. Insect Conserv Divers 5: 331-345.). Environmental filters coupled with spatially structured processes are responsible for the spatial substitution of Cryptocephalinae and Chrysomelinae (Freijeiro & Baselga 2016FREIJEIRO A & BASELGA A. 2016. Spatial and environmental correlates of species richness and turnover patterns in European cryptocephaline and chrysomeline beetles. Zookeys 2016: 81-99.) throughout the European continent. However, these patterns have been observed in studies at wide spatial scales and that present structural gradients of the entire landscape, encompassing wide structural variation and land use.

In fragmented landscapes, where the fragmentation process was preceded by intense habitat loss, the spatial substitution of species could be the result of several factors such as decreased quantity and quality of habitat and the isolation of remaining habitat patches. Habitat loss and fragmentation may lead to the replacement of original fauna, where initially specialist species are replaced by more generalist species (Morante-Filho et al. 2016MORANTE-FILHO JC, ARROYO-RODRÍGUEZ V & FARIA D. 2016. Patterns and predictors of β-diversity in the fragmented Brazilian Atlantic forest: A multiscale analysis of forest specialist and generalist birds. J Anim Ecol 85: 240-250.). In degraded landscapes, more specialist species tend to respond negatively to fragmentation and disturbance (Devictor et al. 2008DEVICTOR V, JULLIARD R & JIGUET F. 2008. Distribution of specialist and generalist species along spatial gradients of habitat disturbance and fragmentation. Oikos 117: 507-514.). Regarding the isolation of areas, this fact can lead to a difference in the composition of local species due to differences in the dispersal capacity of different species. Isolation may also lead to low genetic variability within populations and to the local extinction of some species.

In our study system, fragments of forest remnants are immersed in a matrix dominated largely by agriculture and livestock, and thus constantly suffer pressures and disturbances caused by these production systems. In addition to the pressure generated by the intense use of agricultural pesticides, these fragments present small areas, are isolated from each other and suffer intense and uncontrolled extirpation of their trees (for wood) and fauna (hunting). Such historical factors may have contributed to the species substitution patterns that we found, since they have significantly altered the landscape overall, both structurally and in terms of land use. Thus, local factors more restricted to the habitat patch, such as structural factors (size, shape, insulation, etc.) and wider ones related to landscape characteristics (amount of vegetation cover, number of patches, spatial heterogeneity, different soil uses etc.) may be causing observed patterns. Different landscape patterns can interfere in both ecosystem services (Duarte et al. 2018DUARTE GT, SANTOS PM, CORNELISSEN TG, RIBEIRO MC & PAGLIA AP. 2018. The effects of landscape patterns on ecosystem services: meta-analyses of landscape services. Landsc Ecol 33: 1247-1257.) and insect communities (Rösch et al. 2013RÖSCH V, TSCHARNTKE T, SCHERBER C & BATÁRY P. 2013. Landscape composition, connectivity and fragment size drive effects of grassland fragmentation on insect communities. J Appl Ecol 50: 387-394.).

Based on our results, we speculate that other variables are causing such high species replacement among fragments. These variables could be related to structural factors of the habitat, such as insulation, shape, amount of edge, etc., as well as characteristics related to landscape characteristics such as habitat quantity, soil use, etc. Thus, in order to elucidate the processes that are structuring Chrysomelidae communities in forest fragments, we suggest that variables related to the local structure at the habitat patch level and variables related to broader scales such as composition and heterogeneity of the surrounding landscape should be considered.

ACKNOWLEGMENTS

We thank the Programa de Pós-Graduação em Entomologia e Conservação da Biodiversidade (Universidade Federal da Grande Dourados) to which TST was a member of during the execution of this work; Dr. Jelly Makoto Nakagaki (Universidade Estadual de Mato Grosso do Sul) for providing the equipment and laboratory at the Centro de Pesquisa em Biodiversidade; Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brazil (CAPES) for TST’s scholarship (process no. 88882.182376 / 2018-01), and FVN scholarship (88882.317337/2019-01) and for partially financing this study (Finance code 001).

REFERENCES

  • AMARAL JAM DO, MOTCHI EP, OLIVEIRA H DE, CARVALHO FILHO A, NAIME UJ & DOS SANTOS RD. 2000. Levantamento semidetalhado dos solos do campo experimental de Dourados, da Embrapa Agropecuária Oeste, município de Dourados, MS. Embrapa Agropecuária Oeste, Dourados.
  • ANDREW NR & HUGHES L. 2004. Species diversity and structure of phytophagous beetle assemblages along a latitudinal gradient: predicting the potential impacts of climate change. Ecol Entomol 29: 527-542.
  • BASELGA A. 2010. Partitioning the turnover and nestedness components of beta diversity. Glob Ecol Biogeogr 19: 134-143.
  • BASELGA A, ORME D, VILLEGER S, DE BORTOLI J & LEPRIEUR F. 2018. betapart: Partitioning Beta Diversity into Turnover and Nestdness Components. R package version: 1.5.1. https://CRAN.R-project.org/package=betapart
    » https://CRAN.R-project.org/package=betapart
  • BRASIL LS, OLIVEIRA-JÚNIOR JM, CALVÃO LB, CARVALHO FG, MONTEIRO-JÚNIOR CS, DIAS-SILVA K & JUEN L. 2018. Spatial, biogeographic and environmental predictors of diversity of Amazonian Zygoptera. Insect Conserv Diver 11(2): 174-184.
  • BOUCHARD P, GREBENNIKOV VV, SMITH ABT & DOUGLAS H. 2009. Diversity of Coleoptera. In: Foottit R & Adler P (Eds), Insect Biodiversity: Science and Society. Blackwell Publishing, p. 265-301.
  • DEVICTOR V, JULLIARD R & JIGUET F. 2008. Distribution of specialist and generalist species along spatial gradients of habitat disturbance and fragmentation. Oikos 117: 507-514.
  • DIDHAM RK, GHAZOUL J, STORK NE & DAVIS AJ. 1996. Insects in fragmented forests: a funcional approach. Tree 11: 225-260.
  • DRAY S ET AL. 2019. adespatial: Multivariate Multiscale Spatial Analysis. R package version: 0.3-4. https://CRAN.R-project.org/package=adespatial
    » https://CRAN.R-project.org/package=adespatial
  • DRAY S ET AL. 2012. Community ecology in the age of multivariate spatial analysis. Ecol Monogr 82: 257-275.
  • DUARTE GT, SANTOS PM, CORNELISSEN TG, RIBEIRO MC & PAGLIA AP. 2018. The effects of landscape patterns on ecosystem services: meta-analyses of landscape services. Landsc Ecol 33: 1247-1257.
  • FAHRIG L ET AL. 2019. Is habitat fragmentation bad for biodiversity? Biol Conserv 230: 179-186.
  • FATTORINI S & BASELGA A. 2012. Species richness and turnover patterns in european tenebrionid beetles. Insect Conserv Divers 5: 331-345.
  • FELINKS B, PARDINI R, DIXO M, FOLLNER K, METZGER JP & HENLE K. 2011. Effects of species turnover on reserve site selection in a fragmented landscape. Biodivers Conserv 20: 1057-1072.
  • FERNÁNDEZ-CHACÓN A, STEFANESCU C, GENOVART M, NICHOLS JD, HINES JE, PÁRAMO F, TURCO M & ORO D. 2014. Determinants of extinction-colonization dynamics in Mediterranean butterflies: The role of landscape, climate and local habitat features. J Anim Ecol 83: 276-285.
  • FIETZ CR & FISCH GF. 2008. O Clima da Região de Dourados, MS, 2nd ed., Embrapa Agropecuária Oeste, Documentos, 92, 32 p.
  • FREIJEIRO A & BASELGA A. 2016. Spatial and environmental correlates of species richness and turnover patterns in European cryptocephaline and chrysomeline beetles. Zookeys 2016: 81-99.
  • GAÑÁN I, BASELGA A & NOVOA F. 2008. Diversity patterns in Iberian Calathus (Coleoptera, Carabidae: Harpalinae): species turnover shows a story overlooked by species richness. Environ Entomol 37: 1488-1497.
  • HADDAD NM ET AL. 2015. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci Adv 1: 1-9.
  • HARRISON S, ROSS SJ AND LAWTON JH. 1992. Beta Diversity on Geographic Gradients in Britain. J Anim Ecol 61: 151-158.
  • JOLIVET P. 2015. Together with 30 years of Symposia on Chrysomelidae! Memories and personal reflections on what we know more about leaf beetles. Zookeys 547: 35-61.
  • LEGENDRE P. 2008. Studying beta diversity: ecological variation partitioning by multiple regression and canonical analysis. J Plant Ecol 1: 3-8.
  • LEGENDRE P. 2014. Interpreting the replacement and richness difference components of beta diversity. Glob Ecol Biogeogr 23: 1324-1334.
  • LEGENDRE P, BORCARD D & ROBERTS DW. 2012. Variation partitioning involving orthogonal spatial eigenfunction submodels. Ecology 93: 1234-1240.
  • LENNON JJ, KOLEFF P, GREENWOOD JJD & GASTON KJ. 2001. The geographical structure of British bird distributions : diversity , spatial turnover and scale. J Anim Ecol 70: 966-979.
  • LEPRIEUR F, TEDESCO PA, HUGUENY B, BEAUCHARD O, DÜRR HH, BROSSE S & OBERDORFF T. 2011. Partitioning global patterns of freshwater fish beta diversity reveals contrasting signatures of past climate changes. Ecol Lett 14: 325-334.
  • LINZMEIER AM & RIBEIRO-COSTA CS. 2009. Spatio-temporal dynamics of Alticini (Coleoptera, Chrysomelidae) in a fragment of Araucaria Forest in the state of Parana, Brazil. Rev Bras Entomol 53(2): 294-299.
  • MARTINS SV. 2001. Recuperação de matas ciliares. Aprenda Fácil: Viçosa, 219 p.
  • MORANTE-FILHO JC, ARROYO-RODRÍGUEZ V & FARIA D. 2016. Patterns and predictors of β-diversity in the fragmented Brazilian Atlantic forest: A multiscale analysis of forest specialist and generalist birds. J Anim Ecol 85: 240-250.
  • OKSANEN J ET AL. 2018. vegan: Community Ecology Package. R package version: 2.5-1. https: //cran.r-project.org/package=vegan
    » https: //cran.r-project.org/package=vegan
  • PIMENTA M & DE MARCO PJ. 2015. Leaf Beetle (Chrysomelidae: Coleoptera) Assemblages in a Mosaic of Natural and Altered Areas in the Brazilian Cerrado. Neotrop Entomol 44: 242-255.
  • R CORE TEAM. 2019. R: A language and environment for statistical computing. https://www.r-project.org/
    » https://www.r-project.org/
  • RIZZINI CT. 1997. Tratado de Fitogeografia do Brasil: aspectos ecológicos, sociológicos e florísticos, 2nd ed., Rio de Janeiro: Âmbito Cultural Edições Ltda.
  • RÖSCH V, TSCHARNTKE T, SCHERBER C & BATÁRY P. 2013. Landscape composition, connectivity and fragment size drive effects of grassland fragmentation on insect communities. J Appl Ecol 50: 387-394.
  • SÁNCHEZ-REYES UJ, NIÑO-MALDONADO S, CLARK SM, BARRIENTOS-LOZANO L & ALMAGUER-SIERRA P. 2019. Successional and seasonal changes of leaf beetles and their indicator value in a fragmented low thorn forest of northeastern Mexico (Coleoptera, Chrysomelidae). Zookeys 825: 71-103.
  • SANDOVAL-BECERRA FM, SÁNCHEZ-REYES UJ, CLARK SM, VENEGAS-BARRERA CS, HORTA-VEGA JV & NIÑO-MALDONADO S. 2018. Influence of Habitat Heterogeneity on Structure and Composition of a Chrysomelidae (Coleoptera) Assemblage in a Temperate Forest in Northeast Mexico. Southwest Entomol 43(1): 115-130.
  • SEKERKA L, LINZMEIER AM, MOURA LA, RIBEIRO-COSTA CS, AGRAIN F, CHAMORRO ML, MANFIO D, MORSE GE & REGALIN R. 2017. Chrysomelidae. In: Catálogo Taxonômico da Fauna do Brasil. http://fauna.jbrj.gov.br/fauna/faunadobrasil/115540
    » http://fauna.jbrj.gov.br/fauna/faunadobrasil/115540
  • TELES TS, RIBEIRO DB, RAIZER J & LINZMEIER AM. 2019. Richness of Chrysomelidae (Coleoptera) depends on the area and habitat structure in semideciduous forest remnants. Iheringia Ser Zool 109: 1-8.
  • TOWNES H. 1972. A light-weight malaise trap. Entomol News 83: 239-247.
  • TSCHARNTKE T ET AL. 2012. Landscape moderation of biodiversity patterns and processes - eight hypotheses. Biol Rev 87: 661-685.
  • WHITTAKER RH. 1960. Vegetation of the Sisiyou Mountains, Orgeon and California. Ecol Monogr 30: 279-338
  • WILLIAMS PH. 1996. Mapping variations in the strength and breadth of biogeographic transition zones using species turnover. Proc R Soc B Biol Sci 263: 579-588.

Publication Dates

  • Publication in this collection
    06 Nov 2020
  • Date of issue
    2020

History

  • Received
    1 July 2019
  • Accepted
    28 Oct 2019
Academia Brasileira de Ciências Rua Anfilófio de Carvalho, 29, 3º andar, 20030-060 Rio de Janeiro RJ Brasil, Tel: +55 21 3907-8100 - Rio de Janeiro - RJ - Brazil
E-mail: aabc@abc.org.br