Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Large differences in carbohydrate degradation and transport potential among lichen fungal symbionts

Abstract

Lichen symbioses are thought to be stabilized by the transfer of fixed carbon from a photosynthesizing symbiont to a fungus. In other fungal symbioses, carbohydrate subsidies correlate with reductions in plant cell wall-degrading enzymes, but whether this is true of lichen fungal symbionts (LFSs) is unknown. Here, we predict genes encoding carbohydrate-active enzymes (CAZymes) and sugar transporters in 46 genomes from the Lecanoromycetes, the largest extant clade of LFSs. All LFSs possess a robust CAZyme arsenal including enzymes acting on cellulose and hemicellulose, confirmed by experimental assays. However, the number of genes and predicted functions of CAZymes vary widely, with some fungal symbionts possessing arsenals on par with well-known saprotrophic fungi. These results suggest that stable fungal association with a phototroph does not in itself result in fungal CAZyme loss, and lends support to long-standing hypotheses that some lichens may augment fixed CO2 with carbon from external sources.

Introduction

Stable fungal associations with algae and/or cyanobacteria, usually referred to as lichens, feature prominently in the history of the discovery and study of symbiosis. In describing the pairing of fungi with microbial photosymbionts for the first time, the Swiss botanist Simon Schwendener proposed that lichen fungal symbionts derive nutrition from “assimilates” of their photosynthesizing partners1. Almost a hundred years later, Smith and colleagues revealed these transferred photosynthates to be polyols and glucose in the case of algae and cyanobacteria, respectively2. They and others traced the transfer of algal fixed carbon into fungal cells, where they found it to be converted into mannitol and arabitol3,4. The fungal-photosymbiont relationship is widely interpreted as conferring net independence from external carbohydrates on the resulting lichen thallus. Accordingly, lichen fungal symbionts have been classified as biotrophs5, and the symbiotic outcome, the lichen thallus, as a “photosynthetic carbon autotroph”6 or “composite autotroph”7.

Fungi are assimilative heterotrophs and thus require a robust machinery of enzymes for scavenging and transporting extracellular nutrients, including carbohydrates. In arbuscular mycorrhizal and ectomycorrhizal fungi, the stable supply of glucose from plants is thought to have led to erosion or loss in many families of carbohydrate-active enzymes (CAZymes)8,9, reflecting a common pattern of compensated trait loss in symbioses10. So what happened to CAZymes in lichen fungal symbionts? Multiple lines of indirect evidence have emerged over the last 40 years to suggest the retention of diverse CAZymes in lichen fungal symbionts (LFSs) or their secondary evolutionary derivates, especially specific plant cell wall-degrading enzymes (PCWDEs). First, molecular phylogenetic studies have shown multiple independent origins of putative saprotroph lineages from lichen fungal symbiont ancestors, both ancient11 and recent12,13,14. How these newly evolved lineages acquired the carbohydrate breakdown arsenal they would presumably need for life without an alga has not been explained. Second, some fungi near the symbiont-to-saprotroph transition have been shown to switch between the two lifestyles facultatively, so-called “optional lichens”15,16; in these cases, the fungus appears not to be obligately dependent on the alga for nutrition. Third, many lichen symbioses exhibit anomalous “substrate specificity”, i.e., they are restricted to specific organic substrates and unable to colonize others, suggesting lack of nutritional autonomy17,18,19. Fourth, lichen fungi are capable of growing axenically in vitro on a variety of sugars other than sugar alcohols, including crystalline cellulose and sucrose (reviewed by Fahselt20). Finally, enzymes involved in breakdown of lichen-exogenous polymers, including cellulose and lignin, have been isolated from lichens in nature (reviewed by Beckett et al.21). These phenomena could be dismissed as exceptions, but their distribution across the fungal symbiont tree hints at deeper underlying fungal capabilities, which if combined with phototroph symbiosis could lead to a kind of “hybrid” lifestyle, in which carbohydrates are obtained from multiple sources. The possibility of multiple carbon sources for lichens was even suggested by Schwendener himself in his original 1869 paper, in which he predicted that two tracks of nutrient acquisition would ultimately be proven: one for lichens that have minimal substrate contact, which he predicted to depend mostly on algal assimilates; and one for lichens that closely hug organic substrates such as tree bark or wood1.

Unlike for many fungi, phenotypic carbohydrate use profiles have seldom, if ever, been developed for LFSs. This is in part due to the recalcitrance of most LFSs to culturing and their extreme slow growth, if culturing is successful. Knowledge gaps around unculturable or slow-growing fungi are common, but have been offset in recent years by genome sequencing. Coupled with the development of widely available databases such as CAZy22, it has become possible to infer CAZymes for species for which a genome, but no experimental evidence, is currently available. Comparative genomic overviews of CAZyme repertoires are now available for many symbiotic fungi9, but no survey exists of comparative CAZyme arsenals in LFSs.

Given the common assumption that lichen symbiont complementarity confers collective autotrophy on the symbiosis, and past results inferring they evolved from saprotrophic or non-lichen biotrophic ancestors (reviewed by Spribille et al.23), we hypothesized that LFSs would exhibit functional losses in CAZymes coinciding with the beginning of stable association with phototroph symbionts, similar to what has been found for PCWDEs in arbuscular- and ectomycorrhizal fungi9,24,25. To test this, however, we would need to map CAZymes across a well-sampled phylogeny and reconstruct the ancestral states of common ancestors, many of which are considerably older than the reconstructed origins of e.g., ectomycorrhizal fungi26. Here we map the occurrence of genes encoding CAZymes at two levels: across representative species of the Leotiomyceta group of Ascomycota, including the origin of the Lecanoromycetes, the largest extant lineage of LFSs; and among major groups within Lecanoromycetes, including species representing different ecological substrate specificities (specialists and generalists) as well as major morphological outcomes of the lichens they occur in (crusts, macrolichens). Our survey of 46 lecanoromycete genomes, 29 of which we sequenced for this project, reveals a complex pattern of retention and loss that lends support to Schwendener’s hypothesis of hidden saprotrophy in some lichens and is not unequivocally consistent with CAZyme erosion upon acquisition of phototroph symbionts.

Results

Data set and phylogenomic reconstruction

We assembled a data set of 83 fungal genomes, including 46 from the class Lecanoromycetes (Supplementary Table 1). Our sampling of Ascomycota genomes outside of Lecanoromycetes was informed by two considerations: 1) selected genomes should be representative of a range of CAZyme repertoires, already mapped by Miyauchi et al.9, amongst others; and 2) they should draw primarily from Leotiomyceta, the group that includes the sibling classes of Lecanoromycetes. Because the few published lecanoromycete genomes currently over-represent fungal symbionts of generalist lichens in one taxonomic group within the subclass Lecanoromycetidae—the order Lecanorales—we generated 29 new genomes for this study, with emphasis on capturing both a diversity of substratum specializations as well as diverse lineages of the subclass Ostropomycetidae (sampling explained in greater detail in Methods). Eighteen of the genomes were obtained as metagenome-assembled genomes (MAGs; Supplementary Table 2). Completeness and quality metrics were comparable for genomes derived from culture and MAGs (Supplementary Figs. 1 and 2). Phylogenomic analysis based on 1310 inferred universally present single-copy orthologs (Fig. 1; Supplementary Fig. 3 and 4 and Supplementary Table 7) recovered major clades and sister group relationships found in recent studies, both among class-level clades of sampled Ascomycota27 as well as within the Lecanoromycetes28,29. For each genome, we performed ab initio gene predictions and obtained functional annotations (CAZymes, InterPro IDs, Pfams). To these we assigned activity on the common plant cell wall substrates cellulose, hemicellulose, lignin and pectin following Miyauchi et al.9.

Fig. 1: Distribution and ancestral states of CAZymes and selected sugar transporters across the evolution of Lecanoromycetes and related classes of Ascomycota projected onto a maximum likelihood phylogenomic tree based on 1310 loci.
figure 1

Symbols beside tree tips refer to life history traits and phototrophic partners of LFS under study. Heatmaps with shades of red indicate the number of genes in different CAZyme classes or involved in degrading complex PCW components. Columns from left to right: AA Auxiliary Activities, CBM Carbohydrate binding module, CE Carbohydrate Esterases, GH Glycoside Hydrolases, GT Glycosyl Transferases, PL Polysaccharide Lyases. cell - Number of genes in 35 CAZyme families involved in cellulose and hemicellulose breakdown. pec - Number of genes in 11 CAZyme families involved in pectin breakdown. lign - Number of genes in three CAZyme families involved in lignin modification. Selection of CAZyme sets follows9, 25. PODs: numbers of heme haloperoxidase and DyP peroxidases potentially involved in lignin modification from Redoxibase (see text). Sugar transporters - selected PF00083 transporters. celldex - Number of cellodextrin transporters. malt - number of maltose transporters. total - Total number of other PF00083 transporters. Below the heatmap are ancestral sizes of CAZyme families involved in (hemi-)cellulose, pectin and lignin modification. Colored circles on tree branches indicate significantly expanded CAZyme families. The size of the circles indicates the number of individual CAFE runs (out of 20 total runs) in which a family was found to be significantly expanded. Expanded gene families recovered in less than five runs are omitted here. Exact numbers are given in Supplementary Fig. 12. Source data are provided as a Source Data file.

CAZymes

Many CAZyme families are shared across all sampled genomes. Glycosyl transferases (GTs), which are involved in glycosylation and the synthesis of polysaccharides, differ little across all analyzed genomes and do not exhibit any significant reduction in Lecanoromycetes, suggesting that a core synthetic machinery remains largely unchanged across the evolution of the sampled Ascomycota (Fig. 1). While significant within-group variation exists, the mean number of glycoside hydrolase (GH) genes in lecanoromycete genomes is 40.7% lower compared to other sampled genomes (p = 0.0000; Supplementary Tables 9, 10 and 11). CAZymes with Auxiliary Activities (AA), Carbohydrate Binding Modules (CBM) and Carbohydrate Esterases (CE) are also reduced significantly in Lecanoromycetes (Supplementary Table 10). The bulk of these differences can be attributed to a small number of CAZyme families (Supplementary Fig. 11), most of which are plant cell wall degrading enzymes (PCWDEs). Additionally three CAZyme families widespread in the sampled Ascomycota were not detected in Lecanoromycetes. These include PL4, which contains pectin degrading rhamnogalacturonan endolyases; CBM67, which binds to L-rhamnose and frequently occurs in multi-domain protein with enzymes in GH78 and PL1; and AA13, which contains lytic polysaccharide monooxygenases (LPMOs) involved in starch breakdown. In contrast to these reductions, numerous CAZyme families are not reduced at all in Lecanoromycetes compared to other sampled Ascomycota. Indeed some, including those involved in degradation of endogenous fungal cell wall polysaccharides, such as GH128 and AA5, are even expanded in Lecanoromycetes (Supplementary Fig. 11).

All lecanoromycete genomes possessed genes encoding enzymes predicted to act on plant cell wall compounds, including cellulose and hemicellulose (e.g., GH5 and GH43) and lignin (AA1, AA2 and AA5). Depending on the symbiont configuration, these organic polymers may be produced by the phototroph partner and/or be exogenous to the lichen symbiosis (see Discussion). Principal components analyses (PCA) of CAZyme family numbers in the sampled genomes reveals differences in the amount of variation in CAZyme composition among lecanoromycete and other sampled ascomycete genomes (phylogenetically corrected PCA: Fig. 2; regular PCA: Supplementary Fig. 10). For CAZymes predicted to act on cellulose, hemicellulose and pectin, the number of predicted gene families varies less among lecanoromycete genomes than among comparable classes of Ascomycota (Fig. 2A, C), reflected in tight clustering in Lecanoromycetes versus wide scattering in other ascomycete classes. For lignin-modifying enzymes the variation is similar (Fig. 2B). Within Lecanoromycetes, however, the greatest amount of variation in predicted gene sets is exhibited in the subclass Ostropomycetidae (represented in our sample by 24 genomes), specifically in cellulose/hemicellulose (Fig. 2A) and pectin degradation (Fig. 2C), resulting in scattering in the PCA ordination. By contrast, genomes from the subclasses Acarosporomycetidae, Umbilicariomycetidae and Lecanoromycetidae form a cluster of similar predicted CAZyme sets, resulting in tight clusters in the PCA ordination. These latter three subclasses are represented in our dataset by one, four and 17 genomes, respectively; for convenience results reported for Lecanoromycetidae will refer by extension to results from all three of these subclasses unless otherwise specified.

Fig. 2: Similarity of CAZyme sets involved in the breakdown of different complex plant-based polysaccharides based on phylogenetically corrected Principal Components Analysis.
figure 2

Different colors indicate taxonomic groups; genomes from the subclasses Acarosporomycetidae and Umbilicariomycetidae use the same color codes as Lecanoromycetidae for simplicity. A Similarity of CAZyme families involved in cellulose and hemicellulose breakdown. B Similarity of CAZyme families involved in lignin breakdown, C Similarity of CAZyme families involved in pectin breakdown. Displayed below are representative members of the two subclasses Ostropomycetidae (with light green border) and Lecanoromycetidae (dark green border). D Graphis scripta, E Icmadophila ericetorum, F Mycoblastus sanguinarius, G Peltigera leucophlebia, H Xylographa carneopallida, I Agyrium rufum, J Evernia prunastri, K Cladonia macilenta. Image credits: Agyrium rufum: Paul Cannon (fungi.myspecies.info); Creative Commons: BY-NC 4.0. Peltigera leucophlebia: Jason Hollinger, uploaded by Amada44, CC-BY 2.0, https://commons.wikimedia.org/w/index.php?curid=24213606. Evernia prunastri: by Jason Hollinger, CC-BY 2.0, https://commons.wikimedia.org/w/index.php?curid=50595319. Cladonia macilenta: Bruce McCune & Sunia Yang - Lichen, CC-BY 4.0-NC, https://lichens.twinferntech.net/pnw/species/Cladonia_macilenta.shtml; other images by the authors. Source data are provided as a Source Data file.

Genomes from the Ostropomycetidae possess higher numbers of predicted CAZyme genes that act on cellulose and hemicellulose than those of Lecanoromycetidae, which are reflected in significantly higher numbers of GHs and CEs (p = 0.0000 and p = 0.0024 respectively; Supplementary Table 10). The differences are driven by one lineage in particular, represented by five genomes from the orders Ostropales and Gyalectales and referred to here as the OG clade (Fig. 1), and are much less if the two suborders are compared without the OG clade genomes. The OG clade contains numbers of predicted CAZyme genes for cellulose and hemicellulose modification that are over fourfold more numerous than the lowest lecanoromycete numbers, which are in Cetradonia and Peltigera, and equal or exceed the numbers in well-studied eurotiomycete saprotrophs such as Penicillium. These disparities are accounted for in large part by gene assignments to two CAZyme families, GH5 and GH43. Because both of these are large, heterogeneous families that include multifunctional CAZymes, we mapped putative GH5 and GH43 orthologs from the analyzed genomes against sequences of experimentally validated enzymes and predicted which are secreted (Fig. 3). Extensive gene duplication in Ostropomycetidae is found in gene sequences close to characterized cellulases (GH5 subfamily 5; EC 3.2.1.4) and from both lecanoromycete subclasses in sequences close to characterized 1,3-beta-glucosidases (GH5 subfamily 9; EC 3.2.1.58) and endo-1,4-beta-mannosidases (GH5 subfamily 7; EC 3.2.1.78). Notably, many putative lecanoromycete CAZyme genes from GH5 do not closely cluster with any characterized sequences and form their own clades with sequences from other sampled classes of Ascomycota (Fig. 3A).

Fig. 3: Gene trees of two CAZyme families involved in cellulose (GH5) and hemicellulose (GH43) breakdown. Each tree includes all experimentally characterized sequences combined with all sequences from the 83 genomes studied here.
figure 3

CAZyme Subfamilies are labeled with numbers and gray rectangles. Three columns along tree tips display additional information of corresponding sequences when available. EC - Enzyme Code of experimentally characterized sequences downloaded from cazy.org. Sub Loc - Predicted subcellular location of Enzyme with DeepLoc. Tax - Taxonomic assignment of organisms from which the sequence comes from. For larger subfamilies, functions of characterized sequences based in Enzyme Code numbers are given as colored squares. Sequences used for heterologous expression experiments are marked with red arrows in GH5 subfamily 5. Source data are provided as a Source Data file.

The subclass Ostropomycetidae also possesses a larger proportional representation of genes coding for enzymes predicted to be involved in hemicellulose, specifically xylan, breakdown. Families GH6, GH7, GH11 and GH62, GH67 and GH131 were absent in Lecanoromycetidae and are only present in Ostropomycetidae in the OG clade and the saprotroph Agyrium (Supplementary Fig. 11). Hemicellulose breakdown also involves several CEs of which CE2 is present in the majority of Ostropomycetidae genomes and absent in Lecanoromycetidae (Supplementary Fig. 11). CE2 contains acetylxylan esterases (AXEs;30); acetylxylan is a major component of hemicellulose. An analysis of predicted GH43 gene sequences and their subcellular location compared to those of characterized enzymes shows that most characterized GH43s are in Ostropomycetidae (subfamilies 1, 6, 24, 26 and 36), but as with GH5s, the majority of sequences are not close to any characterized CAZymes (Fig. 3B).

