Next Article in Journal
The Complete Mitochondrial Genomes of Aelia sibirica and A. fieberi (Hemiptera, Pentatomidae), and Phylogenetic Implications
Next Article in Special Issue
A Mitogenome-Based Phylogeny of Pilargidae (Phyllodocida, Polychaeta, Annelida) and Evaluation of the Position of Antonbruunia
Previous Article in Journal
The Forest Refugium of the Bükk Mountains, Hungary—Vegetation Change and Human Impact from the Late Pleistocene
Previous Article in Special Issue
Cosmetopoeia of African Plants in Hair Treatment and Care: Topical Nutrition and the Antidiabetic Connection?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Vegetation in Elevational Diversity Patterns of Tenebrionid Beetles in Central Italy

Department of Life, Health and Environmental Sciences, University of L’Aquila, Via Vetoio, 67100 L’Aquila, Italy
Diversity 2024, 16(2), 110; https://doi.org/10.3390/d16020110
Submission received: 26 December 2023 / Revised: 27 January 2024 / Accepted: 6 February 2024 / Published: 8 February 2024
(This article belongs to the Special Issue Diversity in 2023)

Abstract

:
Vegetation tends to vary in a systematic fashion along elevational gradients, leading to the possibility of recognizing distinct vegetational belts, which are frequently used to describe and interpret elevational variations in biodiversity. However, anthropogenic changes can create landscapes dominated by secondary grasslands in areas formerly occupied by forests, thus altering the natural sequence of vegetation types. The present research illustrates how the distribution of tenebrionid beetles in central Italy is influenced by secondary vegetation. Classical schemes of vegetational belts were modified into a scheme of main vegetation types that include secondary vegetations. Tenebrionid species presence/absence in each vegetation type was then assessed. Species richness tended to decrease with elevation in both natural and secondary vegetations. Geophilous (ground-dwelling) species (which include many endemics) prevailed in natural and secondary grasslands, while xylophilous species (associated with trees) prevailed in the natural forests. Similarities in tenebrionid composition indicated the presence of two main groups: one associated with forests and the other with natural and secondary grasslands. Geophilous species prevailed among tenebrionids with Mediterranean distributions, whereas xylophilous species prevailed among species distributed mainly in Europe and the Palearctic. High values of richness, biogeographical complexity and proportion of endemics make secondary vegetations of high conservation concern.

1. Introduction

Centuries or even millennia of low-intensity agricultural practices (especially mowing) and pastoralism (grazing) have profoundly altered European landscapes, turning natural forest vegetations into open areas represented by seminatural grasslands [1,2,3]. These landscapes, which are important biodiversity reservoirs [3,4,5,6,7], are now menaced by higher human pressures (e.g., intensification of agricultural practices and urbanization, afforestation, intensive grazing) and spontaneous succession towards woody vegetation after the abandonment of traditional land management [3,8,9,10].
Located in the middle of the Mediterranean global biodiversity hotspot [11,12,13,14], Italy occupies a prominent position for European biodiversity conservation due to its complex biogeographical history [15,16,17,18], being home for exceptionally high numbers of plant and animal species, many of which are endemic [19,20,21,22,23]. The Italian territory is also characterized by a huge variety of landscapes, with many resulting from transformations produced over millennia of human activities [24,25,26,27].
About 77% of the Italian territory is classified as montane (35%: areas occupied by relief higher than 600–700 m) or hilly (42%: areas with lower elevations, but not flat) [28]. In particular, the Italian territory is largely dominated by two main chains: the Alps in the north, and the Apennines along the peninsula [29]. Therefore, Italian landscapes are profoundly influenced by the mountainous nature of the territory [24,25,26,27], with vegetation types that vary with elevation [26]. The systematic variation of vegetation along elevational gradients—one of the most universal ecological patterns—leads to the possibility of recognizing distinct vegetational belts that roughly parallel the latitudinal variation of biomes [30,31]. However, millennia of anthropogenic changes in the Apennines, such as agriculture and grazing [24,26,32], created landscapes dominated by grasslands at elevations which would be naturally occupied by forests [3]. Thus, Apennine landscapes offer the intriguing opportunity of investigating how secondary open vegetations influenced elevational patterns of biodiversity.
The aim of the present paper is to investigate how the presence of secondary vegetation influences the distribution of a group of insects, tenebrionid beetles (Coleoptera Tenebrionidae), along the whole elevational gradient of a region (Latium) of Central Apennines. Because of millenary human activities, e.g., [24,25,26,33,34,35,36,37], most of the lowland areas of Latium have been transformed into agricultural fields and grasslands, as have vegetations on mountains, with the current presence of secondary grasslands interspersed within more or less well-preserved natural vegetations through the whole elevational gradient. Previous research showed that tenebrionid diversity in Latium declines with elevation mainly because of the filtering effects of increasingly harsher conditions with increasing elevation [38,39,40], but the role of different types of vegetation in determining elevational patterns in tenebrionid diversity has been not investigated. Under the assumption of environmental filtering processes, we expect that insect species richness should decline with elevation in both natural and secondary vegetations, independently from their ecological preferences and biogeographical origin. However, because of the highly contrasting nature of open and forest vegetations, we expect that species assemblages that can be found in secondary vegetations within elevational belts that should be naturally occupied by forests are not a simplified version (i.e., remnants) of forest assemblages, but very different species assortments arising from natural open habitats occurring in other belts. Tenebrionids are particularly well suited to investigate this issue as they can be classified into two main ecological groups that can be distinctly associated with open and forest vegetations: geophilous (i.e., those associated with soil) and xylophilous (i.e., those associated with trees) [41]. Whereas geophilous species are expected to be more common in open vegetation, xylophilous are expected to be associated mainly with tree vegetations.
To test these hypotheses, we investigated here the elevational patterns of tenebrionid beetles in Latium with reference to the sequence of the main vegetation types along the elevational gradient, taking into account their association with forest and secondary open vegetations. Specifically, we tested the following predictions:
(1)
Prediction 1. If tenebrionid richness varies along the elevational gradient in response to a filtering process [39,40], declining patterns in species richness should be observed in both forest and secondary open vegetations. To test this prediction, richness was calculated separately for natural and secondary vegetations, and the respective elevational patterns compared.
(2)
Prediction 2. If tenebrionids of secondary open vegetations are not an impoverished version of the communities inhabiting former forests, but are composed of species associated with natural open vegetations, we expect a low similarity in species composition between tenebrionid of secondary and natural vegetations even within the same elevational belt, but a high similarity between secondary and open natural vegetations. Additionally, we expect a high relative abundance of geophilous species even at elevations in which xylophilous species should be more frequent. To test this prediction, similarities in species composition among vegetation types and levels of association between ecological groups (geophilous vs. xylophilous) and vegetation types were investigated.
(3)
Prediction 3. Previous research showed that biogeographical composition of tenebrionid beetles in Latium changed with elevation, with species mainly distributed in the Mediterranean basin declining with elevation, being less tolerant to decreasing temperatures [38]. As a result of this filtering process, it can be expected that their proportion declines with elevation in both secondary and natural vegetations. Variation in the incidence of Mediterranean species in different vegetation types was analyzed for testing this prediction.
(4)
Prediction 4. Not only tenebrionid species richness, but also biogeographical diversity declines with increasing elevation [38]. Under the filtering hypothesis, this decline in biogeographical diversity should be observed in both natural and secondary vegetations. To test this prediction, diversity profiles based on the biogeographical composition of tenebrionid assemblages were constructed for the different vegetation types.
(5)
Prediction 5. The proportion of tenebrionid endemic species is known to increase with elevation [38]. This may be explained by the fact that they are almost entirely represented by geophilous tenebrionids, which should be more associated with open habitats that naturally dominate the highest elevations. If this assumption is correct, endemics should also represent a conspicuous component of tenebrionid assemblages of secondary open vegetations at high elevations. For this prediction, frequencies of endemics were compared among vegetation types.

2. Materials and Methods

2.1. Study Area

