Alpine cushion plants inhibit the loss of phylogenetic diversity in severe environments

Positive interactions and interdependence in communities

Facilitative interactions play crucial roles in community organization, and the stress gradient hypothesis (SGH) provides a simple conceptual framework for the context-dependency of competitive and facilitative interactions. The idea is that positive interactions are more common under high physical and consumer stress, where species benefit from stress-tolerant neighbors, than in benign environments. We explore insights from the SGH into ecological generality, niche theory, community assembly, and diversity effects on ecosystem function and discuss how the SGH can inform our understanding of rapid evolution, mutualisms, exotic invasions, and facilitation cascades. We suggest that, with escalating global stresses, the SGH may provide a conceptual template for an interdependent perspective in ecology that can contribute to conservation and restoration efforts.

Cushions serve as conservation refuges for the Himalayan alpine plant diversity: Implications for nature-based environmental management

Globally, rapid climate and land-use changes in alpine environments are posing severe risks to their bountiful biodiversity and ecosystem services. Currently, nature-based solutions are fast-emerging as the preferred approach to address the challenges of environmental sustainability. In alpine environments, cushion plants owing to their unique architecture and adaptability offer a potential nature-based system to plan biodiversity conservation and habitat restoration strategies. Here, we employed an analytical framework to test whether and how the cushion plants facilitate the sustenance of alpine plant diversity in Kashmir Himalaya. We specifically aimed to answer: what are the effects of the cushion plants on the patterns of alpine species richness and phylogenetic diversity, and whether these effects vary across spatial scales (local versus landscape), cushion types, and changing elevation. We randomly selected pairs of cushion and neighbouring non-cushion plots (size 100 m2) across 34 different alpine sites in the study region. Within each plot, we randomly laid three 5 m2 quadrats for vegetation sampling, and sampled a total of 204 quadrats in 68 plots with seven cushion types along elevation ranging from 3100 to 3850 m. Our results revealed positive effects of the cushions by supporting a higher community species richness (SR) and phylogenetic diversity (PD). The effects were consistent both at the local (i.e., quadrat) and landscape (i.e., plot) scales, but varied significantly with the cushion type. Interestingly, SR and PD showed an increasing trend with increase in elevation in cushion communities, thereby supporting stress gradient hypothesis. Along the elevational gradient, the cushion communities showed phylogenetic overdispersion, but clustering by non-cushions. Overall, our study provides empirical evidence to reinforce the role of the cushions as conservation refugia for an imperilled alpine plant diversity in the Himalaya. Looking ahead, we highlight the far-reaching implications of our findings in guiding the nature-based environmental management of alpine ecosystems worldwide.

Degeneration of foundation cushion species induced by ecological constraints can cause massive changes in alpine plant communities

Foundational cushion plants can re-organize community structures and sustain a prominent proportion of alpine biodiversity, but they are sensitive to climate change. The loss of cushion species can have broad consequences for associated biota. The potential plant community changes with the population dynamics of cushion plants remain, however, unclear. Using eight plant communities along a climatic and community successional gradient, we assessed cushion population dynamics, the underlying ecological constraints and hence associated plant community changes in alpine communities dominated by the foundational cushion plant Arenaria polytrichoides. The population dynamics of Arenaria are attributed to ecological constraints at a series of life history stages. Reproductive functions are constrained by increasing associated beneficiary plants; subsequent seedling establishment is constrained by temperature, water and light availability, extreme climate events, and interspecific competition; strong competitive exclusion may accelerate mortality and degeneration of cushion populations. Along with cushion dynamics, species composition, abundance and community structure gradually change. Once cushion plants completely degenerate, previously cushion-dominated communities shift to relatively stable communities that are overwhelmingly dominated by sedges. Climate warming may accelerate the degeneration process of A. polytrichoides. Degeneration of this foundational cushion plant will possibly induce massive changes in alpine plant communities and hence ecosystem functions in alpine ecosystems. The assessment of the population dynamics of foundation species is critical for an effective conservation of alpine biodiversity.

Plant-plant interactions influence post-fire recovery depending on fire history and nurse growth form

Backgrounds

Plant-plant interactions are among the most important factors affecting the natural recovery of vegetation. While the impacts of nurse plants on species composition and biodiversity are well documented, the effects of different nurse’s growth forms on all biodiversity components including taxonomic, functional, and phylogenetic diversity have been less studied and compared, especially for their effects on different times after fire disturbance. This research was focused on comparing the effects of a perennial grass (Elymus hispidens), a perennial herb (Phlomis cancellata), and a high shrub species (Lonicera nummulariifolia) on species composition and the biodiversity components, and how these impacts change across five sites with short-term (1 and 4 years sites), long-term (10 and 20 years sites) times since last fire and a control site where no fire was known in recorded history in semi-arid shrublands of Fereizi Chenaran located in Northeast of Iran. The changes of species composition and taxonomic, functional, and phylogenetic diversity were calculated with respect to the presence/absence of nurse’s growth forms, fire history, and their interactions.

Results

Nurse shrubs affected species composition and all biodiversity components, whereas all indices were reduced when considering Elymus grass as nurse plant. On the other hand, the herb Phlomis enhanced species composition and taxonomic diversity, while it had a negative effect on functional and phylogenetic diversity. Such specific effects of nurse types were mostly observed under long timescales (i.e., 10- and 20-year sites). Interestingly, the relative importance of nurse types and time since the last fire largely explained the variation of species composition and biodiversity components, with larger effects of nurse types on all biodiversity components. However, we found a significant contribution of fire explaining variation of species composition and phylogenetic diversity.

Conclusions

These results indicated nurse plants can affect the post-fire recovery of vegetation by providing specific mechanisms controlling beneficiary relatedness depending on their growth forms and time scales since the last fire. Therefore, these findings suggest perennial plants in the form of nurse species as a useful factor to develop techniques of active restoration in burned ecosystems.

Plant-plant interactions determine natural restoration of plant biodiversity over time, in a degraded mined land

