High Diversity of Type I Polyketide Genes in Bacidia rubella as Revealed by the Comparative Analysis of 23 Lichen Genomes

Genomic features of lichen-associated black fungi

Lichens are mutualistic associations consisting of a primary fungal host, and one to few primary phototrophic symbiont(s), usually a green alga and/or a cyanobacterium. They form complex thallus structures, which provide unique and stable habitats for many other microorganisms. Frequently isolated from lichens are the so-called black fungi, or black yeasts, which are mainly characterized by melanized cell walls and extremophilic lifestyles. It is presently unclear in which ways these fungi interact with other members of the lichen symbiosis. Genomic resources of lichen-associated black fungi are needed to better understand the physiological potential of these fungi and shed light on the complexity of the lichen consortium. Here, we present high-quality genomes of 14 black fungal lineages, isolated from lichens of the rock-dwelling genus Umbilicaria. Nine of the lineages belong to the Eurotiomycetes (Chaetothyriales), four to the Dothideomycetes, and one to the Arthoniomycetes, representing the first genome of a black fungus in this class. The PacBio-based assemblies are highly contiguous (5-42 contigs per genome, mean coverage of 79-502, N50 of 1.0-7.3 mega-base-pair (Mb), Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness generally ≥95.4%). Most contigs are flanked by a telomere sequence, suggesting we achieved near chromosome-level assemblies. Genome sizes range between 26 and 44 Mb. Transcriptome-based annotations yielded ~11,000-18,000 genes per genome. We analyzed genome content with respect to repetitive elements, biosynthetic genes, and effector genes. Each genome contained a polyketide synthase gene related to the dihydroxynaphthalene-melanin pathway. This research provides insights into genome content and metabolic potential of these relatively unknown, but frequently encountered lichen associates.

Are there conserved biosynthetic genes in lichens? Genome-wide assessment of terpene biosynthetic genes suggests ubiquitous distribution of the squalene synthase cluster

Lichen-forming fungi (LFF) are prolific producers of functionally and structurally diverse secondary metabolites, most of which are taxonomically exclusive and play lineage-specific roles. To date, widely distributed, evolutionarily conserved biosynthetic pathways in LFF are not known. However, this idea stems from polyketide derivatives, since most biochemical research on lichens has concentrated on polyketide synthases (PKSs). Here, we present the first systematic identification and comparison of terpene biosynthetic genes of LFF using all the available Lecanoromycete reference genomes and 22 de novo sequenced ones (111 in total, representing 60 genera and 23 families). We implemented genome mining and gene networking approaches to identify and group the biosynthetic gene clusters (BGCs) into networks of similar BGCs. Our large-scale analysis led to the identification of 724 terpene BGCs with varying degrees of pairwise similarity. Most BGCs in the dataset were unique with no similarity to a previously known fungal or bacterial BGC or among each other. Remarkably, we found two BGCs that were widely distributed in LFF. Interestingly, both conserved BGCs contain the same core gene, i.e., putatively a squalene/phytoene synthase (SQS), involved in sterol biosynthesis. This indicates that early gene duplications, followed by gene losses/gains and gene rearrangement are the major evolutionary factors shaping the composition of these widely distributed SQS BGCs across LFF. We provide an in-depth overview of these BGCs, including the transmembrane, conserved, variable and LFF-specific regions. Our study revealed that lichenized fungi do have a highly conserved BGC, providing the first evidence that a biosynthetic gene may constitute essential genes in lichens.

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.

Identification of a biosynthetic gene cluster for a red pigment cristazarin produced by a lichen-forming fungus Cladonia metacorallifera

Lichens are known to produce many novel bioactive metabolites. To date, approximately 1,000 secondary metabolites have been discovered, which are predominantly produced by the lichen mycobionts. However, despite the extensive studies on production of lichen secondary metabolites, little is known about the responsible biosynthetic gene clusters (BGCs). Here, we identified a putative BGC that is implicated in production of a red pigment, cristazarin (a naphthazarin derivative), in Cladonia metacorallifera. Previously, cristazarin was shown to be specifically induced in growth media containing fructose as a sole carbon source. Thus, we performed transcriptome analysis of C. metacorallifera growing on different carbon sources including fructose to identify the BGC for cristazarin. Among 39 polyketide synthase (PKS) genes found in the genome of C. metacorallifera, a non-reducing PKS (coined crz7) was highly expressed in growth media containing either fructose or glucose. The borders of a cristazarin gene cluster were delimited by co-expression patterns of neighboring genes of the crz7. BGCs highly conserved to the cristazarin BGC were also found in C. borealis and C. macilenta, indicating that these related species also have metabolic potentials to produce cristazarin. Phylogenetic analysis revealed that the Crz7 is sister to fungal PKSs that biosynthesize an acetylated tetrahydoxynaphthalene as a precursor of melanin pigment. Based on the phylogenetic placement of the Crz7 and putative functions of its neighboring genes, we proposed a plausible biosynthetic route for cristazarin. In this study, we identified a lichen-specific BGC that is likely involved in the biosynthesis of a naphthazarin derivative, cristazarin, and confirmed that transcriptome profiling under inducing and non-inducing conditions is an effective strategy for linking metabolites of interest to biosynthetic genes.

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.

Discovery and excavation of lichen bioactive natural products

Lichen natural products are a tremendous source of new bioactive chemical entities for drug discovery. The ability to survive in harsh conditions can be directly correlated with the production of some unique lichen metabolites. Despite the potential applications, these unique metabolites have been underutilized by pharmaceutical and agrochemical industries due to their slow growth, low biomass availability, and technical challenges involved in their artificial cultivation. At the same time, DNA sequence data have revealed that the number of encoded biosynthetic gene clusters in a lichen is much higher than in natural products, and the majority of them are silent or poorly expressed. To meet these challenges, the one strain many compounds (OSMAC) strategy, as a comprehensive and powerful tool, has been developed to stimulate the activation of silent or cryptic biosynthetic gene clusters and exploit interesting lichen compounds for industrial applications. Furthermore, the development of molecular network techniques, modern bioinformatics, and genetic tools is opening up a new opportunity for the mining, modification, and production of lichen metabolites, rather than merely using traditional separation and purification techniques to obtain small amounts of chemical compounds. Heterologous expressed lichen-derived biosynthetic gene clusters in a cultivatable host offer a promising means for a sustainable supply of specialized metabolites. In this review, we summarized the known lichen bioactive metabolites and highlighted the application of OSMAC, molecular network, and genome mining-based strategies in lichen-forming fungi for the discovery of new cryptic lichen compounds.

Genomic Analysis of the Xanthoria elegans and Polyketide Synthase Gene Mining Based on the Whole Genome

Xanthoria elegans is a lichen symbiosis, that inhabits extreme environments and can absorb UV-B. We reported the de novo sequencing and assembly of X. elegans genome. The whole genome was approximately 44.63 Mb, with a GC content of 40.69%. Genome assembly generated 207 scaffolds with an N50 length of 563,100 bp, N90 length of 122,672 bp. The genome comprised 9,581 genes, some encoded enzymes involved in the secondary metabolism such as terpene, polyketides. To further understand the UV-B absorbing and adaptability to extreme environments mechanisms of X. elegans, we searched the secondary metabolites genes and gene-cluster from the genome using genome-mining and bioinformatics analysis. The results revealed that 7 NR-PKSs, 12 HR-PKSs and 2 hybrid PKS-PKSs from X. elegans were isolated, they belong to Type I PKS (T1PKS) according to the domain architecture; phylogenetic analysis and BGCs comparison linked the putative products to two NR-PKSs and three HR-PKSs, the putative products of two NR-PKSs were emodin xanthrone (most likely parietin) and mycophelonic acid, the putative products of three HR-PKSs were soppilines, (+)-asperlin and macrolactone brefeldin A, respectively. 5 PKSs from X. elegans build a correlation between the SMs carbon skeleton and PKS genes based on the domain architecture, phylogenetic and BGC comparison. Although the function of 16 PKSs remains unclear, the findings emphasize that the genes from X. elegans represent an unexploited source of novel polyketide and utilization of lichen gene resources.

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.

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.