Diversity of Mountainous Lawn Grasses

How does the diversity of lawn grass communities vary across environmental gradients?

Fescue grass

Abstract

Lawn grass dominates urban and suburban landscapes (Ignatieva et al., 2025), sequestering atmospheric carbon and providing aesthetic grounds for leisure and wellbeing. However, increased urbanisation is driving global declines in cities’ green spaces (HUGSI 2024, pp.8). It is therefore important to study lawn-dominated landscapes and their capacity to respond to environmental change. The soil microbiome contains symbiont predictors of landscape ecological health, aiding with nutrient cycling, disease resistance, and water acquisition (Wagg et al., 2019). Plant biodiversity provides pollinator habitats and may mitigate urban heat (Francoeur et al., 2021). Little research investigates variables affecting lawn grass communities comprehensively. The objectives of this study were to identify bacterial, fungal, and dicot taxons associated with F.rubra, a common lawn grass, investigating changes in abundance and diversity along elevation, temperature, and moisture gradients. A small scale microhabitat was used to explain large scale macrohabitat variation; both had the same broad habitat (F.rubra grasslands) allowing data extrapolation. Genomic root endophyte data was obtained directly from F.rubra in the macrohabitat. Associated dicots were counted in the microhabitat. Results found precipitation significantly affected endophyte diversity, while temperature-precipitation interactions explained significant variation in endophytes and dicot diversity, varying along topographic gradients.

Keywords

Alkane hydroxylase, bioremediation, bioaugmentation, biosurfactants, hydrocarbonoclastic bacteria, rhamnolipids.

Introduction

Mountainous grasslands harbour a diverse range of habitats and biological niches as a result of altitude variability, often dominated by herbaceous and non-woody plants at high elevations (Schirpke et al., 2017). The Norwegian Scandes Mountains stretch to the Arctic Circle, where increased temperatures are melting the ice. The Scandinavian montane grasslands are a tundra ecoregion, regarded as high priority for conservation (WWF, 2000). As rising subarctic temperatures cause these mountains to warm faster than the global average (Lind et al., 2022), the impacts of climate change on their landscape biodiversity are becoming increasingly relevant.

Elevation gradients can explain biodiversity in terms of temperature lapse rate (TLR), whereby temperature decreases with elevation (Rolland, 2003). However the spatial distribution of biodiversity along elevation gradients is also affected by various factors including topography and interspecific interactions (Lomolino, 2001). Therefore, elevation gradients should be assessed alongside biogeographic factors to gain a comprehensive understanding of multivariate effects on regional and landscape biodiversity. For instance, topographic variation in tundra habitats (e.g. aspect, landscape ruggedness, urban planning) can drive localised temperature differences (Kiedrzyński et al., 2024; Aalto et al., 2022). Global warming is also causing higher shifts in precipitation from snowfall to rainfall in the Scandes (Lind et al., 2022), modifying the selective pressures faced by native species. This includes changes to soil moisture, temperature, and nutrient availability which influences the distribution of plant and soil microbiota communities (Zádrapová et al., 2024). The symbiotic endophytes associated with common florae provide essential metabolic resources to hosts, promoting plant growth and immunity (Omomowo et al., 2019). In particular, a group of endophytic fungi called Dark Septate Endophytes (DSE) has been identified in elevations up to 6150m (Kotilínek et al., 2017). The relative abundance, distribution, and diversity of plants and endophytic microbiota in such habitats can inform us of soil quality and ecosystem health. Growing up to altitudes of 3350m (St. John et al., 2012), Red Fescue (Festuca rubra) is a highly resilient cool season lawn grass, widespread across arctic and temperate zones of the Northern Hemisphere (Barkworth et al., 2007; POWO Kew, 2020). F. rubra can thrive in nutrient poor conditions (Shamim et al., 2021), inhabiting a range of climates from moist meadows to drier road verges and pastures. It is known to colonise disturbed habitats by spreading rhizomes (St. John et al., 2012) and outcompete other Festuca species, achieving greater land coverage over time. (Ferschl et al., 2024; Wallace et al., 2017). Drought tolerance exhibited by F. rubra (St John et al., 2012) suggests that diversity of associated bacteria and dicot species decreases with moisture, whereas endophytic fungi require higher moisture for growth and metabolism (Talley et al., 2002). The high prevalence of F. rubra is useful for data extrapolation across different grassland habitats.

Microclimate evaluation is a practical way to explain wider macroclimate heterogeneity (Dorigo et al., 2021). Little research exists regarding the effect of elevation gradients on biodiversity while also considering other correlates. Therefore we have used a microclimate model to examine the variable effects of precipitation, moisture and temperature on landscape and regional diversity. This can inform agricultural practice (Seoane et al., 2006), which is beneficial to local farming economies. Hence, there is both financial and ecological value to understanding complex relationships within montane grassland communities, and responses to environmental shifts under the growing threat of climate change.

Approach

Genomic data of endophytic symbionts associated with Festuca rubra (red fescue) were provided by a Norwegian research team, collected at 12 sites in the Scandes Mountains. Metagenomic sequences on 16S ribosomal RNAs and Internal Transcribed Spacers (ITS) were performed to identify taxonomic groups for bacteria and fungi, respectively. The output was 200 Amplicon Sequence Values (ASVs) for each location and domain. ASV abundances were estimated using a rarefaction analysis to account for varying sequence depths. A factorial design allowed integration of 2 IVs as discrete categories; Temperature (6.5, 8.5, 10.5 °C), and Precipitation (600, 1200, 2000, 2700 mm) using expected values at each site. Data was collected for each combination. For bacterial and fungal diversity, linear models were used to assess correlations between Shannon’s Diversity Index (H’) and temp or precip. Differences between categories, including precip-temp interactions, were analysed using 2-way ANOVA. TukeyHSD was used post-hoc for pairwise comparison. Assumptions of a 2-way ANOVA: (i) observations are independent, (ii) observations are normally distributed in every sample. (iii) variances of the samples are the same. Homoscedasticity (iii) and linearity of the model were checked visually using residuals vs fitted plots, and (ii) using a Shapiro-Wilk test. Taxonomies were then grouped by class instead of phylum to investigate correlations in precip / temp for specific subgroups. Multiple linear regressions were performed to identify these groups. P-values were adjusted for multiple comparisons. Correlations were visually assessed by examining residuals either side of a reference line (y = x).

Diversity of dicot species associated with lawn grass was assessed in a local sloped microclimate (F.rubra grassland, hypotenuse 12m), with elevation quantified as distance from road (m). Long-term records of soil moisture and temperature allowed correlation modelling between elevation (E), soil moisture (M) and temperature (T) to predict M and T at any given location. Random sampling with 0.5m quadrats provided data on the abundance (% cover) of 11 target herb species (Species A-K) using morphological identification. A specimen from each species was analysed using PCR and Sanger sequencing. Electropherogram outputs were analysed using BLAST to identify taxonomies. Regression analysis assessed the predictive models for the York study site. H’ was then calculated for each sample position. Linear regression analysis and model fitting was conducted to assess correlations between response (H’) and explanatory variables (E, M, T). Outliers were removed prior. Linear modelling was additionally used to assess correlations between (i) abundance of individual species and E, and (ii) abundance of species-species to identify relationships indicating shared environmental preferences.

Results

Fig. 1 Fungal diversity
Fig. 1 Effects of precip and temp on Fungal Diversity Data are Shannons’ Diversity Index (H’) where N=8 for each location. Violin plot [2a.] visualises slightly right-skewed data in group precip=2700, so the median was chosen for interaction plot visualisation [2b.] to reduce the effect of outliers. [2c.] represents the median and IQR within each condition. Significant p-values are illustrated. Data was normally distributed (Shapiro-Wilk, W = 0.98222, p = 0.2312).

Results from a 2-way ANOVA [Table 1.] found precip significantly affects fungal diversity (F(3)=4.7265; p<.005). There was insufficient evidence for temp affecting H’, however there were significant temp-precip interactions (F(6)=2.6862; p<.05). This is likely because temp affects evaporation rates and in turn, atmospheric moisture. TukeyHSD post-hocs found H’ was significantly higher when temp=10.5, precip=2000 versus precip=1200 [Fig. 1b]. No such differences were observed in lower temperatures. This suggests root symbiont fungi are more sensitive to changes in precipitation at high temperatures, and that fungal diversity is maximal in warm, moist conditions. Although fungal diversity was significantly affected by changes to precip / temp, we could not identify any specific fungal classes significantly associated with precip or temp. However, ASV abundance of the Exophiala genera was positively associated with temp (F(2,91)=10.16; R^2 = 0.1645; p<.0005).

ANOVA tables for fungi and bacteria Fig. 2 Bacterial diversity
Fig. 2 Effects of precip and temp on Bacterial Diversity Data are Shannons’ Diversity Index (H’) where N=8 for each location. Distribution was slightly left-skewed [3a], so median was chosen for interaction visualisation [3b.] to reduce the effect of outliers. [3c.] represents the median [Q1, Q3] and IQR within each condition. After removal of outliers using upper and lower IQR limits, Shapiro-Wilk found data was not normally distributed (W = 0.96409, p<.05). However, departure from normality was small and residuals were examined to follow a normal distribution, so assumptions were declared as partially met.

Multiple Linear Regression

2-way ANOVA [Table 2] found no significant effects of precip or temp on bacterial diversity alone, however interactions were significant at a 5% level. As normality was only partially met, this finding should be interpreted with a degree of caution. Bacilli (p<.05) and Alphaproteobacteria (p=0.02755189) were significantly correlated with precip and / or temp. Overall, precip and temp explained a significant amount of variation in Bacilli abundance (F(2,91) =6.931; R^2 =0.1131; p<.005 **) and Alphaproteobacteria (F(2,91) =6.678; R^2 =0.1088; p <.005). Bacilli was significantly positively correlated with temp, and significantly inversely correlated with precip [Table 3a.], suggesting a preference for warm, moderately dry conditions. However, visual assessment of the residuals plot for Bacilli suggested a linear model was not viable; an exponential model may have been a better fit. Hence, these correlations should be interpreted with caution as we cannot deduce meaningful relationships.

Alphaproteobacteria abundance shared no significant relationship with temp, however it was significantly inversely correlated with precip [Table 3b.]. Residuals supported the model fit. Subgroup order Rhizobiales was the most abundant, followed by Sphingomonadales. Both orders are frequently observed in Poaceae rhizospheres (Blaine et al., 2017; Jha et al., 2019). Examination of residuals suggested both linear models were accurate fits. Rhizobiales was negatively correlated with precip as seen in Table 3c. (F(2,91) =6.58; R2 =0.1071; p<.005). Although Sphingomonadales was inversely correlated with temp [Table 3d.], the regression model was not significant (F(2,91) =2.043; R2 = 0.02194; p>.05). It is therefore reasonable to assume order Rhizobiales contributes to the trend in precip observed in Alphaproteobacteria. Rhizobiales family Xanthobacteraceae occupied the same significant trend in precip (F(2,91) =5.705, t =-3.360, p<.005, R^2 adj=0.09189), suggesting this family of bacteria has a preference for, or is tolerant of, drier habitats.

Predictive correlations Observed correlations with diversity
Fig. 3. Correlations between environmental variables and dicot diversity. Data are the linear regressions for environmental variables. Predictive correlation models (top row) show the linear and non-linear relationships between explanatory variables, with annotated equations and R2 values. Correlations with diversity (H’) and Elevation (E), Soil Temperature (T) and Soil Moisture (M) are visualised (bottom row) after removal of outliers using Cook’s Distance method, N=77.

Linear microclimatic models of E and M explained 46.48% of variance in plant diversity data (F(1,77) = 68.74; p <.005), T explains 42.93% (F(1,77) = 59.67; p <.005). T was weakly positively correlated with H’, with a gradient of 0.1003, 95% CI [0.1258, 0.0742]. Correlations between H’ and E / M were comparably weaker. E against H’ exhibited a positive gradient of 0.04908, 95% CI [0.03729, 0.06087], while for M this was -0.03108, 95% CI [-0.02362, -0.03855]. This suggests temperature has the greatest effect on dicot diversity. However, given the interaction between E, T and M seen in the predictive correlation models [Fig. 4.], E and M may still contribute to this combined effect (e.g. greater elevation may increase sunlight exposure for photosynthesis, increasing plant metabolism). Such confounding variables may also explain the variation in data.

species-species correlations Species J
Fig. 6. Species J Abundance along an Elevation Gradient Data is % Abundance of Species J (N=84), significantly positively correlated with elevation (S = 51542; Rho = 0.4781642, p<.005).

A moderate positive correlation was identified between Species J (Centaurea nigra, Table 4.) and elevation [Fig. 6.]. The null hypothesis was rejected and the alternative accepted. This suggests growth of C. nigra improves at higher elevations, likely resulting from increased light availability and moderate-low moisture, its preferential conditions (Hill et al., 1999). Significant correlations were also identified between Species D-F, D-J, and C-I [Fig. 5.], indicating shared environmental preferences.

