Next Article in Journal
Treatment of Invasive Candidiasis: A Narrative Review
Next Article in Special Issue
Development of an Improved Carotenoid Extraction Method to Characterize the Carotenoid Composition under Oxidative Stress and Cold Temperature in the Rock Inhabiting Fungus Knufia petricola A95
Previous Article in Journal
Ecological Determinants of Sporotrichosis Etiological Agents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oryzines A & B, Maleidride Congeners from Aspergillus oryzae and Their Putative Biosynthesis

1
University of Bristol, School of Chemistry, Cantock’s Close, Bristol BS8 1TS, UK
2
Leibniz Universität Hannover, Biomolekulares Wirkstoffzentrum und Institut für Organische Chemie, Schneiderberg 38, 30167 Hannover, Germany
3
University of Bristol, School of Chemistry, Cantock’s Close, Bristol BS8 1TS, UK
4
Leibniz Universität Hannover, Biomolekulares Wirkstoffzentrum und Institut für Organische Chemie, Schneiderberg 38, 30167 Hannover, Germany
*
Author to whom correspondence should be addressed.
J. Fungi 2018, 4(3), 96; https://doi.org/10.3390/jof4030096
Submission received: 19 July 2018 / Revised: 8 August 2018 / Accepted: 10 August 2018 / Published: 13 August 2018
(This article belongs to the Special Issue Fungal Polyketides and Other Secondary Metabolites)

Abstract

:
Aspergillus oryzae is traditionally used in East Asia for the production of food and brewing. In addition, it has been developed into a suitable host for the heterologous expression of natural product biosynthetic genes and gene clusters, enabling the functional analysis of the encoded enzymes. A. oryzae shares a 99.5% genome homology with Aspergillus flavus, but their secondary metabolomes differ significantly and various compounds unique to A. oryzae have been reported. While using A. oryzae as a host for heterologous expression experiments we discovered two new metabolites in extracts of A. oryzae M-2-3 with an unusual maleidride backbone, which were named oryzine A and B. Their structures were elucidated by high resolution mass spectrometry (HRMS) and nuclear magnetic resonance (NMR) analysis. Their structural relationships with known maleidrides implied involvement of a citrate synthase (CS) and a polyketide (PKS) or fatty acid synthase (FAS) in their biosynthesis. Analysis of the A. oryzae genome revealed a single putative biosynthetic gene cluster (BGC) consistent with the hypothetical biosynthesis of the oryzines. These findings increase knowledge of the chemical potential of A. oryzae and are the first attempt to link a novel product of this fungus with genomic data.

Graphical Abstract

1. Introduction

Aspergillus oryzae is a filamentous fungus widely used in the food production industry of many East Asian countries to ferment soy beans as basis for soy sauces, miso and fermented bean paste, but also for brewing of various beverages [1]. In addition, enzymes used in food processing, such as α-amylases, aminopeptidases, lactases, lipases or proteases, are isolated in large scale from A. oryzae [2]. Due to its economic importance the fungus has been extensively studied and genetic tools became readily available to manipulate the organism. In combination with its advantageously limited set of native compounds A. oryzae has been successfully developed as a heterologous expression system for fungal secondary metabolites [2,3]. This system has been important in the discovery and functional analysis of various biosynthetic gene clusters (BGC) including xenovulenes, pleuromutilins, maleidrides and aphidicolin [4,5,6,7].
Whole genome sequencing of the RIB40 strain of A. oryzae and other Aspergillus species revealed a 99.5% similarity between the genomes of Aspergillus flavus and A. oryzae. Therefore, it has been suggested that the latter is a domesticated variety of A. flavus, mainly differing in the lack of aflatoxin production [8]. An in-depth chemical comparison between the two species showed striking differences in the metabolites produced. Thus, A. oryzae produces kojic acid, the diterpene indoles aflavinines, the aflatrem precursor 13-desoxypaxilline, the sesquiterpenes parasiticolides and various alkaloids, e.g., ditryptoleucine and oryzamides [9]. Other compounds produced by A. oryzae have been reported recently [10,11,12]. However, identifications of the respective producing strains were solely based on morphological characteristics and ribosomal genes, which are insufficient to clearly differentiate them from A. flavus [13]. Even though the genome sequence of A. oryzae has been available since 2005 [14], only the BGCs of aflatoxin, cyclopiazonic acid, aflatrem, kojic acid, 3, 5-dihydroxybenzoic acid and penicillin have been identified [9,15,16,17]. This is a small number compared with the 32 annotated polyketide synthase (PKS) and 27 annotated non-ribosomal peptide synthetase (NRPS) BGCs of A. oryzae, and leaves much room for new findings [9].
In this study, we report the identification and structure elucidation of two new secondary metabolites of the maleidride family produced by A. oryzae M-2-3 [18]. In addition, a candidate BGC encoded by the genome of A. oryzae is bioinformatically characterized and a putative biosynthetic pathway is presented.

2. Materials and Methods

2.1. General

Optical rotations were determined with an ADP 220 polarimeter at 589 nm and IR spectra were recorded with a Perkin–Elmer FTIR instrument. NMR spectra were recorded on a Varian VNMRS-500 spectrometer (1H 500 MHz, 13C 125 MHz). HR-ESI-MS mass spectra were determined with a Bruker Daltonics Apex IV FT-ICR. HPLC analysis was performed on a Waters 2795HT HPLC system. Detection was achieved by UV between 200 and 400 nm using a Waters 998 diode array detector, and by simultaneous electrospray (ES) mass spectrometry using a Waters ZQ spectrometer detecting between 150 and 600 m/z units. Chromatography (flow rate 1 mL/min) was performed using a Phenomenex Kinetex column (2.6 μ, C18, 100 Å, 4.6 × 100 mm) equipped with a Phenomenex Security Guard precolumn (Luna C5 300 Å). The mobile phase consisted of solvent A (HPLC grade H2O containing 0.05% formic acid) and solvent B (HPLC grade CH3CN containing 0.045% formic acid) subjected to the following gradient: 0 min, 10% B; 10 min, 90% B; 12 min, 90% B; 13 min, 10% B; 15 min, 10% B.

