Climate-specific biosynthetic gene clusters in populations of a lichen-forming fungus

Genome-wide differentiation corresponds to climatic niches in two species of lichen-forming fungi

Lichens can withstand fluctuating environmental conditions such as hydration-desiccation cycles. Many species distribute across climate zones, suggesting population-level adaptations to conditions such as freezing and drought. Here, we aim to understand how climate affects population genomic patterns in lichenized fungi. We analysed population structure along elevational gradients in closely related Umbilicaria phaea (North American; two gradients) and Umbilicaria pustulata (European; three gradients). All gradients showed clear genomic breaks splitting populations into low-elevation (Mediterranean zone) and high-elevation (cold temperate zone). A total of 3301 SNPs in U. phaea and 138 SNPs in U. pustulata were driven to fixation between the two ends of the gradients. The difference between the species is likely due to differences in recombination rate: the sexually reproducing U. phaea has a higher recombination rate than the primarily asexually reproducing U. pustulata. Cline analysis revealed allele frequency transitions along all gradients at approximately 0°C, coinciding with the transition between the Mediterranean and cold temperate zones, suggesting freezing is a strong driver of population differentiation. Genomic scans further confirmed temperature-related selection targets. Both species showed similar differentiation patterns overall, but different selected alleles indicate convergent adaptation to freezing. Our results enrich our knowledge of fungal genomic functions related to temperature and climate, fungal population genomics, and species responses to environmental heterogeneity.

Exploration on cold adaptation of Antarctic lichen via detection of positive selection genes

Lichen as mutualistic symbiosis is the dominant organism in various extreme terrestrial environment on Earth, however, the mechanisms of their adaptation to extreme habitats have not been fully elucidated. In this study, we chose the Antarctic dominant lichen species Usnea aurantiacoatra to generate a high-quality genome, carried out phylogenetic analysis using maximum likelihood and identify genes under positive selection. We performed functional enrichment analysis on the positively selected genes (PSGs) and found that most of the PSGs focused on transmembrane transporter activity and vacuole components. This suggest that the genes related to energy storage and transport in Antarctic U. aurantiacoatra were affected by environmental pressure. Inside of the 86 PSGs screened, two protein interaction networks were identified, which were RNA helicase related proteins and regulator of G-protein signaling related proteins. The regulator of the G-protein signaling gene (UaRGS1) was chosen to perform further verification by the lichen genetic manipulation system Umbilicaria muhlenbergii. Given that the absence of UmRgs1 resulted in elevated lethality to cold shock, the role for UaRgs1 in Antarctic U. aurantiacoatra resistance to cold can be inferred. The investigation of lichen adaptation to extreme environments at the molecular level will be opened up.

An Antarctic lichen isolate (Cladonia borealis) genome reveals potential adaptation to extreme environments

Cladonia borealis is a lichen that inhabits Antarctica's harsh environment. We sequenced the whole genome of a C. borealis culture isolated from a specimen collected in Antarctica using long-read sequencing technology to identify specific genetic elements related to its potential environmental adaptation. The final genome assembly produced 48 scaffolds, the longest being 2.2 Mbp, a 1.6 Mbp N50 contig length, and a 36 Mbp total length. A total of 10,749 protein-coding genes were annotated, containing 33 biosynthetic gene clusters and 102 carbohydrate-active enzymes. A comparative genomics analysis was conducted on six Cladonia species, and the genome of C. borealis exhibited 45 expanded and 50 contracted gene families. We identified that C. borealis has more Copia transposable elements and expanded transporters (ABC transporters and magnesium transporters) compared to other Cladonia species. Our results suggest that these differences contribute to C. borealis' remarkable adaptability in the Antarctic environment. This study also provides a useful resource for the genomic analysis of lichens and genetic insights into the survival of species isolated from Antarctica.

Lichen holobionts show compositional structure along elevation

Holobionts are dynamic ecosystems that may respond to abiotic drivers with compositional changes. Uncovering elevational diversity patterns within these microecosystems can further our understanding of community-environment interactions. Here, we assess how the major components of lichen holobionts-fungal hosts, green algal symbionts, and the bacterial community-collectively respond to an elevational gradient. We analyse populations of two lichen symbioses, Umbilicaria pustulata and U. hispanica, along an elevational gradient spanning 2100 altitudinal metres and covering three major biomes. Our study shows (i) discontinuous genomic variation in fungal hosts with one abrupt genomic differentiation within each of the two host species, (ii) altitudinally structured bacterial communities with pronounced turnover within and between hosts, and (iii) altitude-specific presence of algal symbionts. Alpha diversity of bacterial communities decreased with increasing elevation. A marked turnover in holobiont diversity occurred across two altitudinal belts: at 11°C-13°C average annual temperature (here: 800-1200 m a.s.l.), and at 7°C-9°C average annual temperature (here: 1500-1800 m a.s.l.). The two observed zones mark a clustering of distribution limits and community shifts. The three ensuing altitudinal classes, that is, the most frequent combinations of species in holobionts, approximately correspond to the Mediterranean, cool-temperate, and alpine climate zones. We conclude that multitrophic microecosystems, such as lichen holobionts, respond with concerted compositional changes to climatic factors that also structure communities of macroorganisms, for example, vascular plants.

Biosynthetic gene cluster synteny: Orthologous polyketide synthases in Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus), which together generate a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here, we provide a comparative view of the biosynthetic gene clusters of three lichen mycobionts derived from Hypogymnia physodes, Hypogymnia tubulosa, and Parmelia sulcata. In addition, we present a high-quality PacBio metagenome of Parmelia sulcata, from which we extracted the mycobiont bin containing 214 biosynthetic gene clusters. Most biosynthetic gene clusters in these genomes were associated with T1PKSs, followed by NRPSs and terpenes. This study focused on biosynthetic gene clusters related to polyketide synthesis. Based on ketosynthase homology, we identified nine highly syntenic clusters present in all three species. Among the four clusters belonging to nonreducing PKSs, two are putatively linked to lichen substances derived from orsellinic acid (orcinol depsides and depsidones, e.g., lecanoric acid, physodic acid, lobaric acid), one to compounds derived from methylated forms of orsellinic acid (beta orcinol depsides, e.g., atranorin), and one to melanins. Five clusters with orthologs in all three species are linked to reducing PKSs. Our study contributes to sorting and dereplicating the vast PKS diversity found in lichenized fungi. High-quality sequences of biosynthetic gene clusters of these three common species provide a foundation for further exploration into biotechnological applications and the molecular evolution of lichen substances.

Biosynthetic Potential of Hypogymnia Holobionts: Insights into Secondary Metabolite Pathways

Lichens are symbiotic associations consisting of a photobiont (algae or cyanobacteria) and a mycobiont (fungus). They are known to produce a variety of unique secondary metabolites. To access this biosynthetic potential for biotechnological applications, deeper insights into the biosynthetic pathways and corresponding gene clusters are necessary. Here we provide a comprehensive view of the biosynthetic gene clusters of all organisms comprising a lichen thallus: fungi, green algae, and bacteria. We present two high-quality PacBio metagenomes, in which we identified a total of 460 biosynthetic gene clusters. Lichen mycobionts yielded 73-114 clusters, other lichen associated ascomycetes 8-40, green algae of the genus Trebouxia 14-19, and lichen-associated bacteria 101-105 clusters. The mycobionts contained mainly T1PKSs, followed by NRPSs, and terpenes; Trebouxia reads harbored mainly clusters linked to terpenes, followed by NRPSs and T3PKSs. Other lichen-associated ascomycetes and bacteria contained a mix of diverse biosynthetic gene clusters. In this study, we identified for the first time the biosynthetic gene clusters of entire lichen holobionts. The yet untapped biosynthetic potential of two species of the genus Hypogymnia is made accessible for further research.

How to build a lichen: from metabolite release to symbiotic interplay

Exposing their vegetative bodies to the light, lichens are outstanding amongst other fungal symbioses. Not requiring a pre-established host, 'lichenized fungi' build an entirely new structure together with microbial photosynthetic partners that neither can form alone. The signals involved in the transition of a fungus and a compatible photosynthetic partner from a free-living to a symbiotic state culminating in thallus formation, termed 'lichenization', and in the maintenance of the symbiosis, are poorly understood. Here, we synthesise the puzzle pieces of the scarce knowledge available into an updated concept of signalling involved in lichenization, comprising five main stages: (1) the 'pre-contact stage', (2) the 'contact stage', (3) 'envelopment' of algal cells by the fungus, (4) their 'incorporation' into a pre-thallus and (5) 'differentiation' into a complex thallus. Considering the involvement of extracellularly released metabolites in each phase, we propose that compounds such as fungal lectins and algal cyclic peptides elicit early contact between the symbionts-to-be, whereas phytohormone signalling, antioxidant protection and carbon exchange through sugars and sugar alcohols are of continued importance throughout all stages. In the fully formed lichen thallus, secondary lichen metabolites and mineral nutrition are suggested to stabilize the functionalities of the thallus, including the associated microbiota.

Linking Lichen Metabolites to Genes: Emerging Concepts and Lessons from Molecular Biology and Metagenomics

Lichen secondary metabolites have tremendous pharmaceutical and industrial potential. Although more than 1000 metabolites have been reported from lichens, less than 10 have been linked to the genes coding them. The current biosynthetic research focuses strongly on linking molecules to genes as this is fundamental to adapting the molecule for industrial application. Metagenomic-based gene discovery, which bypasses the challenges associated with culturing an organism, is a promising way forward to link secondary metabolites to genes in non-model, difficult-to-culture organisms. This approach is based on the amalgamation of the knowledge of the evolutionary relationships of the biosynthetic genes, the structure of the target molecule, and the biosynthetic machinery required for its synthesis. So far, metagenomic-based gene discovery is the predominant approach by which lichen metabolites have been linked to their genes. Although the structures of most of the lichen secondary metabolites are well-documented, a comprehensive review of the metabolites linked to their genes, strategies implemented to establish this link, and crucial takeaways from these studies is not available. In this review, I address the following knowledge gaps and, additionally, provide critical insights into the results of these studies, elaborating on the direct and serendipitous lessons that we have learned from them.

Metagenomics Shines Light on the Evolution of "Sunscreen" Pigment Metabolism in the Teloschistales (Lichen-Forming Ascomycota)

Fungi produce a vast number of secondary metabolites that shape their interactions with other organisms and the environment. Characterizing the genes underpinning metabolite synthesis is therefore key to understanding fungal evolution and adaptation. Lichenized fungi represent almost one-third of Ascomycota diversity and boast impressive secondary metabolites repertoires. However, most lichen biosynthetic genes have not been linked to their metabolite products. Here we used metagenomic sequencing to survey gene families associated with production of anthraquinones, UV-protectant secondary metabolites present in various fungi, but especially abundant in a diverse order of lichens, the Teloschistales (class Lecanoromycetes, phylum Ascomycota). We successfully assembled 24 new, high-quality lichenized-fungal genomes de novo and combined them with publicly available Lecanoromycetes genomes from taxa with diverse secondary chemistry to produce a whole-genome tree. Secondary metabolite biosynthetic gene cluster (BGC) analysis showed that whilst lichen BGCs are numerous and highly dissimilar, core enzyme genes are generally conserved across taxa. This suggests metabolite diversification occurs via re-shuffling existing enzyme genes with novel accessory genes rather than BGC gains/losses or de novo gene evolution. We identified putative anthraquinone BGCs in our lichen dataset that appear homologous to anthraquinone clusters from non-lichenized fungi, suggesting these genes were present in the common ancestor of the subphylum Pezizomycotina. Finally, we identified unique transporter genes in Teloschistales anthraquinone BGCs that may explain why these metabolites are so abundant and ubiquitous in these lichens. Our results support the importance of metagenomics for understanding the secondary metabolism of non-model fungi such as lichens.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi

Introduction

Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria.

Methods

The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation.

Results

In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations.

Discussion

We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Fungal Host Affects Photosynthesis in a Lichen Holobiont

Corals and lichens are iconic examples of photosynthetic holobionts, i.e., ecological and evolutionary units resulting from the tightly integrated association of algae and prokaryotic microbiota with animal or fungal hosts, respectively. While the role of the coral host in modulating photosynthesis has been clarified to a large extent in coral holobionts, the role of the fungal host in this regard is far less understood. Here, we address this question by taking advantage of the recent discovery of highly specific fungal-algal pairings corresponding to climatically adapted ecotypes of the lichen-forming genus Umbilicaria. Specifically, we compared chlorophyll a fluorescence kinetics among lichen thalli consisting of different fungal-algal combinations. We show that photosynthetic performance in these lichens is not only driven by algal genotype, but also by fungal host species identity and intra-host genotype. These findings shed new light on the closely intertwined physiological processes of fungal and algal partners in the lichen symbiosis. Indeed, the specific combinations of fungal and algal genotypes within a lichen individual-and the resulting combined functional phenotype-can be regarded as a response to the environment. Our findings suggest that characterizing the genetic composition of both eukaryotic partners is an important complimentary step to understand and predict the lichen holobiont's responses to environmental change.

Genome mining as a biotechnological tool for the discovery of novel biosynthetic genes in lichens

Natural products (NPs) and their derivatives are a major contributor to modern medicine. Historically, microorganisms such as bacteria and fungi have been instrumental in generating drugs and lead compounds because of the ease of culturing and genetically manipulating them. However, the ever-increasing demand for novel drugs highlights the need to bioprospect previously unexplored taxa for their biosynthetic potential. Next-generation sequencing technologies have expanded the range of organisms that can be explored for their biosynthetic content, as these technologies can provide a glimpse of an organism's entire biosynthetic landscape, without the need for cultivation. The entirety of biosynthetic genes can be compared to the genes of known function to identify the gene clusters potentially coding for novel products. In this study, we mine the genomes of nine lichen-forming fungal species of the genus Umbilicaria for biosynthetic genes, and categorize the biosynthetic gene clusters (BGCs) as "associated product structurally known" or "associated product putatively novel". Although lichen-forming fungi have been suggested to be a rich source of NPs, it is not known how their biosynthetic diversity compares to that of bacteria and non-lichenized fungi. We found that 25%-30% of biosynthetic genes are divergent as compared to the global database of BGCs, which comprises 1,200,000 characterized biosynthetic genes from plants, bacteria, and fungi. Out of 217 BGCs, 43 were highly divergant suggesting that they potentially encode structurally and functionally novel NPs. Clusters encoding the putatively novel metabolic diversity comprise polyketide synthases (30), non-ribosomal peptide synthetases (12), and terpenes (1). Our study emphasizes the utility of genomic data in bioprospecting microorganisms for their biosynthetic potential and in advancing the industrial application of unexplored taxa. We highlight the untapped structural metabolic diversity encoded in the lichenized fungal genomes. To the best of our knowledge, this is the first investigation identifying genes coding for NPs with potentially novel properties in lichenized fungi.

A genomic journey in the secondary metabolite diversity of fungal plant and insect pathogens: from functional to population genomics

Fungal pathogens produce a broad array of secondary metabolites (SMs), which allow the fungus to thrive in its natural habitat and gain competitive advantage. Analysis of the genetically encoded blueprints for SM assembly highlighted that only a small portion of the SMs these fungi are capable of producing are known, and even fewer have been investigated for their natural function. Using molecular tools, a lot of progress has been made recently in identifying the blueprint products and linking them to their ecological purpose such as the peptide virulence factor fusaoctaxin A released by Fusarium graminearum during infection of wheat or the F. oxysporum polyketide bikaverin that provides competitive advantage against bacteria in tomato. In addition, population genomics have given particularly important insights into the species-specific plasticity of the SM blueprint arsenal, showcasing the ongoing evolution and adaptation of fungal pathogens. This approach holds promise in inferring roles in pathogenicity of many more fungal SMs.

A Candidate Gene Cluster for the Bioactive Natural Product Gyrophoric Acid in Lichen-Forming Fungi

Natural products of lichen-forming fungi are structurally diverse and have a variety of medicinal properties. Despite this, they have limited implementation in industry mostly because the corresponding genes are unknown for most of their natural products. Here, we implement a long-read sequencing and bioinformatic approach to identify the putative biosynthetic gene cluster of the bioactive natural product gyrophoric acid (GA). Using 15 high-quality genomes representing nine GA-producing species of the lichen-forming fungal genus Umbilicaria, we identify the most likely GA cluster and investigate the cluster gene organization and composition across the nine species. Our results show that GA clusters are promiscuous within Umbilicaria, and only three genes are conserved across species, including the polyketide synthase (PKS) gene. In addition, our results suggest that the same cluster codes for different, but structurally similar compounds, namely, GA, umbilicaric-, and hiascic acid, bringing new evidence that lichen metabolite diversity is also generated through regulatory mechanisms at the molecular level. Ours is the first study to identify the most likely GA cluster and, thus, provides essential information to open new avenues for biotechnological approaches to producing and modifying GA and similar lichen-derived compounds. GA PKS is the first tridepside PKS to be identified.

IMPORTANCE The implementation of natural products in the pharmaceutical industry relies on the possibility of modifying the natural product (NP) pathway to optimize yields and pharmacological effects. Characterization of genes and pathways underlying natural product biosynthesis is a major bottleneck for exploiting the medicinal properties of the natural products. Genome mining is a promising and relatively cost- and time-effective approach to utilize unexplored NP resources for drug discovery. In this study, we identify the most likely gene cluster for the lichen-forming fungal depside gyrophoric acid in nine Umbilicaria species. This compound shows cytotoxic and antiproliferative properties against several cancer cell lines and is also a broad-spectrum antimicrobial agent. This information paves the way for generating GA analogs with modified properties by selective activation/deactivation of genes.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.

Chemical properties of key metabolites determine the global distribution of lichens

In lichen symbioses, fungal secondary metabolites provide UV protection on which lichen algae such as trebouxiophycean green algae-the most prominent group of photobionts in lichen symbioses-sensitively depend. These metabolites differ in their UV absorbance capability and solvability, and thus vary in their propensity of being leached from the lichen body in humid and warm environments, with still unknown implications for the global distribution of lichens. In this study covering more than 10,000 lichenised fungal species, we show that the occurrence of fungal-derived metabolites in combination with their UV absorbance capability and their probability of being leached in warm and humid environments are important eco-evolutionary drivers of global lichen distribution. Fungal-derived UV protection seems to represent an indirect environmental adaptation in which the lichen fungus invests to protect the trebouxiophycean photobiont from high UV radiation in warm and humid climates and, by doing this, secures its carbon source.

Depside and Depsidone Synthesis in Lichenized Fungi Comes into Focus through a Genome-Wide Comparison of the Olivetoric Acid and Physodic Acid Chemotypes of Pseudevernia furfuracea

Primary biosynthetic enzymes involved in the synthesis of lichen polyphenolic compounds depsides and depsidones are non-reducing polyketide synthases (NR-PKSs), and cytochrome P450s. However, for most depsides and depsidones the corresponding PKSs are unknown. Additionally, in non-lichenized fungi specific fatty acid synthases (FASs) provide starters to the PKSs. Yet, the presence of such FASs in lichenized fungi remains to be investigated. Here we implement comparative genomics and metatranscriptomics to identify the most likely PKS and FASs for olivetoric acid and physodic acid biosynthesis, the primary depside and depsidone defining the two chemotypes of the lichen Pseudevernia furfuracea. We propose that the gene cluster PF33-1_006185, found in both chemotypes, is the most likely candidate for the olivetoric acid and physodic acid biosynthesis. This is the first study to identify the gene cluster and the FAS likely responsible for olivetoric acid and physodic acid biosynthesis in a lichenized fungus. Our findings suggest that gene regulation and other epigenetic factors determine whether the mycobiont produces the depside or the depsidone, providing the first direct indication that chemotype diversity in lichens can arise through regulatory and not only through genetic diversity. Combining these results and existing literature, we propose a detailed scheme for depside/depsidone synthesis.