Next Article in Journal
What Role Might Non-Mating Receptors Play in Schizophyllum commune?
Next Article in Special Issue
Production and Activity of Cristazarin in the Lichen-Forming Fungus Cladonia metacorallifera
Previous Article in Journal
Macrofungi as a Nutraceutical Source: Promising Bioactive Compounds and Market Value
Previous Article in Special Issue
Effect of Isolation Conditions on Diversity of Endolichenic Fungal Communities from a Foliose Lichen, Parmotrema tinctorum
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 7
Issue 5
10.3390/jof7050398
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis Identifies a Gene Cluster for the Biosynthesis of Biruloquinone, a Rare Phenanthraquinone, in a Lichen-Forming Fungus Cladonia macilenta

by
Wonyong Kim
1,*,
Min-Hye Jeong
1,
Sung-Hwan Yun
2 and
Jae-Seoun Hur
1,*
1
Korean Lichen Research Institute, Sunchon National University, Suncheon 57922, Korea
2
Department of Medical Sciences, Soonchunhyang University, Asan 31538, Korea
*
Authors to whom correspondence should be addressed.
J. Fungi 2021, 7(5), 398; https://0-doi-org.brum.beds.ac.uk/10.3390/jof7050398
Submission received: 23 April 2021 / Revised: 3 May 2021 / Accepted: 13 May 2021 / Published: 20 May 2021
(This article belongs to the Special Issue Application of Lichen-Forming Fungi for Industrial Use)
Browse Figures
Versions Notes

Abstract

:
Lichens are prolific producers of natural products of polyketide origin. We previously described a culture of lichen-forming fungus (LFF) Cladonia macilenta that produces biruloquinone, a purple pigment that is a phenanthraquinone rarely found in nature. However, there was no genetic information on the biosynthesis of biruloquinone. To identify a biosynthetic gene cluster for biruloquinone, we mined polyketide synthase (PKS) genes from the genome sequence of a LFF isolated from thalli of C. macilenta. The 38 PKS in C. macilenta are highly diverse, many of which form phylogenetic clades with PKS previously characterized in non-lichenized fungi. We compared transcriptional profiles of the 38 PKS genes in two chemotypic variants, one producing biruloquinone and the other producing no appreciable metabolite in vitro. We identified a PKS gene (hereafter PKS21) that was highly upregulated in the LFF that produces biruloquinone. The boundaries of a putative biruloquinone gene cluster were demarcated by co-expression patterns of six clustered genes, including the PKS21. Biruloquinone gene clusters exhibited a high degree of synteny between related species. In this study we identified a novel PKS family responsible for the biosynthesis of biruloquinone through whole-transcriptome analysis.

1. Introduction

More than a thousand polyketide-derived secondary metabolites (SMs) have been isolated from lichen thalli, many of which are exclusively found in lichens [1,2,3]. Several polyketide-derived SMs, such as barbatic acid, didymic acid, rhodocladonic acid, and thamnolic acid, have been extracted from lichen thalli of Cladonia macilenta Hoffm. [4,5]. In addition, a polyketide, biruloquinone, was discovered in an isolated mycobiont of Cladonia macilenta [6]. Biruloquinone is a phenanthraquinone rarely found in nature. It was first discovered in lichen thalli of Parmelia birulae [7], and later in a plant pathogenic fungus Mycosphaerella rubella [8]. Lichen-forming fungi (LFF) occasionally produce novel phenolic compounds that have not been found in lichen thalli grown in nature [9]. Many of these compounds derived from lichen thalli or LFF culture exhibit interesting biological activities, including anti-oxidative, antibacterial, anti-inflammatory, antifungal, and anticancer properties [3,9,10]. Biruloquinone also has pharmaceutical potential, inhibiting acetylcholinesterase, one of the key enzymes involved in Alzheimer’s disease progression [6].
Fungi produce structurally diverse polyketides, such as the health-threatening aflatoxin [11], the cholesterol-lowering drug lovastatin [12,13], and melanin pigments [14,15,16,17] via iterative type I polyketide synthase (PKS). Fungal PKS can be broadly classified into non-reducing PKS (NR-PKS), partially reducing (PR-PKS) and highly reducing PKS (HR-PKS), based on PKS domain architecture. In addition, there are PKS–nonribosomal peptide synthetase (NRPS) hybrid enzymes that comprise PKS domains found in PR-PKS and a module of NRPS. The ketoacyl synthase (KS) domain of PKS has been considered evolutionarily conserved, and thus it has served as a proxy for the similarity of the entire PKS [18]. Given the highly oxygenated tetracyclic core of biruloquinone, it was predicted to be biosynthesized by an NR-PKS [19]. Thioesterase (TE) domains of NR-PKS have been functionally diverged to have different product release mechanisms that may have led to structural diversity of polyketide-derived compounds. Polyketide product release mechanisms by TE domains could be generally differentiated into five distinct modes [20]: (i) macrolactone formation (intramolecular cyclization via O–C bond formation), (ii) depside or melleolide formation (intermolecular cross-coupling via O–C bond formation), (iii) hydrolytic release, (iv) pyrone formation (intramolecular cyclization via O–C bond formation), and (v) Claisen condensation (intramolecular cyclization via C–C bond formation, followed by deacetylation in some cases).
In our previous study, LFF were isolated from thalli of C. macilenta [19]. Among the isolated LFF, two distinct chemotypes were observed: some produced a copious amount of purple pigment, which was structurally identified as biruloquinone [6], and the others produced no such pigment in axenic culture (Figure 1A; Supplementary Figure S1). To investigate genome-encoded metabolic potentials of C. macilenta, we sequenced the genome of one of the isolated LFF with purple pigmentation [21]. The aims of this study were to explore metabolic potentials of C. macilenta by categorizing PKS and to identify a biosynthetic gene cluster (BGC) responsible for the production of biruloquinone through whole-transcriptome analysis.

2. Materials and Methods

2.1. Genome Annotation and Biosynthetic Gene Cluster Identification

