A comprehensive catalogue of polyketide synthase gene clusters in lichenizing fungi

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.

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.

Evolution-Informed Discovery of the Naphthalenone Biosynthetic Pathway in Fungi

Fungi produce a wide diversity of secondary metabolites with interesting biological activities for the health, industrial, and agricultural sectors. While fungal genomes have revealed an unexpectedly high number of biosynthetic pathways that far exceeds the number of known molecules, accessing and characterizing this hidden diversity remain highly challenging. Here, we applied a combined phylogenetic dereplication and comparative genomics strategy to explore eight lichenizing fungi. The determination of the evolutionary relationships of aromatic polyketide pathways resulted in the identification of an uncharacterized biosynthetic pathway that is conserved in distant fungal lineages. The heterologous expression of the homologue from Aspergillus parvulus linked this pathway to naphthalenone compounds, which were detected in cultures when the pathway was expressed. Our unbiased and rational strategy generated evolutionary knowledge that ultimately linked biosynthetic genes to naphthalenone polyketides. Applied to many more genomes, this approach can unlock the full exploitation of the fungal kingdom for molecule discovery.

Beyond Soil-Dwelling Actinobacteria: Fantastic Antibiotics and Where to Find Them

Bacterial secondary metabolites represent an invaluable source of bioactive molecules for the pharmaceutical and agrochemical industries. Although screening campaigns for the discovery of new compounds have traditionally been strongly biased towards the study of soil-dwelling Actinobacteria, the current antibiotic resistance and discovery crisis has brought a considerable amount of attention to the study of previously neglected bacterial sources of secondary metabolites. The development and application of new screening, sequencing, genetic manipulation, cultivation and bioinformatic techniques have revealed several other groups of bacteria as producers of striking chemical novelty. Biosynthetic machineries evolved from independent taxonomic origins and under completely different ecological requirements and selective pressures are responsible for these structural innovations. In this review, we summarize the most important discoveries related to secondary metabolites from alternative bacterial sources, trying to provide the reader with a broad perspective on how technical novelties have facilitated the access to the bacterial metabolic dark matter.

Deletion of a Rare Fungal PKS CgPKS11 Promotes Chaetoglobosin A Biosynthesis, Yet Defers the Growth and Development of Chaetomium globosum

We previously reported that chaetoglobosin A (ChA) exhibits a great potential in the biocontrol of nematodes and pathogenic fungi. To improve the production of ChA, a CRISPR-Cas9 system was created and applied for eliminating potential competitive polyketide products. One of the polyketide synthase encoding genes, Cgpks11, which is putatively involved in the biosynthesis of chaetoglocin A, was disrupted. Cgpks11 deletion led to the overexpression of the CgcheA gene cluster, which is responsible for ChA biosynthesis, and a 1.6-fold increase of ChA. Transcription of pks-1, a melanin PKS, was simultaneously upregulated. Conversely, the transcription of genes for chaetoglocin A biosynthesis, e.g., CHGG_10646 and CHGG_10649, were significantly downregulated. The deletion also led to growth retardation and seriously impaired ascospore development. This study found a novel regulatory means on the biosynthesis of ChA by CgPKS11. CgPKS11 affects chaetoglobosin A biosynthesis, growth, and development in Chaetomium globosum.

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.

The Architecture of Metabolism Maximizes Biosynthetic Diversity in the Largest Class of Fungi

Ecological diversity in fungi is largely defined by metabolic traits, including the ability to produce secondary or “specialized” metabolites (SMs) that mediate interactions with other organisms. Fungal SM pathways are frequently encoded in biosynthetic gene clusters (BGCs), which facilitate the identification and characterization of metabolic pathways. Variation in BGC composition reflects the diversity of their SM products. Recent studies have documented surprising diversity of BGC repertoires among isolates of the same fungal species, yet little is known about how this population-level variation is inherited across macroevolutionary timescales. Here, we applied a novel linkage-based algorithm to reveal previously unexplored dimensions of diversity in BGC composition, distribution, and repertoire across 101 species of Dothideomycetes, which are considered the most phylogenetically diverse class of fungi and known to produce many SMs. We predicted both complementary and overlapping sets of clustered genes compared with existing methods and identified novel gene pairs that associate with known secondary metabolite genes. We found that variation among sets of BGCs in individual genomes is due to nonoverlapping BGC combinations and that several BGCs have biased ecological distributions, consistent with niche-specific selection. We observed that total BGC diversity scales linearly with increasing repertoire size, suggesting that secondary metabolites have little structural redundancy in individual fungi. We project that there is substantial unsampled BGC diversity across specific families of Dothideomycetes, which will provide a roadmap for future sampling efforts. Our approach and findings lend new insight into how BGC diversity is generated and maintained across an entire fungal taxonomic class.

The Mosaic Architecture of NRPS-PKS in the Arbuscular Mycorrhizal Fungus Gigaspora margarita Shows a Domain With Bacterial Signature

As obligate biotrophic symbionts, arbuscular mycorrhizal fungi (AMF) live in association with most land plants. Among them, Gigaspora margarita has been deeply investigated because of its peculiar features, i.e., the presence of an intracellular microbiota with endobacteria and viruses. The genome sequencing of this fungus revealed the presence of some hybrid non-ribosomal peptide synthases-polyketide synthases (NRPS-PKS) that have been rarely identified in AMF. The aim of this study is to describe the architecture of these NRPS-PKS sequences and to understand whether they are present in other fungal taxa related to G. margarita. A phylogenetic analysis shows that the ketoacyl synthase (KS) domain of one G. margarita NRPS-PKS clusters with prokaryotic sequences. Since horizontal gene transfer (HGT) has often been advocated as a relevant evolutionary mechanism for the spread of secondary metabolite genes, we hypothesized that a similar event could have interested the KS domain of the PKS module. The bacterial endosymbiont of G. margarita, Candidatus Glomeribacter gigasporarum (CaGg), was the first candidate as a donor, since it possesses a large biosynthetic cluster involving an NRPS-PKS. However, bioinformatics analyses do not confirm the hypothesis of a direct HGT from the endobacterium to the fungal host: indeed, endobacterial and fungal sequences show a different evolution and potentially different donors. Lastly, by amplifying a NRPS-PKS conserved fragment and mining the sequenced AMF genomes, we demonstrate that, irrespective of the presence of CaGg, G. margarita, and some other related Gigasporaceae possess such a sequence.

IMA Genome - F13: Draft genome sequences of Ambrosiella cleistominuta, Cercospora brassicicola, C. citrullina, Physcia stellaris, and Teratosphaeria pseudoeucalypti

Draft genomes of the fungal species Ambrosiella cleistominuta, Cercospora brassicicola, C. citrullina, Physcia stellaris, and Teratosphaeria pseudoeucalypti are presented. Physcia stellaris is an important lichen forming fungus and Ambrosiella cleistominuta is an ambrosia beetle symbiont. Cercospora brassicicola and C. citrullina are agriculturally relevant plant pathogens that cause leaf-spots in brassicaceous vegetables and cucurbits respectively. Teratosphaeria pseudoeucalypti causes severe leaf blight and defoliation of Eucalyptus trees. These genomes provide a valuable resource for understanding the molecular processes in these economically important fungi.

Advances and Prospects of Phenolic Acids Production, Biorefinery and Analysis

Biotechnological production of phenolic acids is attracting increased interest due to their superior antioxidant activity, as well as other antimicrobial, dietary, and health benefits. As secondary metabolites, primarily found in plants and fungi, they are effective free radical scavengers due to the phenolic group available in their structure. Therefore, phenolic acids are widely utilised by pharmaceutical, food, cosmetic, and chemical industries. A demand for phenolic acids is mostly satisfied by utilising chemically synthesised compounds, with only a low quantity obtained from natural sources. As an alternative to chemical synthesis, environmentally friendly bio-based technologies are necessary for development in large-scale production. One of the most promising sustainable technologies is the utilisation of microbial cell factories for biosynthesis of phenolic acids. In this paper, we perform a systematic comparison of the best known natural sources of phenolic acids. The advances and prospects in the development of microbial cell factories for biosynthesis of these bioactive compounds are discussed in more detail. A special consideration is given to the modern production methods and analytics of phenolic acids.

What Is in Umbilicaria pustulata? A Metagenomic Approach to Reconstruct the Holo-Genome of a Lichen

Lichens are valuable models in symbiosis research and promising sources of biosynthetic genes for biotechnological applications. Most lichenized fungi grow slowly, resist aposymbiotic cultivation, and are poor candidates for experimentation. Obtaining contiguous, high-quality genomes for such symbiotic communities is technically challenging. Here, we present the first assembly of a lichen holo-genome from metagenomic whole-genome shotgun data comprising both PacBio long reads and Illumina short reads. The nuclear genomes of the two primary components of the lichen symbiosis-the fungus Umbilicaria pustulata (33 Mb) and the green alga Trebouxia sp. (53 Mb)-were assembled at contiguities comparable to single-species assemblies. The analysis of the read coverage pattern revealed a relative abundance of fungal to algal nuclei of ∼20:1. Gap-free, circular sequences for all organellar genomes were obtained. The bacterial community is dominated by Acidobacteriaceae and encompasses strains closely related to bacteria isolated from other lichens. Gene set analyses showed no evidence of horizontal gene transfer from algae or bacteria into the fungal genome. Our data suggest a lineage-specific loss of a putative gibberellin-20-oxidase in the fungus, a gene fusion in the fungal mitochondrion, and a relocation of an algal chloroplast gene to the algal nucleus. Major technical obstacles during reconstruction of the holo-genome were coverage differences among individual genomes surpassing three orders of magnitude. Moreover, we show that GC-rich inverted repeats paired with nonrandom sequencing error in PacBio data can result in missing gene predictions. This likely poses a general problem for genome assemblies based on long reads.

Biosynthetic Gene Content of the 'Perfume Lichens' Evernia prunastri and Pseudevernia furfuracea

Lichen-forming fungi produce a vast number of unique natural products with a wide variety of biological activities and human uses. Although lichens have remarkable potential in natural product research and industry, the molecular mechanisms underlying the biosynthesis of lichen metabolites are poorly understood. Here we use genome mining and comparative genomics to assess biosynthetic gene clusters and their putative regulators in the genomes of two lichen-forming fungi, which have substantial commercial value in the perfume industry, Evernia prunastri and Pseudevernia furfuracea. We report a total of 80 biosynthetic gene clusters (polyketide synthases (PKS), non-ribosomal peptide synthetases and terpene synthases) in E. prunastri and 51 in P. furfuracea. We present an in-depth comparison of 11 clusters, which show high homology between the two species. A ketosynthase (KS) phylogeny shows that biosynthetic gene clusters from E. prunastri and P. furfuracea are widespread across the Fungi. The phylogeny includes 15 genomes of lichenized fungi and all fungal PKSs with known functions from the MIBiG database. Phylogenetically closely related KS domains predict not only similar PKS architecture but also similar cluster architecture. Our study highlights the untapped biosynthetic richness of lichen-forming fungi, provides new insights into lichen biosynthetic pathways and facilitates heterologous expression of lichen biosynthetic gene clusters.