Multiple origins of high reciprocal symbiotic specificity at an intercontinental spatial scale among gelatinous lichens (Collemataceae, Lecanoromycetes)

Myrmecia, Not Asterochloris, Is the Main Photobiont of Cladonia subturgida (Cladoniaceae, Lecanoromycetes)

This study explores the diversity of photobionts associated with the Mediterranean lichen-forming fungus Cladonia subturgida. For this purpose, we sequenced the whole ITS rDNA region by Sanger using a metabarcoding method for ITS2. A total of 41 specimens from Greece, Italy, France, Portugal, and Spain were studied. Additionally, two specimens from Spain were used to generate four cultures. Our molecular studies showed that the genus Myrmecia is the main photobiont of C. subturgida throughout its geographic distribution. This result contrasts with previous studies, which indicated that the main photobiont for most Cladonia species is Asterochloris. The identity of Myrmecia was also confirmed by ultrastructural studies of photobionts within the lichen thalli and cultures. Photobiont cells showed a parietal chloroplast lacking a pyrenoid, which characterizes the species in this genus. Phylogenetic analyses indicate hidden diversity within this genus. The results of amplicon sequencing showed the presence of multiple ASVs in 58.3% of the specimens studied.

Cyanobacterial Variability in Lichen Cephalodia

The ecological success of lichens is related to both myco- and photobionts which condition the physiological limits of the lichen symbioses and thus affect their ecological niches and geographic ranges. A particular type of lichen, called cephalolichen, is characterized by housing both green algal and cyanobacterial symbionts-the latter is restricted to special structures called cephalodia. In this type of lichen, questions related to specialization within species or within individuals are still unsolved as different patterns have previously been observed. In order to study the variability at the intrathalline, intraspecific, and interspecific level, cyanobionts from different cephalodia within the same thalli and from different thalli were genetically analysed in three cephalolichen species at two different forests (18 thalli, 90 cephalodia). The results showed variability in the cephalodial Nostoc OTUs in all the studied species, both at the intrathalline and intraspecific levels. The variability of Nostoc OTUs found in different cephalodia of the same thallus suggests low specialization in this relationship. Additionally, differences in OTU diversity in the three studied species and in the two forests were found. The variability observed may confer an increased ecological plasticity and an advantage to colonize or persist under additional or novel habitats or conditions.

Specialization patterns in symbiotic associations: A community perspective over spatial scales

Specialization, contextualized in a resource axis of an organism niche, is a core concept in ecology. In biotic interactions, specialization can be determined by the range of interacting partners. Evolutionary and ecological factors, in combination with the surveyed scale (spatial, temporal, biological, and/or taxonomic), influence the conception of specialization. This study aimed to assess the specialization patterns and drivers in the lichen symbiosis, considering the interaction between the principal fungus (mycobiont) and the associated Nostoc (cyanobiont), from a community perspective considering different spatial scales. Thus, we determined Nostoc phylogroup richness and composition of lichen communities in 11 Nothofagus pumilio forests across a wide latitudinal gradient in Chile. To measure specialization, cyanobiont richness, Simpson's and d' indices were estimated for 37 mycobiont species in these communities. Potential drivers that might shape Nostoc composition and specialization measures along the environmental gradient were analysed. Limitations in lichen distributional ranges due to the availability of their cyanobionts were studied. Turnover patterns of cyanobionts were identified at multiple spatial scales. The results showed that environmental factors shaped the Nostoc composition of these communities, thus limiting cyanobiont availability to establish the symbiotic association. Besides, specialization changed with the spatial scale and with the metric considered. Cyanolichens were more specialized than cephalolichens when considering partner richness and Simpson's index, whereas the d' index was mostly explained by mycobiont identity. Little evidence of lichen distributional ranges due to the distribution of their cyanobionts was found. Thus, lichens with broad distributional ranges either associated with several cyanobionts or with widely distributed cyanobionts. Comparisons between local and regional scales showed a decreasing degree of specialization at larger scales due to an increase in cyanobiont richness. The results support the context dependency of specialization and how its consideration changes with the metric and the spatial scale considered. Subsequently, we suggest considering the entire community and widening the spatial scale studied as it is crucial to understand factors determining specialization.

Phylogeny and Ecology of Trebouxia Photobionts From Bolivian Lichens

In the past few years, new phylogenetic lineages in Trebouxia were detected as a result of molecular approaches. These studies included symbiont selectivity in lichen communities, transects along altitudinal gradients at local and global scales and the photobiont diversity in local populations of lichen-forming fungal species. In most of these studies, phylogenetic and haplotype analyses based on the internal transcribed spacer (ITS) locus have continuously allowed the recognition of new monophyletic lineages, which suggests that still numerous undiscovered Trebouxia lineages can be hidden in lichens from unexplored areas, especially in the tropics. Here, we estimated the biodiversity of photobionts in Bolivian Andean vegetation and assessed their specificity. About 403 lichen samples representing 42 genera, e.g., Haematomma, Heterodermia, Hypotrachyna, Lecanora, Lepra, Leucodermia, Parmotrema, Pertusaria, Polyblastidium, and Usnea, containing Trebouxia photobionts, were analyzed. ITS ribosomal DNA (rDNA) and rbcL markers were used. We obtained Trebouxia sequences from Bolivian samples belonging to already described clades A, C, I, and S. Thirty-nine Trebouxia lineages were distinguished within these clades, while 16 were new. To reveal the structure of the community of Bolivian photobionts and their relationships with mycobionts, the comparative effects of climate, altitude, geographical distances, substrate, and habitat type, as well as functional traits of lichens such as growth forms, propagation mode and secondary metabolites, were analyzed. Furthermore, new Bolivian records were included in analysis on a global scale. In our study, the mycobiont genus or even species are the most important factors correlated with photobiont identity. Moreover, we revealed that the community of Bolivian photobionts is shaped by altitude.

