Skip to main content

Analysis of functional and phylogenetic diversity of bat communities: insights from resource use in the Pantanal and Cerrado

Abstract

Background

Assessing the structure of communities requires analysis of multiple dimensions of biological diversity. Such approaches provide a broader understanding of the ecological and historical factors involved in the formation and maintenance of communities. Metrics such as functional and phylogenetic diversity are crucial as they unveil how communities respond to environmental changes, providing a deeper understanding of their current state and resilience. This goes beyond species richness or community composition, enabling more profound insights into their present condition and resilience. Therefore, analysis in this sense allows accessing differences in the distribution of species, considering different landscape mosaics, and allowing strategic prioritization of conservation initiatives. In this sense, we investigated whether the composition of functional and phylogenetic groups of chiropteran species differs between the Pantanal and Cerrado biomes and various habitat types. We hypothesize that communities with greater conservation integrity are more diverse than communities that suffer anthropogenic influences.

Methods

The surveyed communities include two areas in the Pantanal and two in the Cerrado, with one conservation unit and one that suffers anthropic influence, for each biome. We assessed the composition of bat assemblages by examining functional richness (FRic), functional uniformity (FEve), functional divergence (FDiv), and functional dispersion (FDis). Metrics that reflect phylogenetic richness and divergence were also used, such as phylogenetic diversity (PD), mean pairwise distance (MPD), and mean nearest taxon distance (MNTD).

Results and discussion

Overall, conservation units sustain higher levels of functional and phylogenetic diversity compared to areas with anthropic influence. In the latter, species occupy fewer ecological niches, indicating that bats can persist in degraded environments, albeit the communities support only a portion of the functional groups.

Conclusion

Our results showed, through analyses across multiple dimensions of functional and phylogenetic diversity, that bat populations are impacted by environmental degradation and fragmentation. Communities in more conserved areas exhibited higher functional and phylogenetic richness, indicating a greater number of occupied niches. These findings highlight the importance of exploring measures of biological diversity and their spatial scales to advance our understanding of biodiversity dynamics and optimize conservation planning for both species and communities.

Background

Studies considering the relationships between landscape structure and patterns of habitat use by chiropterans remain limited and predominantly focus on the consequences of fragmentation [1,2,3,4]. Studies on landscape mosaic structure and evolutionary history have been developed to describe geographic gradients in species richness [5,6,7]. However, interest in phylogenetic and functional diversity has greatly increased in recent years, considering different taxonomic groups [6].

The multifaceted roles bats play in resource exploitation are better understood through understanding the structure and composition of their functional groups [8], incorporating various perspectives such as taxonomic, evolutionary, and ecological characteristics. Some classifications, derived from diverse methodologies [9,10,11], often categorize bats into groups such as aerial insectivores, trawling piscivore-insectivores, and gleaning carnivores, among others [11,12,13]. These foraging categories have contributed to a better understanding of the niches occupied by bats, thereby informing us about which species are more susceptible to deforestation [14]. Evaluating groups from a functional perspective has proven to be more effective than narrower assessments [15]. Chiropterans are excellent as a study group due to their wide range of dietary preferences, shelter uses, foraging behaviors, and ecological roles. They provide multiple ecosystem services including pollination, seed dispersal, and insect population control, and they are indicators of toxin accumulation [16, 17].

Metrics of functional and phylogenetic diversity provide valuable insights into community responses to environmental changes, offering a more comprehensive understanding of resilience and functional integrity than species richness or community composition alone [18, 19]. Assessing community phylogenetic diversity, especially by identifying regions rich in evolutionary history, deepens our understanding of structural constraints across different landscapes [20]. This framework suggests that areas with a higher number of phylogenetic branches tend to exhibit greater functional trait diversity, enhancing their potential for evolutionary adaptation to environmental changes [20, 21].

At the core of the concept of functional diversity is the role of functional traits, defined as the characteristics of organisms that influence both their fitness and the functioning of ecosystems [22,23,24]. The accuracy of functional diversity assessments largely depends on the selection of ecologically important traits [22,23,24]. In this context, functional diversity entails comparing trait differences among species and individuals to assess their effects on ecosystem function and their responses to environmental changes [23, 25]. Such analyses of functional diversity help quantify the elements of biodiversity that are essential for ecosystem operation [26].

The mammal fauna of the Pantanal region has been documented in several inventories [4, 27,28,29,30,31,32,33]. However, chiropteran diversity may be underestimated, especially due to differences in species composition between flooded and non-flooded areas. Some species are exclusively found in the floodplain, while the greatest richness of chiropterans is noted in the upland areas adjacent to it [33]. In the Cerrado, chiropteran surveys have identified between 11 and 25 species [34,35,36,37,38,39,40,41,42], which represent approximately 14% of Brazil’s total chiropteran species [43]. This limited diversity is likely influenced by extensive deforestation, environmental degradation, habitat fragmentation, land use for agriculture, and livestock farming, which are the main factors influencing biological diversity in the Pantanal and Cerrado biomes [44,45,46,47].

Previous studies have demonstrated that communities in conserved areas exhibit greater heterogeneity in functional and phylogenetic groups compared to altered areas [48,49,50,51]. Research on the phylogenetic and functional relationships within bat communities remains scarce in these biomes. However, investigations into functional and phylogenetic diversity in the Cerrado have shown greater diversity than previously recorded [52]. Understanding the impacts of anthropogenic actions on species diversity and functional groups is crucial for developing effective landscape management guidelines, considering species persistence and resilience in these areas. In this study, we evaluated this hypothesis, in the context of the Pantanal and Cerrado regions.

Methods

Study area

In the Pantanal, western Brazil, the Reserva Particular do Patrimônio Natural (RPPN) Sesc Pantanal (SescPant) covering 108,000 hectares in Barão de Melgaço, Mato Grosso (16°41’S, 56°24’W), is characterized by a diverse mosaic of savannas and forest formations, including fields with murundus (rounded earthmounds covered by woody vegetation), shrubby areas, seasonally dry forests, and cambará forests (monospecific formation dominated by Vochysia divergens, Vochysiaceae). The cambará forest divides large open areas, characterizing a diversity of habitats within the conservation unit [53, 54]. In contrast, the Santa Lúcia Ranch (RanchPant), in the adjacent region (16°53’S, 55°54’W), left bank of the São Lourenço River, is more altered by human activity and cattle raising. The region comprises wet savanna with forest-islands (murundus), grasslands, and areas of sparse vegetation widely used as pasture. Denser forests along the rivers, interspersed with cambarazal forests and acurizais (forests of Scheelea phalerata palm), show a clear impact of livestock activity on the environmental heterogeneity [13] (Fig. 1).

Fig. 1
figure 1

The main map shows the neotropical region, highlighting the four sampled locations in Brazil. Smaller maps illustrate the communities in the Pantanal: (1) Sesc Pantanal (SescPant) and Santa Lúcia Ranch (RanchPant), Cerrado communities (2) Parque Sesc Serra Azul (SerAzul), and (3) Serra da Mesa (SerMesa). In 1 and 2 the polygons delimit the conservation units

In the Cerrado region, the Parque Sesc Serra Azul (SerAzul), Rosário Oeste, Mato Grosso (14º29’S, 55º44’W), is a conservation unit that spans approximately 5,000 hectares [55]. The region is a mix of dry and floodable pastures, with or without sparse trees, and various types of seasonal forests, including gallery forests and forests with babassu palms (Orbignya phalerata). Human influence in the surrounding areas, as in many other regions of the Cerrado, has resulted in a mosaic of agricultural exploitation interspersed with original vegetation [47, 56, 57]. On the other hand, the Serra da Mesa (SerMesa) region, encompassing the municipalities of Niquelândia, Colinas do Sul, Minaçu, and Uruaçu in Goiás (13°45’ – 14°35’S and 47°50’ – 49°15’W), includes a variety of landscapes ranging from open fields to dense, humid forests along watercourses. The area reflects the typical mesophytic habitats of the Cerrado core, with cerrado sensu stricto and cerradão (woodland savanna), contrasting with gallery forests that usually support a diversity of plant and animal species [56, 57] (Fig. 1).

