Next Article in Journal
Identification and Characterization of 24-Dehydrocholesterol Reductase (DHCR24) in the Two-Spotted Cricket, Gryllus bimaculatus
Next Article in Special Issue
What Are the Most Important Factors Influencing Springtail Tetrodontophora bielanensis?
Previous Article in Journal
Insect Epigenetic Mechanisms Facing Anthropogenic-Derived Contamination, an Overview
Previous Article in Special Issue
Community Structure, Biodiversity and Spatiotemporal Distribution of the Black Flies (Diptera: Simuliidae) Using Malaise Traps on the Highest Mountain in Thailand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity Patterns of Dung Beetles along a Mediterranean Elevational Gradient

1
Department of Life, Health and Environmental Sciences, University of L’Aquila, Via Vetoio, 67100 L’Aquila, Italy
2
cE3c-Centre for Ecology, Evolution and Environmental Changes, Azorean Biodiversity Group, Faculdade de Ciências e Engenharia do Ambiente, Universidade dos Açores, rua Capitão João D’Avila, 9700-042 Angra do Heroísmo, Portugal
*
Author to whom correspondence should be addressed.
Insects 2021, 12(9), 781; https://doi.org/10.3390/insects12090781
Submission received: 30 July 2021 / Revised: 23 August 2021 / Accepted: 28 August 2021 / Published: 31 August 2021
(This article belongs to the Collection Insects in Mountain Ecosystems)

Abstract

:

Simple Summary

Mountains are ideal natural laboratories to study how biodiversity is influenced by environmental characteristics because climate varies rapidly from lowlands to high elevations. Scientists have investigated how the number of species varies with elevation for the most disparate plant and animal groups worldwide. However, species richness is only one aspect of biodiversity that does not consider species abundances. The so-called Hill numbers are a unified family of mathematical indices that express biodiversity in terms of both richness and abundance. We used Hill numbers to investigate how dung beetle diversity varies along an elevational gradient in a Mediterranean mountain. We found that scarabaeids were the most abundant dung beetle group. These insects construct subterranean nests protecting their offspring from desiccation during warm and dry summer climatic conditions. Additionally, in accordance with their preference for open habitats, we found that dung beetles are more abundant and diversified in grasslands than in woodlands. In the woodlands, diversity increased with elevation because of tree thinning, whereas in the grasslands, diversity decreased because of increasingly harsher environmental conditions. This indicates a trade-off in the beetle response to elevation between the positive effects of increasing the availability of suitable habitats and the worsening of environmental conditions.

Abstract

Most studies of biodiversity–elevational patterns do not take species abundance into consideration. Hill numbers are a unified family of indices that use species abundance and allow a complete characterization of species assemblages through diversity profiles. Studies on dung beetle responses to elevation were essentially based on species richness and produced inconsistent results because of the non-distinction between different habitats and the use of gradients dispersed over wide areas. We analyzed dung beetle diversity in a Mediterranean mountain (central Italy) for different habitats (woodlands vs. grasslands) and taxonomic groups (scarabaeids and aphodiids). Scarabaeids were the most abundant. Since scarabaeids are able to construct subterranean nests, this indicates that the warm and dry summer climatic conditions of high elevations favor species capable of protecting their larvae from desiccation. Dung beetles were more abundant and diversified in grasslands than in woodlands, which is consistent with their preference for open habitats. In the woodlands, diversity increased with increasing elevation because of increasing tree thinning, whereas, in the grasslands, diversity decreased with elevation because of increasingly harsher environmental conditions. These results indicate a trade-off in the beetle response to elevation between the positive effects of increasing the availability of more suitable habitats and the decrease of optimal environmental conditions.

1. Introduction

Patterns of variation in biodiversity along elevational gradients have long attracted the interest of naturalists, and they have become increasingly popular in recent years [1,2]. Alpha-diversity (i.e., the diversity of species within a particular area or ecosystem), beta-diversity (i.e., species turnover), functional diversity, animal body size, biogeographical composition, and several other aspects of community structure vary markedly with altitude in response to the associated variation in environmental parameters [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]. For example, climatic conditions change more rapidly with altitude than with latitude (the world average is a drop of about 1 °C for every 150 m above sea level against 1 °C for every 150 km poleward), which makes mountain areas ideal natural laboratories to investigate patterns of variation in biodiversity in response to environmental factors within geographically restricted areas [6,18].
Most of the studies investigating elevational patterns in biodiversity have dealt with species richness (i.e., number of species found at different altitudes) as a measure of alpha-diversity, finding either a monotonic decrease in species richness with increasing elevation or a hump-shaped trend with a mid-elevational peak [13,17,19]. The monotonic decrease is generally interpreted as a consequence of increasing harsher conditions (such as lower temperatures, higher radiation, low partial pressures of both oxygen and carbon dioxide, stronger winds, lower soil nutrients, less stable substrates, and shorter plant growing seasons), decreasing productivity, decreasing available area (because of the conical shape of mountains), and nested species distributions, in which higher assemblages tend to be progressively smaller subsamples comprising the most tolerant species among those already present in the lower ones [13,18]. By contrast, a hump-shaped pattern is frequently interpreted as a result of purely stochastic processes (species ranges tend to overlap in domain centers due to dispersal limitations, resulting in a mid-domain richness peak) [13,18].
Species richness is, however, a problematic measure of the diversity of an assemblage because it depends on sampling intensity and does not take species abundance distribution into consideration [20]. The first problem (i.e., the effect of sampling effort) can be addressed by using rarefaction/extrapolation procedures to obtain richness values that are standardized to a common sampling effort or by calculating asymptotic estimators that are relatively independent of additional sampling efforts [20]. The second problem (i.e., the abundance problem) has led many ecologists to propose measures (indexes) of diversity that combine species richness and the proportion of each species into a single metric [21,22].
Among the vast multitude of diversity measures proposed to take the abundance of the species in an assemblage into account, Hill numbers qD are a mathematically unified family of indices that (1) combine information on species richness, species rarity, and species dominance and (2) are all expressed in the same units (i.e., effective number of species); therefore, they are comparable between each other [20,23]. As the diversity order q can assume any value, it is possible to obtain a complete characterization of an assemblage by constructing a diversity profile plotting qD versus q. For these reasons, Hill numbers are increasingly used in community ecology research [20,24].
In this paper, we used Hill numbers to investigate the elevational patterns of diversity of dung beetles in a Mediterranean montane area. Dung beetles are one of the most popular beetle groups among insect ecologists [25] because of their taxonomic richness (with 7000 known species worldwide [25]), the many essential ecosystems services they provide (such as nutritional cycling [26,27], their soil maintenance characteristics [28,29,30,31,32,33,34], seed dispersal [35], and their ability to protect livestock from dipterans that develop in dung [36]), their responses to environmental changes [37,38,39], and their sensitivity to anthropogenic disturbances [40,41,42,43]. Therefore, it is not surprising that dung beetle assemblages have been the object of various research studies dealing with elevational patterns of diversity.
Most studies on dung beetle diversity in the mountainous areas in Europe, North America, Southeast Asia, South Africa, and South America have highlighted a decrease in the number of species with increasing altitude (e.g., [37,44,45,46,47,48,49,50,51,52,53,54]). In the Western Palearctic, dung beetle richness was found to decline with increasing elevation in the Northern French Alps [44,45], in the Southern French Alps [46], in the Sierra Nevada (Spain) [46,47], in Mediterranean grassland sites in the Southern French Alps [49], and in the Kütahya Province in western Turkey [53]. Hump-shaped patterns were, however, found in the northern French Alps (with the mode within 1500–1800 m along a gradient from about 700 to about 3000 m) [46,48], in the Western Rhodopes Mountains (with the mode at 900 m along a gradient from 200 to 2000 m) [55], and in the Iberian Central System (with the mode at 1000 m along a gradient from 500 to 2200) [55]. No effect of altitude on species richness was observed in a mountainous area near the Cantabrian Sea (namely the western portion of the Picos de Europa, Spain) [56] and in Sierra de Baza (south-eastern Iberian Peninsula) [57].
These studies, however, were characterized by various sources of complexity that make comparisons and interpretations difficult. First, with the exception of a single study [56], all of the available research (e.g., [44,45,53,55,57,58]) used elevational gradients extended over wide areas, sometimes including regions with different ecological settings and biogeographical histories. When gradients are constructed by taking localities dispersed over wide areas, the communities sampled at different elevations might belong to faunas influenced by different biogeographical and evolutionary processes [58].
Second, many studies on the elevational patterns of dung beetle diversity did not consider the differences in the type of habitat with which dung beetles are associated (e.g., [44,45,46,55]), while others investigated a single habitat [49] or compared the response between different habitats [53,56]. Most European dung beetles are generally associated with open habitats, and it has been reported that as vegetation cover increases, species richness and total abundance decrease [59,60,61,62]. For example, Romero-Alcaraz and Ávila [57] found that Mediterranean dung beetle communities are characterized by higher values of abundance and diversity in open habitats (e.g., open scrublands with little vertical development and open pastures) than in closed (forest) habitats. Moreover, responses to habitat type may change with elevation [62]. Even the local climatic setting, which is highlighted by the presence of different vegetation types, influences diversity patterns. For example, in the mountain grasslands of the Southern Alps, Errouissi et al. [49] found that dung beetle diversity (expressed using the Margalef index, which balances the richness by the number of sampled individuals) decreased with increasing elevation for Mediterranean sites but increased in temperate sites. Thus, results from different studies may appear contradictory if they were conducted in different habitats, or are difficult to interpret if data from different habitat types were merged.
Third, available studies were mostly conducted using species richness as a measure of diversity, with only cursory notes on other aspects of diversity (e.g., negative influences of elevation on the Shannon index and the Shannon-based equitability [44,45] as well as on the Margalef index [49] for Alpine dung beetle assemblages; negative influence on the Shannon index in an Afromontane forest [54]; no influence of elevation on the Simpson index in a tropical area in Brazil [52]). To the best of our knowledge, no research used diversity profiles.
In this paper, we applied a diversity profile analysis based on Hill numbers to dung beetle assemblages sampled along an elevational gradient at a narrow geographical scale to overcome the problem of mixing communities belonging to different faunas. Additionally, we conducted separate analyses for open and closed habitats to assess if elevational patterns differed between habitat types.
Specifically, we tested the following two hypotheses:
(1)
In open habitats (grasslands), dung beetle diversity should show a decline in biodiversity with increasing elevation as a consequence of increasing harsher conditions (lower temperatures and higher insulation);
(2)
If most of dung beetles prefer open habitats and if the disadvantages of harsher climatic conditions are overwhelmed by the advantages provided by a more favorable habitat type, dung beetle diversity in closed habitats (forests) should show an opposite pattern (i.e., biodiversity should increase with increasing elevation) because tree cover decreases with elevation.
We tested these hypotheses by both comparing diversity curves at different elevations and by investigating elevational variations for specific Hill numbers corresponding to species richness, Shannon diversity, and Simpson diversity in closed and open habitats.
Because of the influence of the sample size in diversity estimates, we conducted our study by using (1) empirical estimates obtained by standardizing data on the basis of the sampling effort (number of active traps over the entire sampling period at each elevation); (2) estimates based on rarefaction/extrapolation to an equal coverage; and (3) asymptotic estimates.