2.2. Fermentation and Extraction

Three liters of A. oryzae M-2-3 pTAex3-tenS(8309)·egfp [19] (30 flasks × 100 mL) were grown in production medium (20 g/L starch, 10 g/L peptone, 2 g/L NaNO3, 2 g/L KCl, 1 g/L K2HPO4, 0.5 g MgSO4.7H2O, 0.01 g FeSO4.7H2O) for 6 days at 28 °C and shaken at 200 rpm. Cells and media were homogenized using a hand-held electric blender and then acidified to pH 4.0 using 37% aqueous HCl. An equal volume of ethyl acetate was added and the mixture was stirred for 10 min. and then vacuum filtered through Whatman no. 1 filter paper. The filtrate was transferred into a separating funnel and shaken vigorously. The organic layer was washed once with concentrated brine solution and then with deionized water. The organic phase was dried (MgSO4), filtered and evaporated to dryness. The crude extract was dissolved in 10% aqueous methanol and defatted by extraction with hexane. The methanolic layer was evaporated to dryness to yield crude extract (239 mg) and analyzed by liquid chromatography mass spectrometry (LCMS).

2.3. Isolation of Oryzine A 3 and Oryzine B 4

The crude extract was dissolved in HPLC grade methanol and subjected to mass-directed preparative HPLC. Purification of compounds was achieved using a Waters mass-directed auto-purification system comprising of a Waters 2767 autosampler, Waters 2545 pump system, a Phenomenex LUNA column (5 μ, C18, 100 Å, 10 × 250 mm) equipped with a Phenomenex Security Guard precolumn (Luna C5 300 Å) eluted at 16 mL/min. Solvents used were: A, HPLC grade H2O + 0.05% formic acid; B, HPLC grade CH3CN + 0.045% formic acid. The post-column flow was split (100:1) and the minority flow was made up with solvent A to 1 mL/min for simultaneous analysis by diode array detector (Waters 2998; Milford, MA, USA), evaporative light scattering (Waters 2424; Milford, MA, USA) and ESI mass spectrometry in positive and negative modes (Waters Quatro Micro; Milford, MA, USA). The gradient started with 40% to 80% B over 15 min, continued for 0.5 min to 95% and ended with 1 min isocratic conditions. Fractions of 24 consecutive runs were combined to yield 9.5 mg of 3 and 9.0 mg of 4.

2.4. Spectroscopic and Spectrometric Data

The corresponding 1D and 2D NMR spectra of 3 and 4 are depicted in the Supplementary Information (Figures S1–S10).

2.4.1. Oryzine A 3

Light brown viscous oil, [α]22D −16.9 (c = 0.23, MeOH); IR (neat): νmax 2932, 2872, 2342, 1761, 1631, 1355, 1190, 1084 cm−1; 1H NMR (CDCl3, 500 MHz) δ = 0.92 (3H, t, J = 7 Hz, H-1), 1.36 (2 H, m, H-2), 1.37 (1 H, m, H-3a), 1.45 (1 H, m, H-3b), 1.65 (1 H, m, H-4a), 1.81 (1 H, m, H-4b), 2.01 (1 H, ddd, J = 12, 10.5, 12 Hz, H-6a), 2.55 (1 H, ddd, J = 9, 6, 12 Hz, H-6b), 3.65 (1H, t, J = 9, 12 Hz, H-7), 4.43 (1 H, m, H-5), 5.95 (1 H, b, H-10a), 6.54 (1 H, b, H-10b), 13C NMR (CDCl3, 125 MHz) δ = 14.4 (C-1), 22.9 (C-2), 27.8 (C-3), 35.5 (C-4), 36.1 (C-6), 45.2 (C-7), 79.6 (C-5), 131.8 (C-10), 136.1 (C-8), 170.1(C-9), 176.4 (C-11). HRESIMS calculated for C11H17O4: 213.1121; observed 213.1130 [M + H]+.

2.4.2. Oryzine B 4

Light color viscous oil, IR (neat): νmax 2962, 2930, 2873, 1736, 1582, 1454, 1045, 878 cm−1; 1H NMR (CDCl3, 500 MHz) δ = 0.92 (3 H, t, J = 7 Hz, H-1), 1.38 (2 H, m, H-2), 1.46 (2 H, m, H-3), 1.71 (1 H, m, H-4a), 1.79 1 H, m, H-4b), 4.99 (1 H, t, J = 6.2 Hz, H-5), 6.79 (1 H, b, H-10a), 7.19 (1 H, b, H-10b), 7.96 (1 H, b, H-6); 13C NMR (CDCl3, 125 MHz) δ = 14.0 (C-1), 22.6 (C-2), 27.3 (C-3), 33.1 (C-4), 80.8 (C-5), 125.0 (C-7), 128.5 (C-8), 133.6 (C-10), 153.5 (C-6), 169.5 (C-9), 171.7 (C-11). HRESIMS calculated for C11H14O4Na: 233.0784; observed 233.0799 [M + Na]+.

2.5. Oryzine Biosynthetic Gene Cluster Prediction