Biodiversity Patterns and Ecological Preferences of the Photobionts Associated With the Lichen-Forming Genus Parmelia

The worldwide, ecologically relevant lichen-forming genus Parmelia currently includes 41 accepted species, of which the Parmelia sulcata group (PSULgp) and the Parmelia saxatilis group (PSAXgp) have received considerable attention over recent decades; however, phycobiont diversity is poorly known in Parmelia s. lat. Here, we studied the diversity of Trebouxia microalgae associated with 159 thalli collected from 30 locations, including nine Parmelia spp.: P. barrenoae, P. encryptata, P. ernstiae, P. mayi, P. omphalodes, P. saxatilis, P. serrana, P. submontana, and P. sulcata. The mycobionts were studied by carrying out phylogenetic analyses of the nrITS. Microalgae genetic diversity was examined by using both nrITS and LSU rDNA markers. To evaluate putative species boundaries, three DNA species delimitation analyses were performed on Trebouxia and Parmelia. All analyses clustered the mycobionts into two main groups: PSULgp and PSAXgp. Species delimitation identified 13 fungal and 15 algal species-level lineages. To identify patterns in specificity and selectivity, the diversity and abundance of the phycobionts were identified for each Parmelia species. High specificity of each Parmelia group for a given Trebouxia clade was observed; PSULgp associated only with clade I and PSAXgp with clade S. However, the degree of specificity is different within each group, since the PSAXgp mycobionts were less specific and associated with 12 Trebouxia spp., meanwhile those of PSULgp interacted only with three Trebouxia spp. Variation-partitioning analyses were conducted to detect the relative contributions of climate, geography, and symbiotic partner to phycobiont and mycobiont distribution patterns. Both analyses explained unexpectedly high portions of variability (99 and 98%) and revealed strong correlations between the fungal and algal diversity. Network analysis discriminated seven ecological clusters. Even though climatic conditions explained the largest proportion of the variation among these clusters, they seemed to show indifference relative to climatic parameters. However, the cluster formed by P. saxatilis A/P. saxatilis B/Trebouxia sp. 2/Trebouxia sp. S02/Trebouxia sp. 3A was identified to prefer cold-temperate as well as humid summer environments.

Opening the Gap: Rare Lichens With Rare Cyanobionts - Unexpected Cyanobiont Diversity in Cyanobacterial Lichens of the Order Lichinales

The last decades of research led to a change in understanding of lichens that are now seen as self-sustaining micro-ecosystems, harboring diverse microbial organisms in tight but yet not fully understood relationships. Among the diverse interdependencies, the relationship between the myco- and photobiont is the most crucial, determining the shape, and ecophysiological properties of the symbiotic consortium. Roughly 10% of lichens associate with cyanobacteria as their primary photobiont, termed cyanolichens. Up to now, the diversity of cyanobionts of bipartite lichens resolved by modern phylogenetic approaches is restricted to the filamentous and heterocytous genera of the order Nostocales. Unicellular photobionts were placed in the orders Chroococcales, Pleurocapsales, and Chroococcidiopsidales. However, especially the phylogeny and taxonomy of the Chroococcidiopsidales genera remained rather unclear. Here we present new data on the identity and phylogeny of photobionts from cyanolichens of the genera Gonohymenia, Lichinella, Peccania, and Peltula from a broad geographical range. A polyphasic approach was used, combining morphological and cultivation-depending characteristics (microscopy, staining techniques, life cycle observation, baeocyte motility, and nitrogen fixation test) with phylogenetic analyses of the 16S rRNA and 16S–23S ITS gene region. We found an unexpectedly high cyanobiont diversity in the cyanobacterial lichens of the order Lichinales, including two new genera and seven new species, all of which were not previously perceived as lichen symbionts. As a result, we describe the novel unicellular Chroococcidiopsidales genera Pseudocyanosarcina gen. nov. with the species Pseudocyanosarcina phycocyania sp. nov. (from Peltula clavata, Australia) and Compactococcus gen. nov. with the species Compactococcus sarcinoides sp. nov. (from Gonohymenia sp., Australia) and the new Chroococcidiopsidales species Aliterella compacta sp. nov. (from Peltula clavata, Australia), Aliterella gigantea sp. nov. (from Peltula capensis; South Africa), Sinocapsa ellipsoidea sp. nov. (from Peccania cerebriformis, Austria), as well as the two new Nostocales species Komarekiella gloeocapsoidea sp. nov. (from Gonohymenia sp., Czechia) and Komarekiella globosa sp. nov. (from Lichinella cribellifera, Canary Islands, Spain). Our study highlights the role of cyanolichens acting as a key in untangling cyanobacterial taxonomy and diversity. With this study, we hope to stimulate further research on photobionts, especially of rare cyanolichens.

The diversity and community structure of symbiotic cyanobacteria in hornworts inferred from long-read amplicon sequencing

PREMISE

Nitrogen-fixing endosymbioses with cyanobacteria have evolved independently in five very different plant lineages. Expanding knowledge of these symbioses promises to improve the understanding of symbiosis evolution and broaden the toolkit for agricultural engineering to reduce artificial fertilizer use. Here we focused on hornworts, a bryophyte lineage in which all members host cyanobacteria, and investigated factors shaping the diversity of their cyanobiont communities.