Restoration of degraded environments is essential to mitigate adverse impacts of human activities on ecosystems. Plant-plant interactions may provide effective means for restoring degraded arid lands, but little is understood about these impacts. In this regard, we analyzed the effects of two dominant nurse plants (i.e., Artemisia sieberi and Stipa arabica) on taxonomic, functional, and phylogenetic diversity across different ages of land abandonment (i.e., control, recent, and old ages) in a limestone mine site in Iran. In addition, we considered two spatial scales: i) the plot scale (i.e., under 1m2 plots) and ii) the vegetation-patch scale (i.e., under the canopies of nurse plants), to assess nurse plant effects, land abandonment ages, and their relative importance on biodiversity facets by performing Kruskal-Wallis H test and variation partitioning analysis. Our results indicated an increase in taxonomic, functional, and phylogenetic diversity at the plot scale, when considering the presence of nurse plants under old ages of land abandonment. Such significant differences were consistent with the positive effects of Artemisia patches on taxonomic diversity and Stipa patches on functional and phylogenetic diversity. In addition, we found a larger contribution from nurse plants than land abandonment age on biodiversity variation at both spatial scales studied. Therefore, these results indicate the importance of plant-plant interactions in restoring vegetation, with their effects on the presence of beneficiary species and their functional and phylogenetic relatedness depending on the nurse life forms under the stress-gradient hypothesis.

Diversity patterns of cushion plants on the Qinghai-Tibet Plateau: A basic study for future conservation efforts on alpine ecosystems

The Qinghai-Tibet Plateau (QTP) is an important cushion plant hotspot. However, the distribution of cushion plants on the QTP is unknown, as are the factors that drive cushion plant distribution, limiting our understanding of the evolution of cushion species in the region. In this study, we assessed spatial patterns of total cushion plant diversity (including taxonomic and phylogenetic) over the entire QTP and compared patterns of diversity of cushion plants with different typologies (i.e., compact vs. loose). We also examined how these patterns were related to climatic features. Our results indicate that the southern QTP hosts the highest total cushion plant richness, especially in the south-central Hengduan Mountains subregion. The total number of cushion species declines from south to north and from southeast to northwest. Compact cushion plants exhibit similar patterns as the total cushion plant richness, whereas loose cushion plants show random distribution. Cushion plant phylogenetic diversity showed a similar pattern as that of the total cushion plant richness. In addition, cushion plant phylogenetic community structure was clustered in the eastern and southwestern QTP, whereas random or overdispersed in other areas. Climatic features represented by annual energy and water trends, seasonality and extreme environmental factors, had significant effects on cushion plant diversity patterns but limited effects on the phylogenetic community structure, suggesting that climatic features indeed promote the formation of cushion plants. Because cushion plants play vital roles in alpine ecosystems, our findings not only promote our understanding of the evolution and formation of alpine cushion plant diversity but also provide an indispensable foundation for future studies on cushion plant functions and thus alpine ecosystem sustainability in the entire QTP region.

Discriminating ecological processes affecting different dimensions of α- and β-diversity in desert plant communities

Understanding the spatial distribution of plant diversity and its drivers are major challenges in biogeography and conservation biology. Integrating multiple facets of biodiversity (e.g., taxonomic, phylogenetic, and functional biodiversity) may advance our understanding on how community assembly processes drive the distribution of biodiversity. In this study, plant communities in 60 sampling plots in desert ecosystems were investigated. The effects of local environment and spatial factors on the species, functional, and phylogenetic α- and β-diversity (including turnover and nestedness components) of desert plant communities were investigated. The results showed that functional and phylogenetic α-diversity were negatively correlated with species richness, and were significantly positively correlated with each other. Environmental filtering mainly influenced species richness and Rao quadratic entropy; phylogenetic α-diversity was mainly influenced by dispersal limitation. Species and phylogenetic β-diversity were mainly consisted of turnover component. The functional β-diversity and its turnover component were mainly influenced by environmental factors, while dispersal limitation dominantly effected species and phylogenetic β-diversity and their turnover component of species and phylogenetic β-diversity. Soil organic carbon and soil pH significantly influenced different dimensions of α-diversity, and soil moisture, salinity, organic carbon, and total nitrogen significantly influenced different dimensions of α- and β-diversity and their components. Overall, it appeared that the relative influence of environmental and spatial factors on taxonomic, functional, and phylogenetic diversity differed at the α and β scales. Quantifying α- and β-diversity at different biodiversity dimensions can help researchers to more accurately assess patterns of diversity and community assembly.

Diverse phylogenetic neighborhoods enhance community resistance to drought in experimental assemblages

Although the role played by phylogeny in the assembly of plant communities remains as a priority to complete the theory of species coexistence, experimental evidence is lacking. It is still unclear to what extent phylogenetic diversity is a driver or a consequence of species assembly processes. We experimentally explored how phylogenetic diversity can drive the community level responses to drought conditions in annual plant communities. We manipulated the initial phylogenetic diversity of the assemblages and the water availability in a common garden experiment with two irrigation treatments: average natural rainfall and drought, formed with annual plant species of gypsum ecosystems of Central Spain. We recorded plant survival and the numbers of flowering and fruiting plants per species in each assemblage. GLMMs were performed for the proportion of surviving, flowering, fruiting plants per species and for total proportion of surviving species and plants per pot. In water limited conditions, high phylogenetic diversity favored species coexistence over time with higher plant survival and more flowering and fruiting plants per species and more species and plants surviving per pot. Our results agree with the existence of niche complementarity and the convergence of water economy strategies as major mechanisms for promoting species coexistence in plant assemblages in semiarid Mediterranean habitats. Our findings point to high phylogenetic diversity among neighboring plants as a plausible feature underpinning the coexistence of species, because the success of each species in terms of surviving and producing offspring in drought conditions was greater when the initial phylogenetic diversity was higher. Our study is a step forward to understand how phylogenetic relatedness is connected to the mechanisms determining the maintenance of biodiversity.

Pollination success increases with plant diversity in high-Andean communities

Pollinator-mediated plant–plant interactions have traditionally been viewed within the competition paradigm. However, facilitation via pollinator sharing might be the rule rather than the exception in harsh environments. Moreover, plant diversity could be playing a key role in fostering pollinator-mediated facilitation. Yet, the facilitative effect of plant diversity on pollination remains poorly understood, especially under natural conditions. By examining a total of 9371 stigmas of 88 species from nine high-Andean communities in NW Patagonia, we explored the prevalent sign of the relation between conspecific pollen receipt and heterospecific pollen diversity, and assessed whether the incidence of different outcomes varies with altitude and whether pollen receipt relates to plant diversity. Conspecific pollen receipt increased with heterospecific pollen diversity on stigmas. In all communities, species showed either positive or neutral but never negative relations between the number of heterospecific pollen donor species and conspecific pollen receipt. The incidence of species showing positive relations increased with altitude. Finally, stigmas collected from communities with more co-flowering species had richer heterospecific pollen loads and higher abundance of conspecific pollen grains. Our findings suggest that plant diversity enhances pollination success in high-Andean plant communities. This study emphasizes the importance of plant diversity in fostering indirect plant–plant facilitative interactions in alpine environments, which could promote species coexistence and biodiversity maintenance.