The study was conducted in Latium (central Italy, Figure 1). This region extends for 17,242 km2, of which 20% is flat, 54% hilly, and 26% mountainous [28]. In coastal areas, sandy beaches are the main type of habitat with also some remnants of dunes, while rocky shores are restricted to small areas in the southernmost part of the region. Coastal plains are largely occupied by cultivated fields, interspersed with remnants of natural vegetation mainly represented by Mediterranean shrublands and maquis, with rare fragments of hygrophilous woodlands. The mountainous territory of the region includes both pre-Apennine and Apennine mountains. The pre-Apennines include groups of mountains mainly occupied by mixed broadleaved forests, including typical evergreen broadleaves (Quercus ilex L. and Quercus suber L.) and chestnut woods (Castanea sativa L.). The Apennines of Latium are a continuation of the Apennines of Abruzzo, and are mainly characterized by the presence of deciduous oak woods (Quercus cerris L., Quercus pubescens Willd.) and beech woods (Fagus sylvatica L.). For further details on the vegetation of Latium, see, for example, [42,43]. Because of its complex topography (the highest peak is Mount Gorzano, with 2458 m, on the border with the Abruzzo region), the region shows a wide variety of climates, ranging from typically Mediterranean climates along the coasts (with average annual temperatures even around 17 °C, less than 1000 mm of annual precipitation, and a distinct period of summer aridity) to subalpine climates on the highest peaks (with average annual temperatures < 10 °C, about 1500 mm of annual precipitation, and a lack of a period of summer aridity) [44,45].
For the purposes of this research, classical schemes of vegetational belts used for Italian mountains [46,47,48,49,50] were modified into a scheme of main vegetation types that includes secondary open vegetations [51]:
  • Dunes (about sea level). This belt includes the vegetation of mobile and fixed dunes and extends inland to a very limited extent. Because of coastal exploitation, this habitat is now currently highly fragmented and under severe threats, including those posed by climate change [52,53]. It represents a natural open vegetation.
  • Mediterranean shrubland and maquis: holm-oak forests, cork oak forests, coastal pinewoods, Mediterranean maquis with xerophilic species (~0–800 m). This belt is represented by a mosaic of vegetation types, ranging from garrigues to Mediterranean forests. Coastal pinewoods were mostly planted centuries ago and are mixed with natural vegetation. This belt roughly corresponds to the Mediterranean belt of Pignatti [49]. It represents a natural and more or less closed (forest) vegetation.
  • Arid grasslands of lowlands and low hills (~0–800 m). This form of secondary open vegetation is the result of anthropic activities (low intensity agriculture and grazing) within the former belt.
  • Mixed forests (Turkey oaks and chestnuts) (~400–1000 m). This belt, which roughly corresponds to the Samnite belt of Pignatti [49], is characterized by mixed deciduous forests. The presence of large chestnut groves is anthropogenic, but ancient (for example, on the Abani Hills, they date back to the 17th century); currently they are managed as coppices or orchards for fruit production [54,55]. It represents a substantially natural closed vegetation.
  • Grasslands of medium and high hills (~400–1000 m). This form of secondary open vegetation is the result of anthropic activities (especially grazing) within the former belt.
  • Beech forests (~700–1800 m). This belt is characterized by pure or mixed forests of Fagus sylvatica and roughly corresponds to the sub-Atlantic belt of Pignatti [49]. In fact, in peninsular Italy, old-growth primary forest remnants can only be found at high elevation (say, >1600 m), whereas most beech forests are secondary old-growth forests, but occupy areas where they would be the natural vegetation anyway [56,57]. This vegetation represents a substantially natural closed vegetation.
  • Montane grasslands (~700–1800 m). This form of secondary open vegetation is the result of anthropic activities (especially grazing) within the former belt.
  • Pseudo-alpine grasslands (beyond the tree line) (~1600–2000 m). This belt is represented by the grassland vegetation above the tree line and corresponds to the Mediterraneo-altomontane belt of Pignatti [49]. It represents a substantially natural open vegetation.
For the purposes of this paper, vegetations 1, 2, 4, 6, and 8 are considered as natural (either open or forest) vegetations, although they were also subject to human influence, and 3, 5, 7 as secondary open vegetation. Examples of these vegetation types are illustrated in Figure 2.
Figure 1a illustrates the distribution of these vegetations deduced from Corine Land Cover (CLC) categories as given in [58]. CLCs (fourth and five level) falling in the same vegetation type were merged (e.g., CLCs 3.2.3, 3.1.1.1, 3.1.2.1, 3.1.3.1.1, and 3.1.3.2.1 were merged to obtain the distribution of Mediterranean shrubland and maquis) whereas the different types of secondary vegetations all corresponding to the same CLC 3.2.1 were distinguished on the basis of their elevation.

2.2. Data Collection

The elevational distribution of 85 tenebrionid species and subspecies (hereafter species for simplicity, as the dividing line between species and subspecies, as applied to these insects, is arguably arbitrary [59,60,61,62,63,64,65]) from Latium was assessed using a database comprising about 4600 records spanning from 1860 to 2020 [38,39,40,66], plus new records. These data originated from a critical review of available literature and examination of material preserved in public and private collections, for a total of about 29,000 specimens. Species were collected using a variety of methods, including hand searching (on the ground, under stones and fallen trees, on trees, under bark, etc.), pitfall traps, light traps, aerial traps, soil examination (with sand and litter sieving), in both open and forest vegetations from sites distributed through the whole elevational gradient (Figure 1b). In particular, the use of pitfall traps, litter sieving and hand searching on the ground assures that geophilous species were adequately sampled also in forest habitats. Overall, because of the variety of sources used, they form a random sample not affected by biases due to collector preferences for certain areas, habitats, or species. The use of belts based on vegetation settings, instead of elevational limits, ensures that the results are not influenced by species elevational shifts due to climate changes, which were in fact too small to affect overall patterns [67].
Species presence/absence in each vegetation type was based on species distribution, species ecological preferences resulting from literature, e.g., [68,69,70,71] and personal observations conducted in Latium during the last 30 years in all the aforementioned vegetation types. The genus Lagria (traditionally regarded as a member of a separate family, Lagriidae) and the subfamily Alleculinae (traditionally regarded as a distinct family, Alleculidae) were excluded, as their ecological characteristics are very different from that other tenebrionids (both are mainly flower-visiting insects), and faunal and ecological information on these beetles are still scant. Introduced and synanthropic species associated with human food were also not considered. Nomenclature followed Iwan and Löbl [72]. The validity of Leptoderis italicus Solier 2015 [73] and Pachychila italica Ferrer, 2018 [74] should be confirmed, and, for the present, they are listed as Pachychila frioli Solier, 1835 and Leptoderis collaris (Linnaeus, 1767), respectively.
On the basis of their biology, tenebrionids were divided into two broad ecological categories [41]: geophilous (i.e., ground-dwelling species, with larvae living on or in the soil, and adults usually found under stones, or in the sand) and xylophilous (i.e., species associated with plants, especially trees, whose larvae usually develop in rotten wood and adults occupy subcortical spaces). For a few species that can occupy both habitats, the preferred habitat was considered.
Tenebrionids were classified biogeographically by using chorotypes. Chorotypes are groups of generalized distributional patterns based on inductive methods: individual species distributions are mapped and species with similar distributions are grouped together to define chorotypes as recurrent patterns [75,76,77]. Once a chorotype is defined, any additional species exhibiting a similar distribution can be ascribed to that chorotype. Species assigned to the same chorotype can be members of entirely different taxonomic or ecological groups. Analyzing the chorotype composition of local species assemblages might help to identify the ecological and historical reasons that have produced these assemblages, since species with comparable distributions should also have similar macroecological demands [78]. Chorotypes used in this paper followed Vigna Taglianti et al. [79] and were subsequently grouped into fewer groups of wider meaning, as follows: European (Centraleuropean, European, and S-European chorotypes: 32%), Palearctic (Centralasiatic-Europeo-Mediterranean, Holarctic, Sibero-European, Turano-Europeo-Mediterranean, and Turano-Mediterranean chorotypes: 14%), and Mediterranean (Mediterranean, N-African, W-Mediterranean, and E-Mediterranean chorotypes: 38%). Species that occur only within the Italian territory were classified as endemic (16%). All data (species classification as either geophilous or xylophilous, chorotype, general distribution, and occurrence in the various vegetations) are provided in Table S1.

2.3. Data Analysis