2. Materials and Methods

2.1. Study Area

The study was conducted along an elevational gradient of about 1000 m (from about 950 m to 2000 m) in the Central Apennines (Latium, Central Italy) within the “Montagne della Duchessa” Regional Nature Reserve. The reserve has a surface of 35.4 km2 and an altitudinal range of 1450 m (from 950 to 2239). The Montagne della Duchessa mountains are a Cretaceous carbonatic massif with a landscape composed of peaks above 2000 m and high-sloped valleys created by glacial and karstic erosion.
Between 800 and 1200 m, the vegetation is mainly represented by mixed woods, dominated by oaks (Quercus cerris, Quercus pubescens, and Quercus ilex) associated with hornbeams (Ostrya carpinifolia), ash trees (Fraxinus excelsior), and maples (Acer pseudoplatanus, A. campestris, and A. monspessulanum). Open habitats are occupied by grasslands corresponding to the Habitat of Community Interest 6210, defined as ‘semi-natural dry grasslands and scrubland facies on calcareous substrates (Festuco-Brometalia)’, according to the European Commission classification [63,64]. From 1200 to 1500 m, woods are dominated by the hornbeam, and above 1500, the beech (Fagus sylvatica) is the dominant species. The tree line occurs at 1700–1800 m. Above the tree line, vegetation is represented by juniper brushlands and calcareous grasslands. Secondary meadows due to extensive pasture occur along all the gradient, but these are particularly common at the highest elevations.
The large herbivores present in the study area during the sampling period included about 300–400 deer, 80–100 roe deer, and 300 wild boars as sedentary species. Domesticated herbivores included about 400–450 sheep, 300 cows, and 150 horses, which are moved from lowland to high altitude grasslands in the summer months (June–September).

2.2. Sampling Protocol

The elevational gradient was divided into seven intervals (belts) of equal size (150 m, 950–1100 m, 1100–1250 m, 1250–1400 m, etc.). In each belt, dung beetles were sampled using two baited traps: one placed in a grassland and the other placed in a woodland, with the exception of the highest belt, which was above the tree line and therefore entirely occupied by grasslands.
Beetles were collected using ‘bait-surface-grid’ traps [57,65] made up of buckets (diameter: 186 mm, depth:168 mm) sunk in the ground with their rim at soil level. Each trap was filled with 400 mL of water to prevent dung beetles from escaping the trap after falling in. Traps were covered with a 20 mm mesh iron net, and 1 kg of fresh bovine dung was placed on the net. This type of trap is considered to be one of the most efficient methods to sample dung beetles [57,65,66].
Sampling was conducted once a month from June to November 2018. Each trap was active for 48 h; after this time, the traps were emptied to collect the beetles that had fallen into the buckets, which were left on site to be re-used in the following month [52,57,67,68]. The location of the traps therefore remained constant over the entire sampling period. The samples were washed and stored in 70% ethanol. The dung beetles were identified to the species level using keys [69]. All material is preserved in the senior author’s (SF) collection.

2.3. Data Analysis

We performed analyses both by maintaining forest and grassland samples separately and by merging the two habitats to investigate the impact of not distinguishing habitat types in the study of dung beetle assemblages. In all cases, temporal replicates were aggregated prior analyses. To take into account the loss of certain traps in some months and the lack of wooded habitats in the highest belt, for each belt, the number of collected individuals in a given habitat was standardized proportionally to the minimum number of active traps in that habitat. The same procedure was adopted for analyses in which data from the different habitats were merged.
To analyze how dung beetle diversity varied along the elevational gradient, we calculated diversity profiles using Hill numbers (qD) for each dung beetle sample. Hill numbers differ among themselves only by an exponent, the diversity order q, which determines the measure’s sensitivity to the relative abundances of species. To obtain a complete characterization of an assemblage, a diversity profile is constructed by plotting qD versus q from q = 0 to q = 3 or 4 (beyond this it changes very little, whereas q < 0 are not used, as they are dominated by the abundances of rare species and have poor statistical sampling properties) [24]. We constructed diversity profile curves for q varying between 0 and 3 with increments of 0.25.
A first set of analyses was conducted by calculating the Hill numbers on the data standardized as described above (empirical estimates). As Hill numbers are influenced by sample size, we performed a second set of analyses using a rarefaction/extrapolation procedure [20] on raw data. Extrapolation should not exceed doubling [20]; thus, we used the double of the minimum value of sampled individuals as a common value for rarefaction/extrapolation. Finally, we developed a third set of analyses using asymptotic values [20] to deal with the possibility of undersampling. For each index, the associated 95% confidence intervals were also calculated [20]. Analyses were conducted using the R packages SpadeR version 0.1.1 [70] with the function Diversity to compute empirical and asymptotic diversity estimates and iNEXT version 2.0.20 [20,71] with the function iNEXT to compute rarefied and extrapolated diversity estimates. The R package ggplot2 version 3.3.3 [72] was used for the graphics.
Hill numbers include three widely used species diversity measures as special cases: species richness (q = 0), exponential Shannon diversity (q = 1 as the limit case since the equation for Hill numbers is undefined for q = 1), and Simpson diversity (q = 2) [20,23]. Thus, we modelled how these three particular measures varied with elevation by using generalized linear models (GLMs). We used a Poisson distribution and a log-link for species richness (rounded to integers for estimates obtained after rarefaction/extrapolation or as asymptotes), whereas a Gaussian distribution was used for the Shannon and Simpson diversity. The last procedure was also used to model how standardized abundances varied with elevation. GLMs were applied by using the R package stats [73], and model performances were assessed using a Kullback–Leibler-divergence-based R2. R2 estimates were computed using the R package rsq version 2.2 with the function rsq [74].
Analyses were conducted for all dung beetles (geotrupids + aphodiids + scarabaeids) and for aphodiids and scarabaeids separately. Systematics of Scarabaeoidea is much discussed; more traditional authors include all of the main groups as subfamilies of a few families (with the whole Aphodiidae considered as a subfamily—Aphodiinae—of Scarabaeinae), while others tend to identify a higher number of families (with Aphodiidae as a family on their own hence the separation from the Scarabaeidae) [75]. We adopted this second approach, which is the most commonly used by European taxonomists [75] and that was previously followed in various ecological studies (e.g., [45,46,49,56,57]). However, in our case, the choice is merely nomenclatorial and does not affect the results.

3. Results

