Xinhui
Wang
,
Scott A.
Jarmusch
,
Jens C.
Frisvad
and
Thomas O.
Larsen
*
DTU Bioengineering, Technical University of Denmark, DK-2800, Kgs. Lyngby, Denmark. E-mail: tol@bio.dtu.dk
First published on 19th May 2022
Covering: up to the end of 2021
Aspergilli are biosynthetically ‘talented’ micro-organisms and therefore the natural products community has continually been interested in the wealth of biosynthetic gene clusters (BGCs) encoding numerous secondary metabolites related to these fungi. With the rapid increase in sequenced fungal genomes combined with the continuous development of bioinformatics tools such as antiSMASH, linking new structures to unknown BGCs has become much easier when taking retro-biosynthetic considerations into account. On the other hand, in most cases it is not as straightforward to prove proposed biosynthetic pathways due to the lack of implemented genetic tools in a given fungal species. As a result, very few secondary metabolite biosynthetic pathways have been characterized even amongst some of the most well studied Aspergillus spp., section Nigri (black aspergilli). This review will cover all known biosynthetic compound families and their structural diversity known from black aspergilli. We have logically divided this into sub-sections describing major biosynthetic classes (polyketides, non-ribosomal peptides, terpenoids, meroterpenoids and hybrid biosynthesis). Importantly, we will focus the review on metabolites which have been firmly linked to their corresponding BGCs.
To further investigate both inter- and intraspecies genomic variation as well as the chemical potential of species in genus Aspergillus, a large sequencing project including one representative strain of >350 known Aspergillus spp. has been initiated using the Joint Genome Institute (JGI) fungal genome pipeline.4,5 One of the first major outcomes from this initiative was a comprehensive bioinformatic analysis of all detected BGCs from species in section Nigri.6 In line with the above mentioned studies, a major conclusion from this work was that secondary metabolism is unique to individual species. The authors reported a total of 2717 encoding secondary metabolite BGCs from a total of 37 genomes and found that on average, each species contains 8.75 unique BGCs.6 Similar to previous studies,7,8 the authors also identified and linked groups of BGCs, that are functionally closely related and encode the production of the same or very similar secondary metabolites, into 455 biosynthetic gene cluster families.6 This large number is in sharp contrast to the relatively few secondary metabolites that have so far been linked to their corresponding BGCs, indicating the overall potential for discovery in Nigri.
The species in the A. niger group (now called Aspergillus subgenus Circumdati section Nigri) have been treated by Mosseray (1934), and Raper and Fennell (1965).24,25 Similarly, very detailed morphological and physiological taxonomic studies were done by Murakami et al.26–36 with an emphasis on industrial black Kuro koji strains. Later new taxonomies and phylogenies were suggested, based on morphology, physiology, household gene sequences and secondary metabolite profiles.3,6,13,15–17,20,22,37–72
Aspergillus subgenus Circumdati section Nigri presently consists of 5 formal series: ser. Nigri, ser. Carbonarii, ser. Heteromorphi (biseriate species, i.e. species with both metulae and phialides on the globose conidial Aspergillus head (vesicle)), in addition to ser. Japonici (uniseriate species, i.e. species with phialides only on the globose conidial Aspergillus head (vesicle)) and ser. Homomorphi (biseriate).3,46 Ser. Nigri includes 14 species; Ser. Carbonarii includes 4 species; series Heteromorphi includes 2 species; series Japonici includes 15 species, and finally ser. Homomorphi only has one member (Table 1). Many other names have been published and frequently used, especially A. acidus, A. aureus var. pallidus, A. awamori, A. coreanus, A. inuii, A. kawachii and A. nakazawai (all synonyms of A. luchuensis) and A. pulverulentus, A. saitoi and A. saitoi var. kagoshimaensis (synonyms of A. tubingensis), and A. aureus, A. batatae, A. miyakoensis, A. usamii and A. usamii var. shirousamii (all synonyms of A. niger).61,62Aspergillus awamori is a name often used, and was originally a name for a domesticated form of A. niger, later it was regarded as cryptic phylogenetic species in A. niger based on a non-domesticated strain from the plant Welwitschia mirabilis.54 This interpretation of A. awamori was abandoned and the name A. welwitschiae given to this taxon, while the domesticated forms were named A. luchuensis, which is the most commonly used Kuro koji fungus.61,62 The names A. phoenicis and A. ficuum have been rejected in order to conserve the well-known name A. niger.73 The number of accepted species in Aspergillus section Nigri is therefore currently 36. Additionally, many of the species of Aspergillus section Nigri have been genome sequenced (Table 1). Thus there are good opportunities for mining these genomes for BGCs encoding for secondary metabolites.
Series | Species |
---|---|
Nigri | A. brasiliensis, A. chiangmaiensis, A. costaricensis, A. eucalypticola, A. lacticoffeatus, A. luchuensis, A. neoniger, A. niger, A. piperis, A. pseudopiperis, A. pseudotubingensis, A. tubingensis, A. vinaceus, and A. welwitchiae |
Carbonarii | A. carbonarius (=A. fonsecaeus), A. ibericus, A. sclerotiicarbonarius, and A. sclerotioniger |
Heteromorphi | A. ellipticus and A. heteromorphus |
Japonici | A. aculeatinus, A. aculeatus, A. brunneoviolaceus, A. fijiensis, A. floridensis, A. hydei, A. japonicus, A. labruscus, A. oxumae, A. saccharolyticus, A. serratalhadensis, A. siamensis, A. trinidadensis, A. uvarum and A. violaceofuscus |
Homomorphi | A. homomorphus 3,18 |
The species of Aspergillus section Nigri that have been genome sequenced | |
A. niger,6,11,74–80A. luchuensis,75,81–85A. tubingensis,6,75,86,87A. aculeatus,6,88A. aculeatinus in addition to A. brasiliensis, A. brunneoviolaceus, A. carbonarius, A. costaricensis, A. ellipticus, A. eucalypticola, A. fijiensis, A. heteromorphus, A. homomorphus, A. ibericus, A. indologenus, A. japonicus, A. lacticoffeatus, A. neoniger, A. piperis, A. saccharolyticus, A. sclerotiicarbonarius, A. sclerotioniger, A. uvarum, A. vadensis, A. violaceofuscus, and A. welwitschiae6,89 |
Some species in Aspergillus section Nigri can produce mycotoxins, notably ochratoxin A (A. carbonarius, A. lacticoffeatus, A. niger, A. sclerotioniger, A. welwitschiae),40,47,90–96 and fumonisins B2, B4 and B6 (A. lacticoffeatus, A. niger and A. welwitschiae).46,95–100 Secalonic acid, only found in the uniseriate Aspergillus section Nigri species, has also been mentioned as a mycotoxin.101 Oxalic acid is a widespread weak mycotoxin of Aspergillus niger and related species in series Nigri and Carbonarii. It has been linked to kidney failure after growth and in human lungs and transport of calcium oxalate, via the blood to the kidneys.102 Oxalic acid, citric acid and gluconic acid are produced by many species in section Nigri. These natural products should be regarded as secondary metabolites when they are accumulated, secreted and regulated by the general regulator LaeA, in addition to being taxonomically restricted,103–105 but as primary metabolites when they are in flux in the TCA cycle (and therefore not secreted and accumulated and furthermore present in all fungi and many other organisms). Altogether, a large number of secondary metabolites have been reported from Aspergillus niger and similar species,102,106–109 and preliminary genome scale models for Aspergillus niger metabolic networks have been published.110,111
Metabolites | Producing organism | References |
---|---|---|
Fumonisin Py | A. welwitschiae | 112 |
Fumonisin La | A. welwitschiae | 112 |
Ergosterimide B | A. tubingensis YP-2 | 113 |
Demethylincisterol A5 | A. tubingensis YP-2 | 113 |
Sclerolizine | A. sclerotiicarbonarius | 114 |
Emindole SC | A. sclerotiicarbonarius | 114 |
Carbonarins I | A. sclerotiicarbonarius | 114 |
Carbonarins J | A. sclerotiicarbonarius | 114 |
Aculenes A–C | A. aculeatus | 115 |
Okaramine S | A. aculeatus | 115 |
Epi-10,23-dihydro-24,25-dehydroaflavinine | A. aculeatus | 115 |
Acu-dioxomorpholine | A. aculeatus | 115 |
Acucalbistrin A | A. aculeatus | 115 |
Acucalbistrin B | A. aculeatus | 115 |
Dihydrocarolic acid | A. niger | 116 |
Penitricin D | A. niger | 116 |
New diketopiperazine dimer | A. niger | 117 |
Sclerotionigrin A | A. sclerotioniger | 118 |
Sclerotionigrin B | A. sclerotioniger | 118 |
Asperpyrone D | A. tubingensis | 119 |
Methyl protodioscin | A. niger | 120 |
Aurasperone H | A. niger | 121 |
Secalonic acid F and D | A. aculeatus | 122 |
Gibberellic acid | A. niger | 123 |
Aspiperidine oxide | A. indologenus | 124 |
Pre-aurantiamin | A. brunneoviolaceus, A. japonicus and A. aculeatus | 65 |
Homopyrone A | A. homomorphus | 125 |
Homopyrone B | A. homomorphus | 125 |
Name of metabolites | BGC names or key genes | Producer organism | Linkage between metabolite and BGC | Key reference(s) |
---|---|---|---|---|
Naphtho-γ-pyrones | albA (fwnA) | A. niger | The function of albA was confirmed by gene deletion in auxotrophic A. niger strain ATCC 11414 (KB1001) | 126 and 127 |
In another study, the loci involved in spore pigmentation of A. niger were identified. Metabolic profiling of the fawn mutant revealed the production of naphtho-γ-pyrones depends on fwnA (albA) | ||||
aun and bfo (tailoring genes) | A. niger | The tailoring enzymes included in a gene cluster that is separate from PKS core gene albA were identified by gene deletions in A. niger FGSC A1180 (aun) | 128 | |
A. aculeatus | Heterologous expression of aunB and bfoB in A. aculeatus to explore the substrate specificity of the dimerization | |||
Kotanin | ktnS | A. niger | Bioinformatics analysis was used for searching NR-PKS containing BGCs and gene deletion experiments were performed for functional analysis in the NHEJ-deficient A. niger strain FGSC A1180 | 129 |
TAN-1612 | adaA | A. niger | Transcriptional factor was replaced by promoter and the gene cluster was identified using a targeted genome mining approach A. niger ATCC 1015. The function of bifunctional MβL-TE was identified by heterologous pathway reconstitution in Saccharomyces cerevisiae | 130 |
Azanigerones | azaA and azaB | A. niger | The silent gene was activated by overexpression of transcriptional factor (azaR). Two convergent PKSs in A. niger ATCC 1015 were identified by transcriptional analysis and deletion of a key PKS gene (azaB) | 131 |
Nigerpyrones | epaA | A. niger | The epigenetic regulator gcnE was deleted in a highly mutated industrial strain, A. niger FGSC A1279. The function of key gene epaA was confirmed by gene deletion by homologous recombiation and complementary experiments | 132 |
2,4-Dihydroxy-3,5,6-trimethylbenzaldehyde | dtbA | A. niger | The function of dtbA was demonstrated by heterologous expression in A. nidulans. R domain was replaced by a TE release domain in the DtbA to expand the structural diversity | 133 |
Funalenone | phnA | Penicillium herquei | A. niger is a funalenone producer with several homologous genes to those found in P. herquei. The phnA was deleted using homologous recombination and hygromycin as selection marker. The function of phnA was further confirmed by heterologous expression in S. cerevisiae. Point mutation of PhnA was used to investigate the roles of PT and TE/CLC of PhnA | 134 |
Aculinic acid and aculin A | acuA | A. aculeatus | The function of core gene was verified by integration into the IS1 locus in A. nidulans (heterologous expression) and AMA1-based episomal expression and 13C labeling experiments in A. aceleatus | 115 and 135 |
Fumonisins | fum5 | A. niger and A. welwitschiae | The fumonisin biosynthesis have been extensively studied in the different species in Fusarium. The fumonisin production and biosynthesis in black aspergilli were explored mainly with homologous genes analysis | 136–138 |
Homopyrone A and B | AhpA | A. homomorphus | Retro-biosynthesis, genome mining, and gene deletions were used to successfully link to the biosynthetic gene cluster | 125 |
Aculeacin A | EcdA | A. aculeatus | Aculeacin biosynthesis is very closely related to echinocandin B biosynthesis. The gene cluster of aculeacin A was identified by investigating the evolution of chemical diversity in echinocandin lipopeptide metabolites. The proline hydroxylase, Aa-HtyE, was characterized through heterologous expression in E. coli | 139 and 140 |
Malformins | MlfA | A. brasiliensis | The function of core gene mlfA was confirmed by targeted gene deletion in A. brasiliensis and genetic complementation experiments | 107 |
Atrofuranic acid | — | A. niger | The compound was discovered when atromentin synthetases were heterologously expressed in A. niger in order to explore the reaction mechanism of the thioesterase domain | 141 |
Ergosterimide | — | A. niger | This compound was proposed to be a Diels–Alder adduct of ergosteroid and maleimide | 142 |
Taxol | — | A. aculeatinus | The biosynthesis pathway of taxol formed from endophytic fungi was not clear | 143 |
Pyranonigrins E–J | pynA | A. niger | The pyn gene cluster was activated by expressing the transcriptional regulator pynR under control of the arginase (aga) promotor. The biosynthetic pathway was identified by targeted gene deletion in A. niger, heterologous in vivo bioconversion in S. cerevisiae SCKW12 and in vitro assays | 144 and 145 |
Pyranonigrins A, S | pyrA | P. thymicola | The metabolites were first reported from A. niger. The biosynthetic pathway was elucidated by genome mining and heterologous expression in A. nidulans | 146 |
Pyranoviolin A | pyvA | A. violaceofuscus | The BGC was identified via bioinformatics analysis and comparison with the homologous genes to those of pyranonigrins. The PKS–NRPS was confirmed through targeted gene deletion | 147 |
Ochratoxin A | otaA and otaB | A. carbonarius | The NRPS was verified by targeted gene disruption | 148 and 149 |
Acurin A | acrA and acrB | A. aculeatus | The functions of genes were identified through gene deletion experiments | 150 |
Yanuthones | yanA and yanG | A. niger | The PKS was confirmed by heterologous expression in A. nidulans. The biosynthetic pathway was investigated by individually deleting all genes | 151 and 152 |
Aculene A–C | aneB and aneC | A. aculeatus | The biosynthetic pathway was investigated through gene inactivation, feeding experiments and heterologous reconstitution approaches | 153 |
Alkylcitric acids | akcA, akcB and akcD | A. niger | The gene cluster was revealed by genomic data analysis, overexpression of transcriptional factor and targeted gene deletion in A. niger | 154 |
Final product | Gene name | NR-PKS domain organization | |
---|---|---|---|
a The biosynthetic pathway of phenalenone was identified in Penicillium herquei, which can provide the comparative analysis for the biosynthesis of funalenone. b Putatively assigned to a group based on their domain organization and closest characterized homologue. | |||
NR-PKS | Azanigerone A (1) | azaA (group VII) | |
Kotanin (2) | ktnS (group V) | ||
Naphtho-γ-pyrones | AlbA (fwnA) (group III) | ||
TAN-1612 (3) | adaA (group V) | ||
2,4-Dihydroxy-3,5,6-trimethylbenzaldehyde (4), 6-ethyl-2,4-dihydroxy-3,5-dimethylbenzaldehyde (39) | dtbA (group VII) | ||
Nigerpyrone (5) | epaA (group VII) | ||
Phenalenone (42)a | phnA (group III)b | ||
PR-PKS | Aculinic acid (48) | acuA | |
Ochratoxin A (50) | otaA | ||
Yanuthones | yanA | ||
HR-PKS | Fumonisins | fum 1 | |
Azanigerone A (1) | azaB | ||
Homopyrone A and B (70, 71) | ahpA | ||
Pyranonigrins E–J (79–84) | pynA (PKS–NRPS) | ||
Pyranonigrin A (85) | pyrA (PKS–NRPS) | ||
Pyranoviolin A (87) | pyvA (PKS–NRPS) |
Fig. 2 Proposed biosynthetic pathway of DHN-melanin and naphthopyrone dimers in black aspergilli.128 |
In order to explore the naphtho-γ-pyrones biosynthesis in A. niger, albA, a homolog of PKS alb1 originating in A. fumigatus was deleted and led to the abolishment of production of all naphtho-γ-pyrones.126 A contemporary study showed that the naphtho-γ-pyrones polyketides were specifically dependent on fwnA (albA) using UV mutagenesis.127 In addition, both studies showed that albA in A. niger is responsible for the biosynthesis of melanin, which is consistent with the prediction that A. niger DHN-melanin (9) production is similar to the pathway elucidated for A. fumigatus.126,127,158,159 Moreover, a putative ortholog of A. fumigatus ayg1, aygA, was identified in A. niger. When aygA was deleted, melanin production was blocked but naphtho-γ-pyrones were still observed, indicating aygA is not involved in the naphtho-γ-pyrone biosynthetic pathway but only in the biosynthesis of 9 (Fig. 2).126 However, since the downstream biosynthetic genes responsible for naphtho-γ-pyrones are not clustered with the PKS gene albA, researchers were not able to identify the enzymes responsible for methylation and dimerization. It has recently been discovered that dimerization of naphtho-γ-pyrones monomers is catalyzed by cytochrome P450 enzymes in black aspergilli and they are not in the same gene cluster as the polyketide synthase gene.128 To identify the genes involved in naphthopyrone biosynthesis, the dimeric derivatives in five lineages of black aspergilli were re-evaluated by Obermaier and Müller.128 They identified aurasperone B (10) and aurasperone A (11) in A. niger (FGSC A1180), A. tubingensis and A. carbonarius. Bifonsecin B (12) and nigerone (13) were found in A. niger (ATCC 36626) and one strain of A. brasiliensis. The genes encoding tailoring enzymes form a gene cluster that is separate from the polyketide synthase encoding gene. It was therefore proposed that the transcription of these genes is co-regulated with the expression of albA (PKS gene involved in the biosynthesis of YWA1).128 Because YWA1 (7) is also the precursor of the spore pigment melanin, the upregulation most likely coincide with sporulation.158
In order to explore the oxidative enzyme for dimerization from a sporulating culture, publicly available transcriptomes of A. carbonarius obtained from both static and shaken liquid cultures were screened. This identified a CYP enzyme (aunB) clustered with two O-MT genes (aunD, aunE). Other producers of dimeric naphtho-γ-pyrones, such as A. brasiliensis (bfo), A. niger (aun), and A. tubingensis (aun), also had homologous clusters.128 To investigate the function of these enzymes in the biosynthesis of naphthopyrone dimers, an aunB deletion strain was generated in A. niger FGSC A1180 and led to the accumulation of the monomeric compounds fonsecin B (14) and flavasperone (15). This indicated the aunB gene to be responsible for dimerization of the monomeric precursor.128 Additionally, deletion of the O-MT gene aunD resulted in an accumulation of fonsecin (16), which identified that aunD methylates fonsecin. Another methyltransferase gene, aunE, was proposed to be involved in methylation of the hydroxyl group at position 6. To explore the substrate specificity of the CYP enzymes for catalyzing monomeric naphthopyrones into dimers, aunB from A. niger FGSC A1180 and bfoB from A. brasiliensis NRRL 26652 were heterologously expressed in A. aculeatus, which has not been reported to produce dimeric naphthopyrones. AunB was shown to catalyze the regioselective dimerization of fonsecin B (14) and rubrofusarin B (17) to aurasperone B (10) and A (11), respectively, while bfoB catalyze dimerization of bifonsecin B (12) and nigerone (13) (Fig. 2).128 Additional monomeric metabolites, fonsecin (16) and flavasperone (15), were not converted into dimers. Further exploration showed the substrate dependency of the atroposelectivity for the CYP enzymes. (M)-Nigerone was produced by a BfoB catalyzed dimerization of the monomers fonsecin B (14) (enantiomeric excess >95%). However, feeding with rubrofusarin B (17) resulted in an excess of (P)-nigerone. The exact mechanisms leading to atropselectivity is complex and future research will need to explain how the enzymes interact with the various monomeric substrates.128
The BGC responsible for biosynthesis of kotanin was identified using a combination of bioinformatics and gene deletion experiments (Fig. 3);129 searching only for BGCs containing NR-PKS encoding genes as no reduction steps were predicted to be involved in formation of the core dihydroxy coumarin structure from a linear pentaketide. Candidate clusters were further restricted to include only those containing cytochrome P450s, monooxygenases and methyltransferases, all of which would be required for kotanin (2) biosynthesis. KtnS was identified as the NR-PKS responsible for kotanin (2) biosynthesis and further confirmed by targeted gene deletion in the NHEJ-deficient A. niger strain FGSC A1180, resulting in complete loss of 4,7-dihydroxy-5-methylcoumarin (19) production. A functional analysis of the genes involved in the same cluster was also performed by targeted gene deletions, which confirmed the involvement of the O-methyltransferase KtnB in 7-demethylsiderin (20) production and the cytochrome P450 monooxygenase KtnC in kotanin (2) biosynthesis.
Fig. 3 (A) The gene cluster ktn in A. niger and (B) the proposed biosynthetic pathway of P-(+)-kotanin (2).129,163 |
Furthermore, an oxidative phenol coupling catalyzed by the cytochrome P450 monooxygenase KtnC leads to formation of P-(+)-orlandin (18). Docking experiments were performed to gain insights into the mechanism of regio- and stereoselectivity, which revealed that it is likely dependent on substrate orientation in the active site of cytochrome P450 monooxygenase KtnC. Subsequent research also showed that KtnC selectivity catalyzes the intermolecular coupling of 7-demethylsiderin (20), regioselectively to the 8,8′-dimer and stereoselectively to the P-atropisomer.163 Finally, P-(+)-kotanin (2) is then generated following O-methylation of P-(+)-orlandin (Fig. 3).
The PT domain of this compound belongs to group V NR-PKSs, which utilize a trans-acting metallo-β-lactamase-type TE (MβL-TE) instead of a C-terminal release domain; this MβL-TE PT domain would be responsible for catalyzing the C6–C11 regioselective cyclization. However, no MβL-TE or similar enzymes for product release was found in ktn, suggesting a spontaneous (or other) release mechanism of the final polyketide. In addition, the cyclization pattern of the core polyketide structure by ktnS is C2–C7 and not the typical C6–C11 cyclization pattern seen for the group V NR-PKSs. This might be because of the short pentaketide chain or this branch represents a separate group of cyclization encoding PT domains.
Because of the similarities between the ada gene cluster for naphthacenedione 3 and the known apt cluster linked to the biosynthesis of anthraquinone asperthecin, the biosynthesis of these chemically different metabolites were done in parallel via heterologous pathway reconstitution in Saccharomyces cerevisiae.130 This revealed that a novel α-hydroxylation-dependent Claisen cyclization cascade, which involves a flavin-dependent monooxygenase that hydroxylates the α-carbon of an acyl carrier protein-bound polyketide and a bifunctional metallo-β-lactamase-type thioesterase (MβL-TE). The bifunctional MβL-TE catalyzes the cyclization of the fourth ring to afford the naphthacenedione scaffold upon α-hydroxylation, whereas it performs hydrolytic release of an anthracenone product (21) in the absence of α-hydroxylation. The MβL-TE is the first example of a thioesterase in polyketide biosynthesis that catalyzes a Claisen-like condensation without an α/β hydrolase fold. These mechanistic features should be general to the biosynthesis of tetracyclic naphthacenedione compounds in fungi.130
The pathway begins with the production of decaketide (C20) backbone 22 synthesized and cyclized by AdaA, followed by the hydroxylation of the C-2 position by the monoxygenase Ada C to produce 23. Next, the final cyclization in 24 production is mediated by the MβL-TE, AdaB. Finally, AdaD then performs the C-9 O-methylation to complete the biosynthesis of 3 (Fig. 4).
Fig. 4 (A) Organization of the ada gene cluster. (B) Biosynthesis of TAN-1612 (3) in A. niger ATCC 1015.130 |
Azaphilones represent a class of metabolites that consist of a highly-oxygenated bicyclic core and a chiral quaternary center. This family is structurally diverse and includes numerous bioactive compounds.170 Studies have scanned the sequenced genome of A. niger ATCC 1015 for dual PKS gene clusters and located a single gene cluster containing a HR-PKS azaA and a NR-PKS azaB. These appeared to be homologs of afoE and afoG in the A. nidulans afo cluster.131,168,171 Analysis of the nearby genes in the dual PKS cluster revealed the presence of a Zn(II)2Cys6 zinc finger transcription factor azaR. Further investigation revealed that several tailoring enzymes were likely involved in modification of a final polyketide product. Interestingly, most of these genes have a corresponding homolog in the asperfuranone gene cluster.168 The extensive overlap between the aza and afo BGCs clearly indicated that the aza cluster encoded for the production of an asperfuranone-like compound (Fig. 5).131
Fig. 5 The A. niger aza cluster and the proposed biosynthetic pathway of azanigerones from the aza cluster.131 |
RT-PCR transcriptional analysis indicated there were 12 genes in the aza cluster were under control of the pathway-specific regulator azaR.131 After the overexpression of azaR, azanigerone A (1) was obtained as the most abundant product. Several other related compounds were observed at earlier time points, including FK17-P2a (25), azanigerones B (26) and C (27). In addition, azanigerone D (28) was observed to replace azanigerone A as the dominant product at later time points.131
In order to investigate the collaborative mechanism of NR-PKS AzaA and HR-PKS AzaB in azanigerone biosynthesis, an azaB-disruption strain was generated, which upon culturing led to the accumulation of two new compounds, identified as azanigerones E (29) and F (the O4-acetylated version of 29). This suggests a convergence synthesis model, in which AzaA and AzaB biosynthesize two discrete polyketide products that are combined at a later step in the pathway. Importantly, this is the first report of convergent collaboration between a NR-PKS and a HR-PKS in secondary metabolite biosynthesis.131 Furthermore, through enzyme expression and a series of in vitro assays, the critical role of AzaH was confirmed in morphing the benzaldehyde intermediate (FP17-2Pa) into the bicyclic, pyran-containing azaphilone core structure. Altogether, these findings enabled a pathway for biosynthesis of the azanigerones in A. niger to be proposed.131
The biosynthesis of azanigerone A (1) begins with the polyketide assembly by AzaA, forming the hexaketide precursor (31) from successive condensations of five malonyl-CoA units. Then the reduction of the terminal ketone is catalyzed by the ketoreductase AzaE resulting in the formation of 25. Subsequently, the monooxygenase AzaH was demonstrated to be the only enzyme to hydroxylate 25, which facilitates formation of the pyran-ring and generation of 29 after dehydration. In parallel, the 2,4-dimethylhexanoyl chain (31) is synthesized by the HR-PKS AzaB and is combined with 29 to generate 26 in a reaction facilitated by the acyltransferase AzaD. Next, 26 is hydroxylated by the FAD-dependent monooxygenase AzaG or AzaL to yield 27, followed by C–C oxidative cleavage by cytochrome P450 AzaI to obtain 32, and oxidation of the aldehyde to a carboxylic acid by AzaJ, yielding 1. 1 can then undergo addition of NH3 to generate 28, which was observed to replace 1 production at later time points (Fig. 5).
Since the structures of 5, 33 and 34 all contain one phenylacetate and two acetate units, this suggests that they may originate from a polyketide biosynthetic pathway. The PhlB gene (phenylacetyl-CoA ligase) of Penicillium chrysogenum is responsible for the incorporation of phenylacetate,173 a homologous BLAST search was thus performed resulting in the identification of ORF An09g01820 (epaB) in A. niger. Analysis of nearby genes in the epa cluster revealed the presence of a core hybrid NR-PKS enzyme EpaA, an acyl-CoA transferase EpaC, and three genes orf 1–3 encoding an oxidoreductase, a 3-hydroxybenzoate 4-hydro-xylase, and a salicylate hydroxylase, respectively.172
The NR-PKS responsible for the production of 5, 33 and 34 was confirmed to be epaA through gene deletion by homologous recombination, which led to the abolishment of the production of all three metabolites. The complementary strain (epaA-com) was able to produce them all, which further confirmed that the identified BGC is responsible for the synthesis of 5, 33 and 34. In addition, feeding purified 5 and 33 into cell-free extracts (including all biosynthetic enzymes and cofactors for the biosynthesis of 34) of the ΔgcnE mutant did not produce 34, suggesting that 5 and 33 were by-products in the biosynthesis of 34.172
Based on these results, the biosynthetic pathway of 5, 33 and 34 was proposed (Fig. 6).172 The pathway starts with the biosynthesis of phenylacetyl triketide (35) by the NR-PKS EpaA. This prompted the researchers to hypothesize the presence of a phenylacetyl-CoA initiator unit, followed by releasing of a polyketide chain as an aldehyde via the R-domain. Next, the cyclization and dehydration by nonenzymatic rearrangement likely generates 5.172 For the biosynthesis of 33 and 34, an additional methyl group is added through the C-methyltransferase domain before the final chain release, yielding 36. Subsequent, reduction to generate an aldehyde product, oxidation and transamination to produce an intermediate amide (37), followed by enolization, cyclization and water loss to generate 33.172 The researchers proposed that EpaC (An09g01800), an acyl-CoA transferase can additionally catalyze the transfer of 2-methylsuccinyl-CoA, a common intermediate in the ethylmalonyl-CoA pathway, to generate an alternative product 38 and then form 34 by cyclization and water loss.172
Fig. 6 Putative gene cluster and the proposed biosynthetic pathway of nigerpyrone, carbonarone and pestalamide in A. niger.172 Adapted from Wang et al. |
Heterologous expression of dtbA in A. nidulans demonstrated that DtbA catalyzes production of two polyketides, 2,4-dihydroxy-3,5,6- trimethylbenzaldehyde (4) and 6-ethyl-2,4-dihydroxy-3,5-dimethylbenzaldehyde (39) (Fig. 7). In order to further expand the structural diversity of products produced by fungal NR-PKSs, this BGC was engineered with a TE release domain instead of the R-domain. By replacing the DtbA R domain with the AusA (austinol biosynthesis) or ANID_06448 TE domain, DtbAΔR+TE chimeric PKSs were formed allowing the development of two metabolites with carboxylic acids replacing the corresponding aldehydes (40, 41). The TE domains were selected based on phylogenetic relatedness to DtbA, as well as NR-PKS with a similar CMeT domain location. This indicates that the release domains are selective catalysts that co-evolved with their cognate CMeT domains, rather than being simple non-selective and slow hydrolases.133
Fig. 7 (A) Domain architecture, putative starter units, unreleased intermediates, and released polyketides for DtbA in A. niger ATCC 1015. (B) Domain architecture, putative starter units, unreleased intermediates, and released polyketides for hybrid DtbAΔR+TE(AusA or ANID_06448) in A. nidulans.133 |
Fig. 8 (A) The gene cluster phn and (B) the proposed biosynthetic pathway of phenalenone in P. herquei.134 |
The previous study found that the core of phenalenone polyketides is synthesized by the combined functions of a NR-PKS (PhnA) and a flavin-dependent monooxygenase (FMO) (PhnB) that catalyze an unusual aromatic hydroxylation of the NR-PKS product prephenalenone (43). PhnA synthesizes a novel angular naphtha-γ-pyrone product that is oxidatively morphed into the phenalenone structure by a FMO PhnB. Importantly, the combined activities of PhnA and PhnB represent a new strategy of polyketide cyclization in nature. In addition, the neighboring genes encode FMO (PhnB and PhnG), O-MT (PhnC), and prenyltransferase (PrT, PhnF) that are consistent with the structural features of herqueinone.134
To confirm the involvement of phnA in the biosynthesis of herqueinone, the KS region of phnA was deleted using homologous recombination and hygromycin as selection marker, leading to the complete loss of the production of herqueinone. To identify the polyketide product of PhnA and the functions of its cyclization domains, phnA was heterologous expressed in Saccharomyces cerevisiae strain BJ5464-NpgA, yielding the naphtha-γ-pyrone 43. Feeding experiments with purified 43 to the ΔphnA-blocked mutant restored the production of herqueinone, which further confirmed 43 as a key intermediate in the biosynthesis of herqueinone.134
To investigate the role of the PT domain, as well as the TE/CLC domain, the researchers constructed an inactivated TE domain of PhnA either by point mutation or domain deletion (PhnA S2009A and ΔTE), which led to the production of two macrolactones (44, 45) instead of 43. The formation of 44 and 45 in this process requires ketone reduction of the C11 and C13 carbonyl groups, respectively, and subsequent attack by the hydroxyl groups on the C1 thioester. It should be noted that the ketone reduction steps are catalyzed by endogenous yeast ketoreductase and not the domains inside PhnA. Therefore, although the TE/CLC mutation abolished the second aromatic ring derived from the C1–C10 Claisen condensation in 43, whereas first aromatic ring fixed from the C4–C9 aldol condensation remained intact. This confirmed that the PT domain in PhnA is only responsible for directing C4–C9 cyclization while the TE/CLC domain catalyzes C1–C10 cyclization.134
To identify the enzymes responsible for the pyrone opening and C8–C13 cyclization for the formation of 42, targeted gene deletion of selected phn genes was carried out. The inactivation of phnB led to the abolishment of production of herqueinone and the accumulation of 43, which indicated its important role in transforming 43 into a phenalenone structure. Further investigation showed that PhnB was responsible for converting the γ-pyrone 43 into the tricyclic phenalenone as well as C2 hydroxylation.
Thus, following release of the polyketide product from PhnA and rapid formation of the hemiketal 43, the active site of PhnB catalyzed the C2-hydroxylation through deprotonation of the hemiketal C13-OH to yield a ketone intermediate (46). Subsequently, tautomerization of the C3 enol group leads to rapid cyclization to yield the cyclized intermediate, which then aromatizes through dehydration and keto–enol tautomerization to yield the phenalenone product. The remaining enzymes encoded in the gene clusters, including PhnC (O-MT), PhnF (PrT) and PhnG (FMO), have been confirmed to transform prephenalenone to herqueinone.134 Funalenone is the C2-methoxy derivative of phenalenone, and the structural features are reflected in the A. niger gene cluster, which only contains homologues to PhnA, PhnB, PhnC (O-MT), and PhnD (Table 5).
Proposed function of encoded enzyme | A. niger | |
---|---|---|
a In the clusters listed here, each contains a PhnB homolog with >60% sequence identity. b CAK38307 contains two conserved domains, corresponded to PhnB and PhnD, respectively. | ||
PhnA | Nonreducing polyketide synthase (SAT-KS-AT-PT-ACP-ACP-TE/CLC) | CAK38306 |
PhnB | FAD-dependent hydroxylase | CAK38307 |
PhnC | O-Methyltransferase | CAK38305 |
PhnD | Unknown function (Dabb superfamily protein) | CAK38307b |
PhnE | Transcriptional factor | No homolog |
PhnF | Aromatic prenyltransferase | No homolog |
PhnG | FAD dependent pyridine nucleotide-disulfide oxidoreductase | No homolog |
Next, the acu gene cluster, including the acuA (51 synthase encoding gene) and eight surrounding genes acuB–I, was defined (Fig. 9), leading to the proposed biosynthetic pathway for the formation of 48 and 52 from 51 (Fig. 10). The pathway initiates the biosynthesis of 51 by the PKS AcuA, followed by decarboxylation by AcuB to produce 57, similar to what is seen in 49 biosynthesis, followed by oxidation catalyzed by AcuC to produce 58. Next, AcuD and AcuE catalyzes two oxidations of the aromatic ring in the gentisyl alcohol (59) and 60. Inspired by how 51 is used as the starting point in biosynthesis of patulin, the researchers proposed that the following aromatic ring opening is catalyzed through a Baeyer–Villiger type of oxidation, catalyzed by AcuF, with the carboxylic acid at C1 subsequently being reduced to an aldehyde by AcuG, yielding 61. Next, as also proposed for the patulin biosynthetic pathway, a hemiacetal is formed,178 followed by the reduction of the double bond in the side chain by a putative dehydrogenase AcuH. Finally, by non-enzymatic hemiacetal formation, keto–enol tautomerism results in the production of 48, which exists as two diastereomers (both R/S configurations at C1).
Fig. 10 Proposed biosynthetic pathway of aculinic acid and aculin A in A. aculeatus.52 |
In addition, the 55 moiety was proposed from the terpene pathway and recently confirmed it is derived from a 62 precursor through demethylation and will be further described in this review in Section 6.2.153 However, the linking of 48 and 55 is still unclear but it was hypothesized that AcuI is involved in the attachment of proline to 48 to form 52-54. Therefore, 6-MSA is incorporated into 48, which then incorporates into 52–54. Ultimately, further exploration is needed to elucidate the intermediate biosynthetic pathway steps and confirm the function of enzymes involved in acu gene cluster.
Fumonisins, along with host-specific toxins (AAL), are known as sphinganine-analog mycotoxins (SAMT) due to their structural similarity to the early precursors of sphingolipid pathways. Fumonisins consist of a linear polyketide-derived backbone with two methyl, one amine, one to four hydroxyls, and two tricarboxylic ester groups located at different positions along the backbone (Fig. 11). Because this cluster was first discovered and identified in Fusarium species, we first describe the FUM cluster in Fusarium and then the homologous gene cluster fum in Aspergillus species.
A cluster of 17 genes required for fumonisin biosynthesis and regulation was identified in Fusarium using heterologous expression which include: FUM1, FUM6–8, FUM3 (FUM9), FUM10–11, FUM2 (FUM12), FUM13–19, FUM 20–21.138 These genes can be categorized into four different groups: (a) backbone gene encoding the dimethylated polyketide that undergoes condensation with alanine (PKS FUM1 and aminotransferase FUM8), (b) decorating post-PKS steps of biosynthesis (FUM3 (FUM9), FUM2 (FUM12), FUM6, FUM13, FUM7, FUM10, FUM14), (c) regulatory genes or transcription factors (FUM21) and (d) unassigned or non-essential genes (FUM11, FUM19, FUM15-18 and FUM20). In the B series fumonisins that consist of a linear 20-carbon backbone, the biosynthetic pathway starts with linear dimethylated polyketide formation (C3–C20), followed by the polyketide condensation with alanine (C1–C2), and finally by carbonyl reduction, oxygenation and esterification with two propane-1,2,3-tricarboxylic acids.185
The fum cluster in A. niger includes 11 homologue genes to Fusarium spp.: the polyketide synthase (fum1), three hydroxylases (fum3, fum6, fum15), a dehydrogenase (fum7), an aminotransferase (fum8), an acyl-CoA synthase (fum10), a carbonyl reductase (fum13), a condensation-domain protein (fum14), an ABC transporter (fum19), and transcription factor (fum21) (Fig. 12). However, the Aspergillus cluster does not contain FUM11, 2, 16, 17, and 18. The absence of FUM2, which catalyzes C-10 hydroxylation, is consistent with the fact that A. niger produces 64, 65 and 66, but not 63 and FB3 (67) (Fig. 13).137 In addition, the absence of FUM11, 16, 17, and 18 in A. niger also supports the hypothesis that these genes are not essential for fumonisin production.138 In contrast, the BGC in A. niger includes a short-chain of dehydrogenase gene (designated here as sdr1) that is not present in the Fusarium cluster for which the function has yet to be characterized.
Fig. 12 Fumonisin gene cluster in F. verticillioides and A. niger, abbreviated as FUM in Fusarium and fum in Aspergillus. Genes a–d (dark blue arrows) are predicted to encode: a – polyketide synthase; b – major facilitator superfamily transport protein; c – 5-methyltetrahydropteroyltriglutamate-homocysteine methyltransferase/methionine synthase; d – alkaline phytoceramidase.136 |
Fig. 13 The biosynthetic pathway of FB1, B2, B3 and B4. The structures in dotted frame cannot be produced by Aspergillus species because of the absence of FUM12.138 |
In Aspergillus spp., the B series fumonisin biosynthetic pathway starts with formation of a linear dimethylated polyketide by fum1 and polyketide condensation with alanine by fum8, followed by carbonyl reduction by fum13. Two hydroxyls are introduced at C-14 and C-15 by fum 6, followed by esterification with two propane-1,2,3-tricarboxylic acids by fum 7, 10 and 14 to form FB4. Next, the C-5 hydroxylation is catalyzed by fum3 to produce 64.
Additionally, one study investigated the involvement of the fum8 in 64 (biosynthesis in A. niger). The disruption of fum8, coding for an α-oxoamine synthase, resulted in the loss of 64 biosynthesis in A. niger indicating that this gene is involved in 64 biosynthesis in A. niger. The study also showed that 64 is not necessary for vegetative growth and UV tolerance in this fungus.186 Moreover, a study of the black aspergilli from California grapes showed that fumonisin-nonproduction was associated with the absence of six fum genes in A. welwitschiae but not in A. niger.100,136
Recently, two classes of non-aminated fumonisins named FPy and Fla (68 and 69), were isolated from A. welwitschiae, making them the first reported naturally occurring fumonisins lacking an amine.112 These non-aminated metabolites showed significantly less toxicity in the Lemna minor (duckweed) bioassay than the fumonisin B type, providing new insights into the toxicity mechanism of fumonisins. Further exploration with 13C and 15N stable isotope incorporation suggested a post-biosynthetic oxidative deamination of 65. This work shows that certain Aspergillus species can produce an enzyme that naturally detoxifies fumonisins. These enzymes may have many practical applications in food production systems.112
The proposed biosynthetic pathway of 70 and 71 start with the generation of cinnamic acid (Fig. 14), which was catalyzed from phenylalanine by the phenylalanine ammonia-lyase, AhpB. The putative ligase, AhpE, might ligate the cinnamic acid to produce cinnamoyl-CoA, which is the activated starter unit in the biosynthesis. Next, AhpA catalyzed the biosynthesis of the backbone of the homopyrones by the addition of three malonyl-CoA extender units to offer a 15-carbon chain.125 Additionally, the KR and DH in PKS AhpA were predicted to catalyze the reduction at C-9 and the subsequent water loss, respectively. A freestanding thioesterase (TE) named AhpC was proposed to be responsible for the hydrolytic release and the following spontaneous lactonization. After the pyrone formation, the products are proposed to be methylated on C-13 by the O-methyltransferase ahpD.125
Fig. 14 The gene cluster ahp and the proposed biosynthetic pathway of homopyrone A and B in A. homomorphus.125 |
Although the canonical function of C domains is to catalyze the formation of a peptide bond between the growing peptide, the CT domain (a subset of C domains) can be recruited by fungal NRPS to perform the equivalent reaction of a thioesterase domain. To date, all fungal NRPSs that assemble macrocyclic peptides terminate with a condensation-like CT domain that catalyzes the cyclization reaction, and thus, is a near universal macro-cyclization strategy used by fungal NRPSs.188,190 Optional domains include epimerization (E), formylation (F), methylation (M), heterocyclization (Cy), reduction (R), and oxidation (Ox) domains for additional modifications. In this section, aculeacin A (72)/echinocandin B (73) (ecd and hty) and malformins (mlfA) are secondary metabolites produced from section Nigri derived from NRPSs. Additionally, a putative BGC for asperazine biosynthesis was proposed in the A. tubingensis G131 and this cluster shows homology with the fumitremorgin biosynthesis cluster identified in A. fumigatus Af293 and the ditryptophenaline production in A. flavus.191,192 However, genetic experiments should be performed in the future to confirm the involvement of this gene cluster in biosynthesis of asperazine.
Fig. 15 Gene clusters of echinocandin B (73) in A. pachycristatus and aculeacin A (72) in A. aculeatus.139 Genes in dark blue are of unknown function or believed not to participate in the biosynthesis. |
The researchers proposed one possible pathway for the biosynthesis of 73 (Fig. 16): linoleoyl-AMP, produced by EcdI from linoleate, is transferred to the initiation T0 (first thiolation domain) of EcdA. The linoleoyl-S-phosphopantetheinyl-T0 is sequentially extended with L-ornithine, L-threonine, L-proline, L-homotyrosine, L-threonine, and 4R-methyl-L-proline to form the linear hexapeptide. Thereafter, the terminal condensation (CT) performs macrocyclization of the NRPS product and the cyclic scaffold is released from EcdA. In addition, an intriguing structural feature of echinocandins is the incorporation of hydroxylated amino acids. Later research identified five oxygenation events by these three iron dependent enzymes, EcdG, EcdH and EcdK.195 EcdH, a predicted P450 type hemeprotein monoxygenase, has been shown to iteratively hydroxylate C5 and C4 of L-ornithine after the formation of the cyclic peptide, while the ecdG gene encodes α-ketoglutarate (α-KG)-dependent mononuclear, nonheme iron oxygenase that has been shown to hydroxylate free L-homoTyr to generate 3-hydroxyl-L-homoTyr. Furthermore, another α-KG-dependent nonheme iron oxygenase EcdK is involved in hydroxylating C5 of methyl-proline, to produce 4R-methyl-proline (4-Me-Pro).195 However, the last echinocandin oxygenation event is still to be characterized.
Fig. 16 Proposed biosynthetic pathway of echinocandin B (73). Confirmed hydroxylation steps are highlighted in red, while proposed one is in blue.195 |
Moreover, to investigate the evolution of chemical diversity in echinocandin lipopeptide metabolites, the gene cluster responsible for 72 synthesis was characterized via DNA sequences for additional echinocandin-type pathway gene clusters.139 In the gene cluster of 72, the hty genes are located immediately upstream of and are separated by, two unknown genes and a cytochrome P450 from the core NRPS gene. Recently, Zhang et al. characterized the α-ketoglutarate/Fe2+-dependent proline hydroxylase (HtyE) from A. pachycristatus and A. aculeatus in the respective 73 and 72 BGCs through heterologous expression in Escherichia coli.196,197 This study showed that both Ap- and Aa-HtyE converted L-proline to trans-4- and trans-3-hydroxyproline but at different ratios.
Fig. 17 The structures of malformin A2 (74) and malformin C (75), and predicted gene cluster for malformin production in A. brasiliensis.107 |
Since malformins are pentapeptides that contain a disulphide bond, only one gene cluster (mlf) containing NRPS genes and disulphide bond-associated genes were included as search criteria. This search led to the identification of the likely related gene cluster responsible for biosynthesis of malformin in 18 fungal strains, importantly, including previous reported malformin producing species.107,109 The NRPS responsible for malformin biosynthesis was confirmed to be mlfA through targeted gene deletion in A. brasiliensis, which led to a total abolishment of 74 and 75 production. Furthermore, a genetic complementation by a constitutively expressed mlfA in the mlfAΔ strain restored 74 and 75 production, to further confirm the function of mlfA. Based on these results, it was hypothesized that some of the individual NRPS adenylation domains are less selective since multiple malformins disappear after deletion of the NRPS.
When the researchers heterologously expressed atromentin synthetases in A. niger in order to explore the reaction mechanism of the thioesterase domain, they discovered the metabolite atrofuranic acid (76), as well as metabolite X (77), gyroporin (78) and p-hydroxybenzoic acid (decomposed product of atrofuranic acid). However, atromentin synthetases produce atromentin in A. oryzae, which implied that the physiology of A. niger initiates a cross-specificity and cross-chemistry event on the NRPS-like enzyme. When feeding with 2-13C-labeled L-tyrosine, which was incorporated into 76 as a precursor, the biosynthesis pathway for metabolites produced in A. niger could be proposed (Fig. 18). Biosynthesis of 78 was suggested to be proceed via an uncharacterized metabolite, 77, proposed to be derived by a Piancatelli-type rearrangement from 76,202,203 which subsequently by decarboxylation and oxidation leads to formation of 78.
Fig. 18 Proposed biosynthetic pathway of atrofuranic acid (76) in Aspergillus species from section Nigri. A subsequent Piancatelli-type rearrangement and decarboxylation results in the formation of gyroporin (78).141 |
Further exploration showed that this cross-chemistry was not limited to A. niger but also occurred in A. brasiliensis, which contains an intrinsic atromentin synthetase, indicating that modifications of benzoquinone core structures are common to Aspergillus species of the section Nigri. Therefore, characterization of products from NRPS-like enzymes may be hampered by cross-chemistry occurring on metabolite core structures during heterologous expression.141 It appears that using different heterologous expression platforms in parallel to elucidate the true nature of products generated by uncharacterized NRPS-like enzymes is critical.
The BGC which is related to production of 79 (not to be confused with pyranonigrin E reported in another study210) was induced and the metabolites was isolated in A. niger ATCC 1015 as a new PKS–NRPS hybrid secondary metabolite after the pyn gene cluster was activated by the expression of the pathway-specific Zn2Cys6 transcriptional regulator pynR under control of the arginase (aga) promoter.144 Furthermore, cluster activation resulted in induced transcriptions of the pyranonigrin encoding genes pynA, pynB, pynC, pynD, and pynE. PynA encodes a PKS–NRPS hybrid enzyme, PynB encodes a FAD/FMN-containing dehydrogenase, PynC encodes a methyltransferase, PynD encodes a cytochrome P450 oxidase, and PynE encodes an NAD-binding protein.144
In subsequent research, Yamamoto et al. found more pyranonigrins through activating the pyn cluster in high yields by replacing of the pynR promoter with the robust glaA promoter. Subsequently, the complete pathway of pyranonigrins was identified by targeted gene deletion, heterologous in vivo bioconversion and in vitro assays. Furthermore, this study identified three extra genes, including the protease-like protein PynH that performs an exo-methylene formation, the stand-alone thioesterase PynI that is capable of concomitant tetramic acid formation and chain release from the PKS–NRPS hybrid enzyme PynA, and FAD-containing monooxygenase (FMO) PynG that catalyzes an epoxidation-mediated cyclization to form the dihydro-γ-pyrone moiety.145
Altogether, the proposed biosynthetic pathway initiates with biosynthesis of the PKS–NRPS product by PynA, followed by release of the intermediate from PynA by PynI, generating 84. Next, PynC catalyzes N-methylation of 84, resulting in 83. Then PynG catalyzes an epoxidation-mediated cyclization to form the dihydro-γ-pyrone moiety, which is followed by PynD-catalyzed alcohol oxidation and enolization to form 82. PynH then catalyzes the formation of 79, which is dimerized to form 80. Alternatively, PynE can catalyze the conversion of 79 to 81, which can be reversed by PynB.144,171
A very notable discovery in this pathway was PynH, a predicted protease that performs dehydration of a serine side chain to form an exo-methylene in pyranonigrin E. It is highly unusual to find an exo-methylene-forming enzyme that is involved in the modification of a fungal polyketide-nonribosomal peptide hybrid compound. Another important finding is that the stand-alone thioesterase PynI is capable of tetramic acid formation and chain release from the PKS–NRPS hybrid enzyme PynA (Fig. 20).
Fig. 20 Biosynthesis pathway of pyranonigrins E–J (79–84) in A. niger.145 |
In addition, the potent antioxidant pyranonigrin A (85), which was first reported in A. niger,211,212 was found to be biosynthesized by a PKS–NRPS gene cluster pyr in P. thymicola, which is different from the pyn cluster for 79. Its biosynthetic pathway is elucidated in P. thymicola by genome mining and heterologous expression.146
Basic Local Alignment Search Tool (BLAST) search results revealed a PKS–NRPS (PyrA) from P. thymicola IBT 5891 that exhibits moderate sequence similarity to PynA (50% identity/67% similarity) and An18g00520 (61% identity/75% similarity), both of which are from A. niger A1179.213 In addition, An18g00520 was shown to be responsible for the production of tetramic acid type compound, proposed to be the precursor to 86.214 Based on these information, the heterologous expression of pyrA in A. nidulans was carried out, resulting in the production of the product of An18g00520. This indicated the product encoded by pyr cluster shares the same biosynthetic precursor with 86.146 Subsequent heterologous expression of all the genes involved in the entire pyr gene cluster in A. nidulans led to the production of 85. Based on the findings from both in vivo and in vitro studies, only four genes, pyrA–D, are required for the biosynthesis of 85.146
Therefore, the biosynthetic pathway of 85 in P. thymicola initiates the biosynthesis of the PKS–NRPS product by the PyrA, followed by product release catalyzed by either a Dieckmann cyclase (DKC) and/or the hydrolase PyrD. Then, FAD-binding monooxygenase PyrC may catalyze epoxidation followed by subsequent ring closure to form the pyrano[2,3-c]pyrrole core, followed by conversion to 86 and finally the α-carbon of the glycine part of 85 is oxidized by cytochrome P450 PyrB (Fig. 21).
Fig. 21 The proposed biosynthesis of pyranonigrin A (85) in P. thymicola IBT 5891.146 |
The BGC pyv of 87 was identified via bioinformatics analysis and comparison to the homologous genes of pyranonigrins (Fig. 22).144,146,214 This revealed the pyv cluster encodes a PKS–NRPS hybrid PyvA, a cytochrome P450 monooxygenase PyvB, two flavin-dependent monooxygenases (FMOs) PyvC and PyvE, an α/β hydrolase PyvD, a glucose-methanol-choline (GMC) family of oxidoreductase PyvF, and O-methyltransferase PyvH. The PKS–NRPS responsible for 87 was confirmed to be PyvA through targeted gene deletion, which led to complete abolishment of 87 production. Based on the bioinformatic analysis and previous studies on pyranonigrins, the biosynthetic pathway leading to 87 was proposed (Fig. 23).
Fig. 23 Proposed biosynthetic pathway of pyranoviolin A (87) in A. violaceofuscus.147 |
The pathway initiates by the PKS PyvA to synthesize C10 backbone from one acetyl-Co unit and four malonyl-CoA units, followed by the condensation with the thiolation (T) domain-bound glycine activated by the adenylation (A) domain. Then the product release is catalyzed by Dieckmann condensation (DKC) and/or the α/β hydrolase PyvD, yielding the tetramic acid moiety. Next, the FMO PyvC may catalyze the epoxidation of one of the olefins on the polyketide chain, followed by the epoxide ring-opening and the dihydro-γ-pyrone ring formation. Next, the P450 PyvB maybe catalyze two consecutive reactions, in which the dihydro-γ-pyrone is oxidized to γ-pyrone and C-7 is hydroxylated to yield pyranonigrin F (88) [note that two distinct natural products are individually named “pyranonigrin F”].145,215 Finally, the O-methyltransferase PyvH, which is specifically found in the 87 pathway, catalyze the methylation for the C-3 hydroxy group to produce 87.
In Aspergillus species, a PKS required for 50 biosynthesis was identified in A. ochraceus and partially characterized in Aspergillus westerdijkiae.220,221 In addition, in A. carbonarius, the AcOTAnrps gene encoding the NRPS responsible for the ligation of the phenylalanine to the polyketide group, was identified by gene disruption.148 In the same study, the order of reactions in the biosynthetic pathway of 50 in A. carbonarius was clarified demonstrating that the amide bond between the phenylalanine and the polyketide dihydroisocoumarin is catalyzed by the synthetase and that it precedes the chlorination step.148 Furthermore, a flavin-halogenase in the biosynthesis of 50 in A. carbonarius was identified by gene deletion and confirmed that the chlorination step which converts ochratoxin B (OTB, 89) to 50 represents the final stage of the biosynthetic pathway.
Recently, in order to study the regulatory mechanism of 50 biosynthetic pathway, Wang et al. sequenced and assembled the complete genome of A. ochraceus fc-1 and used that for comparative genomic analyses with five other sequenced ochratoxin A-producing fungi (A. carbonarius, A. niger, A. steynii, A. westerdijkiae, and P. nordicum).149 By comparing the genomic analysis of these six 50 producers, a conserved ochratoxin A (50) BGC was identified, including a PKS (otaA), an NRPS (otaB), a cytochrome P450 monooxygenase (otaC), a halogenase (otaD), and a basic leucine zipper (bZIP) transcription factor (otaR1), which are all required for 50 biosynthesis. A GAL4-like Zn2Cys6 binuclear cluster DNA-binding protein (otaR2) and a FAD-dependent oxidoreductase (otaE) were also identified but they are not essential to the 50 production (Fig. 24). Ochratoxin A (50) biosynthesis begins with the formation of 7-methylmellein catalyzed by OtaA utilizing acetyl-CoA and four malonyl-CoA units, followed by the oxidation to OTβ by OtaC and OtaE. Next, OTβ and L-phenylalanine are combined by the NRPS OtaB to form an amide bond to synthesize 89. Finally, 89 is chlorinated by a halogenase OtaD to give 50 (Fig. 25).
Fig. 24 Proposed ochratoxin A (50) BGC family in A. ochraceus, A. niger and A. carbonarius. Genes shown in the same color displayed high similarity with each other. |
Fig. 25 Proposed ochratoxin A (50) biosynthetic pathway, including intermediate metabolites and catalytic enzymes.148 |
Since the structure of 90 is very similar to that of fusarin C, for which the biosynthesis is partly characterized, it was hypothesized that 90 originated from a polyketide and an amino acid fusion. The candidate clusters were further narrowed down to identify the putative PKS and NRPS genes for the production of 90. Intriguingly, bioinformatics analysis indicated that the scaffold for 90 formation potentially could be synthesized by unlinked PKS and NRPS enzymes. In addition, the researchers also found that uncoupled PKS and NRPS encoding genes are conserved in a fungal clade in section Nigri. Moreover, the PKS didn't appear to have a release domain but instead the predicted NRPS contained a reductase domain for the potential release of 90 backbone. Examination of genes encoding the PKS and NRPS responsible for 90 biosynthesis were identified to be acrA and acrB through gene deletion experiment in A. aculeatus, which led to the abolishment of 90 production. Further bioinformatics analyses showed a lack of an epimerase domain is present in the NRPS indicating L-serine is incorporated into 90. In addition to the PKS and NRPS encoding genes of 90, there are six other surrounding genes (acrC–G) including a regulatory transcription factor (acrR). Next, a functional analysis of the genes involved in the formation of 90 was also performed by targeted gene deletions. This confirmed the involvement of the transcription factor acrR and O-methyltransferase (acrG). However, no accumulation of potential intermediates were found in acrC (hydrolase), acrD (oxidoreductase), acrE (CYP450) and acrF (CYP450) deletion strains, which might be because of the unstable intermediate products or a scenario where more than two enzymes interact with each other during the reactions. Altogether, these results enabled Wolff et al. to propose a pathway for biosynthesis of 90.150
The pathway starts with the polyketide assembly by the PKS AcrA forming a heptaketide backbone from successive condensations of six malonyl-CoA units (Fig. 26). Then the fusion of an amino-acid moiety is catalyzed by the NRPS AcrB. The functions of AcrC–F are still unclear, however the researchers hypothesized that the PK-NRP hybrid product release from the enzyme by reductive release to set up the formation of the lactam ring by Aldol condensation, followed by a water loss to generate a double bond followed by reduction. Finally, an O-methyltransferase AcrG is needed for the final step in the biosynthesis of 90.
Fig. 26 The gene cluster in A. aculeatus and the proposed biosynthetic pathway of acurin A (90).150 |
The yan BGC for 49 contains ten genes in the genome of A. niger. The PKS-encoding YanA was heterologously expressed in A. nidulans to confirm its involvement in the biosynthesis of the initial polyketide product 6-MSA (51). The biosynthetic pathway of 49 was further investigated by individually deleting all genes within the yan cluster (Fig. 27) which resulted in the accumulation of numerous yanuthone intermediates and enabled a biosynthetic pathway to be proposed (Fig. 28).
Fig. 27 The yan cluster. The yanA 6-MSA synthase-encoding gene is flanked by nine cluster genes (yanB, yanC yanD, yanE, yanF, yanG, yanH, yanI, and yanR) whose products contain all necessary activities for conversion of 6-MSA into yanuthone D (49).151 |
Fig. 28 Biosynthesis of yanuthone D (49) in A. niger.151,223 |
The pathway starts with the biosynthesis of 6-MSA (51) catalyzed by the PKS YanA, followed by decarboxylation by YanB producing m-cresol (57), and hydroxylation by the cytochrome P450 YanC, generating toluquinol (91). The latter two steps for epoxidation and prenylation are still not very clear since no biosynthetic intermediate has been obtained from the targeted gene deletions. It has been proposed that the farnesylation likely precede epoxidation since the C-prenylation reaction generally requires the π-electron of a prenyl acceptor to attack the allylic cation generated on the prenyl donor. Therefore, 91 may also be accepted by YanG, followed by epoxidation by YanD/YanE generating 7-deacetoxyyanuthone (92), as epoxidized 91 loses a π-electron at the prenylation site due to the epoxide formation.151,223 Next, cytochrome P450 YanH catalyzes the conversion to 22-deacetylyanuthone A (93), which is then converted to yanuthone E (94) by the O-mevalon transferase YanI. Importantly, this is the first time that O-mevalon transferase activity has been molecularly characterized. Finally, the oxidase YanF catalyzes the formation of 49 from 94.171 Deletion of the individual genes in the BGC also led to formation of not only hydroxylated metabolites but also a glycosylated yanuthone.50
In addition to 49 and its analogues with a C7 scaffold defined as class I yanuthones, yanuthone X1 and X2 were also discovered, which instead contain a C6-core scaffold oxygenated at C-16 defined as class II yanuthones and derived from an unknown precursor.108,222
Because the L-proline moiety in 55 was likely introduced by a nonribosomal peptide synthase (NRPS), the researchers accessed the genome sequence from JGI of A. aculeatus ATCC16872 (=BCRC32190) and scanned for all adenylation domain encoding genes and their neighboring genes.224 A candidate BGC, ane, was found containing one single-module NRPS (AneB) with adenylation and thiolation domains as well as a Class I terpene cyclase (AneC). Further investigation revealed several genes nearby aneB and aneC with putative activities consistent with the structure of 55, including a phytanoyl-CoA dioxygenase (aneA), three cytochrome P450 monooxygenases (aneD, aneF and aneG) that could perform oxidative modifications on the terpene scaffold, a α,β-hydrolase (aneE) that may transfer or release L-proline from AneB to the sesquiterpene moiety, and a proline iminopeptidase (aneH) (Fig. 29).
Fig. 29 The gene cluster ane in A. aculeatus and the biosynthetic pathway of aculene A.153 |
To further explore the functions of those enzymes, genetic inactivation on aneC was performed which produced no aculene products and showed no accumulation of any intermediates. Next, intron-free aneC from A. aculeatus was cloned and transformed into Saccharomyces cerevisiae strain BJ5464-NpgA, which showed the production of dauca-4, 7-diene (95) as a sole product. These results confirmed that AneC is a terpene cyclase that produce 95 with an ent-daucane skeleton for the biosynthesis of 55.
Each of the three P450 monooxygenase genes (aneG, aneF and aneD) were then inactivated. The results indicated the involvement of oxidation by AneG, AneD and AneF at C-2, C-10 and C12 respectively. In addition, the researchers hypothesized that the P450s could also play a role in the demethylation process. To further verify the function of the three P450 monooxygenases (ane G, D, F), heterologous co-expression of these enzymes with AneC using a S. cerevisiae strain harboring a single copy of the cytochrome P450 reductase gene from A. terreus (AtCPR) were performed first.224 They confirmed the role of aneF, aneD, aneG when they coexpressed aneF with aneC, aneC/F/D and aneC/G/F. However, the enzyme involved in the second oxidation at C-2 to generate a ketone functionality from a C-2 hydroxy group and decarboxylation at C-12 was still unknown because of the decreasing amounts of downstream metabolites when coexpressing aneC/F/D/G. This might be because of degradation by endogenous enzymes in S. cerevisiae of the aculenes homologues. Therefore, the coexpression of aneC/F/D/G in A. oryzae was performed, which led to the production of aculene D (96) as a major product together with asperaculane G (97). Thus, AneG performs C-H activation at the C-2 position of asperaculane E (98) and catalyzes a double oxidation to install a C-2 ketone group. The involvement of 97 as pathway intermediate was validated, as the feeding of 97 to ΔaneC restored the production of 55. Overall, the researchers proposed that all the three P450s are necessary for this stepwise demethylation process. AneF oxidizes C-12 to a carboxylate, AneG installs a necessary electro-withdrawing carbonyl group on C-2, while AneD introduces a10-hydroxy group to assist the tautomerization.
Next, the modification enzymes involved in the gene cluster was investigated. Deletion of AneA, which encodes an α-ketoglutarate-dependent dioxygenase, led to the abolishment of the production of 55 and aculene C (99), and accumulation of 56 and 96, and a minor component asperaculane B (100) (shunt product). This suggest that the unsaturation is catalyzed by AneA.
Deletion of aneB and aneE both led to the absence of 55 and 56 and the accumulation of 96 and 100. This result suggests that AneB adenylates L-proline and AneE transfers this activated L-proline derivative to 96 to form 56.
Compared to the downstream pathway products, there are several upstream intermediates containing a 14-hydroxy group. To confirm that the 14-hydroxylated compounds are shunt products, the researchers supplied 96 and 100 to aneC-deleted strain. The former restores the production of aculene but the latter does not, which suggest that the metabolites with 14-hydroxy group are shunt products in the biosynthesis of 55.
Overall, the biosynthetic pathway of 55 initiates with the biosynthesis of 95 by AneC, followed by C-12 carboxylation by AneF to produce asperaculane D (101), and C-10 hydroxylation by AneD, yielding asperaculane E (98). Then AneG installs an electron-withdrawing carbonyl group at the C-2 position, which triggers C-12 decarboxylation to yield the nordaucane skeleton in 96. Next, the unsaturation between C5 and C6 are catalyzed by AneA, yielding 55 and 99. Finally, the adenylation is catalyzed by AneB, followed by transferring the activated L-proline derivatives to 96, yielding 56.
Fig. 30 The gene cluster akc in A. niger. The boundaries of the gene cluster has not been defined.154 |
Fig. 31 Predicted biosynthetic pathway of alkylcitric acids.154 |
The proposed alkylcitrate pathway starts with the production of hexylcitric acid (102) catalyzed by a citrate synthase (akcB) and a fatty acid synthase (akcA and akcD). Next, the dihydroxylation of hexylcitric is catalyzed by 2-methylcitrate dehydratase (akcC) to produce hexylaconitic acid A (103). Finally, a decarboxylation reaction is catalyzed by hexylaconitic acid decarboxylase (hadA), yielding hexylitaconic acid B (104) (Fig. 31).154
This journal is © The Royal Society of Chemistry 2023 |