Discussion

The notion that fungi are more sensitive to changes in precipitation at high temperatures is supported by a meta-analysis (Yadav et al., 2025), which found significant seasonal differences in fungal colonisation frequencies; monsoon (80.9%), winter (61.8%) and summer (51.4%). Aligning with our findings, this is likely due to higher resource availability and decreased risk of desiccation in warm, moist versus warm, dry climates. The Exophilia genus, positively correlated with temp, is known to comprise DSE species tolerant of extreme temperatures (Oliveira et al., 2025; de León et al., 2026). However, this does not indicate high temperatures are preferential. An alternative explanation could be decreased competition with other fungi due to increased salinity, an indirect result of increased temperature due to evaporation (Khamidov et al., 2022). Existing literature does not align with our findings on bacteria. Longitudinal research (Zádrapová et al., 2024) has found significantly higher nutrient levels (Al, Fe, P, S) and Proteobacteria abundance in soils exposed to high versus low precipitation. This contrasts our findings of Alphaproteobacteria, indicating that nutrient availability induced by precip contributes to changes in microbial diversity. Alternatively, herbivory gradients may explain the variation in soil microbiota. All 3 Festuca species are major food sources for Scandinavian reindeer herds (Rangifer tarandus) during the feeding season (Granath et al., 2007). Seed dispersal brought about by this mutualism may cause distribution variance, which was not factored in this experiment.

Factorial design limited the parameters under which temp and precip could be examined. Consequently, the macroclimate model fails to predict changes in microbial abundance under extreme conditions. Genomic data also contained taxonomic errors such that 2 ASVs belonging to family Xanthobacteraceae were incorrectly classified as Rhodoplanes, indicating potential for critical Type I and II errors. Moreover, existing literature describes Epichlӧe festucae, a common fungal symbiont found in aerial organs of 65.7% of F.rubra (Pereira et al., 2019; Wallace et al. 2017). As genomic data came purely from root endophytes, no E. festucae ASVs were identified. Investigating aboveground microbiota such as E. festucae may have provided greater insight into how important symbioses are affected along environmental gradients.

Microclimate analysis provided interesting small-scale findings, using predictive correlations. Species J (C. nigra) dominates British NVC community MG5; grassland communities characterized by a preference for well-drained soils (Wallace et al., 2017). This aligns with correlates indicated by the elevation gradient [Fig. 6.]. However, moistures and temperatures exceeding optimum for C.nigra growth still reveal interacting confounders, including CO2 concentration, soil pH, and CO2:NO2 ratios (Quaderi et al., 2013). Both J and F were significantly correlated with D, but not one another, despite being known to share similar elevation and habitat preferences (Hill et al., 1999; ADHB, 2008). It is more likely their coexistence with D is due to shared resources than plant-plant interspecific mutualism (Li et al., 2025). Future investigation should assess effects of environmental variables in J, D, and F on a larger scale to explore these relationships more comprehensively. Elevation only explained 48.46% of microclimatic variance. The remainder is often an indirect result of broader factors (eg. precip influencing nutrient availability, temp influencing light availability). The model was also limited in that it could not uniformly explain and predict topographic variation in the Scandes; although temperature was correlated with elevation in the local microclimate, TLR may have the opposite large-scale impact (Rolland, 2003). The extent of this is particularly relevant as global warming is increasing upper elevation thresholds for montane grassland communities (Lamprecht et al., 2018). Future research should employ factorial design with wider parameters to examine variables independently and identify causal relationships. Herbivory gradients could be monitored via herd tracking. It is challenging and tedious to replicate and predict changes in mountainous landscapes without confounding influence. To overcome this, simulated microclimate modelling could be used (Hungerford et al., 1989) whereby statistical algorithms are developed to represent the actual conditions faced by organisms. The output can then be used to explain and predict large-scale biodiversity responses to climate change.

Overall, this study suggests diversity of rhizospheric bacteria and plants associated with lawn grass decreases with moisture. It revealed that while fungal endophytes prefer higher soil moisture contents, this can be limited by temperature, demonstrating the importance of multivariate analysis to understand changes in micro and macroclimate biodiversity along environmental gradients.

References

Aalto, J., Tyystjärvi, V., Niittynen, P., Kemppinen, J., Rissanen, T., Gregow, H. and Luoto, M. (2022). Microclimate temperature variations from boreal forests to the tundra. Agricultural and Forest Meteorology, 323(109037), pp.1–12. doi:https://doi.org/10.1016/j.agrformet.2022.109037.

AHDB Cereals & Oilseeds (2008). The encyclopaedia of arable weeds. [online] AHDB, Warwickshire: Agriculture and Horticulture Development Board, pp.1–222. Available at:

Bangroo, S.A., Najar, G.R. and Rasool, A. (2017). Effect of altitude and aspect on soil organic carbon and nitrogen stocks in the Himalayan Mawer Forest Range. CATENA, [online] 158, pp.63–68. doi:https://doi.org/10.1016/j.catena.2017.06.017.

Barkworth, M.E. (2007). Manual of Grasses for North America. Ogden: Utah State University Press, pp.1–640.

Blain, N.P., Helgason, B.L. and Germida, J.J. (2017). Endophytic root bacteria associated with the natural vegetation growing at the hydrocarbon-contaminated Bitumount Provincial Historic site. Canadian Journal of Microbiology, 63(6), pp.502–515. doi:https://doi.org/10.1139/cjm-2017-0039.

Choma, A., Komaniecka, I. and Zebracki, K. (2016). Structure, biosynthesis and function of unusual lipids A from nodule-inducing and N 2 -fixing bacteria. Biochimica et Biophysica Acta (BBA) - Molecular and Cell Biology of Lipids, 1862(2), pp.196–209. doi:https://doi.org/10.1016/j.bbalip.2016.11.004.

de León, L.R., González-Abradelo, D., Castillo-Marenco, T., Pérez-Llano, Y., Moreno-Perlin, T., Dávila-Ramos, S., Yarzábal Rodríguez, L.A., Sánchez-Carbente, M. del R., Fernández-Ocaña, A.M., Gostinčar, C., Gunde-Cimerman, N. and Batista-García, R.A. (2026). Exophiala as a polyextremotolerant fungal model: Ecology, adaptive strategies, and the emergence of ‘urban oligotrophic specialists’. Fungal Biology Reviews, 56(100487), pp.1–18. doi:https://doi.org/10.1016/j.fbr.2026.100487.

Dorigo, L., Boscutti, F. and Sigura, M. (2021). Landscape and microhabitat features determine small mammal abundance in forest patches in agricultural landscapes. PeerJ, 9(e12306), pp.1–26. doi:https://doi.org/10.7717/peerj.12306.

Hester, E., Vaksmaa, A., Valè, G., Monaco, S., Jetten, M. and Lüke, C. (2022). Effect of water management on microbial diversity and composition in an Italian rice field system. FEMS Microbiology Ecology, 98. doi: 10.1093/femsec/fiac018.

Hill, M.O., Mountford, J.O., Roy, D.B. and Bunce, R.G.H. (1999). Ellenberg’s indicator values for British plants. ECOFACT Volume 2 Technical Annex, 2a, pp.1–46.

HUGSI (2024). Urban Green Space Report 2024 HUGSI - Husqvarna Urban Green Space Insights present: How green are cities? Husqvarna, pp.1–50.

Hungerford, R.D., Nemani, R.R., Running, S.W. and Coughlan, J.C. (1989). MTCLIM: a mountain microclimate simulation model. Intermountain Research Station, INT-414, pp.1–56. doi:https://doi.org/10.2737/int-rp-414.

Ignatieva, M., Nielsen, S. and Martin, D.J. (2025). Recognising lawns as a part of ‘designed nature’. Pioneering study of lawn’s plant biodiversity in Australian context. Urban Ecosystems, 28(4), pp.1–21. doi:https://doi.org/10.1007/s11252-025-01745-z.

Ferschl, B., Szalai, M.Z., Gere, A., Kocsis, T. and Kotroczó, Z. (2024). Effects of Species-Rich Perennial Inter-Row Cover on Weed Flora and Soil Coverage in an Apple Orchard: A Case Study of Opportunities and Limitations in a Dry Continental Climate. Agronomy, [online] 14(11), p.2716. doi:https://doi.org/10.3390/agronomy14112716 [Journal uses article no. rather than page range].

Francoeur, X.W., Dagenais, D., Paquette, A., Dupras, J. and Messier, C. (2021). Complexifying the urban lawn improves heat mitigation and arthropod biodiversity. Urban Forestry & Urban Greening, 60(127007), p.127007. doi:https://doi.org/10.1016/j.ufug.2021.127007.

Granath, G., Vicari, M., Bazely, D.R., Ball, J.P., Puentes, A. and Tomo Rakocevic (2007). Variation in the abundance of fungal endophytes in fescue grasses along altitudinal and grazing gradients. Ecography, 30(3), pp.422–430. doi:https://doi.org/10.1111/j.0906-7590.2007.05027.x.

Jha, P.N., Gomaa, A.-B., Yanni, Y.G., El-Saadany, A.-E.Y., Stedtfeld, T.M., Stedtfeld, R.D., Gantner, S., Chai, B., Cole, J., Hashsham, S.A. and Dazzo, F.B. (2019). Alterations in the Endophyte-Enriched Root-Associated Microbiome of Rice Receiving Growth-Promoting Treatments of Urea Fertilizer and Rhizobium Biofertilizer. Microbial ecology, 79(2), pp.367–382. doi:https://doi.org/10.1007/s00248-019-01406-7.

Khamidov, M., Ishchanov, J., Hamidov, A., Donmez, C. and Djumaboev, K. (2022). Assessment of Soil Salinity Changes under the Climate Change in the Khorezm Region, Uzbekistan. International Journal of Environmental Research and Public Health, 19(14), pp.1–13. doi:https://doi.org/10.3390/ijerph19148794.

Kiedrzyński, M., Tomczyk, P.P., Zielińska, K.M., Kiedrzyńska, E. and Wąsowicz, P. (2024). Evidence of conservative range in mountain grasses during past climate change: Only contractions or local expansions possible. Global Ecology and Conservation, 51(e02889), pp.1–14. doi:https://doi.org/10.1016/j.gecco.2024.e02889.

Kotilínek, M., Hiiesalu, I., Košnar, J., Šmilauerová, M., Šmilauer, P., Altman, J., Dvorský, M., Kopecký, M. and Doležal, J. (2017). Fungal root symbionts of high-altitude vascular plants in the Himalayas. Scientific Reports, 7(1). doi:https://doi.org/10.1038/s41598-017-06938-x.

Kullman, L. (2018). A Review and Analysis of Factual Change on the Max Rise of the Swedish Scandes Treeline, in Relation to Climate Change over the Past 100 Years. Journal of Ecology & Natural Resources, 2(6), pp.1–16. doi:https://doi.org/10.23880/jenr-16000150.

Lamprecht, A., Semenchuk, P.R., Steinbauer, K., Winkler, M. and Pauli, H. (2018). Climate change leads to accelerated transformation of high-elevation vegetation in the central Alps. New Phytologist, 220(2), pp.447–459. doi:https://doi.org/10.1111/nph.15290.

Li, Z., Wei, J., Du, W., Huang, R., Song, L., Tian, Q. and Zhou, X. (2025). Environmental response strategies for the spatial distribution of seed plants in Gansu. Frontiers in plant science, 16(1526269), pp.1–18. doi:https://doi.org/10.3389/fpls.2025.1526269.

Lind, P., Belušić, D., Toivonen, E., Dobler, A., Pedersen, R.A., Wang, F., Matte, D., Kjellström, E., Landgren, O., Lindstedt, D., Christensen, O. and Jens Hesselbjerg Christensen (2022). Climate change information over Fenno-Scandinavia produced with a convection-permitting climate model. Climate Dynamics, 61(1-2), pp.519–541. doi:https://doi.org/10.1007/s00382-022-06589-3.

Lomolino, MARK.V. (2001). Elevation gradients of species-density: historical and prospective views. Global Ecology and Biogeography, 10(1), pp.3–13. doi:https://doi.org/10.1046/j.1466-822x.2001.00229.x.

Maciá-Vicente, J.G., Glynou, K. and Piepenbring, M. (2016). A new species of Exophiala associated with roots. Mycological Progress, 15(2), pp.1–12. doi:https://doi.org/10.1007/s11557-016-1161-4.

Oliveira, J.A.D and Pereira, O.L. (2025). Opening the black box of an underexplored and megadiverse group: history, taxonomy and functionality of dark septate endophytes (DSE). Fungal Biology Reviews, 54, p.100457. doi:https://doi.org/10.1016/j.fbr.2025.100457.

Omomowo, O.I. and Babalola, O.O. (2019). Bacterial and Fungal Endophytes: Tiny Giants with Immense Beneficial Potential for Plant Growth and Sustainable Agricultural Productivity. Microorganisms, 7(11), pp.1–15. doi:https://doi.org/10.3390/microorganisms7110481.

Pereira, E., Vázquez de Aldana, B.R., San Emeterio, L. and Zabalgogeazcoa, I. (2019). A Survey of Culturable Fungal Endophytes From Festuca rubra subsp. pruinosa, a Grass From Marine Cliffs, Reveals a Core Microbiome. Frontiers in Microbiology, 9(3321), pp.1–14. doi:https://doi.org/10.3389/fmicb.2018.03321.

POWO Kew (2020). Festuca rubra L. | Plants of the World Online | Kew Science. [online] Plants of the World Online. Available at: https://powo.science.kew.org/taxon/urn:lsid:ipni.org:names:30010036-2#publications [Accessed 2 May 2026].

Qaderi, M.M., Lynch, A.L., Godin, V.J. and Reid, D.M. (2013). Single and interactive effects of temperature, carbon dioxide, and watering regime on the invasive weed black knapweed (Centaurea nigra). Écoscience, 20(4), pp.328–338. doi:https://doi.org/10.2980/20-4-3631.

Rolland, C. (2003). Spatial and Seasonal Variations of Air Temperature Lapse Rates in Alpine Regions. Journal of Climate, [online] 16(7), pp.1032–1046. doi:https://doi.org/10.1175/1520-0442(2003)016%3C1032:SASVOA%3E2.0.CO;2.

Shamim, Z., Razzaq, H. and Awan, M.T. (2021). Wild Germplasm for Genetic Improvement in Crop Plants. Academic Press, Elsevier, pp.259–268.

St. John, L., Tilley, D., Hunt, P. and Wright, S. (2012). Plant Guide for Red Fescue (Festuca rubra). Plant Materials Centre, Aberdeen, Idaho: USDA-Natural Resources Conservation Service, pp.1–5.

Schirpke, U., Kohler, M., Leitinger, G., Fontana, V., Tasser, E. and Tappeiner, U. (2017). Future impacts of changing land-use and climate on ecosystem services of mountain grassland and their resilience. Ecosystem Services, 26, pp.79–94. doi:https://doi.org/10.1016/j.ecoser.2017.06.008.

Seoane, J., Justribó, J.H., García, F., Retamar, J., Rabadán, C. and Atienza, J.C. (2006). Habitat-suitability modelling to assess the effects of land-use changes on Dupont’s lark Chersophilus duponti: A case study in the Layna Important Bird Area. Biological Conservation, 128(2), pp.241–252. doi:https://doi.org/10.1016/j.biocon.2005.09.032.

Talley, S.M., Coley, P.D. and Kursar, T.A. (2002). The Effects of Weather on Fungal Abundance and Richness among 25 Communities in the Intermountain West. BMC Ecology, 2(7), pp.1–11. doi:https://doi.org/10.1186/1472-6785-2-7.

Wallace, H.L. and Prosser, M.V. (2017). A review of the National Vegetation Classification for the Calthion group of plant communities in England and Wales / Hilary Wallace and Mike Prosser, Ecological Surveys (Bangor) and Floodplain Meadows Partnership. Natural England Joint Publication ; JP021 ed. United Kingdom: Natural England, pp.1–106.

Wagg, C., Schlaeppi, K., Banerjee, S., Kuramae, E.E. and van der Heijden, M.G.A. (2019). Fungal-bacterial diversity and microbiome complexity predict ecosystem functioning. Nature Communications, 10(1), pp.1–10. doi:https://doi.org/10.1038/s41467-019-12798-y.

WWF (2000). High Arctic tundra | Ecoregions | WWF. [online] World Wildlife Fund Inc. Available at: https://www.worldwildlife.org/ecoregions/na1110 [Accessed 2 May 2026].

Yadav, G. and Meena, M. (2025). Unveiling the hidden culturable endophytic fungal diversity in aerial vegetative parts of Wrightia tinctoria (Roxb.) R.Br. of southern Aravalli hills. Scientific Reports, 15(1), pp.1–16. doi:https://doi.org/10.1038/s41598-025-10980-5.

Zádrapová, D., Chakraborty, A., Žáček, P., Korecký, J., Bhar, A. and Roy, A. (2024). Exploring the Rhizospheric Microbial Communities under Long-Term Precipitation Regime in Norway Spruce Seed Orchard. International Journal of Molecular Sciences, 25. doi: 10.3390/ijms25179658.