Environment and host identity structure communities of green algal symbionts in lichens

Contrasting Symbiotic Patterns in Two Closely Related Lineages of Trimembered Lichens of the Genus Peltigera

Species circumscription is key to the characterization of patterns of specificity in symbiotic systems at a macroevolutionary scale. Here, a worldwide phylogenetic framework was used to assess the biodiversity and symbiotic patterns of association among partners in trimembered lichens from the genus Peltigera, section Chloropeltigera. We sequenced six loci of the main fungal partner and performed species discovery and validation analyses to establish putative species boundaries. Single locus phylogenies were used to establish the identity of both photobionts, Nostoc (cyanobacterium) and Coccomyxa (green alga). Distribution and specificity patterns were compared to the closely related clade, section Peltidea, which includes mainly Peltigera species with trimembered thalli. For section Chloropeltigera, eight fungal species (including five newly delimited putative species) were found in association with nine Nostoc phylogroups and two Coccomyxa species. In contrast, eight fungal species (including three newly delimited putative species) in section Peltidea were found in association with only four Nostoc phylogroups and the same two Coccomyxa species as for section Chloropeltigera. This difference in cyanobiont biodiversity between these two sections can potentially be explained by a significantly higher frequency of sexual reproductive structures in species from section Chloropeltigera compared to section Peltidea. Therefore, horizontal transmission of the cyanobiont might be more prevalent in Chloropeltigera species, while vertical transmission might be more common in Peltidea species. All Peltigera species in section Chloropeltigera are generalists in their association with Nostoc compared to more specialized Peltigera species in section Peltidea. Constrained distributions of Peltigera species that associate strictly with one species of green algae (Coccomyxa subellipsoidea) indicate that the availability of the green alga and the specificity of the interaction might be important factors limiting geographic ranges of trimembered Peltigera, in addition to constraints imposed by their interaction with Nostoc partners and by climatic factors.

The complexity of symbiotic interactions influences the ecological amplitude of the host: A case study in Stereocaulon (lichenized Ascomycota)

Symbiosis plays a fundamental role in nature. Lichens are among the best known, globally distributed symbiotic systems whose ecology is shaped by the requirements of all symbionts forming the holobiont. The widespread lichen-forming fungal genus Stereocaulon provides a suitable model to study the ecology of microscopic green algal symbionts (i.e., phycobionts) within the lichen symbiosis. We analysed 282 Stereocaulon specimens, collected in diverse habitats worldwide, using the algal ITS rDNA and actin gene sequences and fungal ITS rDNA sequences. Phylogenetic analyses revealed a great diversity among the predominant phycobionts. The algal genus Asterochloris (Trebouxiophyceae) was recovered in most sampled thalli, but two additional genera, Vulcanochloris and Chloroidium, were also found. We used variation-partitioning analyses to investigate the effects of climatic conditions, substrate/habitat characteristic, spatial distribution and mycobionts on phycobiont distribution. Based on an analogy, we examined the effects of climate, substrate/habitat, spatial distribution and phycobionts on mycobiont distribution. According to our analyses, the distribution of phycobionts is primarily driven by mycobionts and vice versa. Specificity and selectivity of both partners, as well as their ecological requirements and the width of their niches, vary significantly among the species-level lineages. We demonstrated that species-level lineages, which accept more symbiotic partners, have wider climatic niches, overlapping with the niches of their partners. Furthermore, the survival of lichens on substrates with high concentrations of heavy metals appears to be supported by their association with toxicity-tolerant phycobionts. In general, low specificity towards phycobionts allows the host to associate with ecologically diversified algae, thereby broadening its ecological amplitude.

Comparing Sanger sequencing and high-throughput metabarcoding for inferring photobiont diversity in lichens

The implementation of HTS (high-throughput sequencing) approaches is rapidly changing our understanding of the lichen symbiosis, by uncovering high bacterial and fungal diversity, which is often host-specific. Recently, HTS methods revealed the presence of multiple photobionts inside a single thallus in several lichen species. This differs from Sanger technology, which typically yields a single, unambiguous algal sequence per individual. Here we compared HTS and Sanger methods for estimating the diversity of green algal symbionts within lichen thalli using 240 lichen individuals belonging to two species of lichen-forming fungi. According to HTS data, Sanger technology consistently yielded the most abundant photobiont sequence in the sample. However, if the second most abundant photobiont exceeded 30% of the total HTS reads in a sample, Sanger sequencing generally failed. Our results suggest that most lichen individuals in the two analyzed species, Lasallia hispanica and L. pustulata, indeed contain a single, predominant green algal photobiont. We conclude that Sanger sequencing is a valid approach to detect the dominant photobionts in lichen individuals and populations. We discuss which research areas in lichen ecology and evolution will continue to benefit from Sanger sequencing, and which areas will profit from HTS approaches to assessing symbiont diversity.

Sharing of photobionts in sympatric populations of Thamnolia and Cetraria lichens: evidence from high-throughput sequencing

In this study, we explored the diversity of green algal symbionts (photobionts) in sympatric populations of the cosmopolitan lichen-forming fungi Thamnolia and Cetraria. We sequenced with both Sanger and Ion Torrent High-Throughput Sequencing technologies the photobiont ITS-region of 30 lichen thalli from two islands: Iceland and Öland. While Sanger recovered just one photobiont genotype from each thallus, the Ion Torrent data recovered 10-18 OTUs for each pool of 5 lichen thalli, suggesting that individual lichens can contain heterogeneous photobiont populations. Both methods showed evidence for photobiont sharing between Thamnolia and Cetraria on Iceland. In contrast, our data suggest that on Öland the two mycobionts associate with distinct photobiont communities, with few shared OTUs revealed by Ion Torrent sequencing. Furthermore, by comparing our sequences with public data, we identified closely related photobionts from geographically distant localities. Taken together, we suggest that the photobiont composition in Thamnolia and Cetraria results from both photobiont-mycobiont codispersal and local acquisition during mycobiont establishment and/or lichen growth. We hypothesize that this is a successful strategy for lichens to be flexible in the use of the most adapted photobiont for the environment.

Sequencing genomes from mixed DNA samples - evaluating the metagenome skimming approach in lichenized fungi

The metagenome skimming approach, i.e. low coverage shotgun sequencing of multi-species assemblages and subsequent reconstruction of individual genomes, is increasingly used for in-depth genomic characterization of ecological communities. This approach is a promising tool for reconstructing genomes of facultative symbionts, such as lichen-forming fungi, from metagenomic reads. However, no study has so far tested accuracy and completeness of assemblies based on metagenomic sequences compared to assemblies based on pure culture strains of lichenized fungi. Here we assembled the genomes of Evernia prunastri and Pseudevernia furfuracea based on metagenomic sequences derived from whole lichen thalli. We extracted fungal contigs using two different taxonomic binning methods, and performed gene prediction on the fungal contig subsets. We then assessed quality and completeness of the metagenome-based assemblies using genome assemblies as reference which are based on pure culture strains of the two fungal species. Our comparison showed that we were able to reconstruct fungal genomes from uncultured lichen thalli, and also cover most of the gene space (86-90%). Metagenome skimming will facilitate genome mining, comparative (phylo)genomics, and population genetics of lichen-forming fungi by circumventing the time-consuming, sometimes unfeasible, step of aposymbiotic cultivation.