Genome-wide analysis of biosynthetic gene cluster reveals correlated gene loss with absence of usnic acid in lichen-forming fungi

Abstract

Lichen-forming fungi are known to produce a large number of secondary metabolites. Some metabolites are deposited in the cortical layer of the lichen thallus where they exert important ecological functions, such as UV filtering. The fact that closely related lineages of lichen-forming fungi can differ in cortical chemistry suggests that natural product biosynthesis in lichens can evolve independent from phylogenetic constraints. Usnic acid is one of the major cortical pigments in lichens. Here we used a comparative genomic approach on 46 lichen-forming fungal species of the Lecanoromycetes to elucidate the biosynthetic gene content and evolution of the gene cluster putatively responsible for the biosynthesis of usnic acid. Whole-genome sequences were gathered from taxa belonging to different orders and families of Lecanoromycetes, where Parmeliaceae is the most well-represented taxon, and analyzed with a variety of genomic tools. The highest number of biosynthetic gene clusters was found in Evernia prunastri, Pannoparmelia angustata, and Parmotrema austrosinense, respectively, and lowest in Canoparmelia nairobiensis, Bulbothrix sensibilis, and Hypotrachyna scytodes. We found that all studied species producing usnic acid contain the putative usnic acid biosynthetic gene cluster, whereas the cluster was absent in all genomes of species lacking usnic acid. The absence of the gene cluster was supported by an additional unsuccessful search for ß-ketoacylsynthase, the most conserved domain of the gene cluster, in the genomes of species lacking usnic acid. The domain architecture of this PKS cluster-homologous to the already known usnic acid PKS cluster (MPAS) and CYT450 (MPAO)-varies within the studied species, whereas the gene arrangement is highly similar in closely related taxa. We hypothesize that the ancestor of these lichen-forming fungi contained the putative usnic acid producing PKS cluster and that the gene cluster was lost repeatedly during the evolution of these groups. Our study provides insight into the genomic adaptations to the evolutionary success of these lichen-forming fungal species and sets a baseline for further exploration of biosynthetic gene content and its evolutionary significance.