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.

How do you solve a problem like Letharia? A new look at cryptic species in lichen-forming fungi using Bayesian clustering and SNPs from multilocus sequence data.

The inclusion of molecular data is increasingly an integral part of studies assessing species boundaries. Analyses based on predefined groups may obscure patterns of differentiation, and population assignment tests provide an alternative for identifying population structure and barriers to gene flow. In this study, we apply population assignment tests implemented in the programs STRUCTURE and BAPS to single nucleotide polymorphisms from DNA sequence data generated for three previous studies of the lichenized fungal genus Letharia. Previous molecular work employing a gene genealogical approach circumscribed six species-level lineages within the genus, four putative lineages within the nominal taxon L. columbiana (Nutt.) J.W. Thomson and two sorediate lineages. We show that Bayesian clustering implemented in the program STRUCTURE was generally able to recover the same six putative Letharia lineages. Population assignments were largely consistent across a range of scenarios, including: extensive amounts of missing data, the exclusion of SNPs from variable markers, and inferences based on SNPs from as few as three gene regions. While our study provided additional evidence corroborating the six candidate Letharia species, the equivalence of these genetic clusters with species-level lineages is uncertain due, in part, to limited phylogenetic signal. Furthermore, both the BAPS analysis and the ad hoc ΔK statistic from results of the STRUCTURE analysis suggest that population structure can possibly be captured with fewer genetic groups. Our findings also suggest that uneven sampling across taxa may be responsible for the contrasting inferences of population substructure. Our results consistently supported two distinct sorediate groups, 'L. lupina' and L. vulpina, and subtle morphological differences support this distinction. Similarly, the putative apotheciate species 'L. lucida' was also consistently supported as a distinct genetic cluster. However, additional studies will be required to elucidate the relationships of other L. columbiana s.l. populations with the two sorediate genetic clusters.