METHODS

We sampled hornworts and adjacent soils in upstate New York throughout the hornwort growing season. We included all three sympatric hornwort species in the area, allowing us to directly compare partner selectivity. To profile cyanobacteria communities, we established a metabarcoding protocol targeting rbcL-X with PacBio long reads.

RESULTS

The hornwort cyanobionts detected were phylogenetically diverse, including clades that do not contain other known plant symbionts. We found significant overlap between hornwort cyanobionts and soil cyanobacteria, a pattern not previously reported in other plant-cyanobacteria symbioses. Cyanobiont communities differed between host plants only centimeters apart, but we did not detect an effect of sampling time or host species on the cyanobacterial community structure.

CONCLUSIONS

This study expands the phylogenetic diversity of known symbiotic cyanobacteria. Our analyses suggest that hornwort cyanobionts have a tight connection to the soil background, and we found no evidence that time within growing season, host species, or distance at the scale of meters strongly govern cyanobacteria community assembly. This study provides a critical foundation for further study of the ecology, evolution, and interaction dynamics of plant-cyanobacteria symbiosis.

Long-chain alkenes and alkadienes of eight lichen species collected in Japan

The hydrocarbons of eight lichen species isolated in Japan were analyzed, and diverse mono-, di-, and tri-unsaturated alkenes were detected. The positions of the double bonds of C₁₇ alkadienes (heptadecadiene) and C₁₇-C₂₀ alkenes were determined by mass spectrometry of their dimethyl disulfide adducts. We found that the six lichens containing green algal photobionts were distinguished by the presence of 1,8-heptadecadiene, 6,9-heptadecadiene, and 8- and 7-heptadecenes. On the other hand, 1-octadecene, 4-octadecene, and 5-nonadecene were the major alkene components of the two lichens with cyanobacterial photobionts. These alkadienes and alkenes were present in large quantities in the lichen samples. In particular, 1,8-heptadecadiene accounted for more than 90% of the total alkenes in all four lichens containing it. Our results provide new insights into the origin of C₁₇ alkadienes and C₁₇-C₂₀ alkenes in environmental and geological samples, and these alkenes can potentially be applied as lichen biomarkers.

Global Biodiversity Patterns of the Photobionts Associated with the Genus Cladonia (Lecanorales, Ascomycota)

The diversity of lichen photobionts is not fully known. We studied here the diversity of the photobionts associated with Cladonia, a sub-cosmopolitan genus ecologically important, whose photobionts belong to the green algae genus Asterochloris. The genetic diversity of Asterochloris was screened by using the ITS rDNA and actin type I regions in 223 specimens and 135 species of Cladonia collected all over the world. These data, added to those available in GenBank, were compiled in a dataset of altogether 545 Asterochloris sequences occurring in 172 species of Cladonia. A high diversity of Asterochloris associated with Cladonia was found. The commonest photobiont lineages associated with this genus are A. glomerata, A. italiana, and A. mediterranea. Analyses of partitioned variation were carried out in order to elucidate the relative influence on the photobiont genetic variation of the following factors: mycobiont identity, geographic distribution, climate, and mycobiont phylogeny. The mycobiont identity and climate were found to be the main drivers for the genetic variation of Asterochloris. The geographical distribution of the different Asterochloris lineages was described. Some lineages showed a clear dominance in one or several climatic regions. In addition, the specificity and the selectivity were studied for 18 species of Cladonia. Potentially specialist and generalist species of Cladonia were identified. A correlation was found between the sexual reproduction frequency of the host and the frequency of certain Asterochloris OTUs. Some Asterochloris lineages co-occur with higher frequency than randomly expected in the Cladonia species.

Complex Interaction Networks Among Cyanolichens of a Tropical Biodiversity Hotspot

Interactions within lichen communities include, in addition to close mutualistic associations between the main partners of specific lichen symbioses, also more elusive relationships between members of a wider symbiotic community. Here, we analyze association patterns of cyanolichen symbionts in the tropical montane forests of Taita Hills, southern Kenya, which is part of the Eastern Afromontane biodiversity hotspot. The cyanolichen specimens analyzed represent 74 mycobiont taxa within the order Peltigerales (Ascomycota), associating with 115 different variants of the photobionts genus Nostoc (Cyanobacteria). Our analysis demonstrates wide sharing of photobionts and reveals the presence of several photobiont-mediated lichen guilds. Over half of all mycobionts share photobionts with other fungal species, often from different genera or even families, while some others are strict specialists and exclusively associate with a single photobiont variant. The most extensive symbiont network involves 24 different fungal species from five genera associating with 38 Nostoc photobionts. The Nostoc photobionts belong to two main groups, the Nephroma-type Nostoc and the Collema/Peltigera-type Nostoc, and nearly all mycobionts associate only with variants of one group. Among the mycobionts, species that produce cephalodia and those without symbiotic propagules tend to be most promiscuous in photobiont choice. The extent of photobiont sharing and the structure of interaction networks differ dramatically between the two major photobiont-mediated guilds, being both more prevalent and nested among Nephroma guild fungi and more compartmentalized among Peltigera guild fungi. This presumably reflects differences in the ecological characteristics and/or requirements of the two main groups of photobionts. The same two groups of Nostoc have previously been identified from many lichens in various lichen-rich ecosystems in different parts of the world, indicating that photobiont sharing between fungal species is an integral part of lichen ecology globally. In many cases, symbiotically dispersing lichens can facilitate the dispersal of sexually reproducing species, promoting establishment and adaptation into new and marginal habitats and thus driving evolutionary diversification.