Numerous enzymes can potentially be involved in lignin modification17. In Lecanoromycetes, CAZymes associated with lignin modification are predicted in all genomes, but again the highest numbers are found in Ostropomycetidae (Fig. 1). The numbers of predicted genes in Graphis, Thelotrema and Varicellaria are among the highest in the analyzed genomes. Most of these genes are members of CAZyme family AA1 (laccase, ferroxidase and multicopper oxidases, EC 1.10.3.x), and several are phylogenetically close to experimentally characterized enzymes from EC 1.10.3.2, which included laccase, p-diphenol:oxygen oxidoreductase and ferroxidase (Supplementary Fig. 8). We also found smaller numbers of predicted AA2 which contain Class II peroxidases (EC 1.11.1.x; Supplementary Fig. 9) in Lecanoromycetes, but none of these are phylogenetically close to any experimentally characterized peroxidases (Supplementary Fig. 15). As not all commonly studied enzymes involved in lignin modification are included in the CAZy database, we also mapped ortholog groups for all peroxidases in RedOxibase (31, see Supplementary Fig. 6). Two ascomycete peroxidase groups that have been shown to have activity on lignin molecules, heme haloperoxidases and dye-decolouring peroxidases (DyP), exhibit large differences in predicted gene numbers in Lecanoromycetes (Fig. 1). Heme haloperoxidases occur in three ortholog groups, of which one was recovered in all sampled Lecanoromycetes but with greatly elevated numbers in the OG clade, as well as in Pseudevernia and Mycoblastus (Supplementary Fig. 6). A second heme haloperoxidase ortholog group is found in several Ostropomycetidae, while a third was not recovered in Lecanoromycetes. DyPs likewise occurred in three ortholog groups, one of which was recovered in most Ostropomycetidae and in about half of Lecanoromycetidae, and another of which was recovered only in the wood-dwelling genus Xylographa (Supplementary Fig. 6). The third DyP ortholog group was not detected in any lecanoromycete genome.

The largest differences in predicted PCWDEs among lecanoromycete genomes come from CAZymes acting on pectins. All but three of the Ostropomycetidae genomes possess a few predicted CAZymes involved in pectin degradation (Fig. 1). These include GHs (GH28, 49, 53, GH79, GH108), CEs (CE8, CE12, CE15) and the only polysaccharide lyases (PL; from PL1 and PL3, containing pectate lyases, EC 4.2.2.2) predicted in Lecanoromycetes (Supplementary Fig. 8). Lecanoromycetidae, by contrast, almost completely lack gene predictions for pectin degradative enzymes and completely lack PLs. Similar to CAZyme predictions for cellulose, hemicellulose degradation and lignin modification, the largest numbers of predicted genes involved in pectin degradation come from the OG clade of Ostropomycetidae. In these genomes, predicted gene numbers for CAZymes involved in pectin degradation even exceed those of known model saprotrophs in Eurotiomycetes, though the difference is not significant.

In addition to the PCWDEs outlined above, all sampled Lecanoromycetes possess a gene assigned to GH32, a family which includes invertases involved in the conversion of sucrose to glucose or fructose. Invertases are generally interpreted as an indicator of use of apoplastic plant sucrose. However, lecanoromycete GH32s do not cluster closely with any characterized invertases (Supplementary Fig. 13).

Predictions of gene family contraction and expansion

In order to establish whether CAZyme patterns in Lecanoromycetes are due to gene gain or loss, we reconstructed ancestral gene numbers for each of the three main groups of PCWDEs (cellulose/hemicellulose, pectin, lignin). All three PCWDE groups are reduced already at the split between Eurotiomycetes and Lecanoromycetes (ELS node; Fig. 1), especially enzymes involved in (hemi-)cellulose and pectin breakdown. Much of the signal for this pattern comes from GH5 and GH43 which contain many well-characterized cellulases and hemicellulases, respectively (Fig. 3; Supplementary Fig. 7). Another apparent shedding of CAZymes, especially of those involved in pectin degradation, occurs concomitantly with the origin of Lecanoromycetes, the largest extant clade of LFSs (ALSS and ALSL nodes; Fig. 1). Here, too, the signal comes mainly from GHs. In several cases our ancestral state reconstruction indicates complete losses, despite the prediction of several of the relevant families in Ostropomycetidae genomes (Supplementary Fig. 8).

Gene family expansion analyses of all CAZymes revealed six significantly expanded families along different branches in the phylogeny. The only GHs expanded within Lecanoromycetes are GH43, which are expanded just after the last common ancestor of Ostropomycetidae (Fig. 1; Supplementary Table 8); another, GH18, was expanded in four model runs at the base of Sordariomycetes (Supplementary Fig. 12). Other significantly expanded CAZyme families in Lecanoromycetes include Auxiliary Activity (AA) families AA3, AA7 and AA9, which are variously expanded in Xylographa, the OG clade and Cladonia. Auxiliary Activity families typically act in concert with other CAZymes and the three expanded AA families all contain genes involved in cellulose breakdown. The only other expanded CAZyme family was CE10, now considered to be a family of esterases acting on non-carbohydrate substrates, which was expanded at the last common ancestor of Cladonia species and Umbilicariomycetidae.

Transporters

As an additional line of evidence for potential use of exogenous carbohydrates, we mapped the numbers of predicted sugar transporters across all genomes (Fig. 1). Most predicted transporters were more or less evenly represented across most genomes, but two groups of transporters exhibited distinct presence/absence patterns within Lecanoromycetes. Cellodextrin transporters, close orthologs of experimentally demonstrated cellodextrin transporters from Aspergillus and Penicillium (Supplementary Fig. 5), were predicted in all Ostropomycetidae genomes except Schaereria, and were most numerous in Lignoscripta, Stictis, Ptychographa and Xylographa, all of which are LFSs of lichens with high wood or bark specificity (Fig. 1). Cellodextrin transporters are involved in transmembrane import of cellobiose and other cellodextrins, which are short beta-linked glucose fragments of cellulose. Of the remaining 22 lecanoromycete genomes, all four Umbilicariomycetidae and the one Acarosporomycetidae genome had cellodextrin transporters, while only five of the 17 Lecanoromycetidae genomes did; in all of these, gene numbers were well below the average observed in Ostropomycetidae. We detected a similar pattern for the predicted occurrence of maltose transporters (Fig. 1), involved in transmembrane import of alpha-linked starch breakdown products: 22 of 24 Ostropomycetidae genomes had predicted maltose transporters, but only five lecanoromycete genomes outside of Ostropomycetidae had any, three of them in Umbilicariomycetidae.

Evidence of cellulase functionality

