Next Article in Journal
Mapping QTL for Phenological and Grain-Related Traits in a Mapping Population Derived from High-Zinc-Biofortified Wheat
Next Article in Special Issue
Population and Conservation Status of Buxbaumia viridis (DC.) Moug. & Nestl. in Romania
Previous Article in Journal
Predictive Modeling of Kudzu (Pueraria montana) Habitat in the Great Lakes Basin of the United States
Previous Article in Special Issue
Quo Vadis, Orthotrichum pulchellum? A Journey of Epiphytic Moss across the European Continent
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temperature and Precipitation More Than Tree Cover Affect the Distribution Patterns of Epiphytic Mosses within the Orthotrichaceae Family in China and Adjacent Areas

1
Department of Botany, University of Ostrava, Chittussiho 10, 710 00 Ostrava, Czech Republic
2
Institute of Biology, University of Opole, Oleska 48, 45-052 Opole, Poland
3
Institute of Botany, Faculty of Biology, Jagiellonian University, Gronstajowa 3, 30-387 Kraków, Poland
4
College of Life and Environmental Sciences, Shanghai Normal University, 100 Guilin Road, Shanghai 200234, China
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(1), 222; https://doi.org/10.3390/plants12010222
Submission received: 17 November 2022 / Revised: 17 December 2022 / Accepted: 30 December 2022 / Published: 3 January 2023
(This article belongs to the Special Issue New Knowledge in Bryology 2.0)

Abstract

:
Epiphytes, including vascular and non-vascular, constitute a large part of global plant biodiversity. Distribution of obligatory epiphytic bryophytes results from climate and local habitat conditions. The most important epiphytic bryophytes and at the same time poorly investigated and taxonomically problematic ones belong to the family Orthotrichaceae. Epiphytic mosses are also ideal organisms for species modelling, because of having no roots, they are highly dependent on external environmental conditions. For this purpose, we used the ecological niche modelling approach to define their potential distribution in China and adjacent areas and explore factors that shape this distribution. We used 617 occurrence records of 23 species from six genera within the Orthotrichaceae family. Our results suggest that the distribution of members of the Orthotrichaceae family is predominantly affected by bioclimatic variables, especially bio10 (mean temperature of the warmest quarter), bio15 (precipitation seasonality), bio18 (precipitation of the warmest quarter), bio19 (precipitation of the coldest quarter), bio9 (mean temperature of the driest quarter), and bio2 (mean diurnal range). However, the distribution of particular genera is ruled by a different set of those variables. The distribution of two genera (Leratia and Ulota) is also highly influenced by land cover (especially mixed/other trees), whereas human footprint shows a moderate contribution to models of three genera (Lewinskya, Orthotrichum, Nyholmiella). Based on the occupied climatic niche and distribution patterns, representatives of the studied family are divided into two groups. The ‘western-montane group‘ is characterised by lower temperatures and lower precipitation whereas the ‘eastern-lowland’ group‘ by more humid and warmer conditions.

1. Introduction

The distribution of natural vegetation at large spatial scales is directly impacted by climate [1,2,3]. However, potential drivers of species-richness patterns that have also been considered are environmental factors such as water availability and habitat heterogeneity [4,5,6,7,8]. In the case of obligatory epiphytic bryophytes, patterns of distribution are the results of the relationship between climate and specific environments [9,10,11,12]. Epiphytic bryophytes are generally willing to occupy the bark of trees both in temperate forests and in more humid tropical areas [13]. With approximately 850 widely distributed subcosmopolitan species, the family Orthotrichaceae [14,15,16,17,18,19,20,21,22] is one of the main components of the cryptogam epiphytic communities and occurs on tree trunks and branches, and occasionally they are found also in saxicolous habitats. In Central and East Asia, the family is represented by eight genera, namely Leratia, Lewinskya, Macrocoma, Macromitrium, Nyholmiella, Orthotrichum, Ulota, and Zygodon, which can be found in most biomes, from tundra to wet tropical forests. Representatives of Macromitrium and partly also Ulota are associated with more humid tropical areas and are characterised by resistance to strong competition [23,24]. Whereas the other epiphytic taxa of the family, which are less steady competitors in comparison to other species of bryophytes and lichens, prefer rather sparse forests, bushes, and also alleys or solitairy trees that are exposed to harsh, less humid, or even semi-dry climatic conditions [11,12,25,26,27,28]. Therefore, species within such genera as Macromitrium and Leratia are very often observed in the south-eastern part of China [23,24,29,30], where the climate is warmer and more humid, and the temperature difference between day and night smaller [31], whereas representatives of Lewinskya, Nyholmiella, and Orthotrichum are much more frequently observed in the mountain areas in the central and western part of the country [32,33,34]. This observation is related to the distribution of climatic zones (wet and warm vs. dry and cold) in China and generally in eastern Asia [31,35,36].
Ecological niche modelling (ENM) is a widespread method to determine species distribution patterns. In simple terms, ENM uses computer algorithms to predict the distributions of species across geographic space and time by using occurrence data and environmental data. ENM estimates the relationship between species records at sites and the environmental and/or spatial characteristics of those sites [37]. Despite the frequent use of this method, and despite the fact that mosses are an ideal model organism for species modelling because, having no roots, they are highly dependent on external environmental conditions, the general patterns of distribution of epiphytic mosses from the Orthotrichaceae family are still underexplored.
In this paper, we aim to: (i) investigate macroscale factors affecting the distribution of epiphytic mosses from the family Orthotrichaceae in Central and East Asia; (ii) find the differences in ecological niches occupied by representatives of the family Orthotrichaceae; and (iii) assess the impact of anthropopressure on species distribution patterns. The results can help to define factors influencing the geographical distribution patterns of the Orthotrichaceae family and point out hotspots of the highest epiphyte diversity and at the same time the most important areas for conservation purposes.

2. Materials and Methods

2.1. Characteristics of the Studied Genera

In this study, distribution data of species from six genera within the Orthotrichaceae family (Leratia, Lewinskya, Macromitrium, Nyholmiella, Orthotrichum, and Ulota) were evaluated. Brief descriptions of the genera and their taxonomic distinctness are explained below. The genera Macrocoma and Zygodon, which also occur in the study area, were not included in the modelling, as their taxonomic revision is still ongoing and detailed data on the distribution of particular species was not available.

2.1.1. Leratia Broth. & Paris

The taxonomy of this genus was not completely clear in the past. It includes three very geographically isolated species. Goffinet et al. [15], using molecular data, clearly proved that the Asian and Australasian Leratia obtusifolia, North American L. exigua and New Caledonian L. neocaledonica belong to an isolated and strongly supported clade, to which they applied the generic name Leratia. Only Leratia exigua was used for modelling, as it is the only representative of this genus occurring in China and other parts of the Asian continent [34].
Leratia exigua can be distinguished from similar orthotrichaceous mosses by a combination of features: decurrent leaves with plane margins; superficial stomata and hyaline; and keeled endostome segments [38].

2.1.2. Lewinskya F.Lara, Garilleti & Goffinet

The genus Lewinskya was described as a result of taxonomic and molecular studies that confirmed polyphyly of the broadly conceived traditional genus Orthotrichum s.l. Three segregates have been split from it, including Lewinskya, Nyholmiella, and Pulvigera [19,39,40,41]. Species of the genus Lewinskya differ from the representatives of the genus Orthotrichum mainly in the presence of superficial stomata [19]. Currently, the genus Lewinskya comprises 70 taxa worldwide [41], of which 14 Lewinskya species and two varieties have been recorded in China [34]. Here for the purpose of distribution modelling, we used data for the six most common species, namely L. affinis, L. exigua, L. hookeri, L. speciosa, L. striata, and L. vladikavkana.