New lineages of photobionts in Bolivian lichens expand our knowledge on habitat preferences and distribution of Asterochloris algae

We studied the biodiversity of Asterochloris photobionts found in Bolivian lichens to better understand their global spatial distribution and adaptation strategies in the context of a worldwide phylogeny of the genus. Based on nuclear ITS rDNA, the chloroplast rbcL gene and the actin type I gene we reconstructed a phylogenetic tree that recovered nine new Asterochloris lineages, while 32 Bolivian photobiont samples were assigned to 12 previously recognized Asterochloris lineages. We also show that some previously discovered Asterochloris photobiont species and lineages may occur in a broader spectrum of climatic conditions, and mycobiont species and photobionts may show different preferences along an altitude gradient. To reveal general patterns of of mycobiont specificity towards the photobiont in Asterochloris, we tested the influence of climate, altitude, geographical distance and effects of symbiotic partner (mycobiont) at the species level of three genera of lichen forming fungi: Stereocaulon, Cladonia and Lepraria. Further, we compared the specificity of mycobionts towards Asterochloris photobionts in cosmopolitan, Neotropical, and Pantropical lichen forming fungi. Interestingly, cosmopolitan species showed the lowest specificity to their photobionts, but also the lowest haplotype diversity. Neotropical and Paleotropical mycobionts, however, were more specific.

Functional Traits in Lichen Ecology: A Review of Challenge and Opportunity

Community ecology has experienced a major transition, from a focus on patterns in taxonomic composition, to revealing the processes underlying community assembly through the analysis of species functional traits. The power of the functional trait approach is its generality, predictive capacity such as with respect to environmental change, and, through linkage of response and effect traits, the synthesis of community assembly with ecosystem function and services. Lichens are a potentially rich source of information about how traits govern community structure and function, thereby creating opportunity to better integrate lichens into 'mainstream' ecological studies, while lichen ecology and conservation can also benefit from using the trait approach as an investigative tool. This paper brings together a range of author perspectives to review the use of traits in lichenology, particularly with respect to European ecosystems from the Mediterranean to the Arctic-Alpine. It emphasizes the types of traits that lichenologists have used in their studies, both response and effect, the bundling of traits towards the evolution of life-history strategies, and the critical importance of scale (both spatial and temporal) in functional trait ecology.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

Extensive photobiont sharing in a rapidly radiating cyanolichen clade

Recent studies have uncovered remarkable diversity in Dictyonema s.lat. basidiolichens, here recognized as subtribe Dictyonemateae. This group includes five genera and 148 species, but hundreds more await description. The photobionts of these lichens belong to Rhizonema, a recently resurrected cyanobacterial genus known by a single species. To further investigate photobiont diversity within Dictyonemateae, we generated 765 new cyanobacterial sequences from 635 specimens collected from 18 countries. The ITS barcoding locus supported the recognition of 200 mycobiont (fungal) species among these samples, but the photobiont diversity was comparatively low. Our analyses revealed three main divisions of Rhizonema, with two repeatedly recovered as monophyletic (proposed as new species), and the third mostly paraphyletic. The paraphyletic lineage corresponds to R. interruptum and partnered with mycobionts from all five genera in Dictyonemateae. There was no evidence of photobiont-mycobiont co-speciation, but one of the monophyletic lineages of Rhizonema appears to partner predominantly with one of the two major clades of Cora (mycobiont) with samples collected largely from the northern Andes. Molecular clock estimations indicate the Rhizonema species are much older than the fungal species in the Dictyonemateae, suggesting that these basidiolichens obtained their photobionts from older ascolichen lineages and the photobiont variation in extant lineages of Dictyonemateae is the result of multiple photobiont switches. These results support the hypothesis of lichens representing "fungal farmers," in which diverse mycobiont lineages associate with a substantially lower diversity of photobionts by sharing those photobionts best suited for the lichen symbiosis among multiple and often unrelated mycobiont lineages.

Myco- and photobiont associations in crustose lichens in the McMurdo Dry Valleys (Antarctica) reveal high differentiation along an elevational gradient

Climatically extreme regions such as the polar deserts of the McMurdo Dry Valleys (78° S) in Continental Antarctica are key areas for a better understanding of changes in ecosystems. Therefore, it is particularly important to analyze and communicate current patterns of biodiversity in these sensitive areas, where precipitation mostly occurs in form of snow and liquid water is rare. Humidity provided by dew, clouds, and fog are the main water sources, especially for rock-dwelling crustose lichens as one of the most common vegetation-forming organisms. We investigated the diversity and interaction specificity of myco-/photobiont associations of 232 crustose lichen specimens, collected along an elevational gradient (171–959 m a.s.l.) within the McMurdo Dry Valleys. The mycobiont species and photobiont OTUs were identified by using three markers each (nrITS, mtSSU, RPB1, and nrITS, psbJ-L, COX2). Elevation, positively associated with water availability, turned out to be the key factor explaining most of the distribution patterns of the mycobionts. Pairwise comparisons showed Lecidea cancriformis and Rhizoplaca macleanii to be significantly more common at higher elevations and Carbonea vorticosa and Lecidea polypycnidophora at lower elevations. Lichen photobionts were dominated by the globally distributed Trebouxia OTU, Tr_A02 which occurred at all habitats. Network specialization resulting from myco-/photobiont bipartite network structure varied with elevation and associated abiotic factors. Along an elevational gradient, the spatial distribution, diversity, and genetic variability of the lichen symbionts appear to be mainly influenced by improved water relations at higher altitudes.

