Specialization patterns in symbiotic associations: A community perspective over spatial scales

Specialization, contextualized in a resource axis of an organism niche, is a core concept in ecology. In biotic interactions, specialization can be determined by the range of interacting partners. Evolutionary and ecological factors, in combination with the surveyed scale (spatial, temporal, biological, and/or taxonomic), influence the conception of specialization. This study aimed to assess the specialization patterns and drivers in the lichen symbiosis, considering the interaction between the principal fungus (mycobiont) and the associated Nostoc (cyanobiont), from a community perspective considering different spatial scales. Thus, we determined Nostoc phylogroup richness and composition of lichen communities in 11 Nothofagus pumilio forests across a wide latitudinal gradient in Chile. To measure specialization, cyanobiont richness, Simpson's and d′ indices were estimated for 37 mycobiont species in these communities. Potential drivers that might shape Nostoc composition and specialization measures along the environmental gradient were analysed. Limitations in lichen distributional ranges due to the availability of their cyanobionts were studied. Turnover patterns of cyanobionts were identified at multiple spatial scales. The results showed that environmental factors shaped the Nostoc composition of these communities, thus limiting cyanobiont availability to establish the symbiotic association. Besides, specialization changed with the spatial scale and with the metric considered. Cyanolichens were more specialized than cephalolichens when considering partner richness and Simpson's index, whereas the d′ index was mostly explained by mycobiont identity. Little evidence of lichen distributional ranges due to the distribution of their cyanobionts was found. Thus, lichens with broad distributional ranges either associated with several cyanobionts or with widely distributed cyanobionts. Comparisons between local and regional scales showed a decreasing degree of specialization at larger scales due to an increase in cyanobiont richness. The results support the context dependency of specialization and how its consideration changes with the metric and the spatial scale considered. Subsequently, we suggest considering the entire community and widening the spatial scale studied as it is crucial to understand factors determining specialization.