The sampling periods lasted between 10 and 15 nights and were conducted using mist nets. The number of nets deployed per night ranged from 5 to 7, and each site was sampled for 1 to 3 days, with a sampling duration of six hours following sunset. The nets were checked every 15 to 30 min, depending on the number of animals captured. The nets were primarily set up in clearings or open corridors within forested areas, in locations with flowering or fruiting trees, and near rivers or bodies of water. All data presented in this study were obtained from previous sampling by the authors [4, 13, 55]. Sampling was conducted in four distinct areas over different periods. In the Pantanal region, data from RPPN Sesc Pantanal were collected between 1999 and 2008, covering all months of the year, while collections at Santa Lúcia Ranch took place in November and December 2014 and October 2015. In the Cerrado communities, sampling was conducted at Parque Sesc Serra Azul in August 2012, as well as in March and August 2013. In the Serra da Mesa region, data were gathered in September, October, and December 1995, in addition to February and May 1996. To ensure comparability across communities and adjust for sampling effort variations, we standardized the data by normalizing each value relative to the total captures per community, using the ‘decostand’ function along with the ‘total’ method from the ‘vegan’ package [58].

Sample coverage

The sampling effectiveness was evaluated using Hill numbers, which are a measure of taxonomic species diversity within communities, based on the methodology detailed by Chao et al. [59]. Hill numbers (the effective number of species) are used to quantify species diversity by considering both species richness (the number of species) and their relative abundance. When q = 0, the focus is solely on species richness, meaning the number of distinct species present is counted without regard to how abundant each species is. Using q = 0 is particularly important in areas where the capture effort is unknown (one of the study areas), as the data are based on a substantial faunal survey effort that is not fully available. Although the data offers valuable information, it has limitations. However, it still allows for inferences about the communities and enables comparisons by equalizing sample sizes. The ‘iNEXT’ package [60] was utilized for this analysis, facilitating the computation and visualization of rarefaction and extrapolation sampling curves. This package also provided 95% confidence intervals and allowed us to estimate the asymptotic diversity values. This approach was particularly important given the different sampling efforts employed across the communities, ensuring that our diversity estimates accurately reflect the varying levels of sampling intensity.

Molecular data

The MEGA software [61] was used for aligning the sequences and manually editing alignments. Maximum likelihood gene tree estimations were carried out with IQ-TREE software [62], and the best-fit model was GTR + F + I + G4 chosen according to the BIC identified by ModelFinder [63] included in the IQ-TREE. Branch support was assessed using an approximate likelihood-ratio test based on the Shimodaira-Hasegawa-like procedure (SH-aLRT) [64] and nonparametric bootstrap frequencies. The tree was edited with FigTree software [65]. Overall, the bootstrap and SH-aLRT values were above 80, with only a few clades showing support values under this threshold.

A total of 70 sequences were used, 42 from the present study and 28 obtained from GenBank (Table S1). As an external group, sequences from GenBank of the bats Rousettus aegyptiacus (AB085740), Rhinolophus ferrumequinum (AB085731), and the Canidae Cerdocyon thous (KU253528), with 1,140 base pairs, were included. Neoplatymops mattogrossensis was removed from the phylogenetic diversity analysis due to a lack of samples for comparison. Details of the amplification conditions can be found in Table S2. The list of primers used in the amplification is available in Table S3. Additional information on the deoxyribonucleic acid (DNA) extraction and sequencing methods (cytochrome b gene) is described in detail in the study by de Oliveira et al. [66].

Phylogenetic diversity analyzes

Based on the phylogeny of the rooted maximum likelihood tree, we calculated: phylogenetic diversity (PD) [67], mean pairwise distance (MPD) [68], and mean nearest taxon distance (MNTD) [68].

The PD measure is the sum of the length of all branches in a phylogeny [68], and represents the sum of the evolutionary history in a sample; it is considered the first phylogenetic diversity metric. However, it does not consider species abundance within the sample; Faith’s PD is expected to increase as the number of species in a sample rises.

The MPD measure quantifies the average pairwise phylogenetic distance between all species in an assemblage [22, 68]. MPD is a metric that weights the basal structure of the phylogeny (e.g., relationships between species from different families) [69]. The MNTD measure quantifies the pairwise phylogenetic distance between the closest relatives in an assemblage [22, 68]. It is a terminal metric that weights relationships at the ends of the phylogeny (e.g., species within the same genus) [69]. High values of MPD indicate that phylogenetically grouped. High values of MNTD suggest that closely related species do not co-occur in the community, that is, some species have branches much longer than average. On the other hand, low values suggest co-occurrence.

To eliminate metric redundancy, phylogenetic diversity was quantified using standardized effect size (SES) [22]. To calculate SES, species richness for each site was fixed, and species at the ends of the phylogeny were randomized 1000 times. A significant positive SES value indicates an over-dispersed phylogenetic structure, while significant negative values indicate phylogenetic clustering [69].

Phylogenetic signal

To assess whether functional traits present a phylogenetic signal, we used Blomberg’s K and Significance Tests. These tests quantify the extent to which trait variation is explained by the structure of a given phylogenetic tree. To compare K values across different phylogenies and studies, standardization is required [22, 70]. This value is based on the expected value derived from the evolution of the trait under a Brownian motion model within the observed phylogenetic tree [70].

Characterization of food guild and foraging strategy as functional traits

To better understand bat assemblages, a species-level classification was employed based on diet and foraging behavior [71]. This classification system was based on the analysis of existing literature and included items such as seeds, fruits, nectar, invertebrates, vertebrates, fish, and blood. Each species was then categorized based on its predominant, secondary, and tertiary diets when applicable (Table S4). Some bat species may consume more than one type of resource, either regularly or occasionally, depending on the availability of the predominant resource. These variations may occur seasonally, as documented in various studies [72,73,74].

The foraging strategies have been characterized into five categories: aerial, gleaning, sedentary, nomadic, and trawling (Table S4). Aerial foragers capture prey while flying. Gleaning involves collecting food by picking it up from a solid surface, such as leaves, branches, or the ground. The trawling foraging strategy involves capturing prey on over water or, in some cases, in water. Nomadic foragers are associated with frugivores, and their populations move to different regions based on resource availability, suggesting more frequent shelter renewal. Sedentary foragers, also associated with a frugivorous diet, follow a regularly fixed route, typically when preference is linked to plants with continuous year-round production. The characterization of foraging strategies followed the propositions of Fenton [75], Soriano [76], Schnitzler and Kalko [12], Schnitzler et al. [11], and Denzinger and Schnitzler [77].

In addition, we performed Principal Coordinates Analysis (PCoA) using the ‘ape’ package to classify species into functional groups based on Bray-Curtis dissimilarities, utilizing the ‘ecodist’ package [78]. This analysis groups species based on their traits, thereby providing a better understanding of community composition and the dynamics of diversity.

Functional diversity analysis

Functional diversity was assessed using the widely known metric ‘FD’ [79]. It was calculated by generating a functional dendrogram from a matrix containing species and their associated functional traits, based on the total branch length. This approach yielded indices for functional richness (FRic), functional uniformity (FEve), functional divergence (FDiv), and functional dispersion (FDis) [22]. Additionally, mean pairwise distance (MPD) and mean nearest taxon distance (MNTD) analyses were also performed for functional diversity [68].

Functional richness (FRic) quantifies the volume of niche space occupied by species based on their traits highlighting the range of ecological roles within a community [80, 81]. It is measure as the minimum convex polygon that encompasses all species and is particularly sensitive to species at the edges of the trait distribution [82]. Functional evenness (FEve) assesses how evenly species abundances are distributed across the occupied niche space, reflecting how efficiently a community uses available resources within its ecological limits [81].

Functional divergence (FDiv) and dispersion (FDis) measure niche differentiation and the average distance of species from the functional centroid, respectively, with high FDiv values indicating significant niche differentiation and high FDis values reflecting broader utilization of niche space edges [81].

