Distribution of plant mycorrhizal traits along an elevational gradient does not fully mirror the latitudinal gradient

The influence of mycorrhizal symbiosis on ecosystem processes depends on the mycorrhizal type and status of plants. Early research hypothesized that the proportion of arbuscular mycorrhizal (AM) species decreases and of ectomycorrhizal (ECM) and ericoid mycorrhizal (ERM) species increases along increasing elevations and latitudes. However, there is very scarce information about this pattern along elevation gradients. We aimed to test this hypothesis and to describe the trends in plant mycorrhizal status by examining the Pyrenean mountain range (from 400 to 3400 m asl). The distribution of plant mycorrhizal types: AM, ECM, ERM, and non-mycorrhizal (NM) and status (obligately, OM, or facultatively, FM mycorrhizal plants, FM) were identified based on the Pyrenean Floristic Atlas and analyzed for climatic and edaphic drivers. The proportion of AM plants decreased slightly with elevation, while ECM species peaked at 1000 m asl. The proportion of ERM and NM plant species rose with increasing elevation. The proportion of FM species increased, and OM species decreased with increasing elevation. The change of AM and ECM species, and OM and FM species, along the elevational gradient, corresponds broadly to changes along the latitudinal gradient, driven by a combination of climatic and edaphic factors. Differently, the elevational occurrence of NM plant species is mainly driven only by climatic factors (low temperature) and that of ERM species by only edaphic factors (low pH). Large-scale macroecological studies (≥ 50 km grid cell) well reflect the effects of climate on the distribution of plant mycorrhizal traits, but local data (≤ 1 km grid cell) are needed to understand the effects of soil conditions and land use.