Zhongshu
Song‡
a,
Walid
Bakeer‡
ab,
James W.
Marshall
a,
Ahmed A.
Yakasai
a,
Rozida Mohd
Khalid
a,
Jerome
Collemare
c,
Elizabeth
Skellam
d,
Didier
Tharreau
e,
Marc-Henri
Lebrun
fg,
Colin M.
Lazarus
h,
Andrew M.
Bailey
h,
Thomas J.
Simpson
a and
Russell J.
Cox
*ad
aSchool of Chemistry, University of Bristol, Cantock's Close, Bristol, BS8 1TS, UK
bMicrobiology Department, Faculty of Pharmacy, Beni Suef University, Egypt
cUMR1345, IRHS-INRA, 49071 Beaucouzé Cedex, France
dInstitute for Organic Chemistry, Leibniz University of Hannover, Schneiderberg 1B, 30167, Hannover, Germany. E-mail: russell.cox@oci.uni-hannover.de
eUMR BGPI, CIRAD, Campus International de Baillarguet, 34398 Montpellier Cedex 5, France
fUR 1290 BIOGER-CPP, INRA, Campus AgroParisTech, 78850 Thiverval-Grignon, France
gUMR 5240 MAP, CNRS, UCB, INSA, Bayer CropScience, 69263 Lyon Cedex 09, France
hSchool of Biological Sciences, University of Bristol, 24 Tyndall Avenue, Bristol BS8 1TQ, UK
First published on 1st June 2015
The ACE1 and RAP1 genes from the avirulence signalling gene cluster of the rice blast fungus Magnaporthe oryzae were expressed in Aspergillus oryzae and M. oryzae itself. Expression of ACE1 alone produced a polyenyl pyrone (magnaporthepyrone), which is regioselectively epoxidised and hydrolysed to give different diols, 6 and 7, in the two host organisms. Analysis of the three introns present in ACE1 determined that A. oryzae does not process intron 2 correctly, while M. oryzae processes all introns correctly in both appressoria and mycelia. Co-expression of ACE1 and RAP1 in A. oryzae produced an amide 8 which is similar to the PKS-NRPS derived backbone of the cytochalasans. Biological testing on rice leaves showed that neither the diols 6 and 7, nor amide 8 was responsible for the observed ACE1 mediated avirulence, however, gene cluster analysis suggests that the true avirulence signalling compound may be a tyrosine-derived cytochalasan compound.
ACE1 belongs to an infection-specific secondary metabolite gene cluster (Scheme 1),6 expressed only during appressorium-mediated penetration, and not at any other stage of the M. oryzae life cycle.7 The 12.4 kb ACE1 gene encodes an enzyme consisting of a fungal highly-reducing polyketide synthase (hrPKS)8 fused to a single module of a nonribosomal peptide synthetase (NRPS).4 Similar fungal synthetases have been investigated genetically and biochemically and are known to produce acyl tetramic acids or pyrrolidones including pretenellin A 1,9 prefusarin C 2 (ref. 10) and preaspyridone A 3 (ref. 11 and 12) as well as other potent bioactive compounds including cytochalasin K 4 (ref. 13 and 14) and the important cholesterol biosynthesis inhibitor lovastatin 5 (ref. 15) (Scheme 1). It appears likely, therefore, that ACE1 encodes a biosynthetic protein which makes a small molecule consisting of a polyketide fused to an amino acid. Because ACE1 is under very tight temporal and cell type-specific control in M. oryzae, it has not yet been possible to identify or purify the avirulence compound. Knowledge of the ACE1 metabolite structure would be highly useful for investigating avirulence signalling and developing novel resistant cultivars. The research of our group and that of others has focused on the heterologous expression of other fungal PKS-NRPS genes in Aspergillus oryzae from gene clusters of known function.12,16–19 Here we set out to use heterologous expression of genes from a cryptic pathway to determine its unknown chemical product, since knowledge of the ACE1 compound will enhance efforts to design new compounds of high potential utility in crop protection.
Scheme 1 Relationship between biosynthetic gene clusters and the structures of known compounds in various fungi (only partial gene clusters are shown). The ACE1 gene cluster (MGG_08381 – MGG_08391, MGG_12447)6 is shown and the ACE1/RAP1 gene pair which show significant homology to other genes known to produce bioactive compounds are highlighted. The top panel shows key chemical steps catalysed by the fungal PKS-NRPS TENS, which synthesises pretenellin A 1 in partnership with the trans-acting ER TENC. KS, β-ketoacyl synthase; AT, acyl transferase, DH, dehydratase; ER, enoyl reductase; ERO, defective enoyl reductase; CMeT, C-methyl transferase; KR, β-ketoacyl reductase; ACP, acyl carrier protein; C, condensation; A, adenylation; T, thiolation; DKC, Dieckmann cyclase; R, thiolester reductase. |
Fig. 1 Diagrammatic representation of fungal expression systems used in this work (not to scale). Introns indicated by grey bars. PKS-NRPS functional domains as described in the legend to Scheme 1. P = promoter, T = terminator. (A) cassettes used for expression of ACE1 and RAP1; (B) relative positions of ACE1 introns in cDNA and mRNA; (C) cassettes used for testing intron splicing, I = 1/2/3; (D) cassettes used for promoter testing and coexpression of RAP1. |
This plasmid was then used to transform A. oryzae M-2-3. Transformants were selected on minimal medium and 100 colonies were selected and grown further on both minimal medium and minimal medium plus maltose. Four transformants showed weak fluorescence and these were picked together with 10 non-fluorescent transformants and grown in liquid medium (see ESI†). The resultant organic extracts were examined by LCMS but no significant chemical differences were observed between extracts from any of the transformants and the untransformed A. oryzae M-2-3 strain (data not shown).
Each ACE1 intron was amplified with ca. 20 bp of flanking 5′ and 3′ exon sequences and cloned as in-frame fusions with eGFP downstream of either the native ACE1 promoter (PAce1), the strong constitutive promoter PgpdA from Aspergillus nidulans or the inducible promoter PamyB from A. oryzae. A total of nine plasmids were constructed (Fig. 1C). These constructs were transformed into either A. oryzae or M. oryzae, along with positive and negative controls, selected on the appropriate selection media and observed for fluorescence (Table 1 and Fig. 2). RT-PCR was performed for each intron construct using primers located either side of the introns and the products were sequenced (see ESI† for details).
Construct name | Expression host | Conditions | Fluorescent transformants |
---|---|---|---|
pPamyB·intron1·eGFP·argB | A. oryzae | M | 80% |
pPamyB·intron2·eGFP·argB | A. oryzae | M | 0% |
pPamyB·intron3·eGFP·argB | A. oryzae | M | 90% (Fig. 2A) |
Empty vector | A. oryzae | M | 0% |
eGFP | A. oryzae | M | 87% |
pPgpdA·intron1eGFP·hyg | M. oryzae | M, A | >70% |
pPgpdA·intron2eGFP·hyg | M. oryzae | M, A | >70% (Fig. 2B) |
pPgpdA·intron3eGFP·hyg | M. oryzae | M, A | >70% |
Empty vector | M. oryzae | M, A | 0% |
eGFP | M. oryzae | A | >70% |
pPACE1·intron1eGFP·hyg | M. oryzae | A | >60% (Fig. 2C) |
pPACE1·intron2eGFP·hyg | M. oryzae | A | >60% (Fig. 2D) |
pPACE1·intron3eGFP·hyg | M. oryzae | A | >60% (Fig. 2E) |
Empty vector | M. oryzae | M | 0% |
eGFP | M. oryzae | M, A | >60% |
The results demonstrated that introns 1 and 3 are spliced identically in A. oryzae and M. oryzae whereas incorrect splicing of intron 2 by A. oryzae results in a frameshift and premature stop codon (Fig. S1†).
The new compounds appeared to be identical for the two ACE1 expression constructs (RT = 33.2 min; uv max = 214, 272, 361 nm; m/z = 315 [M]Na+, 293 [M]H+, 275 [M − H2O]H+). Both compounds were purified and shown to be identical by NMR and HRMS (m/z 291.1238 Da [M − H]−, consistent with a molecular formula of C16H20O5).
1H NMR analysis (see ESI†) showed that the sample consisted of a 1:1 mixture of diastereomers each of which contained three methyl groups, a CHX, five contiguous olefinic protons, and two sharp doublets (J = 2.0 Hz) in the olefinic region. COSY correlations were used to determine that the olefinic protons were part of a distinct spin-system with an olefinic methyl group, connected by a vicinal coupling, visible only as a slight broadening of the terminal olefin doublet in the 1D 1H NMR. The COSY also showed that the CHX (which appeared as a complex multiplet) and one of the methyl groups (itself appearing as a pair of doublets) formed another distinct coupling system. An HMBC spectrum determined that the lone methyl group and olefins were all part of the same chain with a formula of C11H17O2, leaving C5H3O3 unaccounted for. The characteristic 13C NMR resonance of δ 163.6 ppm revealed that an ester group was present in the structure accounting for a further CO2, and leaving C4H3O unassigned.
The remaining two unassigned protons showed HMBC correlations to the ester carbonyl and two other quaternary carbons, the chemical shifts of which (δ 162.4 ppm and δ 162.6 ppm) were consistent with enolic carbons. Although these protons were coupled, they did not appear to be on adjacent carbons. An IR spectrum showed a characteristic absorbance at 1685 cm−1 suggesting that the ester group was present as part of a pyrone. Thus the structure was determined as the novel pyrone 6 to which we have assigned the name 12,13-dihydroxymagnaporthepyrone.
The polyketide origin of 6 was confirmed by feeding [1,2-13C2]-acetate and [methyl-13C]-methionine to A. oryzae pPamyB·ACE1[Δintrons1–3]eGFP·argB. Pyrone 6 isolated after the labelled acetate feed showed the presence of distinctive doublets in the 13C NMR indicative of the incorporation of intact 13C2 units (Fig. 4 and ESI†). Matching of the 2JCC values showed intact acetate incorporations consistent with the heptaketide origin of 6. The origin of both methyls from methionine was confirmed by the observation of a 9.4% incorporation of label at C-15 and C-16.
Fig. 4 Structures and biosynthetic origin of the pyrones obtained by expressing ACE1 in A. oryzae M-2-3 and M. oryzae. |
Thus 6 consists entirely of a polyketide and does not contain the expected amino acid component. To rule out NRPS perturbance by the C-terminal eGFP fusion we modified the expression system to remove the 3′-eGFP gene to create pPamyB·ACE1[Δintrons1–3]·argB and transformed A. oryzae M-2-3. Selected transformants were grown in liquid medium and secondary metabolites detected using LCMS; 12,13-dihydroxymagnaporthepyrone 6 was the only new compound detected.
To overcome the possibility that A. oryzae lacks the specific amino acid required by the ACE1 NRPS in M. oryzae we constructed plasmids for homologous expression of ACE1 in M. oryzae Guy11 using the PgpdA constitutive promoter and hygromycin selection. We thus constructed an E. coli shuttle vector based on the pUC18 vector with a Gateway cassette downstream from PgpdA carrying a hygromycin resistance marker. The ACE1 eGFP construct was transferred into this vector by Gateway LR recombination in vitro to give pPgpdA·ACE1eGFP·hyg which was used to transform M. oryzae Guy11. The plasmid pPgpdA·eGFP·hyg was used as a control. Hygromyin resistant transformants were selected for each construct and grown on complete medium plates and examined for fluorescence. For pPgpdA·eGFP·hyg, 18/20 transformants were fluorescent while only 5/100 of the pPgpdA·ACE1eGFP·hyg transformants were fluorescent.
The five M. oryzae transformants expressing ACE1eGFP were then grown in complete liquid medium for 7 days then extracted and examined by LCMS as described in the ESI.† All five transformants produced a new compound which was not produced by the controls. This compound had a different retention time to 12,13 dihydroxymagnaporthepyrone 6 produced in A. oryzae (RT = 31.6 min, Fig. 3B) and a different UV spectrum (uvmax = 210 nm, 248 nm, 340 nm) although the low resolution ESI MS spectrum was identical to 6, as was the molecular formula predicted from the HRMS. Thus it appeared that this compound was isomeric to 12,13-dihydroxymagnaporthepyrone 6. 1H NMR showed it to consist of a single stereoisomer, but compared to 6 it had only three contiguous olefinic CH resonances, one of which was coupled to a CHO methine which was in turn coupled to another CHO methine. The methyl doublet of 6 was replaced by a methyl doublet coupling to an olefinic CH. These data, supported by 13C and 2D spectra showed this compound to be 10,11-dihydroxymagnaporthepyrone 7.
For A. oryzae we created a plasmid carrying both the ACE1 and RAP1 genes in which ACE1 lacking all introns was expressed under the control of PamyB, and RAP1 was expressed from the strong constitutive alcohol dehydrogenase promoter Padh from A. oryzae (pPamyB·ACE1[Δintrons 1–3]·argB·Padh·RAP1, Fig. 1A). A. oryzae M-2-3 protoplasts were transformed with pPamyB·ACE1[Δintrons 1–3]·argB·Padh·RAP1 and transformants were selected on minimal medium.
Twenty seven of the resulting A. oryzae transformants were selected and transferred to new agar plates. Ten of these transformants were grown in liquid medium containing starch for seven days, then extracted in the usual way. LCMS analysis of the extracts showed that of the A. oryzae transformants, five produced a new peak compared to the controls in the neutral organic extract, and in lesser amounts, in the acidic and mycelial extracts. The new compound had different retention times to compounds 6 and 7 (RT = 22.5 min, uv max = 270, 280 nm, m/z 481.4 [M − H]− (Fig. 5)).
A single A. oryzae transformant was grown at large scale (800 mL) to yield 65.5 mg of crude extract. The new compound was purified by HPLC and analysed by NMR and HRMS (m/z 482.3265 Da [M]H+, consistent with a molecular formula of C30H43NO4). 1H NMR analysis showed that the compound contained four methyl groups, six olefinic protons and an aromatic group indicative of tyrosine. 13C NMR, COSY, HMBC and HSQC determined the structure as 8 (Fig. 6, ESI†).
Fig. 7 Biological testing of 12,13-dihydroxymagnaporthepyrone 6 on both resistant (IR64 and C101LAC) and susceptible (Maratelli) rice cultivars. Testing of 7 gave similar results. 8 Had no effect. |
Further investigation using RT-PCR and subsequent sequence analysis clearly showed that ACE1 intron-2 is not spliced correctly in A. oryzae, resulting in an mRNA with a frame-shift and premature stop codon. Intron-2 is located centrally in the PKS portion of ACE1, upstream of the sequence encoding the key acyl carrier protein (ACP) domain. Any truncated ACE1 protein resulting from the translation of this incorrectly spliced mRNA would not be capable of producing a polyketide product.
Removal of the ACE1 introns resulted in transcripts which were translated correctly in A. oryzae. Indeed more than 75% of transformants produced full-length ACE1-eGFP protein fusions as evidenced by strong fluorescence of transformed hyphae, and the production of novel metabolites. Expression of both full-length and NRPS-deleted ACE1 intronless constructs resulted in the production of the same compound, 12,13-dihydroxymagnaporthepyrone 6 which is exclusively of polyketide origin as shown by isotopic labelling. Expression of the same clones in M. oryzae resulted in the production of the closely related 10,11-dihydroxymagnaporthepyrone 7. Similar compounds were observed by Vederas and coworkers when the lovasatin nonaketide synthase (LNKS) was expressed without its cognate ER encoded by lovC.15
These results can be interpreted if the PKS portion of ACE1 produces the polyunsaturated pyrone 9 (which we name magnaporthepyrone, Scheme 2). This compound is then epoxidised by different monooxygenases in either A. oryzae or M. oryzae and spontaneous hydrolysis of the epoxides results in the observed diols 6 and 7. Support for this hypothesis comes from the expression of tenS in A. oryzae where we observed the presence of prototenellin C 10 (ref. 22) which displays exactly the same chemical motif as 12,13-dihydroxymagnaporthepyrone 6 – thus A. oryzae must possess a monooxidase selective for the terminal methylbutenyl motif of such precursors. Likewise, M. oryzae is known to produce diol compounds such as pyriculol 11 (ref. 23) which must also derive from epoxidation of a polyunsaturated polyketide followed by hydrolysis. Therefore, M. oryzae must have a monooxygenase with selectivity for epoxidation of mid-chain polyenes.
Scheme 2 Deduced biosynthesis of ACE1 products in the presence and absence of RAP1. X = ACE1 ACP domain. |
The biosynthesis of 6 and 7 do not require an enoyl reductase step, consistent with the observation that the ACE1 ER domain is predicted to be inactive. Coexpression of ACE1 with RAP1, which encodes a trans-acting ER, changed the biosynthetic pathway. In the heterologous host A. oryzae, where we constructed a single vector to carry both genes, the amide 8 was produced and compounds derived from 9 were no longer observed.
Amide 8 is closely related to pyrone 9 and its biosynthesis can be rationalised by the intervention of the RAP1 ER at the pentaketide stage of biosynthesis (Scheme 2). Pentaketide intermediate 12 is recognised by the ACE1 KS domain and extended to 13, but in the absence of prior enoyl reduction the next required CMeT step is prevented. Continued chain extension by an unselective KS then gives the tricarbonyl 14 which can spontaneously offload itself from the synthase as the pyrone 9. Alternatively, in the presence of the RAP1 ER, pentaketide 12 is reduced giving a substrate which can be extended by the KS and then correctly methylated to give 15 before following three more extension and processing cycles to give the fully extended nonaketide 16. Thus the ACE1 CMeT domain appears to display an additional level of substrate selectivity in blocking correct chain processing in the absence of 10,11 reduction. Crucially, however, the KS is unselective and continues to extend the chain to the tricarbonyl 14 which can be off-loaded as the pyrone 9. Thus the lack of selectivity by the KS effectively removes incorrectly processed chains from the PKS, preventing its blockage.
In many other investigated PKS-NRPS systems the NRPS then attaches the amino acid to the fully extended polyketide and either reductively cleaves it,24 or performs a non-reductive Dieckmann cyclisation reaction to release an acyl tetramic acid.9 In this case it appears that tyrosine is used as the amino acid and attached to the polyketide in the same way as observed in the cases of pretenellin A 1 and preaspyridone 3. After amide formation, however, the fully elaborated product is probably reductively released. This could be via 2 electron reduction to form an aldehyde (as in the classic case of fungal lysine biosynthesis25) which could be further reduced to a primary alcohol by an adventitious A. oryzae enzyme, or via double (i.e. 4 electron) reduction catalysed by the NRPS R-domain itself as observed during the biosynthesis of myxochelin.26
Consistent with this hypothesis is the observation that the ACE1 terminal-domain has an intact conserved NADPH binding motif (GXSXXG) and the catalytic triad Ser-Tyr-Lys more closely resembles ‘reducing’ R-domains than DKC cyclising domains (see ESI†).24 Qiao et al. predicted a reductive release for ccsA, the PKS-NRPS gene involved in cytochalasin K 4 biosynthesis in Aspergillus clavatus,27 based on the same rationale, and the co-expression of ccsA with ccsC by Oikawa and co-workers led to the production of the amide 17,14 also at the alcohol oxidation state.
The Oikawa compound 17 differs from 8 in being an octaketide rather than a nonaketide and being evidently constructed from phenylalanine rather than tyrosine (Scheme 3). A further difference is that 8 is unsaturated between carbons 8 and 9 of the polyketide, whereas 17 is fully saturated at the corresponding position. Genetic experiments have already shown that ccsA and ccsC are involved in cytochalasin K 4 biosynthesis, most likely via ketocytochalasin 18.27 Although 17 is unlikely to be a direct precursor of 4, it is probably a shunt metabolite from the corresponding aldehyde. Vederas, Tang and coworkers have also recently elucidated that ccsB from the cytochalasin K biosynthetic pathway is an FAD-dependent oxygenase which catalyses (among other reactions) the creation of the distinctive carbonate moiety of 4.28
The structural similarities between the ccsA/C product 17 and the ACE1/RAP1 product 8 strongly suggest that 8 is also a shunt metabolite of a cytochalasan-like biosynthetic pathway.13 This is supported by the observation that 8 does not show the expected biological activity as an avirulence signalling compound.
To our knowledge, only two other gene clusters have been correlated with cytochalasan biosynthesis in fungi. These are the gene clusters for chaetoglobosin A 19 biosynthesis in Penicillium expansum and Chaetomium globosum.29,30 Chaetoglobosin A 19 is biosynthesised from tryptophan and a nonaketide likely to be almost identical to that of amide 8, except for the unsaturation between carbons 4 and 5. However, several chaetoglobosins are also known which are fully saturated at this position (e.g. chaetoglobosin C 20)31 indicating that amide 8 could include the correct polyketide precursor. Analysis of the tailoring genes in the ACE1 cluster show that they have more similarity to the chaetoglobosin biosynthetic genes than those in the ccs cytochalasin gene cluster (Table S2†).32
The nonaketide nature of the amide 8 and the fact that tyrosine is incorporated, distinguishes the RAP1/ACE1 product from the ccsA/ccsC compound and mean that the likely ultimate product of the ACE1 pathway will be a member of the tyrosine-derived cytochalasans, an example of which is phomopsichalasin (also known as diaporthichalasin) 21.33,34 In addition, the differences between the ACE1 and ccs gene clusters, such as the lack of a ccsB homolog, plus additional oxidoreductases similar to those found in the chaetoglobosin A 19 gene cluster, indicate that the tailoring of the polyketide backbone is likely to be more similar to 19, which could arise from a polyketide identical or very similar to that of 8.
Thus we have shown that heterologous expression is an effective tool to probe cryptic gene clusters in fungi, which enables important metabolic relationships to be unravelled. It is now clear that avirulence signalling between the rice blast fungus Magnaporthe oryzae and rice is mediated by a natural product bearing structural similarities to the cytochalasans. Current and future work in our group is focused on exploring the enzymatic activities of the proteins encoded by other genes from the ACE1 biosynthetic cluster to identify the bioactive compound responsible for avirulence signalling.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c4sc03707c |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2015 |