Phylogenetic distance and resource availability mediate direction and strength of plant interactions in a competition experiment

Phylogenetic ecology uses evolutionary history to improve understanding of plant interactions. Phylogenetic distance can mediate plant interactions such as competition (e.g., via limiting similarity) and facilitation (e.g., via niche complementarity), influencing community assembly patterns. Previous research has found evidence both for and against a relationship between phylogenetic distance and the strength of plant interactions, and has found that other factors, such as trait differences, may be more influential. In addition to phylogenetic distance and species' traits, environmental conditions can also influence competition, with facilitative interactions-particularly among distantly related species-potentially becoming more pronounced under stressful, resource-limited conditions. We tested the prediction that greater phylogenetic distance is associated with decreased competition in a greenhouse experiment using plant species of the North American tallgrass prairie. We calculated the Relative Interaction Index for 81 species pairs using plant height, leaf length, and biomass as indicators of performance. We found that phylogenetic distance alone did not significantly affect competition. However, the interaction between phylogenetic distance and stressful conditions (sandier soils with low nutrient availability and water retention vs. resource-rich potting soil) altered plant traits and competition. Under stressful conditions, more distantly related species competed more strongly, leading to smaller plants. Conversely, under benign conditions more distantly related species pairs competed less and were larger. These results were contrary to our expectations that distant relatives would compete less under stressful conditions. Our experiment provides evidence that, while relatedness alone may not drive competition, phylogenetic distance can nonetheless be influential through interactions with environmental conditions.

Evolutionary assembly of flowering plants into sky islands

Alpine floras (plants in the vegetation belts above the climatic treelines) experience cold climates, particularly in temperate mountains during winter, but they are generally species-rich. Yet, biogeographers have not determined whether these floras represent evolutionarily independent but convergent assemblages drawn from their regional floras, or whether they originated from particular clades pre-adapted to harsh conditions. Here, we analyse the evolutionary relationships of angiosperm (flowering plant) species in 63 alpine floras worldwide (~7,000 species) in comparison with their regional floras (~94,000 species) and with the entire global flora. We find that each of the alpine floras represents an assemblage of more closely related species than their respective regional floras. The degree of phylogenetic clustering of species in alpine floras in tropical mountains exceeds that in temperate mountains. However, in relation to the global flora, temperate alpine floras are phylogenetically closely related subsets of floras that colonized cold temperate areas during interglacial periods. We conclude that alpine floras include a few dominant families that have evolved tolerance to low temperature, and that evolutionary niche conservatism explains their phylogenetic clustering, compared with species in their regional species pools.

Phylogenetic structure of understorey annual and perennial plant species reveals opposing responses to aridity in a Mediterranean biodiversity hotspot

Aridity is a critical driver of the diversity and composition of plant communities. However, how aridity influences the phylogenetic structure of functional groups (i.e. annual and perennial species) is far less understood than its effects on species richness. As perennials have to endure stressful conditions during the summer drought, as opposed to annuals that avoid it, they may be subjected to stronger environmental filtering. In contrast, annuals may be more susceptible to interannual climatic variability. Here we studied the phylogenetic structure of the annual and perennial components of understorey plant communities, along a regional aridity gradient in Mediterranean drylands. Specifically, we asked: (1) How do species richness (S) and phylogenetic structure (PS) of annuals and perennials in plant communities respond to aridity? (2) What is the contribution of other climatic and topo-edaphic variables in predicting S and PS for both components? (3) How does the taxonomic and phylogenetic turnover of annuals and perennials vary with spatial and environmental distances? We assessed annuals' and perennials' species richness, the phylogenetic structure at deep and shallow phylogenetic levels, and taxonomic and phylogenetic turnover along spatial and environmental distances. We found no relationship between annuals' richness and aridity, whereas perennials' richness showed a unimodal pattern. The phylogenetic structure of annuals and perennials showed contrasting responses to aridity and negatively correlated with topo-edaphic variables. We found phylogenetic clustering at intermediate-to-higher aridity levels for annuals, and at lower aridity levels for perennials. Both taxonomic and phylogenetic turnover in annuals and perennials correlated with the environmental distance rather than with spatial distance between communities, suggesting adaptation to local factors. Overall, our results show a decoupling in the response of the phylogenetic structure of annual and perennial components of plant communities to aridity in Mediterranean drylands. Our findings have significant implications for land management strategies under climate change.

Ecophylogenetics redux

Species' evolutionary histories shape their present-day ecologies, but the integration of phylogenetic approaches in ecology has had a contentious history. The field of ecophylogenetics promised to reveal the process of community assembly from simple indices of phylogenetic pairwise distances - communities shaped by environmental filtering were composed of closely related species, whereas communities shaped by competition were composed of less closely related species. However, the mapping of ecology onto phylogeny proved to be not so straightforward, and the field remains mired in controversy. Nonetheless, ecophylogenetic methods provided important advances across ecology. For example the phylogenetic distances between species is a strong predictor of pest and pathogen sharing, and can thus inform models of species invasion, coexistence and the disease dilution/amplification effect of biodiversity. The phylogenetic structure of communities may also provide information on niche space occupancy, helping interpret patterns of facilitation, succession and ecosystem functioning - with relevance for conservation and restoration - and the dynamics among species within foodwebs and metacommunities. I suggest leveraging advances in our understanding of the process of evolution on phylogenetic trees would allow the field to progress further, while maintaining the essence of the original vision that proved so seductive.

Effect of aridity on species assembly in gypsum drylands: a response mediated by the soil affinity of species