Dietary traits were quantified by determining the proportion of various food types consumed by each species when they used multiple food sources. Foraging strategy traits were coded in binary terms. Quantitative traits, such as size and weight, were measured in millimeters and grams. Due to variability in measurement scales, standardization was necessary to mitigate scale effects, as stated in Pla et al. [83].

To ensure the accurate calculation of indices, all trait values were standardized to a range of zero to one using the ‘decostand’ function along with the ‘range’ method from the ‘vegan’ package [58]. This conversion ensured that the minimum value on the original scale becomes zero and the maximum value one, making all traits continuous [81, 83]. The ‘FD’ package was used to calculate functional diversity metrics using the ‘dbFD’ function for distance-based indices [84]. The 3D graphs were generated using the ‘plotly’ package [85]. All diversity analyses were conducted on the R platform [86].

Results

The collected samples varied considerably in size and composition, reflecting the ecological complexity of the analyzed environments. In the Pantanal, the SescPant site recorded 721 captures, including 40 species and 25 genera, while the RanchPant community had 659 captures, representing 29 species and 17 genera. In the Cerrado, the conservation area community SerAzul documented 562 captures with 39 species and 27 genera, and the SerMesa community recorded 490 captures, encompassing 33 species and 27 genera.

The rarefaction curves illustrate different levels of bat species richness among the four communities evaluated (Fig. 2). Within the Pantanal, the reserve community (SescPant) exhibits an extrapolation curve that has not reached an asymptote, indicating the potential for higher species richness than been recorded so far. In contrast, the curve for the ranch community (RanchPant), a more anthropogenically impacted area suggests lower species richness, with sampling approaching completion. Despite the proximity of the ranch and the reserve, which are separated by the banks of the São Lourenço River (with sampling points up to 20 km apart), the differences in bat communities between these areas are striking.

In the Cerrado region, the communities of Parque Sesc Serra Azul (SerAzul) and Serra da Mesa (SerMesa) show comparable species richness (Fig. 2), with both curves indicating that the species richness surpasses the sampled data. After rarefaction, there are no expressive differences in bat species richness across all studied communities in the Pantanal and Cerrado regions. However, within the Pantanal, the contrast between the conservation-oriented SescPant area and the livestock-driven RanchPant, is notably pronounced.

The effectiveness of the sampling effort was assessed by calculating the sample coverage at each community, which was particularly important as the sampling effort in some communities was smaller, requiring us to ensure equivalent effort per unit of area and time. In the Pantanal communities, sample coverage was 98% for SescPant and 99% for RanchPant. In the Cerrado communities, the sample coverage was 98.17% for SerMesa and 97.87% for SerAzul, indicating near-complete sampling for both communities.

Fig. 2
figure 2

Rarefaction curves, represented by solid lines for sample size and dotted lines for extrapolation, include 95% confidence intervals shown as shaded areas. These curves assess the number of individuals captured across various communities, which are categorized based on their conservation status: conserved areas (SescPant, SerAzul) and human-influenced areas (RanchPant, SerMesa). The communities are differentiated based on diversity, specifically q = 0 (species richness)

In the Principal Coordinates Analysis (PCoA), a clustering pattern emerged based on species’ foraging strategies. The foraging strategy delineated the groups, with some sharing similarities due to certain species consuming multiple types of food (Fig. 3).

Fig. 3
figure 3

Principal Coordinates Analysis (PCoA) was performed using Bray-Curtis dissimilarities to classify species into six groups representing foraging strategies: Gleaning (Nomadic and Sedentary), Trawling, and Aerial. Diets were categorized as (1) Frugivore, (2) Piscivore-Insectivore, (3) Insectivore, (4) Sanguivore, (5) Granivore, (6) Omnivore, (7) Nectarivore, (8) Carnivore-Insectivore, and (9) Frugivore-Insectivore

The phylogenetic diversity (PD, Fig. 4a) observed in the SerAzul community is high, with positive sesPD values suggesting that this community has higher PD values than expected for its observed species richness. Conversely, negative sesPD values were observed for the other three communities (SescPant, RanchPant, SerMesa), indicating that their phylogenetic diversity is lower than expected given their species richness.

The MPD metric, which assesses the average phylogenetic distance among all individuals within a community (Fig. 4b), was high in the SescPant and SerAzul communities. This indicates the greatest phylogenetic distance among the species within these communities. Additionally, the positive sesMPD for SescPant suggests that this community has higher MPD values than expected (Table 1; Fig. 4b), showing exceptional phylogenetic diversity when contrasted with the others. The communities with the greatest average distance to the nearest taxon (MNTD) were the two from the Cerrado (SerMesa, SerAzul). However, the positive standardized effect size for MNTD (sesMNTD) in SerAzul indicates that this community has a higher phylogenetic diversity between species than observed (Table 1; Fig. 4c).

Table 1 Phylogenetic and functional diversity indices showing phylogenetic diversity (PD), mean pairwise distance (MPD), mean nearest taxon distance (MNTD) and standardized effect size (sesPD, sesMPD, sesMNTD), calculated for the Pantanal (SescPant, RanchPant) and Cerrado (SerMesa, SerAzul) communities. In bold are values statistically significant, with a p-value lower than 0.05

Mean pairwise distance (MPD) and mean nearest taxon distance (MNTD) (Table 1) for phylogenetic diversity (Fig. 4) and functional diversity (Fig. 5) showed higher diversity in the conservation unit communities (SescPant, SerAzul). This finding is supported by the FRic values, which showed the highest functional richness in these communities (Fig. 6).

Fig. 4
figure 4

Standardized effect size: (a) sesPD, (b) sesMPD, and (c) sesMNTD, along with 5% and 95% quantiles (dashed lines) of the simulated null communities of Faith’s phylogenetic diversity (PD), mean pairwise distance (MPD) and mean nearest taxon distance (MNTD) assessed for the communities sampled in Pantanal (SescPant, RanchPant) and Cerrado (SerMesa, SerAzul)

Blomberg’s K values were tested for all functional traits to evaluate the intensity of the phylogenetic signal of these characteristics in the phylogeny of these communities. The results suggest a phylogenetic signal with a high degree of variation for some phylogenetic traits, such as the insectivore, sanguivore, frugivore, and nectarivore diets, and the foraging strategies gleaning nomadic, gleaning, gleaning sedentary, aerial-trawling, aerial and trawling. Blomberg’s K values show a non-significant or absent phylogenetic signal for the carnivore, piscivore-insectivore, and granivore diets, the same occurring for the forearm and weight traits (Table 2).

Table 2 Blomberg’s K values tested for all functional traits. K represents the intensity of the phylogenetic signal for each functional trait in the phylogeny of the communities. Bold text indicates a significant statistical value with a p-value lower than 0.05

In terms of functional diversity, the high MPD (Fig. 5a) and low MNTD (Fig. 5b) values indicate that all communities exhibit substantial functional differentiation between species pairs, thereby indicating high functional diversity at the species level. However, conservation units exhibit even greater functional diversity compared to communities subjected to more anthropogenic influence.

In the Pantanal, the community experiencing the highest level of anthropogenic influence (RanchPant), exhibited the lowest functional and phylogenetic diversity. The low MPD and MNTD values, which were significant (p) in the context of phylogenetic diversity (Fig. 4), suggest that it is the most phylogenetically impoverished community. Moreover, the significant values of the standardized effect size suggest that the diversity measures may be overestimated, meaning that the phylogenetic diversity values reported by the metrics are higher than expected. Furthermore, the values of functional evenness (FEve), functional divergence (FDiv), and functional dispersion (FDis, Table 3) were also low in this community, indicating few niche spaces occupied, many species closely related functionally, and low dispersal of functional traits.

Table 3 Values of the functional diversity metrics assessed include species functional richness (FRic), functional evenness (FEve), functional divergence (FDiv), and functional dispersion (FDis). These metrics were evaluated for communities sampled in the Pantanal (SescPant, RanchPant) and Cerrado (SerMesa, SerAzul). The observed species richness (N) was recorded for each community

Conversely, the conservation units’ communities in both sampled biomes (SescPant, SerAzul) displayed the greatest richness (FRic), increased functional divergence (FDiv), and greatest functional dispersion (FDis, Fig. 6). These metrics indicate that a wide variety of niches are occupied and traits are broadly distributed within these communities (Table 3).

The community in the Cerrado with the most anthropogenic influence (SerMesa) had the lowest functional richness (FRic) and the second-lowest functional dispersion (FDis, Fig. 6), indicating that few niches are occupied with restricted functional dispersion. However, this community showed the greatest functional uniformity (FEve) and elevated functional divergence (FDiv, Fig. 6), indicating that functional traits are evenly distributed across species, which fulfill a diverse array of ecological functions. These metrics, which represent MPD and MNTD, suggest that most species are closely related phylogenetically. However, phylogenetic divergence is greater among the basal taxa, while phylogenetic distances are smaller between the terminal taxa (Fig. 4b and c).

Fig. 5
figure 5

Standardized effect size: (a) sesMPD, and (b) sesMNTD, along with 5% and 95% quantiles (dashed lines) of the simulated null communities of functional diversity (FD), mean pairwise distance (MPD), and mean nearest taxon distance (MNTD) assessed for the communities sampled in Pantanal (SescPant, RanchPant) and Cerrado (SerMesa, SerAzul)

Fig. 6
figure 6

Functional trait diversity plotted in a 3D graph, assessed for the sampled communities in the Pantanal (SescPant, RanchPant) and Cerrado (SerMesa, SerAzul). Spatial Functional Richness (FRic), Functional Evenness (FEve), Functional Divergence (FDiv), and Functional Dispersion (FDis) were evaluated

Discussion

We tested the hypothesis that conservation areas, which maintain habitat integrity in the Pantanal and Cerrado regions, support bat communities with greater functional and phylogenetic heterogeneity. Our findings validated this hypothesis, as we observed a positive association between habitat integrity and functional and phylogenetic diversity. This association was adjusted for species richness using a standardized effect size (SES). Previous studies have shown that vegetation homogeneity can act as an environmental filter, leading to the dominance of functional and phylogenetic clustering in bat communities [87,88,89]. In contrast, the heterogeneity of vegetation types found in more conserved areas increases phylogenetic diversity [89]. These results were consistent with our conclusions and support the importance of conservation efforts in maintaining habitat integrity and promoting biodiversity.

Land use for ranching, as exemplified by RanchPant in the Pantanal, has led to a marked decline in both functional and phylogenetic diversity within bat communities, as shown by the lower values of phylogenetic diversity indices (Table 1) and reduced functional richness (Table 3). This decline is indicative of a shrinking ecological niche space, which is particularly evident when compared to more conserved areas [6]. Previous studies have highlighted that biodiversity loss has been aggravated by substantial habitat reduction and fragmentation at the landscape level, one of the primary drivers of these observed patterns [87, 90, 91]. In RanchPant, anthropogenic influence has promoted conservative habitat filtering, favoring frugivorous species that are closely related phylogenetically and share similar functional traits, adapted to degradation and environmental fragmentation [6, 92]. This selection has resulted in bat communities that share advantageous traits for survival in degraded conditions, but with impoverished biological diversity compared to less anthropized areas [4, 28].

Despite the pressures of habitat alteration, the present study reveals the prevalence of frugivores and nectarivores in all analyzed communities. This finding corroborates the results of previous studies, which indicate that in contexts of intensified land use, there is an expansion in the abundance of functional traits associated with plant-based diets [93]. This association suggests a potential adaptation of these animals specialized in fruits and nectar to environmental transformations induced by human activity, possibly facilitated by the availability of viable food resources, even in profoundly modified landscapes.

Bat species play crucial roles as seed dispersers and pollinators, especially in the Neotropical region, where many plant species rely on frugivores for dispersal services [94]. Additionally, bats contribute significantly to the pollination of numerous tropical crops, supporting a variety of economically and ecologically important plant species [95]. Therefore, conserved areas are crucial for maintaining the resilience of key species and functional groups that underpin the stability of these biomes. Meanwhile, even degraded areas continue to support constant regeneration, particularly in regions affected by livestock-induced fragmentation [96]. This underscores the importance of bats in both protected and anthropized areas in supporting essential ecological functions.

Phylogenetic signal analyses indicate limited variation in the evolution of traits such as forearm length, body mass, and diet (carnivore, granivore and piscivore-insectivore), with a low phylogenetic signal observed when testing the value of K. This suggests that only a small portion of the variation in these functional traits can generally be attributed to adaptive selection occurring after the divergence of a common ancestor [97]. The bat superfamily Noctilionoidea exhibits a remarkable ecological breadth, with divergent echolocation, specialized flight strategies linked to insectivory, and parallel evolution of diverse plant-based diets, such as nectar, pollen, and fruits, evolving from the ancestral insectivorous condition [98]. Frugivorous, nectarivorous, and insectivorous conditions appear to evolve before carnivorous and hematophagous ones in bat phylogeny. The insectivorous condition within the order is considered the oldest, with all other guilds evolving from this initial condition [98]. However, the lack of data on closely related community traits may reduce the ability to detect a phylogenetic signal for functional traits, especially those that have evolved recently [99].

Conclusions

Our study focuses on the biodiversity of bat communities in the Pantanal and Cerrado regions and how it is affected by local-scale land use. We found that degraded areas have lower functional and phylogenetic diversity, which is consistent with previous studies confirming the sensitivity of bats to local-scale deforestation [51, 87, 100, 101]. However, there is still a lack of research on the mechanisms behind these changes and the impact of habitat complexity and disturbance on bat diversity.

Therefore, our results emphasize the importance of integrating both functional and phylogenetic dimensions into conservation strategies to ensure the preservation of species, their evolutionary histories, and the functional roles they play, which are essential for sustaining ecosystem functions and services. They also reinforce that the effects of different forms of landscape use should not only be viewed from the perspective of species arrangements (losses or gains) but also under the traits of their components, with potential implications for the functionality of ecological systems, particularly through the recognized contributions of bats in regeneration processes. Within each bat community in the Pantanal and Cerrado, functional and phylogenetic diversity is shaped by environmental filtering and limiting similarity pressures. In protected areas, environmental filtering supports greater functional and phylogenetic diversity, while in anthropized areas, limiting similarity pressures and the effects of degradation tend to reduce this diversity.

In summary, our research uncovered strong spatial patterns in the functional diversity of bat communities across two biomes, when contrasting community structures in more intact and anthropized areas. This suggests that anthropogenic pressures not only reduce species richness but also homogenize the functional attributes of communities, potentially compromising ecosystem resilience and functionality. Such insights are critical for developing and implementing conservation programs aimed at preserving biodiversity.

Data availability

All data generated or analyzed during this study are included in this published article.

Abbreviations

FRic:

Functional richness

FEve:

Functional uniformity

FDiv:

Functional divergence

FDis:

Functional dispersion

PD:

Phylogenetic diversity

MPD:

Mean pairwise distance

MNTD:

Mean nearest taxon distance

SES:

Standardized effect size

sesPD:

Standardized effect size phylogenetic diversity

sesMPD:

Standardized effect size mean pairwise distance

sesMNTD:

Standardized effect size mean nearest taxon distance

PCoA:

Principal Coordinates Analysis

DNA:

Deoxyribonucleic acid

SescPant:

Sesc Pantanal

RanchPant:

Santa Lúcia Ranch

SerAzul:

Parque Sesc Serra Azul

SerMesa:

Serra da Mesa

RPPN:

Reserva Particular do Patrimônio Natural

