Lichen Distribution Patterns in the Ecoregions of Italy

: An outline of the main distribution patterns of lichens in the ecoregions of Italy, accounting for their climatic, geographic, and environmental features, is still missing. On the basis of a GIS-based analysis, we summarized: (1) the main features (e.g., surface, climate, landscape, topographic heterogeneity, bedrock, eutrophication) of the 9 ecoregions adopted in ITALIC, the information system on Italian lichens, and (2) the patterns of richness, functional traits, and ecological requirements of lichens in the ecoregions. Our GIS-based analysis describes for the ﬁrst time the main features of the 9 ecoregions adopted in ITALIC, highlighting di ﬀ erences which could explain the main lichen patterns. Overall, the exploration of the Italian lichen biota is still a work in progress, some regions being still underexplored, especially in the South, with new taxa being reported every year. Our research could provide a baseline for further advancements in the understanding of species richness and community composition of Italian lichens, at a regional scale.


Introduction
Climatic, geographic, and environmental conditions are quite heterogeneous across Italy, resulting in substantial ecoregional diversity. This is reflected in distinctive lichen distribution patterns that likely result in strong differences of local community compositions [1]. A well-known example is described by Nimis and Tretiach [2], concerning the relevant difference in lichen community composition between the eastern and western parts of the Italian Peninsula.
In his first checklist of Italian lichens, Nimis [3] provided a first outline of the potential distribution of lichens within 9 ecoregions (Alpine, Oromediteranean, Subalpine, Montane, Dry and Humid Submediterranean, Padanian, Dry and Humid Mediterranean), which were originally delimited on the basis of altitudinal and bio-climatic constraints. For each ecoregion, a commonness-rarity score was provided for each taxon (for details see [4]), and recently patterns of commonness and rarity of Italian lichens were analyzed, indicating that rarity could depend on bioclimate and reduced habitat availability, while commonness may mainly reflect anthropogenic disturbance [5].
The distribution data included in the first checklist [3] were converted into the first release of ITALIC (http://italic.units.it) [4], the information system on Italian lichens. This online database, which is now available in its version 6.0, provides access to a wealth of taxonomic, ecologic, and distributional data on the lichens known to occur in the country. The system is freely accessible online, without registration. The data are reusable under a Creative Commons CC BY license. Currently, ITALIC [6] is based upon the latest national checklist [7], and it is being constantly updated by a flow of new data deriving from

Ecoregions and OGUs
The distribution of lichen taxa is reported as presence/absence in the administrative regions [5]. A commonness-rarity score [4], ranging from "extremely rare" to "extremely common", has been calculated for each infrageneric taxon, for each of the 9 ecoregions of the country (Figure 1) as originally delimited by Nimis [3]. Each administrative region includes several OGUs, deriving from the overlap of administrative and ecoregion borders. The delimitation of the 9 ecoregions mainly follows the altitudinal limits of the potential vegetation types, which can be very different from northern to southern, and from eastern to western Italy, so that altitude delimiting each ecoregion can differ among administrative regions (Table 1).
A GIS map was generated from the digital elevation model (DEM) of the country, obtained by the national geoportal (http://www.pcn.minambiente.it), with a resolution of 20 m. The original map was split into administrative regions, and inside each region the ecoregions were delimited by means of the r.reclass, r.to.vect, and v.patch of the software GRASS GIS 7.4 (https://grass.osgeo.org/), thus producing a GIS map of a total of 97 OGUs ( Figure 2). Subsequently, the ecoregions present in each region were extracted from each of the 20 raster maps using the assigned ranges of elevation (Table 1) using the function GRASS GIS r.reclass. The ecoregions occurring each administrative region were transformed into vectors (r.to.vect of GRASS GIS) and joined into a single map of all Italy (v.patch of GRASS GIS). The subdivision between "humid" and "dry" ecoregions is based on Nimis and Tretiach [2], in which the concept of "Thyrrenian Italy" and the influence of humid western winds to lichen communities were discussed. The "Padanian" ecoregion is a sub-Mediterranean portion of the country, including especially the whole Po plain, where high anthropization strongly affects lichen diversity [4]. Each of the 97 OGUs is made of a number of polygons ranging from 1 to n. Since computation complexity increases with the number of polygons, it was decided to take into account only the 20 larger polygons for each OGU. This allowed to drastically reduce computation complexity while retaining almost the whole surface of each OGU. In fact, discarding the smaller polygons led to the loss of 0.6% only of the total surface of the country. The number of ecoregions for each administrative region, together with their surface in percentage, are reported in Table 2.   Figure 1, and the text for details), and the elevation range are provided.

Functional Traits and Ecological Requirements
Information on species functional traits and ecological requirements were retrieved from ITALIC (see [4] for further details).
Reproductive strategy (mostly sexually; mostly asexual, by soredia or soredia-like structures; mostly asexual by isidia; mostly asexual, by thallus fragmentation) (B) Ecological requirements Substrata (epiphytic on bark or wood; foliicolous on evergreen leaves; saxicolous; terricolous) pH of the substratum (indicator values ranging from 1, very acid, to 5, strongly basic) Eutrophication (indicator values ranging from 1, no eutrophication, to 5, high eutrophication) Light (indicator values ranging from 1, very low solar radiation, to 5, very high solar radiation)

Environmental Data
All environmental data were retrieved from freely accessible online sources (see details below). Italy has a surface of about 300,000 square km, of which 35% is made of mountains, 42% of hills, and 23% of plains. Given the wide altitudinal and longitudinal gradients, the country hosts a wide array of climatic conditions. Landscape analyses are based upon the vectorial maps of Corine Land Cover, IV level, 2012 (https: //land.copernicus.eu/pan-european/corine-land-cover).
Landscape data were extracted for each OGU by using software v.overlay of GRASS GIS, within the environment Rstudio (4.0.2), and are: Forest coverage: percentage of woodland cover (code 3.1) Natural and semi-natural landscapes: percentage of surface occupied by patches of woodlands and semi-natural environments (code 3) Urbanized areas: percentage of surface occupied by patches of artificial environments (code 1) Diversity of the landscape: number of different Corinne IV level environments. Geological data were extracted from the Carta Geologica d'Italia, obtained by the National Geoportal (http://www.pcn.minambiente.it). Since this map has a low resolution at a 1:500,000 scale (http://www.isprambiente.gov.it/it/cartografia) we have produced a simplified map, relevant to lichen ecology, by partitioning areas with a prevalence of siliceous or carbonatic rocks. For each OGU, the percentages of both types of substrata were extracted.
Topographic data were extracted from the national DEM (http://www.pcn.minambiente.it) by means of the Zonal statistics function of QGIS. The data are: Minimum and maximum elevation. Elevation diversity, expressed as the standard deviation of the elevation in an OGU.
Latitude of the centroid (WGS84), taken from the wider polygon of each OGU. As far as pollutants are concerned, two deposition maps were used, both obtained from EMEP (European Monitoring and Evaluation Programme, http://www.emep.int/mscw/), which were used as proxies of eutrophication. Their resolution is 10 square kilometers.

Data Analysis
The ecological and lichenological features of the 9 ecoregions were descriptively analyzed. For each variable, box plot charts were produced, in order to describe trends in environmental, climatic, and landscape variables. Once it was verified that the values of variables were assuming a non-normal distribution (Shapiro-Wilk test), the Kruskal-Wallis test [11] was applied to test differences. When a p-value < 0.05 was obtained, the Wilcoxon test was applied to each of the possible pairs [12]) to check for differences between ecoregions.
Species-area relationships across the 9 ecoregions were explored by means of Arrhenhius power function [13]: S = c·Az Where S is the number of taxa, A the area, and c and z two fitted parameters. The former represents the number of taxa expected per area unit, the latter is the angular coefficient of the straight line in a log-log space. The model was produced using the function "SSArrhenius" from the vegan package [14], The graph was produced by means of the ggplot2 package [15]. All analyses were conducted in R version 4.0.2 [16].

Environmental Features of the Ecoregions
The Dry Submediterranean (D) is the widest ecoregion (   Landscape diversity (Figure 4c) is higher in lowland ecoregions and decreases with altitude, while the contrary applies to the coverage of natural and semi-natural habitats (Figure 4d), which has a minimum in the Padanian ecoregion. At high altitudes, the surface occupied by natural and semi-natural habitats is always higher than 85%, while in the lowlands it is always lower than 20%.
Woodland coverage is higher in the Montane (60%), and Subalpine (40%) ecoregions (Figure 4e), while the higher percentages of disturbed areas (Figure 4f), are in the Padanian and the Humid and Dry Mediterranean ecoregions. The Oromediterranean ecoregion has the highest geographic variability in the percent of woodland cover, being scattered in several administrative regions, and on different mountain peaks.   Average temperature ( • C), and precipitation (mm), as expected, have an opposite trend from high to low altitudes OGUs (Figure 4g,h). Low-altitude ecoregions encompass a broader latitudinal span than high-altitude ecoregions (Figure 5a), the Alpine and Subalpine ecoregions being limited to the Alps and Central Apennines, while the Oromediterranean ecoregion is limited to the southernmost high mountains. The Padanian ecoregion is limited to the Po Plain, where temperature seasonality (Figure 5b) is higher, indicating a relevant climatic continentality.

Lichen Richness in the Ecoregions
Species richness is positively correlated to total surface, topographic heterogeneity, coverage of natural and semi-natural habitats, and of woodlands, and is negatively correlated to the amount of urbanized, heavily disturbed areas ( Figure 6). Thus, the Montane ecoregion has the highest species richness, the Padanian ecoregion the lowest (Figure 7).

Lichen Richness in the Ecoregions
Species richness is positively correlated to total surface, topographic heterogeneity, coverage of natural and semi-natural habitats, and of woodlands, and is negatively correlated to the amount of urbanized, heavily disturbed areas ( Figure 6). Thus, the Montane ecoregion has the highest species richness, the Padanian ecoregion the lowest (Figure 7).
The ecoregions with the highest woodland coverage (Montane and Subalpine, Figure 4e), are also those with the highest species richness, hosting more than 50% of all taxa known to occur in the country. These two ecoregions, however, strongly differ in surface (Figure 3), as well as in latitudinal range (Figure 5a). The Subalpine ecoregion is restricted to the northern part of the country and has a very small surface (less than 2%), while the Montane ecoregion occurs in all administrative regions, and has a wider surface (ca. 13%). The Dry and Humid Submediterranean, Alpine, and Humid Mediterranean ecoregions host between 35% and 50% of the taxa known to occur in Italy. On the contrary, heavily disturbed areas have the lowest lichen diversity (Padanian and Dry Mediterranean ecoregions, Figure 4f), with less than 500 taxa each. The Oromediterranean ecoregion as well has a low number of taxa, but it occupies a very small surface, less than 0.5% of the country ( Table 2). The OGUs in the Humid Mediterranean ecoregion, while heavily anthropized, host a relatively high lichen diversity, comparable to that of the two Submediterranean ecoregions. Average latitude (Figure 5a) of the OGUs of the 9 ecoregions does not reflect species richness (Figure 7). The ecoregions with the highest woodland coverage (Montane and Subalpine, Figure 4e), are also those with the highest species richness, hosting more than 50% of all taxa known to occur in the country. These two ecoregions, however, strongly differ in surface (Figure 3), as well as in latitudinal range (Figure 5a). The Subalpine ecoregion is restricted to the northern part of the country and has a very small surface (less than 2%), while the Montane ecoregion occurs in all administrative regions, and has a wider surface (ca. 13%). The Dry and Humid Submediterranean, Alpine, and Humid Mediterranean ecoregions host between 35% and 50% of the taxa known to occur in Italy. On the contrary, heavily disturbed areas have the lowest lichen diversity (Padanian and Dry Mediterranean ecoregions, Figure 4f), with less than 500 taxa each. The Oromediterranean ecoregion as well has a low number of taxa, but it occupies a very small surface, less than 0.5% of the country ( Table 2). The OGUs in the Humid Mediterranean ecoregion, while heavily anthropized, host a relatively high lichen diversity, comparable to that of the two Submediterranean ecoregions. Average latitude (Figure 5a) of the OGUs of the 9 ecoregions does not reflect species richness (Figure 7).
The relationship between number of taxa and surface of ecoregions is reported in Figure 8. Species richness in 5 out of 9 ecoregions (Oromediterranean, Alpine, Humid and Dry submediterranean, and Humid Mediterranean) is well described by the model (significance level of the model: alpha = 0.05). Two ecoregions (Subalpine, Montane) have a larger number of taxa than expected, while two (Padanian and Dry Mediterranean) host far fewer taxa than expected.  Table 3 Reports the incidence of taxa with oceanic, suboceanic, and subcontinental distribution patterns in the ecoregions. At lower elevations (Mediterranean and Submediterranean ecoregions), the basic difference between the western and eastern sides of the Italian Peninsula is confirmed, with a higher incidence of oceanic-suboceanic species with subtropical affinities along the western side, which is subject to the influence of humid, maritime air masses. On the contrary, very low values are reached in the Alpine Subalpine, Oromediterranean, and Dry Mediterranean ecoregions. Subcontinental species have the highest incidence in the Dry Mediterranean ecoregion, mainly due to the occurrence of steppic or desert species.

Lichen Traits and Ecological Requirements in the Ecoregions
As far as substrates are concerned (Figure 9a), epiphytic taxa are obviously scarce in the highaltitude northern ecoregions (Alpine and Subalpine), saxicolous and terricolous species being the  Table 3 Reports the incidence of taxa with oceanic, suboceanic, and subcontinental distribution patterns in the ecoregions. At lower elevations (Mediterranean and Submediterranean ecoregions), the basic difference between the western and eastern sides of the Italian Peninsula is confirmed, with a higher incidence of oceanic-suboceanic species with subtropical affinities along the western side, which is subject to the influence of humid, maritime air masses. On the contrary, very low values are reached in the Alpine Subalpine, Oromediterranean, and Dry Mediterranean ecoregions. Subcontinental species have the highest incidence in the Dry Mediterranean ecoregion, mainly due to the occurrence of steppic or desert species. Table 3. Incidence (absolute numbers and percentages) of taxa with oceanic, suboceanic, and subcontinental phytoclimatic range in the lichen flora of the 9 ecoregions.

Lichen Traits and Ecological Requirements in the Ecoregions
As far as substrates are concerned (Figure 9a), epiphytic taxa are obviously scarce in the high-altitude northern ecoregions (Alpine and Subalpine), saxicolous and terricolous species being the most abundant, the latter reaching their maximum, up to 20% of the total, due to the presence of Alpine tundras. Saxicolous taxa have the maximum in the Dry Mediterranean ecoregion, which has a scanty woodland cover. Epiphytes have the maximum (almost 50%) in the Montane, the Submediterranean, and the Humid Mediterranean ecoregions. Epiphyllous taxa are a very minor component of the Italian lichen flora (30 taxa), which occur almost exclusively in the Humid Mediterranean ecoregion, which is well in accordance with their mainly tropical affinities. Except for the type of photobiont, species traits do not show distinctive patterns across ecoregions. Crustose taxa are 65-70% of the whole lichen biota (Figure 9b). Endolithic taxa, however, occur mostly in ecoregions with a small coverage of woodlands

Discussion
Our GIS-based analysis allowed us to describe, for the first time, the main features of the 9 ecoregions adopted in ITALIC, highlighting the main differences that could explain the main patterns of lichen richness, functional traits, and ecological requirements. While the 9 ecoregions are bioclimatically well-delimited, they are also part of a continuum, and often share similar features, which gradually change with altitude and latitude. The ecological factors used for delimiting the ecoregions are those which are known to considerably influence lichen species richness and diversity, shaping the lichen biota of the country.
The areas of the OGUs belonging to each ecoregion are strongly affected by altitude, with generally smaller OGUs at higher elevations. Theoretical models of species-area relationship (SAR) evidence that the increase in area is often correlated with an increase in habitat heterogeneity, thus leading to higher richness [17][18][19][20]. However, in Italy, the area occupied by natural and semi-natural habitats strongly

Discussion
Our GIS-based analysis allowed us to describe, for the first time, the main features of the 9 ecoregions adopted in ITALIC, highlighting the main differences that could explain the main patterns of lichen richness, functional traits, and ecological requirements. While the 9 ecoregions are bioclimatically well-delimited, they are also part of a continuum, and often share similar features, which gradually change with altitude and latitude. The ecological factors used for delimiting the ecoregions are those which are known to considerably influence lichen species richness and diversity, shaping the lichen biota of the country.
The areas of the OGUs belonging to each ecoregion are strongly affected by altitude, with generally smaller OGUs at higher elevations. Theoretical models of species-area relationship (SAR) evidence that the increase in area is often correlated with an increase in habitat heterogeneity, thus leading to higher richness [17][18][19][20]. However, in Italy, the area occupied by natural and semi-natural habitats strongly increases with altitude, which positively affects lichen richness. Thus, the combined effect of reduced area and increase of "naturality" weakens the positive correlation between area and species richness. In particular, the Subalpine ecoregion, restricted to the northern regions, and with a small surface, has a high species richness, comparable to that of the much broader Montane ecoregion. In particular, larch and stone-pine forests of the Subalpine ecoregion were found to be a hotspot of lichen diversity of particular conservation concern [21][22][23]. As far as the Montane ecoregion is concerned, the high species richness is probably related to its wide latitudinal extent, since it ranges throughout the country, spanning a relevant biogeographic gradient [24,25]. This should be reflected in heterogeneous community composition across sites [1], a pattern which would, however, require further investigation.
In the Padanian and Dry Mediterranean ecoregions, species richness is negatively impacted by high anthropization [26], even where the landscape is heterogeneous. The Padanian ecoregion is the most disturbed, virtually lacking natural or semi-natural habitats that were replaced by intensive agriculture, industrial areas, and urban areas. Under these circumstances, depositions of both nitrogen oxides and ammonia, as well as habitat destruction and fragmentation directly impact on species richness [27,28]. The interplay between these drivers caused the historical depauperation of a previously rich and specialized lichen biota that was replaced by a few generalist species which are homogenizing the lichen communities [29,30].
Species richness is also associated with topographic heterogeneity, likely due to its effect on microand meso-climatic conditions. This is particularly evident in the high species richness of the Alpine, Subalpine, and Oromediterranean ecoregions, even if they account for a very small surface (ca 6.5%). This pattern is consistent with the outstanding richness of the lichen biota of the Alps [5]. However, besides environmental factors, this situation may also reflect an exploration bias, since high-elevation areas have a long-lasting tradition of lichenological surveys.
Our analysis also revealed some robust patterns of species functional traits and ecological features across the 9 ecoregions. For example, lichens with Trentepohlia photobionts have their maximum in Submediterranean and Humid Mediterranean ecoregions, which is in accordance with their affinity for warm-humid conditions [31][32][33]. At low temperatures, these taxa are less competitive than other lichens [34], since trentepohlioid algae are more prone to freezing [35]. While a sexual reproductive strategy prevails in all ecoregions, thallus fragmentation, a particular type of asexual reproduction, seems to be clearly linked to high altitudes, as already observed by [36]. This may reflect an adaptation to cold conditions [37], and to the presence of lichen-grazing mammals in arctic-subarctic ecosystems dominated by tundras with many species of Cladonia, in which this strategy is often adopted.
Substrate availability is a major driver of lichen distribution in the different ecoregions. Epiphytic lichens are scarce or absent in the Alpine ad Oromediterranean ecoregions, which by definition are devoid of trees. Terricolous species are most abundant in sites where agriculture and urbanization are limited, or absent. Saxicolous species are most abundant in ecoregions with low coverage of woodlands, as the Alpine, Oromediterranean, and Dry Mediterraneans. This does not apply to the Padanian ecoregion, where the prevalence of sedimentary substrates, heavy anthropization, and intensive agricultural practices reduce the availability of saxicolous habitats.
Substrate pH is another relevant factor which contributes to lichen species richness [38]. Epiphytic species growing on acidic bark prevail in the communities of the Subalpine and Montane ecoregions, where the incidence of coniferous forests (Corine Land Cover code 3.2.2) is higher. The Subalpine ecoregion is characterized by woods dominated by Larix decidua and Picea abies, while Abies alba and Fagus sylvatica, whose bark is acidic as well, are characteristic of the Montane ecoregion [39]. Acidophilic saxicolous and terricolous lichens reach a maximum at high altitudes. Nimis and Schiavon [40] in a study of epiphytic lichen vegetation along the western coast of peninsular Italy, showed that this is basically different from that of the eastern coast. The western side of the Italian Peninsula has a mild-humid climate under the influence of Tyrrhenian maritime winds, while the eastern, Adriatic side, located on the lee-side of the Apennines, is subject to cold-dry air masses from the Balkans during winter. This basic bioclimatic difference results in important compositional differences between the humid and dry Mediterranean/Submediterranean ecoregions, the former hosting a much higher number of species with oceanic/suboceanic distribution patterns and tropical/subtropical affinities.
Overall, the exploration of the Italian lichen biota is still a work in progress [3]. Some regions are still underexplored, especially in the South, new taxa being reported every year (e.g., [41]). Our research could provide a baseline for further advancements in the understanding of species richness and community composition of Italian lichens at a regional scale.