In total, we collected 12,324 dung beetles belonging to 42 species (Table 1; Table S1). Aphodiids represented 66% of the collected species and 18% of the collected individuals; scarabaeids accounted for 27% of the collected species and 81% of the collected individuals; geotrupids were 7% of the collected species and less than 2% of the collected individuals. Thus, geotrupids were the less represented taxon in terms of both richness and abundance. Aphodiids were more diversified, but they were less abundant than scarabaeids. The proportional abundance of aphodiids was lower than that of scarabaeids at all elevations except for the first belt, where aphodiids accounted for about 67% of total dung beetles, whereas scarabaeids were 25%. The proportional abundance of aphodiids tended to decline with elevation, whereas that of scarabaeids tended to increase, although the relationships were marginally not significative (Spearman rank correlation tests: Aphodiidae rs = −0.750, Scarabaeidae rs = 0.750; p = 0.052 in both cases). Dung beetles were much more abundant in grasslands than in woodlands (average values of standardized abundances for belts 1–6 ± standard error: 1483.972 ± 636.049 in grasslands, 112.639 ± 24.667 in woodlands; Wilcoxon test for paired data: p = 0.028).
Diversity curves constructed with different approaches (empirical, rarefied/extrapolated, and asymptotic values) showed similar patterns, with substantial differences only occurring for q = 0, which clearly indicates the impact of sampling size on species richness estimates (Figure 1).
Using empirical values, woodland diversity curves clearly indicated a direct relationship between diversity and elevation (Figure 1A) for q ≥ 1.5; however, the confidence intervals were largely overlapping (Figure S1A). By contrast, grassland diversity curves (Figure 1B) showed distinct patterns of decreasing diversity with increasing elevation. For q = 0, diversity values were characterized by large uncertainties (Figure S1B) and did not show any clear relationship with elevation. However, with increasing q-values, confidence intervals became narrower, and diversity values showed distinct elevational patterns. The first belt (950–1100 m) showed the highest diversity followed by the second belt (1100–1250 m). Belts 3–5 (1250–1400, 1400–1550, 1550–1700 m, respectively) showed intermediate values with largely overlapping confidence intervals. The seventh belt (1850–2000 m) showed a lower diversity, and the sixth belt (1700–1850 m) was that with the lowest diversity. When the two habitats were grouped, the emerging pattern was similar to that obtained for the grasslands (Figure 1C and Figure S1C). The first two belts (950–1100 m and 1100–1250 m, respectively) showed overlapping confidence intervals, both with clearly higher diversity values at increasing q-values. The third and fourth belts (1250–1400 and 1400–1550 m, respectively) presented largely overlapping confidence intervals, with diversity values only slightly lower than those of the first and second belts. The fifth belt (1550–1700 m) had a substantially lower diversity, and the sixth and seventh belts (1700–1850 and 1850–2000 m, respectively) showed the lowest diversity.
The use of a rarefaction/extrapolation approach (Figure 1D–F and Figure S1D–F) produced similar results, but with clearer patterns. Again, woodland diversity curves (Figure 1D and Figure S1D) indicated a clear positive relationship between diversity and elevation for q ≥ 1.5 although confidence intervals overlapped by a large amount. By contrast, in the grassland habitat (Figure 1E and Figure S1E), diversity tended to decrease with elevation, even for q = 0. When the two habitats were grouped (Figure 1F and Figure S1F), the emerging pattern was similar to that obtained for the grasslands but with the first and second belts showing more similar diversity values. Curves of asymptotic values (Figure 1G–I, Figure S1G–I) differed from the previous ones, especially for the exaggerated estimates of diversity at q = 0, with extremely large confidence intervals (Figure 1G and Figure S1G).
We focused on the diversity values obtained for q = 0 (species richness), q = 1 (exponential Shannon diversity), and q = 2 (Simpson diversity) to model how they varied with elevation (Figure 2, Figures S2 and S3, Table 2).
In the woodlands, species richness did not show any elevational trend (Figure 2A, Figures S2A and S3A), whereas both the Shannon (Figure 2B, Figures S2B and S3B) and Simpson (Figure 2C, Figures S2C and S3C) diversity increased with elevation (see also Table 2). By contrast, in the grasslands, all three indices decreased with elevation when using values obtained with rarefaction/extrapolation (Figure S2D–F, Table 2).
Using empirical (Figure 2E,F) and asymptotic (Figure S3E,F) values, the Shannon and Simpson indexes had the same behavior, whereas richness exhibited no significant relationship with elevation using empirical estimates (Figure 2D) and a positive relationship using asymptotic values (Figure S3D) (see also Table 2). When the two habitats were merged, both the Shannon and Simpson diversity decreased with elevation (Figure 2H,I, Figures S2H,I and S3H,I), whereas richness showed a significant decrease only when calculated using rarefaction/extrapolation (Figure S2G) (see also Table 2). Abundance increased significantly with altitude for grasslands and for the two habitats merged, but it did not increase for the woodlands (Figure 3, Table 2).
When the aphodiids and scarabaeids were analyzed separately, diversity profiles did not show any clear pattern, with largely overlapping confidence intervals occurring (Figure 4 and Figures S4–S8). Only the curves for the 3rd, 6th, and 7 h belts in the aphodiids (for grasslands and the two habitats merged) appeared well separated (lower diversity) from the others.
Relationships between aphodiid diversity and elevation (Figure 5, Figures S9 and S10, Table 3) followed the patterns obtained with all dung beetle species, but the decline in the Shannon and Simpson indices for the aggregated habitats were not significant because of the opposite trends found in the woodlands and in the grasslands. In the scarabaeids (Figure 6, Figures S11 and S12, Table 4), all of the trends were negative (except for richness in woodlands with empirical and asymptotic estimates) but not significant (except for Shannon estimates in aggregated habitats, which were marginally not significant with empirical values and marginally significant with asymptotic values).
Aphodiid abundance declined with elevation in woodlands (Figure 7A, Table 3) but increased in grasslands (Figure 7B) and in the two habitats when merged (Figure 7C), whereas scarabaeids increased in both habitats and globally (Figure 7D–F, Table 4).

4. Discussion

The most diversified dung beetle group in the study area was the aphodiids, with 27 species, followed by the scarabaeids, with 12 species; the geotrupids were present with 3 species. This situation is representative of the whole Italian fauna, which includes 144 aphodiids, 52 scarabaeids, and 19 geotrupids [69]. Although less diversified, the scarabaeids were much more abundant than the aphodiids. Dung beetle communities in northern Europe are almost completely dominated by aphodiids, whereas scarabaeids dominate in tropical and Mediterranean ecosystems [76]. This pattern can be related to differences in excrement utilization between aphodiids and scarabaeids. According to excrement utilization, dung beetles can be roughly classified into three main categories: dwellers, tunnellers, and rollers [25,76,77], but this classification may be imprecise [78]. Scarabaeids are typically tunnellers (i.e., they dig chambers more or less directly underneath the pat for feeding or breeding) or rollers (i.e., they form a ball of dung that can be rolled away from the pat and buried for feeding or breeding), whereas most aphodiids are dwellers (i.e., they feed in the dung pat as adults and lay eggs within or under the dung mass where they undergo larval development and do not construct nests) [25,46,76,77], but there are aphodiids that construct more or less elaborate pedotrophic nests [78]. In our samples, the only aphodiid species that could not be classified as dwellers were Colobopterus erraticus, Acrossus rufipes, and Bodilopsis rufa [78].
In a cold and humid climate, excrement remains humid for long periods, dung removal is absent or slow, and thus competition among dung beetles is considered reduced; by contrast, in tropical and Mediterranean areas, excrement tends to dry and to be removed quickly, which possibly intensifies competition, favoring species able to protect the offspring through subterranean pedotrophic nests [76,79,80]. This may explain why temperate communities are dominated by dwellers (aphodiids), whereas tropical and Mediterranean communities are dominated by rollers and tunnellers (scarabaeids). The higher abundance of scarabaeids in our study area is consistent with this general pattern and is in line with previous findings in other Mediterranean mountain areas [58].
In Mediterranean areas, climate tends to be more temperate at high altitudes; thus, one might expect a predominance of aphodiids at the highest elevations because of the presence of colder climates. For example, Martin-Piera et al. [57] found that high elevation communities in the Iberian mountains were dominated by aphodiids and geotrupids, and Barbero et al. [81] found that aphodiids were numerically dominant in an Alpine valley (Sagone Valley, Western Alps). In our study, we found that aphodiids were only more abundant than scarabaeids at the lowest elevation, possibly because of the presence of more humid conditions allowed by the tree cover. The higher abundance of scarabaeids at the highest elevations suggests that conditions remain less favorable to aphodiids even with increasing elevation. This can be explained by the fact that average lower temperatures at higher elevations do not prevent excrement from a rapid desiccation, which is probably because of the drying effects of increased windiness and insulation. Thus, even high-altitude Mediterranean dung beetles must cope with the rapid desiccation of the excrement, especially in the summer, when precipitation may be completely absent. In these hot, dry, and insulated environments, scarabaeids can be favored because the use of subterranean pedotrophic nests may reduce the risk of larval death by desiccation.
On the whole, dung beetles were much more abundant in grasslands than in woodlands, which is consistent with the overall preferences of these beetles for open habitats [57,62] and the less abundant presence of potential food in the woodlands. Although deer and roe deer were relatively common in forests, sheep, cows, and horses, which were frequently observed in large numbers in open areas through the entire elevational gradient, were much less frequently observed in the forests, where their pats were also much less abundant. Because of the higher dung beetle abundance in grasslands, the results obtained with the two habitats merged mirrored those obtained for the grasslands alone, but they were less clear because of the confounding effects of the opposite patterns found in the woodlands. The physical structure of habitats can be an important determining factor in the composition and distribution of dung beetle assemblages [49,54,56,57]. Our results highlight the importance of distinguishing open assemblages from the forest ones.
In most cases, previous studies on the altitudinal distribution of dung beetle communities showed a decrease in species with increasing elevation (e.g., [48,50,51,53]), but mid-elevation peaks have been also documented in two Mediterranean areas [55], and in some cases, no relationship was found [57]. In our study, species richness decreased distinctly with elevation in the grassland habitat, whereas no clear pattern was found for the woodland assemblages. Mid-elevation peaks recorded in Mediterranean dung beetles were observed at around 800–900 m in the Iberian Central System and less clearly at around 900–1100 in the Western Rhodopes; above these altitudes, dung beetle richness decreased monotonically [55]. The absence of a mid-elevation peak in our study can be attributed to the fact that our elevational gradient started from 950 m, which is beyond these mid elevations. The difference in the species–elevation patterns observed in grasslands (decreasing number of species with increasing elevation) and woodlands (no relationship) are consistent with trends observed in a Western Anatolian mountain area, where species richness decreased with elevation in open habitats but not in forest habitats [53].
Woodland and grassland assemblages showed different diversity profiles. In general, the diversity curves of woodland assemblages were below those of grassland assemblages. A lower diversity of dung beetle assemblages in wooded habitats seems to be a common situation in the Western Palearctic region, where most of the dung beetles prefer open habitats [57,62], and the presence of open habitats have been increasing by millenary activities of grazing and deforestation [82]. Additionally, the woodland and grassland assemblages showed opposite responses to increasing elevation. In the woodlands, diversity tended to increase with elevation, while in the grasslands, diversity tended to decrease with elevation. Although these opposite patterns can be observed along almost all q-value ranges, they become clearer at increasing q-values, especially for q ≥ 1, which indicates the importance of the contribution of the most dominant species in generating the elevational increase (woodlands) or decrease (grasslands) in diversity. Examination of the elevational variation in the Shannon (q = 1) and Simpson (q = 2) diversity showed an increase in diversity for the woodland assemblages and a decrease for the grassland assemblages. These contrasting patterns can be interpreted as a result of a trade-off in the dung beetle response to elevation between the positive influence of the increasing availability of open habitats, which become more frequent with increasing elevations (although the tree-line in the study area is at around 1700–1800 m, woodland cover tends to decrease with elevation along the whole gradient), and the decrease of optimal environmental conditions (because of higher temperature, lower water availability, and higher insulation). A recent study on Afromontane dung beetle assemblages, where dung beetle abundance, richness, and Shannon diversity declined with elevation, showed positive correlations between mean ground cover and dung beetle abundance and species richness, whereas negative correlations were found between mean canopy cover and dung beetle abundance and species richness [54]. Although from a completely different biome, these results support our conclusion that in the woodlands, diversity increased with increasing elevation because of decreasing tree cover.
The decline of diversity with elevation in grasslands is likely due to increasingly hostile environmental conditions [53] and cannot be related to the absence of large herbivores at higher elevations in the study area, because during the summer months the cattle are moved from lowland to mountain pastures, where they graze at elevations between 1400 and 2100 m. However, further studies should be performed to reveal how environmental conditions act on species diversity patterns along the gradient.
The strong decrease in diversity with increasing elevation in the grassland assemblages was contrasted by an increase in abundances, especially in the scarabaeids. This increase in abundance was mostly due to the very high abundance values attained by the scarabaeid Onthophagus fracticornis, which accounted for more than 90% of the total dung beetles collected in the grasslands at 1700–1850 m and some 80% at 1850–2000 m. Among the aphodiids, Bodilopsis rufa was particularly abundant (some 14% of total captured individuals) in the grassland habitat of the highest belt. These two species are eurytopic dung beetles with wide elevational ranges, and the ability of B. rufa to dig pedotrophic nests [78] may make this species a superior competitor compared to other aphodiids. The very high abundance of dung beetles observed in the 1700–1850 m can be associated with the distribution of grazing activities in the study area, where, during our sampling, cattle were very numerous and frequent in the pastures located at 1700–1800 m, allowing a larger availability of trophic sources to exist here.
Romero-Alcaraz and Ávila [57] found a positive correlation between dung beetle diversity (expressed as average Shannon diversity per trap) and elevation, which seems to be in contrast with our general trends. In fact, the lack of a significant correlation between richness and elevation and the presence of a positive correlation between the Shannon diversity and elevation found in Romero-Alcaraz and Ávila’s study paralleled our results for the woodlands and can be explained by the fact that their sampling sites were mostly placed in woodlands and shrublands. In general, we found that with increasing elevation, assemblages tend to become less diverse because of the increasing dominance of only a few species, which is consistent with a decrease in equitability found in Alpine dung beetle assemblages along an elevational gradient (1750 to 2760 m) in the Maurienne Valley (France) [45]. Thus, the main effect of elevation was that of decreasing dung beetle diversity because environmental conditions become less favorable; however, the abundance of cattle at a high altitude allows a few species to reach very high values of abundance.