The genome assembly of C. macilenta [21] was annotated using the GenSAS (v.6.0) annotation pipeline [22]. Default settings were used unless otherwise noted. In brief, low-complexity regions and repeats were masked using RepeatModeler (v1.0.11) and RepeatMasker (v4.0.7) (www.repeatmasker.org, accessed on 5 February 2020), setting the DNA source to ‘Fungi’. A masked consensus sequence was generated, on which ab initio gene prediction was performed using the following tools: (i) Augustus (v3.3.1) [23], selecting As. nidulans as a trained organism; (ii) GeneMark-ES (v4.33) [24]; (iii) Genscan (v1.0) [25], using a parameter setting for Human and other vertebrates; and (iv) GlimmerM (v2.5.1) [26], selecting Aspergillus as a trained organism. For homology-based predictions, the NCBI reference transcript and protein databases for Fungi were searched, using (v) BLAST+ (v2.7.1) [27] and (vi) DIAMOND (v0.9.22) [28], respectively. For the consensus gene model prediction using EVidenceModeler (v06-25-2012) [29], the above-mentioned standalone gene predictions were weighted as follows: (i) five, (ii) ten, (iii) one, (iv) one, (v) five, and (vi) five. A total of 10,705 open reading frames were predicted in the current genome annotation. For genome mining, the genome assembly and annotation files of C. macilenta were processed by the antiSMASH program (v5.0+) [30], with a parameter setting: ‘--minimal’.

2.2. Phylogenetic Analysis

KS domain sequences were extracted from 38 PKS (19 NR-PKS and 19 R-PKS) identified in the C. macilenta genome sequence, and 111 PKS (62 NR-PKS and 49 R-PKS) that have been linked to known compounds in other filamentous fungi (Supplementary Table S1), using the online tool NaPDoS [31]. A total of 149 KS domain sequences were aligned using MAFFT (v7.310) [32] with the ‘auto’ setting, and the resulting multiple sequence alignment was trimmed for poorly aligned regions using Gblocks (v0.91b) [33] with the parameter: ‘-b4 = 5’. A maximum likelihood tree was constructed using the RAxML program (v8.2) [34] and annotated by iTOL (v5.7) [35]. Nodal supports were evaluated by 1000 bootstrap replications.

2.3. RNA-Seq Experiment

LFF isolated from thalli of C. macilenta (accession number: KoLRI003786) were obtained from the Korean Lichen and Allied Bioresource Center at Sunchon National University (Suncheon, Korea). A LFF producing biruloquinone (purple strain, accession number KoLRI021765) and a LFF producing no appreciable metabolite (white strain, accession number KoLRI021766) were cultured on malt extract agar media (BD Biosciences, Baltimore, MD, USA) for revitalization. The colonies of purple and white strains were harvested and homogenized to be mycelial solutions in sterilized water, using a bead beater. Approximately 100 mg of shredded mycelia were inoculated into malt extract broth (MEB) (BD Biosciences, Baltimore, MD, USA). The purple and white strains were grown in 100 mL of MEB at 18 °C in an orbital shaker (150 rpm). Fifteen days after incubation, production of biruloquinone was confirmed in the purple strain by a high-performance liquid chromatography (HPLC) analysis (Supplementary Figure S1). Chromatographic separation was achieved on a Prominence Modular LC-20A HPLC instrument (Shimadzu, Kyoto, Japan), as previously described [19]. After chemical profiling, mycelia were harvested, ground to a fine powder in liquid nitrogen, and subjected to total RNA extraction using an easy-spin total RNA extraction kit (iNtRON Biotechnology, Seoul, Korea). cDNA libraries were constructed, using the TruSeq RNA library preparation kit (San Diego, CA, USA), and sequenced on the HiSeq2000 platform at Macrogen Inc. (Seoul, Korea). The quality of raw reads (paired-end, 100 bp) was assessed with the FastQC program (v0.11.3; www.bioinformatics.babraham.ac.uk/projects/fastqc, accessed on 10 November 2019). Then, Illumina adapters were trimmed, and poor-quality reads with a mean quality lower than 15 in 4-bp sliding windows from 3′ ends of reads were trimmed, and trimmed reads shorter than 36 bp were further filtered out using the Trimmomatic program [36]. Filtered reads were mapped to the genome sequence of C. macilenta using the HISAT2 program (v2.1.0) [37]. Gene expression levels in reads per kilobase per million mapped reads (RPKM) values were computed and normalized by effective library size estimated by trimmed mean of M values, using the edgeR R package (v3.26.8) [38].

2.4. Identification of Syntenic Gene Clusters

Genome assemblies of Cladonia species were downloaded from the Joint Genome Institute (JGI) or the National Center for Biotechnology Information (NCBI): C. borealis (NCBI accession: JAFEKC000000000), C. grayi (JGI accession: Cgr/DA2myc/ss v2.0), C. metacorallifera (NCBI accession: GCA_000482085.2), C. rangiferina (NCBI accession: GCA_006146055.1), and C. uncialis (NCBI accession: GCA_002927785.1). The genome assemblies were annotated using the GenSAS pipeline, and BGCs were mined from the five related Cladonia species, using antiSMASH, as described above. To search for homologous BGCs in the six Cladonia spp., we used the BiG-SCAPE program [39] that is useful for investigating the conservation and variation of BGCs in related species. Based on the Jaccard index of domain types, domain sequence similarity, and domain adjacency index, the BiG-SCAPE program calculates a similarity matrix between pairwise combinations of clusters where smaller values indicate greater BGC similarity [39]. A cutoff value of 0.5 was used to identify homologous gene clusters.

2.5. Data Availability

The sequences of the biruloquinone BGC and its flanking genes in C. macilenta (Cma_00864–Cma_00876) have been submitted to GenBank and may be found under accession number MW980728–MW980740. The RNA-seq data for transcriptome analysis of C. macilenta LFF have been deposited in NCBI SRA (PRJNA723447).

3. Results

3.1. The Biosynthesis of Biruloquinone

On the basis of the chemical structure, biruloquinone is proposed to be biosynthesized from an octaketide that is likely biosynthesized by an NR-PKS (Figure 1B) [19]. A proposed folding mechanism of an octaketide product in the catalytic cavity of a yet-unknown NR-PKS suggested that a pyrone ring can be formed when the polyketide intermediate is released by a hydrolytic activity of the TE domain of NR-PKSs, as in the case of the cercosporin biosynthetic pathway [40]. The resulting tetracyclic aromatic ring needs to be further modified by tailoring enzymes to yield biruloquinone. The most parsimonious routes of biruloquinone formation require oxidation of two carbon units and methylation of a hydroxyl group (Figure 1B).