To test Prediction 1, elevational variation in species richness was investigated for native and secondary vegetations separately.
To test Prediction 2, similarity in species composition between vegetation types was investigated by using two coefficients of dissimilarity (β-diversity): the Dice–Sørensen coefficient (βsor, which expresses the total β-diversity) and the Simpson coefficient (βsim, which expresses the pure turnover component, i.e., the compositional differences after removing the effect of nestedness, that is the compositional change caused by ordered species loss) [80,81]). Distance matrices of βsor and βsim were then subject to non-metric multidimensional scaling (NMDS) with Procrustes analysis. NMDS is a multivariate ordination technique that is very useful to disclose multiple biotic relationships [17,18,82]. In NMDS the first two dimensions were retained for visualization, in which the axis with the highest variance was standardized between 0 and 1, and the other axis was rescaled according to the first one. Additionally, chi-square tests were used to test differences in the proportional distribution of geophilous and xylophilous species among vegetation types at different elevations.
To test Prediction 3, chi-square tests were used to test differences in the proportional distribution of main biogeographical groups among vegetation types. Chi-square tests were also used to investigate the association between ecological and biogeographical groups and association between endemics and vegetation types.
To test Prediction 4, tenebrionid assemblages were compared by constructing diversity profiles on the basis of the frequencies of chorotypes. To this end, we used Hill numbers qD [83,84,85,86,87]. Hill numbers are a mathematically unified family of diversity indices which differ among themselves only by the exponent q (which is the diversity order): as q can assume any value, it is possible to calculate, for a given assemblage of species, a diversity profile plotting qD versus q [85]. Although Hill numbers are typically used to express diversity with reference to the number of individuals belonging to the different species, they can be applied to any category (see, for example, their use for chemical diversity [87]), and in the case of the present research they were calculated using the number of species belonging to the different chorotypes, thus expressing the biogeographical diversity. Hill numbers are gaining increasing favor in community ecology [84,85,86], since they are all expressed in units of effective numbers, that is, the number of equally abundant categories (chorotypes in our case) required to obtain the same value of diversity. Therefore, they are perfectly comparable [84,86]. Hill numbers include three widely used diversity measures as special cases: richness (for q = 0), the exponential version of Shannon diversity (when q = 1 as the limiting case, since the equation for Hill numbers is undefined for q = 1) and Simpson diversity (q = 2) [84,86]. Since the parameter q controls the sensitivity of the measure to the relative abundances of the categories, for q = 0 the measure strongly reflects the contribution of rare categories; for q = 1, categories are weighed in proportion to their abundance; for q > 1, more weight is put on abundant categories, and the index can be roughly read as the effective number of the dominant categories in the assemblage, being more influenced by the most abundant categories [84,85,87]. We constructed diversity profiles by plotting qD versus q from q = 0 to q = 4 (as beyond q = 3 or 4 there is little variation, whereas values of diversity with q < 0 are not used, as they are dominated by the abundances of rare species and have poor statistical sampling properties) [85].
All calculations were performed in R 4.1.3 software [88]. Chi-square tests were performed with the stats package. Analyses based on Hill numbers were conducted using the package SpadeR version 0.1.1 [89]. β-diversity and NMDS analyses were performed using vegan 2.6–2 [90] and recluster 2.9 [91] packages. Figure 1a,b were constructed in ArcGis Pro 3.1.3 [92]. Figure 1b was realized with the function Hillshade using data from TINITALY [93].

3. Results

Species richness varied among vegetations, with the highest richness recorded in the Mediterranean shrubland and maquis vegetation, where about 45% of the total fauna occurred (Figure 3). In general, except for the first belt (represented by the sandy shore vegetation) species richness decreased with elevation, both for natural vegetation (Mediterranean shrubland and maquis: 38 species > Mixed forests: 31 species > Beech forests: 15 species > Pseudo-alpine grasslands: 4 species) and for secondary vegetation (Arid grasslands of lowlands and low hills: 28 species > Grasslands of medium and high hills: 11 species > Montane grasslands: 8 species), as expected according to Prediction 1.
The study fauna included much more geophilous species than xylophilous species (geophilous species: 55 species, 65%; xylophilous species: 30 species, 35%; χ2 = 7.353, df = 1, p = 0.007). The proportion of geophilous and xylophilous species varied among vegetation types (χ2 = 64.021, df = 7, p < 0.000001), with geophilous species being prevalent in coastal, secondary, and natural high elevation vegetations, and xylophilous species in natural forests (Figure 4). These results indicate that geophilous species are associated with open vegetations (both secondary and natural), as expected according to Prediction 2.
The study fauna included 14 species (about 16%) endemic to Italy. The remaining, non-endemic species (71 species) belong to 12 chorotypes (Table S1). If chorotypes are grouped into groups of broader meaning, it appears that vegetations of lower elevations were characterized by relatively high proportions of Mediterranean species (37–68%), whereas at higher elevations European and Palearctic species accounted for 70–100% of the species (χ2 = 58.460, df = 21, p < 0.00001) (Figure 5), thus supporting Prediction 3.
Diversity profiles based on chorotypes (Figure 6) indicated that: (1) grasslands of medium and high hills, montane grasslands, and pseudo-alpine grasslands are the vegetation types with the lowest chorotypological diversity for the whole range of q-values; (2) Mediterranean shrubland and maquis, arid grasslands of lowlands and hills, mixed forests, and beech forests are the vegetation types with the highest chorotypological diversity for the whole range of q-values; (3) dunes occupy an intermediate position. Overall, diversity profiles indicate no difference at low elevations between natural and secondary vegetations (Mediterranean shrubland and maquis, and arid grasslands of lowlands and low hills), while a clear difference exists in diversity for medium to high belts (between mixed and beech forests on one hand, and grasslands of medium and high hills on the other hand). These results support Prediction 4. A clear decline in diversity is seen in the open habitats independently of the naturalness (with the exception of dunes, that are less diversified than arid grasslands of lowlands and low hills), as expected according to Prediction 1.
Proportion of geophilous and xylophilous species varied among the main biogeographical groups (χ2 = 18.273, df = 3, p = 0.0004), with a disproportionate presence of geophilous species among endemics (93%, thus supporting Prediction 5) and Mediterranean species (81%); by contrast, xylophilous species were only slightly more abundant than geophilous ones among European (59%) and Palearctic species (58%).
Proportion of endemics varied among vegetation types (Figure 7, χ2 = 29.741, df = 7, p = 0.0001), with endemic species being highly associated with some open habitats such as montane grasslands (63%) and pseudo-alpine grasslands (75%) (Prediction 5).
Relationships between vegetation types outlined by overall similarities (βsor) in tenebrionid composition indicated the presence of two main groups (Figure 8a): one including forest vegetations (vegetations 2, 4, and 6) and the other including natural and secondary grasslands (vegetations 5, 7, and 8); dunes (vegetation 1) are well distinct and arid pastures of lowlands and low hills (vegetation 3) occupy an intermediate position. When only the pure turnover is considered (βsim, Figure 8b), arid pastures of lowlands and low hills appear to be more similar to vegetation 5 (grasslands of medium and high hills), which suggests that the higher distance obtained with βsor was due to the influence of nestedness on β-diversity pairwise relationships. In both ordinations, mixed forests (vegetation 4) occupy an intermediate position between Mediterranean shrubland and maquis (vegetation 2) and beech forests (vegetation 6), which indicates that the fauna found in mixed forests comprises species of both Mediterranean shrubland and maquis and beech forests. Specifically, about 23% of the species recorded in the mixed forests were also found in both the Mediterranean shrubland and maquis and the beech forests, about 45% were shared with the Mediterranean shrubland and maquis (but not with the beech forest), and about 16% with the beech forests (but not with the Mediterranean shrubland and maquis). Similarly, grasslands of medium and high hills (vegetation 5) showed an intermediate position between arid grasslands of lowlands and low hills (vegetation 3) and montane grasslands (vegetation 7). Specifically, about 27% of the species recorded in the grasslands of medium and high hills were also found in both the arid grasslands of lowlands and low hills and the montane grasslands, about 55% were shared with the arid grasslands of lowlands and low hills (but not with the montane grasslands), and about 18% with the montane grasslands (but not with the arid grasslands of lowlands and low hills). All species of the pseudo-alpine grasslands can be found in the montane grasslands. Thus, the communities of mixed forests share a large fraction (about 84%) of their species with those of the maquis and/or the beech forests, whereas those of the grasslands of medium and high hills share all of their species with those of the arid grasslands of lowlands and low hills and/or montane grasslands. These results support Prediction 2.