Previous studies found that plant communities on infertile soils are relatively resistant to climatic variation due to stress tolerance adaptations. However, the species assemblies in gypsum soil habitats require further investigation. Thus, we considered the following questions. (1) Do harsher arid conditions determine the characteristics of the species that form plant assemblages? (2) Is the selection of the species that assemble in arid conditions mediated by their ability to grow on gypsum soils? (3) Is the selection of species that assemble in harsher conditions related to phylogenetically conserved functional traits? Perennial plant communities were analysed in 89 gypsum-soil sites along a 400 km climate gradient from the central to southeastern Iberian Peninsula. Each local assemblage was analysed in 30 × 30 m plots and described based on taxonomic, functional (soil plant affinity) and phylogenetic parameters. The mean maximum temperatures in the hottest month, mean annual precipitation and their interaction terms were used as surrogates for the aridity conditions in generalized linear models. In the hottest locations, the gypsophily range narrowed and the mean gypsophily increased at the community level, thereby suggesting the filtering of species and the dominance of soil specialists in the actual plant assemblies. Drier sites had higher taxonomic diversity. The species that formed the perennial communities were close in evolutionary terms at the two ends of the aridity gradient. The mean maximum temperatures in the hottest month had the main abiotic filtering effect on perennial plant communities, which was mediated by the ability of species to grow on gypsum soils, and thus gypsum specialists dominated the species assemblies in the hottest locations. In contrast, the perennial communities on gypsum soils were relatively resistant to changes in precipitation. Our findings suggest that the warmer environmental conditions predicted by global change models will favour gypsum specialists over generalists.

Spatial patterns and interspecific relationships of two dominant cushion plants at three elevations on the Kunlun Mountain, China

One of the most important ecological processes is the formation of interspecific relationships in relation to spatial patterns among alpine cushion plants in extreme environmental habitats. However, such relationships remain poorly understood. Here, we examined the spatial patterns of alpine cushion plants along an altitudinal gradient of environmental severity and the interspecific relationship between two cushion species (Thylacospermum caespitosum and Androsace tangulashanensis) on the eastern Kunlun Mountain of China. Our results showed that the two species were highly aggregated within a distance of 2.5-5 m at the mid (S2) altitude, whereas they were randomly distributed at the low (S1) and high (S3) altitudes. A positive spatial interaction between the two species was observed over shorter distances at the mid (S2) altitude, and the spatial patterns were related to the size of individuals of the two species. Moreover, the impact of A. tangulashanensis on T. caespitosum (RII_{T. caespitosum}) was negative in all the study plots, and a positive impact of T. caespitosum on A. tangulashanensis (RII_{A. tangulashanensis}) was only observed at the mid (S2) altitude. Together, these results demonstrated that the spatial patterns of these two cushions varied with environmental severity, since the outcome of the interactions were different, to some extent, at the three altitudes. Plant size is the main factor affecting the spatial correlation and interspecific relationship between two cushions. Therefore, its potential influence should be considered when discussing interspecific relationships among cushions and their community construction at small scales in alpine ecosystems.

Assessment of Cu, Zn, Mn, and Fe enrichment in Mt. Kenya soils: evidence for atmospheric deposition and contamination

Mountains are the preferred sites for studying long-range atmospheric transportation and deposition of heavy metals, due to their isolation and steep temperature decrease that favors cold trapping and condensation of particulate forms of heavy metals. Any enrichment of heavy metals in mountains is presumed to primarily occur through atmospheric deposition. In this particular study, we assessed the status of 27 subsurface soils collected along two elevation gradients of Mt. Kenya using enrichment factors (EFs) as the ecological risk assessments. The collected soils were analyzed for total organic carbon, zinc (Zn), iron (Fe), manganese (Mn), and copper (Cu). The mean concentration of Mn, Fe, Zn, and Cu was 0.376 mg/kg, 47.6 mg/kg, 12.3 mg/kg, and 4.88 mg/kg in Chogoria and 0.560 mg/kg, 113 mg/kg, 12.7 mg/kg, and 2.70 mg/kg in Naro Moru respectively. These concentrations were below the US-EPA maximum permissible levels for soils, implying that the levels recorded had low toxicity. Meanwhile, the mean enrichment factors for Mn, Cu, and Zn were 0.447, 131, and 78.8 in Chogoria and 0.463, 38.9, and 53.0 in Naro Moru respectively. This implied that Zn and Cu in Chogoria sites were extremely enriched, while in Naro Moru, enrichment levels ranged from significant to extreme. However, Mn was found to have minimal enrichment in all the sites. Lower montane forest and bamboo zone recorded relatively high enrichment due to distance from source of pollution. Ericaceous zone also had high mean enrichment due to influence of wind which favors higher deposition at mid-elevations.

When facilitation meets clonal integration in forest canopies

Few studies have explored how - within the same system - clonality and positive plant-plant interactions might interact to regulate plant community composition. Canopy-dwelling epiphytes in species-rich forests provide an ideal system for studying this because many epiphytic vascular plants undertake clonal growth and because vascular epiphytes colonize canopy habitats after the formation of nonvascular epiphyte (i.e. bryophyte and lichen) mats. We investigated how clonal integration of seven dominant vascular epiphytes influenced inter-specific interactions between vascular epiphytes and nonvascular epiphytes in a subtropical montane moist forest in southwest China. Both clonal integration and environmental buffering from nonvascular epiphytes increased survival and growth of vascular epiphytes. The benefits of clonal integration for vascular epiphytes were higher when nonvascular epiphytes were removed. Similarly, facilitation from nonvascular epiphytes played a more important role when clonal integration of vascular epiphytes was eliminated. Overall, clonal integration had greater benefits than inter-specific facilitation. This study provides novel evidence for interactive effects of clonality and facilitation between vascular and nonvascular species, and has implications for our understanding of a wide range of ecosystems where both high levels of clonality and facilitation are expected to occur.

Is Phylogeny More Useful than Functional Traits for Assessing Diversity Patterns Under Community Assembly Processes?

Phylogenetic and functional diversities and their relationship are important for understanding community assembly, which relates to forest sustainability. Thus, both diversities have been used in ecological studies evaluating community responses to environmental changes. However, it is unclear whether these diversity measures can uncover the actual community assembly processes. Herein, we examined their utility to assess such assembly processes by analyzing similarities in phylogenetic, functional, and taxonomic α- and β-diversities along an elevational gradient. Additionally, we examined the relationships among environment, phylogeny, and functional traits within the community. Based on our results, we evaluated whether phylogenetic or functional diversity could better reveal the actual community assembly processes. We found that taxonomic, phylogenetic, and functional α-diversities were correlated with one another. Although the functional α-diversity showed a linear correlation with the elevational gradient, taxonomic and phylogenetic α-diversities showed unimodal patterns. Both phylogenetic and functional β-diversities correlated with taxonomic β-diversity, but there was no significant relationship between the former. Overall, our results evidenced that phylogenetic diversity and taxonomic diversity showed similar patterns, whereas functional diversity showed a relatively independent pattern, which may be due to limitations in the functional trait dimensions used in the present study. Although it is difficult to unravel whether the environment shapes phylogeny or functional traits within a community, phylogenetic diversity is a good proxy for assessing the assembly processes, whereas functional diversity may improve knowledge on the community by maximizing information about the functional trait dimensions.

The spatial frequency of climatic conditions affects niche composition and functional diversity of species assemblages: the case of Angiosperms

Climatic conditions vary in spatial frequency globally. Spatially rare climatic conditions provide fewer suitable environments than common ones and should impose constraints on the types of species present locally and regionally. We used data on 467 North American angiosperms to test the effects of the spatial frequency of climatic conditions on ecological niche specialisation and functional diversity. We predicted that rare climates should favour generalist species that are able to inhabit a broader range of climatic conditions. Our results show that climate frequency filters species that differ in niche breadths and rare environments host species combinations with greater functional diversity. The proposed analytical approaches and hypotheses can be adapted to investigate different aspects of ecological assemblies and their biodiversity. We discuss different mechanisms regarding how spatial frequency of environments can affect niche composition and functional diversity. These should be useful while developing theoretical frameworks for generating a deeper understanding of its underpinnings.

The dark side of facilitation: native shrubs facilitate exotic annuals more strongly than native annuals

Positive interactions enhance biodiversity and ecosystem function, but can also exacerbate biological invasions. Facilitation of exotic invaders by exotic foundation species (invasional meltdown) has been studied extensively, but facilitation of exotic invaders by native foundation species has attracted less attention. Specifically, very few studies have examined the extent that native foundation species facilitate native and exotic competitors. Understanding the processes that mediate interactions between native and exotic species can help explain, predict, and improve management of biological invasions. Here, we examined the effects of native foundation shrubs on the relative abundance of the annual plant community – including native and exotic taxa – from 2015–2018 in a desert ecosystem at Carrizo Plain National Monument, California, USA (elevation: 723 m). Shrub effects varied by year and by the identity of annual species, but shrubs consistently enhanced the abundance of the annual plant community and facilitated both native (n=17 species) and exotic (n=4 species) taxa. However, at the provenance level, exotic annuals were facilitated 2.75 times stronger in abundance than native annuals, and exotic annuals were always more abundant than natives both near and away from shrubs. Our study reaffirms facilitation as an important process in the organisation of plant communities and confirms that both native and exotic species can form positive associations with native foundation species. However, facilitation by native foundation species can exacerbate biological invasions by increasing the local abundance of exotic invaders. Thus, the force of facilitation can have a dark side relevant to ecosystem function and management.

Functionally distinct assembly of vascular plants colonizing alpine cushions suggests their vulnerability to climate change

BACKGROUND AND AIMS

Alpine cushion plants can initially facilitate other species during ecological succession, but later on can be negatively affected by their development, especially when beneficiaries possess traits allowing them to overrun their host. This can be reinforced by accelerated warming favouring competitively strong species over cold-adapted cushion specialists. However, little empirical research has addressed the trait-based mechanisms of these interactions. The ecological strategies of plants colonizing the cushion plant Thylacospermum caespitosum (Caryophyllaceae), a dominant pioneer of subnival zones, were studied in the Western Himalayas.

METHODS

To assess whether the cushion colonizers are phylogenetically and functionally distinct, 1668 vegetation samples were collected, both in open ground outside the cushions and inside their live and dead canopies, in two mountain ranges, Karakoram and Little Tibet. More than 50 plant traits related to growth, biomass allocation and resource acquisition were measured for target species, and the phylogenetic relationships of these species were studied [or determined].

KEY RESULTS

Species-based trait-environment analysis with phylogenetic correction showed that in both mountain ranges Thylacospermum colonizers are phylogenetically diverse but functionally similar and are functionally different from species preferring bare soil outside cushions. Successful colonizers are fast-growing, clonal graminoids and forbs, penetrating the cushion by rhizomes and stolons. They have higher root-to-shoot ratios, leaf nitrogen and phosphorus concentrations, and soil moisture and nutrient demands, sharing the syndrome of competitive species with broad elevation ranges typical of the late stages of primary succession. In contrast, the species from open ground have traits typical of stress-tolerant specialists from high and dry environments.

CONCLUSION

Species colonizing tight cushions of T. caespitosum are competitively strong graminoids and herbaceous perennials from alpine grasslands. Since climate change in the Himalayas favours these species, highly specialized subnival cushion plants may face intense competition and a greater risk of decline in the future.

Plant-plant interactions influence phylogenetic diversity at multiple spatial scales in a semi-arid mountain rangeland

Molecular phylogenies are increasingly used to understand how biotic interactions and environment shape phylogenetic community structure (PCS). However, we do not understand the effects of plant-plant interactions and environment on PCS and phylogenetic diversity across spatial scales, particularly in rangelands. Here, we ask: (1) do plant-plant interactions and environment affect PCS and phylogenetic diversity differently across the three spatial scales of the patch, the community, and the habitat? (2) What are the impacts of dominant cushion-nurse plants on the phylogenetic structure of plant communities? We assessed the PCS of semi-arid plant communities along an elevation gradient at the patch, community and habitat scales. Then, we assessed co-occurrence patterns along two sample slopes. Our results indicated important roles for biotic interactions and environmental filtering in determining phylogenetic diversity, with biotic interactions, in particular, having a stronger tendency to increase phylogenetic diversity. This is most likely due to the asymmetrical effects of nurse plants across the three spatial scales on our two different slopes. The impact of biotic interactions caused non-random phylogenetic patterns in more severe environments. In conclusion, biotic interactions influence phylogenetic diversity by altering PCS across aspects and along elevation gradients.

Differences in the fungal communities nursed by two genetic groups of the alpine cushion plant, Silene acaulis

Foundation plants shape the composition of local biotic communities and abiotic environments, but the impact of a plant's intraspecific variations on these processes is poorly understood. We examined these links in the alpine cushion moss campion (Silene acaulis) on two neighboring mountain ranges in the French Alps. Genotyping of cushion plants revealed two genetic clusters matching known subspecies. The exscapa subspecies was found on both limestone and granite, while the longiscapa one was only found on limestone. Even on similar limestone bedrock, cushion soils from the two S. acaulis subspecies deeply differed in their impact on soil abiotic conditions. They further strikingly differed from each other and from the surrounding bare soils in fungal community composition. Plant genotype variations accounted for a large part of the fungal composition variability in cushion soils, even when considering geography or soil chemistry, and particularly for the dominant molecular operational taxonomic units (MOTUs). Both saprophytic and biotrophic fungal taxa were related to the MOTUs recurrently associated with a single plant genetic cluster. Moreover, the putative phytopathogens were abundant, and within the same genus (Cladosporium) or species (Pyrenopeziza brassicae), MOTUs showing specificity for each plant subspecies were found. Our study highlights the combined influences of bedrock and plant genotype on fungal recruitment into cushion soils and suggests the coexistence of two mechanisms, an indirect selection resulting from the colonization of an engineered soil by free-living saprobes and a direct selection resulting from direct plant-fungi interactions.

Spatial scale-dependent phylogenetic signal in species distributions along geographic and elevation gradients in a mountainous rangeland

The mechanisms determining community phylogenetic structure range from local ecological mechanisms to broad biogeographical processes. How these community assembly processes determine phylogenetic structure and patterns in rangeland communities across multiple spatial scales is still poorly understood. We sought to determine whether the structure of herbaceous and shrub assemblages along local environmental gradients (elevation) and broad geography (latitude) exhibited phylogenetic signal at different spatial scales, across 2,500 ha of a mountainous rangeland. We analyzed species distribution and phylogenetic data at two spatial scales: the community level (1 m2 sample units obtained by stratified random sampling) and the habitat level (plant assemblages identified categorically based on environmental and geographical variables). We found significant phylogenetic signal in structure and pattern at both spatial scales, along local elevational, and latitudinal gradients. Moreover, beta diversity was affected by different environmental variables in herbaceous and shrub species distributions across different spatial scales. Our results highlight the relative importance of local ecological mechanisms, including niche-based deterministic processes (environmental filtering and species interactions) as well as those of biogeographical processes, such as stochastic dispersal limitation and habitat specialization in plant assemblages of mountainous rangeland.

Seasonal comparison of bacterial communities in rhizosphere of alpine cushion plants in the Himalayan Hengduan Mountains

Positive associations between alpine cushion plants and other species have been extensively studied. However, almost all studies have focused on the associations between macrofauna. Studies that have investigated positive associations between alpine cushion plants and rhizospheric microbes have been limited to the vegetation growing season. Here, we asked whether the positive effects that alpine cushion plants confer on rhizospheric microbe communities vary with seasons. We assessed seasonal variations in the bacterial diversity and composition in rhizosphere of two alpine cushion plants and surrounding bare ground by employing a high throughput sequencing method targeting the V3 region of bacterial 16S rRNA genes. Soil properties of the rhizosphere and the bare ground were also examined. We found that cushion rhizospheres harbored significantly more C, N, S, ammonia nitrogen, and soil moisture than the bare ground. Soil properties in cushion rhizospheres were not notably different, except for soil pH. Bacterial diversities within the same microhabitats did not vary significantly with seasons. We concluded that alpine cushion plants had positive effects on the rhizospheric bacterial communities, even though the strength of the effect varied in different cushion species. Cushion species and the soil sulfur content were probably the major factors driving the spatial distribution and structure of soil bacterial communities in the alpine communities dominated by cushion plants.

Thylacospermum caespitosum population structure and cushion species community diversity along an altitudinal gradient

As alpine plants, cushion species are particularly susceptible to environment changes. Thus, understanding population structure and community diversity variation of cushion plants along elevational gradients is crucial for estimating their response to predicted climate changes. In this study, Thylacospermum caespitosum populations from three elevations (low, medium, and high) in three climate zones of China (the Kunlun, Qilian, and Tianshan Mountains) were selected to evaluate the effect of elevation on the structure of T. caespitosum populations and species diversity of cushion communities. Results showed that elevation substantially influenced T. caespitosum populations (size structure, density, and death rate), as well as richness (α-diversity) and microhabitat species pool (species pool) of cushion communities. In the low elevations, T. caespitosum populations were in decline due to a lower ratio of small plants and higher mortality compared with populations at medium and high elevations. The α-diversity and species pool in cushion communities were significantly increased with decreased elevation, but the importance value of T. caespitosum decreased accordingly. Moreover, there was a significant positive correlation between elevation and relative importance value (the importance of one species in the community) of T. caespitosum (r = 0.883; P < 0.01). Elevation was significantly negatively correlated with the mortality rate of T. caespitosum (r = - 0.855; P < 0.01), α-diversity (r = - 0.933; P < 0.001), and species pool (r = - 0.885; P < 0.01). The declining characters of T. caespitosum population structure were obvious in low elevation populations. This decline may directly or indirectly relate to environmental change. Effects of elevation can provide an early indication of range contractions and population declines of cushion species with future climate warming. We call for more mechanistic studies of climate change impacts on cushion populations, particularly in alpine systems near the snow line.

Local trampling disturbance effects on alpine plant populations and communities: Negative implications for climate change vulnerability

Global change is modifying species communities from local to landscape scales, with alterations in the abiotic and biotic determinants of geographic range limits causing species range shifts along both latitudinal and elevational gradients. An important but often overlooked component of global change is the effect of anthropogenic disturbance, and how it interacts with the effects of climate to affect both species and communities, as well as interspecies interactions, such as facilitation and competition. We examined the effects of frequent human trampling disturbances on alpine plant communities in Switzerland, focusing on the elevational range of the widely distributed cushion plant Silene acaulis and the interactions of this facilitator species with other plants. Examining size distributions and densities, we found that disturbance appears to favor individual Silene growth at middle elevations. However, it has negative effects at the population level, as evidenced by a reduction in population density and reproductive indices. Disturbance synergistically interacts with the effects of elevation to reduce species richness at low and high elevations, an effect not mitigated by Silene. In fact, we find predominantly competitive interactions, both by Silene on its hosted and neighboring species and by neighboring (but not hosted) species on Silene. Our results indicate that disturbance can be beneficial for Silene individual performance, potentially through changes in its neighboring species community. However, possible reduced recruitment in disturbed areas could eventually lead to population declines. While other studies have shown that light to moderate disturbances can maintain high species diversity, our results emphasize that heavier disturbance reduces species richness, diversity, as well as percent cover, and adversely affects cushion plants and that these effects are not substantially reduced by plant-plant interactions. Heavily disturbed alpine systems could therefore be at greater risk for upward encroachment of lower elevation species in a warming world.

Frequency and intensity of facilitation reveal opposing patterns along a stress gradient

Disentangling the different processes structuring ecological communities is a long-standing challenge. In species-rich ecosystems, most emphasis has so far been given to environmental filtering and competition processes, while facilitative interactions between species remain insufficiently studied. Here, we propose an analysis framework that not only allows for identifying pairs of facilitating and facilitated species, but also estimates the strength of facilitation and its variation along environmental gradients. Our framework combines the analysis of both co-occurrence and co-abundance patterns using a moving window approach along environmental gradients to control for potentially confounding effects of environmental filtering in the co-abundance analysis. We first validate our new approach against community assembly simulations, and exemplify its potential on a large 1,134 plant community plots dataset. Our results generally show that facilitation intensity was strongest under cold stress, whereas the proportion of facilitating and facilitated species was higher under drought stress. Moreover, the functional distance between individual facilitated species and their facilitating species significantly changed along the temperature-moisture gradient, and seemed to influence facilitation intensity, although no general positive or general negative trend was discernible among species. The main advantages of our robust framework are as follows: It enables detecting facilitating and facilitated species in species-rich systems, and it allows identifying the directionality and intensity of facilitation in species pairs as well as its variation across long environmental gradients. It thus opens numerous opportunities for incorporating functional (and phylogenetic) information in the analysis of facilitation patterns. Our case study indicated high complexity in facilitative interactions across the stress gradient and revealed new evidence that facilitation, similarly to competition, can operate between functionally similar and dissimilar species. Extending the analyses to other taxa and ecosystems will foster our understanding how complex interspecific interactions promote biodiversity.

Nurse species and indirect facilitation through grazing drive plant community functional traits in tropical alpine peatlands

Facilitation among plants mediated by grazers occurs when an unpalatable plant extends its protection against grazing to another plant. This type of indirect facilitation impacts species coexistence and ecosystem functioning in a large array of ecosystems worldwide. It has nonetheless generally been understudied so far in comparison with the role played by direct facilitation among plants. We aimed at providing original data on indirect facilitation at the community scale to determine the extent to which indirect facilitation mediated by grazers can shape plant communities. Such experimental data are expected to contribute to refining the conceptual framework on plant–plant–herbivore interactions in stressful environments. We set up a 2‐year grazing exclusion experiment in tropical alpine peatlands in Bolivia. Those ecosystems depend entirely on a few, structuring cushion‐forming plants (hereafter referred to as “nurse” species), in which associated plant communities develop. Fences have been set over two nurse species with different strategies to cope with grazing (direct vs. indirect defenses), which are expected to lead to different intensities of indirect facilitation for the associated communities. We collected functional traits which are known to vary according to grazing pressure (LDMC, leaf thickness, and maximum height), on both the nurse and their associated plant communities in grazed (and therefore indirect facilitation as well) and ungrazed conditions. We found that the effect of indirectly facilitated on the associated plant communities depended on the functional trait considered. Indirect facilitation decreased the effects of grazing on species relative abundance, mean LDMC, and the convergence of the maximum height distribution of the associated communities, but did not affect mean height or cover. The identity of the nurse species and grazing jointly affected the structure of the associated plant community through indirect facilitation. Our results together with the existing literature suggest that the “grazer–nurse–beneficiary” interaction module can be more complex than expected when evaluated in the field.

Riders in the sky (islands): using a mega-phylogenetic approach to understand plant species distribution and coexistence at the altitudinal limits of angiosperm plant life

AIM

Plants occurring on high-alpine summits are generally expected to persist due to adaptations to extreme selective forces caused by the harshest climates where angiosperm life is known to thrive. We assessed the relative effects of this strong environmental filter and of other historical and stochastic factors driving plant community structure in very high-alpine conditions (up to 4,000m).

LOCATION

European Alps, Écrins National Park, France.

METHODS

Using species occurrence data collected from floristic surveys on 15 summits (2,791 m - 4,102 m a.s.l.) throughout the Écrins range, along with existing molecular sequence data obtained from GenBank, we used a mega-phylogenetic approach to evaluate the phylogenetic structure of high-alpine plant species assemblages. We used three nested species pools and two null models to address the importance of species-specific and species-neutral processes for driving coexistence.

RESULTS

Compared to the entire species pool of the study region, alpine summits exhibited a strong signal of phylogenetic clustering. Restricting statistical sampling to environmentally and historically defined species pools reduced the significance of this pattern. However, we could not reject a model that explicitly incorporates neutral colonization and local extinction in shaping community structure for dominant plant orders. Between summits, phylogenetic turnover was generally lower than expected. Environmental drivers did not explain overall phylogenetic patterns, but we found significant geographic and climatic structure in phylogenetic diversity at finer taxonomic scales.

MAIN CONCLUSIONS

Although we found evidence for strong phylogenetic clustering within alpine summits, we were not able to reject models of species-neutral processes to explain patterns of floristic diversity. Our results suggest that plant community structure in high-alpine regions can also be shaped by neutral processes, and not through the sole action of environmental selection as traditionally assumed for harsh and stressful environments.

Revisiting phylogenetic signal; strong or negligible impacts of polytomies and branch length information?

Background

Inaccurate estimates of phylogenetic signal may mislead interpretations of many ecological and evolutionary processes, and hence understanding where potential sources of uncertainty may lay has become a priority for comparative studies. Importantly, the sensitivity of phylogenetic signal indices and their associated statistical tests to incompletely resolved phylogenies and suboptimal branch-length information has been only partially investigated.

Methods

Here, we use simulations of trait evolution along phylogenetic trees to assess whether incompletely resolved phylogenies (polytomic chronograms) and phylogenies with suboptimal branch-length information (pseudo-chronograms) could produce directional biases in significance tests (p-values) associated with Blomberg et al.’s K and Pagel’s lambda (λ) statistics, two of the most widely used indices to measure and test phylogenetic signal. Specifically, we conducted pairwise comparisons between the p-values resulted from the use of “true” chronograms and their degraded counterparts (i.e. polytomic chronograms and pseudo-chronograms), and computed the frequency with which the null hypothesis of no phylogenetic signal was accepted using “true” chronograms but rejected when using their degraded counterparts (type I bias) and vice versa (type II bias).

Results

We found that the use of polytomic chronograms in combination with Blomberg et al.’s K resulted in both, clearly inflated estimates of phylogenetic signal and moderate levels of type I and II biases. More importantly, pseudo-chronograms led to high rates of type I biases. In contrast, Pagel’s λ was strongly robust to either incompletely resolved phylogenies and suboptimal branch-length information.

Conclusions

Our results suggest that pseudo-chronograms can lead to strong overestimation of phylogenetic signal when using Blomberg et al.’s K (i.e. high rates of type I biases), while polytomies may be a minor concern given other sources of uncertainty. In contrast, Pagel’s λ seems strongly robust to either incompletely resolved phylogenies and suboptimal branch-length information. Hence, Pagel’s λ may be a more appropriate alternative over Blomberg et al.’s K to measure and test phylogenetic signal in most ecologically relevant traits when phylogenetic information is incomplete.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0898-y) contains supplementary material, which is available to authorized users.

Evolutionary conservatism explains increasing relatedness of plant communities along a flooding gradient

Abiotic filters have been found either to increase or reduce evolutionary relatedness in plant communities, making it difficult to generalize responses of this major feature of biodiversity to future environmental change. Here, we hypothesized that the responses of phylogenetic structure to environmental change ultimately depend on how species have evolved traits for tolerating the resulting abiotic changes. Working within ephemeral wetlands, we tested whether species were increasingly related as flooding duration intensified. We also identified the mechanisms underlying increased relatedness by measuring root aerenchyma volume (RAV), a trait which promotes waterlogging tolerance. We found that species-specific responses to flooding explained most of the variation in occurrence for 63 vascular plant species across 5170 plots. For a subset of 22 species, we attributed these responses to variation in RAV. Large RAV specifically increased occurrence when flooding lasted for longer time periods, because large RAV reduced above-ground biomass loss. As large RAV was evolutionarily conserved within obligate wetland species, communities were more phylogenetically related as flooding increased. Our study shows how reconstructing the evolutionary history of traits that influence the responses of species to environmental change can help to predict future patterns in phylogenetic structure.

Functional Plant Types Drive Plant Interactions in a Mediterranean Mountain Range

Shrubs have positive (facilitation) and negative (competition) effects on understory plants, the net interaction effect being modulated by abiotic conditions. Overall shrubs influence to great extent the structure of plant communities where they have significant presence. Interactions in a plant community are quite diverse but little is known about their variability and effects at community level. Here we checked the effects of co-occurring shrub species from different functional types on a focal understory species, determining mechanisms driving interaction outcome, and tested whether effects measured on the focal species were a proxy for effects measured at the community level. Growth, physiological, and reproductive traits of Euphorbia nicaeensis, our focal species, were recorded on individuals growing in association with four dominant shrub species and in adjacent open areas. We also recorded community composition and environmental conditions in each microhabitat. Shrubs provided environmental conditions for plant growth, which contrasted with open areas, including moister soil, greater N content, higher air temperatures, and lower radiation. Shrub-associated individuals showed lower reproductive effort and greater allocation to growth, while most physiological traits remained unaffected. Euphorbia individuals were bigger and had more leaf N under N-fixing than under non-fixing species. Soil moisture was also higher under N-fixing shrubs; therefore soil conditions in the understory may counter reduced light conditions. There was a significant effect of species identity and functional types in the outcome of plant interactions with consistent effects at individual and community levels. The contrasting allocation strategies to reproduction and growth in Euphorbia plants, either associated or not with shrubs, showed high phenotypic plasticity and evidence its ability to cope with contrasting environmental conditions.

Plant-plant interactions as a mechanism structuring plant diversity in a Mediterranean semi-arid ecosystem

Plant-plant interactions are among the fundamental processes that shape structure and functioning of arid and semi-arid plant communities. Despite the large amount of studies that have assessed the relationship between plant-plant interactions (i.e., facilitation and competition) and diversity, often researchers forget a third kind of interaction, known as allelopathy. We examined the effect of plant-plant interactions of three dominant species: the perennial grass Lygeum spartum, the allelopathic dwarf shrub Artemisia herba-alba, and the nurse shrub Salsola vermiculata, on plant diversity and species composition in a semi-arid ecosystem in NE Spain. Specifically, we quantified the interaction outcome (IO) based on species co-occurrence, we analyzed diversity by calculation of the individual species-area relationship (ISAR), and compositional changes by calculation of the Chao-Jaccard similarity index. We found that S. vermiculata had more positive IO values than L. spartum, and A. herba-alba had values between them. Lygeum spartum and A. herba-alba acted as diversity repellers, whereas S. vermiculata acted as a diversity accumulator. As aridity increased, A. herba-alba transitioned from diversity repeller to neutral and S. vermiculata transitioned from neutral to diversity accumulator, while L. spartum remained as diversity repeller. Artemisia herba-alba had more perennial grass species in its local neighborhood than expected by the null model, suggesting some tolerance of this group to its "chemical neighbor". Consequently, species that coexist with A. herba-alba were very similar among different A. herba-alba individuals. Our findings highlight the role of the nurse shrub S. vermiculata as ecosystem engineer, creating and maintaining patches of diversity, as well as the complex mechanism that an allelopathic plant may have on diversity and species assemblage. Further research is needed to determine the relative importance of allelopathy and competition in the overall interference of allelopathic plants.

Recreational trails as a source of negative impacts on the persistence of keystone species and facilitation

Hiking trails, which are among the most common forms of infrastructure created for nature-based tourism, can alter key ecological processes. Trails can damage plants that facilitate the establishment and growth of other species leading to changes in community and functional composition. This can be a particular concern in harsh alpine ecosystems where plant communities are often dominated by one or two keystone species that provide shelter to a suite of beneficiary species. We analysed how a hiking trail affects interspecific facilitation by a dominant trampling-sensitive nurse shrub in the highest National Park in Australia. First we assessed the effects of the trail on the abundance, size and density of the nurse shrub at different distances from the trail. We then compared species richness and composition between areas in, and out, of the nurse shrub's canopy at different distances from the trail. To better understand why some species may benefit from facilitation and any effects of the trail on the quality of facilitation we compared functional composition between quadrats using community trait weighted means calculated by combining plant composition with species functional traits (canopy height, leaf area, % dry weight of leaves and specific leaf area). The abundance, size and density of nurse shrubs was lower on the trail edges than further away, particularly on the leeward edge, where there was more bare ground and less shrub cover. There were differences in species richness, cover, composition and functional composition in and outside the nurse shrub canopy. The shrubs appeared to facilitate species with more competitive, but less stress tolerant traits (e.g. taller plants with leaves that were larger, had high specific leaf area and low dry matter content). However, despite reductions in nurse shrubs near the trail, where they do exist, they appear to provide the same 'quality' of facilitation as nurse shrubs further away. However, longer-term effects may be occurring as the loss of nurse shrubs alters the wind profile of the ridgeline and therefore succession. The use of a steel mesh walkway along the trail may facilitate the regeneration of nurse shrubs and other plants that require protection from wind. Our results highlight the importance of diversifying recreation ecology research to assess how trails affect important ecological processes.