3.2. Phylogenetic Placement of C. macilenta PKS

To search for candidate PKSs responsible for the biosynthesis of biruloquinone, we mined the genome of C. macilenta, using antiSMASH. A total of 38 PKS genes were identified: 19 NR-PKSs and 19 R-PKSs (Supplementary Table S1). Since PKS genes are similar in their deduced amino acid sequences, and PKS domain architecture tends to afford end products similar in chemical structure, we analyzed the 38 C. macilenta PKS in an evolutionary framework, with 62 NR-PKSs and 49 R-PKSs that have been characterized in other filamentous fungi (Supplementary Table S1). Phylogenetic analysis of KS domain sequences revealed two major PKS clades: R-PKS (Figure 2A) and NR-PKS (Figure 2B). The R-PKS clade includes three different PKS types differing in PKS domain architecture: HR-PKS, PR-PKS, and PKS-NRPS hybrid enzymes (Figure 2A). The NR-PKS clade has been grouped into seven [41] or eight major groups [42] by protein sequence similarity and PKS domain architecture. We also identified the eight previously described groups of NR-PKS (I–VIII) (Figure 2B). The PKS genes identified in C. macilenta are highly diverse, forming phylogenetic clades with many different PKSs previously characterized in non-lichenized fungi (Figure 2). For candidate PKSs for the biosynthesis of biruloquinone, we focused on the NR-PKS groups II, III, IV, and V. NR-PKSs in these groups are known to afford polycyclic compounds, such as aflatoxin and melanin, while the NR-PKS groups I, VI, VII, and VIII are known to produce monocyclic compounds, such as orsellinic acid [41,42]. In C. macilenta, eleven PKSs (PKS13, PKS15, PKS16, PKS17, PKS18, PKS21, PKS22, PKS30, PKS31, PKS33, PKS34) were predicted to be involved in biosynthesis of polycyclic compounds.

3.3. Transcriptome Analysis of C. macilenta PKS

To narrow down candidate PKS genes responsible for the biosynthesis of biruloquinone, we performed RNA-seq experiments and compared PKS gene expression levels between the purple strain and the white strain in axenic culture [19]. Among the 38 PKS genes, we observed that an NR-PKS (PKS21) was highly upregulated in the purple strain, while the other PKS genes exhibited similar expression levels in the purple and white strains (Figure 3A). The expression of the PKS21 was 195-fold greater in the purple strain, suggesting that the PKS21 is likely involved in the biosynthesis of biruloquinone. The PKS21 belongs to the NR-PKS group IV and is most closely related to CTB1 (Figure 2B), an NR-PKS involved in the biosynthesis of cercosporin in Cercospora nicotianae [40]. The TE domain of CTB1 is shown to directly catalyze pyrone formation during the ring cyclization process of the cercosporin intermediate [40]. A pyrone formation is also required for the biosynthesis of the biruloquinone intermediate when it is released from PKS enzyme (Figure 1B). The close phylogenetic relationship between CTB1 and the PKS21 further strengthened the putative assignment of the PKS21 to the biosynthesis of biruloquinone.

3.4. Identification of Biruloquinone Biosynthesis Gene Clusters

As the biosynthesis of biruloquinone requires O-methylation at a hydroxyl group and oxidation at two different carbon units (Figure 1B), we examined flanking regions of the PKS21 gene to search for related biosynthetic enzymes. Expectedly, we found two genes, encoding O-methyltransferase (gene ID: Cma_00869) and FAD-dependent monooxygenase (Cma_00870) (Table 1), which may be required for tailoring the biruloquinone intermediate. In addition, there were two major facilitator superfamily (MFS) genes (Cma_00868 and Cma_00872) and a gene encoding a GAL4-type transcription factor (Cma_00867), surrounding the PKS21 gene (Table 1). To determine the borders of a biruloquinone BGC in C. macilenta, we compared expression levels of genes neighboring the PKS21 between the purple strain and white strain. Visualization of RNA-seq reads mapped on the genome of C. macilenta indicated that the six genes related to secondary metabolism (brq1–brg6; Table 1) were coordinately upregulated in the purple strain, while the six genes remained nearly silent in the white strain (Figure 3B). Expression levels of other genes of hypothetical function flanking the six genes (brq1–brg6) were basal and similar between the purple and white strains (Figure 3B,C). Taken together, we were able to delimit the borders of the biruloquinone BGC in C. macilenta.
To investigate whether other Cladonia species also possess a homologous BGC for the biosynthesis of biruloquinone, we mined genome sequences of five Cladonia species, C. borealis, C. grayi, C. metacorallifera, C. rangiferina, and C. uncialis. The BiG-SCAPE program [39] identified a homologous BGC in C. borealis and C. metacorallifera. However, the genomes of Cladonia grayi, C. rangiferina, and C. uncialis did not carry a homologous BGC. The biruloquinone BGCs identified in C. borealis and C. metacorallifera were highly syntenic to the one in C. macilenta, including all the six core genes (brq1–brg6) (Figure 4), suggesting that C. borealis and C. metacorallifera also have genetic potential to produce biruloquinone.

4. Discussion