4. Discussion

Overall, the considered fauna represents about 25% of the entire Italian fauna (about 334 species/subspecies if Lagria and Alleculinae are excluded; estimate based on [73,74,94,95,96,97,98,99,100]). This high diversity can be related to the relatively large surface of the region, its position in the center of the Italian peninsula and the variety of landscapes [66,67,94].
The landscapes of Latium region (central Italy) have been profoundly modified by human presence, leading to a mosaic of natural, seminatural and artificial landscapes, with open areas of secondary vegetation interspersed within belts naturally characterized by forest vegetation [24,25,26,27]. Given the different preferences of animals for biotopes with different types of vegetation [68,69], species associated with open areas might infiltrate belts where forest is the expected vegetation if anthropogenic transformations produced open areas. Results obtained in this research are in accordance with the hypothesis that species richness declines with elevation in both natural and secondary vegetations because of environmental filtering processes, and this pattern is independent from species’ ecological preferences and biogeographical origin. Additionally, as hypothesized on the basis of species preferences for open and tree vegetations, species assemblages of open vegetations within elevational belts that should be naturally occupied by forests are not subsets of those of forest assemblages, but are more similar to those of natural open habitats occurring in other belts.
Consistent with results obtained in previous studies in which there was no distinction among vegetation types [38,39,40], the current study shows that, in general, tenebrionid species richness decreases upwards, in both natural and secondary vegetation (Prediction 1), which suggests that the drivers of such elevational decline in tenebrionid diversity (possibly represented by increasing harsher conditions) are more important and likely different from changes in vegetation. However, it is important to note that the lowest belt, represented by sandy shore vegetation, has fewer species than the following belt, occupied by the Mediterranean maquis. This is clearly a reflection of the harsh conditions of coastal environments (e.g., high temperatures, especially on the ground in summer, low water availability, low productivity, low vegetation diversity), where tenebrionid communities are mainly represented by few, highly adapted species [101,102]. Although this decline in diversity might be influenced by the fact that, with increasing elevation, available area also decreases, previous research indicated that decrease in species richness of tenebrionids along the elevational gradient of Latium can be observed also after controlling for available area [39,40].
As expected, ecological preferences influenced species association with vegetation (Prediction 2): geophilous (ground-dwelling) and xylophilous (arboreal) species showed different proportions among vegetation types, with geophilous species prevailing in open vegetations, and xylophilous species in natural forests, independently from their elevation. Despite the higher proportion of the geophilous species in the study fauna, secondary open vegetations have fewer species than the corresponding natural (forest) vegetation found at the same elevations. These results stress the conservation importance of forests for tenebrionid beetles, especially for saproxylic species, some of which are imperiled [71,103]. Replacement of forests with secondary open habitats would impact dramatically on the diversity of tenebrionids at elevations that would be naturally occupied by forests. Additionally, the tenebrionids inhabiting secondary (open) vegetations are not subsets of those found in the corresponding natural (forest) vegetation, as shown by their low similarity in species composition (Prediction 2). This indicates that the assemblages of secondary vegetations are not impoverished versions of those associated with the original forest vegetation. If, on the one hand, this implies that secondary vegetation types cannot retain original communities of elevational belts naturally occupied by forests, on the other hand it suggests that they might serve as a habitat extension for some species found in natural vegetation at a different elevations, which supports the importance of secondary vegetation for certain species [68,69]. Moreover, endemics represent a conspicuous proportion (more than 50% of the species) of the tenebrionid assemblages of high-elevation open habitats, included those occupied by secondary vegetation. In general, the study area is home to a large number of endemics (at least 14 taxa), which account for about 10% of the total number of species/subspecies endemic to Italy (about 137 species/subspecies; estimate based on [74,94,95,96,97,98,99,100]). These results further support the importance of secondary (seminatural) vegetations for biodiversity conservation in the Apennines [68,69]. Interestingly, endemics occur almost exclusively among geophilous species, which can be explained by the low dispersal capabilities of most ground-dwelling tenebrionids, which are flightless. In general, the study area shows a predominance of geophilous species, which reflects both their prevalence in the whole Italian fauna (more than 75%) [41] and the prevalence of open areas in the study area (with less of 40% of the surface occupied by forests) [104].
Biogeographically, the study fauna is characterized by a relatively high proportion of species mainly distributed in the Mediterranean basin. This obviously reflects the geographical position of the study area, in the middle of the Mediterranean basin. However, Mediterranean species represent a conspicuous fraction (~40–70%) of the species assemblages in the three lowest vegetations (Prediction 3), whereas they are a minority of species in the mixed forests, in the grasslands of medium and high hills, and in the beech forests, being finally absent in the montane and pseudo-alpine grasslands. This reduction in the Mediterranean species is associated with an increase in European and Palearctic species. Differently from Mediterranean species, which are typically thermophilic species, these non-Mediterranean species are associated with colder climates, and their presence in the study area (in many cases limited to middle-high elevations) increases its biogeographical diversity, especially at low and intermediate elevations (Prediction 4). This is especially illustrated by the high biogeographical diversity (diversity profiles based on Hill numbers) shown by mixed forests and beech forests. Although arid grasslands of lowlands and hills are an anthropogenic vegetation, they host a biogeographical diversity even higher than that of Mediterranean shrubland and maquis, when the role of rare chorotypes is emphasized (Hill numbers with low q-values), thus highlighting the importance of secondary vegetation for these insects. Overall, vegetations at intermediate elevations (beech forests, mixed forests, and arid grasslands of lowlands and hills) showed higher biogeographical diversity values, compared to those of higher (montane and pseudo-alpine grasslands) and lower (dunes and possibly Mediterranean shrubland and maquis) elevations (Prediction 4).
We found a strong association between chorotypes and lifestyle, which explain the high proportion of endemic species at high elevations. Geophilous species prevail among Mediterranean and endemic species, whereas xylophilous species are mainly associated with European and Palearctic distributions. As high elevations are dominated by open vegetation, this leads to a predominance of geophilous species, and hence favors a higher proportion of endemics (Prediction 5). By contrast, European and Palearctic species are mainly associated with tree vegetations, especially those that predominate in areas with colder climates, being therefore prevalent in the mixed and beech forests. These results further stress the importance of forest vegetation for tenebrionids.
Finally, as expected according to Prediction 2, tenebrionid assemblages of open vegetations, both natural and secondary, appear to be relatively similar in species composition, being distinctly separated from that of forest vegetation, both in terms of overall composition (Dice–Sørensen coefficient) and pure turnover (Simpson coefficient). In particular: (1) dunes are profoundly different in tenebrionid composition from all other vegetations; (2) Mediterranean shrubland, mixed forests and beech forests clustered apart; (3) grasslands of middle elevation, montane grasslands and pseudo-alpine grasslands form another group; (4) arid grasslands of low elevations show an intermediate position between the Mediterranean shrubland and the grasslands of middle elevations. These results demonstrate the importance of vegetation types in constraining tenebrionid distributions across the elevational gradient.

5. Conclusions

Grasslands of secondary vegetation may represent important habitats for tenebrionid beetles in central Italy. They are characterized by high levels of species richness, biogeographical diversity and endemism, but are not reservoirs of species that are associated with the forest vegetation that would naturally occur at certain elevations. Tenebrionid assemblages in these habitats are not impoverished versions of those recorded from the forests naturally occurring at the same elevations, but form very different species assortments. Because of increasing human pressures on these seminatural habitats, it is important that their importance is recognized and valued in conservation programs, to ensure they are not subject to further degradation processes. However, this should not divert conservation efforts from protecting natural vegetation. If, on one hand, secondary vegetations can play a pivotal role for biological conservation in the Mediterranean biodiversity hotspots, it is also important that natural vegetations are adequately preserved, especially for the conservation of xylophilous species.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d16020110/s1, Table S1: Tenebrionid species distribution (0 = absence, 1 = presence) in different vegetation types in Latium and their ecological and biogeographical characterization. Chorotypes: CEU: Centraleuropean, END: Endemic to Italy; EUR: European, SEU: S-European, MED: Mediterranean; NAF: N-African; WME: W-Mediterranean; EME: E-Mediterranean; CEM: Centralasiatic-Europeo-Mediterranean; OLA: Holarctic; SIE: Sibero-European; TEM: Turano-Europeo-Mediterranean; TUM: Turano-Mediterranean.

Funding

This research was partially funded by the University of L’Aquila, Department of Life, Health and Environmental Sciences, FFO 2023 “Struttura di comunità in ecosistemi di sky island”.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data supporting reported results are provided as Table S1.

Acknowledgments

I am grateful to Maria Laura Carranza and Letizia Di Biase for discussions about vegetation types, to Davide Bergamaschi for his help with Figure 1, and to three anonymous reviewers for their comments on a previous version of this paper. Thanks are due to Letizia Di Biase and Andrea Di Giulio for allowing me to use their photos in Figure 2.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Poschlod, P.; WallisDeVries, M.F. The historical and socioeconomic perspective of calcareous grasslands—Lessons from the distant and recent past. Biol. Conserv. 2002, 104, 361–376. [Google Scholar] [CrossRef]
  2. Hejcman, M.; Hejcmanová, P.; Pavlů, V.; Beneš, J. Origin and history of grasslands in Central Europe—A review. Grass Forage Sci. 2013, 68, 345–363. [Google Scholar] [CrossRef]
  3. Giarrizzo, E.; Burrascano, S.; Chiti, T.; de Bello, F.; Lepš, J.; Zavattero, L.; Blasi, C. Re-visiting historical semi-natural grasslands in the Apennines to assess patterns of changes in species composition and functional traits. Appl. Veg. Sci. 2017, 20, 247–258. [Google Scholar] [CrossRef]
  4. Hobohm, C.; Bruchmann, I. Endemische Gefäßpflanzen und ihre Habitate in Europa: Plädoyer für den Schutz der Grasland-Ökosysteme. Berichte Reinhold-Tüxen-Ges. 2009, 21, 142–161. [Google Scholar]
  5. WallisDeVries, M.F.; van Swaay, C.A.M. Grasslands as habitats for butterflies in Europe. In Grasslands in Europe of High Nature Value; Veen, P., Jefferson, R., de Smidt, J., van der Straaten, J., Eds.; KNNV Publishing: Zeist, The Netherlands, 2009; pp. 27–34. [Google Scholar]
  6. Wilson, J.B.; Peet, R.K.; Dengler, J.; Pärtel, M. Plant Species Richness: The World Records. J. Veg. Sci. 2012, 23, 796–802. [Google Scholar] [CrossRef]
  7. Habel, J.C.; Dengler, J.; Janišová, M.; Török, P.; Wellstein, C.; Wiezik, M. European grassland ecosystems: Threatened hotspots of biodiversity. Biodivers. Conserv. 2013, 22, 2131–2138. [Google Scholar] [CrossRef]
  8. Öckinger, E.; Eriksson, A.K.; Smith, H.G. Effects of grassland abandonment, restoration and management on butterflies and vascular plants. Biol. Conserv. 2006, 133, 291–300. [Google Scholar] [CrossRef]
  9. Assini, S.; Filipponi, F.; Zucca, F. Land cover changes in an abandoned agricultural land in the Northern Apennine (Italy) between 1954 and 2008: Spatio-temporal dynamics. Plant Biosyst. 2014, 149, 807–817. [Google Scholar] [CrossRef]
  10. Bonari, G.; Fajmon, K.; Malenovský, I.; Zelený, D.; Holuša, J.; Jongepierová, I.; Kočárek, P.; Konvička, O.; Uřičář, J.; Chytrý, M. Management of semi-natural grasslands benefiting both plant and insect diversity: The importance of heterogeneity and tradition. Agric. Ecosyst. Environ. 2017, 246, 243–252. [Google Scholar] [CrossRef]
  11. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.B.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  12. Médail, F.; Myers, N. Mediterranean Basin. In Hotspots Revisited: Earth’s Biologically Richest and Most Endangered Terrestrial Ecoregions; Mittermeier, R.A., Gil, P.R., Hoffman, M., Pilgrim, J., Brooks, T., Mittermeier, C.G., Lamoreux, J., Da Fonseca, G.A.B., Eds.; CEMEX: Agrupación Sierra Madre, Mexico, 2004; pp. 144–147. [Google Scholar]
  13. Mittermeier, R.A.; Turner, W.R.; Larsen, F.W.; Brooks, T.M.; Gascon, C. Global biodiversity conservation: The critical role of hotspots. In Biodiversity Hotspots: Distribution and Protection of Conservation Priority Areas; Zachos, F.E., Habel, J.C., Eds.; Springer: Heidelberg, Germany, 2011; pp. 2–22. [Google Scholar]
  14. Birdlife International. Ecosystem Profile. Mediterranean Basin Biodiversity Hotspot. Available online: https://www.cepf.net/sites/default/files/mediterranean-basin-2017-ecosystem-profile-english_0.pdf (accessed on 25 November 2023).
  15. Fattorini, S.; Ulrich, W. Spatial distributions of European Tenebrionidae point to multiple postglacial colonization trajectories. Biol. J. Linn. Soc. 2012, 105, 318–329. [Google Scholar] [CrossRef]
  16. Schmitt, T.; Fritz, U.; Delfino, M.; Ulrich, W.; Habel, J.C. Biogeography of Italy revisited: Genetic lineages confirm major phylogeographic patterns and a pre-Pleistocene origin of its biota. Front. Zool. 2021, 18, 34. [Google Scholar] [CrossRef] [PubMed]
  17. Fattorini, S. Odonate diversity patterns in Italy disclose intricate colonization pathways. Biology 2022, 11, 886. [Google Scholar] [CrossRef] [PubMed]
  18. Fattorini, S. Biogeographical patterns of earwigs in Italy. Insects 2023, 14, 235. [Google Scholar] [CrossRef] [PubMed]
  19. Stoch, F. How many endemic species? Species richness assessment and conservation priorities in Italy. Belg. J. Entomol. 2000, 2, 125–133. [Google Scholar]
  20. Ruffo, S.; Stoch, F. (Eds.) Checklist and Distribution of the Italian Fauna. 10,000 Terrestrial and Freshwater Species. 2 Serie, Sez. Scienze della Vita, 2nd ed.; Memorie del Museo Civico di Storia Naturale di Verona: Verona, Italy, 2007; Volume 17. [Google Scholar]
  21. Blasi, C.; Biondi, E. La Flora in Italia; Ministero dell’Ambiente e della Tutela del Territorio e del Mare, Sapienza Università Editrice: Rome, Italy, 2017; pp. 1–704. [Google Scholar]
  22. Blasi, C.; Boitani, L.; La Posta, S.; Manes, F.; Marchetti, M. Stato della Biodiversità in Italia; Palombi Editori: Roma, Italy, 2005; pp. 1–468. [Google Scholar]
  23. Stinca, A.; Musarella, C.M.; Rosati, L.; Laface, V.L.A.; Licht, W.; Fanfarillo, E.; Wagensommer, R.P.; Galasso, G.; Fascetti, S.; Esposito, A.; et al. Italian vascular flora: New findings, updates and exploration of floristic similarities between regions. Diversity 2021, 13, 600. [Google Scholar] [CrossRef]
  24. Sestini, A. Conosci l’Italia. Volume VII. Il Paesaggio; Touring Club Italiano: Milano, Italy, 1963; pp. 1–232. [Google Scholar]
  25. Bonapace, U. (Ed.) I Paesaggi Umani; Touring Club Italiano: Milano, Italy, 1977; pp. 1–224. [Google Scholar]
  26. Pignatti, S. Ecologia del Paesaggio; UTET: Torino, Italy, 1994; pp. 1–228. [Google Scholar]
  27. Touring Club Italiano. Il Paesaggio Italiano: Idee Contributi Immagini; Touring Club Italiano: Milano, Italy, 2000; pp. 1–224. [Google Scholar]
  28. ISTAT. Annuario Statistico Italiano; Istituto Nazionale di Statistica: Roma, Italy, 2023; pp. 1–857. Available online: https://www.istat.it/storage/ASI/2022/ASI_2022.pdf (accessed on 10 December 2023).
  29. Caveri, L.; Da Polenza, A.; Salvatori, F.; Smiraglia, C. (Eds.) Montagne d’Italia; Touring Club Italiano: Milano, Italy, 2002; pp. 1–320. [Google Scholar]
  30. 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]
  31. 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]
  32. Ortisi, A. Boschi e legno in Europa. In Storia dell’Agricoltura in Europa; ETAS-BNL: Roma, Italy, 1980; pp. 64–115. [Google Scholar]
  33. Kelly, M.G.; Huntley, B. An 11,000-year record of vegetation and environment from Lago di Martignano, Latium, Italy. J. Quat. Sci. 1991, 6, 209–224. [Google Scholar] [CrossRef]
  34. Sadori, L.; Giardini, M.; Giraudi, C.; Mazzini, I. The plant landscape of the imperial harbour of Rome. J. Archaeol. Sci. 2010, 37, 3294–3305. [Google Scholar] [CrossRef]
  35. Sadori, L.; Mercuri, A.M.; Mariotti Lippi, M. Reconstructing past cultural landscape and human impact using pollen and plant macroremains. Plant Biosyst. 2010, 144, 940–951. [Google Scholar] [CrossRef]
  36. Fanfarillo, E.; Latini, M.; Nicolella, G.; Abbate, G. Development of a new GIS-based method to detect high natural value farmlands: A case study in Central Italy. Ann. Bot. 2017, 7, 25–31. [Google Scholar]
  37. Mensing, S.A.; Schoolman, E.M.; Tunno, I.; Noble, P.J.; Sagnotti, L.; Florindo, F.; Piovesan, G. Historical ecology reveals landscape transformation coincident with cultural development in central Italy since the Roman Period. Sci. Rep. 2018, 8, 2138. [Google Scholar] [CrossRef]
  38. Fattorini, S. Variation in zoogeographical composition along an elevational gradient: The tenebrionid beetles of Latium (Central Italy). Entomologia 2013, 1, 33–40. [Google Scholar] [CrossRef]
  39. 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]
  40. 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]
  41. Fattorini, S.; Mantoni, C.; Audisio, P.; Biondi, M. Taxonomic variation in levels of endemism: A case study of Italian tenebrionid beetles. Insect. Conserv. Divers. 2019, 12, 351–361. [Google Scholar] [CrossRef]
  42. Lucchese, F. Atlante della Flora Vascolare del Lazio, Cartografia, Ecologia e Biogeografia. Vol. 1. Parte Generale e Flora Alloctona; Regione Lazio, Direzione Capitale Naturale, Parchi e Aree Protette: Rome, Italy, 2017; pp. 1–356. [Google Scholar]
  43. Lucchese, F. Atlante della Flora Vascolare del Lazio, Cartografia, Ecologia e Biogeografia. Vol. 2. La Flora di Maggiore Interesse Conservazionistico; Regione Lazio, Direzione Capitale Naturale, Parchi e Aree Protette: Rome, Italy, 2018; p. 400. [Google Scholar]
  44. Blasi, C. Fitoclimatologia del Lazio; Università “La Sapienza”-Roma, Regione Lazio: Roma, Italy, 1994; pp. 1–58. [Google Scholar]
  45. Salvati, L.; Perini, L.; Bajocco, S.; Sabbi, A. Climate aridity and landuse change: A regional-scale analysis. Geogr. Res. 2012, 50, 193–203. [Google Scholar] [CrossRef]
  46. Fenaroli, L.; Giacomini, V. Conosci l’Italia. Vol. II. La Flora; Touring Club Italiano: Milano, Italy, 1958; pp. 1–276. [Google Scholar]
  47. Fenaroli, L. Note Illustrative della Carta della Vegetazione Reale d’Italia; Ministero dell’Agricoltura e delle Foreste: Roma, Italy, 1970; pp. 1–128. [Google Scholar]
  48. Tomaselli, R. Note Illustrative della Carta della Vegetazione Naturale Potenziale d’Italia; Ministero dell’Agricoltura e delle Foreste: Roma, Italy, 1970; pp. 1–64. [Google Scholar]
  49. Pignatti, S. I piani di vegetazione in Italia. Giorn. Bot. Ital. 1979, 113, 411–428. [Google Scholar] [CrossRef]
  50. Ubaldi, D. Le fasce della vegetazione italiana su basi fitosociologiche. Giorn. Bot. ital. 1989, 123 (Suppl. 1), 106. [Google Scholar]
  51. Carpaneto, G.M. Note sulla distribuzione geografica ed ecologica dei Coleotteri Scarabaeoidea Laparosticti nell’Italia Appenninica (I Contributo). Boll. Ass. Romana. Entomol. 1986, 40, 31–58. [Google Scholar]
  52. Stanisci, A.; Acosta, A.; Ercole, S.; Blasi, C. Plant communities on coastal dunes in Lazio (Italy). Ann. Bot. 2004, 2, 39–48. [Google Scholar]
  53. Prisco, I.; Carboni, M.; Acosta, A.T.R. The fate of threatened coastal dune habitat in Italy under climate change scenarios. PLoS ONE 2013, 8, e68850. [Google Scholar] [CrossRef] [PubMed]
  54. De Felici, S.; Vigna Taglianti, A. I Coleotteri Carabidi dei Colli Albani (Coleoptera, Carabidae). Boll. Ass. Romana Entomol. 1994, 49, 3–96. [Google Scholar]
  55. Angelini, A.; Mattioli, W.; Merlini, P.; Corona, P.; Portoghesi, L. Empirical modelling of chestnut coppice yield for Cimini and Vicani mountains (Central Italy). Ann. Silvic. Res. 2013, 37, 7–12. [Google Scholar]
  56. Ziaco, E.; Alessandrini, A.; Blasi, S.; Di Filippo, A.; Dennis, S.; Piovesan, G. Communicating old-growth forest through an educational trail. Biodivers. Conserv. 2012, 21, 131–144. [Google Scholar] [CrossRef]
  57. Piovesan, G.; Alessandrini, A.; Baliva, M.; Chiti, T.; D’Andrea, E.; De Cinti, B.; Di Filippo, A.; Hermanin, L.; Lauteri, M.; Scarascia Mugnozza, G.; et al. Structural patterns, growth processes, carbon stocks in an Italian network of old-growth beech forests. Ital. J. For. Mt. Environ. 2010, 65, 557–590. [Google Scholar] [CrossRef]
  58. Direzione Regionale Capitale Naturale, Parchi e Aree Protette. Carta delle formazioni naturali e seminaturali al IV e V livello Corine Land Cover. Available online: http://dati.lazio.it/catalog/dataset/cus-lazio-approfondimento-delle-formazioni-naturali-e-seminaturali-iv-e-v-livello-corine-land-cover (accessed on 25 January 2024).
  59. Trichas, A. The genus Dendarus Latreille, 1829 (Coleoptera, Tenebrionidae: Dendarini) in Greece (A systematic account of the genus with description of a new species and four new systematic combinations). In Advances in Arachnology and Developmental Biology; Makarov, S.E., Dimitrijević, R.N., Eds.; SASA, Belgrade and UNESCO MAB Serbia: Belgrade, Serbia, 2008; pp. 417–462. [Google Scholar]
  60. Ferrer, J. Contribución al conocimiento de los Asidini iberobaleares. Segunda nota. Las Alphasida (Glabrasida) del grupo Tricostatae Escalera 1922 (Coleoptera, Tenebrionidae, Pimeliinae). Boletín Soc. Entomológica Aragonesa A 2008, 43, 61–73. [Google Scholar]
  61. Condamine, F.L.; Soldati, L.; Rasplus, J.Y.; Kergoat, G.J. New insights on systematics and phylogenetics of Mediterranean Blaps species (Coleoptera: Tenebrionidae: Blaptini), assessed through morphology and dense taxon sampling. Syst. Entomol. 2011, 36, 340–361. [Google Scholar] [CrossRef]
  62. Ferrer, J. Contribución al conocimiento del Género Phylan Dejean, 1821, y descripción de una specie nueva del género Heliopates Dejean, 1834 (Coleoptera, Tenebrionidae, Pedinini). Boletín Soc. Entomológica Aragonesa 2011, 49, 75–82. [Google Scholar]
  63. Pons, J.; Bruvo, B.; Petitpierre, E.; Plohl, M.; Ugarkovic, D.; Juan, C. Complex structural features of satellite DNA sequences in the genus Pimelia (Coleoptera: Tenebrionidae): Random differential amplification from a common ‘satellite DNA library’. Heredity 2004, 9, 418–427. [Google Scholar] [CrossRef]
  64. Soldati, F.; Soldati, L. Species delimitation using morphological and molecular tools in the Asida (Polasida) jurinei Solier, 1836 species complex. Preliminary results. (Coleoptera: Tenebrionidae: Tentyrinae). Cah. Sci. Muséum Lyon 2006, 10, 111–116. [Google Scholar]
  65. Stroscio, S.; Baviera, C.; Frati, F.; Lo Paro, G.; Nardi, F. Deep genetic divergence in the darkling beetle Pimelia rugulosa (Coleoptera, Tenebrionidae) reflects Plio-Pleistocenic paleogeographic history of Sicily. J. Zool. Syst. Evol. Res. 2011, 49, 196–203. [Google Scholar] [CrossRef]
  66. Fattorini, S. Regional insect inventories require long time, extensive spatial sampling and good will. PLoS ONE 2013, 8, e62118. [Google Scholar] [CrossRef]
  67. Fattorini, S.; Salvati, L. Tenebrionid beetles as proxy indicators of climate aridity in a Mediterranean area. Ecol. Indic. 2014, 38, 256–261. [Google Scholar] [CrossRef]
  68. Fattorini, S. Use of insect rarity for biotope prioritisation: The tenebrionid beetles of the Central Apennines (Italy). J. Insect Conserv. 2010, 14, 367–378. [Google Scholar] [CrossRef]
  69. Fattorini, S. Biotope prioritisation in the Central Apennines (Italy): Species rarity and cross-taxon congruence. Biodivers. Conserv. 2010, 19, 3413–3429. [Google Scholar] [CrossRef]
  70. Fattorini, S. I Coleotteri Tenebrionidi di Roma (Coleoptera Tenebrionidae). Fragm. Entomol. 2013, 45, 87–142. [Google Scholar] [CrossRef]
  71. Fattorini, S.; Maltzeff, P.; Salvati, L. Use of insect distribution across landscape-soil units to assess conservation priorities in a Mediterranean coastal reserve: The tenebrionid beetles of Castelporziano (Central Italy). Rend. Lincei-Sci. Fis. 2015, 26, 353–366. [Google Scholar] [CrossRef]
  72. Iwan, D.; Löbl, I. Catalogue of Palaearctic Coleoptera. Vol. 5. Revised and Updated Second Edition. Tenebrionoidea; Brill: Leiden, The Netherlands, 2020; p. 946. [Google Scholar]
  73. Ferrer, J. Revisión del género Leptoderis Billberg, 1820 y comentarios sobre el origen, composición, anatomía y necrofagia de la tribu Elenophorini (Coleoptera, Tenebrionidae). Boln. SEA 2015, 57, 19–38. [Google Scholar]
  74. Ferrer, J. Contribución al conocimiento del género Pachychila Eschscholtz, 1831 y descripción de dos especies nuevas de la Península Ibérica y una de Italia (Coleoptera, Tenebrionidae, Pimeliinae). Boln. SEA 2018, 62, 35–54. [Google Scholar]
  75. Fattorini, S. On the concept of chorotype. J. Biogeogr. 2015, 42, 2246–2251. [Google Scholar] [CrossRef]
  76. Fattorini, S. A history of chorological categories. Hist. Philos. Life Sci. 2016, 38, 12. [Google Scholar] [CrossRef] [PubMed]
  77. Gatto, C.A.F.R.; Cohn-Haft, M. Spatial Congruence Analysis (SCAN): A method for detecting biogeographical patterns based on species range congruences. PLoS ONE 2021, 16, e0245818. [Google Scholar] [CrossRef] [PubMed]
  78. Olivero, J.; Real, R.; Márquez, A.L. Fuzzy chorotypes as a conceptual tool to improve insight into biogeographic patterns. Syst. Biol. 2011, 60, 645–660. [Google Scholar] [CrossRef] [PubMed]
  79. Vigna Taglianti, A.; Audisio, P.A.; Biondi, M.; Bologna, M.A.; Carpaneto, G.M.; De Biase, A.; Fattorini, S.; Piattella, E.; Sindaco, R.; Venchi, A.; et al. A proposal for a chorotype classification of the Near East fauna, in the framework of the Western Palearctic region. Biogeographia 1999, 20, 31–59. [Google Scholar] [CrossRef]
  80. Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 2010, 19, 134–143. [Google Scholar] [CrossRef]
  81. Baselga, A.; Leprieur, F. Comparing methods to separate components of beta diversity. Methods Ecol. Evol. 2015, 6, 1069–1079. [Google Scholar] [CrossRef]
  82. Moulpied, M.; Smith, C.H.; Robertson, C.R.; Johnson, N.A.; Lopez, R.; Randklev, C.R. Biogeography of freshwater mussels (Bivalvia: Unionida) in Texas and implications on conservation biology. Divers. Distrib. 2022, 28, 1458–1474. [Google Scholar] [CrossRef]
  83. Jost, L. The relation between evenness and diversity. Diversity 2010, 2, 207–232. [Google Scholar] [CrossRef]
  84. 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]
  85. 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]
  86. 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]
  87. Petrén, H.; Köllner, T.G.; Junker, R.R. Quantifying chemodiversity considering biochemical and structural properties of compounds with the R package CHEMODIV. New Phytol. 2023, 237, 2478–2492. [Google Scholar] [CrossRef] [PubMed]
  88. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 15 March 2022).
  89. 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).
  90. Oksanen, J.; Simpson, G.L.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Solymos, P.; Stevens, M.H.H.; Szoecs, E. Vegan: Community Ecology Package. Available online: https://CRAN.R-project.org/package=vegan (accessed on 20 April 2022).
  91. Dapporto, L.; Ramazzotti, M.; Fattorini, S.; Vila, R.; Talavera, G.; Dennis, R.L.H. recluster: Ordination Methods for the Analysis of Beta-Diversity Indices. Available online: https://CRAN.R-project.org/package=recluster (accessed on 15 March 2022).
  92. ESRI. ArcGis Desktop: Version 10.4; Environmental Systems Research Institute (ESRI): Redlands, CA, USA, 2016. [Google Scholar]
  93. Tarquini, S.; Isola, I.; Favalli, M.; Battistini, A.; Dotta, G. TINITALY, a Digital Elevation Model of Italy with a 10 Meters Cell Size, Version 1.1. Istituto Nazionale di Geofisica e Vulcanologia (INGV). 2023. Available online: https://tinitaly.pi.ingv.it/ (accessed on 3 October 2023).
  94. Fattorini, S. Tenebrionid beetle distributional patterns in Italy: Multiple colonisation trajectories in a biogeographical crossroad. Insect Conserv. Divers. 2014, 7, 144–160. [Google Scholar] [CrossRef]
  95. Aliquò, V.; Soldati, F. Updating the CD-rom on Coleoptera Tenebrionidae of Italy and the check-list of the same family. Biodivers. J. 2014, 5, 429–442. [Google Scholar]
  96. Leo, P.; Fancello, L. Observaciones sobre Akis bacarozzo (Schrank, 1786) y Akis tuberculata Kraatz, 1865, y nota sinonímica (Coleoptera, Tenebrionidae). Rev. Gaditana De Entomol. 2019, 10, 31–46. [Google Scholar]
  97. Leo, P.; Fancello, L. Clitobius subplumbeus (Fairmaire, 1856) bona species: Nuovo inquadramento sistematico delle popolazioni mediterranee attribuite a Clitobius ovatus (Erichson, 1843) (Coleoptera, Tenebrionidae). Ann. Mus. Civ. Stor. Nat. “Giacomo Doria” 2020, 112, 155–166. [Google Scholar]
  98. Ponel, P.; Lo Cascio, P.; Leo, P. A new species of Catomus Allard, 1876 (Coleoptera: Tenebrionidae: Helopini) from the Aeolian archipelago (Sicily, Italy). Zootaxa 2020, 4743, 2. [Google Scholar] [CrossRef] [PubMed]
  99. Kamiński, M.J.; Bouchard, P.; Leo, P.; Iwan, D. On the type species of the genus Dendarus Dejean, 1821 (Coleoptera: Tenebrionidae). Zootaxa 2021, 5048, 575–580. [Google Scholar] [CrossRef]
  100. Leo, P.; Pezzi, G.; Riolo, F. Prima segnalazione per l’Italia di Cryphaeus cornutus (Fischer von Waldheim, 1823) e della tribù Toxicini con osservazioni sulla distribuzione della specie in Europa (Coleoptera Tenebrionidae Toxicini). Boll. Soc. Entomol. Ital. 2021, 153, 123–128. [Google Scholar] [CrossRef]
  101. Fattorini, S. Ecology and conservation of tenebrionid beetles in Mediterranean coastal areas. In Insect Ecology and Conservation; Fattorini, S., Ed.; Research Signpost: Trivandrum, Kerala, India, 2008; pp. 165–297. [Google Scholar]
  102. Fattorini, S. Adaptations of tenebrionid beetles to Mediterranean sand dune environments and the impact of climate change (Coleoptera: Tenebrionidae). Fragm. Entomol. 2023, 55, 1–20. [Google Scholar]
  103. Carpaneto, G.M.; Baviera, C.; Biscaccianti, A.B.; Brandmayr, P.; Mazzei, A.; Mason, F.; Battistoni, A.; Teofili, C.; Rondinini, C.; Fattorini, S.; et al. A Red List of Italian Saproxylic Beetles: Taxonomic overview, ecological features and conservation issues (Coleoptera). Fragm. Entomol. 2015, 47, 53–126. [Google Scholar] [CrossRef]
  104. ISTAT. Rapporto sul Territorio 2020. Ambiente, Economia e Società; Istituto Nazionale di Statistica: Roma, Italy, 2020; pp. 1–201. Available online: https://www.istat.it/storage/rapporti-tematici/territorio2020/Rapportoterritorio2020.pdf (accessed on 10 December 2023).