The prediction of the BGC responsible for oryzine production is based on the genome information available for A. oryzae RIB40 (BioProject: PRJNA28175). The producing strain of the oryzines (M-2-3) is an auxotrophic mutant derived from A. oryzae RIB203, which is genetically very close to RIB40 (https://www.nrib.go.jp/data/asp/strain.html). Due to annotation gaps a part of the genome was re-annotated. Gene prediction was accomplished by applying AUGUSTUS version 3.3. [20] on the A. oryzae RIB40 chromosome 8 (RefSeq: NC_036442.1) using the precomputed gene model for this species. For comparison purpose, in addition, GeneMark-ES version 4.3.3. [21] was applied by using the parameters for self-training (–ES) and fungal genomes (–fungus). Clustered genes in A. oryzae related to the byssochlamic acid BGC from Byssochlamys fulva were identified with the MultiGeneBlast architecture search [22] using bfpks1, bfl2 and bfl3 [6] as templates. Homologies of the putative oryzine BGC with related BGCs from A. flavus NRRL 3357 and B. fulva were analyzed with the Artemis comparison tool (ACT) [23].

3. Results

As part of a larger project involving the use of A. oryzae M-2-3 [18] as a host for the heterologous expression of engineered fungal polyketide synthase genes, we were interested in the biosynthesis of pretenellin A 1. This compound is produced by the tenellin synthetase (TENS) which is a hybrid polyketide synthase non-ribosomal peptide synthetase (PKS-NRPS) which works with a trans-acting enoyl reductase, encoded by tenC or dmbC [24,25]. We engineered a system containing a truncated tenS gene, lacking the sequence encoding the NRPS, together with dmbC [19]. We predicted that this would produce the tenellin polyketide 2 with nominal mass 210 Da. Comparison of transformed (Figure 1C,D) vs WT A. oryzae M-2-3 (Figure 1A,B) revealed the presence of a new compound (RT 30.3 min) in the transformed strain with the expected nominal mass (Figure 1D, [M + H]+ 211.1). This was accompanied by a compound (RT 29.0 min, [M + H]+ 213.1) which is present in both WT and transformed strains, although at significantly lower titer in the WT strain (<1%). Both compounds were purified from extracts of the transformed A. oryzae strain by mass-directed HPLC. Full structure elucidation (vide infra) revealed these two compounds to be unrelated to pretenellin A 1, showing that the original experiment had failed.

3.1. Structure Elucidation

Oryzine A 3 was obtained as light brown viscous oil. The HRESIMS data of 3 indicated a molecular formula of C11H17O4 (calc. 213.1121, obs. 213.1130, [M + H]+). The carbon spectra (125 MHz) revealed the presence of eleven carbons, which included two carbonyl groups at δC 176.4 (C-11) and 170.1 (C-9), four methylene groups in the alkane region of the spectra and one methyl carbon at δC 14.4 (C-1, Figure 2). The HSQC (heteronuclear single quantum coherence spectroscopy) spectrum revealed 8 protonated carbons. The 1H NMR spectrum showed two distinct doublets in the alkene region at δH 5.95 and δH 6.54 which were assigned to two geminal methylene protons attached to C-10. One of these methylene protons shows HMBC (heteronuclear multiple bond correlation) correlation to an alkene carbon at δC 136.1 which was assigned C-8. Both geminal methylene protons (H-10) show HMBC correlations to C-9 and a methine at δC 45.2 (C-7). The triplet at δH 3.65 was assigned to a methine proton (H-7) which showed an HMBC correlation to C-11, C-9 and an alkene carbon at δC 136.1 (C-8). H-7 shows COSY (correlation spectroscopy) coupling to two geminal methylene protons at δH 2.01 and 2.55 (H-6). These, in turn, show COSY and HMBC connections to a methine proton at δH 4.43 (C-5) attached to an oxygen atom. The chemical shift of δC 79.6 for C-5 supports its link to an oxygen atom. This data suggested the presence of a pyrone with a carboxylic acid at C-7. The methine protons at δH 4.43 (H-5) and methylene protons at δH 2.01 and 2.55 (H-6) show COSY connections to two other geminal protons at δH 1.65 and 1.81 (H-4). Two multiplets at δH 1.37 and 1.46 were assigned to geminal methylene protons (3-CH2) which displays a COSY coupling to one of the methylenes at δH 1.8 (H-4). The broad signal at δH 1.35–1.39 was assigned to methylene protons (2-CH2) showing correlations to δC 27.8 (C-3) in the HMBC spectrum and to the methyl at δH 0.92 (H-1). Both HMBC and COSY correlations show an aliphatic chain (C1-C4) attached to the pyran at C-5 methine.
The relative orientation of the protons at the two stereocentres C-5 and C-7 was established using the coupling constants of H-5 and H-7 with their adjacent methylene protons at δH 2.01 (H-6a) and δH 2.55 (H-6b, Figure 2). A J value of 6 Hz between H-5 and H-6b indicated this pair of protons to be gauche, while the J value of 10 Hz between H-5 and H-6a suggesting this pair as anti and thus axial/axial. The methine proton H-7 displayed a J value of 12 Hz with H-6a and 9 Hz with H-6b. These values suggest that H-7 is also axial.
The structure of 4 was assigned using HRMS and NMR analysis and by comparison to 3. The chemical formula is C11H14O4 (HRMS, observed 233.0799; calculated 233.0784 for [M + Na]+). The 13C NMR displayed the presence of 11 carbons including two carbonyls at δC 169.5 (C-9) and 171.7 (C-11), a methine carbon at δC 80.8 (C-5) attached to an oxygen atom, a methyl group at δC 14.0 (C-1), three methylene groups and four olefinic carbons at δC 125.0 (C-7), 128.5 (C-8), 133.6 (C-10) and 153.5 (C-6, Figure 2). The HSQC spectrum showed seven protonated carbons. In the 1H NMR, the broad doublet signal downfield at δH 7.96 was assigned to the methine (H-6) attached to the olefinic carbon at δC 153.5 (C-6) in the HSQC spectrum. In the HMBC spectrum, H-6 shows linkage to methine C-5 and to the olefinic carbon C-7. This confirmed the presence of an olefin between C-6 and C-7. H-6 also shows HMBC correlations with the olefin carbons C-8 and C-11. The two broad resonances at δH 6.79 and 7.19 were assigned to geminal methylene protons attached to C-10 showing correlations to the methine protons H-5 and H-6 in the COSY, and HMBC correlations to the olefin carbon C-7 and to the C-9 carbonyl. The triplet signal at δH 4.99 was assigned to the H-5 methine, displaying HMBC correlations with the methylene groups at δC 27.3 (C-3) and 33.1 (C-4) and to olefin carbons C-7 and C-6. The two multiplet signals at δH 1.71 and 1.79 were assigned to geminal methylene protons (H-4) correlated to C-2 and C-3 in the HMBC. The methyl signal at δH 0.92 (H-1) showed HMBC correlations to the methylene groups at C-2 and C-3. This compound was thus assigned the same carbon skelton as oryzine A 3, but unsaturated between C-6 and C-7.

3.2. Identification of A Putative Oryzine Biosynthetic Gene Cluster

The structures of the oryzines indicate a biosynthetic relationship to the maleidrides [6,26,27]. In general, the biosynthesis of maleidrides is characterized by condensation of an acyl CoA thiolester of variable chain-length with oxaloacetic acid, catalyzed by a citrate synthase (CS). In addition, a 2-methyl citrate dehydratase (2MCD) is required to catalyze subsequent steps [6]. Therefore, the A. oryzae genome was analyzed for a BGC containing genes encoding those enzymes. MultiGeneBlast [22] analysis (Figure S11) indicated the presence of a single BGC that contains genes encoding a CS and a 2MCD. The respective genes had no PKS in their proximity but were instead surrounded by α and β subunits of a fungal FAS. To verify that this gene cluster is the only possible candidate, a manual search of the A. oryzae genome for CS genes was also conducted. A total of five homologues of CS were identified on the A. oryzae genome by using the CS (Bfl2) of the byssochlamic acid BGC (ANF07286) from Byssochlamys fulva as template. The genes in the proximity to these five targets were analyzed in detail. Four of these genes (amino acid sequence identity with Bfl2 between 55.5% and 11.2%) were not close to any recognized types of biosynthetic genes. The fifth gene (37.7% amino acid sequence identity) was the one predicted by MultiGeneBlast. In addition to the 2MCD and FAS subunits various other biosynthetic genes including oxygenases, lactonases and a decarboxylase were found in this region of the genome. Hence, this BGC was considered as the most likely candidate for oryzine biosynthesis (Figure 3).

3.3. Characteristics of the Putative Oryzine BGC

The putative oryzine BGC has a size of around 47 kb and contains 18 predicted genes (Table 1, Table S1). The two largest genes (oryfasA, oryfasB) encode the α- and β-subunits of a fungal fatty acid synthase (FAS). Additional copies of FAS genes can be found within the genome of A. oryzae. OryfasA is located at the downstream part of the BGC next to the citrate synthase (oryE) and a putative dehydrogenase (oryD). In contrast, oryfasB is found at the upstream end of the BGC flanked by a 2MCD (oryR) and a gene encoding a protein with high similarity to sterigmatocystin P450 monooxygenase (oryQ). The short amino acid sequence of OryQ (170 aa) implies a truncation, probably hampering its functionality. Between the two FAS encoding genes additional biosynthetic genes can be identified. OryG is homologous to the α -ketoglutarate-dependent taurine dioxygenase (P37610, 31% identity) and OryH as well as OryL contain a domain that belongs to the lactonase superfamily. OryM encodes a protein related to aconitate decarboxylases, whereas oryP is most likely encodes an acyl-CoA ligase. In addition, multiple putative transporters (oryC, oryF, oryN) and a putative transcription factor (oryP) can be found. The genes oryJ and oryK cannot be assigned to any known biosynthetic family.

3.4. Proposed Biosynthesis of the Oryzines

With the set of identified biosynthetic genes in hand, a putative biosynthetic pathway can be proposed (Figure 4). The two subunits of the fungal FAS (OryfasA, OryfasB) probably form octenoic acid 5. This fatty acid is most likely activated by the acyl-CoA ligase OryP to give octenyl CoA 6 before the CS (OryE) acts to catalyze condensation with oxaloacetate to form tricarboxylic acid 7. At this stage, two pathways appear feasible. In the first, decarboxylation and concomitant dehydration (OryM), followed by tautomerization gives the diene 9. Reduction of this could give the known metabolite piliformic acid 10 from the fungus Xylaria mali which has been shown to be derived from FAS-derived octanoic acid and oxaloacetate by feeding studies [28]. On the pathway to 3 and 4, however, hydroxylation of the diene 9 (OryG) and lactonisation (OryH or OryL) could give oryzine B 4 directly. Final enoyl reduction (OryD) would then give Oryzine A 3. An alternative pathway branching from 7 could involve epoxidation to 11 (OryG) followed by ring closure (OryH or OryL). The resulting diol 12 could be a substrate for the MCDH to give unsaturated alcohol 13. Finally, concerted decarboxylation and dehydration (OryM) could give the observed 4. We do not yet have enough experimental evidence to distinguish these possibilities, but the first pathway is preferred due to the proximity to the known metabolite piliformic acid 10.

4. Discussion

It is clear that our original genetic engineering experiment in A. oryzae M-2-3 did not give the anticipated result. However, the transformed organism produced two compounds, oryzines A 3 and B 4, in titers significantly higher than untransformed strains of A. oryzae. As various other transformants did not show increased production of oryzines it can be speculated that the regulation of gene expression within the oryzine BGC was influenced by the random insertion of the heterologous expressed tenS gene into the A. oryzae genome. Related compounds such as piliformic acid 10 are known to be derived from fatty acids and oxaloacetate [28]; while the furofuranones 4-epi-ethiosiloide 14, sporothriolide 15 and discosiolide 16, are also thought to be derived from fatty acids and oxaloacetate [29,30]. These compounds differ from 3 and 10 by having differing carbon-chain lengths, and by a different lactonisation pattern. It is therefore likely that the oryzines, piliformic acid 10 and the furofuranones 1416 are assembled by similar biosynthetic machinery. However, to date, none of these compounds have been linked to any BGC. The presence of a fatty acyl and an oxaloacetyl moiety is also known in the maleidride class of secondary metabolites where a key citrate synthase-type enzyme is responsible for their formation. We thus suspected that the oryzines are a class of hitherto unrecognised maleidrides and the biosynthesis should also require a citrate synthase step.
Based on the BGC prediction, the first steps of oryzine biosynthesis are likely to be similar to those of the elucidated byssochlamic acid (BA 17) biosynthesis (Figure 5) [6]. The main difference is the mode of production of the alkenoic acid: in case of BA 17 this is produced by a polyketide synthase; but in the case of 3 and 4 the alkanoic acid is probably produced by a fatty acid synthase. FAS components of fungal BGC are known, but their essential role for the biosynthesis of secondary metabolites has only been partially demonstrated. In case of the aflatoxins, sterigmatocystin and dothistromin the FAS subunits provide the initial hexanoyl starter unit for the PKS, which forms norsolorinic acid 18 [31,32] (Figure 5). A homology search of FAS genes in the genome of A. oryzae indicated the presence of four additional copies of both subunits. Only one of these gene pairs had no biosynthetic genes in their proximity, therefore the respective proteins are likely part of primary metabolism. The second FAS copy belongs to the previously identified aflatoxin BGC [31]. The third homologous pair is nested inside an NRPS gene cluster that shows high similarity to the aspercryptin BGC from A. nidulans which is known to require octanoic acid during its biosynthesis [33]. The last FAS is part of an azaphilone-type BGC similar to those identified from various Monascus and Talaromyces (Penicillium) species [34]. The high diversity of FAS containing secondary metabolite gene clusters in A. oryzae is striking and highlights the importance of such proteins for natural product biosynthesis.
The BA 17 system requires a hydrolase encoded by bfl1, presumably to release the polyketide from the PKS at the end of its synthesis [6]. However in the case of fungal FAS a dedicated acyl transferase, which is an integral component, releases the fatty acid at the end of biosynthesis. Interestingly it has been shown that these AT domains can be engineered to release short fatty acids [35]. Consistent with this idea is the fact that the proposed oryzine cluster does not contain a gene homologous to bfl1. Similarly, the BA BGC contains genes which encode proteins responsible for the dimerisation (to form nonadrides such as 17) which are not present in the proposed oryzine BGC; and oryzine biosynthesis requires oxidative steps and thus the oryzine BGC has encoded oxygenases not present in the BA BGC (Figure S12).
It is likely that maleic acid anhydride intermediates such as 19 (Figure 5) are not formed during oryzine biosynthesis due to the presence of an aconitate decarboxylase, which is not present in the other known maleidride gene clusters. The loss of a carboxyl functionality at C-10 would prevent anhydride formation, which explains the structural differences between the oryzines and the maleic acid anhydrides. In other systems one could imagine oxidations and lactonisations with different regioselectivities leading to the furofuranones 1416.
Even though only a limited set of biosynthetic steps is necessary to explain the assembly of the oryzines, the BGC seems rather large. However, in order to correctly define the borders of the cluster transcription studies should be conducted. It is also questionable if all predicted genes encode functional proteins. A BGC comparison with A. flavus revealed the presence of a homologous BGC with >99.5% similarity (Figure S12), highlighting the potential of A. flavus to produce oryzines as well. Whether these compounds have been overlooked in the past by researchers, or are not produced under the tested conditions or whether small mutations in the A. flavus BGC abolished production of 3 and 4 or closely related compounds must also be evaluated in the future.

5. Conclusions

The oryzines isolated from Aspergillus oryzae constitute a novel class of maleidrides, which are structurally related to furofuranones 1416. In our previous work we have proven the links between fungal gene clusters by either knockout or heterologous expression in engineered strains of A. oryzae. Here, however, the heterologous expression strategy is obviated by the presence of the genes of interest in the expression host itself; while knockout strategies were frustrated by the low titer of 3 and 4 in untransformed A. oryzae. However, comparison of the BGCs of A. oryzae with the byssochlamic acid gene cluster from B. fulva revealed a single candidate cluster containing FAS, citrate synthase, dehydratase, decarboxylase and oxygenase genes, and no other likely BGCs in the A. oryzae genome. The presence of a FAS and a CS in a single cluster and their link to a final product is shown for the first time, adding to our knowledge about fungal biosynthetic pathways. These observations also allow the development of a firm hypotheses for the biosynthesis (and hence the composition of biosynthetic gene clusters) of the furofuranones and related compounds such as piliformic acid 10. Further progress, however, will have to await full genome sequencing of the organisms responsible for the biosynthesis of these compounds. Future work will involve the development of new methods to understand such clusters.

Supplementary Materials

The following are available online at https://www.mdpi.com/2309-608X/4/3/96/s1: Figures S1–S10, 1D and 2D NMR data for 3 and 4; Figure S11, MultiGeneBlast (architecture search) of the byssochlamic acid biosynthetic genes bfpks1, bfl2 and bfl3 against the genomes of Aspergillus oryzae RIB40 and A. flavus NRRL3357; Figure S12, Artemis comparison tool (ACT) analysis between the putative oryzine biosynthetic gene cluster (BGC) of A. flavus (A) and A. oryzae (B) and the byssochlamic acid BGC of Byssochlamys fulva (C); Table S1, Predicted genes and corresponding nucleotide as well as protein sequences of the re-annotated oryzine BGC from A. oryzae RIB40.

Author Contributions

Conceptualization, R.J.C. and T.J.S.; Investigation, Z.W. and E.K.; Writing-Original Draft Preparation, R.J.C. and E.K.; Supervision, R.J.C. and T.J.S.; Project Administration, R.J.C. and T.J.S.; Funding Acquisition, R.J.C. and Z.W.

Funding

This research was funded by EPSRC (Engineering and Physical Sciences research council) grant number EP/F066104/1 and DFG (Deutsche Forschungsgemeinschaft) grant number CO 1328/4-1.

Acknowledgments

The authors thank Jenny Rose and Walid Bakeer for technical assistance. Zahida Wasil Thanks the Government of Pakistan (Higher Education Commission) for funding. We are grateful to Daniel Wibberg for the re-annotation of the A. oryzae genomic data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Park, H.S.; Jun, S.C.; Han, K.H.; Hong, S.B.; Yu, J.H. Diversity, Application and synthetic biology of industrially important aspergillus fungi. Adv. Appl. Microbiol. 2017, 100, 161–202. [Google Scholar] [CrossRef] [PubMed]
  2. Pariza, M.W.; Johnson, E.A. Evaluating the safety of microbial enzyme preparations used in food processing: Update for a new century. Regul. Toxicol. Pharmacol. 2001, 33, 173–186. [Google Scholar] [CrossRef] [PubMed]
  3. He, Y.; Wang, B.; Chen, W.; Cox, R.J.; He, J.; Chen, F. Recent advances in reconstructing microbial secondary metabolites biosynthesis in Aspergillus spp. Biotechnol. Adv. 2018, 36, 739–783. [Google Scholar] [CrossRef] [PubMed]
  4. Schor, R.; Schotte, C.; Wibberg, D.; Kalinowski, J.; Cox, R.J. Three previously unrecognised classes of biosynthetic enzymes revealed during the production of xenovulene A. Nat. Commun. 2018, 9, 1963. [Google Scholar] [CrossRef] [PubMed]
  5. Bailey, A.M.; Alberti, F.; Kilaru, S.; Collins, C.M.; De Mattos-Shipley, K.; Hartley, A.J.; Hayes, P.; Griffin, A.; Lazarus, C.M.; Cox, R.J.; et al. Identification and manipulation of the pleuromutilin gene cluster from Clitopilus passeckerianus for increased rapid antibiotic production. Sci. Rep. 2016, 6, 1–11. [Google Scholar] [CrossRef] [PubMed]
  6. Williams, K.; Szwalbe, A.J.; Mulholland, N.P.; Vincent, J.L.; Bailey, A.M.; Willis, C.L.; Simpson, T.J.; Cox, R.J. Heterologous production of fungal maleidrides reveals the cryptic cyclization involved in their biosynthesis. Angew. Chem.-Int. Ed. 2016, 55, 6784–6788. [Google Scholar] [CrossRef] [PubMed]
  7. Fujii, R.; Minami, A.; Tsukagoshi, T.; Sato, N.; Sahara, T.; Ohgiya, S.; Gomi, K.; Oikawa, H. Total biosynthesis of diterpene aphidicolin, a specific inhibitor of dna polymerase α: Heterologous expression of four biosynthetic genes in Aspergillus oryzae. Biosci. Biotechnol. Biochem. 2011, 75, 1813–1817. [Google Scholar] [CrossRef] [PubMed]
  8. Machida, M.; Yamada, O.; Gomi, K. Genomics of Aspergillus oryzae: Learning from the history of koji mold and exploration of its future. DNA Res. 2008, 15, 173–183. [Google Scholar] [CrossRef] [PubMed]
  9. Rank, C.; Klejnstrup, M.L.; Petersen, L.M.; Kildgaard, S.; Frisvad, J.C.; Held Gotfredsen, C.; Ostenfeld Larsen, T. Comparative chemistry of Aspergillus oryzae (RIB40) and A. flavus (NRRL 3357). Metabolites 2012, 2, 39–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Shaaban, M.; El-Metwally, M.M.; Nasr, H. A new diketopiperazine alkaloid from Aspergillus oryzae. Nat. Prod. Res. 2014, 28, 86–94. [Google Scholar] [CrossRef] [PubMed]
  11. Ren, R.; Chen, C.-J.; Hu, S.-S.; Ge, H.-M.; Zhu, W.-Y.; Tan, R.-X.; Jiao, R.-H. Drimane Sesquiterpenoids from the Aspergillus oryzae QXPC-4. Chem. Biodivers. 2015, 12, 371–379. [Google Scholar] [CrossRef] [PubMed]
  12. Liu, L.; Bao, L.; Wang, L.; Ma, K.; Han, J.; Yang, Y.; Liu, R.; Ren, J.; Yin, W.; Wang, W.; et al. Asperorydines A-M: Prenylated Tryptophan-Derived Alkaloids with Neurotrophic Effects from Aspergillus oryzae. J. Org. Chem. 2018, 83, 812–822. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, C.Z.; Liou, G.Y.; Yuan, G.F. Comparison of Aspergillus flavus and Aspergillus oryzae by amplified fragment length polymorphism. Bot. Bull. Acad. Sin. 2004, 45, 61–68. [Google Scholar]
  14. Machida, M.; Asai, K.; Sano, M.; Tanaka, T.; Kumagai, T.; Terai, G.; Kusumoto, K.I.; Arima, T.; Akita, O.; Kashiwagi, Y.; et al. Genome sequencing and analysis of Aspergillus oryzae. Nature 2005, 438, 1157–1161. [Google Scholar] [CrossRef] [PubMed]
  15. Seshime, Y.; Juvvadi, P.R.; Kitamoto, K.; Ebizuka, Y.; Nonaka, T.; Fujii, I. Aspergillus oryzae type III polyketide synthase CsyA is involved in the biosynthesis of 3,5-dihydroxybenzoic acid. Bioorganic Med. Chem. Lett. 2010, 20, 4785–4788. [Google Scholar] [CrossRef] [PubMed]
  16. Marui, J.; Ohashi-Kunihiro, S.; Ando, T.; Nishimura, M.; Koike, H.; Machida, M. Penicillin biosynthesis in Aspergillus oryzae and its overproduction by genetic engineering. J. Biosci. Bioeng. 2010, 110, 8–11. [Google Scholar] [CrossRef] [PubMed]
  17. Marui, J.; Yamane, N.; Ohashi-Kunihiro, S.; Ando, T.; Terabayashi, Y.; Sano, M.; Ohashi, S.; Ohshima, E.; Tachibana, K.; Higa, Y.; et al. Kojic acid biosynthesis in Aspergillus oryzae is regulated by a Zn(II)2Cys6transcriptional activator and induced by kojic acid at the transcriptional level. J. Biosci. Bioeng. 2011, 112, 40–43. [Google Scholar] [CrossRef] [PubMed]
  18. Gomi, K.; Iimura, Y.; Hara, S. Integrative Transformation of Aspergillus oryzae with a Plasmid Containing the Aspergillus nidulans argB Gene. Agric. Biol. Chem. 1987, 51, 2549–2555. [Google Scholar] [CrossRef]
  19. Bakeer, W. Overexpression of Secondary Metabolism Genes from Magnaporthe grisea and Beauveria bassiana. Ph.D. Thesis, University of Bristol, Bristol, UK, 2011. [Google Scholar]
  20. Stanke, M.; Steinkamp, R.; Waack, S.; Morgenstern, B. Augustus: A web server for gene finding in eukaryotes. Nucleic Acids Res. 2004, 32, W309–W312. [Google Scholar] [CrossRef] [PubMed]
  21. 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]
  22. Medema, M.; Takano, E.; Breitling, R. Detecting Sequence Homology at the Gene Cluster Level with MultiGeneBlast. Mol. Biol. Evol. 2013, 30, 1218–1223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Carver, T.J.; Rutherford, K.M.; Berriman, M.; Rajandream, M.A.; Barrell, B.G.; Parkhill, J. ACT: The Artemis comparison tool. Bioinformatics 2005, 21, 3422–3423. [Google Scholar] [CrossRef] [PubMed]
  24. Fisch, K.M.; Bakeer, W.; Yakasai, A.A.; Song, Z.; Pedrick, J.; Wasil, Z.; Bailey, A.M.; Lazarus, C.M.; Simpson, T.J.; Cox, R.J. Rational Domain Swaps Decipher Programming in Fungal Highly Reducing Polyketide Synthases and Resurrect an Extinct Metabolite. J. Am. Chem. Soc. 2011, 133, 16635–16641. [Google Scholar] [CrossRef] [PubMed]
  25. Heneghan, M.N.; Yakasai, A.A.; Williams, K.; Kadir, K.A.; Wasil, Z.; Bakeer, W.; Fisch, K.M.; Bailey, A.M.; Simpson, T.J.; Cox, R.J.; et al. The programming role of trans-acting enoyl reductases during the biosynthesis of highly reduced fungal polyketides. Chem. Sci. 2011, 2, 972–979. [Google Scholar] [CrossRef]
  26. Szwalbe, A.; Williams, K.; O’Flynn, D.; Bailey, A.; Mulholland, N.; Vincent, J.; Willis, C.; Cox, R.; Simpson, T. Novel nonadride, heptadride and maleic acid metabolites from the byssochlamic acid producer Byssochlamys fulva IMI 40021–an insight into the biosynthesis of maleidrides. Chem. Commun. 2015, 51, 17088–17091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Williams, K.; Szwalbe, A.; Dickson, C.; Desson, T.; Mulholland, N.; Vincent, J.; Clough, J.; Bailey, A.; Butts, C.; Willis, C.; et al. Genetic and chemical characterisation of the cornexistin pathway provides further insight into maleidride biosynthesis. Chem. Commun. 2017, 53, 7965–7968. [Google Scholar] [CrossRef] [PubMed]
  28. Chesters, N.C.J.E.; O’Hagan, D. Biosynthesis of the fungal metabolite, piliformic acid (2-hexylidene-3-methylsuccinic acid). J. Chem. Soc. Perkin. Trans. 1 1997, 827–834. [Google Scholar] [CrossRef]
  29. Krohn, K.; Ludewig, K.; Aust, H.-J.; Dräger, S.; Schulz, B. Biologically active metabolites from fungi. 3. Sporothriolide, discosiolide, and 4-epi-ethisolide New furofurandiones from Sporothrix sp., Discosia sp., and Pezicula livida. J. Antibiot. (Tokyo) 1994, 47, 113–118. [Google Scholar] [CrossRef] [PubMed]
  30. Surup, F.; Kuhnert, E.; Lehmann, E.; Heitkämper, S.; Hyde, K.D.; Fournier, J.; Stadler, M. Sporothriolide derivatives as chemotaxonomic markers for Hypoxylon monticulosum. Mycology 2014, 5, 110–119. [Google Scholar] [CrossRef] [PubMed]
  31. Hitchman, T.S.; Schmidt, E.W.; Trail, F.; Rarick, M.D.; Linz, J.E.; Townsend, C.A. Hexanoate synthase, a specialized type I fatty acid synthase in aflatoxin B1biosynthesis. Bioorg. Chem. 2001, 29, 293–307. [Google Scholar] [CrossRef] [PubMed]
  32. Mahanti, N.; Bhatnagar, D.; Cary, J.W.; Joubran, J.; Linz, J.E. Structure and function of fas-1A, a gene encoding a putative fatty acid synthetase directly involved in aflatoxin biosynthesis in Aspergillus parasiticus. Appl. Environ. Microbiol. 1996, 62, 191–195. [Google Scholar] [PubMed]
  33. Chiang, Y.; Ahuja, M.; Oakley, C.; Entwistle, R.; Asokan, A.; Zutz, C.; Wang, C.; Oakley, B. Development of Genetic Dereplication Strains in Aspergillus nidulans Results in the Discovery of Aspercryptin. Angew. Chem. Int. Ed. 2015, 55, 1662–1665. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, W.; Chen, R.; Liu, Q.; He, Y.; He, K.; Ding, X.; Kang, L.; Guo, X.; Xie, N.; Zhou, Y.; et al. Orange, red, yellow: Biosynthesis of azaphilone pigments in Monascus fungi. Chem. Sci. 2017, 8, 4917–4925. [Google Scholar] [CrossRef] [PubMed]
  35. Zhu, Z.; Zhou, Y.; Krivoruchko, A.; Grininger, M.; Zhao, Z.; Nielsen, J. Expanding the product portfolio of fungal type I fatty acid synthases. Nat. Chem. Biol. 2017, 13, 360–362. [Google Scholar] [CrossRef] [PubMed]