5. Conclusions

We found that dung beetle assemblages along a Mediterranean elevational gradient were dominated by scarabaeids, a group of beetles that protect their offspring through subterranean pedotrophic nests. This indicates that warm and dry summer climatic conditions at high elevations favor species that are able to protect their larvae from desiccation. In accordance with our prediction based on the overall preferences of these insects for open habitats, dung beetles were much more abundant and diversified in grasslands than in woodlands. Dung beetle diversity variation with elevation followed opposite patterns in the woodlands and grasslands. In the woodlands, diversity increased with increasing elevation because of decreasing tree cover, which causes the landscape to be more open and hence more favorable to dung beetles. By contrast, in the grasslands, diversity decreased with elevation because of increasingly harsher environmental conditions. These results indicate a trade-off in the beetle response to elevation between the positive effects of the increasing availability of more suitable habitats and the decrease of optimal environmental conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12090781/s1. Table S1: Number of dung beetles collected in two habitats (woodlands and grasslands along an elevational gradient in central Italy. Figure S1. Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy with confidence intervals, Figure S2: Relationship between dung beetle species richness, exponential Shannon diversity, and Simpson diversity and altitude along an elevational gradient in central Italy. Curves were constructed using values based on rarefaction/extrapolation for the two habitats, both separately and merged, Figure S3: Relationship between dung beetle species richness, exponential Shannon diversity, and Simpson diversity and altitude along an elevational gradient in central Italy. Curves were constructed using asymptotic values for the two habitats, both separately and merged, Figure S4: Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy. Curves were constructed for aphodiids and scarabaeids using empirical values, Figure S5: Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy with confidence intervals. Curves were constructed for aphodiids and scarabaeids using empirical values. Figure S5. Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy constructed for aphodiids and scarabaeids using rarefied/extrapolated values, Figure S6: Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy constructed for aphodiids and scarabaeids using rarefied/extrapolated values with confidence intervals, Figure S7: Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy constructed for aphodiids and scarabaeids using asymptotic values, Figure S8: Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy constructed for aphodiids and scarabaeids using asymptotic values with confidence intervals, Figure S9: Relationship between aphodiid species richness, exponential Shannon diversity, and Simpson diversity and altitude along an elevational gradient in central Italy. Curves were constructed using values based on rarefaction/extrapolation for the two habitats, both separately and merged, Figure S10: Relationship between aphodiid species richness, exponential Shannon diversity, and Simpson diversity and altitude along an elevational gradient in central Italy. Curves were constructed using asymptotic values for the two habitats, both separately and merged, Figure S11: Relationship between scarabaeid species richness, exponential Shannon diversity, and Simpson diversity and altitude along an elevational gradient in central Italy. Curves were constructed using values based on rarefaction/extrapolation for the two habitats, both separately and merged, Figure S12: Relationship between scarabaeid species richness, exponential Shannon diversity, and Simpson diversity and altitude along an elevational gradient in central Italy. Curves were constructed using asymptotic values for the two habitats, both separately and merged.

Author Contributions

Conceptualization, S.F.; methodology, S.F. and N.T.; doftware, S.F. and N.T.; validation, C.M.; formal analysis, S.F., N.T., C.M. and S.D.P.; investigation, C.M. and E.P.; resources, S.F.; data curation, C.M., E.P. and S.D.P.; writing—original draft preparation, S.F.; writing—review and editing, S.F.; visualization, N.T. and S.F.; project administration, S.F.; funding acquisition, S.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the department grant Bando Ricerca FFO 2021 ‘Pattern altimetrici della diversità animale, vegetale e microbica’ to Simone Fattorini.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are openly available as Supplementary Materials (Table S1).

Acknowledgments

Riccardo Poloni and Stefano Ziani assisted us with beetle identification. The staff of the Riserva Naturale Montagne della Duchessa allowed us to conduct this research. A special thanks to Emanuela Peria and Giampiero Di Clemente for helping us during the fieldwork and for the information about large herbivores. We are grateful to the three anonymous reviewers for their comments on a previous version of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Peters, M.; Hemp, A.; Appelhans, T.; Behler, C.; Classen, A.; Detsch, F.; Ensslin, A.; Ferger, S.W.; Frederiksen, S.B.; Gebert, F.; et al. Predictors of elevational biodiversity gradients change from single taxa to the multi-taxa community level. Nat. Commun. 2016, 7, 13736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fattorini, S.; Di Biase, L.; Chiarucci, A. Recognizing and interpreting vegetational belts: New wine in the old bottles of a von Humboldt’s legacy. J. Biogeogr. 2019, 46, 1643–1651. [Google Scholar] [CrossRef]
  3. Di Biase, L.; Fattorini, S.; Cutini, M.; Bricca, A. The Role of Inter- and Intraspecific Variations in Grassland Plant Functional Traits along an Elevational Gradient in a Mediterranean Mountain Area. Plants 2021, 10, 359. [Google Scholar] [CrossRef]
  4. Kikvidze, Z.; Pugnaire, F.I.; Brooker, R.W.; Choler, P.; Lortie, C.J.; Michalet, R.; Callaway, R.M. Linking patterns and processes in alpine plant communities: A global study. Ecology 2005, 86, 1395–1400. [Google Scholar] [CrossRef] [Green Version]
  5. Körner, C. Alpine Plant Life: Functional Plant Ecology of High Mountain Ecosystems, 2nd ed.; Springer: Heidelberg, Germany, 1999; p. 349. [Google Scholar]
  6. Körner, C. The use of “altitude” in ecological research. Trends Ecol. Evol. 2007, 22, 569–574. [Google Scholar] [CrossRef]
  7. Körner, C. Alpine Treelines—Functional Ecology of the Global High Elevation Tree Limits, 1st ed.; Springer: Basel, Switzerland, 2012; p. 220. [Google Scholar]
  8. Le Roux, P.C.; McGeoch, M.A. Interaction intensity and importance along two stress gradients: Adding shape to the stress-gradient hypothesis. Oecologia 2010, 162, 733–745. [Google Scholar] [CrossRef]
  9. McCain, C.M.; Grytnes, J.A. Elevational gradients in species richness. In Encyclopedia of Life Sciences (eLS); John Wiley & Sons: Chichester, UK, 2010; pp. 1–10. [Google Scholar]
  10. Hoiss, B.; Krauss, J.; Potts, S.G.; Roberts, S.; Steffan-Dewenter, I. Altitude acts as an environmental filter on phylogenetic composition, traits and diversity in bee communities. Proc. R. Soc. Lond. 2012, 279, 4447–4456. [Google Scholar] [CrossRef] [Green Version]
  11. Sanders, N.J. Elevational gradients in ant species richness: Area, geometry, and Rapoport’s rule. Ecography 2002, 25, 25–32. [Google Scholar] [CrossRef]
  12. Sanders, N.J.; Rahbek, C. The patterns and causes of elevational diversity gradients. Ecography 2012, 35, 1–3. [Google Scholar] [CrossRef]
  13. Fattorini, S. Disentangling the effects of available area, mid-domain constraints, and species environmental tolerance on the altitudinal distribution of tenebrionid beetles in a Mediterranean area. Biodivers. Conserv. 2014, 23, 2545–2560. [Google Scholar] [CrossRef]
  14. Luo, Y.; Liu, J.; Tan, S.; Cadotte, M.W.; Xu, K.; Gao, L.; Li, D. Trait variation and functional diversity maintenance of understory herbaceous species coexisting along an elevational gradient in Yulong Mountain, Southwest China. Plant Divers. 2016, 38, 303–311. [Google Scholar] [CrossRef] [PubMed]
  15. Camacho, L.; Avilés, L. Decreasing predator density and activity explain declining predation of insect prey along elevational gradients. Am. Nat. 2019, 194, 334–343. [Google Scholar] [CrossRef]
  16. Lazarina, M.; Charalampopoulos, A.; Psaralexi, M.; Krigas, N.; Michailidou, D.E.; Kallimanis, A.S.; Sgardelis, S.P. Diversity patterns of different life forms of plants along an elevational gradient in Crete, Greece. Diversity 2019, 11, 200. [Google Scholar] [CrossRef] [Green Version]
  17. Fattorini, S.; Mantoni, C.; Di Biase, L.; Pace, L. Mountain biodiversity and sustainable development. In Encyclopedia of the UN Sustainable Development Goals. Life on Land; Leal Filho, W., Azul, A., Brandli, L., Özuyar, P., Wall, T., Eds.; Springer: Cham, Switzerland, 2020; pp. 1–31. [Google Scholar]
  18. Fattorini, S.; Mantoni, C.; Di Biase, L.; Strona, G.; Pace, L.; Biondi, M. Elevational patterns of generic diversity in the tenebrionid beetles (Coleoptera Tenebrionidae) of Latium (Central Italy). Diversity 2020, 12, 47. [Google Scholar] [CrossRef] [Green Version]
  19. Moradi, H.; Fattorini, S.; Oldeland, J. Influence of elevation on the species-area relationship. J. Biogeogr. 2020, 46, 304–315. [Google Scholar]
  20. Chao, A.; Gotelli, N.J.; Hsieh, T.C.; Sander, E.L.; Ma, K.H.; Colwell, R.K.; Ellison, A.M. Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies. Ecol. Monogr. 2014, 84, 45–67. [Google Scholar] [CrossRef] [Green Version]
  21. Magurran, A.E. Ecological Diversity and Its Measurement; Princeton University Press: Princeton, NJ, USA, 1988; p. 192. [Google Scholar]
  22. Magurran, A.E. Measuring Biological Diversity; Wiley Blackwell: Oxford, UK, 2004; p. 264. [Google Scholar]
  23. Chao, A.; Jost, L. Coverage-based rarefaction and extrapolation: Standardizing samples by completeness rather than size. Ecology 2012, 93, 2533–2547. [Google Scholar] [CrossRef] [PubMed]
  24. Chiu, C.H.; Chao, A. Distance-based functional diversity measures and their decomposition: A framework based on Hill numbers. PLoS ONE 2020, 9, e100014. [Google Scholar] [CrossRef] [Green Version]
  25. Hanski, I. The dung insect community. In Dung Beetle Ecology; Hanski, I., Cambefort, Y., Eds.; Princeton University Press: Princeton, NJ, USA, 1991; pp. 5–21. [Google Scholar]
  26. Gillard, P. Coprophagous beetles in pasture ecosystems. J. Aust. Inst. Agric. Sci. 1967, 33, 30–34. [Google Scholar]
  27. Yokoyama, K.; Kai, H.; Koga, T.; Kawaguchi, S. Effect of dung beetle, Onthophagus lenzii H. on nitrogen transformation in cow dung and dung balls. Soil Sci. Plant Nutr. 1991, 37, 341–345. [Google Scholar] [CrossRef]
  28. Macqueen, A.; Beirne, B.P. Effects of cattle dung and dung beetle activity on growth of beardless wheatgrass in British Columbia. Can. J. Plant Sci. 1975, 55, 961–967. [Google Scholar] [CrossRef] [Green Version]
  29. Kabir, S.M.H.; Howlader, J.A.; Begum, J. Effect of dung beetle activities on the growth and yield of wheat plants. Bangladesh J. Agric. Res. 1985, 10, 49–55. [Google Scholar]
  30. Galbiati, C.; Bensi, C.; Conceiçaõ, C.H.C.; Florcovski, J.L.; Calafiori, M.H. Estudo comparativo entre besouros do esterco, Dichotomius analypticus (Mann, 1829) e Onthophagus gazella (F.), sobre as pastagens, em condiçoes brasileiras. Ecossistema 1995, 20, 109–118. [Google Scholar]
  31. Bertone, M.A. Dung Beetles (Coleoptera: Scarabaeidae and Geotrupidae) in North Carolina Pastures and Their Implications for Pasture Improvement; North Carolina State University: Raleigh, NC, USA, 2004; p. 159. Available online: https://repository.lib.ncsu.edu/bitstream/handle/1840.16/1952/etd.pdf?sequence=1&isAllowed=y (accessed on 30 May 2021).
  32. Bang, H.S.; Lee, J.H.; Kwon, O.S.; Na, Y.E.; Jang, Y.S.; Kim, W.H. Effects of paracoprid dung beetles (Coleoptera: Scarabaeidae) on the growth of pasture herbage and on the underlying soil. Appl. Soil Ecol. 2005, 29, 165–171. [Google Scholar] [CrossRef]
  33. Lastro, E. Dung Beetles (Coleoptera: Scarabaeidae and Geotrupidae) in North Carolina Pasture Ecosystem; North Carolina State University: Raleigh, NC, USA, 2006; p. 121. Available online: https://repository.lib.ncsu.edu/bitstream/handle/1840.16/1498/etd.pdf?sequence=1&isAllowed=y (accessed on 30 May 2021).
  34. Yamada, D.; Imura, O.; Shi, K.; Shibuya, T. Effect of tunneler dung beetles on cattle dung decomposition, soil nutrients and herbage growth. Graesslia 2007, 53, 121–129. [Google Scholar] [CrossRef]
  35. Andresen, E.; Feer, F. The role of dung beetles as secondary seed dispersers and their effect on plant regeneration in tropical rainforests. In Seed Fate: Predation, Dispersal and Seedling Establishment; Forget, P.-M., Lambert, J.E., Hulme, P.E., Vander Wall, S.B., Eds.; CAB International: Wallingford, UK, 2005; pp. 331–349. [Google Scholar]
  36. Bishop, A.L.; McKenzie, H.J.; Spohr, L.J.; Barchia, I.M. Interactions between dung beetles (Coleoptera: Scarabaeidae) and the arbovirus vector Culicoides brevitarsis Kieffer (Diptera: Ceratopogonidae). Aust. J. Entomol. 2005, 44, 89–96. [Google Scholar] [CrossRef]
  37. Hanski, I.; Krikken, J. Dung beetles in tropical forest in Southeast Asia. In Dung Beetle Ecology; Hanski, I., Cambefort, Y., Eds.; Princeton University Press: Princeton, NJ, USA, 1991; pp. 179–197. [Google Scholar]
  38. Davis, A.L.V. Compositional differences between dung beetle (Coleoptera: Scarabaeidae s. str.) assemblages in winter and summer rainfall climates. Afric. Entomol. 1994, 2, 45–51. [Google Scholar]
  39. Davis, A.J.; Huijbregts, H.; Krikken, J. The role of local and regional processes in shaping dung beetle communities in tropical forest plantations in Borneo. Glob. Ecol. Biogeogr. 2000, 9, 281–292. [Google Scholar] [CrossRef]
  40. Feinsinger, P. Designing Field Studies for Biodiversity Conservation; The Nature Conservancy and Island Press: Washington, DC, USA, 2001. [Google Scholar]
  41. de Andrade, R.B.D.; Barlow, J.; Louzada, J.; Vaz-de-Mello, F.Z.; Souza, M.; Silveira, J.M.; Cochrane, M.A. Quantifying responses of dung beetles to fire disturbance in tropical forests: The importance of trapping method and seasonality. PLoS ONE 2011, 6, e26208. [Google Scholar]
  42. Gardner, T.A.; Hernández, M.I.; Barlow, J.; Peres, C.A. Understanding the biodiversity consequences of habitat change: The value of secondary and plantation forests for neotropical dung beetles. J. Appl. Ecol. 2008, 45, 883–893. [Google Scholar] [CrossRef]
  43. Nichols, E.S.; Gardner, T.A. Dung beetles as a candidate study taxon in applied biodiversity conservation research. In Ecology and Evolution of Dung Beetles; Simmons, L.W., Ridsdill-Smith, T.J., Eds.; Blackwell Publishing: Oxford, UK, 2011; pp. 267–291. [Google Scholar]
  44. Lumaret, J.P.; Stiernet, N. Montane dung beetles. In Dung Beetle Ecology; Hanski, I., Cambefort, Y., Eds.; Princeton University Press: Princeton, NJ, USA, 1991; pp. 242–254. [Google Scholar]
  45. Lumaret, J.P.; Stiernet, N. Biogeography of dung beetle communities in the western and central Alps (Coleoptera, Scarabaeoidea). Biogeographia 1992, 16, 425–436. [Google Scholar]
  46. Jay-Robert, P.; Lobo, J.M.; Lumaret, J.P. Altitudinal turnover and species richness variation in European montane dung beetle assemblages. Arct. Antarct. Alp. Res. 1997, 29, 196–205. [Google Scholar] [CrossRef]
  47. Avila, J.M.; Pascual, F. Contribución al conocimiento de los escarabeidos coprófagos (Coleoptera, Scarabaeoidea) de Sierra Nevada: III. Distribución altitudinal y temporal. Boll. Mus. Reg. Sci. Nat. Torino 1988, 6, 217–240. [Google Scholar]
  48. Lobo, J.M.; Halffter, G. Biogeographical and ecological factors affecting the altitudinal variation of mountainous communities of coprophagous beetles (Coleoptera: Scarabaeoidea): A comparative study. Ann. Entomol. Soc. Am. 2000, 93, 115–126. [Google Scholar] [CrossRef]
  49. Errouissi, F.; Jay-Robert, P.; Lumaret, J.P.; Piau, O. Composition and structure of dung beetle (Coleoptera: Aphodiidae, Geotrupidae, Scarabaeidae) assemblages in mountain grasslands of the Southern Alps. Ann. Entomol. Soc. Am. 2004, 97, 701–709. [Google Scholar] [CrossRef]
  50. Escobar, F.; Halffter, G.; Arellano, L. From forest to pasture: An evaluation of the influence of environment and biogeography on the structure of beetle (Scarabaeinae) assemblages along three altitudinal gradients in the Neotropical region. Ecography 2007, 30, 193–208. [Google Scholar] [CrossRef]
  51. Herzog, S.K.; Hamel-Leigue, A.C.; Larsen, T.H.; Mann, D.J.; Soria-Auza, R.W.; Gill, B.D.; Edmonds, W.D.; Spector, S. Elevational distribution and conservation biogeography of Phanaeine dung beetles (Coleoptera: Scarabaeinae) in Bolivia. PLoS ONE 2013, 8, e64963. [Google Scholar] [CrossRef] [Green Version]
  52. Nunes, C.A.; Braga, R.F.; Figueira, J.E.C.; Neves, F.D.S.; Fernandes, G.W. Dung beetles along a tropical altitudinal gradient: Environmental filtering on taxonomic and functional diversity. PLoS ONE 2016, 11, e0157442. [Google Scholar] [CrossRef]
  53. Şenyüz, Y.; Lobo, J.M.; Dindar, K. Altitudinal gradient in species richness and composition of dung beetles (Coleoptera: Scarabaeidae) in an eastern Euro-Mediterranean locality: Functional, seasonal and habitat influences. Eur. J. Entomol. 2019, 116, 309–319. [Google Scholar] [CrossRef]
  54. Stanbrook, R.; Wheater, C.P.; Harris, W.E.; Jones, M. Habitat type and altitude work in tandem to drive the community structure of dung beetles in Afromontane forest. J. Insect Conserv. 2021, 25, 159–173. [Google Scholar] [CrossRef]
  55. Lobo, J.M.; Guéorguiev, B.; Chehlarov, E. Convergences and divergences between two European mountain dung beetle assemblages (Coleoptera, Scarabaeoidea). Anim. Biodivers. Conserv. 2007, 30.1, 83–96. [Google Scholar]
  56. Menéndez, R.; Gutiérrez, D. Altitudinal effects on habitat selection of dung beetles (Scarabaeoidea Aphodiidae) in the northern Iberian peninsula. Ecography 1996, 19, 313–317. [Google Scholar] [CrossRef]
  57. Romero-Alcaraz, E.; Ávila, J.M. Effect of elevation and type of habitat on the abundance and diversity of Scarabaeoid dung beetle (Scarabaeoidea) assemblages in a Mediterranean area from Southern Iberian Peninsula. Zool. Stud. 2000, 39, 351–359. [Google Scholar]
  58. Martín-Piera, F.; Veiga, C.M.; Lobo, J.M. Ecology and biogeography of dung-beetle communities (Coleoptera, Scarabaeoidea) in an Iberian mountain range. J. Biogeogr. 1992, 19, 677–691. [Google Scholar] [CrossRef]
  59. Hanski, I. The community of coprophagous beetles (Coleoptera. Scarabaeidae and Hydrophilidae) in northern Europe. Ann. Zool. Fenn. 1980, 46, 57–73. [Google Scholar]
  60. Baz, A. Selección de macrohábitat por algunas especies y análisis de una comunidad de escarabeidos coprófagos (Coleoptera) del macizo de Ayllón (Sistema Central, España). Ann. Soc. Entomol. Fr. 1988, 24, 203–210. [Google Scholar]
  61. Lumaret, J.P.; Kirk, A.A. South temperate dung beetles. In Dung Beetle Ecology; Hanski, I., Cambefort, Y., Eds.; Princeton University Press: Oxford, UK, 1991; pp. 97–115. [Google Scholar]
  62. Palusci, E.; Mantoni, C.; Strona, G.; Fattorini, S. Wildfire does not affect dung beetle diversity of high-altitude Mediterranean habitats. Int. J. Wildland Fire 2021, 30, 636–642. [Google Scholar]
  63. European Commission. Council Directive 92/43/EEC of 21 May 1992 on the Conservation of Natural Habitats and of Wild Fauna and Flora. 1992. Available online: http://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:31992L0043 (accessed on 30 May 2021).
  64. European Commission. Interpretation Manual of European Union Habitats; Eur 27; European Commission, DG Environment: Strasbourg, France, 2007; p. 142. Available online: http://ec.europa.eu/environment/nature/legislation/habitatsdirective/docs/2007_07_im.pdf (accessed on 30 May 2021).
  65. Lobo, J.M.; Martín–Piera, F.; Veiga, C.M. Las trampas pitfall con cebo, sus posibilidades en el estudio de las comunidades coprófagas de Scarabaeidae (Col.). I. Características determinantes de su capacidad de captura. Rev. Ecol. Biol. Sol. 1988, 25, 77–100. [Google Scholar]
  66. Labidi, I.; Errouissi, F.; Nouira, S. Spatial and temporal variation in species composition, diversity, and structure of Mediterranean dung beetle assemblages (Coleoptera: Scarabaeidae) across a bioclimatic gradient. Environ. Entomol. 2012, 41, 785–801. [Google Scholar] [CrossRef]
  67. Louzada, J.; Lima, A.P.; Matavelli, R.; Zambaldi, L.; Barlow, J. Community structure of dung beetles in Amazonian savannas: Role of fire disturbance, vegetation and landscape structure. Landsc. Ecol. 2010, 25, 631–641. [Google Scholar] [CrossRef]
  68. Feer, F. Variations in dung beetles assemblages (Coleoptera: Scarabaeidae) within two rain forest habitats in French Guiana. Rev. Biol. Trop. 2013, 61, 753–768. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Ballerio, A.; Rey, A.; Uliana, M.; Rastelli, M.; Rastelli, S.; Romano, M.; Colacurcio, L. Coleotteri Scarabeoidei d’Italia. Available online: http://www.societaentomologicaitaliana.it/Coleotteri%20Scarabeoidea%20d’Italia%202014/index.htm (accessed on 30 May 2021).
  70. Chao, A.; Ma, K.H.; Hsieh, T.C.; Chiu, C.H. SpadeR: Species-Richness Prediction and Diversity Estimation with R; Version 0.1.1; 2016; Available online: https://CRAN.R-project.org/package=SpadeR (accessed on 30 May 2021).
  71. Hsieh, T.C.; Ma, K.H.; Chao, A. iNEXT: iNterpolation and EXTrapolation for Species Diversity; R Package Version 2.0.20; 2020; Available online: https://cran.r-project.org/web/packages/iNEXT/index.html (accessed on 30 May 2021).
  72. Wickham, H. ggplot2: Elegant Graphics for Data Analysis, 2nd ed.; Springer: New York, NY, USA, 2016; p. 260. [Google Scholar]
  73. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 30 May 2021).
  74. Zhang, D. rsq: R-Squared and Related Measures; R Package Version 2.2; 2021; Available online: https://CRAN.R-project.org/package=rsq (accessed on 30 May 2021).
  75. Carpaneto, G.M.; Piattella, E.; Valerio, L. Insecta Coleoptera Scarabeoidea. In Checklist and Distribution of the Italian Fauna. 10,000 Terrestrial and Inland Water Species; Ruffo, S., Stoch, F., Eds.; Memorie del Museo Civico di Storia Naturale di Verona: Verona, Italy, 2006; Volume 17, pp. 193–196. [Google Scholar]
  76. Finn, J.A.; Gittings, T. A review of competition in north temperate dung beetle communities. Ecol. Entomol. 2003, 28, 1–13. [Google Scholar] [CrossRef]
  77. Halffter, G.; Edmonds, W.D. The Nesting Behaviour of Dung Beetles (Scarabaeinae): An Ecological and Evolutive Approach; Instituto de Ecologìa: México City, México, 1982; 176p. [Google Scholar]
  78. Tonelli, M. Some considerations on the terminology applied to dung beetle functional groups. Ecol. Entomol. 2021, 46, 772–776. [Google Scholar] [CrossRef]
  79. Hanski, I.; Cambefort, Y. Competition in dung beetles. In Dung Beetle Ecology; Hanski, I., Cambefort, Y., Eds.; Princeton University Press: Oxford, UK, 1991; pp. 305–329. [Google Scholar]
  80. Carpaneto, G.M.; Piattella, E. Competizione interspecifica e andamento stagionale di una comunità coprofaga in un’area verde urbana di Roma (Coleoptera, Scarabaeoidea). Boll. Ass. Romana Entomol. 1990, 44, 67–79. [Google Scholar]
  81. Barbero, E.; Palestrini, C.; Rolando, A. Le comunità di Scarabaeoidea della Val Sangone: Considerazioni ecologiche e zoogeografiche. Biogeographia 1992, 16, 437–449. [Google Scholar] [CrossRef]
  82. Kadiri, N.; Lobo, J.M.; Lumaret, J.-P. Conséquences de l’interaction entre préférences pour l’habitat et quantité de ressources trophiques sur les communautés d’insectes coprophages (Coleoptera: Scarabaeoidea). Acta Oecol. 1997, 18, 107–119. [Google Scholar] [CrossRef]