References

  1. Aguirre LF. Structure of a neotropical savanna bat community. J Mammal. 2002;83(3):775–84. https://doiorg.publicaciones.saludcastillayleon.es/10.1644/1545-1542(2002)0830775:SOANSB2.0.CO;2.

    Article  Google Scholar 

  2. Bobrowiec PED, Gribel R. Effects of different secondary vegetation types on bat community composition in Central Amazonia, Brazil. Anim Conserv. 2010;13(2):204–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1469-1795.2009.00322.x.

    Article  Google Scholar 

  3. Ávila-Gómez ES, Moreno CE, García-Morales R, Zuria I, Sánchez-Rojas G, Briones-Salas M. Deforestation thresholds for phyllostomid bat populations in tropical landscapes in the Huasteca region, Mexico. Trop Conserv Sci. 2015;8(3):646–61. https://doiorg.publicaciones.saludcastillayleon.es/10.1177/194008291500800305.

    Article  Google Scholar 

  4. de Oliveira MB, Brandão ML, Cordeiro JLP, de Oliveira LFB, Peracchi AL. Composition and frequency of capture of bat assemblages in a landscape mosaic in Northern Pantanal, Mato Grosso, Brazil. Mastozool Neotrop. 2020;27(1):39–52. https://doiorg.publicaciones.saludcastillayleon.es/10.31687/saremMN.20.27.1.0.25.

    Article  Google Scholar 

  5. Brum FT, Kindel A, Hartz SM, Duarte LD. Spatial and phylogenetic structure drive frugivory in Tyrannidae birds across the range of Brazilian Araucaria forests. Oikos. 2012;121(6):899–906. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1600-0706.2011.19978.x.

    Article  Google Scholar 

  6. Cosset CC, Edwards DP. The effects of restoring logged tropical forests on avian phylogenetic and functional diversity. Ecol Appl. 2017;27(6):1932–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/eap.1578.

    Article  PubMed  Google Scholar 

  7. Grimshaw JR, Higgins CL. Environmental correlates of phylogenetic structure in Mexican bat communities. J Mammal. 2017;98(6):1657–66. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/jmammal/gyx101.

    Article  Google Scholar 

  8. Stevens RD, Cox SB, Strauss RE, Willig MR. Patterns of functional diversity across an extensive environmental gradient: vertebrate consumers, hidden treatments and latitudinal trends. Ecol Lett. 2003;6(12):1099–108. https://doiorg.publicaciones.saludcastillayleon.es/10.1046/j.1461-0248.2003.00541.x.

    Article  Google Scholar 

  9. Fleming TH. Opportunism versus specialization: the evolution of feeding strategies in frugivorous bats. In: Estrada A, Fleming TH, editors. Frugivores and seed dispersal. Dr W. Junk, Springer Dordrecht;; 1986. pp. 105–18.

  10. Fenton MB. The foraging behaviour and ecology of animal-eating bats. Can J Zool. 1990;68(3):411–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1139/z90-061.

    Article  Google Scholar 

  11. Schnitzler HU, Moss CF, Denzinger A. From spatial orientation to food acquisition in echolocating bats. Trends Ecol Evol. 2003;18(8):386–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/S0169-5347(03)00185-X.

    Article  Google Scholar 

  12. Schnitzler HU, Kalko EK. Echolocation by insect-eating bats: we define four distinct functional groups of bats and find differences in signal structure that correlate with the typical echolocation tasks faced by each group. Biosci. 2001;51(7):557–69. https://doiorg.publicaciones.saludcastillayleon.es/10.1641/0006-3568(2001)0515b0557:EBIEB5d2.0.CO;2.

    Article  Google Scholar 

  13. de Oliveira MB, Brandão ML, Cordeiro JLP, de Oliveira LFB, Peracchi AL. Trophic structure and foraging strategies in a bat community in northern Pantanal. Brazil Mammalia. 2022;87(1):41–51. https://doiorg.publicaciones.saludcastillayleon.es/10.1515/mammalia-2021-0196.

    Article  Google Scholar 

  14. Medellín RA, Equihua M, Amin MA. Bat diversity and abundance as indicators of disturbance in neotropical rainforests. Conserv Biol. 2000;14(6):1666–75. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1523-1739.2000.99068.x.

    Article  PubMed  Google Scholar 

  15. Gonçalves F, Fischer E, Dirzo R. Forest conversion to cattle ranching differentially affects taxonomic and functional groups of neotropical bats. Biol Conserv. 2017;210(Part A):343–8. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.biocon.2017.04.021.

    Article  Google Scholar 

  16. Aguiar L, Antonini Y. Diet of two sympatric insectivores bats (Chiroptera: Vespertilionidae in the Cerrado of Central Brazil. Rev Bras Zool. 2008;25(1):28–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S0101-81752008000100005.

    Article  Google Scholar 

  17. Jones G, Jacobs DS, Kunz TH, Willig MR, Racey PA. Carpe noctem: the importance of bats as bioindicators. Endanger Species Res. 2009;8(1–2):93–115. https://doiorg.publicaciones.saludcastillayleon.es/10.3354/esr00182.

    Article  Google Scholar 

  18. Violle C, Navas ML, Vile D, Kazakou E, Fortunel C, Hummel I, Garnier E. Let the concept of trait be functional! Oikos. 2007;116(5):882–92. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.2007.0030-1299.15559.x.

    Article  Google Scholar 

  19. Ricotta C, Godefroid S, Celesti-Grapow L. Common species have lower taxonomic diversity evidence from the urban floras of Brussels and Rome. Divers Distrib. 2008;14(3):530–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1472-4642.2008.00467.x.

    Article  Google Scholar 

  20. Forest F, Grenyer R, Rouget M, Davies JT, Cowling RM, Faith DP, Balmford A, Manning JC, Procheş S, van der Bank M, Reeves G, Hedderson TAJ, Savolainen V. Preserving the evolutionary potential of floras in biodiversity hotspots. Nature. 2007;445(7129):757. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nature05587.

    Article  CAS  PubMed  Google Scholar 

  21. Safi K, Cianciaruso MV, Loyola RD, Brito D, Armour-Marshall K, Diniz-Filho JAF. Understanding global patterns of mammalian functional and phylogenetic diversity. Philos Trans R Soc Lond B Biol Sci. 2011;366(1577):2536–44. https://doiorg.publicaciones.saludcastillayleon.es/10.1098/rstb.2011.0024.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Swenson NG. Functional and phylogenetic ecology in R. New York: Springer; 2014. p. 212.

    Book  Google Scholar 

  23. Petchey OL, Gaston KJ. Functional diversity: back to basics and looking forward. Ecol Lett. 2006;9(6):741–58. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1461-0248.2006.00924.x.

    Article  PubMed  Google Scholar 

  24. Cianciaruso MV, Batalha MA, Gaston KJ, Petchey OL. Including intraspecific variability in functional diversity. Ecol. 2009;90(1):81–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1890/07-1864.1.

    Article  CAS  Google Scholar 

  25. Cadotte MW, Dinnage R, Tilman D. Phylogenetic diversity promotes ecosystem stability. Ecol. 2012;93(sp8):S223–33. https://doiorg.publicaciones.saludcastillayleon.es/10.1890/11-0426.1.

    Article  Google Scholar 

  26. Tilman D. Functional diversity. In: Levin SA, editor. Encyclopedia of Biodiversity. San Diego: Academic Press (Elsevier Science); 2001. pp. 109–20.

    Chapter  Google Scholar 

  27. Rodrigues FHG, Medri IM, Tomas WM, Mourão GM. Revisão do conhecimento sobre ocorrência e distribuição de mamíferos do Pantanal. Corumbá: Embrapa Pantanal; 2002. p. 41.

    Google Scholar 

  28. Alho CJR. Biodiversity of the Pantanal: response to seasonal flooding regime and to environmental degradation. Braz J Biol. 2008;68(4):957–66. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S1519-69842008000500005.

    Article  CAS  PubMed  Google Scholar 

  29. Alho CJR, Fischer E, Oliveira-Pissini LF, Santos CF. Bat-species richness in the Pantanal floodplain and its surrounding uplands. Braz J Biol. 2011;71(1):311–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S1519-69842011000200010.

    Article  CAS  PubMed  Google Scholar 

  30. Cáceres NC, Carmignotto AP, Fischer E, Santos CF. Mammals from Mato Grosso do sul, Brazil. Check List. 2008;4(3):321–35. https://doiorg.publicaciones.saludcastillayleon.es/10.15560/4.3.321.

    Article  Google Scholar 

  31. Oliveira AKM, Conte CO, Oliveira-Pissini LF. Diversidade Da Quiropterofauna do Instituto De Pesquisas do Pantanal E Entorno, Aquidauana, Mato Grosso do sul, Brasil. Rev Bras Biocienc. 2011;9:96–102.

    Google Scholar 

  32. Oliveira SL, Souza LAS, Silva HK, Faria KC. Spatial configuration of the occurrence of bat species (Mammalia: Chiroptera) in eastern Mato Grosso, Brazil. Biota Neotrop. 2015;15:e20140122. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/1676-06032014012214.

    Article  Google Scholar 

  33. Fischer E, Silveira M, Munin RL, Camargo G, Santos CF, Pereira MJR, Fischer W, Eriksson A. Bats in the dry and wet Pantanal. Hystrix. 2018;29(1):11. https://doiorg.publicaciones.saludcastillayleon.es/10.4404/hystrix-00019-2017.

    Article  Google Scholar 

  34. Willig MR. Composition, microgeographic variation, and sexual dimorphism in Caatingas and Cerrado bat communities from Northeastern Brazil. Bull Carnegie Mus Nat Hist. 1983;23:1–131.

    Article  Google Scholar 

  35. Gonçalves E, Gregorin R. Quirópteros Da Estação Ecológica Da Serra Das Araras, Mato Grosso, Brasil, com o primeiro registro de Artibeus gnomus E A. Anderseni para o Cerrado. Lundiana: Int J Biodivers. 2004;5(1):143–9. https://doiorg.publicaciones.saludcastillayleon.es/10.35699/2675-5327.2004.22016.

    Article  Google Scholar 

  36. Bordignon M. Diversidade De Morcegos (Mammalia, Chiroptera) do Complexo Aporé-Sucuriú, Mato Grosso do sul, Brasil. Rev Bras Zool. 2006;23(4):1002–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S0101-81752006000400004.

    Article  Google Scholar 

  37. Camargo G, Fischer E, Gonçalves F, Fernandes G, Ferreira S. Morcegos do Parque Nacional Da Serra Da Bodoquena, Mato Grosso do sul, Brasil. Chiropt Neotrop. 2009;15(1):417–24.

    Google Scholar 

  38. Bezerra AMR, Marinho-Filho J. Bats of the Paranã River Valley, Tocantins and Goiás states, Central Brazil. Zootaxa. 2010;2725(1):41–56. https://doiorg.publicaciones.saludcastillayleon.es/10.11646/zootaxa.2725.1.3.

    Article  Google Scholar 

  39. Ferreira CMM, Fischer E, Pulchério-Leite A. Fauna De morcegos em remanescentes urbanos de Cerrado em Campo Grande, Mato Grosso do sul. Biota Neotrop. 2010;10(3):155–60. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S1676-06032010000300017.

    Article  Google Scholar 

  40. Zortéa M, Alho CJ. Bat diversity of a Cerrado habitat in central Brazil. Biodivers Conserv. 2008;17(4):791–805. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10531-008-9318-3.

    Article  Google Scholar 

  41. Zortéa M, Melo FR, Carvalho JC, Rocha ZD. Morcegos Da Bacia do rio Corumbá, Goiás. Chiropt Neotrop. 2010;16(1):611–7.

    Google Scholar 

  42. Sousa RFD, Venere PC, Faria KDC. Bats in forest remnants of the Cerrado savanna of eastern Mato Grosso, Brazil. Biota Neotrop. 2013;13:236–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S1676-06032013000200023.

    Article  Google Scholar 

  43. Garbino GST, Gregorin R, Lima IP, Loureiro L, Moras LM, Moratelli R, Nogueira MR, Pavan AC, Tavares VC, do Nascimento MC, Peracchi AL. Updated checklist of Brazilian bats: versão 2020. Comitê da Lista de Morcegos do Brasil—CLMB. Sociedade Brasileira para o Estudo de Quirópteros (Sbeq). 2020. https://www.sbeq.net/lista-de-especies. Accessed 7 Feb 2022.

  44. Ponce VM, Cunha CN. Vegetated earthmounds in Tropical savannas of Central Brazil: a synthesis: with special reference to the Pantanal do Mato Grosso. J Biogeogr. 1993;20:219–25. https://doiorg.publicaciones.saludcastillayleon.es/10.2307/2845673.

    Article  Google Scholar 

  45. Harris MB, Tomas W, Mourão G, Silva CJ, Guimarães E, Sonoda F, Fachim E. Safeguarding the Pantanal wetlands: threats and conservation initiatives. Conserv Biol. 2005;19(3):714–20. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1523-1739.2005.00708.x.

    Article  Google Scholar 

  46. Klink CA, Machado RB. Conservation of the Brazilian cerrado. Conserv Biol. 2005;19(3):707–13. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1523-1739.2005.00702.x.

    Article  Google Scholar 

  47. Strassburg BB, Brooks T, Feltran-Barbieri R, Iribarrem A, Crouzeilles R, Loyola R, Latawiec AE, Oliveira Filho FJB, Scaramuzza CAM, Scarano FR, Soares-Filho B, Balmford A. Moment of truth for the Cerrado hotspot. Nat Ecol Evol. 2017;1(4):1–3. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41559-017-0099.

    Article  Google Scholar 

  48. Laliberté E, Wells JA, DeClerck F, Metcalfe DJ, Catterall CP, Queiroz C, et al. Land-use intensification reduces functional redundancy and response diversity in plant communities. Ecol Lett. 2010;13(1):76–86. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1461-0248.2009.01403.x.

    Article  PubMed  Google Scholar 

  49. Flynn DF, Gogol-Prokurat M, Nogeire T, Molinari N, Richers BT, Lin BB, Simpson N, Mayfield MM, DeClerck F. Loss of functional diversity under land use intensification across multiple taxa. Ecol Lett. 2009;12(1):22–33. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1461-0248.2008.01255.x.

    Article  PubMed  Google Scholar 

  50. Mayfield MM, Bonser SP, Morgan JW, Aubin I, McNamara S, Vesk PA. What does species richness tell us about functional trait diversity? Predictions and evidence for responses of species and functional trait diversity to land-use change. Glob Ecol Biogeogr. 2010;19(4):423–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1466-8238.2010.00532.x.

    Article  Google Scholar 

  51. Carvalho WD, Meyer CF, Xavier BDS, Mustin K, Castro IJD, Silvestre SM, Pathek DB, Capaverde UD Jr, Hilário R, Toledo JJ. Consequences of replacing native savannahs with acacia plantations for the taxonomic, functional, and phylogenetic α-and β-diversity of bats in the Northern Brazilian Amazon. Front Ecol Evol. 2020;8:609214.

    Article  Google Scholar 

  52. Pereira MJR, Fonseca C, Aguiar LM. Loss of multiple dimensions of bat diversity under land-use intensification in the Brazilian Cerrado. Hystrix. 2018;29(1):25. https://doiorg.publicaciones.saludcastillayleon.es/10.4404/hystrix-00020-2017.

    Article  Google Scholar 

  53. Cordeiro JLP. Estrutura e heterogeneidade da paisagem de uma unidade de conservação no nordeste do Pantanal (RPPN SESC Pantanal), Mato Grosso, Brasil: efeitos sobre a distribuição e densidade de antas (Tapirus terrestris) e de cervos-do-pantanal (Blastocerus dichotomus). Doctoral thesis, Universidade Federal do Rio Grande do Sul; 2004.

  54. Oliveira LFB, Cordeiro JLP, Hasenack H. Padrões E tendências espaço-temporais na estrutura de uma paisagem antropizada no norte do Pantanal. In: Peres CA, Barlow J, Gardner TA, Vieira ICG, editors. Conservação Da Biodiversidade em paisagens antropizadas do Brasil. Curitiba: Editora UFPR; 2013. pp. 231–52.

    Google Scholar 

  55. Louzada NSV, do Monte Lima AC, Pessôa LM, Cordeiro JLP, Oliveira LFB. New records of phyllostomid bats for the state of Mato Grosso and for the cerrado of Midwestern Brazil (Mammalia: Chiroptera). Check List. 2015;11(3):1644–1644. https://doiorg.publicaciones.saludcastillayleon.es/10.15560/11.3.1644.

    Article  Google Scholar 

  56. Pennington RT, Ratter JA. Neotropical savannas and seasonally dry forests: plant diversity, biogeography, and conservation. Boca Raton: CRC; 2006. p. 508.

    Book  Google Scholar 

  57. Ribeiro JF, Walter BMT. As principais fitofisionomias do Bioma Cerrado. In: Sano SM, Almeida SP, Ribeiro JF, editors. Cerrado: Ecologia E Flora. Brasília: Embrapa Cerrados; 2008. pp. 151–212.

    Google Scholar 

  58. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’hara RB et al. Package ‘vegan’. Community Ecology Package, version 2.6-4. 2018. Chrome-extension://efaidnbmnnnibpcajpcglclefindmkaj/ https://cran.ism.ac.jp/web/packages/vegan/vegan.pdf. Accessed Sept 2022.

  59. Chao A, Gotelli NJ, Hsieh TC, Sander EL, Ma KH, Colwell RK, Ellison AM. Rarefaction and extrapolation with Hill numbers: a framework for sampling and estimation in species diversity studies. Ecol Monogr. 2014;84(1):45–67. https://doiorg.publicaciones.saludcastillayleon.es/10.1890/13-0133.1.

    Article  Google Scholar 

  60. Hsieh TC, Ma KH, Chao A, Hsieh MT. Package ‘iNEXT’. 2016. http://chao.stat.nthu.edutw/blog/software-download. Accessed Feb 2022.

  61. Tamura K, Stecher G, Kumar S. MEGA 11: Molecular Evolutionary Genetics Analysis Version 11. Mol Biol Evol. 2021;38(7):3022–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/molbev/msab120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Trifinopoulos J, Nguyen LT, von Haeseler A, Minh BQ. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016;44(W1):W232–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/nar/gkw256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Kalyaanamoorthy S, Minh BQ, Wong TK, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/nmeth.4285.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21. https://doiorg.publicaciones.saludcastillayleon.es/10.1093/sysbio/syq010.

    Article  CAS  PubMed  Google Scholar 

  65. Rambaut A, Drummond A. FigTree: tree figure drawing tool, version v1.4.4. Edinburgh: Institute of Evolutionary Biology, University of Edinburgh. 2018. http://tree.bio.ed.ac.uk/software/figtree/. Accessed May 2023.

  66. de Oliveira MB, Povill C, de Oliveira LFB, Bonvicino CR. Phylogeography and potential distribution of Sturnira lilium and S. giannae (Chiroptera: Stenodermatinae) with range extension for S. Giannae in the Cerrado and Pantanal biomes. Zool. 2024;41:e23061. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S1984-4689.v41.e23061.

    Article  Google Scholar 

  67. Faith DP. Conservation evaluation and phylogenetic diversity. Biol Conserv. 1992;61(1):1–10. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/0006-3207(92)91201-3.

    Article  Google Scholar 

  68. Webb CO. Exploring the phylogenetic structure of ecological communities: an example for rain forest trees. Am Nat. 2000;156(2):145–55. https://doiorg.publicaciones.saludcastillayleon.es/10.1086/303378.

    Article  PubMed  Google Scholar 

  69. Webb CO, Ackerly DD, McPeek MA, Donoghue MJ. Phylogenies and community ecology. Annu Rev Ecol Syst. 2002;33(1):475–505. https://doiorg.publicaciones.saludcastillayleon.es/10.1146/annurev.ecolsys.33.010802.150448.

    Article  Google Scholar 

  70. Blomberg SP, Garland T Jr, Ives AR. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evol. 2003;57(4):717–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.0014-3820.2003.tb00285.x.

    Article  Google Scholar 

  71. Wilman H, Belmaker J, Simpson J, de la Rosa C, Rivadeneira MM, Jetz W. EltonTraits 1.0: species-level foraging attributes of the world’s birds and mammals: ecological archives E095‐178. Ecol. 2014;95(7):2027–2027. https://doiorg.publicaciones.saludcastillayleon.es/10.1890/13-1917.1.

    Article  Google Scholar 

  72. Findley JS. Bats: a community perspective. New York: Cambridge University Press; 1993. p. 167.

    Google Scholar 

  73. Bonato V, Facure KG, Uieda W. Food habits of bats of subfamily Vampyrinae in Brazil. J Mammal. 2004;85(4):708–13. https://doiorg.publicaciones.saludcastillayleon.es/10.1644/BWG-121.

    Article  Google Scholar 

  74. Andrade TY, Thies W, Rogeri PK, Kalko EK, Mello MA. Hierarchical fruit selection by neotropical leaf-nosed bats (Chiroptera: Phyllostomidae). J Mammal. 2013;94(5):1094–101. https://doiorg.publicaciones.saludcastillayleon.es/10.1644/12-MAMM-A-244.1.

    Article  Google Scholar 

  75. Fenton MB, Whitaker JO Jr, Vonhof MJ, Waterman JM, Pedro WA, Aguiar L, Baumgarten JE, Bouchard S, Faria DM, Portfors CV, Rautenbach NIL, Scully W, Zortéa M. The diet of bats from Southeastern Brazil: the relation to echolocation and foraging behaviour. Rev Bras Zool. 1999;16(4):1081–5. https://doiorg.publicaciones.saludcastillayleon.es/10.1590/S0101-81751999000400017.

    Article  Google Scholar 

  76. Soriano PJ. Functional structure of bat communities in tropical rainforests and andean cloud forests. Ecotropicos. 2000;13(1):1–20.

    Google Scholar 

  77. Denzinger A, Schnitzler HU. Bat guilds, a concept to classify the highly diverse foraging and echolocation behaviors of microchiropteran bats. Front Physiol. 2013;4:164. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fphys.2013.00164.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Goslee SC, Urban DL. The ecodist package for dissimilarity-based analysis of ecological data. J Stat Softw. 2007;22:1–19. https://doiorg.publicaciones.saludcastillayleon.es/10.18637/jss.v022.i07.

    Article  Google Scholar 

  79. Petchey OL, Gaston KJ. Functional diversity (FD), species richness, and community composition. Ecol Lett. 2002;5(3):402–11. https://doiorg.publicaciones.saludcastillayleon.es/10.1046/j.1461-0248.2002.00339.x.

    Article  Google Scholar 

  80. Cornwell WK, Schwilk DW, Ackerly DD. A trait-based test for habitat filtering: convex hull volume. Ecol. 2006;87(6):1465–71. https://doiorg.publicaciones.saludcastillayleon.es/10.1890/0012-9658(2006)875b1465:ATTFHF5d2.0.CO;2.

    Article  Google Scholar 

  81. Villéger S, Mason NWH, Mouillot D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecol. 2008;89(8):2290–301. https://doiorg.publicaciones.saludcastillayleon.es/10.1890/07-1206.1.

    Article  Google Scholar 

  82. Fontana S, Petchey OL, Pomati F. Individual-level trait diversity concepts and indices to comprehensively describe community change in multidimensional trait space. Funct Ecol. 2016;30(5):808–18. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/1365-2435.12551.

    Article  Google Scholar 

  83. Pla L, Casanoves F, Di Rienzo J. Quantifying functional biodiversity. Berlim: Springer Science & Business Media; 2012. p. 98.

    Book  Google Scholar 

  84. Laliberté E, Legendre P, Shipley B, Laliberté ME. Package ‘fd’. Measuring functional diversity from multiple traits, and other tools for functional ecology, version 1.0–12. 2014. chrome-extension://efaidnbmnnnibpcajpcglclefindmkaj/https://citeseerx.ist.psu.edu/documentrepid=rep1type=pdfdoi=c0123be0f0737359fa9303aece03fd500d2499b1. Accessed Sept 2022.

  85. Sievert C, Parmer C, Hocking T, Chamberlain S, Ram K, Corvellec M, Despouy P, Brüggemann S. Package ‘plotly’. R Foundation for Statistical Computing, Vienna. 2021. https://cran.r-project.org/package=plotly. Accessed Jan 2022.

  86. R Core Team. R: A language and environment for statistical computing. Austria: R Foundation for Statistical Computing. 2023. https://www.R-project.org/. Accessed Oct 2023.

  87. Farneda FZ, Rocha R, Aninta SG, López-Baucells A, Sampaio EM, Palmeirim JM, Bobrowiec PED, Dambros CS, Meyer CF. Bat phylogenetic responses to regenerating amazonian forests. J Appl Ecol. 2022;59(8):1986–96. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/1365-2664.14041.

    Article  Google Scholar 

  88. Bogoni JA, Carvalho-Rocha V, da Silva PG. Spatial and land‐use determinants of bat species richness, functional diversity, and site uniqueness throughout the largest tropical country, Brazil. Mammal Rev. 2022;52(2):267–83. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/mam.12279.

    Article  Google Scholar 

  89. Silva DC, Oliveira HFM, Domingos FMCB. Cerrado bat community assembly is determined by both present-day and historical factors. J Biogeogr. 2024;00:1–16. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/jbi.14798.

    Article  Google Scholar 

  90. Newbold T, Hudson LN, Arnell AP, Contu S, de Palma A, Ferrier Hill SLL, Hoskins AJ, Lysenko I, Phillips HRP, Burton VJ, Chng CWT, Emerson S, Gao D, Pask-Hale G, Hutton J, Jung M, Sanchez-Ortiz K, Simmons BI, Whitmee S, Zhang H, Scharlemann JPW, Purvis A. Has land use pushed terrestrial biodiversity beyond the planetary boundary? A global assessment. Science. 2016;353(6296):288–91. https://doiorg.publicaciones.saludcastillayleon.es/10.1126/science.aaf2201.

    Article  CAS  PubMed  Google Scholar 

  91. Chisholm RA, Lim F, Yeoh YS, Seah WW, Condit R, Rosindell J. Species–area relationships and biodiversity loss in fragmented landscapes. Ecol Lett. 2018;21(6):804–13. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/ele.12943.

    Article  PubMed  PubMed Central  Google Scholar 

  92. Cadotte MW, Tucker CM. Should environmental filtering be abandoned? Trends Ecol Evol. 2017;32(6):429–37. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.tree.2017.03.004.

    Article  PubMed  Google Scholar 

  93. Ausprey IJ, Newell FL, Robinson SK. Functional response traits and altered ecological niches drive the disassembly of cloud forest bird communities in tropical montane countrysides. J Anim Ecol. 2022;91(11):2314–28. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/1365-2656.13816.

    Article  PubMed  Google Scholar 

  94. Terborgh J, Pitman N, Silman M, Schichter H, Núñez VP. Maintenance of tree diversity in tropical forests. In: Levey DJ, Silva WR, Galetti M, editors. Seed dispersal and frugivory: ecology, evolution and conservation. Third International Symposium-Workshop on Frugivores and Seed Dispersal, São Pedro, Brazil. Wallingford UK: CABI Publishing; 2002. pp. 1–17.

  95. Kunz TH, Braun de Torrez E, Bauer D, Lobova T, Fleming TH. Ecosystem services provided by bats. Ann N Y Acad Sci. 2011;1223(1):1–38. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1749-6632.2011.06004.x.

    Article  PubMed  Google Scholar 

  96. Kelm DH, Wiesner KR, Helversen OV. Effects of artificial roosts for frugivorous bats on seed dispersal in a Neotropical forest pasture mosaic. Conserv Biol. 2008;22(3):733–41. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1523-1739.2008.00925.x.

    Article  PubMed  Google Scholar 

  97. Safi K, Meiri S, Jones KE, Smith FA, Lyons K. Evolution of body size in bats. In: Smith FA, Lyons SK, editors. Animal body size: linking pattern and process across Space, Time, and Taxonomic Group. Chicago and London: University of Chicago Press; 2013. pp. 95–115.

    Chapter  Google Scholar 

  98. Schondube JE, Herrera-M L, del Rio CM. Diet and the evolution of digestion and renal function in phyllostomid bats. Zool. 2001;104(1):59–73. https://doiorg.publicaciones.saludcastillayleon.es/10.1078/0944-2006-00007.

    Article  CAS  Google Scholar 

  99. Srivastava DS, Cadotte MW, MacDonald AAM, Marushia RG, Mirotchnick N. Phylogenetic diversity and the functioning of ecosystems. Ecol Lett. 2012;15(7):637–48. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1461-0248.2012.01795.x.

    Article  PubMed  Google Scholar 

  100. Cisneros LM, Fagan ME, Willig MR. Environmental and spatial drivers of taxonomic, functional, and phylogenetic characteristics of bat communities in human-modified landscapes. PeerJ. 2016;4:e2551. https://doiorg.publicaciones.saludcastillayleon.es/10.7717/peerj.2551.

    Article  PubMed  PubMed Central  Google Scholar 

  101. Frank HK, Frishkoff LO, Mendenhall CD, Daily GC, Hadly EA. Phylogeny, traits, and biodiversity of a neotropical bat assemblage: close relatives show similar responses to local deforestation. Am Nat. 2017;190(2):200–12. https://doiorg.publicaciones.saludcastillayleon.es/10.1086/692534.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Martha L Brandão, José LP Cordeiro, Fabricio Escarlate-Tavares, and Fernanda Trierveiler for their efforts in the field, and Gustavo Staut for his support and permission to work at the Santa Lucia ranch, Brazilian Pantanal. We thank Sesc Pantanal for their support of the project and assistance with field sampling. The project received partial support from the Conselho Nacional de Desenvolvimento Científico (CNPq, process 400713/2013-6) for Luiz FB de Oliveira. Laboratory analyses were funded by grants awarded to Cibele R Bonvicino (CNPq, process 311712/2021-5) and by the Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ, E-26/210.925/2019). This research was supported by a PhD scholarship awarded to Marcione B. de Oliveira from the CNPq (process 162273/2018-5). We thank the two anonymous reviewers for their valuable comments on the manuscript.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization (MBO, LFBO, CRB); data curation and analysis (MBO); funding acquisition (LFBO, CRB); investigation (MBO, LFBO, CRB); writing (MBO) – original draft; writing – review and editing (MBO, LFBO, CRB).

Corresponding author

Correspondence to Marcione B. de Oliveira.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

40693_2025_137_MOESM1_ESM.docx

Supplementary Material 1: Table S1. List of specimens used in phylogenetic diversity analyses, with taxon, GenBank accession number or voucher number, and font. *Due to the lack of data on Promops nasutus in the present study, a congener, Promops centralis, was utilized. Table S2. Standardized PCR amplification cycles for the cytochrome b gene of the analyzed chiropteran genera. Some genera had two standardized cycles. Table S3. Primer list employed in the amplification (A) and sequencing (S) steps of the cytochrome b gene. Table S4. Compilation of functional traits for the sampled bats. Mean forearm length (mm) and body mass (g) were calculated based only on adult individuals of species with more than ten captures* in the present study. For species with few individuals, data from the literature were sourced from Velazco et al. (2010), Moratelli et al. (2013), Reis et al. (2017), and Velazco and Patterson (2019). A species-level classification based on diet and foraging behavior was employed based on Findley (1993), Fenton (1999), Soriano (2000), Schnitzler and Kalko (2001), and Schnitzler et al. (2003).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

de Oliveira, M.B., de Oliveira, L.F.B. & Bonvicino, C.R. Analysis of functional and phylogenetic diversity of bat communities: insights from resource use in the Pantanal and Cerrado. Rev. Chil. de Hist. Nat. 98, 2 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40693-025-00137-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40693-025-00137-z

Keywords