Figure 1. LCMS analysis of extracts of A. oryzae. (A) WT A. oryzae M-2-3, extracted ion chromatogram (EIC+) for m/z 213.1 (oryzine B 4, [M + H]+); (B) WT A. oryzae M-2-3, EIC+ for m/z 211.1 (oryzine A 3, [M + H]+); (C) A. oryzae M-2-3 pTAex3-tenS(8309)·egfp EIC+ for m/z 213.1 (oryzine B 4, [M + H]+); (D) A. oryzae M-2-3 pTAex3-tenS(8309)·egfp EIC+ for m/z 211.1 (oryzine A 3, [M + H]+). * unrelated compound.
Figure 1. LCMS analysis of extracts of A. oryzae. (A) WT A. oryzae M-2-3, extracted ion chromatogram (EIC+) for m/z 213.1 (oryzine B 4, [M + H]+); (B) WT A. oryzae M-2-3, EIC+ for m/z 211.1 (oryzine A 3, [M + H]+); (C) A. oryzae M-2-3 pTAex3-tenS(8309)·egfp EIC+ for m/z 213.1 (oryzine B 4, [M + H]+); (D) A. oryzae M-2-3 pTAex3-tenS(8309)·egfp EIC+ for m/z 211.1 (oryzine A 3, [M + H]+). * unrelated compound.
Jof 04 00096 g001
Figure 2. Chemical structures of pretenellin A 1 and the predicted product 2 of the tenellin PKS. Actual structures of oryzine A 3 and B 4.
Figure 2. Chemical structures of pretenellin A 1 and the predicted product 2 of the tenellin PKS. Actual structures of oryzine A 3 and B 4.
Jof 04 00096 g002
Figure 3. Putative biosynthetic gene cluster of the oryzines.
Figure 3. Putative biosynthetic gene cluster of the oryzines.
Jof 04 00096 g003
Figure 4. Proposed biosynthetic pathway for oryzine A 3 and B 4.
Figure 4. Proposed biosynthetic pathway for oryzine A 3 and B 4.
Jof 04 00096 g004
Figure 5. Compounds biosynthetically related to 3.
Figure 5. Compounds biosynthetically related to 3.
Jof 04 00096 g005
Table 1. Predicted genes of the oryzine BGC and their proposed function.
Table 1. Predicted genes of the oryzine BGC and their proposed function.
GeneGenBank Acc. No.Proposed FunctionClosest UniProtKB Hit (Query Coverage/Identity)
oryCXP_023094066TransporterSugar transporter STL1, Saccharomyces cerevisiae, P39932 (94%/33%)
oryD-DehydrogenaseHistidinol dehydrogenase 1, Bacillus clausii, Q5WIU9 (99%/56%)
oryfasAXP_001827206Fatty acid synthase subunit alphaFatty acid synthase subunit alpha, Aspergillus nidulans, Q5B7V0 (99%/46%)
oryEXP_001827205Citrate synthaseCitrate synthase, Dictyostelium discoideum, Q86AV6 (97%/38%)
oryFXP_001827204TransporterEfflux pump VrtL, Penicillium aethiopicum, D7PHY8 (95%/39%)
oryGXP_001827203Alpha-ketoglut-
-arate dependent dioxygenase
Alpha-ketoglutarate-dependent sulfonate dioxygenase, Saccharomyces cerevisiae, Q12358 (88%/32%)
oryHXP_023094065LactonaseGluconolactonase, Zymomonas mobilis, Q01578 (52%/28%)
oryJXP_023094064-Lysine--tRNA ligase, Bacillus cereus, B7JK84 (24%/30%)
oryKXP_001827200-Ankyrin repeat and KH domain-containing protein mask, Drosophila melanogaster, Q9VCA8 (36%/27%)
oryLXP_023094061LactonaseGluconolactonase, Zymomonas mobilis, Q01578 (66%/25%)
oryMXP_001827196Aconitate decarboxylaseCis-aconitate decarboxylase, Aspergillus terreus, Q0C8L3 (98%/55%)
oryNXP_001827195TransporterCitrinin biosynthesis cluster MFS transporter Mrr1, Monascus ruber, A0A161CLJ6 (93%/48%)
oryOXP_023094060Transcription factorUncharacterized transcriptional regulatory protein C530.05, Schizosaccharomyces pombe, O59741 (24%/40%)
oryP-Acyl-CoA ligaseAcyl-CoA synthetase family member 3, Homo sapiens, Q4G176 (92%/23%)
oryQ-P450 monooxygenase (truncated)Probable sterigmatocystin biosynthesis P450 monooxygenase StcF, Aspergillus nidulans, Q12609 (82%/49%)
oryfasBXP_001827193Fatty acid synthase subunit betaFatty acid synthase subunit beta, Yarrowia lipolytica, P34229 (99%/40%)
oryRXP_0018271922-Methylcitrate dehydratase2-methylcitrate dehydratase, Escherichia coli, P77243 (96%/49%)

Share and Cite

MDPI and ACS Style

Wasil, Z.; Kuhnert, E.; Simpson, T.J.; Cox, R.J. Oryzines A & B, Maleidride Congeners from Aspergillus oryzae and Their Putative Biosynthesis. J. Fungi 2018, 4, 96. https://doi.org/10.3390/jof4030096

AMA Style

Wasil Z, Kuhnert E, Simpson TJ, Cox RJ. Oryzines A & B, Maleidride Congeners from Aspergillus oryzae and Their Putative Biosynthesis. Journal of Fungi. 2018; 4(3):96. https://doi.org/10.3390/jof4030096

Chicago/Turabian Style

Wasil, Zahida, Eric Kuhnert, Thomas J. Simpson, and Russell J. Cox. 2018. "Oryzines A & B, Maleidride Congeners from Aspergillus oryzae and Their Putative Biosynthesis" Journal of Fungi 4, no. 3: 96. https://doi.org/10.3390/jof4030096

APA Style

Wasil, Z., Kuhnert, E., Simpson, T. J., & Cox, R. J. (2018). Oryzines A & B, Maleidride Congeners from Aspergillus oryzae and Their Putative Biosynthesis. Journal of Fungi, 4(3), 96. https://doi.org/10.3390/jof4030096

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