To validate the functionality of putative cellulases found in lichens, we selected X. bjoerkii genes MMC18_000518 (MCJ1387675.1 and MMC18_004565 (MCJ1391700.1) after aligning sequences from all Xylographa species recovered in GH5 Subfamily 5 (Fig. 4d), hereafter called cellulase A and B, for further analysis based on their similarity with Cel5A cellulase sequences from Trichoderma reesei, for which there has been substantial structural and functional characterization. The cellulase domain of each gene was expressed in Escherichia coli as a C-terminal fusion with the maltose binding protein, which enhances protein expression and stability and provides a convenient handle for purification. The partially purified proteins were active on both cellulose (β-1,4-linked glucose) and barley β-glucan (alternating β-1,4/β-1,3-linked glucose) but not xylan (β-1,4-linked xylose containing side branches of α-arabinofuranose and α-glucuronic acids) (Fig. 4a–c). Both exhibited a pH optimum of 5, similar to other Cel5A enzymes but differ slightly in their temperature optima. While both are active at 20 °C, only cellulase A is also active at higher temperatures. Closely related orthologs of cellulase A and cellulase B were recovered in all lecanoromycete genomes. Lecanoromycetidae mostly only had a single ortholog with the exception of Lobaria, Pseudocyphellaria, Pseudevernia and Sticta which each have two orthologs (Fig. 4d). In contrast, the majority of Ostropomycetidae genomes have two or more orthologs, and the highest numbers can be found in Xylographa species and Stictis which each have five (Fig. 4d).

Fig. 4: Enzymatic activity of two putative cellulases A and B from Xylographa bjoerkii and orthologs of these genes.
figure 4

Activity of cellulases A and B at different pH conditions (a), at different temperatures (b), and with different substrates (c). Y axis represents mean absorbance (+/− standard error of mean) at 595 nm after adjusting for the dilution factor. Experiments were performed in technical and biological triplicates. d Number of orthologs of cellulase A and B in the 72 genomes that possess them. Different colors indicate taxonomic groups; genomes from the subclasses Acarosporomycetidae and Umbilicariomycetidae use the same color codes as Lecanoromycetidae for simplicity. Source data are provided as a Source Data file.

Discussion

Our analysis of CAZyme and sugar transporter genes paints a picture of a lecanoromycete arsenal that is larger, more diverse and more consistently present than expected. Our data show a reduction in mean numbers of lecanoromycete degradative CAZyme genes relative to other sampled Ascomycota, but also reveal significant differences within Lecanoromycetes. Genes for the breakdown of cellulose, hemicellulose and pectin are disproportionately enriched in the subclass Ostropomycetidae, and some genomes, notably in the OG clade, possess overall CAZyme numbers and functionality similar to those of well-known model saprotrophs such as Aspergillus and Penicillium. Furthermore, LFSs associated with the same genus of algal phototroph, Trebouxia, can retain multiple genes coding for pectin degradation as well as cellodextrin transporters (in Lambiella, Loxospora, Ptychographa and Xylographa), or they can largely lack these genes (in Ramalina). This implies that association with Trebouxia does not in itself result in gene loss, and the exact nature of the fungal-phototroph interaction, and other aspects of symbiosis biology, may need to be reconsidered. Understanding when and where in the fungal mycelium the CAZymes are deployed will be critical to determining their biological significance.

The largest differences in CAZyme gene numbers among the sampled lecanoromycete genomes represent PCWDEs. This raises a central question: what are the targets of LFS PCWDEs? The two most obvious candidates are the phototroph itself, on the one hand, and lichen-exogenous polysaccharides, such as wood and tree bark, or even non-plant polymers, on the other. The first possibility — targeting of the alga — echoes suggestions in studies of ectomycorrhizal fungi, where PCWDEs have been postulated to play a role in “cell softening” or “remodeling” of root tissues in their vascular plant symbionts upon contact initiation32. Cellulases and other CAZymes have been postulated to be involved in haustorial penetration of the algal cell wall33,34 and degradation of algal cell walls in fresh lichen growth tips35. They could also play a role in digestion of dead algal cell walls, especially if algal populations turn over during the lifespan of the thallus, as has long been suspected36. Transcriptomic studies of isolates, co-cultures and natural lichens could provide evidence of this. In two studies that reported CAZyme differential expression in symbiont coculture experiments, one detected upregulation of multifunctional GHs that could also be involved in fungal cell wall modification (GH2 and GH12), but none of the core lecanoromycete cellulases or hemicellulases we report here37, while another reported upregulation of gene models assigned to GH5, GH43 and GH4538 (in their Additional Files 6.3 and 6.4). Most common algal symbionts are thought to contain cellulose in their cell walls but no pectin39; it is unclear if any lichen algal symbiont possesses pectin in its cell walls.

The second possibility, that LFS PCWDEs target lichen-exogenous polysaccharides, is supported by two lines of evidence. The first is experimental. Both cellulases and polygalacturonases, which are active on pectins, have been detected in incubated whole lichens both in the presence40,41 and absence42 of cellulose-containing algal photosymbionts. Cellulase production furthermore has been reported to vary depending on the species of tree it is growing on43. Some lichens have also been experimentally shown to possess an enzymatic and extracellular redox cycling toolbox consistent with lignin modification21,44,45. The second line of evidence relates to the ecological substratum use and relationship to known saprotrophs. The enrichment of PCWDEs we found tracks closely with the secondary origin of non-lichen saprotrophs and, among lichens, those with wood-obligate ecology. Several secondary origins of putative fungal saprotrophs occur within the Lecanoromycetes, some of which include fungi involved in “optional lichens”, in which the fungus can occur either as an LFS or saprotroph46. The most speciose of these saprotroph groups, the Ostropales, arose within our “OG clade” that possesses the largest PCWDE enrichment of all sampled lecanoromycete genomes. We also found cellulase gene family expansions in LFSs of wood-obligate lichens such as Ptychographa and Xylographa, and our experimental evidence for cellulase activity derives from one of these genomes (Xylographa bjoerkii).

Enzymatic targeting of algal cell walls or lichen-exogenous polymers are not mutually exclusive possibilities. Nor do all exogenous substrates necessarily have to be plant cell walls, as lichen symbionts also interact with a wide range of other fungi, notably fungal parasites and bacteria47. Determining where and when these genes are expressed in nature will require consideration of the life cycle and spatial extent of LFS mycelia, both of which are still poorly understood. The mycelia of sexually reproducing LFSs go through an aposymbiotic stage of unknown length, during which it has been suggested they may be saprotrophic48. Lichens at high latitudes or under snowpack may go through seasonal fluctuations in fixed carbon input, which could be augmented by other sources21. Even after symbiosis is established, parts of the mycelium can be free of phototroph cells, and exhibit deviations in carbohydrates that suggest other metabolic processes are in effect than in the phototroph-associated mycelium49. Many lichens include a phototroph-free “prothallus”, in which fungal hyphae radiate beyond the zone of phototroph cells50. Others possess a “hypothallus”, a phototroph-free cushion of mycelium in direct contact with the substrate49. Macrolichens are often anchored onto their substrates by phototroph-free “holdfasts”, “rhizines” or mycelial pegs which have traditionally been interpreted as having an exclusively structural, stabilizing function, but can extend as mycelial networks into xylem51 or living moss mats42. Parts of the mycelium with phototroph cells may also differ based on position relative to growth tips, which are thought to include larger proportions of living cells52, and entire strata can be phototroph-free. We currently lack data on the role CAZymes could play in nutrient scavenging or cell wall remodeling in these phototroph-free sectors of the symbiont mycelium.

In addition to the differences and enrichment patterns we have analyzed here, the lecanoromycete genomes we surveyed are rich in genes that do not cluster with sequences from experimentally characterized CAZymes. Much work remains to be done to characterize their substrates and activity. These include two clades of GH5 cellulases (Fig. 3), numerous sequences assigned to GH43 hemicellulases but without closely clustering orthologs (Fig. 3), putative GH32 invertases (Supplementary Fig. 13), and numerous sequences assigned to AA5 (Supplementary Fig. 17), GH30 (Supplementary Fig. 38), GH72 (Supplementary Fig. 53), and GH79 (Supplementary Fig. 55). The need for more experimental work is especially true of putative lignin-modifying enzymes. Lichen fungi have been shown to possess elements of a lignin modification toolbox including a brown rot-like system of hydroxyl radical production in the presence of a quinone and chelated Fe44,53, as well as a variety of metal-containing oxidoreductases including laccase and heme peroxidases that are induced by starvation and treatment with soluble cellulose and lignin breakdown products21,44. Specifically, experimental work has shown that LsaPOX, a putative heme peroxidase from the fungal symbiont of the lichen Leptogium saturninum, can metabolize the β-O-4 lignin model dimer adlerol into veratraldehyde, a reaction also performed by white-rot fungi45 (a published oligopeptide sequence of LsaPOX is closest to 4206_AniHalPrx05 in Redoxibase, with a close hit in Cladonia grayi, corresponding to the haloperoxidase_.haem. column in Supplementary Fig. 6). Our survey shows that laccase and peroxidase ortholog groups are widespread in Lecanoromycetes and especially enriched in species on wood, with numerous novel sequence groups. However, none of these alone should be interpreted as prima facie evidence of lignin modification in nature, as these enzymes can also be involved in other aspects of lichen biology, such as degradation of secondary metabolites44.

Did stable phototroph association coincide with the loss of PCWDEs in LFSs? Looking only at the reductions in mean CAZyme gene numbers, the answer would appear to be positive, but the occurrence of LFSs in the OG clade (four of which in our sample are LFSs and one of which, Cyanodermella, is an endophyte) — with CAZyme arsenals as large or larger than those of many saprotrophs — shows that LFSs do not necessarily lack CAZyme arsenals. Despite their minority representation in our data set, CAZyme-rich LFSs may in fact be numerous in nature: the Ostropales and Gyalectales which make up the OG clade include over 3200 named LFS species, almost 17% of named LFSs54, in addition to approximately 500 named non-LFS fungi, mostly saprotrophs. For phototroph acquisition to have no evolutionary consequences for the CAZyme arsenal, symbiont-derived photosynthates would likely have to provide a function other than as a substrate for growth and respiration. In fact this has already long been postulated, in the form of a role for polyols as compatible solutes55. Some have even suggested this may be the main role for transferred photosynthates56,57. A further indication that photosynthates are not solely nutritional subsidies could be the retention of invertases in Lecanoromycetes. In ectomycorrhizal fungi, the loss of invertases is thought to limit their ability to access plant sucrose and reinforce their dependence on plant-derived glucose9. Though the orthologs we found did not exhibit high similarity to characterized sequences, the prediction of an invertase is consistent with the ability to culture LFSs on sucrose58 and experimental evidence of invertases in lichen fungi20,59.

A strong inference about whether LFS CAZyme gene reduction coincided with the onset of stable phototroph association requires greater certainty about the ancestral states along the ascomycete backbone, and the lecanoromycete phylogenetic backbone in particular. Our ancestral state reconstruction suggests a gradual loss of GHs and PLs since the last common ancestor of Lecanoromycetes and Eurotiomycetes, and only pectin degradation genes are inferred to have been abruptly reduced. This reconstruction is sensitive to taxon undersampling and may be a conservative estimate. The large CAZyme arsenal of the OG clade, with dozens of unlinked genes and numerous PLs, is unlikely to have been acquired by lateral gene transfer. If the PL-rich OG clade CAZyme arsenal is in fact ancestral, this would imply that CAZyme loss is not an automatic consequence of stable association with phototroph symbionts, but rather of subsequent events or adaptations. If so, this implies that the OG clade CAZyme arsenal was lost no fewer than seven times in our tree. Support for a scenario of frequent mass gene reduction comes from the fact that one such loss appears to have happened within the OG clade itself: Gomphillus exhibits marked gene number reductions compared to its closest sampled relatives. Our data do not currently allow the hypothesis that the OG clade CAZyme arsenal is ancestral to be rejected. If it is not, it becomes more likely that CAZyme loss in LFSs is driven by additional processes.

Lecanoromycetes, the focus of our study, are the largest extant group of LFSs. A fuller picture of CAZyme arsenal retention during lichen symbioses and especially in pairing with specific phototroph partners will emerge as more genomic data become available from other fungal classes. Two additional LFS genomes were included in the outgroups of our dataset, Endocarpon (Eurotiomycetes) and Arthonia (Arthoniomycetes). Both genomes echo the lack of predicted cellobiose transporters found within numerous Lecanoromycetes, but differ in their predicted pectin degradation capacity, with Arthonia possessing gene numbers similar to those found in the Ostropomycetidae. Arthonia also possesses higher numbers of non-CAZyme peroxidases, similar to the OG clade. If the OG clade CAZyme arsenals are related to saprotrophy, similar CAZyme arsenals could be expected in lichen fungal symbionts e.g., from the Dothideomycetes, which have long been postulated to have dual lifestyles as mutualists and saprotrophs60. Within Lecanoromycetes, the most striking apparent functional losses, both in terms of CAZyme genes as well as in cellodextrin and cellobiose transporters, occurred outside of the Ostropomycetidae, especially in the subclass Lecanoromycetidae, which largely lacks these two types of transporters. LFSs in Lecanoromycetidae, like those in Ostropomycetidae, are considered obligate symbionts of algae or cyanobacteria. They share many of the same algal symbionts, especially Trebouxia, and have no known lifecycle differences. However, they differ in general thallus architecture. Ostropomycetidae almost exclusively form crustose thalli in which thallus-to-substrate surface area contact is maximized. Lecanoromycetidae, and to some extent Acarosporomycetidae and Umbilicariomycetidae, include LFSs involved both in crustose thalli as well as so-called macrolichens, in which the thallus often becomes greatly expanded into leaf-, hair- or shrub-like forms, in which surface contact is minimized. Lecanoromycetidae macrolichens include some of the largest lichens by biomass, and in theory these would require more carbon, not less. So what additional adaptational processes could drive CAZyme loss?

Two basic biological traits surrounding the crust-to-macrolichen transition deserve more attention with respect to their potential effect on sugar uptake. In the first, the type of fungal-algal contact differs, from intracellular haustoria in many crust lichens to so-called intraparietal haustoria, which do not breach the cell wall, in some crusts and all macrolichens61. The exact consequence of this is unknown, but intracellular haustoria are characteristic of pathogenic fungi and may require a greater variety of CAZymes to penetrate the algal cell wall. In the second, macrolichens owe their architecture to a well-developed exopolysaccharide gel, termed cortex, which forms a rigid structural scaffold considered a prerequisite to macrolichen formation62. The polysaccharide composition of this layer varies widely across lichen symbioses63 and is poorly characterized and mapped across the phylogenetic tree. Additional evolution in cortex layering happened in many lichens involving Lecanoromycetidae64, though in how many remains unclear. The cortex, which is hydrophilic, also mediates water retention63 and has been shown to operate like a sponge for passive uptake of dissolved nutrients65 and glucose66. One of the largest epiphytic macrolichens has been experimentally demonstrated to take up tree-derived glucose67. Whether this is specifically facilitated by the cortex remains to be tested, but as the cortex mediates the passage of environmental molecules for many macrolichens, it seems likely. If capture of simple sugars is found to be a general function of the cortex in different environments, it could be expected to have significant evolutionary consequences for maintenance of CAZyme genes used for more costly carbohydrates.

The existence of a robust carbohydrate breakdown machinery across a large swathe of LFS evolution calls into question the assumption, underlying decades of ecophysiological work on lichens, that the lichen carbon economy is exclusively the sum of algal CO2 fixation. The finding that some LFSs have CAZyme arsenals on par with saprotrophs lends support to Schwendener’s hypothesis that different types of lichens exist: those with fungal symbionts that depend mostly on their phototroph, and those with fungal symbionts that augment their carbon assimilation from external sources1. The longstanding assumption that LFSs solely utilize fixed carbon must now be weighed against competing hypotheses, including A) that some or many LFSs build their mycelia from non-algal carbon sources, including absorbed monosaccharides and complex polysaccharides, before or during symbiosis; and B) that many different models of carbon acquisition—fixed, seasonal, facultative, scavenged and/or absorbed—may exist under the umbrella of what we currently call lichen symbiosis.

