Honglin Li1,2
Peng Luo- 1CAS Key Laboratory of Mountain Ecological Restoration and Bioresource Utilization & Ecological Restoration and Biodiversity Conservation Key Laboratory of Sichuan Province, Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, China
- 2College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
A comprehensive understanding of the effects of mountain roads on plant diversity is critical to finding the most effective solutions for managing this particular driver. Little is known, however, about the simultaneous effects that road have on the multiple facets of biodiversity, although roads are considered to be one of the major disturbances in the Qionglai mountain range. In this study, we analyzed the impact of roads on the multiple facets of plant diversity (taxonomic, functional and phylogenetic diversity) in the study area using Hill numbers by comparing plant diversity between roadside and interior plots at the landscape scale, then, we used linear mixed models to analyze the effect of mountain roads on the multiple facets of plant diversity along an elevational gradient. The results showed that the roadside plots lacked 29.45% of the total number of species with particular functional traits (such as a relatively high specific leaf area (SLA), a relatively low leaf dry matter content (LDMC) and relatively old clades) and exclusively contained 14.62% of the total number of species. Compared with the interior community, the taxonomic, functional and phylogenetic diversity of roadside community decreased by no more than 26.78%, 24.90% and 16.62%, respectively. Taxonomic and functional diversity of dominant and common species showed greater changes to road disturbances, while rare species showed the greatest change in phylogenetic diversity. Taxonomic homogenization of roadside communities was accompanied by functional and phylogenetic homogenization. Additionally, the impact of roads on these three facets of plant diversity showed the characteristics of peak clipping along the elevation gradient. Our findings highlight the negative impact of roads on the taxonomic, functional and phylogenetic diversity of the Qionglai mountain range, as roads promote communities that are more similar in taxonomic, functional, and phylogenetic composition, and to a greater extent contributed to compositional evenness. These effects tend to be functionally and phylogenetically non-random, and species in some clades or with some functional traits are at higher risk of loss. Our results are important for the conservation and management of nature reserves, especially for local governments aiming to create new infrastructure to connect natural mountainous areas.
Introduction
Although covering a relatively small fraction of Earth’s area (13%~25%, Ives, 1997), mountainous regions support a significant one-third of all terrestrial biological diversity due to their topographic and environmental complexity (Tang et al., 2006; Perrigo et al., 2020; Brighenti et al., 2021). A common feature of many mountains is that they contain road networks that connect lowlands with high-elevation sites. Economic development will mandate the construction of many new roads in the coming decades (Seto et al., 2012; Xie et al., 2021), and these roads will extend into once-inaccessible mountainous areas; this will cause these regions to face trade-offs between biodiversity conservation and economic growth (Miller et al., 2021). Mountain roads can drive various changes in the global mountain community composition and structure via habitat losses (Laurance et al., 2009), fragmentation (Delgado et al., 2007) and pollution from combustion of fossil fuels by the on-road vehicles (Karmakar et al., 2021), all of which affect biodiversity and impact the structure and functioning of ecosystems (Millennium Ecosystem Assessment, 2005; De Pauw et al., 2021).
During the past two decades, research into the effects of roads on plant diversity in mountainous regions has gained momentum (Pauchard and Alaback, 2004; Barbosa et al., 2010; Seipel et al., 2012). Previous studies have focused on the impact of roads on plant communities by assisting the dispersal of non-native species or native generalists, and the results suggest that mountain roads and the colonization of generalists lead to taxonomic homogenization of plant community (Arévalo et al., 2010; Haider et al., 2018). Most of what we currently know about the influence of roads on plant diversity across the globe has come from studies of plant species richness and composition (Lembrechts et al., 2017; Dar et al., 2018; Rashid et al., 2021). However, no studies have investigated how the functional traits and phylogenetic clades respond to environmental filtering and biological competition induced by mountain roads, thus limiting our ability to devise practical strategies for conserving plant diversity.
The Qionglai mountain range, located on the easternmost edge of the Hengduan Mountains in Sichuan Province, is home to many temperate vascular plant species, some of which are ancient relic and endemic species (Myers et al., 2000; Boufford, 2014). This range was recognized as a World Heritage Site in 2006 due to its unique and valuable plant diversity (Boufford, 2014). Climate change and anthropogenic disturbances are identified as the two major reasons driving recent changes in plant diversity and community composition in this region (Li et al., 2013; Qiu et al., 2014; Salick et al., 2019; Yin et al., 2020). Conservation efforts have focused on managing the human activities that affect plant diversity (Kang, 2021; State Council Information Office of China (SCIO), 2015). With the official establishment of the Giant Panda National Park (GPNP) in 2020, clear policies and guidelines for human activities such as enforcement, tourism and infrastructure development guidelines are needed within this region. A comprehensive study of the impact of roads on the plant diversity would provide managers with enough information for them to better predict the gains and losses associated with various decisions when planning new roads or closing existing roads.
Roads are gathering points of human activities, carrying material and cultural exchanges between lowlands and high-elevation sites. In this study area, Yin et al. (2020) found that road was one of the major human activities by quantifying the relative contributions of climate change and human activities to net primary productivity and Kang (2021) concluded by reviewing the impact of anthropogenic disturbance on the habitat of giant pandas that roads were the major factor affecting the habitat of giant pandas. In addition, the State Forestry Administration also believed that roads are a major human disturbance in the study area (State Forestry Administration (SFA), 2015). Nonetheless, it is unclear whether functional and phylogenetic diversity change as a result of road disturbances, whether such changes are in the same direction as taxonomic diversity (i.e., decreased diversity) and, if so, with what intensity these changes occur.
Moreover, biodiversity includes not only the number of species but also aspects such as identity, dominance, evenness and rarity (Hill et al., 2016; Hillebrand et al., 2018). Thus, measuring plant species richness alone may not be sufficient to represent changes in plant diversity (Berdugo et al., 2019). Within communities, species composition changes may have effects on ecosystem functions that are equal to, or stronger than, the effects of changes in species richness (Shi et al., 2016). Thus, assessing road-driven shifts in species abundance is paramount for understanding how these compositional changes affect community and ecosystem functioning in response to human disturbances.
In this study, we investigated whether and how mountain roads affected the taxonomic, functional and phylogenetic diversity of plant communities in the Qionglai mountain range, and whether the responses of plant communities to roads are dependent on the location of the plant communities on the elevation gradient. Specifically, we aimed to answer the following questions: (i) Have there been changes in the three facets of plant diversity due to mountain roads by comparing the pattern of species diversity between roadside plots and interior plots, and if so, to what extent? (ii) Are the shifts in plant diversity synchronous and monodirectional; i.e., does plant diversity uniformly increase or decline along the elevation gradient? For the first question, we expect that species conserved from interior plots or spread along roads may be closely related (even sister species) or functionally similar species. Therefore, roads lead to decreases in overall plant diversity indices by inducing changes in species richness and abundance (Haddad et al., 2015; Razafindratsima et al., 2018). Furthermore, we predict that changes in functional and phylogenetic diversity will be smaller than changes in taxonomic diversity due to functional and phylogenetic redundancy. In relation to our second question, we hypothesize that the pattern of the effects of roads on the functional and phylogenetic diversity indices would be consistent with the effect on taxonomic diversity, i.e., the effects would vary with the elevation gradient. Based on our previous results (Li et al., 2022), the expected negative impacts of roads on diversity may be greater at moderate and high elevations than at low elevations. Answering the abovementioned questions and filling this knowledge gap are important, as changes in the three facets of plant diversity may have different implications for the functioning of ecosystems and their responses to human activities.
Methods
Study area
The Qionglai mountain range is located in the upper Yangtze River of northwest Sichuan Province (Figure 1), China (102°16′ E-104°10′ E, 29°49′ N-31°31′ N). The area comprises a major portion of Giant Panda National Park. The mountain range is situated in the transitional zone between a subtropical monsoon climate zone and a continental plateau climate zone. Its highest peak, Siguniang Mountain, has an elevation of approximately 6250 m a.s.l. The landscape along the elevation gradient is dominated by arid shrubs, subtropical evergreen and deciduous broadleaf forests, coniferous forests, subalpine coniferous forests, scrub meadows, alpine talus vegetation, and permanent snow belts (Sichuan Flora Editorial Committee, 1981). Of these, arid shrubs and permanent snow belts were not included in this study (Li et al., 2022).
Figure 1 Location map of the study area.
G350, G318, G213 and G351 are four major roads in the Qionglai mountain range (Xu et al., 2006; Li et al., 2022). The total length of major roads in the mountain range is about 799.84 km, and the road density is about 0.05 km/km2. Only part of three major roads, G350, G318 and G351, are located in the national park (Figure 1) and these road segments were selected for this study. The density of major roads in the national park is about 0.02 km/km2. All road segments selected in this study comply with the third-level highway technical standards in the Technical Standards of Highway Engineering (JTG B01-2014) issued by the Ministry of Transport of China. These road segments have been in place for more than 50 years, with an average daily traffic volume of more than 1,500 vehicles. Other roads in the area are low-traffic country trails. That is to say, the selected road segments have the same properties, and they collectively represent the overall road conditions within the national park in the Qionglai mountain range.
Floristic data
In this work, field surveys were conducted along the selected segments of three national roads (G318, G350 and G351) in the Qionglai mountain range. For each of the three roads, the lowest and highest sampling sites were first determined. The lowest sampling site for each road was the location that can be sampled with the greatest avoidance of interference from other human activities. The highest sampling site of each road was more than 1,000 meters from the tunnel entrance or high passes. The sampling range of this study includes four vegetation types. To ensure adequate replicate sampling for each vegetation type, referring to the research by Zhang et al. (2009), a sampling site was systematically located at 100 m a.s.l. elevation intervals from the lowest sample site to the highest sample site (see Figure 2). In order to avoid interference from other human activities to the greatest extent, all selected sampling plots were required to be free of farmland or people gathering points within 1,000 m. In addition, sampling sites avoided pasture paths and sites with large amounts of livestock manure. In order to meet the sampling conditions described above, coupled with the difficulties of topography and access, 38 sampling sites were finally established in this study [see also Li et al. (2022)]. Four vegetation types including subtropical evergreen and deciduous broadleaf forests, coniferous forests, subalpine coniferous forests, scrub meadows included 8, 6, 15 and 9 sampling sites, respectively.
Figure 2 Schematic diagram of sampling site and plot setting.Data analysis.
At each sampling site, two 2 m*50 m plots parallel to the road were established, which were called a pair of plots (interior and roadside plots) (see Figure 2). One plot was parallel and next to the roadside was collectively referred to as the roadside plot (representing the communities affected by road). According to the previous study results on the impact of roads on adjacent ecosystems (Hansen and Clevenger, 2005; Arévalo et al., 2010; Deljouei et al., 2018), and considering the accessibility of plots during sampling, another one plot was placed parallel to the roadside plot and was more than 50 m away from the roadside (representing the natural communities away from the road). In the study, we established a total of 38 plot pairs and 76 plots. We recorded the species richness of vascular plants and the abundance of each species in each plot. In the field, we taxonomically identified the vascular plants based on the Flora of China (http://www.eFloras.org). Take photos and collect specimens of uncertain or unknown plants and bring them back to the laboratory to continue searching for identification information. Standardized nomenclature for all species according to the Plant List (http://www.theplantlist.org/).
Functional traits
For each of the 978 species in the 76 plots, we recorded 12 functional traits (Table S3). These heterogeneous characteristics collectively comprise the functional facet of plant diversity required for the calculation of Hill numbers. Data on 3 quantitative functional traits, specific leaf area (SLA), leaf dry matter content (LDMC), and leaf thickness (Lth), were obtained from field measurements. The SLA was calculated as the ratio of the fresh leaf area to the leaf dry matter content. An HPG3110 scanner (Yaxin-1241, Beijing Yaxin Liyi Technology Co., Ltd., Beijing) was used to measure the fresh leaf area. The LDMC was obtained by storing the collected leaves in silica gel and taking them back to the indoor laboratory, and then weighing the leaves after oven-drying at 60°C for at least 48 hours. The Lth was obtained by measuring the thickness of fresh leaves with Digital Vernier caliper with 0.01 mm accuracy (PD-151, Prokit’s Industries Co., Ltd., China). All traits were measured on five to 30 healthy mature individual plants, depending on the trait and species, following Pérez-Harguindeguy et al. (2013). The data of the remaining 9 qualitative functional traits were obtained through the TRY database (Kattge et al., 2020; www.try-db.org) and Flora of China (English edition; www.efloras.org). A functional distance matrix was then built based on the calculation of the Gower distance (Gower, 1971) to allow for the mixing of categorical, ordinal, and continuous data.
To analyze the specific functional compositions of the roadside and interior communities, we first used the functcomp function of the R package FD (Laliberté and Legendre, 2010) to calculate the community-weighted means (hereafter referred to as CWMs) and non-weighted means (hereafter referred to as CMs) of all functional traits based on the abundances of species in each community. Nonmetric multidimensional scaling (NMDS) based on these community-type-level CWM and CM values of the 76 plots was performed using the metaMDS function in the vegan package (Oksanen et al., 2020) to identify the functional determine differences among all plots.
Phylogenetic tree construction
To create a phylogenetic tree for each plot type, we pruned our species list from the mega-tree GBOTB.extended.tre, which currently represents the largest available phylogenetic treatment of vascular plants (Zanne et al., 2014; Smith and Brown, 2018). This was done by using the V.PhyloMaker package with the phylo.maker function and using scenario 3 method to add plant species to the phylogenetic tree (Jin and Qian, 2019).
Taxonomic, functional and phylogenetic facets of plant diversity
The three plant diversity facets were converted to Hill numbers expressed as the effective numbers of species (Jost, 2007; Chao et al., 2014; Chiu et al., 2014). We calculated all plant diversity index facets at the plot level (n = 38 for each plot type), using the “hillR” package (Li, 2018). Then, we examined the pairwise correlations of the multifaceted diversity indices using Pearson’s correlation coefficient.
When q = 0, the diversity index equals the species richness, giving disproportionate weight to rare plants (Jost, 2006; Tuomisto, 2010; Arroyo-Rodríguez et al., 2013; Vega-Álvarez et al., 2019). When q = 1, the diversity index is equivalent to Shannon’s diversity index, and each species is weighted by its abundance in the community (Jost, 2007). Finally, for q = 2 (representing the inverse Simpson concentration), dominant species are favored while rare species are discounted (Jost, 2010; Chao et al., 2012). The functional and phylogenetic counterparts of taxonomic diversity for q = [0, 1, 2] are presented in Chao et al. (2014).
Statistical analysis
To explore the overall compositional differences between the roadside and interior communities, we also evaluated the taxonomic, functional and phylogenetic dissimilarities using normalized dissimilarity measures based on the Hill numbers derived (Chao et al., 2014; Chiu et al., 2014) at the local (i.e., comparing each individual plot to other plots surveyed within the same plot type) and regional (i.e., comparing each individual plot to the pool of plots surveyed within the same plot type) scales. Then, we used linear mixed models to analyze the effects of mountain roads on three facets of plant diversity (taxonomic, functional and phylogenetic diversity) for Hill numbers q = 0, q = 1, and q = 2 along the elevation gradient. The plot type (a factor with two levels: roadside vs. interior plot), elevation and its quadratic term, and interaction between the covariate and factor were used as fixed terms in the models. Due to the nested structure of the data, linear mixed models with the site (38 levels corresponding to 38 roadside and interior plot pairs) set as a random effect (random intercept) nested within roads (three levels corresponding to three roads) were used (thus, 76 plots were nested in 38 sites nested in three roads) (Li et al., 2022). We built all possible fixed factor models and performed model selection based on the Akaike information criterion corrected for small sample sizes (AICc). We tested all possible models with the dredge function from the package MuMIn (Bartoń, 2022). The mixed effects models were fitted with the function lme in the nlme package (Pinheiro et al., 2021).
To reduce the model selection bias effects on the regression coefficient estimates in these selected subsets (Burnham & Anderson, 2002), we calculated the linear mixed model estimate values and their standard errors (SEs) for the entire set of plausible models (ΔAICc< 2) and for marginal() and conditional ( ) R2 values. Following the methodology proposed by Nakagawa and Schielzeth (2013), values accounted for the proportion of variance explained by the fixed factors of the model, while the values yielded the proportion of variance explained by the fixed and random model factors. All statistical analyses were run in the statistical program R 4.1.2 (R Core Team, 2021; https://www.r-project.org/).
Results
Three facets of plant diversity
Overall, 978 species belonging to 395 genera and 113 families were collected in this study. We found 177 species whose abundances reached 0.1% of the total abundance (Table S1) and 98 species whose abundances represented less than 0.001% of the total abundance (Table S2). In roadside plots, 689 plant species (70.45%) were found, in contrast to the 835 species (85.38%) found growing in interior communities. The 143 species (14.62%) occurring exclusively in roadside communities belonged to 14 families of angiosperms and 1 family of gymnosperms. Plants uniquely occurring interior plots included 288 species (29.44%), 70 genera and 40 families (10 families of pteridophytes, 29 families of angiosperms and 1 gymnosperm family). Comparisons of the different facets of plant diversity revealed that lower taxonomic and functional diversity index values were obtained for all q orders (q= [0, 1, 2]) compared to the phylogenetic diversity values (Table 1).
Table 1 Descriptive statistics of taxonomic, functional and phylogenetic Hill numbers for q = [0, 1, 2].
In addition, the three facets of plant diversity showed the same arrangement from relatively high to low values in the two categories of plant plots: interior plots > roadside plots (Figure S1). The taxonomic and functional diversity levels showed the strongest responses in dominant species, whereas phylogenetic diversity showed the strongest effects in rare species (Table 1). The Pearson correlation results derived for the overall plant diversity indices showed that the correlation coefficient between the taxonomic level (qTD), functional level (qFD) and phylogenetic level (qPD) declined with increasing q-values (Figure S2). Positive correlations were obtained between the taxonomic level (qTD) and functional level (qFD) (ρ > 0.963). The strongest correlations (ρ > 0.988) were found between the taxonomic level and functional level for q = 0 (0TD vs. 0FD) in the interior community, whereas the weakest relationship was found between the functional diversity index at q = 0 and the phylogenetic diversity index at q = 2 in the roadside plots (0FD vs. 2PD, ρ = 0.46).
NMDS analyses of CWM and CM
According to the NMDS, the functional identity patterns were dissimilar between roadside and interior plots (Figure 3). The stress values of the NMDS analyses were 0.118 and 0.132 (<0.2) for CM and CWM, respectively, indicating a good representation of data in the two-dimensional NMDS graphs. The NMDS diagram of CM showed that a clear difference in functional compositions between the roadside and interior plots (Figure 3A). Regarding the CWMs, roadside communities were dominated by species with fast leaf economics (such as relatively high SLA and low LDMC values), and interior communities were dominated by species with slow leaf economics, although there was no clear discrimination between roadside and interior plots in the NMDS diagram (Figure 3B).
Figure 3 NMDS of 76 plots based on their CWM (A) and CM (B) values.
Compositional differences between the roadside and interior vegetation communities
Both the local and regional mean dissimilarity values were lower in the roadside communities than in the interior communities for all diversity indices and q orders. The greatest differences were observed for the taxonomic and functional diversity levels at q=0 at the regional scale and at q= 2 at the local scale, whereas the lowest differences were found for phylogenetic dissimilarity at q=2 at both scales (Figure 4). In general, roads increased the community similarity for all facets of plant diversity in the roadside plots, and this effect was stronger for common and dominant species at both scales (Figure 4).
Figure 4 Boxplot comparing mean dissimilarity measures for interior and roadside plots for taxonomic [qTD], functional [qFD] and phylogenetic [qPD]) diversity for q = [0, 1, 2]. Error bars represent 95% confidence intervals. ***P < 0.001, *P < 0.05, n.s. not significant difference between roadside and interior communities.
Elevational patterns of plant diversity
Regarding our second research question, our findings revealed a differential effect of road disturbance on plant diversity along the elevation gradient. At least three “best” models were obtained for each facet of diversity (Table 2) through the model selection process. The model-averaged estimates of plant diversity and their conditional standard errors were listed in Table S4. The contributions of random effects to explaining the variabilities in all plant diversity indices were less than those of fixed effects, as determined by calculating the differences between and (Table 2).
Table 2 The entire set of plausible models (ΔAICc< 2) for plant diversity indices (taxonomic [qTD], functional [qFD], phylogenetic [qPD]) based on AICc model selection.
Most of the best models included the interaction between the plot type and the elevation and/or its quadratic term (Table 2), especially when modeling the phylogenetic diversity index. By comparing the plant diversity between the roadside and interior plots, in terms of the taxonomic and functional diversity indices, the negative effects of roads were greater at moderate and high elevations (>3500 m a.s.l) than at low elevations, while the effects of roads on phylogenetic diversity were greater at low and moderate elevations (<2800 m a.s.l) than at high elevations (Figure 5).
Figure 5 Model predictions (lines) with 95% confidence intervals for each q order (q = [0, 1, 2]) and plant diversity (taxonomic, functional and phylogenetic levels) along elevation gradient. Red solid lines represent predictions in interior plots; black doted lines are used for predictions in the roadside plots.
Discussion
We investigated the taxonomic, functional and phylogenetic responses of dominant, common and rare species to mountain roads, accounting for shifts in the magnitude of road effects along an elevation gradient. Our results show that dominant species had greater changes to road disturbances than common and rare species at the taxonomic and functional levels, whereas the opposite was true at the phylogenetic level. Furthermore, the greatest effect of roads was reflected in mid- and high-elevation areas at taxonomic and functional levels, while at the phylogenetic level, the greatest road impacts occurred in low- and mid-elevation areas. This study provides fundamental information regarding the effects of roads on plant diversity and may help to forecast future community compositions in response to new road construction in similar areas.
Roadside communities favor and exclude some species, functional traits and phylogenetic clades
The environmental heterogeneities induced by road construction may act at different spatial scales, thus modifying the compositions of mountain plant communities (Haider et al., 2018; Neher et al., 2013). We found that 14.62% of the total number of species grew exclusively in roadside communities; this finding confirmed that roads act as corridors for the dispersal of plant species (Lembrechts et al., 2017). Most species that appeared only in roadside plots were annual/perennial species that underwent early flowering, had small LDMCs, low SLA values, and anemo- and epizoochory dispersal modes, such as Cirsium arvense, Ajuga ciliate, Oxalis corniculata, Trifolium pratense, and Erigeron acris. These species are capable of rapid resource uptake, growth and tissue turnover, allowing them to adapt to the often-disturbed and nutrient-rich environments along roads (Balázs et al., 2016; Craven et al., 2018).
In addition, analysis of CWMs of all functional traits revealed no clear difference in functional composition between roadside and interior communities (Figure 3). This result is not surprising since abundance-based metrics attenuate the role of rare species in determining overall functional traits (Norden et al., 2017). At the same time, roads drove the absence of species that were present in interior plots. Most of the species that were absent in roadside plots were gymnosperms, such as Abies fargesii var. faxoniana, Picea brachytyla, Picea likiangensis var. rubescens, Pinus thunbergii, and ferns with relatively high SLA values, such as Adiantum pedatum and Athyrium costulalisorum. The gymnosperm and fern clades are older (Jin and Qian, 2019) than the Asteraceae species favored in the roadside communities (Alexander et al., 2016), suggesting that the changes in functional traits in this study can be attributed to the position of the species in a phylogeny (Munguía-Rosas et al., 2014).
Mountain roads reduce three facets of plant diversity, but by different magnitudes
A total of 32 nonnative species were found in this survey, including 4 invasive species (Li et al., 2022). Nonnative species that disperse via roads often have more similar functional traits (McDougall et al., 2018), so their introduction would lead to a certain decline in the functional diversity of roadside communities. However, considering the relatively small proportion of nonnative species richness in the study area, we believe that the high radiation, high nutrient and high disturbance environment of the roadside community (Johnston and Johnston, 2004; Delgado et al., 2007; De Pauw et al., 2021) led to the disappearance of many native species from the roadside plots, which was the main reason for the decrease of functional diversity. In other words, the communities formed following road construction have lower functional diversity levels than other communities because the environmental filters imposed by the construction of roads select species that are able to cope with these filters and excluded some species with unique functional traits (Laliberté et al., 2014; Li et al., 2018). Some of the older lineages present in the interior communities (such as ferns and gymnosperms) were absent from the roadside communities, thus resulting in the decreased phylogenetic diversity measured in the roadside plots (Arroyo-Rodríguez et al., 2012; Munguía-Rosas et al., 2014; De Pauw et al., 2021).
Furthermore, the effect of roads on the three facets of plant diversity can be explained not only by the decreased species richness but also by the changes that occur in species abundance. Although roads affect overall plant diversity, the multi-comparison of these plant diversity facets showed that road disturbances excluded rare species and reduced the abundance of common and dominant species (Figure 3). Thus, our results suggest that roadside plant communities may be more homogenization in the three facets of plant diversity than interior communities. This is consistent with previous findings regarding species richness in mountainous areas (Haider et al., 2018; Li et al., 2022). We detected those changes in taxonomic diversity would exceed changes in functional and phylogenetic diversity, signaling functional and phylogenetic resilience (Jarzyna and Jetz, 2018).
The use of phylogenetic and functional diversity versus taxonomic diversity provides different, complementary information on community assembly mechanisms in plant communities across the Qionglai mountain range. Phylogenetic niche conservatism holds that that closely related species tend to be ecologically similar (Losos, 2008; Wiens et al., 2010). However, phylogenetic niche conservatism was not common in the present study, suggesting that the plant diversity metrics used in the study were decoupled from one another (Grigoropoulou et al., 2022). The results of the correlation analysis also confirmed this result (Figure S2). The observed mismatches between phylogenetic and functional diversity in the study might have also been due to nonrandom filtering of species with unfit traits (Hacala et al., 2021; Monge-González et al., 2021). On the other hand, the taxonomic and functional plant diversity indices showed relatively strong positive correlations, suggesting that the community assemblages of mountain plants are ruled by functional differentiation (Freitas et al., 2021; Vega-Álvarez et al., 2019). We argue that roads in the Qionglai mountain range limit the fitness of some species and expand the range of distribution of other species due to the combination of direct disturbances and indirect microhabitat amelioration.
Mismatched effects of roads on taxonomic, functional and phylogenetic diversity along the elevation gradient
Here, we describe the differential shifts in taxonomic, functional and phylogenetic diversity along an elevation gradient due to roads in the Qionglai mountain range. These results are congruent with previous studies in which increasing taxonomic (Zhang et al., 2009; Li et al., 2022), functional (Zhang et al., 2014) and phylogenetic (Manish and Pandit, 2018; Zu et al., 2019) diversity levels were reported along elevation gradients and can be explained by moderate elevations being the center region for repeated plant migrations between low and high elevations in the alternating processes of glaciation and deglaciation that occurred along the elevation gradient during the Quaternary (Zhong et al., 1979; Zhang et al., 2009). Thus, plants had repeated opportunities to come into contact, hybridize, and become isolated, and for different populations of the same species to become physically separated and genetically divergent, as they moved into new, isolated niches (Liu et al., 2006; Zu et al., 2019).
To adapt to the complex hydrothermal conditions in mid-elevation regions, plants have also evolved diverse functional properties (Dıíaz and Cabido, 2001). However, the abundant species at higher elevations belonged to specific clades, and traits related to harsh environmental filtering at higher elevation evolved independently of phylogeny (Kitagawa et al., 2018). Species that are able to overcome harsh environments and are common/dominant at high elevations have closer evolutionary distances, thus resulting in the relatively low phylogenetic diversity levels of common and dominant species at high elevations (Webb et al., 2002; Qian, 2014; Kitagawa et al., 2018).
When simultaneously considering the effects of roads and elevation gradients, we revealed that the generally negative effects of roads on the three plant diversity facets were more pronounced at moderate and high elevations than at low elevations for taxonomic and functional diversity. In contrast, for phylogenetic diversity, the effects of mountain roads were smaller at high elevations, especially for common and dominant species. The impacts of roads were greater in places with greater diversity, likely because these places had more floral taxa, functional traits and phylogenetic clades that were not adapted to roadside disturbance (Zu et al., 2019). Therefore, planning new roads in a place with large numbers of plant species with unique traits and phylogenetic clades should be a last-resort decision.
Conclusions and implications for management and conservation
Our studies fill gaps in knowledge regarding the coordination of biodiversity changes in functional, taxonomic and phylogenetic facets of plant communities in response to mountain roads. The results show that mountain roads reduced plant diversity by reducing the abundance of individuals in the community and, to a lesser extent, reducing species richness. Furthermore, we found that the responses of taxonomic and functional diversity to mountain roads were consistent, but the responses of phylogenetic diversity to mountain roads were inconsistent with those of the two other facets of plant diversity. Overall, the multiple facets of plant diversity decrease due to the presence of mountain roads, and that such decreases are not necessarily synchronous or equal along the elevation gradient.
These findings highlight that the importance of considering the multiple facets of plant diversity to understand the effects of roads on plant communities to help management agencies develop effective and rational protection management measures for the national park, since attempts to maximize taxonomic or functional diversity may inadvertently lead to the erosion of phylogenetic diversity. Furthermore, the extent to which montane plant communities may be susceptible to roads depending on their current position along the elevation gradient. When road construction in this area is unavoidable, the road construction plan can be guided according to the facets of diversity with the greatest loss caused by roads, combined with the elevation of the plant community.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
HL: Conceptualization, Methodology, Investigation, Writing-original draft, Writing-review & editing. PL: Writing-review & editing, Funding acquisition. HY: Writing-review & editing. CL, WX, HJ, YC and YH: collected the trait data and conducted the fieldwork. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Demonstration of Monitoring and Protection of Important Species Habitat (grant number 2016YFC0503305); Key Technology and Demonstration for Biodiversity Conservation in Giant Panda National Park (grant number 2018SZDX0036); and the Biodiversity Survey and Assessment Project of the Ministry of Ecology and Environment, People’s Republic of China (grant number 2019HJ2096001006).
Acknowledgments
Thanks to Dr. Andre F. Clewell for his insightful and helpful comments and edits on this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.985673/full#supplementary-material
References
Alexander, J. M., Lembrechts, J. J., Cavieres, L. A., Daehler, C., Haider, S., Kueffer, C., et al. (2016). Plant invasions into mountains and alpine ecosystems: current status and future challenges. Alpine Bot. 126 (2), 89–103. doi: 10.1007/s00035-016-0172-8
Arévalo, J. R., Otto, R., Escudero, C., Fernández-Lugo, S., Arteaga, M., Delgado, J. D., et al. (2010). Do anthropogenic corridors homogenize plant communities at a local scale? a case studied in tenerife (Canary islands). Plant Ecol. 209 (1), 23–35. doi: 10.1007/s11258-009-9716-y
Arroyo-Rodríguez, V., Cavender-Bares, J., Escobar, F., Melo, F. P. L., Tabarelli, M., Santos, B. A. (2012). Maintenance of tree phylogenetic diversity in a highly fragmented rain forest. J. Ecol. 100 (3), 702–711. doi: 10.1111/j.1365-2745.2011.01952.x
Arroyo-Rodríguez, V., Rös, M., Escobar, F., Melo, F. P. L., Santos, B. A., Tabarelli, M., et al. (2013). Plant β-diversity in fragmented rain forests: testing floristic homogenization and differentiation hypotheses. J. Ecol. 101 (6), 1449–1458. doi: 10.1111/1365-2745.12153
Balázs, D., Bernadett, H., Tóthmérész, B. (2016). Grassland vegetation in urban habitats - testing ecological theories. TUEXENIA 36, 379–393. doi: 10.14471/2016.36.017
Barbosa, N. P. U., Wilson Fernandes, G., Carneiro, M. A. A., Júnior, L. A. C. (2010). Distribution of non-native invasive species and soil properties in proximity to paved roads and unpaved roads in a quartzitic mountainous grassland of southeastern Brazil (rupestrian fields). Biol. Invasions. 12 (11), 3745–3755. doi: 10.1007/s10530-010-9767-y
Bartoń, K. (2022) MuMIn: Multi-model inference. Available at: https://CRAN.R-project.org/package=MuMIn.
Berdugo, M., Maestre, F. T., Kéfi, S., Gross, N., Le Bagousse-Pinguet, Y., Soliveres, S. (2019). Aridity preferences alter the relative importance of abiotic and biotic drivers on plant species abundance in global drylands. J. Ecol. 107 (1), 190–202. doi: 10.1111/1365-2745.13006
Boufford, D. E. (2014). Biodiversity hotspot: China ’ s hengduan mountains Vol. 72 (Arnold Arboretum, Harvard University), 24–35. Available at: https://www.biodiversitylibrary.org/part/253573. Arnoldia. Jamaica Plain, Mass.
Brighenti, S., Hotaling, S., Finn, D. S., Fountain, A. G., Hayashi, M., Herbst, D., et al. (2021). Rock glaciers and related cold rocky landforms: Overlooked climate refugia for mountain biodiversity. Global Change Biol. 27 (8), 1504–1517. doi: 10.1111/gcb.15510
Burnham, K. P., Anderson, D. R. (2002). Model selection and multimodel inference: A practical information-theoretic approach (2nd ed.) (New York, NY: Springer-Verlag).
Chao, A., Chiu, C.-H., Hsieh, T. C. (2012). Proposing a resolution to debates on diversity partitioning. Ecology 93 (9), 2037–2051. doi: 10.1890/11-1817.1
Chao, A., Chiu, C.-H., Jost, L. (2014). Unifying species diversity, phylogenetic diversity, functional diversity, and related similarity and differentiation measures through hill numbers. Annu. Rev. Ecol. Evol. Syst. 45 (1), 297–324. doi: 10.1146/annurev-ecolsys-120213-091540
Chiu, C.-H., Jost, L., Chao, A. (2014). Phylogenetic beta diversity, similarity, and differentiation measures based on hill numbers. Ecol. Monographs. 84 (1), 21–44. doi: 10.1890/12-0960.1
Craven, D., Eisenhauer, N., Pearse, W. D., Hautier, Y., Isbell, F., Roscher, C., et al. (2018). Multiple facets of biodiversity drive the diversity–stability relationship. Nat. Ecol. Evol. 2 (10), 1579–1587. doi: 10.1038/s41559-018-0647-7
Dar, P., Reshi, Z., Shah, A. (2018). Altitudinal distribution of native and alien plant species along roadsides in Kashmir himalaya, India. Trop. Ecol. 59, 45–55. doi: 10.1111/j.1442-9993.2010.02134.x.
Delgado, J. D., Arroyo, N. L., Arévalo, J. R., Fernández-Palacios, J. M. (2007). Edge effects of roads on temperature, light, canopy cover, and canopy height in laurel and pine forests (Tenerife, canary islands). Landscape Urban Planning. 81 (4), 328–340. doi: 10.1016/j.landurbplan.2007.01.005
Deljouei, A., Sadeghi, S. M. M., Abdi, E., Bernhardt-Römerman, M., Pascoe, E. L., Marcantonio, M. (2018). The impact of road disturbance on vegetation and soil properties in a beech stand, hyrcanian forest. Eur. J. For. Res. 137 (6), 759–770. doi: 10.1007/s10342-018-1138-8
De Pauw, K., Meeussen, C., Govaert, S., Sanczuk, P., Vanneste, T., Bernhardt-Römermann, M., et al. (2021). Taxonomic, phylogenetic and functional diversity of understorey plants respond differently to environmental conditions in European forest edges. J. Ecol. 109 (7), 2629–2648. doi: 10.1111/1365-2745.13671
Dıíaz, S., Cabido, M. (2001). Vive la différence: plant functional diversity matters to ecosystem processes. Trends Ecol. Evol. 16 (11), 646–655. doi: 10.1016/S0169-5347(01)02283-2
Freitas, C., Brum, F., Silva, C., Maracahipes, L., Carlucci, M., Collevatti, R., Bacon, C. (2021). Incongruent Spatial Distribution of Taxonomic, Phylogenetic, and Functional Diversity in Neotropical Cocosoid Palms. Front. for Glob. Change 4. doi: 10.3389/ffgc.2021.739468.
Gower, J. C. (1971). A general coefficient of similarity and some of its properties. Biometrics 27, 857. doi: 10.2307/2528823
Grigoropoulou, A., Schmidt-Kloiber, A., Múrria, C. (2022). Incongruent latitudinal patterns of taxonomic, phylogenetic and functional diversity reveal different drivers of caddisfly community assembly across spatial scales. Global Ecol. Biogeogr. 31 (5), 1006–1020. doi: 10.1111/geb.13479
Hacala, A., Lafage, D., Prinzing, A., Sawtschuk, J., Pétillon, J. (2021). Drivers of taxonomic, functional and phylogenetic diversities in dominant ground-dwelling arthropods of coastal heathlands. Oecologia 197 (2), 511–522. doi: 10.1007/s00442-021-05032-4
Haddad, N. M., Brudvig, L. A., Clobert, J., Davies, K. F., Gonzalez, A., Holt, R. D., et al. (2015). Habitat fragmentation and its lasting impact on earth's ecosystems. Sci. Advances. 1 (2), e1500052. doi: 10.1126/sciadv.1500052
Haider, S., Kueffer, C., Bruelheide, H., Seipel, T., Alexander, J. M., Rew, L. J., et al. (2018). Mountain roads and non-native species modify elevational patterns of plant diversity. Global Ecol. Biogeogr. 27 (6), 667–678. doi: 10.1111/geb.12727
Hansen, M. J., Clevenger, A. P. (2005). The influence of disturbance and habitat on the presence of non-native plant species along transport corridors. Biol. Conserv. 125 (2), 249–259. doi: 10.1016/j.biocon.2005.03.024
Hillebrand, H., Blasius, B., Borer, E. T., Chase, J. M., Downing, J. A., Eriksson, B. K., et al. (2018). Biodiversity change is uncoupled from species richness trends: Consequences for conservation and monitoring. J. Appl. Ecol. 55 (1), 169–184. doi: 10.1111/1365-2664.12959
Hill, S. L. L., Harfoot, M., Purvis, A., Purves, D. W., Collen, B., Newbold, T., et al. (2016). Reconciling biodiversity indicators to guide understanding and action. Conserv. Letters. 9 (6), 405–412. doi: 10.1111/conl.12291
Jarzyna, M.A., Jetz, W (2018). Taxonomic and functional diversity change is scale dependent. Nat. Commun. 9, 2565. doi: 10.1038/s41467-018-04889-z.
Jin, Y., Qian, H. (2019). V.PhyloMaker: an r package that can generate very large phylogenies for vascular plants. Ecography 42 (8), 1353–1359. doi: 10.1111/ecog.04434
Johnston, F. M., Johnston, S. W. (2004). Impacts of road disturbance on soil properties and on exotic plant occurrence in subalpine areas of the Australian Alps. Arctic Antarctic Alpine Res. 36 (2), 201–207. doi: 10.1657/1523-0430(2004)036[0201:IORDOS]2.0.CO;2
Jost, L. (2006). Entropy and diversity. Oikos 113 (2), 363–375. doi: 10.1111/j.2006.0030-1299.14714.x
Jost, L. (2007). Partitioning diversity into independent alpha and beta components. Ecology 88 (10), 2427–2439. doi: 10.1890/06-1736.1
Jost, L. (2010). The relation between evenness and diversity. Diversity 2 (2), 207–232. doi: 10.3390/d2020207
Kang, D. (2021). A review of the impacts of four identified major human disturbances on the habitat and habitat use of wild giant pandas from 2015 to 2020. Sci. Total Environ. 763, 142975. doi: 10.1016/j.scitotenv.2020.142975
Karmakar, N., Ballav, J., Hazra, A. (2021). Roadside vegetation in and around the megacity of kolkata along an urbanisation gradient. Acta Botanica Hungarica. 63, 77–86. doi: 10.1556/034.63.2021.1-2.6
Kattge, J., Bönisch, G., Díaz, S., Lavorel, S., Prentice, C., Leadley, P., et al. (2020). TRY plant trait database-enhanced coverage and open access. Glob. Change Biol. 26, 119–188. doi: 10.1111/gcb.14904
Kitagawa, R., Koide, D., Mori, A. S. (2018). Different trends in phylogenetic and functional structure of plant communities along an elevation gradient. Ecol. Res. 33 (6), 1233–1243. doi: 10.1007/s11284-018-1638-5
Laliberté, E., Legendre, P. (2010). A distance-based framework for measuring functional diversity from multiple traits. Ecology 91 (1), 299–305. doi: 10.1890/08-2244.1
Laliberté, E., Zemunik, G., Turner, B. L. (2014). Environmental filtering explains variation in plant diversity along resource gradients. Science 345 (6204), 1602–1605. doi: 10.1126/science.1256330
Laurance, W. F., Goosem, M., Laurance, S. G. W. (2009). Impacts of roads and linear clearings on tropical forests. Trends Ecol. evol. 24 (12), 659–669. doi: 10.1016/j.tree.2009.06.009
Lembrechts, J. J., Alexander, J. M., Cavieres, L. A., Haider, S., Lenoir, J., Kueffer, C., et al. (2017). Mountain roads shift native and non-native plant species' ranges. Ecography 40 (3), 353–364. doi: 10.1111/ecog.02200
Li, D. (2018). hillR: taxonomic, functional, and phylogenetic diversity and similarity through hill numbers. J. Open Source Software 3, 1041. doi: 10.21105/joss.01041
Li, S., Dong, S., Zhang, X., Liu, S., Shi, J., Gao, X., et al. (2018). Evolutionary history and functional traits determine the spatial pattern of multifaceted plant diversity in a typical temperate desert disturbed by an expressway. Sci. Total Environ. 635, 972–983. doi: 10.1016/j.scitotenv.2018.04.137
Li, H., Luo, P., Yang, H., Li, T., Luo, C., Wu, S., et al. (2022). Effect of road corridors on plant diversity in the qionglai mountain range, China. Ecol. Indicators. 134, 108504. doi: 10.1016/j.ecolind.2021.108504
Li, X., Tian, H., Wang, Y., Li, R., Song, Z., Zhang, F., et al. (2013). Vulnerability of 208 endemic or endangered species in China to the effects of climate change. Regional Environ. Change. 13 (4), 843–852. doi: 10.1007/s10113-012-0344-z
Liu, J.-Q., Wang, Y.-J., Wang, A.-L., Hideaki, O., Abbott, R. J. (2006). Radiation and diversification within the ligularia–Cremanthodium–Parasenecio complex (Asteraceae) triggered by uplift of the qinghai-Tibetan plateau. Mol. Phylogenet. Evol. 38 (1), 31–49. doi: 10.1016/j.ympev.2005.09.010
Losos, J. B. (2008). Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species. Ecol. Letters. 11 (10), 995–1003. doi: 10.1111/j.1461-0248.2008.01229.x
Manish, K., Pandit, M. K. (2018). Phylogenetic diversity, structure and diversification patterns of endemic plants along the elevational gradient in the Eastern himalaya. Plant Ecol. Divers. 11 (4), 501–513. doi: 10.1080/17550874.2018.1534147
McDougall, K. L., Lembrechts, J., Rew, L. J., Haider, S., Cavieres, L. A., Kueffer, C., et al. (2018). Running off the road: Roadside non-native plants invading mountain vegetation. Biol. Invasions 20, 3461–3473. doi: 10.1007/s10530018-1787-z
Millennium Ecosystem Assessment (2005). Ecosystems and human well-being: Biodiversity synthesis (Washington, DC: World Resources Institute).
Miller, S. C., Wiethase, J. H., Motove Etingue, A., Franklin, E., Fero, M., Wolfe, J. D., et al. (2021). Interactive effects of elevation and newly paved road on avian community composition in a scientific reserve, bioko island, equatorial Guinea. Biotropica 53 (6), 1646–1663. doi: 10.1111/btp.13014
Monge-González, M. L., Guerrero-Ramírez, N., Krömer, T., Kreft, H., Craven, D. (2021). Functional diversity and redundancy of tropical forests shift with elevation and forest-use intensity. J. Appl. Ecol. 58 (9), 1827–1837. doi: 10.1111/1365-2664.13955
Munguía-Rosas, M. A., Jurado-Dzib, S. G., Mezeta-Cob, C. R., Montiel, S., Rojas, A., Pech-Canché, J. M. (2014). Continuous forest has greater taxonomic, functional and phylogenetic plant diversity than an adjacent naturally fragmented forest. J. Trop. Ecol. 30 (4), 323–333. doi: 10.1017/S0266467414000194
Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., Kent, J. (2000). Biodiversity hotspots for conservation priorities. Nature 403 (6772), 853–858. doi: 10.1038/35002501
Nakagawa, S., Schielzeth, H. (2013). A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142. doi: 10.1111/j.2041-210x.2012.00261.x
Neher, D. A., Asmussen, D, Lovell, S. T. (2013). Roads in northern hardwood forests affect adjacent plant communities and soil chemistry in proportion to the maintained roadside area. Sci.Total Environ. 449, 320–327. doi: 10.1016/j.scitotenv.2013.01.062
Norden, N., Boukili, V., Chao, A., Ma, K. H., Letcher, S. G., Chazdon, R. L. (2017). Opposing mechanisms affect taxonomic convergence between tree assemblages during tropical forest succession. Ecol. Letters. 20 (11), 1448–1458. doi: 10.1111/ele.12852
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020) Vegan: Community ecology package. Available at: https://CRAN.R-project.org/package=vegan.
Pauchard, A., Alaback, P. B. (2004). Influence of elevation, land use, and landscape context on patterns of alien plant invasions along roadsides in protected areas of south-central Chile. Conserv. Biol. 18 (1), 238–248. doi: 10.1111/j.1523-1739.2004.00300.x
Pérez-Harguindeguy, N., Díaz, S., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., et al. (2013). New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot. 61 (3), 167–234. doi: 10.1071/BT12225
Perrigo, A., Hoorn, C., Antonelli, A. (2020). Why mountains matter for biodiversity. J. Biogeogr. 47 (2), 315–325. doi: 10.1111/jbi.13731
Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D., R Core Team (2021) Nlme: Linear and nonlinear mixed effects models. Available at: https://cran.r-project.org/web/packages/nlme/index.html.
Qian, H. (2014). Contrasting relationships between clade age and temperature along latitudinal versus elevational gradients for woody angiosperms in forests of south America. J. Vegetation Science. 25 (5), 1208–1215. doi: 10.1111/jvs.12175
Qiu, C., Shen, Z., Peng, P., Mao, L., Zhang, X. (2014). How does contemporary climate versus climate change velocity affect endemic plant species richness in China? Chin. Sci. Bull. 59, 34, 4660–4667. doi: 10.1007/s11434-014-0640-8
Rashid, I., Haq, S. M., Lembrechts, J. J., Khuroo, A. A., Pauchard, A., Dukes, J. S. (2021). Railways redistribute plant species in mountain landscapes. J. Appl. Ecol. 58 (9), 1967–1980. doi: 10.1111/1365-2664.13961
Razafindratsima, O. H., Brown, K. A., Carvalho, F., Johnson, S. E., Wright, P. C., Dunham, A. E. (2018). Edge effects on components of diversity and above-ground biomass in a tropical rainforest. J. Appl. Ecol. 55 (2), 977–985. doi: 10.1111/1365-2664.12985
R Core Team (2021). R: A language and environment for statistical computing (Vienna: R Foundation for Statistical Computing). Available at: https://www.R-project.org/.
Salick, J., Fang, Z., Hart, R. (2019). Rapid changes in eastern Himalayan alpine flora with climate change. Am. J. Bot. 106 (4), 520–530. doi: 10.1002/ajb2.1263
Seipel, T., Kueffer, C., Rew, L. J., Daehler, C. C., Pauchard, A., Naylor, B. J., et al. (2012). Processes at multiple scales affect richness and similarity of non-native plant species in mountains around the world. Global Ecol. Biogeogr. 21 (2), 236–246. doi: 10.1111/j.1466-8238.2011.00664.x
Seto, K. C., Güneralp, B., Hutyra, L. R. (2012). Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proc. Natl. Acad. Sci. 109 (40), 16083–16088. doi: 10.1073/pnas.1211658109
Shi, Z., Xu, X., Souza, L., Wilcox, K., Jiang, L., Liang, J., et al. (2016). Dual mechanisms regulate ecosystem stability under decade-long warming and hay harvest. Nat. Commun. 7 (1), 11973. doi: 10.1038/ncomms11973
Smith, S. A., Brown, J. W. (2018). Constructing a broadly inclusive seed plant phylogeny. Am. J. Bot. 105 (3), 302–314. doi: 10.1002/ajb2.1019
State Council Information Office of China (SCIO) (2015). Available at: http://www.scio.gov.cn/xwfbh/gbwxwfbh/fbh/Document/1395514/1395514.htm.
State Forestry Administration (SFA) (2015). The 4rd national survey report on giant panda in China (Beijing: Science Press).
Tang, Z., Wang, Z., Zheng, C., Fang, J. (2006). Biodiversity in china's mountains. Front. Ecol. Environ. 4 (7), 347–352. doi: 10.1890/1540-9295(2006)004[0347:BICM]2.0.CO;2
Tuomisto, H. (2010). A diversity of beta diversities: straightening up a concept gone awry. part 1. defining beta diversity as a function of alpha and gamma diversity. Ecography 33 (1), 2–22. doi: 10.1111/j.1600-0587.2009.05880.x
Vega-Álvarez, J., García-Rodríguez, J. A., Cayuela, L. (2019). Facilitation beyond species richness. J. Ecol. 107 (2), 722–734. doi: 10.1111/1365-2745.13072
Webb, C. O., Ackerly, D. D., McPeek, M. A., Donoghue, M. J. (2002). Phylogenies and community ecology. Annu. Rev. Ecol. Syst. 33, 475–505. doi: 10.1146/annurev.ecolsys.33.010802.150448
Wiens, J. J., Ackerly, D. D., Allen, A. P., Anacker, B. L., Buckley, L. B., Cornell, H. V., et al. (2010). Niche conservatism as an emerging principle in ecology and conservation biology. Ecol. Letters. 13 (10), 1310–1324. doi: 10.1111/j.1461-0248.2010.01515.x
Xie, S., Wang, X., Yang, T., Huang, B., Wang, W., Lu, F., et al. (2021). Effect of highways on breeding birds: Example of hulunbeier, China. Global Ecol. Conserv. 27, e01554. doi: 10.1016/j.gecco.2021.e01554
Xu, W., Ouyang, Z., Viña, A., Zheng, H., Liu, J., Xiao, Y. (2006). Designing a conservation plan for protecting the habitat for giant pandas in the qionglai mountain range, China. Diversity Distrib. 12 (5), 610–619. doi: 10.1111/j.1366-9516.2006.00236.x
Yin, L., Dai, E., Zheng, D., Wang, Y., Ma, L., Tong, M. (2020). What drives the vegetation dynamics in the hengduan mountain region, southwest China: Climate change or human activity? Ecol. Indicators. 112, 106013. doi: 10.1016/j.ecolind.2019.106013
Zanne, A. E., Tank, D. C., Cornwell, W. K., Eastman, J. M., Smith, S. A., FitzJohn, R. G., et al. (2014). Three keys to the radiation of angiosperms into freezing environments. Nature 506 (7486), 89–92. doi: 10.1038/nature12872
Zhang, J.-T., Li, M., Nie, E. (2014). Pattern of functional diversity along an altitudinal gradient in the baihua mountain reserve of Beijing, China. Braz. J. Bot. 37 (1), 37–45. doi: 10.1007/s40415-014-0047-4
Zhang, D.-C., Zhang, Y.-H., Boufford, D. E., Sun, H. (2009). Elevational patterns of species richness and endemism for some important taxa in the hengduan mountains, southwestern China. Biodivers. Conserv. 18 (3), 699–716. doi: 10.1007/s10531-008-9534-x
Zhong, Z. C., Qin, W. C., Xu, M. Q. (1979). Discussion on the historical evolution of vegetation geography in sichuan. J. Southwest Normal University: Natural Sci. Edition 1, 4–16. doi: 10.13718/j.cnki.xsxb.1979.01.001
Keywords: hill numbers, roads, mountain plant assemblages, functional traits, phylogenetic diversity, biodiversity hotspots
Citation: Li H, Luo P, Yang H, Luo C, Xie W, Jia H, Cheng Y and Huang Y (2022) Assessing the effect of roads on mountain plant diversity beyond species richness. Front. Plant Sci. 13:985673. doi: 10.3389/fpls.2022.985673
Received: 04 July 2022; Accepted: 12 September 2022;
Published: 26 September 2022.
Edited by:
Mani Shrestha, University of Bayreuth, GermanyCopyright © 2022 Li, Luo, Yang, Luo, Xie, Jia, Cheng and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Luo, luopeng@cib.ac.cn; Hao Yang, yanghao@cib.ac.cn