2.1.3. Macromitrium Brid

Among the epiphytic bryophytes, Macromitrium is considered one of the largest moss genera in the world with the majority of its species distributed in tropical and subtropical regions [42,43]. According to Vitt and Ramsay [44], there are 350 species worldwide but only 128 can be morphologically differentiated [42]. In China, there are 17 species and one variety for which occurrence was recently confirmed [38]; however, for the purpose of distribution modelling we used data of the two most common species, M. cavaleriei and M. gymnostomum.

2.1.4. Nyholmiella Holmen & Warncke

The taxonomic position of the genus Nyholmiella has been uncertain in the past and it has been alternately classified as either part of the widely conceived genus Orthotrichum or as a separate genus [45,46]. Sawicki et al. [17] proposed the first systematic consequences since 2004 to reflect phylogenetic reconstruction, and reinstated the name Nyholmiella. Currently, the genus Nyholmiella comprises 3 taxa worldwide [17,47]. Only N. obtusifolia was used for distribution modelling, as it is the only representative of this genus occurring in China and adjacent areas [34]. Historically, the species N. gymnostoma was also reported from the territory of China [48], but after a taxonomic revision it was proved to be N. obtusifolia, and consequently, N. gymnostoma has been excluded from the list of Chinese bryophytes [34].

2.1.5. Orthotrichum Hedw

The genus Orthotrichum is the only orthotrichaceaous moss—of those which occurred in the study area—with immersed stomata [39]. According to Lara et al. [41], the genus comprises 103 taxa worldwide, of which the Northern Hemisphere hosts 83%. In the course of a new taxonomical revision, 29 species and two varieties have been reported from China [34]. Considering that the genus is the most numerous group of all orthotrichaceous mosses in China and the surrounding Asian countries, we used 14 representatives for distribution modelling, namely Orthotrichum alpestre, O. anomalum, O. callistomum, O. consobrinum, O. crenulatum, O. crispifolium, O. dasymitrium, O. griffithii, O. pallens, O. pamiricum, O. pumilum, O. scanicum, O. sordidum, and O. urnigerum.

2.1.6. Ulota D.Mohr

Recent taxonomic studies have shown that the genus Ulota is polyphyletic [19,21,22]. First, based on molecular studies, Ulota phyllantha was displaced into the new genus Plenogemma [19] and latterly U. bellii into the genus Rehubryum [22]. According to Wang and Jia [49], the genus Ulota s. str. comprises about 50 taxa worldwide.
Species of Ulota are distinguished from other orthotrichaceous mosses based on a combination of characters—crisped to slightly twisted leaves (when dry) and leaf base with a well-differentiated border formed by one or a few rows of quadrate to rectangular hyaline cells [39]. The occurrence of eight species is recently known from China [29,38]. Distribution data of six of them, namely U. crispa, U. gigantospora, U. gymnostoma, U. perbreviseta, U. robusta, and U. yunnanensis were used here for modelling.

2.2. Occurrence Data

The data on the distribution of species that were used for modelling comes from our own collections, a query of the literature data, and also from the revision of the moss specimens housed in major European and Chinese herbaria. Field research in twelve provinces of China, including Anhui, Guizhou, Henan, Hubei, Jiangsu, Shaanxi, Shanxi, Shanghai, Sichuan, Xinjiang, Yunnan, and Zhejiang were carried out in 2015–2019. In addition, the authors carried out several bryological expeditions to Central Asia (Tajikistan, Kirgizstan, Kazakhstan, and Afghanistan), where the Orthotrichaceae family was studied in detail in the years 2011–2015. A complete revision of herbarium items of epiphytic bryophytes was also carried out in the herbarium of the Botanical Institute of the Tajikistan Academy of Sciences in Dushanbe.
During the study, the herbarium collections from the following Chinese and European herbaria have been revised (acronyms of the herbaria according to Index herbariorum but for some Chinese herbaria they are unauthorised and introduced only for the purposes of the present article; they are distinguished by not using bold letters):
C—University of Copenhagen, Denmark;
E—Royal Botanic Garden Edinburgh, Scotland, U.K;
GuZU—Guizhou University, China;
H—University of Helsinky, Finland;
HSNU—Herbaria of East China Normal University, China;
HuBU—Hubei University, China;
IFP—Institute of Applied Ecology, Academia Sinica, Shenyang, China;
InMU—Inner Mongolia University, China;
KRA—Herbarium of the Inst. of Botany, Jagiellonian University, Kraków, Poland;
KRAM—Polish Academy of Sciences, Krakow, Poland;
KUN—Herbarium of Kunming Inst. of Botany, the Chinese Academy of Sci., China;
OSTR—University of Ostrava, Czech Republic;
PC—Muséum National d’Historie Naturelle, Paris, France;
PE—Institute of Botany, Chinese Academy of Science, Beijing, China;
ShM—Shanghai Museum, China;
ShU—Shanghai Normal University, China;
TAD—Botanical Inst. of the Tajikistan Academy of Sciences in Dushanbe, Tajikistan;
W—Naturhistorisches Museum, Wien, Austria;
XJU—Xinjiang University, China.
Because the Orthotrichaceae family is considered to be a taxonomically problematic group, whose determination is not elementary, we decided to use mainly distribution data from our own field research (33.8%), as well as revised herbarium specimens (64.4%), supplemented only in a small percentage by records from online databases (1.62%) and distribution data from published papers (0.18%). To prevent possible misidentification of taxa, all available individuals were morphologically determined or re-determined in the BRYOLAB laboratory (part of the Department of Biology and Ecology, at the University of Ostrava).
In summary, 2717 geographic occurrence records were collected from the area of Central Asia (e.g., Afghanistan, Kazakhstan, Kyrgyzstan, and Tajikistan), East and Southeast Asia (e.g., China, Japan, Mongolia, Taiwan, and Vietnam). To avoid problems resulting from inaccurate descriptions of localities we excluded from analysis data with missing geographic coordinates as well as data collected before 1960 [50]. To remove sampling bias, we first deleted duplicate records and then filtered the data through a grid with a WGS84 format and resolution of cells 0.1 × 0.1 -degree, as described in Hijmans et al. [51]. For filtering data, we used a sample function from dismo package in R software [52]. Preparation of occurrence data was performed using the open source software QGIS [53]. As a result, we obtained 617 occurrence records for the whole family Orthotrichaceae, including Leratia (22), Lewinskya (312), Macromitrium (110), Nyholmiella (53), Orthotrichum (354), and Ulota (19).
In addition to data analysis for individual genera, we also evaluated two ecological groups of mosses connected with different distribution patterns: and western-montane epiphytic group (479 records) and eastern-lowland epiphytic group (191 records). Based on distribution records we included Lewinskya affinis, L. hookeri, L. speciosa, L. vladikavkana, Nyholmiella obtusifolia, Orthotrichum alpestre, O. anomalum, O. callistomum, O. crenulatum, O. dasymitrium, O. pallens, O. pamiricum, O. scanicum, O. sordidum, O. urnigerum, and Ulota robusta in the western-montane epiphytic group, whereas Leratia exigua, Macromitrium gymnostomum, M. cavaleriei, O. consobrinum, O. crispifolium, O. griffithii, and Ulota gymnostoma were included in the eastern-lowland epiphytic group. In the case of species with wide distributions, the species were assigned to a particular group depending on their frequency (the majority of localities).
Nomenclature follows Lara et al. [41] and Sawicki et al. [19]. Duplicates of all herbarium specimens together with our own materials collected during field research in Central Asia and China are deposited in OSTR and KRAM.

2.3. Ecological Niche Modelling

For modelling the potential distribution of epiphytic mosses from the family Ortotrichaceae we used MaxEnt software 3.3.3, which is a generative species distribution modelling approach recommended for applications involving presence-only datasets and environmental predictors [54]. For our purpose, 25 variables of continuous environmental predictors, including 19 bioclimatic variables from WorldClim v.2.1 [55], 5 global 1-km consensus land cover v.2.2 [56], and human footprint v.3 [57] were initially used (Table 1). All variables were evaluated with distribution data for the entire Orthotrichaceae family. All environmental data were derived with a spatial resolution of 30 arcsec (approx. 1 km). Because one of the aims of the research was to create a general model of the distribution of epiphytic mosses of the Orthotrichaceae family in the region of Asia, the background area was defined by the geographic boundaries of Afghanistan, Bhutan, Cambodia, China, India, Iran, Japan, Kazakhstan, Kyrgyzstan, Laos, Mongolia, Myanmar, Nepal, North Korea, Pakistan, Philippines, South Korea, Taiwan, Tajikistan, Thailand, Turkmenistan, Uzbekistan, and Vietnam [58]. The number of background points was set to 10,000 (default value) [59].
Environmental covariates were clipped according to the mask and resampled to the same resolution and projection as the bioclimatic variables using the extract by mask function in ArcMap 10.5 software [60]. For each layer, the WGS-84 geographical coordinate system was chosen.
To avoid overfitting of the model we built a correlation matrix (Pearson’s correlation coefficient) and removed highly correlated variables (r > 0.7) [61], by generating a Pearson correlation matrix in Statistica software [62]. To choose which of the strongly correlated variables to remove, we performed a jackknife test of variable importance [63] and compared the percent contribution of the variables. We chose variables that showed higher positive gain values and higher percent contribution or those whose impact was easier to explain from a species biology point of view. For the final models, we used 15 variables, including 9 bioclimatic variables, human footprint, and 5 landcover classes (Table 1).
We ran models with default values (a maximum of 500 iterations, convergence threshold 0.00001, and five auto feature classes). We decided to choose a logistic format, as it is currently considered easier and potentially more accurate for interpretation than the cumulative and raw approaches [59,63]. The model was calibrated using 75% of the occurrence records and tested on the remaining 25%. We performed 20 replicates using the subsample replicated run type and then averaged the results. To provide a different random test/train partition in each replicate, we used the ‘random seed’ option. Because the distribution data came from different sources with different levels of accuracy, we decided to use a 10 percentile training presence as a threshold rule, which eliminates the most outlying data. We evaluated the final model using the area under the curve (AUC) values, which is a threshold-independent evaluation analysis generated by Maxent [48]. We applied the AUC classification according to Swets 1988: <0.6 invalid; 0.6–0.7 weak; 0.7–0.8 fair; 0.8–0.9 good; and 0.9–1.0 excellent. The feature selection was based on a number of occurrence data [59]. The regularization multipliers parameter was set to a default value 1.0.
To derive areas that could serve as potentially suitable for analysed representatives of the family Orthotrichaceae we used maximum training sensitivity plus a specificity-threshold value that was generated for each genus and group independently. The potential probability of occurrence was visualised on the map with a 5-degree scale with breakpoints: 0.2 (very low); 0.4 (low); 0.6 (medium); 0.8 (high); and 1.0 (very high). Probability values were reclassified by using ArcMap 10.5 software [60]. To calculate the overlap between ENMs we first transformed them to binary format using ArcMap 10.5 software [60]. Then we used SDM Toolbox v2.0 [64] to calculate the area of overlap. We depicted predicted ‘western-montane’ (area predicted as potentially available), ‘eastern-lowland’ (area predicted as potentially available), and ‘both’ (area predicted as potentially available for both studied groups) areas in the species distribution.
To compare realised ecological climatic niches of two studied groups of epiphytic mosses (‘western-montane’ vs. ‘eastern-lowland’) we conducted principal component analysis (PCA), which is a method widely used for analyzing ecological niches [65,66]. We conducted PCA on the basis of the correlation matrix using Statistica software [62]. We analysed 7 bioclimatic variables. In PCA the localities were grouped with no a priori assumptions, then marked with the symbol on the scatter plot corresponding to a particular ecoregion.

3. Results

3.1. Distribution Patterns and Ecological Preferences of Particular Genera within the Orthotrichaceae Family

MaxEnt models show excellent performance in predicting current suitable conditions for all studied genera, with an average AUC exceeding 0.9 (Table 2). In general, the bioclimatic variables connected with temperature and precipitation have the highest contributions to the models, especially bio10 and bio15, but also bio18, bio19, bio9 and bio2. However, different bioclimatic variables show the highest contributions to the models of particular genera (Table 2). Only one of the land cover variables (mixed/other trees; cons4) plays a significant role for two genera (Leratia and Ulota), whereas human footprint shows a moderate contribution to models of three genera—Lewinskya, Orthotrichum, and Nyholmiella (Table 2).

3.1.1. Leratia

The predicted current distribution model (Figure 1) shows a high affinity of the species distribution with the southeastern and eastern parts of the studied area. These are primarily territories with a humid and warmer climate (e.g., SE China), or coastal regions with a mild oceanic climate (Korea, Japan).
The highest percent contributions are cons4 (48.6%) and bio19 (16.9%) (Table 2). Resulting plots of environmental conditions cons4 and bio 19 show that the probability of occurrence of genus Leratia increases with an increase in the percentage of land cover by mixed trees (up to 30% coverage) and when the precipitation in the coldest quarter does not exceed 500 mm.

3.1.2. Lewinskya

The predicted current distribution model (Figure 2) shows that the highest probability of distribution of the studied species lies in several, more or less isolated, regions of Asia. Within China, these areas are mainly found in the mountainous regions of Yunnan and Sichuan provinces. Similarly, suitable conditions are provided by the mountain massifs of Central Asia (e.g., Pamir, Tian-Shan, and Altai). Mountain valleys in the north of Iran (Elburz Mts.) or, on the contrary, in the east part of China (Changbai Mts.) also show a high affinity for species distribution.
The highest percent contributions are bio10 (40.8%) and bio15 (23.1%) (Table 2). Occurrence probability is the highest with ca. 10 up to 15 °C for mean temperature of the warmest quarter and ca. 40–80% of precipitation seasonality.

3.1.3. Macromitrium

This distribution model for Macromitrium species (Figure 3) shows that the distribution probability of the studied species is very similar to that of the genus Leratia (Figure 1). The highest affinity for these species is predicted for the southeastern provinces of China, which are characterised by a warmer and more humid climate, as well as to the coastal regions in the east of China and Japan, where an oceanic climate prevails.
The highest percent contributions are bio2 (31.2%) and bio18 (27.3%) (Table 2). The probability of occurrence decreases as the mean diurnal temperature range increases (the highest with values lower than 10). In the case of bio 18, the highest probability of occurrence is when the precipitation of the warmest quarter ranges from 300 up to 1200 mm.

3.1.4. Nyholmiella

The predicted current distribution model (Figure 4) shows that the highest probability of distribution of the studied species is mainly in the mountain ranges of northern Iran, central Asia (Tian-Shan and Altai Mts.) and especially in mountains within central Chinese provinces, e.g., Sichuan, Shaanxi, and Gansu.
The highest percent contributions are bio10 (39.7%) and bio15 (14.5%) (Table 2). Accordingly, the response curve with the highest probability of occurrence is ca. 11–16 °C of the mean temperature of the warmest quarter and ca. 40–60% of precipitation seasonality.

3.1.5. Orthotrichum

The predicted current distribution model (Figure 5) shows that the highest probability of distribution of the studied species is in mountain ranges, e.g., Elburz Mts. (Iran), Kungey Alatau Mts. (Kyrgyzstan), Tian Shan Mts. and Altay Mts. (Xinjiang, China), Yunnan–Guizhou Plateau (China), Daba Mts. (along the Gansu-Sichuan border, China), Central Mts. and East Coast Mts. (east part of Taiwan), and the easternmost portion of the Japanese island of Hokkaidō (area of volcanoes).
The highest percent contributions are bio15 (34.6%) and bio10 (33.9%) (Table 2). Occurrence probability is the highest around 40–70% of precipitation seasonality, and ca. 10–15 °C mean temperature of the warmest quarter.

3.1.6. Ulota

The predicted current distribution model (Figure 6) shows that the highest probability of distribution of the studied species is in the Himalayan range passing through India, Nepal, and Bhutan, the Yunnan–Guizhou Plateau (the north part of Yunnan province, China), south part of Sichuan province (China), south-west part of Hubei province (China), and the central mountain range (Taiwan).
The highest percent contributions are cons4 (25.5%) and bio10 (23.3%) (Table 2). Models show that the best conditions for the genus Ulota are when the area covered with mixed trees is higher than 30% and the mean temperature of the warmest quarter is ca. 7–16 °C.

3.2. Potential Distribution of Ecological Groups within Epiphytic Bryophytes

Based on climatic differences within the vast study area, we defined two ecologically distinct groups whose occurrence reflects these contrasting conditions, and prepared models of their potential distribution.
MaxEnt models produced very accurate results in predicting current suitable conditions for both studied ecological groups (western-montane and eastern-lowland epiphytes), with an average AUC exceeding 0.9 (Table 2). It was proved, that bioclimatic variables connected with temperature and precipitation have the highest contributions to the models, especially bio10, bio15, bio18, and bio9. However, the highest impact on the models of a particular group has a different set of bioclimatic variables (Table 2).

Potential Distribution and Niche Comparisons of ‘western-montane’ and ‘eastern-lowland’ Epiphytic Groups of Species

The model shows the distribution patterns roughly congruent with the observed distributions for both studied groups (Figure 7A and Figure 8A,B). However, potentially suitable areas are much wider than those actually occupied, especially in the case of the western-montane group of epiphytic species (Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7A and Figure 8A,B). The area of potentially suitable habitats is almost three times wider for the western-montane epiphytic group than for the eastern-lowland one (3,805,764 km2 vs. 1,617,601 km2, respectively). In general, the eastern-lowland group occupies more south-east regions; however, both studied groups may also occupy the same geographical space (ranges overlap in an area of 1,658,565 km2) (Figure 7A and Figure 8A,B). For both western-montane and eastern-lowland groups of epiphytic mosses, the highest contributions to the model are bioclimatic variables, bio10 and bio15 for the western-montane group as well as bio9 and bio18 for the eastern-lowland group (Table 2, Figure 7B). Based on the current niche model occurrence probability, the western-montane epiphytic group has an optimum from 9 to 16 °C for mean temperature of the warmest quarter (bio10) and 40–80% of precipitation seasonality (bio15). The response curves based on the current niche model for the eastern-lowland epiphytic group indicate that occurrence probability is highest from about 250 mm for precipitation of the warmest quarter (bio 18) and not more than 12 °C for mean temperature of the driest quarter (bio9).
The principal component analysis (PCA), performed on seven bioclimatic variables, shows two somewhat overlapping groups of points corresponding to western-montane and eastern-lowland epiphytic groups of species (Figure 7C). The first two principal components explain over 69% of the total variance: 47.87 and 21.33% for the first and second axis, respectively. In the climate space, both studied groups differed in the first component (Figure 7C), in which the three variables (bio2, bio9, bio19), and to a lesser extent also the next two (bio10, bio18) showed high correlation. In general, species classified to the eastern-lowland epiphytic group occupy warmer and wettest regions than the second examined group (Figure 7C).

4. Discussion

The results give a useful overview of suitable habitat conditions for epiphytic bryophytes in China and adjacent areas. The potential distribution of the studied mosses largely reflects the topography as well as the distribution of climatic and vegetation zones in the studied area [31,67]. Within the Orthotrichaceae family, there can be observed two species-specific groups characterised by a difference in ecological niches linked to the climatic zones observed in eastern Asia [31,68]. The first one named here eastern-lowland group is related to a more humid, temperate, and subtropical climate, while the second one named western-montane group is related to a dryer and colder climate. Occurrence of the epiphytic mosses from the Orthotrichaceae family is predominantly influenced by climate (temperature and precipitation), which is not surprising because due to their physiology, mosses are dependent on air humidity, water availability, and drought frequency [10,69]. Similar conclusions, but on a smaller spatial scale, have been presented by Ma et al. [23] and Lou et al. [30] in modelling the habitat preferences of selected Macromitrium species in China. The occurrence probability of these species increased especially with increasing annual precipitation.
The occurrence of epiphytic mosses is highly related to the tree species identity as well as forest continuity [10,27,28,70]. However, in our epiphytic mosses, the tree cover (especially consensus prevalence of mixed to other trees) affects especially the distribution of Leratia and Ulota, for which occurrence probability increases along with increases in the tree cover.
Anthropopression seems to have less importance on the distribution of the studied group of mosses than we initially expected. Human footprint has a relatively small or medium impact on the occurrence of the whole family Orthotrichaceae. However, the probability of occurrence of some genera (e.g., Lewinskya and Nyholmiella) increases to a certain degree with the severity of this factor, which is associated with the planting of trees in villages and cities, as well as along roads [25,71]. Little is known about the factors that determine the distribution of epiphytic mosses in urban landscapes; however, in the last few decades, an increase in epiphytic species richness and the number of their localities in city centers were observed [72,73]. Recently a significant heat-island effect on bryophyte species growing on the ground or stones in historical Japanese moss gardens was reported by Oishi et al. [74]. Studies by Żołnierz et al. [75] also showed a significant but weak effect of this phenomenon on epiphytic mosses growing in the center of a large European city.
The results of our modelling showed a clear link between some genera and similar environmental conditions. A very similar affinity to habitat conditions was clearly demonstrated in the genera Leratia, Macromitrium, and Ulota. The potential distribution of the species of these genera was predicted mainly in the south and southeastern provinces of China, which are characterised by a subtropical monsoon climate [68,76]. Moreover, Yunan and Guizhou provinces are located in one of the warmest regions in China with an average daily high temperature of 29°C and high precipitation but also have a few truly tropical and sultry months [68,76]. Quite an opposite pattern in comparison to the previous three genera is clearly visible in the models of Lewinskya, Nyholmiella and Orthotrichum. These epiphytes prefer mainly temperate areas with continental climate often in medium-high to high mountain massifs, which are found both in west, south, and east China, as well as in other adjacent countries. These areas have a significantly cooler climate and the amount of precipitation fluctuates depending on the relief of the mountains. Slopes on the windward side collect more precipitation and, on the contrary, the leeward sides are in a rain shadow. In combination with sparser to open forest vegetation, where the spores of the epiphytic mosses can spread very effectively, these places are an ideal habitat for the occurrence of the above-mentioned genera.
One potential application of the distribution models is the forecasting of future species distribution and predicted species response to climate change [10]. The species growing in humid climates (e.g., in tropical cloud forests) seem to be more sensitive to climate change than species growing in regions characterised by recurring seasonal drought (e.g., species from seasonal lowland forests) [13]. According to Wierzcholska et al. [10], loss of mature forest ecosystems is a more important threat to woodland-specialist epiphytic bryophytes than climate change, which results from limited dispersal capabilities and specific habitat requirements. However, habitat fragmentation seems to have more disastrous effects on moisture-loving species, than on drought-resistant ones [13,77]. Thus, we can assume that representatives of the Orthotrichaceae family connected with dryer ecoregions (e.g., Lewinskya affinis, L. hookeri, L. speciosa, L. vladikavkana, Nyholmiella obtusifolia, Orthotrichum alpestre, O. anomalum, O. callistomum, O. crenulatum, O. dasymitrium, O. pallens, O. pamiricum, O. scanicum, O. sordidum, O. urnigerum, and Ulota robusta) will be more resistant to future climate warming than those from humid ones (e.g., Leratia exigua, Macromitrium gymnostomum, M. cavaleriei, O. consobrinum, O. crispifolium, O. griffithii, and Ulota gymnostoma). However, this phenomenon would require further research on the potential future distribution of the studied species.
Most of the Orthotrichaceae genera used in the modelling showed a strong affinity to the transition zone (ecotone), which is located vertically in the mountain areas (e.g., in the Himalayas or Hengduan Shan) and horizontally at the border of two different climate regions. This zone crosses the area of China diagonally from south to northeast, through the provinces of Yunnan, Sichuan, Shaanxi, and Heibei [31,35,36,67]. Changes that occur in population or community structures at the boundary of two or more habitats are called the ‘edge effect’ [78,79]. A consequence of the edge effect in the transition zone in China is an increase in the species diversity of epiphytic bryophytes. The conditions on the border of both climate zones favour the occurrence of both studied ecological groups (western-montane and eastern-lowland epiphytic bryophytes). Moreover, the phenomenon of a high degree of species migration within the ecotone zone has been observed [80,81], which can affect the spatial pattern of epiphytic bryophytes. The boundary of the climatic zones is therefore a clearly defined hotspot of the epiphytic species (belonging to the Orthotrichaceae family) richness in China and important areas for nature conservation purposes.

Author Contributions

Conceptualization, V.P., L.F., E.K. and M.N.; methodology, L.F., E.K., M.N. and V.P.; formal analysis, V.P., L.F., E.K., M.N., L.F. and S.G.; investigation, V.P., L.F., E.K. and M.N.; writing-original draft preparation, L.F., V.P., E.K., M.N. and S.G; writing-review and editing, L.F., V.P., E.K., M.N. and S.G.; visualization, E.K., V.P. and L.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by research projects of University of Ostrava No. SGS14/PřF/2022.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All authors agree with MDPI Research Data Policies.

Acknowledgments

We are very grateful to all the reviewers for their valuable advice and suggestions to improve the manuscript. The authors are indebted to John Campbell (Scotland) for kindly checking the English.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Overpeck, J.; Rind, D.; Goldberg, R. Climate-induced changes in forest disturbance and vegetation. Nature 1990, 343, 51–53. [Google Scholar] [CrossRef]
  2. Pearson, R.G.; Dawson, T.P. Predicting the impacts of climate change on the distribution of species: Are bioclimate envelope models useful? Glob. Ecol. Biogeogr. 2003, 12, 361–371. [Google Scholar] [CrossRef] [Green Version]
  3. Thuiller, W.; Lavergne, S.; Roquet, C.; Boulangeat, I.; Lafourcade, B.; Araujo, M.B. Consequences of climate change on the tree of life in Europe. Nature 2011, 470, 531–534. [Google Scholar] [CrossRef] [PubMed]
  4. Allen, A.P.; Brown, J.H.; Gillooly, J.F. Global biodiversity, biochemical kinetics, and the energetic-equivalence rule. Science 2002, 297, 1545–1548. [Google Scholar] [CrossRef]
  5. Currie, D.; Mittelbach, G.; Cornell, H.; Field, R.; Guégan, J.-F.; Hawkins, B.; Kaufman, D.; Kerr, J.; Oberdorff, T.; O’Brien, E.; et al. Predictions and tests of climate-based hypotheses broad-scale variation in taxonomic richness. Ecol. Lett. 2004, 7, 1121–1134. [Google Scholar] [CrossRef]
  6. Stein, A.; Gerstner, K.; Kreft, H. Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales. Ecol. Lett. 2014, 17, 866–880. [Google Scholar] [CrossRef]
  7. Rana, S.K.; Price, T.D.; Qian, H. Plant species richness across the Himalaya driven by evolutionary history and current climate. Ecosphere 2019, 10, e02945. [Google Scholar] [CrossRef] [Green Version]
  8. Vintsek, L.; Klichowska, E.; Nowak, A.; Nobis, M. Genetic differentiation, demographic history and distribution models of high alpine endemic vicariants outline the response of species to predicted climate changes in a Central Asian biodiversity hotspot. Ecol. Indic. 2022, 144, e109419. [Google Scholar] [CrossRef]
  9. Gignac, L.D. Bryophytes as Indicators of Climate Change. Bryologist 2001, 104, 410–420. [Google Scholar] [CrossRef]
  10. Wierzcholska, S.; Dyderski, M.K.; Jagodziński, A.M. Potential distribution of an epiphytic bryophyte depends on climate and forest continuity. Glob. Planet. Change 2020, 193, e103270. [Google Scholar] [CrossRef]
  11. Draper, I.; Lara, F.; Albertos, B.; Garilleti, R.; Mazimpaka, V. Epiphytic bryoflora of the Atlas and Antiatlas Mountains, including a synthesis of the distribution of epiphytic bryophytes in Morocco. J. Bryol. 2006, 28, 312–330. [Google Scholar] [CrossRef]
  12. Mazimpaka, V.; Medina, N.G.; Draper, I.; Lara, F. Epiphytic bryophyte flora in dry environments from the Western Mediterranean: The special case of Sierra Alhamilla (Almería, South-eastern Spain). Plant Biosyst. 2009, 143, 113–125. [Google Scholar] [CrossRef]
  13. Zotz, G.; Bader, M. Sampling vascular epiphyte diversity-Species richness and community structure. Ecotropica 2011, 17, 103–112. [Google Scholar]
  14. Goffinet, B.; Bayer, R.J.; Vitt, D.H. Circumscription and phylogeny of the Orthotrichales (Bryopsida) inferred from rbcL sequence analyses. Am. J. Bot. 1998, 85, 1324–1337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Goffinet, B.; Shaw, A.J.; Cox, C.J.; Wickett, N.J.; Boles, S. Phylogenetic inferences in the Orthotrichoideae (Orthotrichaceae: Bryophyta) based on variation in four loci from all genomes. Monogr. Syst. Bot. Mo. Bot. Gard. 2004, 98, 270–289. [Google Scholar]
  16. Sawicki, J.; Plášek, V.; Szczecińska, M. Preliminary studies on the phylogeny of Orthotrichum (Bryophyta) inferred from nuclear ITS sequences. Ann. Bot. Fenn. 2009, 46, 507–515. [Google Scholar] [CrossRef]
  17. Sawicki, J.; Plášek, V.; Szczecińska, M. Molecular studies resolve Nyholmiella (Orthotrichaceae) as a separate genus. J. Syst. Evol. 2010, 48, 183–194. [Google Scholar] [CrossRef]
  18. Sawicki, J.; Plášek, V.; Szczecińska, M. Molecular data do not support the current division of Orthotrichum (Bryophyta) species with immersed stomata. J. Syst. Evol. 2012, 50, 12–24. [Google Scholar] [CrossRef]
  19. Sawicki, J.; Plášek, V.; Ochyra, R.; Szczecińska, M.; Ślipiko, M.; Myszczyński, K.; Kulik, T. Mitogenomic analyses support the recent division of the genus Orthotrichum (Orthotrichaceae, Bryophyta). Sci. Rep. 2017, 7, e4408. [Google Scholar] [CrossRef] [Green Version]
  20. Plášek, V.; Komínková, Z.; Číhal, L.; Guo, S.-L. A noteworthy disjunction of the epiphytic moss Lewinskya graphiomitria (Orthotrichaceae) to China. Acta Soc. Bot. Pol. 2020, 83, 45–53. [Google Scholar] [CrossRef]
  21. Draper, I.; Garilleti, R.; Calleja, J.A.; Flagmeier, M.; Mazimpaka, V.; Vigalondo, B.; Lara, F. Insights into the evolutionary history of the Subfamily Orthotrichoideae (Orthotrichaceae, Bryophyta): New and former supra-specific taxa so far obscured by prevailing homoplasy. Front. Plant Sci. 2021, 12, 427. [Google Scholar] [CrossRef] [PubMed]
  22. Draper, I.; Villaverde, T.; Garilleti, R.; Burleigh, J.G.; McDaniel, S.F.; Mazimpaka, V.; Calleja, J.A.; Lara, F. An NGS-Based Phylogeny of Orthotricheae (Orthotrichaceae, Bryophyta) With the Proposal of the New Genus Rehubryum from Zealandia. Front. Plant Sci. 2022, 13, 882960. [Google Scholar] [CrossRef] [PubMed]
  23. Ma, Y.; Li, D.; Yu, J.; Guo, S. Geographical distribution patterns of Macromitrium and Orthotrichum in China and their relationship with climatic factors. Biodivers. Sci. 2013, 21, 177–184. [Google Scholar] [CrossRef]
  24. Shen, Y.; Yu, J.; Guo, S.L. Macromitrium and Orthotrichum distribution patterns under different climate warming scenarios in China. Acta Ecol. Sin. 2015, 35, 6449–6459. [Google Scholar]
  25. Nowak, A.; Plášek, V.; Nobis, M.; Nowak, S. Epiphytic communities of open habitats in the western Tian-Shan Mts (Middle Asia: Kyrgyzstan). Cryptogam. Bryol. 2016, 37, 415–433. [Google Scholar] [CrossRef]
  26. Plášek, V.; Ochyra, R. Orthotrichum alpestre, a New Addition to the Moss Flora of Poland, with Notes on O. schimperi (Orthotrichaceae: Bryophyta). Acta Soc. Bot. Pol. 2020, 89, e89310. [Google Scholar] [CrossRef]
  27. Medina, N.G.; Albertos, B.; Lara, F.; Mazimpaka, V.; Garilleti, R.; Draper, D.; Hortal, J. Species richness of epiphytic bryophytes: Drivers across scales on the edge of the Mediterranean. Ecography 2014, 37, 80–93. [Google Scholar] [CrossRef]
  28. Medina, N.G.; Bowker, M.A.; Hortal, J.; Mazimpaka, V.; Lara, F. Shifts in the importance of the species pool and environmental controls of epiphytic bryophyte richness across multiple scales. Oecologia 2018, 186, 805–816. [Google Scholar] [CrossRef]
  29. Wang, Q.H.; Jia, Y. A taxonomic revision of the Asian species of Ulota Mohr (Orthotrichaceae). Bryologist 2012, 115, 412–443. [Google Scholar] [CrossRef]
  30. Lou, Y.X.; Si, H.; Guo, S.L. Using macroclimatic models to estimate the distribution ranges of taxonomically challenging taxa, an example with Macromitrium cavaleriei Cardot & Ther. (Orthotrichaceae). J. Bryol. 2014, 36, 271–279. [Google Scholar] [CrossRef]
  31. Beck, H.E.; Zimmermann, N.E.; McVicar, T.R.; Vergopolan, N.; Berg, A.; Wood, E.F. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 2018, 5, 180214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Skoupá, Z.; Ochyra, R.; Guo, S.L.; Sulayman, M.; Plášek, V. Distributional novelties for Lewinskya, Nyholmiella and Orthotrichum (Orthotrichaceae) in China. Herzogia 2017, 30, 58–73. [Google Scholar] [CrossRef]
  33. Skoupá, Z.; Ochyra, R.; Guo, S.L.; Sulayman, M.; Plášek, V. Three remarkable additions of Orthotrichum species (Orthotrichaceae) to the moss flora of China. Herzogia 2018, 31, 88–100. [Google Scholar] [CrossRef]
  34. Plášek, V.; Komínková, Z.; Ochyra, R.; Fialová, L.; Guo, S.; Sulayman, M. A Synopsis of Orthotrichum s. lato (Bryophyta, Orthotrichaceae) in China, with Distribution Maps and a Key to Determination. Plants 2021, 10, 499. [Google Scholar] [CrossRef]
  35. Shi, Y.; Wang, G. Changes in building climate zones over China based on high-resolution regional climate projections. Environ. Res. Lett. 2020, 15, e114045. [Google Scholar] [CrossRef]
  36. Lin, G.; Fu, J.; Jiang, D.; Wang, J.; Wang, Q.; Dong, D. Spatial variation of the relationship between PM2. 5 concentrations and meteorological parameters in China. BioMed Res. Int. 2015, 2015, e84618. [Google Scholar] [CrossRef]
  37. Peterson, A.; Soberón, J. Species Distribution Modeling and Ecological Niche Modeling: Getting the Concepts Right. Nat. Conserv. 2012, 10, 1–6. [Google Scholar] [CrossRef]
  38. Jia, Y.; He, S.; Guo, S.L. Orthotrichaceae. In Moss Flora of China: Erpodiaceae–Climaciaceae; Wu, P.-C., Crosby, M.R., He, S., Eds.; Missouri Botanical Garden Press: St. Louis, MO, USA, 2011; pp. 22–117. [Google Scholar]
  39. Plášek, V.; Sawicki, J.; Ochyra, R.; Szczecińska, M.; Kulik, T. New taxonomical arrangement of the traditionally conceived genera Orthotrichum and Ulota (Orthotrichaceae, Bryophyta). Acta Musei Sil. Sci. Nat. 2015, 64, 169–174. [Google Scholar] [CrossRef] [Green Version]
  40. Plášek, V.; Sawicki, J.; Ochyra, R. Are Orthotrichum and Dorcadion (Orthotrichaceae, Bryophyta) heterotypic generic names? Acta Musei Sil. Sci. Nat. 2016, 65, 193–202. [Google Scholar] [CrossRef] [Green Version]
  41. Lara, F.; Garilleti, R.; Goffinet, B.; Draper, I.; Medina, R.; Vigalondo, B.; Mazimpaka, V. Lewinskya, a new genus to accommodate the phaneroporous and monoicous taxa of Orthotrichum (Bryophyta, Orthotrichaceae). Cryptogam. Bryol. 2016, 37, 361–382. [Google Scholar] [CrossRef]
  42. Crosby, M.R.; Magill, R.E.; Allen, B.; He, S. A Checklist of the Mosses. Missouri Botanical Garden, St. Louis, Missouri. Available online: http://www.mobot.org/MOBOT/tropicos/most/checklist.shtml (accessed on 4 September 2022).
  43. Frey, W.; Stech, M. Marchantiophyta, Bryophyta, Anthocerotophyta. In Syllabus of Plant Families. A. Engler´s Syllabus der Pflanzenfamilien; Frey, W., Ed.; Gebr. Borntraeger: Stuttgart, Germany, 2009; pp. 1–419. [Google Scholar]
  44. Vitt, D.H.; Ramsay, H.P. The Macromitrium complex in Australasia (Orthotrichaceae, Bryopsida). Part I. Taxonomy and phylogenetic relationships. J. Hattori Bot. Lab. 1985, 59, 325–451. [Google Scholar]
  45. Vitt, D.H. The infrageneric evolution, phylogeny, and taxonomy of the genus Orthotrichum (Musci) in North America. Nova Hedwig. 1971, 21, 683–711. [Google Scholar]
  46. Lewinsky, J. The genus Orthotrichum. Morphological studies and evolutionary remarks. J. Hattori Bot. Lab. 1977, 3, 31–61. [Google Scholar]
  47. Otnyukova, T.N. A new species of Orthotrichum (Orthotrichaceae, Musci) from Tuva Republic (South Siberia). Arctoa 2001, 10, 155–156. [Google Scholar] [CrossRef] [Green Version]
  48. Gao, C.; Lai, M.-J. Illustrations of Bryophytes of China; SMC Publishing Inc.: Taipei, China, 2003; pp. 1–1313. [Google Scholar]
  49. Wang, Q.; Jia, Y. A taxonomic revision of Ulota Mohr (Orthotrichaceae) in South and Central America. Plant Divers. 2016, 38, 65–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Zizka, A.; Azevedo, J.; Leme, E.; Neves, B.; Ferreira da Costa, A.; Caceres, D.; Zizka, G. Biogeography and conservation status of the pineapple family (Bromeliaceae). Divers. Distrib. 2019, 18, 183–195. [Google Scholar] [CrossRef]
  51. Hijmans, R.J.; Phillips, S.; Leathwick, J.; Elith, J. DISMO: Species Distribution Modeling. R Package Version 1.1-4. 2017. Available online: https://CRAN.R-project.org/package=dismo (accessed on 4 September 2022).
  52. R Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. 2020. Available online: http://www.R-project.org (accessed on 4 September 2022).
  53. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2020. Available online: http://qgis.osgeo.org (accessed on 4 September 2022).
  54. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  55. Fick, S.E.; Hijmans, R.J. Worldclim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climat. 2017, 32, 4302–4315. [Google Scholar] [CrossRef]
  56. Tuanmu, M.N.; Jetz, W. A global 1-km consensus land-cover product for biodiversity and ecosystem modeling. Global Ecol. Biogeogr. 2014, 23, 1031–1045. [Google Scholar] [CrossRef]
  57. Venter, O.; Sanderson, E.W.; Magrach, A.; Allan, J.R.; Beher, J.; Jones, K.R.; Possingham, H.P.; Laurance, W.F.; Wood, P.; Fekete, B.M.; et al. Last of the Wild Project Version 3 (LWP-3): 2009 Human Footprint; NASA Socioeconomic Data and Applications Center (SEDAC): New York, NY, USA, 2018. [Google Scholar] [CrossRef]
  58. Rodríguez-Castañeda, G.; Hof, A.R.; Jansson, R.; Harding, L.E. Predicting the Fate of Biodiversity Using Species’ Distribution Models: Enhancing Model Comparability and Repeatability. PLoS ONE 2012, 7, e44402. [Google Scholar] [CrossRef] [Green Version]
  59. Phillips, S.J.; Dudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  60. ESRI, Inc. Cell Size and Resampling in Analysis. Available online: https://desktop.arcgis.com/en/arcmap/latest/extensions/spatial-analyst/performing-analysis/cell-size-and-resampling-in-analysis.html (accessed on 4 September 2022).
  61. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; García, M.J.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 2013, 36, 27–46. [Google Scholar] [CrossRef]
  62. TIBCO, STATISTICA (Data Analysis Software System). 2017. Available online: https://docs.tibco.com/products/tibco-statistica/archive (accessed on 4 September 2022).
  63. Baldwin, R.A. Use of maximum entropy modelling in wildlife research. Entropy 2009, 11, 854–866. [Google Scholar] [CrossRef]
  64. Brown, J.; Bennett, J.; French, C. SDMtoolbox 2.0: The next generation Python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. PeerJ 2017, 5, e4095. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Janžekovič, F.; Novak, T. PCA—A Powerful Method for Analyze Ecological Niches. In Principal Component Analysis—Multidisciplinary Applications; Sanguansat, P., Ed.; IntechOpen: London, UK, 2012. [Google Scholar] [CrossRef]
  66. Escoriza, D.; Ben Hassine, J.; Boix, D. Factors regulating the invasive success of an alien frog: A comparison of the ecology of the native and alien populations. Hydrobiologia 2014, 30, 127–138. [Google Scholar] [CrossRef]
  67. Hsin-Shih, C.; Dian-an, Y. A study on climate-vegetation interaction in China: The ecological model for global change. Coenoses 1993, 8, 105–119. [Google Scholar]
  68. China Meteorological Administration (CHMA). Climate. CHMA: Beijing, China. 2022. Available online: https://www.cma.gov.cn/en2014/climate/shtml (accessed on 4 September 2022).
  69. Proctor, M.C.F. Climate responses and limits of bryophytes: Comparisons and contracts with vascular plants. In Bryophyte Ecology and Climate Change; Slack, N., Tuba, Z., Stark, L.R., Eds.; Cambridge University Press: Cambridge, UK, 2011; pp. 35–54. [Google Scholar]
  70. Ódor, P.; Király, I.; Tinya, F.; Bortignon, F.; Nascimbene, J. Patterns and drivers of species composition of epiphytic bryophytes and lichens in managed temperate forests. For. Ecol. Manag. 2013, 306, 256–265. [Google Scholar] [CrossRef] [Green Version]
  71. Plášek, V.; Sawicki, J.; Číhal, L. Orthotrichum pamiricum (Bryophyta), a new epiphytic moss species from Pamir Mountains in Central Asia. Turk. J. Bot. 2014, 38, 754–762. [Google Scholar] [CrossRef]
  72. Fudali, E. Recent tendiences in distribution of epiphytic bryophytes in urban areas: A Wroclaw case study (south-west Poland). Polish Bot. J. 2012, 57, 231–241. [Google Scholar]
  73. Stebel, A.; Fojcik, B. Changes in the epiphytic bryophyte flora in Katowice city (Poland). Cryptogam. Bryol. 2016, 37, 399–414. [Google Scholar] [CrossRef]
  74. Oishi, Y. Urban heat island effects on moss gardens in Kyoto, Japan. Landsc. Ecol. Eng. 2019, 15, 177–184. [Google Scholar] [CrossRef]
  75. Żołnierz, L.; Fudali, E.; Szymanowski, M. Epiphytic Bryophytes in an Urban Landscape: Which Factors Determine Their Distribution, Species Richness, and Diversity? A Case Study in Wroclaw, Poland. Int. J. Environ. Res. Public. Health 2022, 19, 6274. [Google Scholar] [CrossRef] [PubMed]
  76. Zheng, J.; Yin, Y.; Li, B. A new scheme for climate regionalization in China. Acta Geogr. Sin. 2010, 65, 3–13. [Google Scholar]
  77. Číhal, L.; Kaláb, O.; Plášek, V. Modeling the distribution of rare and interesting moss species of the family Orthotrichaceae (Bryophyta) in Tajikistan and Kyrgyzstan. Acta Soc. Bot. Pol. 2017, 86, e3543. [Google Scholar] [CrossRef] [Green Version]
  78. Lourenço, G.M.; Soares, G.R.; Santos, T.P.; Dáttilo, W.; Freitas, A.V.; Ribeiro, S.P. Equal but Different: Natural Ecotones Are Dissimilar to Anthropic Edges. PLoS ONE 2019, 14, e0213008. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Malanson, G.P.; Resler, L.M.; Tomback, D.F. Ecotone Response to Climatic Variability Depends on Stress Gradient Interactions. Clim. Chang. Resp. 2017, 4, 1. [Google Scholar] [CrossRef] [Green Version]
  80. Holland, M.M.; Risser, P.G.; Naiman, R.J. Ecotones: The Role of Landscape Boundaries in the Management and Restoration of Changing Environments; Chapman and Hall: New York, NY, USA, 1991; pp. 1–165. [Google Scholar]
  81. Kent, M.; Moyeed, R.A.; Reid, C.L.; Pakeman, R.; Weaver, R. Geostatistics, spatial rate of change analysis and boundary detection in plant ecology and biogeography. Prog. Phys. Geogr. 2006, 30, 201–231. [Google Scholar] [CrossRef]
Figure 1. Model of the potential distribution of the genus Leratia in China and adjacent areas.
Figure 1. Model of the potential distribution of the genus Leratia in China and adjacent areas.
Plants 12 00222 g001
Figure 2. Model of the potential distribution of the genus Lewinskya in China and adjacent areas.
Figure 2. Model of the potential distribution of the genus Lewinskya in China and adjacent areas.
Plants 12 00222 g002
Figure 3. Model of the potential distribution of the genus Macromitrium in China and adjacent areas.
Figure 3. Model of the potential distribution of the genus Macromitrium in China and adjacent areas.
Plants 12 00222 g003
Figure 4. Model of the potential distribution of the genus Nyholmiella in China and adjacent areas.
Figure 4. Model of the potential distribution of the genus Nyholmiella in China and adjacent areas.
Plants 12 00222 g004
Figure 5. Model of the potential distribution of the genus Orthotrichum in China and adjacent areas.
Figure 5. Model of the potential distribution of the genus Orthotrichum in China and adjacent areas.
Plants 12 00222 g005
Figure 6. Model of the potential distribution of the genus Ulota in China and adjacent areas.
Figure 6. Model of the potential distribution of the genus Ulota in China and adjacent areas.
Plants 12 00222 g006
Figure 7. Comparison of spatial distributions and ecological niches of two studied epiphytic groups: (A) map of predicted potential distributions; (B) percent contribution of environmental variables used for ENMs; (C) principal component analysis (PCA) plot representing the realised niches occupied by two studied groups based on climatic variables. The ellipses represent the 95% confidence intervals; (D) correlations between variables and the first two principal components. For variables abbreviation see Table 1.
Figure 7. Comparison of spatial distributions and ecological niches of two studied epiphytic groups: (A) map of predicted potential distributions; (B) percent contribution of environmental variables used for ENMs; (C) principal component analysis (PCA) plot representing the realised niches occupied by two studied groups based on climatic variables. The ellipses represent the 95% confidence intervals; (D) correlations between variables and the first two principal components. For variables abbreviation see Table 1.
Plants 12 00222 g007
Figure 8. Models of potential distributions of two studied groups: (A) western-montane epiphytic group; (B) eastern-lowland epiphytic group.
Figure 8. Models of potential distributions of two studied groups: (A) western-montane epiphytic group; (B) eastern-lowland epiphytic group.
Plants 12 00222 g008
Table 1. A list of environmental layers. The uncorrelated variables used for modelling are marked in bold.
Table 1. A list of environmental layers. The uncorrelated variables used for modelling are marked in bold.
AbbreviationName/ExplanationUnit
BIO1annual mean temperaturedegrees Celsius × 10
BIO2mean diurnal range (mean of monthly max temp—min temp)degrees Celsius × 10
BIO3isothermality (bio2/bio7) (×100)degrees Celsius × 10
BIO4temperature seasonality (standard deviation ×100)degrees Celsius × 10
BIO5max temperature of warmest monthdegrees Celsius × 10
BIO6min temperature of coldest monthdegrees Celsius × 10
BIO7temperature annual range (bio5-bio6)degrees Celsius × 10
BIO8mean temperature of wettest quarterdegrees Celsius × 10
BIO9mean temperature of driest quarterdegrees Celsius × 10
BIO10mean temperature of warmest quarterdegrees Celsius × 10
BIO11mean temperature of coldest quarterdegrees Celsius × 10
BIO12annual precipitationmillimeters
BIO13precipitation of wettest monthmillimeters
BIO14precipitation of driest monthmillimeters
BIO15precipitation seasonality (coefficient of variation)percentage
BIO16precipitation of wettest quartermillimeters
BIO17precipitation of driest quartermillimeters
BIO18precipitation of warmest quartermillimeters
BIO19precipitation of coldest quartermillimeters
hfphuman footprintpercentage
cons1evergreen/deciduous needle leaf treespercentage
cons2evergreen broadleaf treespercentage
cons3deciduous broadleaf treespercentage
cons4mixed/other treespercentage
cons5shrubspercentage
Table 2. Percentage contribution of variables used for distribution modelling and performance of models measured using AUC of six genera and two studied epiphytic groups. The three highest values of the variables for each model are marked in bold.
Table 2. Percentage contribution of variables used for distribution modelling and performance of models measured using AUC of six genera and two studied epiphytic groups. The three highest values of the variables for each model are marked in bold.
CodeLeratiaLewinskyaMacromitriumNyholmiellaOrthotrichumUlotawestern-montaneeastern-lowland
bio102.6040.81.439.733.99.841.32.6
bio153.3023.13.514.534.61.825.38.5
bio1811.302.327.37.95.115.84.134.7
bio1916.901.71.141.52.81.50.9
bio24.405.631.24.73.15.84.19.7
bio40.206.19.43.95.85.46.117.8
bio95.801.115.80.30.723.30.918.6
cons13.407.40.74.51.24.55.60.6
cons21.400.33.20.20.21.60w.10.5
cons30.100.20.53.113.10.20.2
cons448.6011.12.82.525.50.53.2
cons51.800.81.62.10.60.10.70.5
hfp0.209.63.312.5100.59.72.1
AUC0.9280.9280.9500.9250.9080.9250.9030.944
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

Fialová, L.; Plášek, V.; Klichowska, E.; Guo, S.; Nobis, M. Temperature and Precipitation More Than Tree Cover Affect the Distribution Patterns of Epiphytic Mosses within the Orthotrichaceae Family in China and Adjacent Areas. Plants 2023, 12, 222. https://doi.org/10.3390/plants12010222

AMA Style

Fialová L, Plášek V, Klichowska E, Guo S, Nobis M. Temperature and Precipitation More Than Tree Cover Affect the Distribution Patterns of Epiphytic Mosses within the Orthotrichaceae Family in China and Adjacent Areas. Plants. 2023; 12(1):222. https://doi.org/10.3390/plants12010222

Chicago/Turabian Style

Fialová, Lucie, Vítězslav Plášek, Ewelina Klichowska, Shuiliang Guo, and Marcin Nobis. 2023. "Temperature and Precipitation More Than Tree Cover Affect the Distribution Patterns of Epiphytic Mosses within the Orthotrichaceae Family in China and Adjacent Areas" Plants 12, no. 1: 222. https://doi.org/10.3390/plants12010222

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