Methods

Extended Materials and Methods as well as Supplementary Figures and Tables are available as Supplementary Information. All used software is listed in Supplementary Table 3.

Used genomes

We built a data set consisting of 83 fungal genomes from the phylum Ascomycota, including 46 genomes from the Lecanoromycetes, twenty-nine of which were newly generated for this study, and 37 genomes from the related classes Eurotiomycetes, Dothideomycetes, Arthoniomycetes and Sordariomycetes. Within Lecanoromycetes, the acquisition of new genomes was targeted to include a representation of lineages including 1) LFSs of macrolichens (Cladonia, Evernia, Lobaria, Peltigera, Pseudevernia, Pseudocyphellaria, Ramalina, Sticta, Umbilicaria); 2) LFS lineages involved in crust-forming lichens for which carbon acquisition from the substrate has been postulated, including wood specialists (Lecidea scabridula, Lignoscripta, Ptychographa, Puttea, Xylographa18), bark specialists (Graphis, Loxospora, Schaereria, Stictis, Thelotrema, Varicellaria), specialists of decaying plant matter (Gomphillus, Icmadophila), mineral soil (Dibaeis), rock (Acarospora, Trapelia) and a lichen-on-lichen “parasite” (Lambiella); 3) LFS lineages from crust-forming lichens that behave as ecological generalists (Gyalolechia, Mycoblastus, Toensbergia); and 4) a lecanoromycete saprotroph (Agyrium) and endophyte (Cyanodermella68). Our sampling simultaneously represents a cross-section of lecanoromycete evolution, with 24 genomes from the subclass Ostropomycetidae, 17 from the subclass Lecanoromycetidae, and four and one each from the species-poor subclasses Acarosporomycetidae and Umbilicariomycetidae, respectively. Ten genomes were derived from cultured samples and 18 are metagenome-assembled genomes (MAGs) assembled and binned according to Tagirdzhanova et al.69. Culture-derived genomes differed little from MAGs in estimated completeness and gene numbers (Supplementary Figs. 1 and 2).

Genome assembly and filtering

Raw sequence data was inspected with FastQC 0.11.7 and trimmed with trimmomatic 0.29 to remove adapter remnants and low-quality reads. Trimming parameters are given in Supplementary Table 4. Genomes were assembled using SPAdes 3.12.0 or Abyss 2.0.2. Assembly parameters are provided in Supplementary Table 5. To extract LFS contigs, draft assemblies were filtered using blobtools 1.1.1 or CONCOCT 1.2 and assembly completeness assessed with Quast 4.6.3 and BUSCO 3.0.2. We additionally extracted mitochondrial contigs by blasting mitochondrial genes downloaded from NCBI against each de novo assembly. We then discarded contigs with blast hits of e < 1e-03 and alignment length >500 bp.

Gene calling and functional annotation

We used funannotate 1.8.3 to perform gene-calling and functional annotation for all used genomes in the same way. This reduces potential biases introduced by different gene-calling and annotation methods. De novo sequenced genomes were repeat-masked using RepeatModeler and RepeatMasker. We then used Augustus 3.3, snap 2013_11_29, GeneMark-ES 4.62 and GlimmerHMM 3.0.4 and tRNA-Scan 2.0.5. Functional annotations for all predicted protein sequences were inferred using Interproscan-5.48–83.0, HMMer3 searches against dbCAN (v9) and Pfam (v33.1) as well as eggnog-scanner searches against EggNOG (v4.5.1) databases.

Phylogenomics

Phylogenomic trees were calculated using the phylociraptor pipeline (https://github.com/reslp/phylociraptor). We ran BUSCO 3.0.2 on each genome assembly to identify single copy-orthologous, combined amino-acid sequences of each BUSCO gene from all genomes. Only genes which were present in >50% of genomes were aligned using mafft 7.464 and trimmed using trimal 1.4.1. We calculated single-gene trees using iqtree 2.0.7 and a species tree using ASTRAL 5.7.1. We created a concatenated alignment from all alignments, estimated the best substitution model for each and calculated a tree based on a partitioned analysis of the concatenated alignment using IQ-Tree 2.0.7. We used the concatenated phylogeny to generate an ultrametric tree using r8s 1.81. This tree was used for subsequent analyses. We used custom python and R scripts to plot phylogenomic trees.

Selection of CAZyme groups

We selected sets of CAZyme families involved in (hemi-)cellulose, pectin and lignin degradation based on previous studies9,25 (Supplementary Table 6). For the lignin set we additionally identified class II peroxidases based on similarity of all Ascomycota class II peroxidases downloaded from RedOxiBase (30; accessed Jul. 14, 2021; http://peroxibase.toulouse.inra.fr/) using Orthofinder 2.5.2.

Ancestral state reconstruction of CAZyme families and similarity of CAZyme sets

We reconstructed the ancestral size of each CAZyme family using CAZyme counts from our genome annotations and our ultrametric phylogenomic tree in R. We used anc.ML from phytools (v0.7–70) under an Ornstein-Uhlenbeck model of trait evolution. We then visualize ancestral states with custom R scripts.

We calculated phylogenetically corrected PCAs in R as implemented in the phytools function phyl.pca using maximum-likelihood optimization and regular PCAs using the prcomp function (Supplementary Fig. 10). We used log-transformed gene-counts of CAZyme families and our ultrametric phylogeny as input for all PCAs. The PCAs were visualized using custom R scripts.

Gene family expansion analysis

We analyzed gene family expansions of CAZyme families, based on CAZyme family count data from funannotate annotations and our ultrametric phylogenomic tree using CAFE 5 (git commit 08d27a1). We ran CAFE twenty times while accounting for different gene-birth rates and error-models and summarized the results using custom R and python scripts.

Additional characterization of CAZymes

CAZyme families known to contain important degradative enzymes (e.g., cellulases, hemicellulases) were further characterized using Saccharis v1 (git commit 9a748be; Supplementary Figs. 1461). For each CAZyme family, we downloaded all characterized genes from cazy.org (accessed: March 3, 2021) including additional information such as taxonomy, accession numbers and CAZyme (sub-)family assignments. Using Saccharis, we created MUSCLE 3.8.31 alignments for each gene family of all (from genomes used in this study plus characterized genes from cazy.org) genes and created maximum-likelihood trees using Fasttree 2.1.10. For all genes included in our phylogenetic reconstructions, we also predicted subcellular locations using DeepLoc 1.0. We only considered predictions for subcellular locations if the predicted probability was >70%. We used custom R and python scripts to visualize the trees.

Heterologous expression of putative cellulases

Sequences assigned to GH5 subfamily 5, with confirmed cellulolytic activity, from obligately wood inhabiting LFS Xylographa species were aligned to the characterized and crystalized cellulase domain of Trichoderma reesei (PDB: 3QR3; https://www.rcsb.org/structure/3QR3). Based on sequence similarity, we selected two candidate cellulase sequences from Xylographa bjoerkii for testing cellulolytic activity.

Enzyme activity of the two candidate cellulases A and B was tested by combining in a 1.5 mL microcentrifuge tube 50 µL enzyme; 100 µL buffer A at either pH 4, 5, or 6; and 50 µL either AZCL-HE-cellulose, AZCL-β-glucan, or AZCL-xylan. Tubes were incubated at 4, 20, 37, 50, or 60 °C for 48 h. To measure activity, samples were centrifuged at 13,000 rpm for 5 min to settle any debris, then 100 µL supernatant was removed to a 96-well flat-bottom plate. Absorbance at 595 nm was measured by plate reader and blanked with a sample containing water instead of enzyme.

Identification of sugar- and sugar-alcohol transporters

To identify putative sugar- and sugar-alcohol transporters we used Orthofinder 2.5.2 on all sequences from all genomes with Pfam annotation PF00083 (Sugar_tr; http://pfam.xfam.org/family/sugar_tr; accessed August 5, 2021) combined with characterized sugar transporter sequences from the PF00083 seed set. Additionally, we added characterized cellodextrin (MH648002.1 (NCBI; from Aspergillus niger), S8AIR7 (UniProtKB; from Penicillium oxalicum)) and sugar alcohol transporters (AAX98668.1; from Ambrosiozyma monospora, CAR65543.1, CAG86001.1; from Debaryomyces hansenii, NP_010036.1; from Saccharomyces cerevisiae). We used the presence of characterized sequences in the inferred orthogroups to identify orthologs in each genome included in this study. The number of orthologs of different transporters per genome were visualized using custom R scripts.

Identification of class II peroxidases

First, we used diamond 0.9.22 to search all Ascomycota class II peroxidases downloaded from RedOxiBase (accessed Jul. 14, 2021) against the predicted proteins in all 83 genomes studied here. We then extracted all predicted protein sequences which had a diamond hit to any of the downloaded sequences. Similar to the identification of sugar transporters we used Orthofinder 2.5.2 to classify putative class II peroxidases present in our genomes based on the presence of downloaded genes in individual orthogroups. The number of orthologs per genome for different peroxidases were visualized using custom R scripts.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

De novo generated genome assemblies and corresponding functional annotations are deposited at NCBI under BioProject PRJNA795879. Accession numbers for all genomes used in this study are provided in Supplementary Data 1. Used PFAM and RedoxiBase sequences, annotations of previously sequenced genomes, alignments and phylogenomic trees are available for download in a public Dryad repository (https://doi.org/10.5061/dryad.3xsj3txjb). Source Data are provided with this paper.

Code availability

The analysis workflow used to acquire the results in this paper, including all custom python and R scripts, is available on Github70 (https://github.com/reslp/LFS-cazy-comparative; https://doi.org/10.5281/zenodo.6453144).

References

  1. Schwendener, S. Die Algentypen der Flechtengonidien. (Universitaetsbuchdruckerei, Basel, 1869).

    MATH  Google Scholar 

  2. Smith, D. C., Muscatine, L. & Lewis, D. Carbohydrate movement from autotrophs to heterotrophs in parasitic and mutualistic symbiosis. Biol. Rev. 44, 17–90. (1969).

    CAS  PubMed  Article  Google Scholar 

  3. Drew, E. A. & Smith, D. C. Studies in the physiology of lichens. VII. The physiology of the Nostoc symbiont of Peltigera polydactyla compared with cultured and free-living forms. N. Phytologist 66, 379–388 (1967).

    Article  Google Scholar 

  4. Richardson, D. L., Hill, D. J. & Smith, D. Lichen Physiology - XI. The role of the alga in determining the pattern of carbohydrate movement between lichen symbionts. N. Phytologist 67, 469–486 (1968).

    CAS  Article  Google Scholar 

  5. Smith, D. C. Mechanisms of nutrient movement between the lichen symbionts. In Cellular interactions in symbiosis and parasitism (eds Cook, C. B., Pappas, P. W. & Rudolph, E. D.). 197–227 (Ohio State University Press, 1980).

  6. Palmqvist, K., Dahlman, L., Jonsson, A. & Nash, T. H. III. The carbon economy of lichens. In Lichen Biology. 2nd edn. (ed. Nash, T.H. III) 182–215 (Cambridge University Press, 2008).

  7. Scott, G. D. Evolutionary aspects of symbiosis. In: The Lichens (Ahmadjian, V. & Hale, M. E.) 581–598 (Academic Press, Inc., New York and London, 1973).

  8. Wipf, D., Krajinski, F., van Tuinen, D., Recorbet, G. & Courty, P.-E. Trading on the arbuscular mycorrhiza market: from arbuscules to common mycorrhizal networks. N. Phytologist 223, 1127–1142 (2019).

    CAS  Article  Google Scholar 

  9. Miyauchi, S. et al. Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits. Nat. Commun. 11, 5125 (2020).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. Ellers, J., Kiers, E. T., Currie, C. R., McDonald, B. R. & Visser, B. Ecological interactions drive evolutionary loss of traits. Ecol. Lett. 15, 1071–1082 (2012).

    PubMed  Article  Google Scholar 

  11. Lutzoni, F., Pagel, M. & Reeb, V. Major fungal lineages are derived from lichen symbiotic ancestors. Nature 411, 937–940 (2001).

    ADS  CAS  PubMed  Article  Google Scholar 

  12. Lumbsch, H. T., Schmitt, I., Döring, H. & Wedin, M. ITS sequence data suggest variability of ascus types and support ontogenetic characters as phylogenetic discriminators in the Agyriales (Ascomycota). Mycological Res. 105, 265–274 (2001).

    CAS  Article  Google Scholar 

  13. Lumbsch, H. T., Schmitt, I., Lücking, R., Wiklund, E. & Wedin, M. The phylogenetic placement of Ostropales within Lecanoromycetes (Ascomycota) revisited. Mycological Res. 111, 257–267 (2007).

    Article  Google Scholar 

  14. Thiyagaraja, V. et al. The evolution of life modes in Stictidaceae, with three novel taxa. J. Fungi 7(2), 105 (2021).

    CAS  Article  Google Scholar 

  15. Wedin, M., Döring, H. & Gilenstam, G. Saprotrophy and lichenization as options for the same fungal species on different substrata: environmental plasticity and fungal lifestyles in the Stictis–Conotrema complex. N. Phytologist 164, 459–465 (2004).

    Article  Google Scholar 

  16. Muggia, L., Baloch, E., Stabentheiner, E., Grube, M. & Wedin, M. Photobiont association and genetic diversity of the optionally lichenized fungus Schizoxylon albescens. FEMS Microbiol. Ecol. 75, 255–272 (2011).

    CAS  PubMed  Article  Google Scholar 

  17. Brodo, I. M. Substrate ecology. In The Lichens. (eds Ahmadjian, V.) 401–441 (Academic Press, Inc, 1973).

  18. Spribille, T., Thor, G., Bunnell, F. L., Goward, T. & Björk, C. R. Lichens on dead wood: species-substrate relationships in the epiphytic lichen floras of the Pacific Northwest and Fennoscandia. Ecography 31, 741–750 (2008).

    Article  Google Scholar 

  19. Resl, P., Fernández-Mendoza, F., Mayrhofer, H. & Spribille, T. The evolution of fungal substrate specificity in a widespread group of crustose lichens. Proc. R. Acad. Sci. B 285, 20180640 (2018).

    Google Scholar 

  20. Fahselt, D. Carbon metabolism in lichens. Symbiosis 17, 127–182 (1994).

    CAS  Google Scholar 

  21. Beckett, R. P., Zavarzina, A. G. & Liers, C. Oxidoreductases and cellulases in lichens: possible roles in lichen biology and soil organic matter turnover. Fungal Biol. 117, 431–443 (2013).

    CAS  PubMed  Article  Google Scholar 

  22. Cantarel, B. L. et al. The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238 (2009).

    CAS  PubMed  Article  Google Scholar 

  23. Spribille, T., Resl, P., Stanton, D. E. & Tagirdzhanova, G. Evolutionary biology of lichen symbioses. New Phytologist https://doi.org/10.1111/nph.18048 (2022).

  24. Kobayashi, Y., et al. The genome of Rhizophagus clarus HR1 reveals a common genetic basis for auxotrophy among arbuscular mycorrhizal fungi. BMC Genom. 1–11. https://doi.org/10.1186/s12864-018-4853-0 (2018).

  25. Kohler, A. et al. Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. Nat. Genet. 47, 410–415 (2015).

    CAS  PubMed  Article  Google Scholar 

  26. Lutzoni, F. et al. Contemporaneous radiations of fungi and plants linked to symbiosis. Nat. Commun. 9, 5451 (2018).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. Li, Y. et al. A genome-scale phylogeny of the kingdom Fungi. Curr. Biol. 31, 1653–1665 (2021).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. Miadlikowska, J. et al. A multigene phylogenetic synthesis for the class Lecanoromycetes (Ascomycota): 1307 fungi representing 1139 infrageneric taxa, 317 genera and 66 families. Mol. Phylogenetics Evol. 79, 132–168 (2014).

    Article  Google Scholar 

  29. Nelsen, M. P., Lücking, R., Boyce, C. K., Lumbsch, H. T. & Ree, R. H. The macroevolutionary dynamics of symbiotic and phenotypic diversification in lichens. Proc. Natl Acad. Sci. 117, 21495–21503 (2020).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. Nakamura, M. A., Nascimento, A. S. & Polikarpov, I. Structural diversity of carbohydrate esterases. Biotechnol. Res. Innov. 1, 35–51 (2017).

    Article  Google Scholar 

  31. Savelli, B. et al. RedoxiBase: a database for ROS homeostasis regulated proteins. Redox Biol. 26, 101247 (2019).

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. Plett, J. M. & Martin, F. M. Know your enemy, embrace your friend: using omics to understand how plants respond differently to pathogenic and mutualistic microorganisms. Plant J. 93, 729–746 (2018).

    CAS  PubMed  Article  Google Scholar 

  33. Yagüe, E. & Estévez, M. P. Cellobiose induces extracellular β-1,4-glucanase and β-glucosidase in the epiphytic lichen Evernia prunastri. Plant Physiol. Biochem. 28, 203–207 (1990).

    Google Scholar 

  34. Yagüe, E. & Estévez, M. P. β-1,4-glucanase (cellulase) location in the symbionts of Evernia prunastri. Lichenologist 21, 147–151 (1989).

    Article  Google Scholar 

  35. Honegger, R. The symbiotic phenotype of lichen-forming ascomycetes and their endo- and epibionts. In Fungal Associations, 2nd edn. The Mycota IX, (ed. Hock, B.) 287–338. (Springer-Verlag, 2012).

  36. Elenkin, A. K voprosu o „vnutrennem‘ saprofitizme” („endosaprofitizm”) u lishaynikov’ [On the question of “inner saprotrophy” (“endosaprophytism”) in lichens]. Bull. du Jard. Imp.érial Botanique de. St. -Pétersbourg 2, 65–84 (1902).

    Google Scholar 

  37. Kono, M., Kon, Y., Ohmura, Y., Satta, Y. & Terai, Y. In vitro resynthesis of lichenization reveals the genetic background of symbiosis-specific fungal-algal interaction in Usnea hakonensis. BMC Genom. 21, 671 (2020).

    CAS  Article  Google Scholar 

  38. Armaleo, D. et al. The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata. BMC Genom. 20, 605 (2019).

  39. Domozych, D. S. et al. The cell walls of green algae: a journey through evolution and diversity. Front. Plant Sci. 3, 82 (2012).

  40. Yagüe, E., Orus, M. I. & Estévez, M. P. Extracellular polysaccharidases synthesized by the epiphytic lichen Evernia prunastri (L.) Ach. Planta 160, 212–216 (1984).

    PubMed  Article  Google Scholar 

  41. Guerrero, J., Yagüe, E. & Estévez, M. P. Cellulase production by lichens. J. Plant Physiol. 140, 508–510 (1992).

    CAS  Article  Google Scholar 

  42. de los Ríos, A., Ramírez, R. & Estévez, P. Production of several isoforms of ß-1,4-glucanase by the cyanolichen Peltigera canina. Physiologia Plant. 100, 159–164 (1997).

    Article  Google Scholar 

  43. Tolpysheva, T. Y. Dependence of the cellulosolytic activity in the lichen Hypogymnia physodes (L.) Nyl. on the substrate. Mosc. Univ. Biol. Sci. Bull. 62, 65–68 (2007).

    Article  Google Scholar 

  44. Beckett, R. P. et al. Role of laccases and peroxidases in saprotrophic activities in the lichen Usnea undulata. Fungal Ecol. 14, 71–78 (2015).

    Article  Google Scholar 

  45. Liers, C., Ullrich, R., Hofrichter, M., Minibayeva, F. V. & Beckett, R. P. A heme peroxidase of the ascomyceteous lichen Leptogium saturninum oxidizes high-redox potential substrates. Fungal Genet. Biol. 48, 1139–1145 (2011).

    CAS  PubMed  Article  Google Scholar 

  46. Baloch, E., Lücking, R., Lumbsch, H. T. & Wedin, M. Major clades and phylogenetic relationships between lichenized and non-lichenized lineages in Ostropales (Ascomycota: Lecanoromycetes). Taxon 59, 1483–1494 (2010).

    Article  Google Scholar 

  47. Grimm, M. et al. The lichens’ microbiota, still a mystery. Front. Microbiol. 12, 623839 (2021).

    PubMed  PubMed Central  Article  Google Scholar 

  48. Frey, E. Die Spezifität der Flechtengonidien. Alte und neue Probleme. Ber. der schweizerischen botanischen Ges. 41(2), 180–198 (1932).

    Google Scholar 

  49. Armstrong, R. A. & Smith, S. N. Carbohydrates in the hypothallus and areolae of the crustose lichen Rhizocarpon geographicum (L.) DC. Symbiosis 49, 95 (2009).

    CAS  Article  Google Scholar 

  50. Hammer, S. Prothallus structure in. Cladonia. Bryologist 99, 212–217 (1996).

    Article  Google Scholar 

  51. Orus, M. I. & Ascaso, C. Localización de hifas liquénicas en los tejidos conductores de Quercus rotundifolia Lam. Collectanea Botanica 13, 325–338 (1982).

    Google Scholar 

  52. Greenhalgh, G. N. & Anglesea, D. The distribution of algal cells in lichen thalli. Lichenologist 11, 283–292 (1979).

    Article  Google Scholar 

  53. Moyo, C. E., Beckett, R. P., Trifonova, T. V. & Minibayeva, F. V. Extracellular redox cycling and hydroxyl radical production occurs widely in lichenized Ascomycetes. Fungal Biol. 121, 582–588 (2017).

    CAS  PubMed  Article  Google Scholar 

  54. Lücking, R., Hodkinson, B. P. & Leavitt, S. D. The 2016 classification of lichenized fungi in the Ascomycota and Basidiomycota – Approaching one thousand genera. Bryologist 119, 361–416 (2017).

    Article  Google Scholar 

  55. Farrar, J. F. Physiological buffering. In CRC Handbook of Lichenology. Volume II. (ed. Galun, M.) 101–106 (CRC Press, Boca Raton, FL, USA, 1988).

  56. Smith, D. C. Is a lichen a good model of biological interactions in nutrient-poor environments? In Strategies of microbial life in extreme environments. Dahlem Konferenzen, Life Sciences Research Report 13, (ed. Shilo, M.) 291–303 (Verlag Chemie, 1978).

  57. Smith, D. C. The symbiotic condition. Symbiosis 14, 3–15 (1992).

    CAS  Google Scholar 

  58. Yamamoto, Y., Miura, Y., Higuchi, M., Kinoshita, Y. & Yoshimura, I. Using lichen tissue cultures in modern biology. Bryologist 96, 384–393 (1993).

    Article  Google Scholar 

  59. Parrent, J. L. et al. Friend or foe? Evolutionary history of glycoside hydrolase family 32 genes encoding for sucrolytic activity in fungi and its implications for plant-fungal symbioses. BMC Evol. Biol. 9, 148 (2009).

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  60. Johnson, G. T. Contributions to the study of the Trypetheliaceae. Ann. Mo. Botanical Gard. 27, 1–50 (1940).

    Article  Google Scholar 

  61. Honegger, R. Ultrastructural studies in lichens. I. Haustorial types and their frequencies in a range of lichens with trebouxioid photobionts. N. Phytologist 103, 785–795 (1986).

    Article  Google Scholar 

  62. Poelt, J. Die Entstehung einer Strauchflechte aus einem Formenkreis krustiger Verwandter. Flora 183, 65–72 (1989).

    Article  Google Scholar 

  63. Spribille, T. et al. 3D biofilms: in search of the polysaccharides holding together lichen symbioses. FEMS Microbiol. Lett. https://doi.org/10.1093/femsle/fnaa023 (2020).

  64. Hale, M. E. Fine structure of the cortex in the lichen family Parmeliaceae viewed with the scanning electron microscope. Smithson. Contr Bot. 10, 1–92 (1973).

    Article  Google Scholar 

  65. Honegger, R. Functional aspects of the lichen symbiosis. Ann. Rev. Plant Physiol. Plant Mol. Biol. 42, 553–578 (1991).

    CAS  Article  Google Scholar 

  66. Tapper, R. Glucose uptake by Trebouxia and associated fungal symbiont in the lichen symbiosis. FEMS Microbiol. Lett. 10, 103–106 (1981).

    CAS  Article  Google Scholar 

  67. Campbell, J., Bengtson, J., Fredeen, A. L., Coxson, D. S. & Prescott, C. E. Does exogenous carbon extend the realized niche of canopy lichens? Evidence from sub-boreal forests in British Columbia. Ecology 94, 1186–1195 (2013).

    PubMed  Article  Google Scholar 

  68. Jahn, L. et al. Cyanodermella asteris sp. nov. from the inflorescence of Aster tataricus. Mycotaxon 132, 107–123 (2017).

    Article  Google Scholar 

  69. Tagirdzhanova, G. et al. Predicted input of uncultured fungal symbionts to a lichen symbiosis from metagenome-assembled genomes. Genome Biol. Evol. 13, evab047 (2021).

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  70. Resl, P. reslp/LFS-cazy-comparative: LFS-cazy-comparative-initial-release-V1.0.0 (v1.0.0). Zenodo. https://doi.org/10.5281/zenodo.6453144 (2022).

Download references

Acknowledgements

The early phases of this project were funded by the Austrian Science Fund (FWF grant P25237, “Evolution of Substrate Specificity in Lichens”) and carried out at the Institute of Plant Sciences (Uni Graz). PR would like to acknowledge the Theodor Körner Funds (Vienna, Austria) and Network of Biological Systematics Austria (NOBIS; Vienna, Austria) for funding that enabled genome sequencing of several species. PR would also like to thank Christoph Hahn and Fernando Fernández-Mendoza for many methodological discussions. SW received funding from the Icelandic Research Fund IRF (174307-051), from Deutsche Forschungsgemeinschaft DFG (WE 6443/1-1) and from LMU Munich (startup funds). MW received funding from the Swedish Research Council, grant VR 2016-03589. MW and MK would like to acknowledge support from the NRM Department of Bioinformatics and Genetics, the National Genomics Infrastructure in Stockholm, the SNIC/Uppsala Multidisciplinary Center for Advanced Computational Science and the UPPMAX computational infrastructure, and Linda Phillips for assistance with material. TS would like to acknowledge an NSERC Discovery Grant, a Canada Research Chair in Symbiosis, and the generosity of Susan Dalby and Eskild Petersen, who provided a place to work on this manuscript. Special thanks go to Sophie Dang at the Molecular Biology Service Unit, University of Alberta Department of Biological Sciences, for help in data acquisition, and members of the U of A Lichen Evolution Lab for reading and commenting on the manuscript. We are also grateful to Sigrun Kraker, Theodora Kopun, Tanja Ernst and Andrea Brandl for laboratory assistance.

Author information

Authors and Affiliations

Authors

Contributions

T.S., P.R. and L.M.W. designed this study. Genomes used for this study were generated in the labs of T.S., H.M., J.P.M., S.W. and M.W. G.T. performed metagenome filtering; S.F.R. assembled genomes of Sticta, Pseudocyphellaria and Lobaria, P.M. for Xylographa species; S.F.B., M.K. and L.M. organized cultures and genome sequencing for Stictis; S.W., H.G. and O.S.A. contributed the Peltigera genome; P.R. assembled all other genomes. A.R.B. performed heterologous expression of Xylographa cellulases. P.R. performed gene-calling, functional annotation and comparative genomic analyses. T.S. and P.R. wrote the manuscript with contributions from all authors.

Corresponding author

Correspondence to Toby Spribille.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Source data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Resl, P., Bujold, A.R., Tagirdzhanova, G. et al. Large differences in carbohydrate degradation and transport potential among lichen fungal symbionts. Nat Commun 13, 2634 (2022). https://doi.org/10.1038/s41467-022-30218-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41467-022-30218-6

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing