1. Introduction
A large part of terrestrial ecosystems is impacted due to the conversion of natural habitats into anthropic habitats [
1]. The conversion of natural habitats into anthropogenic ones changes the composition (e.g., types and proportions of land uses) and configuration (e.g., spatial arrangement of landscape elements) of the landscape, causing the loss of native habitats and fragmentation of the remnants [
2]. Thus, habitat loss and fragmentation are the main causes of defaunation, including vertebrates of larger body size such as medium and large-sized mammals [
3,
4,
5,
6].
Habitat amount is commonly evaluated to explain the species richness in forest fragments [
7,
8,
9]. However, the matrix around the forest remnants also influences the occurrence of species and ecological dynamics [
10,
11,
12,
13]. The effects that the matrix can have on species are mainly related to species movement and dispersion, availability of resources and changes in abiotic characteristics. Therefore, the magnitude of the impact that the types and quantities of matrices in the landscape exert on populations depends on characteristics that influence the life history of the species (e.g., ability to move, diet, habitat requirement) [
14]. In this way, anthropogenic matrices can have ambivalent effects on the impacts of habitat loss and habitat fragmentation on biodiversity.
Habitat fragmentation per se changes the original landscape configuration, increasing the number of forest fragments of different sizes and shapes, dispersed in the landscape and under greater edge effect [
15,
16]. The effects of fragmentation on biodiversity have been discussed [
12,
16,
17,
18,
19] and most of the effects are considered as negative [
18,
20,
21], but also positive or neutral [
16,
17]. However, many studies evaluated the effect of habitat fragmentation without distinguishing it from the effect of habitat loss and without considering an appropriate scale [
17]. In this way, fragmentation must be accessed from metrics related to landscape configuration since it is characterized by the division of habitat, independent of habitat loss [
15,
17].
One of the challenges of studies at the landscape level is to define the best scale at which the landscape composition and configuration variables should be collected [
22,
23,
24]. This is not a trivial question, as species’ responses to landscape variables are shaped by ecological processes - such as metapopulation dynamics, dispersion, inter and intra-specific interactions - that act at different scales [
16,
24,
25,
26]. A solution to this challenge is to define the scale of effect of each landscape variable, that is, the scale in which the landscape variable has the greatest effect under the evaluated response variable and consider the multiscale approach, which each variable is evaluated in their respective effect scale [
22,
24].
The Atlantic Forest and the Cerrado are Brazilian domains with high species endemism and are highly threatened due to, mainly, conversion of habitat to agriculture, livestock and urbanization and, therefore, they are considered biodiversity hotspots [
27,
28,
29,
30]. The forest cover of the contemporary Atlantic Forest is restricted to only 12% of its original cover, with approximately 84% of the remainder represented by forest fragments smaller than 50 hectares [
31]. While approximately 50% of the native vegetation cover in the Cerrado has been transformed into agriculture and pasture [
30]. Consequently, the biological populations present in the forest remnants are under the influence of the landscape characteristics to persist in these fragmented landscapes.
Our main goal was to compare how landscape configuration and composition affect medium and large-sized mammal communities in agricultural ecosystems. We also evaluated the scale of effect of each landscape variable for both mammal richness and composition, and we compared the influence of matrix types on the mammal communities. We hypothesized that landscape composition variables affect more the mammal communities than the landscape configuration ones. Also, we expected that the scales of effect vary among variables and the more different the matrix to the native habitat, the greater the influence of the matrix on the mammal communities.
4. Discussion
The scales of effect varied among the landscape and response variables, indicating that the landscape variables affect large and medium-sized mammal species in different aspects, highlighting the importance of evaluating different measures to understand the effect of landscape structure on biodiversity. We found that the landscape composition (forest, pasture and coffee covers) influences more the richness and composition of medium and large-sized mammal than the landscape configuration (number of forest fragments and mean distance among forest fragments) and the pasture in the landscape has a positive influence on species richness while the forest cover influences species composition.
The landscape composition and configuration influence ecological processes and can have different effects depending on the community structure and the spatial scale [
25,
26]. Different variables affect different dimensions of the ecological requirements of organisms leading to scale-dependent responses. As, for example, ecological processes such as metapopulation dynamics can shape ecological responses on larger spatial scales than processes such as competition, predation, resource availability, and anthropogenic impacts such as hunting [
43,
44]. In addition, different species can have different responses from the same variable at the same scale. That is, the quantity of a given agricultural matrix can be a barrier for the movement of species from the interior of the forest [
10,
45,
46]. However, for generalist species, some agriculture matrix can function as an alternative foraging site, or even displacement corridor [
47,
48]. Thus, the scale of effect resulting from the response of species communities is the combination of the variation in the influence that each independent variable has on each species.
In our study, the pasture cover was the only variable that explained the variation in species richness, with a high explanation support (33%) for a positive relationship between the pasture cover and species richness. Although this result is contrary to the expected positive relation between richness and native habitat cover in agricultural landscapes [
7,
49,
50], agroecosystems can favor the presence of generalist species from open areas and/or typical of the Cerrado. Studies have shown that landscapes altered by anthropic processes tend to favor the occurrence of species with low dependence on intact ecosystems or forest interior [
51,
52]. More than 50% of the species registered in our study have wide geographical distribution, varied diet and high movement capacity [
53]. That is, species that can use different types of habitats and that can be favored in modified landscapes [
54]. In addition, some studies have also found no relationship between the richness of medium and large-sized mammals and the amount of native habitat [
8,
21,
55,
56] or others have recorded a high species richness in landscapes with strong urban [
58,
59] or agricultural influence [
8,
49,
50]. In this way, species with a wide ecological niche can be favored by these modified landscapes and, probably, maintain their populations in agricultural landscapes [
60].
On the other hand, the positive relationship between species richness and the pasture cover in the landscape may be due to the use of the almost exclusive remaining fragments by the species, since the proportion of pasture is inversely proportional to that of forest. In addition, we found, for example, that the landscape with the highest number of species (landscape P (e), Supplementary Table 4) has a high pasture cover and a low frequency of individuals of the sampled species. Consequently, the sampled fragment can be used by species present in the landscape as a foraging site, reproduction or as stepping stones, mainly due to the reduced forest cover in the landscape [
61,
62]. In general, species with large home range, as is the case of most medium and large-sized mammals sampled, which persist in modified landscapes depend on the remaining fragments present in the landscape [
63], highlighting the importance of keeping small fragments in agricultural landscapes [
17].
Our results also show that despite forest cover does not explain the variation in species richness, it influences species composition, explaining part of the variation between the sampled communities. Some authors [
64,
65] developed and improved the metapopulation theory based on the analogy to island biogeography theory [
66], where the fragment size and isolation in a terrestrial landscape correspond to the size and isolation of oceanic islands. Since then, studies on habitat fragmentation include variables to represent this analogy and explain the relationship between species-area in continental environments [
67,
68]. However, the matrices around the fragments affect the biodiversity differently than the ocean, since fragmented landscapes are not binary like oceanic islands [
45,
69]. Nevertheless, Fahrig [
16] argued that both effects of fragment size and isolation are shaped by the effect of the sampling area. That is, in a landscape of a certain size, the amount of forest habitat around the sampled forest fragment would be a better predictor for species richness [
16] and other ecological responses [
9,
70] than fragment size and isolation [
16]. In addition, Fahrig [
16] adds that the matrix effect is secondary to the amount of habitat and that the landscape configuration should have little or no effect on the species. Although we have not corroborated this hypothesis for species richness, we found that the amount of native habitat (i.e., forest cover) was the variable with the greatest effect on the mammal composition, followed, and with minor importance, by the coffee cover and finally by mean distance among forest fragments.
As landscape composition has a greater effect on medium and large-sized mammals than landscape configuration, more heterogeneous and structurally more complex agricultural landscapes can favor greater connectivity between the fragments and increase the diversity of species [
71]. In addition, the role and management of the matrix in order to contribute to the conservation of biodiversity is one of the main strategies for the permanence of species in modified landscapes [
14]. However, matrix management involves economic and social aspects in addition to ecological ones [
72] and, therefore, an interdisciplinary approach is necessary for conservation measures in agroecosystems to be applied and effective.
Finally, exploring different responses to biodiversity can be important to understand different aspects of landscape influence on species. As well, the scale in which variables of landscape composition and configuration are measured can interfere with the species-landscape relationship. Therefore, we emphasize the importance of future studies on the relationship between landscape structure and biodiversity to consider different biodiversity responses and a multiscale approach with a wide range of scales. Also, with the advancement of geoprocessing tools and development of free software, such studies are currently more viable. Considering landscape composition and configuration variables on an appropriate spatial scale allows us to expand our knowledge about the relationship between landscape and biodiversity and thus improve management and conservation measures in agricultural ecosystems.