Cophylogenetic patterns in algal symbionts correlate with repeated symbiont switches during diversification and geographic expansion of lichen-forming fungi in the genus Sticta (Ascomycota, Peltigeraceae)

Species in the fungal genus Sticta form symbiotic associations primarily with either green algae or cyanobacteria, but tripartite associations or photosymbiodemes involving both types of photobionts occur in some species. Sticta is known to associate with green algae in the genus Symbiochloris. However, previous studies have shown that algae from other genera, such as Heveochlorella, may also be suitable partners for Sticta. We examined the diversity of green algal partners in the genus Sticta and assessed the patterns of association between the host fungus and its algal symbiont. We used multi-locus sequence data from multiple individuals collected in Australia, Cuba, Madagascar, Mauritius, New Zealand, Reunion and South America to infer phylogenies for fungal and algal partners and performed tests of congruence to assess coevolution between the partners. In addition, event-based methods were implemented to examine which cophylogenetic processes have led to the observed association patterns in Sticta and its green algal symbionts. Our results show that in addition to Symbiochloris, Sticta associates with green algae from the genera Chloroidium, Coccomyxa, Elliptochloris and Heveochlorella, the latter being the most common algal symbiont associated with Sticta in this study. Geography plays a strong role in shaping fungal-algal association patterns in Sticta as mycobionts associate with different algal lineages in different geographic locations. While fungal and algal phylogenies were mostly congruent, event-based methods did not find any evidence for cospeciation between the partners. Instead, the association patterns observed in Sticta and associated algae, were largely explained by other cophylogenetic events such as host-switches, losses of symbiont and failure of the symbiont to diverge with its host. Our results also show that tripartite associations with green algae evolved multiple times in Sticta.

Reproductive and dispersal strategies shape the diversity of mycobiont-photobiont association in Cladonia lichens

Ecological preferences, partner compatibility, or partner availability are known to be important factors shaping obligate and intimate lichen symbioses. We considered a complex of Cladonia species, traditionally differentiated by the extent of sexual reproduction and the type of vegetative propagules, to assess if the reproductive and dispersal strategies affect mycobiont-photobiont association patterns. In total 85 lichen thalli from 72 European localities were studied, two genetic markers for both Cladonia mycobionts and Asterochloris photobionts were analyzed. Variance partitioning analysis by multiple regression on distance matrices was performed to describe and partition variance in photobiont genetic diversity. Asexually reproducing Cladonia in our study were found to be strongly specific to their photobionts, associating with only two closely related Asterochloris species. In contrast, sexually reproducing lichens associated with seven unrelated Asterochloris lineages, thus being photobiont generalists. The reproductive mode had the largest explanatory power, explaining 44% of the total photobiont variability. Reproductive and dispersal strategies are the key factors shaping photobiont diversity in this group of Cladonia lichens. A strict photobiont specialisation observed in two studied species may steer both evolutionary flexibility and responses to ecological changes of these organisms, and considerably limit their distribution ranges.

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.

Bioclimatic factors at an intrabiome scale are more limiting than cyanobiont availability for the lichen-forming genus Peltigera

PREMISE OF THE STUDY

Factors shaping spatiotemporal patterns of associations in mutualistic systems are poorly understood. We used the lichen-forming fungi Peltigera and their cyanobacterial partners Nostoc to investigate the spatial structure of this symbiosis at an intrabiome scale and to identify potential factors shaping these associations.

METHODS

Ninety-three thalli were sampled in Québec, Canada, along a south-north and an east-west transect of ~1300 km each. We identified the two main partners (Peltigera species and Nostoc phylogroups) using molecular markers and modeled the effects of environmental variables and partner occurrence on Peltigera-Nostoc distributions.

KEY RESULTS

Peltigera species showed a high degree of specialization toward cyanobionts, whereas two Nostoc phylogroups dominated both transects by associating with several Peltigera species. Peltigera species had narrower ranges than these two main cyanobionts. Distributions of three Peltigera species were highly associated with precipitation and temperature variables, which was not detected for Nostoc phylogroups at this spatial scale.

CONCLUSIONS

For these cyanolichens, factors driving patterns of symbiotic associations are scale dependent. Contrary to global-scale findings, generalist Peltigera species were not more widespread within the boreal biome than specialists. Nostoc availability was not the only driver of Peltigera species' geographic ranges; environmental factors also contributed to their intrabiome distributions. Climatic conditions (especially precipitation) limited the range of some Peltigera species more than the range of their cyanobacterial partners at an intrabiome (boreal) scale.

Strong specificity and network modularity at a very fine phylogenetic scale in the lichen genus Peltigera

Identifying the drivers and evolutionary consequences of species interactions is a major goal of community ecology. Network-based analyses can provide mathematical tools to detect non-random patterns of interactions, and potentially help predicting the consequences of such patterns on evolutionary dynamics of symbiotic systems. Here, we characterize the structure of a lichen network at a very fine phylogenetic scale, by identifying the photosynthetic partners (i.e., cyanobacteria of the genus Nostoc) of lichenized fungi belonging to a monophyletic section of a single genus (i.e., section Polydactylon of the genus Peltigera), worldwide. Even at such a fine phylogenetic scale, we found that interactions were highly modular and anti-nested, indicating strong preferences in interactions. When considering local Peltigera communities, i.e., datasets at small spatial scales with only a slightly broader phylogenetic range, interactions remained modular but were asymmetric, with generalist Nostoc partners interacting with specialized Peltigera species. This asymmetry was not detected with our global spatial scale dataset. We discuss these results in the light of lichen community assembly, and explore how such interaction patterns may influence coevolution in lichens and the evolutionary stability of the mutualism in general.

Macroevolution of Specificity in Cyanolichens of the Genus Peltigera Section Polydactylon (Lecanoromycetes, Ascomycota)

Patterns of specificity among symbiotic partners are key to a comprehensive understanding of the evolution of symbiotic systems. Specificity of mutualistic partners, within a widespread monophyletic group for which all species are sampled has rarely been explored. Here, we assess the level of specificity between the cosmopolitan lichen-forming fungus (mycobiont) from the genus Peltigera, section Polydactylon, and its cyanobacterial partner Nostoc (cyanobiont). The mycobiont and cyanobiont phylogenies are inferred from five nuclear loci and the rbcLX region, respectively. These sequences were obtained from 206 lichen thalli, representing ca. 40 closely related Peltigera species sampled worldwide, doubling the number of known species in this group. We found a broad spectrum of specificity for both partners ranging from strict specialists to generalists. Overall, mycobionts are more specialized than cyanobionts by associating mostly with one or a few Nostoc phylogroups, whereas most cyanobionts associate frequently with several Peltigera species. Specialist mycobionts are older than generalists, supporting the hypothesis that specialization of mycobionts to one or few cyanobionts, is favored through time in geographic areas where species have been established for long periods of time. The relatively recent colonization of a new geographic area (Central and South America) by members of section Polydactylon is associated with a switch to a generalist pattern of association and an increased diversification rate by the fungal partner, suggesting that switches to generalism are rare events that are advantageous in new environments. We detected higher genetic diversity in generalist mycobionts. We also found that Peltigera species specialized on a single Nostoc phylogroup have narrower geographical distributions compared with generalist species.

Scaling up discovery of hidden diversity in fungi: impacts of barcoding approaches

The fungal kingdom is a hyperdiverse group of multicellular eukaryotes with profound impacts on human society and ecosystem function. The challenge of documenting and describing fungal diversity is exacerbated by their typically cryptic nature, their ability to produce seemingly unrelated morphologies from a single individual and their similarity in appearance to distantly related taxa. This multiplicity of hurdles resulted in the early adoption of DNA-based comparisons to study fungal diversity, including linking curated DNA sequence data to expertly identified voucher specimens. DNA-barcoding approaches in fungi were first applied in specimen-based studies for identification and discovery of taxonomic diversity, but are now widely deployed for community characterization based on sequencing of environmental samples. Collectively, fungal barcoding approaches have yielded important advances across biological scales and research applications, from taxonomic, ecological, industrial and health perspectives. A major outstanding issue is the growing problem of 'sequences without names' that are somewhat uncoupled from the traditional framework of fungal classification based on morphology and preserved specimens. This review summarizes some of the most significant impacts of fungal barcoding, its limitations, and progress towards the challenge of effective utilization of the exponentially growing volume of data gathered from high-throughput sequencing technologies.This article is part of the themed issue 'From DNA barcodes to biomes'.

Ecological Specialization of Two Photobiont-Specific Maritime Cyanolichen Species of the Genus Lichina

All fungi in the class Lichinomycetes are lichen-forming and exclusively associate with cyanobacteria. Two closely related maritime species of the genus Lichina (L. confinis and L. pygmaea) show similar distribution ranges in the Northeast Atlantic, commonly co-occurring at the same rocky shores but occupying different littoral zones. By means of 16S rRNA and phycocyanin operon markers we studied a) the phylogenetic relationships of cyanobionts associated with these species, b) the match of divergence times between both symbionts, and c) whether Lichina species differ in photobiont association and in how geography and ecology affect selectivity. The cyanobionts studied are closely related to both marine and freshwater strains of the genus Rivularia. We found evidence of a high specificity to particular cyanobiont lineages in both species: Lichina pygmaea and L. confinis incorporate specific lineages of Rivularia that do not overlap at the haplotype nor the OTU levels. Dating divergences of the fungal and cyanobacterial partners revealed an asynchronous origin of both lineages. Within each fungal species, selectivity varied across the studied area, influenced by environmental conditions (both atmospheric and marine), although patterns were highly correlated between both lichen taxa. Ecological speciation due to the differential association of photobionts to each littoral zone is suspected to have occurred in marine Lichina.

A multigene phylogenetic synthesis for the class Lecanoromycetes (Ascomycota): 1307 fungi representing 1139 infrageneric taxa, 317 genera and 66 families

The Lecanoromycetes is the largest class of lichenized Fungi, and one of the most species-rich classes in the kingdom. Here we provide a multigene phylogenetic synthesis (using three ribosomal RNA-coding and two protein-coding genes) of the Lecanoromycetes based on 642 newly generated and 3329 publicly available sequences representing 1139 taxa, 317 genera, 66 families, 17 orders and five subclasses (four currently recognized: Acarosporomycetidae, Lecanoromycetidae, Ostropomycetidae, Umbilicariomycetidae; and one provisionarily recognized, 'Candelariomycetidae'). Maximum likelihood phylogenetic analyses on four multigene datasets assembled using a cumulative supermatrix approach with a progressively higher number of species and missing data (5-gene, 5+4-gene, 5+4+3-gene and 5+4+3+2-gene datasets) show that the current classification includes non-monophyletic taxa at various ranks, which need to be recircumscribed and require revisionary treatments based on denser taxon sampling and more loci. Two newly circumscribed orders (Arctomiales and Hymeneliales in the Ostropomycetidae) and three families (Ramboldiaceae and Psilolechiaceae in the Lecanorales, and Strangosporaceae in the Lecanoromycetes inc. sed.) are introduced. The potential resurrection of the families Eigleraceae and Lopadiaceae is considered here to alleviate phylogenetic and classification disparities. An overview of the photobionts associated with the main fungal lineages in the Lecanoromycetes based on available published records is provided. A revised schematic classification at the family level in the phylogenetic context of widely accepted and newly revealed relationships across Lecanoromycetes is included. The cumulative addition of taxa with an increasing amount of missing data (i.e., a cumulative supermatrix approach, starting with taxa for which sequences were available for all five targeted genes and ending with the addition of taxa for which only two genes have been sequenced) revealed relatively stable relationships for many families and orders. However, the increasing number of taxa without the addition of more loci also resulted in an expected substantial loss of phylogenetic resolving power and support (especially for deep phylogenetic relationships), potentially including the misplacements of several taxa. Future phylogenetic analyses should include additional single copy protein-coding markers in order to improve the tree of the Lecanoromycetes. As part of this study, a new module ("Hypha") of the freely available Mesquite software was developed to compare and display the internodal support values derived from this cumulative supermatrix approach.

Molecular phylogeny and symbiotic selectivity of the green algal genus Dictyochloropsis s.l. (Trebouxiophyceae): a polyphyletic and widespread group forming photobiont-mediated guilds in the lichen family Lobariaceae

Dictyochloropsis s.l. is an ecologically important, common but little-studied genus of green algae. Here, we examined the diversity and host selectivity of algae attributed to this genus at both species-to-species and species-to-community levels. We conducted a molecular investigation of 15 cultured strains and several lichen photobionts, using 18S rRNA, rbcL and ITS sequence data. We further used seven alga-specific microsatellite markers to study algal sharing among fungi of the family Lobariaceae in two populations in Madeira and Taiwan (454 lichens). We found that the genus Dictyochloropsis s.l. is polyphyletic. Dictyochloropsis clade 1 comprises only free-living algae whereas Dictyochloropsis clade 2 includes lichenized algae as well as free-living algae. Fungal selectivity towards algae belonging to Dictyochloropsis clade 2 is high. Selectivity varies geographically, with photobionts being restricted to a single region. Finally, we showed that Dictyochloropsis clade 2 individuals are shared among different fungal hosts in communities of lichens of the Lobariaceae. As for other green algal lineages, there is a high amount of cryptic diversity in Dictyochloropsis. Furthermore, co-evolution between Dictyochloropsis clade 2 algae and representatives of the Lobariaceae is manifested at the community level, with several unrelated fungal species being horizontally connected by shared photobiont clones.

Do photobiont switch and cephalodia emancipation act as evolutionary drivers in the lichen symbiosis? A case study in the Pannariaceae (Peltigerales)

Lichen symbioses in the Pannariaceae associate an ascomycete and either cyanobacteria alone (usually Nostoc; bipartite thalli) or green algae and cyanobacteria (cyanobacteria being located in dedicated structures called cephalodia; tripartite thalli) as photosynthetic partners (photobionts). In bipartite thalli, cyanobacteria can either be restricted to a well-delimited layer within the thallus ('pannarioid' thalli) or spread over the thallus that becomes gelatinous when wet ('collematoid' thalli). We studied the collematoid genera Kroswia and Physma and an undescribed tripartite species along with representatives of the pannarioid genera Fuscopannaria, Pannaria and Parmeliella. Molecular inferences from 4 loci for the fungus and 1 locus for the photobiont and statistical analyses within a phylogenetic framework support the following: (a) several switches from pannarioid to collematoid thalli occured and are correlated with photobiont switches; the collematoid genus Kroswia is nested within the pannarioid genus Fuscopannaria and the collematoid genus Physma is sister to the pannarioid Parmeliella mariana group; (b) Nostoc associated with collematoid thalli in the Pannariaceae are related to that of the Collemataceae (which contains only collematoid thalli), and never associated with pannarioid thalli; Nostoc associated with pannarioid thalli also associate in other families with similar morphology; (c) ancestors of several lineages in the Pannariaceae developed tripartite thalli, bipartite thalli probably resulting from cephalodia emancipation from tripartite thalli which eventually evolved and diverged, as suggested by the same Nostoc present in the collematoid genus Physma and in the cephalodia of a closely related tripartite species; Photobiont switches and cephalodia emancipation followed by divergence are thus suspected to act as evolutionary drivers in the family Pannariaceae.

High photobiont diversity in the common European soil crust lichen Psora decipiens

The genetic diversity of green algal photobionts (chlorobionts) in soil crust forming lichens was studied as part of the SCIN-project (Soil Crust InterNational). A total of 64 lichen samples were collected from four different sites along latitudinal and altitudinal gradients in Europe (Tabernas/Spain; Hochtor-Großglockner/Austria; Gynge Alvar/Sweden; Ruine Homburg/Germany). The dominant lichen species at all four sites was Psora decipiens, often occurring with Buellia elegans, Fulgensia bracteata, F. fulgens and Peltigera rufescens. Genetic identification of chlorobionts was carried out using the nuclear marker (nrITS) and a chloroplast marker (psbL-J). We found P. decipiens to be associated with several different species of Trebouxia and Asterochloris, although previously described to only have Asterochloris sp. The phylogenetic analyses revealed a high chlorobiont diversity with 12 well supported clades, including Trebouxia asymmetrica, T. jamesii, T. impressa and other, as yet taxonomically unidentified clades (Trebouxia sp. URa1-4, T. sp. URa6, T. sp. URa7-13). Additionally, five clades of Asterochloris were identified (A. magna, A. sp. URa14 -17). Most of the chlorobiont species appeared to be cosmopolitan, but five clades were unevenly distributed between the sampling sites with only Trebouxia being found in the warm and dry Spanish habitats and combinations of Trebouxia and Asterochloris in the cooler and more humid habitats. The wide range of chlorobiont species might contribute to the observed domination of P. decipiens at all four research sites of the SCIN project which range from a desert in Spain to an alpine site in the Alps of Austria.

Assessing population structure and host specialization in lichenized cyanobacteria

Coevolutionary theory predicts that the distribution of obligately symbiotic organisms will be determined by the dispersal ability and ecological range of both partners. We examined this prediction for lichen-forming fungi that form obligate symbioses with cyanobacteria. We compared genotypes of both partners of 250 lichens collected at multiple spatial scales in British Columbia, Canada. Multilocus sequence data collected from a subset of 128 of the specimens were used to determine the degree of recombination within the cyanobacterial populations. We found that six distinct clusters of cyanobacterial genotypes are distributed throughout the known global phylogeny of the genus Nostoc, and that each appears to be evolving clonally. Fungal specialization is high, with each species associating with either one or two of the cyanobacterial clusters, while cyanobacterial specialization varies, with clusters associating with between one and 12 different fungal species. Specialization also varies geographically, with some combinations restricted to a single site despite the availability of both partners elsewhere. Photobiont association patterns are determined by a combination of genetically based specificity, spatial population structure, and ecological factors and cannot be easily predicted by photobiont dispersal syndromes.

LEC-2, a highly variable lectin in the lichen Peltigera membranacea

Lectins are a diverse group of carbohydrate binding proteins often involved in cellular interactions. A lectin gene, lec-2, was identified in the mycobiont of the lichen Peltigera membranacea. Sequencing of lec-2 open reading frames from 21 individual samples showed an unexpectedly high level of polymorphism in the deduced protein (LEC-2), which was sorted into nine haplotypes based on amino acid sequence. Calculations showed that the rates of nonsynonymous versus synonymous nucleotide substitutions deviated significantly from the null hypothesis of neutrality, indicating strong positive selection. Molecular modeling revealed that most amino acid replacements were around the putative carbohydrate-binding pocket, indicating changes in ligand binding. Lectins have been thought to be involved in the recognition of photobiont partners in lichen symbioses, and the hypothesis that positive selection of LEC-2 is driven by variation in the Nostoc photobiont partner was tested by comparing mycobiont LEC-2 haplotypes and photobiont genotypes, as represented by the rbcLX region. It was not possible to pair up the two types of marker sequences without conflicts, suggesting that positive selection of LEC-2 was not due to variation in photobiont partners.

Multiple-scale environmental modulation of lichen reproduction

It is necessary to understand how environmental changes affect plant fitness to predict survival of a species, but this knowledge is scarce for lichens and complicated by their formation of sexual and asexual reproductive structures. Are the presence and number of reproductive structures in Lobaria pulmonaria, a threatened lichen, dependent on thallus size, and is their formation sequential? Does any size-dependence and sequential formation vary along a climate gradient? Generalized linear mixed models were used to explore the effect of environmental predictors on the size and presence/abundance of each reproductive structure and to determine the probability of a given-sized thallus to develop any reproductive structure. The largest individuals are more likely to develop reproductive structures, and the lichen uses a mixed strategy of early asexual reproduction and late sexual. Macro and microclimatic variables also influenced reproductive capacity. Relationships among climate conditions and lichen size and reproductive capacity can compromise the future viability of the species in the most southern populations of Europe.

Geographic mosaic of symbiont selectivity in a genus of epiphytic cyanolichens

In symbiotic systems, patterns of symbiont diversity and selectivity are crucial for the understanding of fundamental ecological processes such as dispersal and establishment. The lichen genus Nephroma (Peltigerales, Ascomycota) has a nearly cosmopolitan distribution and is thus an attractive model for the study of symbiotic interactions over a wide range of spatial scales. In this study, we analyze the genetic diversity of Nephroma mycobionts and their associated Nostoc photobionts within a global framework. The study is based on Internal Transcribed Spacer (ITS) sequences of fungal symbionts and tRNA(L) (eu) (UAA) intron sequences of cyanobacterial symbionts. The full data set includes 271 Nephroma and 358 Nostoc sequences, with over 150 sequence pairs known to originate from the same lichen thalli. Our results show that all bipartite Nephroma species associate with one group of Nostoc different from Nostoc typically found in tripartite Nephroma species. This conserved association appears to have been inherited from the common ancestor of all extant species. While specific associations between some symbiont genotypes can be observed over vast distances, both symbionts tend to show genetic differentiation over wide geographic scales. Most bipartite Nephroma species share their Nostoc symbionts with one or more other fungal taxa, and no fungal species associates solely with a single Nostoc genotype, supporting the concept of functional lichen guilds. Symbiont selectivity patterns within these lichens are best described as a geographic mosaic, with higher selectivity locally than globally. This may reflect specific habitat preferences of particular symbiont combinations, but also the influence of founder effects.