Keystone mutualisms, such as corals, lichens or mycorrhizae, sustain fundamental ecosystem functions. Range dynamics of these symbioses are, however, inherently difficult to predict because host species may switch between different symbiont partners in different environments, thereby altering the range of the mutualism as a functional unit. Biogeographic models of mutualisms thus have to consider both the ecological amplitudes of various symbiont partners and the abiotic conditions that trigger symbiont replacement. To address this challenge, we here investigate ‘symbiont turnover zones'––defined as demarcated regions where symbiont replacement is most likely to occur, as indicated by overlapping abundances of symbiont ecotypes. Mapping the distribution of algal symbionts from two species of lichen-forming fungi along four independent altitudinal gradients, we detected an abrupt and consistent β-diversity turnover suggesting parallel niche partitioning. Modelling contrasting environmental response functions obtained from latitudinal distributions of algal ecotypes consistently predicted a confined altitudinal turnover zone. In all gradients this symbiont turnover zone is characterized by approximately 12°C average annual temperature and approximately 5°C mean temperature of the coldest quarter, marking the transition from Mediterranean to cool temperate bioregions. Integrating the conditions of symbiont turnover into biogeographic models of mutualisms is an important step towards a comprehensive understanding of biodiversity dynamics under ongoing environmental change.

Holobionts are dynamic ecosystems that may respond to abiotic drivers with compositional changes. Uncovering elevational diversity patterns within these microecosystems can further our understanding of community‐environment interactions. Here, we assess how the major components of lichen holobionts—fungal hosts, green algal symbionts, and the bacterial community—collectively respond to an elevational gradient. We analyse populations of two lichen symbioses, Umbilicaria pustulata and U. hispanica, along an elevational gradient spanning 2100 altitudinal metres and covering three major biomes. Our study shows (i) discontinuous genomic variation in fungal hosts with one abrupt genomic differentiation within each of the two host species, (ii) altitudinally structured bacterial communities with pronounced turnover within and between hosts, and (iii) altitude‐specific presence of algal symbionts. Alpha diversity of bacterial communities decreased with increasing elevation. A marked turnover in holobiont diversity occurred across two altitudinal belts: at 11°C–13°C average annual temperature (here: 800–1200 m a.s.l.), and at 7°C–9°C average annual temperature (here: 1500–1800 m a.s.l.). The two observed zones mark a clustering of distribution limits and community shifts. The three ensuing altitudinal classes, that is, the most frequent combinations of species in holobionts, approximately correspond to the Mediterranean, cool‐temperate, and alpine climate zones. We conclude that multitrophic microecosystems, such as lichen holobionts, respond with concerted compositional changes to climatic factors that also structure communities of macroorganisms, for example, vascular plants.

Molecular methods are routinely used to estimate the effective size of populations (N_e). However, underlying model assumptions are frequently violated to an unknown extent. Although simulations can detect sources of bias and help to adjust sampling strategies and analyses methods, additional information from empirical data can also be used to calibrate methods and improve molecular N_e estimation methods. Here, we take advantage of long‐term genetic and ecological monitoring data of the grey wolf (Canis lupus) in Germany, and detailed population genetic studies in Poland, Spain and Portugal to improve N_e estimation strategies in this species, and species with similar life history traits. We first calculated N_e from average lifetime reproductive success and detailed census data from the German population, which served as a baseline to compare to molecular estimates based on linkage disequilibrium and sibship frequency. This yielded a robust N_e/N_c estimation that we used to calibrate molecular estimates of German, Polish and Iberian wolf populations. The linkage disequilibrium method was strongly influenced by spatial genetic structure, much more than the sibship frequency method. When N_e was estimated in local neighbourhoods, both methods yielded comparable results. Estimates of the metapopulation effective size seemed to correspond generally well with the sum of the estimates of local neighbourhoods. Overall, we found that the number of packs is a good proxy of the effective population size. Using this as a rule of thumb, we evaluated for all European wolf populations the N_e 500 indicator and concluded that half of the European wolf populations do not yet fulfil this criterion.

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host’s ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.