Draft Genome Sequence of Lichen-Forming Fungus Cladonia metacorallifera Strain KoLRI002260

Genomic Resources for the First Federally Endangered Lichen: The Florida Perforate Cladonia (Cladonia perforata)

Thirty years after its designation as a federally endangered species, the Florida Perforate Cladonia (FPC) remains imperiled in isolated populations in the Florida scrub in the southeastern USA. For threatened and endangered species, such as FPC, reference genomes provide critical insight into genomic diversity, local adaptations, landscape-level genetics, and phylogenomics. Using high-throughput sequencing, we assemble the first draft nuclear and mitochondrial genomes for the FPC mycobiont-Cladonia perforata. We also assess genetic diversity within and among populations in southeastern Florida using genome-scale data and investigate diversity across the entire nuclear ribosomal cistron, including the standard DNA barcoding marker for fungi. The draft nuclear genome spanned 33.6 Mb, and the complete, circular mitochondrial genome was 59 Kb. We also generated the first chloroplast genome, to our knowledge, for the photobiont genus associated with FPC, an undescribed Asterochloris species. We inferred the presence of multiple, distinct mycobiont parental genotypes (genets) occurring at local scales in southeastern Florida, and strikingly, no genets were shared among even the closest sample sites. All sampled thalli shared identical mitochondrial genomes, while the nuclear ribosomal cistron showed limited variability-highlighting the genetic resolution provided by nuclear genome-scale datasets. The genomic resources generated here provide critical resources for informed conservation efforts for the FPC.

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.

A comparative genomic analysis of lichen-forming fungi reveals new insights into fungal lifestyles

Lichen-forming fungi are mutualistic symbionts of green algae or cyanobacteria. We report the comparative analysis of six genomes of lichen-forming fungi in classes Eurotiomycetes and Lecanoromycetes to identify genomic information related to their symbiotic lifestyle. The lichen-forming fungi exhibited genome reduction via the loss of dispensable genes encoding plant-cell-wall-degrading enzymes, sugar transporters, and transcription factors. The loss of these genes reflects the symbiotic biology of lichens, such as the absence of pectin in the algal cell wall and obtaining specific sugars from photosynthetic partners. The lichens also gained many lineage- and species-specific genes, including those encoding small secreted proteins. These genes are primarily induced during the early stage of lichen symbiosis, indicating their significant roles in the establishment of lichen symbiosis.Our findings provide comprehensive genomic information for six lichen-forming fungi and novel insights into lichen biology and the evolution of symbiosis.

Production and Activity of Cristazarin in the Lichen-Forming Fungus Cladonia metacorallifera

Lichens are a natural source of bioactive compounds. Cladonia metacorallifera var. reagens KoLRI002260 is a rare lichen known to produce phenolic compounds, such as rhodocladonic, thamnolic, and didymic acids. However, these metabolites have not been detected in isolated mycobionts. We investigated the effects of six carbon sources on metabolite biosynthesis in the C. metacorallifera mycobiont. Red pigments appeared only in Lilly and Barnett's media with fructose at 15 °C after 3 weeks of culture and decreased after 6 weeks. We purified these red pigments using preparative-scale high performance liquid chromatography and analyzed them via nuclear magnetic resonance. Results indicated that 1% fructose-induced cristazarin and 6-methylcristazarin production under light conditions. In total, 27 out of 30 putative polyketide synthase genes were differentially expressed after 3 weeks of culture, implying that these genes may be required for cristazarin production in C. metacorallifera. Moreover, the white collar genes Cmwc-1 and Cmwc-2 were highly upregulated at all times under light conditions, indicating a possible correlation between cristazarin production and gene expression. The cancer cell lines AGS, CT26, and B16F1 were sensitive to cristazarin, with IC50 values of 18.2, 26.1, and 30.9 μg/mL, respectively, which highlights the value of cristazarin. Overall, our results suggest that 1% fructose under light conditions is required for cristazarin production by C. metacorallifera mycobionts, and cristazarin could be a good bioactive compound.

Microbial Communities of Cladonia Lichens and Their Biosynthetic Gene Clusters Potentially Encoding Natural Products

Lichens have been widely used in traditional medicine, especially by indigenous communities worldwide. However, their slow growth and difficulties in the isolation of lichen symbionts and associated microbes have hindered the pharmaceutical utilisation of lichen-produced compounds. Advances in high-throughput sequencing techniques now permit detailed investigations of the complex microbial communities formed by fungi, green algae, cyanobacteria, and other bacteria within the lichen thalli. Here, we used amplicon sequencing, shotgun metagenomics, and in silico metabolomics together with compound extractions to study reindeer lichens collected from Southern Finland. Our aim was to evaluate the potential of Cladonia species as sources of novel natural products. We compared the predicted biosynthetic pathways of lichen compounds from isolated genome-sequenced lichen fungi and our environmental samples. Potential biosynthetic genes could then be further used to produce secondary metabolites in more tractable hosts. Furthermore, we detected multiple compounds by metabolite analyses, which revealed connections between the identified biosynthetic gene clusters and their products. Taken together, our results contribute to metagenomic data studies from complex lichen-symbiotic communities and provide valuable new information for use in further biochemical and pharmacological studies.

Towards a Systems Biology Approach to Understanding the Lichen Symbiosis: Opportunities and Challenges of Implementing Network Modelling

Lichen associations, a classic model for successful and sustainable interactions between micro-organisms, have been studied for many years. However, there are significant gaps in our understanding about how the lichen symbiosis operates at the molecular level. This review addresses opportunities for expanding current knowledge on signalling and metabolic interplays in the lichen symbiosis using the tools and approaches of systems biology, particularly network modelling. The largely unexplored nature of symbiont recognition and metabolic interdependency in lichens could benefit from applying a holistic approach to understand underlying molecular mechanisms and processes. Together with ‘omics’ approaches, the application of signalling and metabolic network modelling could provide predictive means to gain insights into lichen signalling and metabolic pathways. First, we review the major signalling and recognition modalities in the lichen symbioses studied to date, and then describe how modelling signalling networks could enhance our understanding of symbiont recognition, particularly leveraging omics techniques. Next, we highlight the current state of knowledge on lichen metabolism. We also discuss metabolic network modelling as a tool to simulate flux distribution in lichen metabolic pathways and to analyse the co-dependence between symbionts. This is especially important given the growing number of lichen genomes now available and improved computational tools for reconstructing such models. We highlight the benefits and possible bottlenecks for implementing different types of network models as applied to the study of lichens.

An Overview of Genomics, Phylogenomics and Proteomics Approaches in Ascomycota

Fungi are among the most successful eukaryotes on Earth: they have evolved strategies to survive in the most diverse environments and stressful conditions and have been selected and exploited for multiple aims by humans. The characteristic features intrinsic of Fungi have required evolutionary changes and adaptations at deep molecular levels. Omics approaches, nowadays including genomics, metagenomics, phylogenomics, transcriptomics, metabolomics, and proteomics have enormously advanced the way to understand fungal diversity at diverse taxonomic levels, under changeable conditions and in still under-investigated environments. These approaches can be applied both on environmental communities and on individual organisms, either in nature or in axenic culture and have led the traditional morphology-based fungal systematic to increasingly implement molecular-based approaches. The advent of next-generation sequencing technologies was key to boost advances in fungal genomics and proteomics research. Much effort has also been directed towards the development of methodologies for optimal genomic DNA and protein extraction and separation. To date, the amount of proteomics investigations in Ascomycetes exceeds those carried out in any other fungal group. This is primarily due to the preponderance of their involvement in plant and animal diseases and multiple industrial applications, and therefore the need to understand the biological basis of the infectious process to develop mechanisms for biologic control, as well as to detect key proteins with roles in stress survival. Here we chose to present an overview as much comprehensive as possible of the major advances, mainly of the past decade, in the fields of genomics (including phylogenomics) and proteomics of Ascomycota, focusing particularly on those reporting on opportunistic pathogenic, extremophilic, polyextremotolerant and lichenized fungi. We also present a review of the mostly used genome sequencing technologies and methods for DNA sequence and protein analyses applied so far for fungi.

Complete, high-quality genomes from long-read metagenomic sequencing of two wolf lichen thalli reveals enigmatic genome architecture

Fungal genomes display incredible levels of complexity and diversity, and are exceptional study systems for genome evolution. Here we used the Oxford Nanopore MinION sequencing platform to generate high-quality fungal genomes from complex metagenomic samples of lichen thalli. We sequenced two wolf lichens using one flow cell per sample, generating 17.1 Gbps for Letharia lupina and 14.3 Gbps for Letharia columbiana. The resulting L. lupina genome is one of the most contiguous lichen genomes available to date, with 49.2 Mbp contained on 31 contigs. The L. columbiana genome, while less contiguous, is still relatively high quality, with 52.3 Mbp on a total of 161 contigs. Each thallus for both species contained multiple distinct haplotypes, a phenomenon that has rarely been empirically demonstrated. The Oxford Nanopore sequencing technologies are robust and effective when applied to complex symbioses, and have the potential to fundamentally transform our understanding of fungal genetics.

Understanding transport processes in lichen, Azolla-cyanobacteria, ectomycorrhiza, endomycorrhiza, and rhizobia-legume symbiotic interactions

Intimate interactions between photosynthetic and non-photosynthetic organisms require the orchestrated transfer of ions and metabolites between species. We review recent progress in identifying and characterizing the transport proteins involved in five mutualistic symbiotic interactions: lichens, Azolla–cyanobacteria, ectomycorrhiza, endomycorrhiza, and rhizobia–legumes. This review focuses on transporters for nitrogen and carbon and other solutes exchanged in the interactions. Their predicted functions are evaluated on the basis of their transport mechanism and prevailing transmembrane gradients of H ⁺ and transported substrates. The symbiotic interactions are presented in the assumed order from oldest to most recently evolved.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata

BACKGROUND

Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution.

RESULTS

A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting.

CONCLUSIONS

The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

A comprehensive catalogue of polyketide synthase gene clusters in lichenizing fungi

Lichens are fungi that form symbiotic partnerships with algae. Although lichens produce diverse polyketides, difficulties in establishing and maintaining lichen cultures have prohibited detailed studies of their biosynthetic pathways. Creative, albeit non-definitive, methods have been developed to assign function to biosynthetic gene clusters in lieu of techniques such as gene knockout and heterologous expressions that are commonly applied to easily cultivatable organisms. We review a total of 81 completely sequenced polyketide synthase (PKS) genes from lichenizing fungi, comprising to our best efforts all complete and reported PKS genes in lichenizing fungi to date. This review provides an overview of the approaches used to locate and sequence PKS genes in lichen genomes, current approaches to assign function to lichen PKS gene clusters, and what polyketides are proposed to be biosynthesized by these PKS. We conclude with remarks on prospects for genomics-based natural products discovery in lichens. We hope that this review will serve as a guide to ongoing research efforts on polyketide biosynthesis in lichenizing fungi.

Identification of a putative polyketide synthase gene involved in usnic acid biosynthesis in the lichen Nephromopsis pallescens

Usnic acid is a unique polyketide produced by lichens. To characterize usnic acid biosynthesis, the transcriptome of the usnic-acid-producing lichen-forming fungus Nephromopsis pallescens was sequenced using Illumina NextSeq technology. Seven complete non-reducing polyketide synthase genes and nine highly-reducing polyketide synthase genes were obtained through transcriptome analysis. Gene expression results obtained by qPCR and usnic acid detection with LCMS-IT-TOF showed that Nppks7 is probably involved in usnic acid biosynthesis in N. pallescens. Nppks7 is a non-reducing polyketide synthase with a MeT domain that also possesses beta-ketoacyl-ACP synthase, acyl transferase, product template, acyl carrier protein, C-methyltransferase, and Claisen cyclase domains. Phylogenetic analysis shows that Nppks7and other polyketide synthases from lichens form a unique monophyletic clade. Taken together, our data indicate that Nppks7 is a novel PKS in N. pallescens that is likely involved in usnic acid biosynthesis.

Draft Genome Sequence of the Lichen-Forming Fungus Ramalina intermedia Strain YAF0013

Here, we report a draft genome sequence of Ramalina intermedia strain YAF0013. The functional annotation of R. intermedia provides important information related to its ability to produce secondary metabolites. The genome sequence reported here builds the basis for further genome mining.

Modification and functional adaptation of the MBF1 gene family in the lichenized fungus Endocarpon pusillum under environmental stress

The multiprotein-bridging factor 1 (MBF1) gene family is well known in archaea, non-lichenized fungi, plants, and animals, and contains stress tolerance-related genes. Here, we identified four unique mbf1 genes in the lichenized fungi Endocarpon spp. A phylogenetic analysis based on protein sequences showed the translated MBF1 proteins of the newly isolated mbf1 genes formed a monophyletic clade different from other lichen-forming fungi and Ascomycota groups in general, which may reflect the evolution of the biological functions of MBF1s. In contrast to the lack of function reported in yeast, we determined that lysine¹¹⁴ in the deduced Endocarpon pusillum MBF1 protein (EpMBF1) had a specific function that was triggered by environmental stress. Further, the Endocarpon-specific C-terminus of EpMBF1 was found to participate in stress tolerance. Epmbf1 was induced by a number of abiotic stresses in E. pusillum and transgenic yeast, and its stress-resistant ability was stronger than that of the yeast mbf1. These findings highlight the evolution and function of EpMBF1 and provide new insights into the co-evolution hypothesis of MBF1 and TATA-box-binding proteins.

Lichenized Fungi and the Evolution of Symbiotic Organization

Lichen symbioses comprise a fascinating relationship between algae and fungi. The lichen symbiotic lifestyle evolved early in the evolution of ascomycetes and is also known from a few basidiomycetes. The ascomycete lineages have diversified in the lichenized stage to give rise to a tremendous variety of morphologies. Their thalli are often internally complex and stratified for optimized integration of algal and fungal metabolisms. Thalli are frequently colonized by specific nonlichenized fungi and occasionally also by other lichens. Microscopy has revealed various ways these fungi interact with their hosts. Besides the morphologically recognizable diversity of the lichen mycobionts and lichenicolous (lichen-inhabiting) fungi, many other microorganisms including other fungi and bacterial communities are now detected in lichens by culture-dependent and culture-independent approaches. The application of multi-omics approaches, refined microscopic techniques, and physiological studies has added to our knowledge of lichens, not only about the taxa involved in the lichen interactions, but also about their functions.

Identification of the Sfp-Type PPTase EppA from the Lichenized Fungus Evernia prunastri

In the last decades, natural products from lichens have gained more interest for pharmaceutical application due to the broad range of their biological activity. However, isolation of the compounds of interest directly from the lichen is neither feasible nor sustainable due to slow growth of many lichens. In order to develop a pipeline for heterologous expression of lichen biosynthesis gene clusters and thus the sustainable production of their bioactive compounds we have identified and characterized the phosphopantheteinyl transferase (PPTase) EppA from the lichen Evernia prunastri. The Sfp-type PPTase EppA was functionally characterized through heterologous expression in E. coli using the production of the blue pigment indigoidine as readout and by complementation of a lys5 deletion in S. cerevisiae.