Environmental drivers and cryptic biodiversity hotspots define endophytes in Earth's largest terrestrial biome

Understanding how symbiotic associations differ across environmental gradients is key to predicting the fate of symbioses as environments change, and it is vital for detecting global reservoirs of symbiont biodiversity in a changing world.^1,2,3 However, sampling of symbiotic partners at the full-biome scale is difficult and rare. As Earth's largest terrestrial biome, boreal forests influence carbon dynamics and climate regulation at a planetary scale. Plants and lichens in this biome host the highest known phylogenetic diversity of fungal endophytes, which occur within healthy photosynthetic tissues and can influence hosts’ resilience to stress.^4,5 We examined how communities of endophytes are structured across the climate gradient of the boreal biome, focusing on the dominant plant and lichen species occurring across the entire south-to-north span of the boreal zone in eastern North America. Although often invoked for understanding the distribution of biodiversity, neither a latitudinal gradient nor mid-domain effect^5,6,7 can explain variation in endophyte diversity at this trans-biome scale. Instead, analyses considering shifts in forest characteristics, Picea biomass and age, and nutrients in host tissues from 46° to 58° N reveal strong and distinctive signatures of climate in defining endophyte assemblages in each host lineage. Host breadth of endophytes varies with climate factors, and biodiversity hotspots can be identified at plant-community transitions across the boreal zone at a global scale. Placed against a backdrop of global circumboreal sampling,⁴ our study reveals the sensitivity of endophytic fungi, their reservoirs of biodiversity, and their important symbiotic associations, to climate.