Despite huge genome-encoded metabolic potentials of lichen-forming fungi [43,44], little is known about biosynthetic genes responsible for lichen metabolite production. Ascribing specific biosynthetic genes to lichen metabolites has been slow, due to a paucity of genetic information and lack of molecular tools for manipulating LFF recalcitrant to genetic transformation. To circumvent these issues, positive correlations between expression levels of biosynthetic genes and production of metabolites of interest have been used as circumstantial evidence for linking biosynthetic genes to lichen metabolites [45,46,47].
We previously described CmaPKS1 as a candidate PKS for the biosynthesis of biruloquinone by analyzing the phylogenetic relationship and PKS domain architecture of CmaPKS1 [19]. After our reclassification of Cladonia PKS genes, CmaPKS1 was henceforth renamed PKS15, forming a phylogenetic clade with GlPKS1 and NsPKS1 that are known to be involved in melanin production in non-lichenized fungi [16,17]. Thus, this putative assignment was erroneous, and, in this study, we provided layers of evidence that the PKS21 (brq5) is responsible for the biosynthesis of biruloquinone: (i) the PKS21 was the sole PKS upregulated in the LFF-producing biruloquinone; (ii) the PKS21 was closely-related to CTB1 that directly catalyzes pyrone formation, an activity required for the biosynthesis of biruloquinone and cercosporin; and (iii) the putative biruloquinone gene cluster (brq1–brg6) included tailoring enzymes required for the biosynthesis of biruloquinone, such as O-methyltransferase and FAD-dependent monooxygenase.
The boundaries of the putative biruloquinone gene cluster were determined by the co-expression of six genes (brq1–brg6). The gene cluster included two MFS-type transporter genes, brq2 and brq6. The brq2 exhibited a dramatic increase (256-fold) in its expression level in the LFF-producing biruloquinone, while the expression of the brq6 exhibited only three-fold greater in the LFF-producing biruloquinone. Although the precise function of the two transporters remains to be investigated, the brq2 may play a major role in pumping the purple pigment, biruloquinone, out to culture media. Fungal biosynthetic gene clusters typically harbor one or more pathway-specific transcription factors that regulate clustered genes in a coordinated manner [48]. The gene cluster also included a gene encoding GAL4-type transcription factor (brq1) that can be found in many biosynthetic gene clusters in filamentous fungi. However, the brq1 lacks a DNA-binding domain and contains only the middle homology region commonly found in typical Zn(II)2Cys6 zinc cluster TFs, such as yeast GAL4. This might be due to misannotation of the brq1 in our current genome assembly, but there was a report that a GAL4-type TF only with the middle homology region regulates a fungal biosynthetic gene cluster [49]. The brq1 exhibited only two-fold increase in expression in the purple strain, and this slight upregulation seems to be sufficient to induce the other clustered genes in the purple strain.
It is currently unknown how chemotypic variation is introduced to produce biruloquinone in C. macilenta LFF isolated from the same lichen thallus. Production of unexpected secondary metabolites (not observed in natural lichens) in axenic culture of LFF have been frequently observed [6,50,51,52,53]. These chemotypic variations may be attributable to phenotypic instability of isolated mycobionts. A significant proportion of biosynthetic gene clusters in fungi remain silenced in normal growth conditions and are subject to tight controls by epigenetic regulators (e.g., histone modifiers) [54,55,56,57,58,59,60,61,62,63]. Thus, it is conceivable that the biruloquinone BGC was derepressed by dysfunctionally regulated histone modifiers due to genomic instability of isolated mycobionts, which otherwise remains silenced during the symbiosis.
In this study, we connected a biosynthetic gene cluster to biruloquinone with high likelihood. Homologous biruloquinone BGCs were also found in related Cladonia species (subsection Erythrocarpae [64]). Given the observations that biruloquinone is produced by different fungal taxa, such as C. fruticulosa (subclade Graciles in the genus Cladonia), Parmelia birulae, and Mycosphaerella rubella [7,8,65], production of rare phenanthraquinone-derived metabolites is more widespread than previously thought. Ascribing the PKS21 (brq5) to biruloquinone helps elucidate the evolutionary origin of this unique PKS family in fungi.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/jof7050398/s1, Supplementary Figure S1: Detection of biruloquinone in lichen-forming fungi (LFF) isolated from Cladonia macilenta, Supplementary Table S1: Known PKS characterized in non-lichenized fungi and PKS identified in Cladonia macilenta.

Author Contributions

Conceptualization, W.K., M.-H.J. and J.-S.H.; investigation, W.K., M.-H.J. and S.-H.Y.; writing—original draft preparation, W.K.; writing—review and editing, W.K. and J.-S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Korea Research Fellowship Program funded by the Ministry of Science, ICT and Future Planning (2018H1D3A1A01074888), and by The Korean National Research Resource Center Program (2017M3A9B8069471), through the National Research Foundation of Korea.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Jung-Ae Ryu and Seol-Hwa Jang for technical support in mycobiont isolation and culture.

Conflicts of Interest

Authors have no potential conflict of interest.

References

  1. Calcott, M.J.; Ackerley, D.F.; Knight, A.; Keyzers, R.A.; Owen, J.G. Secondary metabolism in the lichen symbiosis. Chem. Soc. Rev. 2018, 47, 1730–1760. [Google Scholar] [CrossRef] [PubMed]
  2. Stocker-Wörgötter, E. Metabolic diversity of lichen-forming ascomycetous Fungi: Culturing, polyketide and shikimate metabolite production, and PKS genes. Nat. Prod. Rep. 2008, 25, 188–200. [Google Scholar] [CrossRef] [PubMed]
  3. Stocker-Wörgötter, E.; Cordeiro, L.M.C.; Iacomini, M. Accumulation of potential pharmaceutically relevant lichen metabolites in lichens and cultured lichen symbionts. In Studies in Natural Products Chemistry; Elsevier: Amsterdam, The Netherlands, 2013; Volume 39, pp. 337–380. [Google Scholar]
  4. Culberson, C.F. Supplement to “Chemical and Botanical Guide to Lichen Products”. Bryologist 1970, 73, 177. [Google Scholar] [CrossRef]
  5. Stenroos, S.; Ahti, T. Phylogeny of the genus Cladonia s.lat. (Cladoniaceae, Ascomycetes) inferred from molecular, morphological, and chemical data. Cladistics 2002, 18, 237–278. [Google Scholar] [CrossRef]
  6. Luo, H.; Li, C.; Kim, J.C.; Liu, Y.; Jung, J.S.; Koh, Y.J.; Hur, J.-S. Biruloquinone, an acetylcholinesterase inhibitor produced by lichen-forming Fungus Cladonia macilenta. J. Microbiol. Biotechnol. 2013, 23, 161–166. [Google Scholar] [CrossRef]
  7. Krivoshchekova, O.E.; Stepanenko, L.S.; Mishchenko, N.P.; Denisenko, V.A.; Maksimov, O.B. A study of aromatic metabolites of lichens of the family Parmeliaceae. II. Pigments. Chem. Nat. Compd. 1983, 19, 270–274. [Google Scholar] [CrossRef]
  8. Arnone, A.; Nasini, G.; de Pava, O.V. A reinvestigation of the structure of biruloquinone, a 9,10-phenanthrenequinone isolated from Mycosphaerella rubella. Phytochemistry 1991, 30, 2729–2731. [Google Scholar] [CrossRef]
  9. Yamamoto, Y.; Hara, K.; Kawakami, H.; Komine, M. Lichen substances and their biological activities. In Recent Advances in Lichenology; Upreti, D.K., Divakar, P.K., Shukla, V., Bajpai, R., Eds.; Springer India: New Delhi, India, 2015; pp. 181–199. [Google Scholar]
  10. Müller, K. Pharmaceutically relevant metabolites from lichens. Appl. Microbiol. Biotechnol. 2001, 56, 9–16. [Google Scholar] [CrossRef]
  11. Yu, J.; Chang, P.-K.; Ehrlich, K.C.; Cary, J.W.; Bhatnagar, D.; Cleveland, T.E.; Payne, G.A.; Linz, J.E.; Woloshuk, C.P.; Bennett, J.W. Clustered pathway genes in aflatoxin biosynthesis. Appl. Environ. Microbiol. 2004, 70, 1253–1262. [Google Scholar] [CrossRef] [Green Version]
  12. Hendrickson, L.; Davis, C.R.; Roach, C.; Nguyen, D.K.; Aldrich, T.; McAda, P.C.; Reeves, C.D. Lovastatin biosynthesis in Aspergillus terreus: Characterization of blocked mutants, enzyme activities and a multifunctional polyketide synthase gene. Chem. Biol. 1999, 6, 429–439. [Google Scholar] [CrossRef] [Green Version]
  13. Meehan, M.J.; Xie, X.; Zhao, X.; Xu, W.; Tang, Y.; Dorrestein, P.C. FT-ICR-MS characterization of intermediates in the biosynthesis of the α-methylbutyrate side chain of lovastatin by the 277 KDa polyketide synthase LovF. Biochemistry 2011, 50, 287–299. [Google Scholar] [CrossRef] [Green Version]
  14. Moriwaki, A.; Kihara, J.; Kobayashi, T.; Tokunaga, T.; Arase, S.; Honda, Y. Insertional mutagenesis and characterization of a polyketide synthase gene (PKS1) required for melanin biosynthesis in Bipolaris oryzae. FEMS Microbiol. Lett. 2004, 238, 1–8. [Google Scholar]
  15. Cheng, Q.; Kinney, K.A.; Whitman, C.P.; Szaniszlo, P.J. Characterization of two polyketide synthase genes in Exophiala lecanii-Corni, a melanized fungus with bioremediation potential. Bioorg. Chem. 2004, 32, 92–108. [Google Scholar] [CrossRef]
  16. Fulton, T.R.; Ibrahim, N.; Losada, M.C.; Grzegorski, D.; Tkacz, J.S. A melanin polyketide synthase (PKS) gene from Nodulisporium sp. that shows homology to the Pks1 gene of Colletotrichum lagenarium. Mol. Gen. Genet. 1999, 262, 714–720. [Google Scholar] [CrossRef]
  17. Zhang, A.; Lu, P.; Dahl-Roshak, A.M.; Paress, P.S.; Kennedy, S.; Tkacz, J.S.; An, Z. Efficient disruption of a polyketide synthase gene (Pks1) required for melanin synthesis through Agrobacterium-mediated transformation of Glarea lozoyensis. Mol. Gen. Genom. 2003, 268, 645–655. [Google Scholar] [CrossRef]
  18. Kroken, S.; Glass, N.L.; Taylor, J.W.; Yoder, O.C.; Turgeon, B.G. Phylogenomic analysis of type I polyketide synthase genes in pathogenic and saprobic ascomycetes. Proc. Natl. Acad. Sci. USA 2003, 100, 15670–15675. [Google Scholar] [CrossRef] [Green Version]
  19. Jeong, M.-H.; Kim, J.A.; Yu, N.H.; Jung, J.S.; Hong, S.G.; Cheong, Y.H.; Hur, J.-S. Isolation and characterization of a non-reducing polyketide synthase gene in Cladonia macilenta. Mycoscience 2015, 56, 49–57. [Google Scholar] [CrossRef]
  20. Liu, L.; Zhang, Z.; Shao, C.-L.; Wang, C.-Y. Analysis of the sequences, structures, and functions of product-releasing enzyme domains in fungal polyketide synthases. Front. Microbiol. 2017, 8, 1685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Park, S.-Y.; Choi, J.; Kim, J.A.; Jeong, M.-H.; Kim, S.; Lee, Y.-H.; Hur, J.-S. Draft genome sequence of Cladonia macilenta KoLRI003786, a lichen-forming fungus producing biruloquinone. Genome Announc. 2013, 1, e00695-13. [Google Scholar] [CrossRef] [Green Version]
  22. Humann, J.L.; Lee, T.; Ficklin, S.; Main, D. Structural and functional annotation of eukaryotic genomes with GenSAS. Methods Mol. Biol. 2019, 1962, 29–51. [Google Scholar]
  23. Stanke, M.; Waack, S. Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics 2003, 19, ii215–ii225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ter-Hovhannisyan, V.; Lomsadze, A.; Chernoff, Y.O.; Borodovsky, M. Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res. 2008, 18, 1979–1990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Richmond, T. Identification of complete gene structures in genomic DNA. Genome Biol. 2000, 1, reports222. [Google Scholar] [CrossRef]
  26. Delcher, A.L.; Harmon, D.; Kasif, S.; White, O.; Salzberg, S.L. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999, 27, 4636–4641. [Google Scholar] [CrossRef]
  27. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [Green Version]
  28. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef]
  29. Haas, B.J.; Salzberg, S.L.; Zhu, W.; Pertea, M.; Allen, J.E.; Orvis, J.; White, O.; Buell, C.R.; Wortman, J.R. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008, 9, R7. [Google Scholar] [CrossRef] [Green Version]
  30. Blin, K.; Shaw, S.; Steinke, K.; Villebro, R.; Ziemert, N.; Lee, S.Y.; Medema, M.H.; Weber, T. AntiSMASH 5.0: Updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 2019, 47, W81–W87. [Google Scholar] [CrossRef] [Green Version]
  31. Ziemert, N.; Podell, S.; Penn, K.; Badger, J.H.; Allen, E.; Jensen, P.R. The natural product domain seeker NaPDoS: A phylogeny based bioinformatic tool to classify secondary metabolite gene diversity. PLoS ONE 2012, 7, e34064. [Google Scholar] [CrossRef] [Green Version]
  32. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  33. Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Stamatakis, A. RAxML Version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed]
  35. Letunic, I.; Bork, P. Interactive Tree of Life (ITOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, W256–W259. [Google Scholar] [CrossRef] [Green Version]
  36. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  37. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-Genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef]
  38. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  39. Navarro-Muñoz, J.C.; Selem-Mojica, N.; Mullowney, M.W.; Kautsar, S.A.; Tryon, J.H.; Parkinson, E.I.; De Los Santos, E.L.C.; Yeong, M.; Cruz-Morales, P.; Abubucker, S.; et al. A computational framework to explore large-scale biosynthetic diversity. Nat. Chem. Biol. 2020, 16, 60–68. [Google Scholar] [CrossRef]
  40. Newman, A.G.; Vagstad, A.L.; Belecki, K.; Scheerer, J.R.; Townsend, C.A. Analysis of the cercosporin polyketide synthase CTB1 reveals a new fungal thioesterase function. Chem. Commun. 2012, 48, 11772–11774. [Google Scholar] [CrossRef] [Green Version]
  41. Ahuja, M.; Chiang, Y.-M.; Chang, S.-L.; Praseuth, M.B.; Entwistle, R.; Sanchez, J.F.; Lo, H.-C.; Yeh, H.-H.; Oakley, B.R.; Wang, C.C.C. Illuminating the diversity of aromatic polyketide synthases in Aspergillus nidulans. J. Am. Chem. Soc. 2012, 134, 8212–8221. [Google Scholar] [CrossRef] [Green Version]
  42. Liu, L.; Zhang, Z.; Shao, C.-L.; Wang, J.-L.; Bai, H.; Wang, C.-Y. Bioinformatical analysis of the sequences, structures and functions of fungal polyketide synthase product template domains. Sci. Rep. 2015, 5, 10463. [Google Scholar] [CrossRef] [Green Version]
  43. Calchera, A.; Dal Grande, F.; Bode, H.B.; Schmitt, I. Biosynthetic gene content of the ‘perfume lichens’ Evernia prunastri and Pseudevernia furfuracea. Molecules 2019, 24, 203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bertrand, R.L.; Abdel-Hameed, M.; Sorensen, J.L. Lichen biosynthetic gene clusters. Part I. Genome sequencing reveals a rich biosynthetic potential. J. Nat. Prod. 2018, 81, 723–731. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, Y.; Geng, C.; Yuan, X.; Hua, M.; Tian, F.; Li, C. Identification of a putative polyketide synthase gene involved in usnic acid biosynthesis in the lichen Nephromopsis pallescens. PLoS ONE 2018, 13, e0199110. [Google Scholar] [CrossRef] [PubMed]
  46. Armaleo, D.; Sun, X.; Culberson, C. Insights from the first putative biosynthetic gene cluster for a lichen depside and depsidone. Mycologia 2011, 103, 741–754. [Google Scholar] [CrossRef] [Green Version]
  47. Elshobary, M.E.; Becker, M.G.; Kalichuk, J.L.; Chan, A.C.; Belmonte, M.F.; Piercey-Normore, M.D. Tissue-specific localization of polyketide synthase and other associated genes in the lichen, Cladonia rangiferina, using laser microdissection. Phytochemistry 2018, 156, 142–150. [Google Scholar] [CrossRef]
  48. Yin, W.; Keller, N.P. Transcriptional regulatory elements in fungal secondary metabolism. J. Microbiol. 2011, 49, 329–339. [Google Scholar] [CrossRef] [Green Version]
  49. Kim, W.; Park, J.-J.; Gang, D.R.; Peever, T.L.; Chen, W. A novel type pathway-specific regulator and dynamic genome environments of a solanapyrone biosynthesis gene cluster in the fungus Ascochyta rabiei. Eukaryot. Cell 2015, 14, 1102–1113. [Google Scholar] [CrossRef] [Green Version]
  50. Molina, M.C.; Crespo, A.; Vicente, C.; Elix, J.A. Differences in the composition of phenolics and fatty acids of cultured mycobiont and thallus of Physconia distorta. Plant Physiol. Biochem. 2003, 41, 175–180. [Google Scholar] [CrossRef]
  51. Brunauer, G.; Hager, A.; Grube, M.; Türk, R.; Stocker-Wörgötter, E. Alterations in secondary metabolism of aposymbiotically grown mycobionts of Xanthoria elegans and cultured resynthesis stages. Plant Physiol. Biochem. 2007, 45, 146–151. [Google Scholar] [CrossRef]
  52. Ejiri, H.; Sankawa, U.; Shibata, S. Graciliformin and its acetates in Cladonia graciliformis. Phytochemistry 1975, 14, 277–279. [Google Scholar] [CrossRef]
  53. Yamamoto, Y.; Matsubara, H.; Kinoshita, Y.; Kinoshita, K.; Koyama, K.; Takahashi, K.; Ahmadjiam, V.; Kurokawa, T.; Yoshimura, I. Naphthazarin derivatives from cultures of the lichen Cladonia cristatella. Phytochemistry 1996, 43, 1239–1242. [Google Scholar] [CrossRef]
  54. Dallery, J.; Adelin, É.; Le Goff, G.; Pigné, S.; Auger, A.; Ouazzani, J.; O’Connell, R.J. H3K4 trimethylation by CclA regulates pathogenicity and the production of three families of terpenoid secondary metabolites in Colletotrichum higginsianum. Mol. Plant Pathol. 2019, 20, 831–842. [Google Scholar] [CrossRef] [Green Version]
  55. Bachleitner, S.; Sørensen, J.L.; Gacek-Matthews, A.; Sulyok, M.; Studt, L.; Strauss, J. Evidence of a demethylase-independent role for the H3K4-specific histone demethylases in Aspergillus nidulans and Fusarium graminearum secondary metabolism. Front. Microbiol. 2019, 10, 1759. [Google Scholar] [CrossRef] [Green Version]
  56. Studt, L.; Janevska, S.; Arndt, B.; Boedi, S.; Sulyok, M.; Humpf, H.-U.; Tudzynski, B.; Strauss, J. Lack of the COMPASS component Ccl1 reduces H3K4 trimethylation levels and affects transcription of secondary metabolite genes in two plant–pathogenic Fusarium species. Front. Microbiol. 2017, 7, 2144. [Google Scholar] [CrossRef] [Green Version]
  57. Studt, L.; Rösler, S.M.; Burkhardt, I.; Arndt, B.; Freitag, M.; Humpf, H.-U.; Dickschat, J.S.; Tudzynski, B. Knock-down of the methyltransferase Kmt6 relieves H3K27me3 and results in induction of cryptic and otherwise silent secondary metabolite gene clusters in Fusarium fujikuroi: Relief of H3K27me3 induces silent SM gene clusters. Environ. Microbiol. 2016, 18, 4037–4054. [Google Scholar] [CrossRef] [Green Version]
  58. Gacek-Matthews, A.; Berger, H.; Sasaki, T.; Wittstein, K.; Gruber, C.; Lewis, Z.A.; Strauss, J. KdmB, a Jumonji histone H3 demethylase, regulates genome-wide H3K4 trimethylation and is required for normal induction of secondary metabolism in Aspergillus nidulans. PLoS Genet. 2016, 12, e1006222. [Google Scholar] [CrossRef]
  59. Gacek-Matthews, A.; Noble, L.M.; Gruber, C.; Berger, H.; Sulyok, M.; Marcos, A.T.; Strauss, J.; Andrianopoulos, A. KdmA, a histone H3 demethylase with bipartite function, differentially regulates primary and secondary metabolism in Aspergillus nidulans: Histone H3K9 demethylase function in growth. Mol. Microbiol. 2015, 96, 839–860. [Google Scholar] [CrossRef] [Green Version]
  60. Soukup, A.A.; Chiang, Y.-M.; Bok, J.W.; Reyes-Dominguez, Y.; Oakley, B.R.; Wang, C.C.C.; Strauss, J.; Keller, N.P. Overexpression of the Aspergillus nidulans histone 4 acetyltransferase EsaA increases activation of secondary metabolite production: H4K12 acetylation indicates SM production potential. Mol. Microbiol. 2012, 86, 314–330. [Google Scholar] [CrossRef] [Green Version]
  61. Reyes-Dominguez, Y.; Bok, J.W.; Berger, H.; Shwab, E.K.; Basheer, A.; Gallmetzer, A.; Scazzocchio, C.; Keller, N.; Strauss, J. Heterochromatic marks are associated with the repression of secondary metabolism clusters in Aspergillus nidulans: Heterochromatin regulation of secondary metabolism. Mol. Microbiol. 2010, 76, 1376–1386. [Google Scholar] [CrossRef] [Green Version]
  62. Bok, J.W.; Chiang, Y.-M.; Szewczyk, E.; Reyes-Dominguez, Y.; Davidson, A.D.; Sanchez, J.F.; Lo, H.-C.; Watanabe, K.; Strauss, J.; Oakley, B.R.; et al. Chromatin-level regulation of biosynthetic gene clusters. Nat. Chem. Biol. 2009, 5, 462–464. [Google Scholar] [CrossRef] [Green Version]
  63. Shwab, E.K.; Bok, J.W.; Tribus, M.; Galehr, J.; Graessle, S.; Keller, N.P. Histone deacetylase activity regulates chemical diversity in Aspergillus. Eukaryot. Cell 2007, 6, 1656–1664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Stenroos, S.; Pino-Bodas, R.; Hyvönen, J.; Lumbsch, H.T.; Ahti, T. Phylogeny of the family Cladoniaceae (Lecanoromycetes, Ascomycota) based on sequences of multiple loci. Cladistics 2019, 35, 351–384. [Google Scholar] [CrossRef]
  65. Arai, Y.; Kinoshita, K.; Usuniwa, Y.; Yamamoto, Y.; Koyama, K.; Takahashi, K. Search for the substances of cultured lichen mycobiont—Focusing on Cladonia fruticlosa. Lichenology 2012, 10, 199. [Google Scholar]
Jof 07 00398 g001 550
Figure 1. Production of biruloquinone. (A) Accumulation of biruloquinone (purple pigment) in an axenic culture of Cladonia macilenta LFF (right panel). Left panel shows an axenic culture of C. macilenta LFF devoid of biruloquinone. (B) A proposed polyketide pathway for the biosynthesis of biruloquinone.
Figure 1. Production of biruloquinone. (A) Accumulation of biruloquinone (purple pigment) in an axenic culture of Cladonia macilenta LFF (right panel). Left panel shows an axenic culture of C. macilenta LFF devoid of biruloquinone. (B) A proposed polyketide pathway for the biosynthesis of biruloquinone.
Jof 07 00398 g001
Jof 07 00398 g002 550
Figure 2. Genealogy of fungal PKS. (A) Phylogeny of reducing type (R)-PKS, including highly reducing (HR)-PKS, partially reducing (PR)-PKS, and PKS-non-ribosomal peptide synthetase (NRPS) hybrid enzyme. (B) Phylogeny of non-reducing type (NR)-PKS. Roman numerals in red indicate eight phylogenetic groups of fungal NR-PKS. Shaded boxes indicate phylogenetic clades (bootstrap values greater than 90) including a C. macilenta PKS and one or more PKSs characterized in other filamentous fungi. End products biosynthesized by PKSs in other filamentous fungi are denoted in the phylogenetic clades (see Supplementary Table S1). The numbers at the internal nodes indicate bootstrap values greater than 70 from 1000 bootstrap replications. Branch lengths are proportional to the inferred amount of evolutionary change, and the scale represents 1.0 amino acid sequence substitutions per site.
Figure 2. Genealogy of fungal PKS. (A) Phylogeny of reducing type (R)-PKS, including highly reducing (HR)-PKS, partially reducing (PR)-PKS, and PKS-non-ribosomal peptide synthetase (NRPS) hybrid enzyme. (B) Phylogeny of non-reducing type (NR)-PKS. Roman numerals in red indicate eight phylogenetic groups of fungal NR-PKS. Shaded boxes indicate phylogenetic clades (bootstrap values greater than 90) including a C. macilenta PKS and one or more PKSs characterized in other filamentous fungi. End products biosynthesized by PKSs in other filamentous fungi are denoted in the phylogenetic clades (see Supplementary Table S1). The numbers at the internal nodes indicate bootstrap values greater than 70 from 1000 bootstrap replications. Branch lengths are proportional to the inferred amount of evolutionary change, and the scale represents 1.0 amino acid sequence substitutions per site.
Jof 07 00398 g002
Jof 07 00398 g003 550
Figure 3. Identification of biruloquinone PKS and associated biosynthetic gene cluster. (A) Expression values (Log2-transformed (RPKM+1)) are shown as heat maps for 38 PKS genes in a LFF that lacks biruloquinone production (white) and a LFF that produces biruloquinone (purple) in axenic culture. (B) Per-base coverage of RNA-seq reads was plotted for a genomic locus surrounding the PKS21 (brq5; gene ID: Cma_00871). RNA-seq reads mapped on the C. macilenta reference genome were subsampled to 60 million reads for visual comparison of expression levels between white (lower panel) and purple (upper panel). Arrows on the x-axis indicate genes (brq1–brg6) and flanking genes. The numbers on the y-axis are per-base coverage. (C) Log2-transformed fold change of RPKM values in purple compared to white. RPKM, reads per kilobase per million mapped reads.
Figure 3. Identification of biruloquinone PKS and associated biosynthetic gene cluster. (A) Expression values (Log2-transformed (RPKM+1)) are shown as heat maps for 38 PKS genes in a LFF that lacks biruloquinone production (white) and a LFF that produces biruloquinone (purple) in axenic culture. (B) Per-base coverage of RNA-seq reads was plotted for a genomic locus surrounding the PKS21 (brq5; gene ID: Cma_00871). RNA-seq reads mapped on the C. macilenta reference genome were subsampled to 60 million reads for visual comparison of expression levels between white (lower panel) and purple (upper panel). Arrows on the x-axis indicate genes (brq1–brg6) and flanking genes. The numbers on the y-axis are per-base coverage. (C) Log2-transformed fold change of RPKM values in purple compared to white. RPKM, reads per kilobase per million mapped reads.
Jof 07 00398 g003
Jof 07 00398 g004 550
Figure 4. Synteny of biruloquinone BGCs in Cladonia species. Arrows indicate open reading frames (ORFs) found in biruloquinone BGC (brq1–brg6) loci in the three Cladonia species. Pfam domains related to secondary metabolism are color-coded in predicted ORFs (see Table 1).
Figure 4. Synteny of biruloquinone BGCs in Cladonia species. Arrows indicate open reading frames (ORFs) found in biruloquinone BGC (brq1–brg6) loci in the three Cladonia species. Pfam domains related to secondary metabolism are color-coded in predicted ORFs (see Table 1).
Jof 07 00398 g004
Table
Table 1. The biruloquinone biosynthetic gene cluster in Cladonia macilenta.
Table 1. The biruloquinone biosynthetic gene cluster in Cladonia macilenta.
ORF a
(Cma)		Size
(aa)		BLASTP Homolog
(Accession)		%
Identity
(Coverage)		Conserved DomainE-Value
00864325extradiol dioxygenase (XP_037155234)71 (100)Memo-like protein (Pfam01875)6 × 10−77
00865393hypothetical protein (CAF9912135)26 (31)(not detected)-
00866737hypothetical protein (CAF9931462)60 (96)(not detected)-
00867699GAL4-type transcription factor (PMD28964), brq156 (80)fungal-specific transcription factor (Pfam04082)1 × 10−14
00868452MFS transporter (OCK96312), brq275 (99)Major Facilitator Superfamily (Pfam07690)1 × 10−21
00869425O-methyltransferase (KAA6412585), brq374 (99)O-methyltransferase (Pfam00891)4 × 10−14
00870430FAD-dependent monooxygenase (KAA6409270), brq472 (99)NAD(P)-binding Rossmann-like domain (Pfam13450)2 × 10−7
008711971polyketide synthase (QIX11496), brq594 (100)ketoacyl synthase (Pfam00109)8 × 10−88
00872590MFS transporter (KAA6410135), brq671 (96)Major Facilitator Superfamily (Pfam07690)3 × 10−32
00873167hypothetical protein (XP_018191177)30 (60)(not detected)-
00874191hypothetical protein (CAF9930055)86 (100)CS domain (Pfam04969)5 × 10−19
00875159hypothetical protein (XP_037155314)87 (98)Polysaccharide biosynthesis (Pfam04669)7 × 10−38
00876970ABC efflux pump (SLM35520)50 (99)ABC transporter transmembrane region (Pfam00664)3 × 10−35
a ORF: open reading frames; genes for the biruloquinone biosynthetic gene cluster are highlighted in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, W.; Jeong, M.-H.; Yun, S.-H.; Hur, J.-S. Transcriptome Analysis Identifies a Gene Cluster for the Biosynthesis of Biruloquinone, a Rare Phenanthraquinone, in a Lichen-Forming Fungus Cladonia macilenta. J. Fungi 2021, 7, 398. https://0-doi-org.brum.beds.ac.uk/10.3390/jof7050398

AMA Style

Kim W, Jeong M-H, Yun S-H, Hur J-S. Transcriptome Analysis Identifies a Gene Cluster for the Biosynthesis of Biruloquinone, a Rare Phenanthraquinone, in a Lichen-Forming Fungus Cladonia macilenta. Journal of Fungi. 2021; 7(5):398. https://0-doi-org.brum.beds.ac.uk/10.3390/jof7050398

Chicago/Turabian Style

Kim, Wonyong, Min-Hye Jeong, Sung-Hwan Yun, and Jae-Seoun Hur. 2021. "Transcriptome Analysis Identifies a Gene Cluster for the Biosynthesis of Biruloquinone, a Rare Phenanthraquinone, in a Lichen-Forming Fungus Cladonia macilenta" Journal of Fungi 7, no. 5: 398. https://0-doi-org.brum.beds.ac.uk/10.3390/jof7050398

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop