Fungal specificity and selectivity for algae play a major role in determining lichen partnerships across diverse ecogeographic regions in the lichen-forming family Parmeliaceae (Ascomycota).

Microbial symbionts are instrumental to the ecological and long-term evolutionary success of their hosts, and the central role of symbiotic interactions is increasingly recognized across the vast majority of life. Lichens provide an iconic group for investigating patterns in species interactions; however, relationships among lichen symbionts are often masked by uncertain species boundaries or an inability to reliably identify symbionts. The species-rich lichen-forming fungal family Parmeliaceae provides a diverse group for assessing patterns of interactions of algal symbionts, and our study addresses patterns of lichen symbiont interactions at the largest geographic and taxonomic scales attempted to date. We analysed a total of 2356 algal internal transcribed spacer (ITS) region sequences collected from lichens representing ten mycobiont genera in Parmeliaceae, two genera in Lecanoraceae and 26 cultured Trebouxia strains. Algal ITS sequences were grouped into operational taxonomic units (OTUs); we attempted to validate the evolutionary independence of a subset of the inferred OTUs using chloroplast and mitochondrial loci. We explored the patterns of symbiont interactions in these lichens based on ecogeographic distributions and mycobiont taxonomy. We found high levels of undescribed diversity in Trebouxia, broad distributions across distinct ecoregions for many photobiont OTUs and varying levels of mycobiont selectivity and specificity towards the photobiont. Based on these results, we conclude that fungal specificity and selectivity for algal partners play a major role in determining lichen partnerships, potentially superseding ecology, at least at the ecogeographic scale investigated here. To facilitate effective communication and consistency across future studies, we propose a provisional naming system for Trebouxia photobionts and provide representative sequences for each OTU circumscribed in this study.