Figure 1. Study area (Latium, central Italy). (a): Distribution of main natural and secondary vegetation types. The inset shows the location of Latium (in red) within the Italian territory. (b): Distribution of tenebrionid records (each dot is a locality for which one or more records are available).
Figure 1. Study area (Latium, central Italy). (a): Distribution of main natural and secondary vegetation types. The inset shows the location of Latium (in red) within the Italian territory. (b): Distribution of tenebrionid records (each dot is a locality for which one or more records are available).
Diversity 16 00110 g001
Figure 2. Examples of natural and secondary vegetations occurring along the elevational gradient in the Apennines of central Italy. (a): Dunes; (b): Mediterranean shrubland and maquis; (c): Arid grasslands of lowlands and low hills; (d): Mixed forests; (e): Grasslands of medium and high hills; (f): Beech forests; (g): Montane grasslands; (h): Pseudo-alpine grasslands. Photos (a,b,f) by S. Fattorini; photos (ce,h) by L. Di Biase; photo (g) by A. Di Giulio.
Figure 2. Examples of natural and secondary vegetations occurring along the elevational gradient in the Apennines of central Italy. (a): Dunes; (b): Mediterranean shrubland and maquis; (c): Arid grasslands of lowlands and low hills; (d): Mixed forests; (e): Grasslands of medium and high hills; (f): Beech forests; (g): Montane grasslands; (h): Pseudo-alpine grasslands. Photos (a,b,f) by S. Fattorini; photos (ce,h) by L. Di Biase; photo (g) by A. Di Giulio.
Diversity 16 00110 g002
Figure 3. Number of tenebrionid species recorded from different vegetations along the elevational gradient in Latium (central Italy). Grid patterns identify secondary vegetations.
Figure 3. Number of tenebrionid species recorded from different vegetations along the elevational gradient in Latium (central Italy). Grid patterns identify secondary vegetations.
Diversity 16 00110 g003
Figure 4. Percentage of geophilous and xylophilous tenebrionid species in different vegetations along the elevational gradient in Latium (central Italy). Grid patterns identify secondary vegetations.
Figure 4. Percentage of geophilous and xylophilous tenebrionid species in different vegetations along the elevational gradient in Latium (central Italy). Grid patterns identify secondary vegetations.
Diversity 16 00110 g004
Figure 5. Biogeographical composition of tenebrionid species in different vegetations in Latium (Central Italy). END: Endemic; EUR: European, PAL: Palearctic; MED: Mediterranean. Grid patterns identify secondary vegetations.
Figure 5. Biogeographical composition of tenebrionid species in different vegetations in Latium (Central Italy). END: Endemic; EUR: European, PAL: Palearctic; MED: Mediterranean. Grid patterns identify secondary vegetations.
Diversity 16 00110 g005
Figure 6. Biogeographical diversity profiles for tenebrionid species in different vegetations along the elevational gradient in Latium (central Italy). Black lines indicate natural vegetations, grey lines indicate secondary vegetations.
Figure 6. Biogeographical diversity profiles for tenebrionid species in different vegetations along the elevational gradient in Latium (central Italy). Black lines indicate natural vegetations, grey lines indicate secondary vegetations.
Diversity 16 00110 g006
Figure 7. Percentages of endemic (END) and non-endemic (Non-END) tenebrionid species from different vegetations in Latium (central Italy). Grid patterns identify secondary vegetations.
Figure 7. Percentages of endemic (END) and non-endemic (Non-END) tenebrionid species from different vegetations in Latium (central Italy). Grid patterns identify secondary vegetations.
Diversity 16 00110 g007
Figure 8. Non-metric multidimensional scaling plots showing β-diversity distances between vegetation types for tenebrionid species composition along the elevational gradient in Latium (central Italy). (a) Results obtained using Dice–Sørensen coefficient (Stress: <0.0001); (b) Results obtained using Simpson coefficient (Stress: 0.017.). Numbers identify vegetation types as follows: 1. Dunes; 2. Mediterranean shrubland and maquis; 3. Arid grasslands of lowlands and low hills; 4. Mixed forests (Turkey oaks and chestnuts); 5. Grasslands of medium and high hills; 6. Beech forests; 7. Montane grasslands; 8. Pseudo-alpine grasslands (beyond the tree line). Black symbols: natural vegetations; white symbols: secondary vegetations. Squares: open vegetations; circles: forest vegetations.
Figure 8. Non-metric multidimensional scaling plots showing β-diversity distances between vegetation types for tenebrionid species composition along the elevational gradient in Latium (central Italy). (a) Results obtained using Dice–Sørensen coefficient (Stress: <0.0001); (b) Results obtained using Simpson coefficient (Stress: 0.017.). Numbers identify vegetation types as follows: 1. Dunes; 2. Mediterranean shrubland and maquis; 3. Arid grasslands of lowlands and low hills; 4. Mixed forests (Turkey oaks and chestnuts); 5. Grasslands of medium and high hills; 6. Beech forests; 7. Montane grasslands; 8. Pseudo-alpine grasslands (beyond the tree line). Black symbols: natural vegetations; white symbols: secondary vegetations. Squares: open vegetations; circles: forest vegetations.
Diversity 16 00110 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fattorini, S. The Role of Vegetation in Elevational Diversity Patterns of Tenebrionid Beetles in Central Italy. Diversity 2024, 16, 110. https://doi.org/10.3390/d16020110

AMA Style

Fattorini S. The Role of Vegetation in Elevational Diversity Patterns of Tenebrionid Beetles in Central Italy. Diversity. 2024; 16(2):110. https://doi.org/10.3390/d16020110

Chicago/Turabian Style

Fattorini, Simone. 2024. "The Role of Vegetation in Elevational Diversity Patterns of Tenebrionid Beetles in Central Italy" Diversity 16, no. 2: 110. https://doi.org/10.3390/d16020110

APA Style

Fattorini, S. (2024). The Role of Vegetation in Elevational Diversity Patterns of Tenebrionid Beetles in Central Italy. Diversity, 16(2), 110. https://doi.org/10.3390/d16020110

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