Figure 1. Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy. Curves were constructed using empirical (AC), rarefied/extrapolated (DF), and asymptotic (GI) values for the two habitats (woodlands and grasslands) separately (woodlands: A,D,G; grasslands: B,E,H) and merged (C,F,I). Dung beetles were sampled in seven altitudinal intervals (belts) of 150 m extent: Belt 1: 950–1100 m, Belt 2: 1100–1250 m, Belt 3: 1250–1400 m, Belt 4: 1400–1550, Belt 5: 1550–1700, Belt 6: 1700–1850, Belt 7: 1850–2000. For ease of reading, 95% confidence intervals are not shown here, but they are reported in Figure S1.
Figure 1. Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy. Curves were constructed using empirical (AC), rarefied/extrapolated (DF), and asymptotic (GI) values for the two habitats (woodlands and grasslands) separately (woodlands: A,D,G; grasslands: B,E,H) and merged (C,F,I). Dung beetles were sampled in seven altitudinal intervals (belts) of 150 m extent: Belt 1: 950–1100 m, Belt 2: 1100–1250 m, Belt 3: 1250–1400 m, Belt 4: 1400–1550, Belt 5: 1550–1700, Belt 6: 1700–1850, Belt 7: 1850–2000. For ease of reading, 95% confidence intervals are not shown here, but they are reported in Figure S1.
Insects 12 00781 g001
Figure 2. Relationship of dung beetle species richness (A,D,G), exponential Shannon diversity (B,E,H), and Simpson diversity (C,F,I) with altitude along an elevational gradient in central Italy. Curves were constructed using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: AC; grasslands: DF) and merged (GI). General linear models were used for fitting. Dashed areas are 95% confidence intervals.
Figure 2. Relationship of dung beetle species richness (A,D,G), exponential Shannon diversity (B,E,H), and Simpson diversity (C,F,I) with altitude along an elevational gradient in central Italy. Curves were constructed using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: AC; grasslands: DF) and merged (GI). General linear models were used for fitting. Dashed areas are 95% confidence intervals.
Insects 12 00781 g002
Figure 3. Relationship between dung beetle species abundance and altitude along an elevational gradient in central Italy. Curves were constructed using standardized abundance values for the two habitats (woodlands and grasslands) separately (woodlands: (A), grasslands: (B) and merged (C)). General linear models were used for fitting. Dashed areas are 95% confidence intervals. Confidence intervals in panels (B,C) are not distinguishable because of their narrowness.
Figure 3. Relationship between dung beetle species abundance and altitude along an elevational gradient in central Italy. Curves were constructed using standardized abundance values for the two habitats (woodlands and grasslands) separately (woodlands: (A), grasslands: (B) and merged (C)). General linear models were used for fitting. Dashed areas are 95% confidence intervals. Confidence intervals in panels (B,C) are not distinguishable because of their narrowness.
Insects 12 00781 g003
Figure 4. Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy. Curves were constructed for aphodiids (AC) and scarabaeids (DF) using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: (A,D); grasslands: (B,E) and merged (C,F)). Dung beetles were sampled in seven altitudinal intervals (belts) of 150 m extent: Belt 1: 950–1100 m, Belt 2: 1100–1250 m, Belt 3: 1250–1400 m, Belt 4: 1400–1550, Belt 5: 1550–1700, Belt 6: 1700–1850, Belt 7: 1850–2000. For ease of reading, 95% confidence intervals are not shown here, but they are reported in Figure S4.
Figure 4. Diversity curves using Hill numbers for dung beetle assemblages along an elevational gradient in central Italy. Curves were constructed for aphodiids (AC) and scarabaeids (DF) using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: (A,D); grasslands: (B,E) and merged (C,F)). Dung beetles were sampled in seven altitudinal intervals (belts) of 150 m extent: Belt 1: 950–1100 m, Belt 2: 1100–1250 m, Belt 3: 1250–1400 m, Belt 4: 1400–1550, Belt 5: 1550–1700, Belt 6: 1700–1850, Belt 7: 1850–2000. For ease of reading, 95% confidence intervals are not shown here, but they are reported in Figure S4.
Insects 12 00781 g004
Figure 5. Relationships between aphodiid species richness (A,D,G), exponential Shannon diversity (B,E,H), and Simpson diversity (C,F,I) and altitude along an elevational gradient in central Italy. Curves were constructed using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: AC; grasslands: DF) and merged (GI). General linear models were used for fitting. Dashed areas are 95% confidence intervals.
Figure 5. Relationships between aphodiid species richness (A,D,G), exponential Shannon diversity (B,E,H), and Simpson diversity (C,F,I) and altitude along an elevational gradient in central Italy. Curves were constructed using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: AC; grasslands: DF) and merged (GI). General linear models were used for fitting. Dashed areas are 95% confidence intervals.
Insects 12 00781 g005
Figure 6. Relationship of scarabaeid species richness (A,D,G), exponential Shannon diversity (B,E,H), and Simpson diversity (C,F,I) and altitude along an elevational gradient in central Italy. Curves were constructed using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: AC; grasslands: DF) and merged (GI). General linear models were used for fitting. Dashed areas are 95% confidence intervals.
Figure 6. Relationship of scarabaeid species richness (A,D,G), exponential Shannon diversity (B,E,H), and Simpson diversity (C,F,I) and altitude along an elevational gradient in central Italy. Curves were constructed using empirical values for the two habitats (woodlands and grasslands) separately (woodlands: AC; grasslands: DF) and merged (GI). General linear models were used for fitting. Dashed areas are 95% confidence intervals.
Insects 12 00781 g006
Figure 7. Relationship between dung beetle species abundance and altitude along an elevational gradient in central Italy. Curves were constructed separately for aphodiids (AC) and scarabaeids (DF) using standardized abundance values for the two habitats separately (woodlands: A,D and grasslands: B,E) and merged (C,F). General linear models were used for fitting. Dashed areas are 95% confidence intervals. Confidence intervals in panels E and F are not distinguishable because of their narrowness.
Figure 7. Relationship between dung beetle species abundance and altitude along an elevational gradient in central Italy. Curves were constructed separately for aphodiids (AC) and scarabaeids (DF) using standardized abundance values for the two habitats separately (woodlands: A,D and grasslands: B,E) and merged (C,F). General linear models were used for fitting. Dashed areas are 95% confidence intervals. Confidence intervals in panels E and F are not distinguishable because of their narrowness.
Insects 12 00781 g007
Table 1. Number of individuals captured for each species of dung beetle collected along an elevational gradient in central Italy.
Table 1. Number of individuals captured for each species of dung beetle collected along an elevational gradient in central Italy.
TaxonomyIndividuals in WoodlandsIndividuals in GrasslandsTotal Number of Collected Individuals
Geotrupidae
Geotrupes spiniger Marsham, 1802172643
Sericotrupes niger (Marsham, 1802)34107141
Trypocopris vernalis (Linnaeus, 1758)02121
Aphodiidae
Acrossus depressus (Kugelann, 1792)101
Acrossus rufipes (Linnaeus, 1758)401151
Amidorus obscurus (Fabricius, 1792)145
Amidorus thermicola (Sturm, 1800)033
Aphodius coniugatus (Panzer, 1795)01313
Aphodius fimetarius (Linnaeus, 1758)04848
Bodiloides ictericus(Laicharting, 1781)011
Bodilopsis rufa (Moll, 1782)48673721
Chilothorax conspurcatus (Linnaeus, 1758)011
Chilothorax paykulli (Bedel, 1907)011
Colobopterus erraticus (Linnaeus, 1758)03535
Coprimorphus scrutator (Herbst, 1789)011
Esymus pusillus pusillus (Herbst, 1789)171027
Euheptaulacus carinatus (Germar, 1824)3300303
Labarrus lividus (Olivier, 1789)01515
Limarus zenkeri (Germar, 1813)255141396
Melinopterus consputus (Creutzer, 1799)134
Melinopterus prodromus (Brahm, 1790)022
Nimbus contaminatus (Herbst, 1783)09393
Nimbus johnsoni (Baraud, 1976)17955234
Nimbus obliteratus (Panzer, 1823)212445
Otophorus haemorrhoidalis (Linnaeus, 1758)011
Phalacronotus biguttatus (Germar, 1824)011
Planolinoides borealis (Gyllenhal, 1827)101
Planolinus fasciatus (Olivier, 1789)5862120
Sigorus porcus (Fabricius, 1792)35457
Trichonotulus scrofa (Fabricius, 1787)011
Scarabaeidae
Copris umbilicatus Abeille de Perrin, 1901022
Euoniticellus fulvus (Goeze, 1777)055
Euonthophagus gibbosus (Scriba, 1790)044
Onthophagus coenobita (Herbst, 1783)34043
Onthophagus fracticornis (Preyssler, 1790)4174837524
Onthophagus illyricus (Scopoli, 1763)077
Onthophagus joannae Goljan, 195310980990
Onthophagus lemur (Fabricius, 1781)03232
Onthophagus taurus (Schreber, 1759)167
Onthophagus vacca (Linnaeus, 1767)011
Onthophagus verticicornis (Laicharting, 1781)4710651112
Sisyphus schaefferi (Linnaeus, 1758)1210211
Total73111,38812,324
Table 2. Results of generalized linear models for the relationships of diversity (AC) and abundance (D) of dung beetles with altitude along an elevational gradient in central Italy. Diversity (expressed as species richness, Shannon diversity, and Simpson diversity) was calculated using empirical (A), rarefied/interpolated (B), and asymptotic (C) values. Analyses were conducted for the two habitats (woodlands and grasslands), both separately and merged. SE = standard error, p = probability, R2 = Kullback–Leibler-divergence-based R2.
Table 2. Results of generalized linear models for the relationships of diversity (AC) and abundance (D) of dung beetles with altitude along an elevational gradient in central Italy. Diversity (expressed as species richness, Shannon diversity, and Simpson diversity) was calculated using empirical (A), rarefied/interpolated (B), and asymptotic (C) values. Analyses were conducted for the two habitats (woodlands and grasslands), both separately and merged. SE = standard error, p = probability, R2 = Kullback–Leibler-divergence-based R2.
WoodlandsGrasslandsHabitats Merged
EstimateSEpR2EstimateSEpR2EstimateSEpR2
(A) Empirical diversity values
Richness0.00010.00050.86010.00820.00020.00030.50950.12190.00020.00030.50300.0900
Shannon0.00400.00100.01670.7965−0.00900.00160.00270.8591−0.00790.00090.00030.9438
Simpson0.00400.0002<0.00010.9942−0.00590.00140.00700.7945−0.00480.00060.00040.9346
(B) Rarefied/interpolated
diversity values
Richness0.00020.00050.76240.0491−0.00120.00040.00380.9321−0.00050.00030.08280.5295
Shannon0.00410.00070.00520.8852−0.00820.00150.00250.8625−0.00810.00090.00020.9464
Simpson0.00410.00030.00020.9756−0.00560.00120.00630.8036−0.00490.00060.00040.9358
(C) Asymptotic diversity values
Richness0.00050.00040.21030.12070.00070.00020.00200.80710.00040.00020.04200.1734
Shannon0.00460.00110.01350.8167−0.00940.00170.00260.8598−0.00820.00090.00030.9440
Simpson0.00440.00030.00020.9777−0.00610.00140.00740.7907−0.00490.00060.00040.9359
(D) Abundance−0.00110.0001<0.00010.28210.0026<0.0001<0.00010.74310.0023<0.0001<0.00010.7174
Table 3. Results of generalized linear models for the relationships between diversity (AC) and abundance (D) of aphodiid dung beetles and altitude along an elevational gradient in central Italy. Diversity (expressed as species richness, Shannon diversity, and Simpson diversity) was calculated using empirical (A), rarefied/interpolated (B), and asymptotic (C) values. Analyses were conducted for the two habitats (woodlands and grasslands), both separately and merged. SE = standard error, p = probability, R2 = Kullback–Leibler-divergence-based R2.
Table 3. Results of generalized linear models for the relationships between diversity (AC) and abundance (D) of aphodiid dung beetles and altitude along an elevational gradient in central Italy. Diversity (expressed as species richness, Shannon diversity, and Simpson diversity) was calculated using empirical (A), rarefied/interpolated (B), and asymptotic (C) values. Analyses were conducted for the two habitats (woodlands and grasslands), both separately and merged. SE = standard error, p = probability, R2 = Kullback–Leibler-divergence-based R2.
WoodlandsGrasslandsHabitats Merged
EstimateSEpR2EstimateSEpR2EstimateSEpR2
(A) Empirical diversity values
Richness0.00010.00060.85990.01050.00030.00040.35140.32690.00030.00030.35900.2248
Shannon0.00280.00080.02720.7432−0.00430.00150.03320.6294−0.00060.00180.76900.0189
Simpson0.00270.00050.00670.8693−0.00370.00120.02870.6494−0.00010.00120.93800.0013
(B) Rarefied/interpolated
diversity values
Richness0.00020.00060.70410.066−0.00040.00040.30300.5069<0.00010.00030.95900.0007
Shannon0.00290.00070.01470.8091−0.00440.00140.02500.6668−0.00070.00180.70430.0313
Simpson0.00280.00060.00800.8573−0.00370.00120.02670.6584−0.00010.00120.92610.0019
(C) Asymptotic diversity values
Richness0.00040.00060.45660.09720.00100.00030.00220.75360.00070.00030.01020.3756
Shannon0.00310.00090.02140.7711−0.00470.00160.02860.6497−0.00070.00190.73470.0250
Simpson0.00290.00060.00820.8560−0.00400.00130.02650.6594−0.00010.00120.92900.0018
(D) Abundance−0.00160.0002<0.00010.44440.00200.0001<0.00010.72630.00080.0001<0.00010.3090
Table 4. Results of generalized linear models for the relationships between diversity (AC) and abundance (D) of scarabaeid dung beetles and altitude along an elevational gradient in central Italy. Diversity (expressed as species richness, Shannon diversity and Simpson diversity) was calculated using empirical (A), rarefied/interpolated (B), and asymptotic (C) values. Analyses were conducted for the two habitats (woodlands and grasslands), both separately and merged. SE = standard error, P = probability, R2 = Kullback–Leibler-divergence-based R2.
Table 4. Results of generalized linear models for the relationships between diversity (AC) and abundance (D) of scarabaeid dung beetles and altitude along an elevational gradient in central Italy. Diversity (expressed as species richness, Shannon diversity and Simpson diversity) was calculated using empirical (A), rarefied/interpolated (B), and asymptotic (C) values. Analyses were conducted for the two habitats (woodlands and grasslands), both separately and merged. SE = standard error, P = probability, R2 = Kullback–Leibler-divergence-based R2.
WoodlandsGrasslandsHabitats Merged
EstimateSEpR2EstimateSEpR2EstimateSEpR2
(A) Empirical diversity values
Richness<−0.00010.00100.97700.0036−0.00020.00050.61360.0492−0.00020.00050.61360.0492
Shannon−0.00180.00120.22420.4373−0.00230.00100.08070.4881−0.00240.00090.05100.5660
Simpson−0.00170.00110.21380.4519−0.00150.00090.15390.3607−0.00160.00080.11210.4259
(B) Rarefied/interpolated
diversity values
Richness−0.00110.00120.36850.3867−0.00080.00090.32050.6141−0.00050.00060.40460.0716
Shannon−0.00220.00130.18410.4962−0.00110.00080.23270.2693−0.00090.00150.56610.0701
Simpson−0.00210.00120.18520.4945−0.00080.00070.29250.2167−0.00050.00110.64750.0451
(C) Asymptotic diversity values
Richness<−0.00010.00100.97700.0036−0.00020.00040.69120.0145−0.00010.00040.84110.0036
Shannon−0.00220.00150.24360.4111−0.00230.00100.07210.5083−0.00250.00090.04650.5805
Simpson−0.00260.00160.18840.4890−0.00150.00090.14900.3670−0.00160.00080.10810.4331
(D) Abundance0.00190.0004<0.00010.84190.0027<0.0001<0.00010.68380.0027<0.0001<0.00010.6825
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mantoni, C.; Tsafack, N.; Palusci, E.; Di Pietro, S.; Fattorini, S. Diversity Patterns of Dung Beetles along a Mediterranean Elevational Gradient. Insects 2021, 12, 781. https://doi.org/10.3390/insects12090781

AMA Style

Mantoni C, Tsafack N, Palusci E, Di Pietro S, Fattorini S. Diversity Patterns of Dung Beetles along a Mediterranean Elevational Gradient. Insects. 2021; 12(9):781. https://doi.org/10.3390/insects12090781

Chicago/Turabian Style

Mantoni, Cristina, Noelline Tsafack, Ettore Palusci, Stefano Di Pietro, and Simone Fattorini. 2021. "Diversity Patterns of Dung Beetles along a Mediterranean Elevational Gradient" Insects 12, no. 9: 781. https://doi.org/10.3390/insects12090781

APA Style

Mantoni, C., Tsafack, N., Palusci, E., Di Pietro, S., & Fattorini, S. (2021). Diversity Patterns of Dung Beetles along a Mediterranean Elevational Gradient. Insects, 12(9), 781. https://doi.org/10.3390/insects12090781

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop