Original Research ARTICLE
Perturbations in small molecule synthesis uncovers an iron-responsive secondary metabolite network in Aspergillus fumigatus
- 1Department of Medical Microbiology and Immunology, University of Wisconsin-Madison, Madison, WI, USA
- 2Division of Molecular Biology/Biocenter, Innsbruck Medical University, Innsbruck, Austria
- 3Boyce Thompson Institute and Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, USA
- 4The J. Craig Venter Institute, Rockville, MD, USA
- 5Department of Bacteriology, University of Wisconsin-Madison, Madison, WI, USA
Iron plays a critical role in survival and virulence of the opportunistic pathogen Aspergillus fumigatus. Two transcription factors, the GATA-factor SreA and the bZip-factor HapX oppositely monitor iron homeostasis with HapX activating iron acquisition pathways (e.g., siderophores) and shutting down iron consumptive pathways (and SreA) during iron starvation conditions whereas SreA negatively regulates HapX and corresponding pathways during iron sufficiency. Recently the non-ribosomal peptide, hexadehydroastechrome (HAS; a tryptophan-derived iron (III)-complex), has been found important in A. fumigatus virulence. We found that HAS overproduction caused an iron starvation phenotype, from alteration of siderophore pools to regulation of iron homeostasis gene expression including sreA. Moreover, we uncovered an iron dependent secondary metabolism network where both SreA and HapX oppositely regulate multiple other secondary metabolites including HAS. This circuitry links iron-acquisition and consumption pathways with secondary metabolism—thus placing HAS as part of a metabolic feedback circuitry designed to balance iron pools in the fungus and presenting iron availability as one environmental trigger of secondary metabolism.
The ubiquitous soil-dwelling filamentous fungus Aspergillus fumigatus can cause the life-threating disease invasive aspergillosis in immunocompromised individuals. It is regarded as the most common airborne fungal pathogen of humans (Tekaia and Latge, 2005). A. fumigatus has the ability to produce hundreds of bioactive small molecules, so called secondary metabolites (Frisvad et al., 2009). With few exceptions, genes required for their production are generally clustered in fungal genomes, each gene in a specific cluster being subject to common regulatory patterns often triggered by largely unknown environmental cues (Hoffmeister and Keller, 2007; Brakhage, 2013). Some of these secondary metabolites have been shown to contribute to the fungus' virulence (Jahn et al., 2002; Bok et al., 2006; Heinekamp et al., 2012; Berthier et al., 2013; Yin et al., 2013).
One compound recently associated with enhanced virulence is hexadehydroastechrome (HAS), an iron-binding, non-ribosomal peptide-derived molecule (Yin et al., 2013). For the majority of organisms, including fungi, iron is an essential cofactor in several enzymatic reactions and serves as a catalyst in electron transport. However, over-abundance of iron can lead to formation of reactive oxygen species that are highly damaging to cellular components (Halliwell and Gutteridge, 1984). As iron availability is essential for survival of microorganism, the mammalian defense machinery includes iron-withholding systems in order to fight microbial infections (Weinberg, 1999; Weiss, 2002; Fluckinger et al., 2004; Ganz, 2009; Leal et al., 2013). Therefore, controlling iron access during infection is a major determinant of successful microbial infection of the mammalian host. Iron deficiency is also known to serve as a regulatory cue for other virulence determinants in microorganisms (Litwin and Calderwood, 1993; Weinberg, 1999; Oglesby-Sherrouse et al., 2014) and has been shown to induce the production of the virulence factor ribotoxin AspF1 in A. fumigatus (Schrettl et al., 2010a).
A. fumigatus harbors two high-affinity iron uptake systems, siderophore-assisted iron uptake and reductive iron uptake (Schrettl et al., 2004a; Schrettl and Haas, 2011). Siderophores are non-ribosomal peptide-derived ferric iron chelators (Haas, 2003). Aspergillus spp. produces the extracellular siderophores fusarinine C (FC) and triacetylfusarinine C (TAFC) and the intracellular iron storage siderophores ferricrocin (FC) and hydroxyl ferricrocin (HFC) (Oberegger et al., 2001; Schrettl et al., 2007; Wallner et al., 2009). The first devoted step in extra- and intracellular siderophore biosynthesis is catalyzed by the L-ornithine N5-oxygenase SidA (Schrettl et al., 2004a, 2007). The non-proteinogenic amino acid produced by SidA then bifurcates into two independent non-ribosomal peptide biosynthesis pathways devoted to produce both groups of siderophores (Schrettl et al., 2004a). In contrast to the reductive iron uptake system, siderophore production was shown to be an essential virulence determinant of A. fumigatus and other fungal species (Eisendle et al., 2003; Schrettl et al., 2004a, 2007; Oide et al., 2006; Greenshields et al., 2007; Gauthier et al., 2010; Hwang et al., 2012; Leal et al., 2013).
As iron exhibits a Janus-faced role in the physiology of an organism, cellular iron homeostasis is tightly regulated in order to ensure iron availability on one side and prevent toxic iron excess on the other side. In the model organism A. nidulans, efficient control of iron homeostasis is executed by the two transcription factors, HapX and SreA that are interconnected in a negative feedback loop (Haas et al., 1999; Hortschansky et al., 2007). HapX is a bZIP protein that is thought to directly sense iron deficiency and execute repression of iron-dependent pathways and sreA expression by protein-protein interaction with the CCAAT-binding complex under iron-limiting conditions (Hortschansky et al., 2007). SreA is a DNA-binding GATA-type transcription factor, which represses hapX expression and iron acquisition pathways under iron-replete conditions (Haas et al., 1999). Both factors are conserved in A. fumigatus and other fungi executing similar roles in iron regulation (Voisard et al., 1993; An et al., 1997; Zhou et al., 1998; Hwang et al., 2008, 2012; Schrettl et al., 2008, 2010a; Gauthier et al., 2010; Wiemann et al., 2012; Leal et al., 2013). In accordance with the requirement of siderophores for establishment of full virulence, deletion of hapX was shown to attenuate virulence in A. fumigatus and Fusarium oxysporum (Schrettl et al., 2010a; Leal et al., 2013; Lopez-Berges et al., 2013).
Considering the importance of iron pools in the fungus, coupled with the iron chelating properties of the virulence factor HAS, we set out to characterize the role of HAS in iron-dependent regulation in A. fumigatus. We show that on the one hand HAS and, unexpectedly several other secondary metabolites, are subject to iron-dependent regulation by SreA and HapX and on the other hand, over-production of HAS changes HapX/SreA-target gene expression resulting in perturbations in iron-acquisition and consumption pathways. Taken together, this work implicates iron gradients as important environmental cues regulating secondary metabolite synthesis.
Materials and Methods
Strains and Culture Conditions
A. fumigatus strains used in this study are listed in Table S1. Strains were maintained as glycerol stocks and activated on solid glucose minimal medium (GMM) at 37°C with appropriate supplements (Shimizu and Keller, 2001). For pyrG auxotrophs, the growth medium was supplemented with 5 mM uridine and uracil. Conidia were harvested in 0.01% Tween 80 and enumerated using a hemacytometer. For RNA-seq analysis strains Af293, TWY32.1, and TWY24.121 were inoculated into 50 mL of liquid GMM at 5 x 106 conidia/mL in duplicate and grown at 25°C and 250 rpm for 96 h in ambient light conditions. The mycelium was harvested and lyophilized before RNA extraction. For biomass, siderophore and iron-dependent gene expression analysis, strains CEA17, TWY37.2, TJW109.3, TWY25.5, and TWY28.3 were grown in 50 mL Aspergillus minimal media (AMM) according to (Pontecorvo et al., 1953) containing 20 mM glutamine and iron concentrations as indicated in the text according to Schrettl et al. (2008). Strains were grown in triplicates (duplicates for gene expression analysis) for 24 h at 37°C, 250 rpm and ambient light conditions with an initial spore concentration of 5 × 106 conidia/mL. For analysis of secondary metabolites and cluster gene expression, strains ATCC46645, Δ hapX, and Δ sreA were inoculated into 50 mL of GMM containing different iron concentrations as indicated in the text at 5 × 106 conidia/mL. Strains were grown in triplicates (duplicates for gene expression analysis) for 5 days (3 days for gene expression analysis) at 25°C and 250 rpm at ambient light conditions. For growth assays strains TJW55.2, TWY32.1, TWY38.6, TWY24.121, TWY35.1, and TW36.1 indicated amount of conidia were inoculated in 2 μL on solidified (Noble Agar, Difco™, BD, USA) GMM containing indicated iron concentrations and FEC, respectively, and incubated for 3 days at 37°C in the dark.
DNA Isolation, Genetic Manipulations and Southern Blot Analysis
For DNA isolation, A. fumigatus strains were grown for 24 h at 37°C in steady state liquid GMM, supplemented with 1 mM FeSO4. DNA isolation was performed as described by Green and Sambrook (2012). The A. fumigatus mutant strains were constructed using a double-joint fusion PCR approach (Szewczyk et al., 2006) and all primers used in the study are listed in Table S2. Briefly, approximately 1 kb fragments flanking the targeted deletion region were amplified by PCR from A. fumigatus strain Af293 or CEA17 genomic DNA and the A. parasiticus pyrG marker gene was amplified from the plasmid pJW24 (Calvo et al., 2004). The three fragments were subjected to fusion PCR to generate deletion cassettes. The A. parasiticus pyrG marker gene was amplified from the plasmid pJW24 using the primer pair pyrG_prom_F/pyrG_term_R (Table S2). The primers “gene”-5R and “gene”-3F contain complement sequences to the primers pyrG_prom_F and pyrG_term_R at their 5'-region, respectively (Table S2). The fusion construct was created by PCR containing 5′ and 3′ gene flanks and the pyrG gene fragment functioning as templates and primers simultaneously. The final PCR fusion product was amplified using primer pairs “gene”-5F/“gene”-3R and the previously PCR-generated fusion construct as template. Transformation was performed as described by Palmer et al. (2008). For selection of Δ sidA transformants, 1 mM FeSO4 was added to the selection media. For multiplex diagnostic PCR, primer pair “gene”-F/“gene”-R were used to identify transformants that lost the respective gene locus and primer pair gpd_int-F/gpd_int-R as internal control (Table S2). Single integration of the transformation construct was confirmed by Southern analysis as described by Green and Sambrook (2012) using P32-labeled probes created by amplification of the respective deletion construct using primer pair “gene”-5F/“gene”-3R (Table S2 and Figure S3).
RNA Extraction and Northern Analysis
Mycelia were harvested by filtration through Miracloth (Calbiochem). Total RNA was extracted with TRIzol reagent (Invitrogen) from freeze-dried mycelia, following the manufacturer's protocol. Northern analysis was performed as described by Green and Sambrook (2012). Probes for northern analysis were constructed at regions internal to the gene of interest using primers listed in Table S2 (“gene”-F/“gene”-R) and labeled with dCTP αP32.
RNA Sample Preparation and Illumina Sequencing (RNASeq)
To identify transcriptionally active genes, extracted RNA samples were subject to DNase treatment using the RNeasy kit (Qiagen), and sequencing libraries were generated using the ScriptSeq kit v2 (Epicenter) following manufacturer's directions. All libraries were sequenced (7 sample per lane) using the Illumina HiSeq2000 instrument (www.illumina.com) on a 2 × 100 bp paired-end run. Transcripts were assembled and expression levels were estimated using normalized reads per kilobase per million mapped reads (RPKM) (Mortazavi et al., 2008) as calculated by the TopHat, Bowtie, and Cufflinks packages (Trapnell et al., 2009, 2010, 2012; Langmead et al., 2010) and CLC Genomics Workbench (CLC Bio). After trimming for quality (>Q30), reads were mapped at 90% length, 90% identity and reads that mapped to more than one location in the genome were excluded. A minimum of 4 reads in each condition were required to analyze gene expression, over 88% of genes had over 4 reads in every condition. Differentially expressed genes were determined in duplicate data by identifying genes that were at least two-fold difference using a T-test cutoff of p < 0.05 (Bullard et al., 2010) and by using SAM with a false discovery rate of 0 as implemented in MeV (Saeed et al., 2006). Replicates of each condition were highly reproducible (R2 > 0.99) (Dataset S1). For gene enrichment analysis, genes differentially regulated were analyzed for enrichment using the Munich Information Center for Protein Sequences (MIPS) Functional Catalog (FunCat) (Ruepp et al., 2004) at http://www.helmholtz-muenchen.de/en/ibis/resourcesservices/index.html (Dataset S2).
Siderophore and Secondary Metabolite Analysis
Analysis of siderophores was carried out by reversed phase high performance liquid chromatography (HPLC) as described previously (Oberegger et al., 2001). Moreover, to quantify extracellular or intracellular siderophores, culture supernatants or cellular extracts were saturated with FeSO4 and siderophores were extracted with 0.2 volumes of phenol. The phenol phase was separated and subsequent to addition of five volumes of diethylether and 1 volume of water, the siderophore concentration of the aqueous phase was measured photometrically using a molar extinction factor of 2996 M−1 cm−1 at 440 nm.
For analysis of pseurotin A, fumitremorgin C and fumagillin supernatant of fungal cultures was lyophilized and re-dissolved in 5 mL double destilled (dd) H2O 1% (v/v) formic acid (FA). From the water crude, 800 μL were partitioned with 800 μL ethyl acetate for three times. The combined organic layers of each sample were evaporated in vacuo, and re-dissolved in 500 μL of 20% (v/v) acetonitrile (ACN) 1% FA (v/v) and 50 μL were examined by HPLC photo diode array (PDA) analysis. The samples were separated on a ZORBAX Eclipse XDB-C18 column (Agilent, 4.6 mm by 150 mm with a 5 μm particle size) by using a binary gradient of 1% (v/v) FA as solvent A and 1% FA in ACN as solvent B using a Flexar Binary Liquid Chromatography (LC) Pump (PerkinElmer) coupled to a Flexar LC Autosampler (Perkin Elmer) and a Flexar PDA Plus Detector (PerkinElmer). The binary gradient started with an isocratic step at 80% A for 2 min followed by a linear gradient to 40% A in 10 min and an additional linear gradient to 100% B in 0.5 min. After each run the column was washed for 5 min using 100% B and was equilibrated for 4 min using 80% A. The flow rate was set to 1.5 mL/min. Identification and relative quantification of secondary metabolites was performed using Chromera Manager (PerkinElmer) by comparison to standards for fumagillin, pseurotin A and fumitremorgin C (Cayman Chemicals, MI, USA) (Dataset S3). Areas under the curve were normalized to dry weights of each sample. For terezine D analysis liquid fungal cultures including fungal tissue and media were frozen using a dry ice acetone bath, and lyophilized. The lyophilized residues were extracted with 75 mL of 10% methanol in ethyl acetate for 3.5 h with vigorously stirring. Extracts were filtered over cotton, evaporated to dryness, and stored in 8 mL vials. Preparation for HPLC-mass spectrometry (MS) Analysis: Crude extracts were suspended in 0.1 mL of methanol and centrifuged to remove insoluble materials, and the supernatant was subjected to HPLC-mass spectrometry analysis using an Agilent 1100 series HPLC connected to a Quattro II low resolution mass-spectrometer (Micromass/Waters) operated in electrospray positive ionization (ESI+). Data acquisition and processing for the HPLC-MS was controlled by Waters MassLynx software. An Agilent Zorbax Eclipse XDB-C18 column (4.6 × 250 mm, 5 μm particle diameter) was used with 0.1% acetic acid in 50/50 methanol/acetonitrile (organic phase) and 0.1 % acetic acid in water (aqueous phase) as solvents at a flow rate of 1.0 mL/min. A solvent gradient scheme was used, starting at 5% organic for 3 min, followed by a linear increase to 100% organic over 35 min, holding at 100% organic for 15 min, then decreasing back to 5% organic for 1 min and holding at 5% organic for the final 6 min, a total of 1 h. Relative abundance of terezine D was determined by integrating their extracted ion chromatograms using MassLynx software (Dataset S3).
Statistical analysis was performed by using GraphPad Prism software using analysis of variance (ANOVA) followed by Tukey's test to show significant differences.
HAS does not Contribute to Siderophore-Assisted Iron Uptake
Since HAS was shown to chelate ferric iron, similar to the intra- and extra-cellular siderophores produced by A. fumigatus, we investigated whether HAS could restore growth defects in an A. fumigatus sidA deletion mutant deficient of siderophore production. Deletion mutants of sidA in A. fumigatus were previously shown to exhibit strong growth defects under iron deficient conditions and to be avirulent in animal models of invasive aspergillosis (Schrettl et al., 2004a, 2007; Hissen et al., 2005; Knox et al., 2014). Single (ΔsidA) and double (OE::hasA/ΔsidA) mutants of sidA with or without an overexpression hasA allele were created as described in Material and Methods (Figure S1). HasA is the has cluster transcription factor positively activating other members of the cluster and HAS biosynthesis (Yin et al., 2013).
Confirming previous results (Schrettl et al., 2004a), sidA deletion mutants showed significant growth defects under iron depleted conditions compared to the wild type (WT), the hasA deletion mutant (ΔhasA), the overexpression hasA mutant (OE::hasA), as well as a mutant overexpressing hasA in a hasD deletion background (OE::hasA/ΔhasD) where hasD encodes the non-ribosomal peptide synthetase required for HAS biosynthesis (Figure 1). We found that OE::hasA in the ΔsidA background did not rescue growth (Figure 1), thereby ruling out the hypothesis that HAS could function as an extra-cellular siderophore assisting in iron uptake. Growth of ΔsidA and OE::hasA/ΔsidA could be restored to wild-type like levels when 2 μM of the Schizosaccharomyces pombe siderophore ferrichrome (Schrettl et al., 2004b) was supplemented to iron deplete conditions. Under extreme iron excess (10 mM), growth rates of the WT and the OE::hasA strains were decreased to similar levels (Figure 3), indicating that HAS, similar to FC and HFC, does not participate in iron detoxification in contrast to the vacuolar iron transporter CccA (Gsaller et al., 2012) (Figure 1).
Figure 1. Growth phenotypes of the wild type, ΔhasA, OE::hasA, OE::hasA/ΔhasD, ΔsidA, and OE::hasA/ΔsidA strains on different iron conditions using solidified media. Indicated amount of spores were spotted on solidified GMM with various iron concentrations and 2 μM ferrichrome (FEC) as indicated. Strains were grown for 3 days at 37°C in the dark.
HAS and Astechrome Bind Iron in vivo
In order to confirm the ability of HAS and its precursor astechrome to bind iron as reported previously (Arai et al., 1981; Yin et al., 2013), we cultivated the WT, OE::hasA, OE::hasA/ΔhasD, and OE::hasA/ΔhasG strains under increasing iron conditions. HasG is a FAD binding protein required to convert astechrome to HAS and both compounds chelate iron (Arai et al., 1981; Yin et al., 2013). OE::hasA and OE::hasA/ΔhasG exhibited a reddish mycelial coloration particularly during iron excess in contrast to the wild type and OE::hasA/ΔhasD (Figure 2).
Figure 2. Growth phenotypes of wild type, OE::hasA, OE::hasA/ΔhasD, and OE::hasA/ΔhasG strains in different liquid iron conditions. Indicated strains were incubated in AMM supplemented with indicated iron concentrations for 1 day at 37°C at 200 rpm.
Overproduction of HAS Causes Growth Defects and Increased Siderophore Production Under Iron-Replete Conditions
Since HAS could be excluded to function as an extracellular siderophore, we investigated if it could function in iron detoxification under iron-replete and iron excess conditions. Therefore, we assessed biomass of the WT, ΔhasA, OE::hasA, and OE::hasA/ΔhasD strains. During iron-replete (30 μM), iron deficient (1 μM) and no iron (0 μM) conditions OE::hasA showed decreased biomass production compared to the other strains tested (Figure 3). In case of OE::hasA this growth defect could be rescued by addition of 1 mM iron to the media (Figure 3). Taken together these data suggest that HAS causes growth deficiency due to its iron-chelating ability.
Figure 3. Mycelial dry weights of wild type, ΔhasA, OE::hasA and OE::hasA/ΔhasG under different iron conditions. Indicated strains were grown in liquid 50 mL AMM for 1 day at 37°C at 250 rpm. Harvested mycelia was lyophilized before weighing. Error bars indicate standard deviations for triplicates of each strain. Asterisks indicate significance as calculated by Tukey's test following ANOVA: **p < 0.05.
In line with the observed biomass reduction, production of intra- and extra-cellular siderophores was increased under iron-replete conditions (30 μM) in strains overproducing HAS (OE::hasA) compared to the WT or OE::hasA/ΔhasD (Figure 4). At iron-replete conditions (30 μM iron), OE::hasA produces more extracellular and intracellular siderophores than WT or the other has mutants thus supporting the hypothesis that HAS causes the fungus to experience iron-deficiency in this condition. However, under iron starvation (0 μM) no increase in siderophore production could be observed in OE::hasA compared to the WT, most likely due to the fact that siderophore production is already fully induced under those conditions in all strains. In iron excess conditions, siderophore production is similar across all strains (Figure 4).
Figure 4. Production of extra and intracellular siderophores in the wild type, ΔhasA, OE::hasA, and OE::hasA/ΔhasD strains in different iron conditions. Indicated strains were grown in liquid AMM supplemented with indicated iron concentrations at 37°C for 1 day at 200 rpm. Siderophore concentration was measured as described in Material and Methods and normalized to the amount produced by the WT and plotted in log scale. Error bars indicate standard deviations for triplicates of each strain. Asterisks indicate significance as calculated by Tukey's test following ANOVA: *p < 0.0005.
HAS and Astechrome Overproduction Increases Expression of Genes Involved in Iron Acquisition and Represses Genes Involved in Iron-Dependent Pathways
In A. fumigatus, iron starvation induces the expression of genes involved in siderophore biosynthesis (e.g., sidA) and reductive iron uptake (e.g., ftrA). Moreover, iron starvation represses genes involved in iron dependent pathways such as cycA (Schrettl et al., 2010a). During iron-replete conditions (30 μM), OE::hasA and OE::hasA/ΔhasG increase expression of sidA and ftrA but repress cycA expression compared to the WT, ΔhasA and OE::hasA/ΔhasD, respectively (Figure 5A). These data indicate that HAS or astechrome accumulation causes iron starvation, which is in agreement with the reduced biomass formation and increased siderophore production under this condition (Figures 3, 4). Since expression of genes involved in iron acquisition and iron-utilization are repressed by SreA and HapX, respectively, in A. fumigatus (Schrettl et al., 2008, 2010a), we tested expression of sreA in iron deplete conditions and conditions exposed to iron sufficiency for a short time. As previously described (Schrettl et al., 2008), shifting to iron-replete conditions caused induction of sreA expression compared to iron-deplete conditions in the WT. This induction is less pronounced in OE::hasA and OE::hasA/ΔhasG compared to the other strains tested (Figure 5B). These data are in line with the growth and siderophore production data indicating that HAS overexpression decreases iron bioavailability.
Figure 5. Expression analysis for select HAS target genes for the wild type, ΔhasA, OE::hasA, OE::hasA/ΔhasD, and OE::hasA/ΔhasG under different iron conditions. (A) Strains were grown for 24 h as described in Material and Methods under the indicated iron conditions. (B) Strains were grown for 24 h in −Fe media as described in Material and Methods. One set of cultures was supplemented with 10 μM FeSO4 for 30 min before harvesting. Northern blot hybridization of P32 labeled DNA probes to indicated target genes. Ethidium bromide-stained rRNA visualized as control.
HAS and Other Secondary Metabolites are Controlled by Ambient Iron Conditions in a HapX- and SreA-Dependent Manner
As our cumulative data demonstrate a role of HAS in iron homeostasis of A. fumigatus we were curious if the has gene cluster itself is subject to HapX- or SreA-mediated iron-dependent expression. Therefore, we grew the WT, ΔhapX, and ΔsreA under iron deplete (0 μM), replete (30 μM) and excess (1 mM) conditions and assessed HAS production and gene expression by northern blot analysis. Expression of select has genes was repressed by iron depleted conditions in a HapX-independent manner (Figure 6A). Gene expression of hasB and concomitant terezine D production, the first stable intermediate of the HAS pathway, was induced by medium-high iron concentrations and SreA-deficiency, which increases the cellular iron accumulation (Figures 6A,B).
Figure 6. Analysis of HAS biosynthesis in the wild type, ΔhapX and ΔsreA mutants under different iron conditions. (A) Strains were grown for 3 days at 25°C at 250 rpm in GMM supplemented with indicated iron concentrations. Northern blot hybridization of P32 labeled DNA probes to indicated target genes. Ethidium bromide-stained rRNA visualized as control. (B) Strains were grown at 25°C for 5 days at 250 rpm in GMM supplemented with indicated iron concentrations. Terezine D was measured as described in Material and Methods and plotted as percent of the amount produced in the wild type at 0, 0.03 and 1 mM Fe, respectively, in log scale. Error bars indicate standard deviations for triplicates of each strain. Asterisks indicate significance as calculated by Tukey's test following ANOVA: **p < 0.0005.
Since we observed an iron-dependent expression of the has cluster and HAS production, we investigated expression of select genes from the gene clusters responsible for fumitremorgin C, pseurotin A and fumagillin production, as these secondary metabolites are known to be produced under the same wild-type conditions and therefore convenient to measure (Wiemann et al., 2013). Similar to the observed regulation for the has cluster, select genes from the pso and fma supercluster (Wiemann et al., 2013) were co-activated followed by increased product formation with increasing iron concentration in the WT (Figures 7, 8).
Figure 7. Expression analysis for select secondary metabolite cluster genes for the wild type, ΔhapX and ΔsreA mutants under different iron conditions. Strains were grown for 3 days at 25°C for 3 days at 250 rpm. Northern blot hybridization of P32 labeled DNA probes to indicated target genes. Ethidium bromide-stained rRNA visualized as control.
Figure 8. Secondary metabolite analysis of wild type, ΔhapX and ΔsreA mutants in an iron-dependent manner normalized to iron replete conditions. The strains were grown for 5 days in GMM supplemented with indicated concentrations of iron for 5 days at 25°C. Secondary metabolites were analyzed as described in material and Methods and plotted as percent of the amount produced at 0.03 mM for each metabolite in each strain individually. Error bars indicate standard deviations for triplicates of each strain. Asterisks indicate significance as calculated by Tukey's test following ANOVA: **p < 0.0005 and *p < 0.001.
Both SreA and HapX loss perturbed metabolite production, reflected in transcriptional regulation in the sreA but not hapX mutant. In contrast to SreA negative regulation of HAS, fumagillin and pseurotin A production and gene expression was decreased in the ΔsreA mutant with increasing iron concentrations (Figures 7, 8). Production of fumagillin and pseurotin A, as well as fumitremorgin, was greatly reduced in ΔsreA in high-iron conditions (Figure 9). All four secondary metabolites were differentially produced in the ΔhapX mutant depending on iron concentration (Figure 8) with production of fumagillin and pseurotin significantly increased in ΔhapX compared to the WT under no-iron and iron-replete conditions (Figure 9).
Figure 9. Secondary metabolite analysis of wild type, ΔhapX and ΔsreA mutants in an iron-dependent manner normalized to the wild type. The strains were grown for 5 days in GMM supplemented with indicated concentrations of iron for 5 days at 25°C. Secondary metabolites were analyzed as described in Material and Methods and plotted as percent of the wild-type amount produced by ΔhapX and ΔsreA, respectively, for each iron condition. Error bars indicate standard deviations for triplicates of each strain. Asterisks indicate significance as calculated by Tukey's test following ANOVA: ***p < 0.0001 and *p < 0.05.
A curious observation was made for fumitremorgin C where its production followed much of the same trend as fumagillin and pseurotin A in both ΔsreA ΔhapX mutants (Figures 8, 9) but, unlike fumagillin and pseurotin, stayed unaffected by increasing iron concentration in WT (Figure 8). Interestingly—and in contrast to fumagillin and pseurotin A cluster genes - expression of ftmH was highest under iron-deficient conditions, despite the lack of product formation under those conditions (Figure 7). Notably, fumitremorgin C production is significantly increased in the ΔhapX mutant compared to the WT under no-iron but not iron excess conditions, indicating that HapX functions as a repressor under these conditions (Figure 9).
HAS has Global Influence on Genes Involved in Secondary Metabolism and Iron Acquisition and Consumption
The above data strongly suggested HAS to be a player in an iron acquisition and consumption circuitry involving secondary metabolite production. In order to more thoroughly interrogate these hypotheses, we performed RNA-seq analysis comparing overexpression hasA (e.g., OE::hasA = high HAS) to WT (low HAS) or to OE::hasA/ΔhasD (no HAS). Compared to the WT, 583 and 435 genes were up- and down-regulated in the OE::hasA strain, respectively. When comparing the OE::hasA strain to the strain OE::hasA/ΔhasD, 767 were up- and 394 genes were down-regulated in the strain able to produce HAS (Figure 10; Dataset S1). A common set of 404 and 291 genes to up- and down-regulated, respectively, was identified between the two comparisons (Figure 10; Dataset S1). Enrichment analysis of these common gene sets of the OE::hasA strain compared to both the WT and the OE::hasA/ΔhasD strains, shows that “secondary metabolism” is a category to be both positively and negatively affected (Figures S1, S2; Dataset S2). Analyzing genes from this category more closely, based on characterized and predicted (Khaldi et al., 2010; Medema et al., 2011; Inglis et al., 2013) gene clusters, provided an accurate view on the genes affected. As expected, genes belonging to the HAS cluster were found to be up-regulated (Table 1).
Figure 10. RNA seq data analysis of gene sets differentially regulated comparing OE::hasA to WT and OW::hasA to OE::hasA/ΔhasD. (A) log(2)-log(2) scatter plot and distribution of genes comparing OE::hasA to WT and OW::hasA to OE::hasA/ΔhasD. Genes from the gliotoxin gene cluster are represented in green. (B) Number of genes at least two-fold differentially regulated comparing the indicated mutants in Venn diagrams. Numbers in overlapping areas resemble common set of genes in the indicated comparison.
Table 1. Secondary metabolite cluster genes differentially regulated comparing OE::hasA vs. WT and OE::hasA vs. OE::hasA/ΔhasD.
Additionally, genes from several predicted but uncharacterized clusters also are positively affected: five and eight genes from two uncharacterized NRPS-like gene clusters on chromosome III (AFUA_3g02670; 02680, 02685, 02690, and 02760) and V (AFUA_5g09960; 10040, 10090, 10160, 10180, 10200, 10210, and 10220), respectively, as well as one gene belonging to an uncharacterized PKS cluster on chromosome I (AFUA_1g17690), four genes belonging to an uncharacterized PKS gene cluster on chromosome III (AFUA_3g01400; 01460; 01500, and 01560) and three genes belonging to an uncharacterized PKS cluster on chromosome IV (AFUA_4g14530, 14570, and 14580) (Table 1). Interestingly, the uncharacterized PKS-encoding gene on chromosome I (AFUA_1g17740) as well as one of the uncharacterized PKS-encoding gene on chromosome III (AfuA_3g14700) and additional five genes from the same cluster (AFUA_3g14710, 14740, 14750; 14760, and 14770) were only found to be up-regulated comparing the OE::hasA strain to OE::hasA/ΔhasD, in which neither HAS nor any intermediate is produced (Table 1). Similarly, three (fgaPT1/easL, fgaDH/easD, and fgaP450-2/easM) of the 11 characterized genes from the fumigaclavine gene cluster (Unsold and Li, 2006; Wallwey et al., 2010; Robinson and Panaccione, 2012; Wallwey et al., 2012) as well as one uncharacterized NRPS-encoding gene [AFUA_3g13730; nrps6/pesG; (Cramer et al., 2006b; Stack et al., 2007)] were only found to be up-regulated when comparing OE::hasA to OE::hasA/ΔhasD (Table 1).
Of the secondary metabolite gene clusters found to be down-regulated in the common gene set, all of the 13 gliotoxin cluster genes (gliZ, gliI, gliJ gliP, gliC, gliM, gliG, gliK, gliA, gliN, gliF, gliT, and gliH) (Gardiner et al., 2005; Balibar and Walsh, 2006; Schrettl et al., 2010b; Davis et al., 2011; Forseth et al., 2011; Scharf et al., 2012a,b) (Figure 10), as well as the adjacently located NRPS-encoding gene nrps9/pesJ [AFUA_6g09610; (Cramer et al., 2006b; Stack et al., 2007)], seven (ftmA, ftmB, ftmC, ftmD, ftmE, ftmG, and ftmH) of the eight fumitremorgin cluster genes (Maiya et al., 2006; Grundmann et al., 2008; Kato et al., 2009; Steffan et al., 2009; Wiemann et al., 2013), and three (pyr6, pyr7, and pyr8) of the nine pyripyropene cluster genes (Itoh et al., 2012) were identified (Table 1). Additionally, one gene (AFUA_1g10370) adjacent to the NRPS-encoding gene nrps1/pes1/pesB [AFUA_1g10380; (Cramer et al., 2006b; Reeves et al., 2006; O'Hanlon et al., 2012)] and the uncharacterized PKS-encoding gene (AFUA_3g02570) were also negatively affected (Table 1).
Overall, as noted above, the effect of hasA overexpression was dampened in the absence of hasD, with the fold differences in the secondary metabolite clusters often being smaller and the number of differentially expressed genes fewer. These results suggest that for the majority of gene clusters HAS itself, and not just HasA, plays a role in regulation of secondary metabolite biosynthetic genes.
Another category that we found to be enriched specifically among the genes positively affected comparing OE::hasA to the WT and OE::hasA to OE::hasA/ΔhasD belong to the functional group “heavy metal ion transport” (Figure S1; Table S2). This included the metalloreductase-encoding gene fre2/freB, the putative ferric chelate reductase-encoding gene fre7, the putative zinc transporter fetD/fet4, and two putative siderophore transporters (AFUA_3g13670 and AFUA_7g04730) (Table 2) all previously identified to be induced by iron or zinc starvation (Schrettl et al., 2008; Yasmin et al., 2009; Blatzer et al., 2011a,b). Two of the genes in this category were also induced by exposure to neutrophils (Sugui et al., 2008), one of them encoding a putative high affinity copper transporter (AFUA_2g03730) and the other one encoding a putative metalloreductase (AFUA_6g02820) (Table 2). These results support the hypothesis that HAS is a component of an iron-centered regulon, that when over-produced induced iron starvation responsive genes.
Table 2. Genes significantly up-regulated comparing OE::hasA vs. WT and OE::hasA vs. OE::hasA/ΔhasD of the category “heavy metal ion transport.”
On the other hand, several genes repressed by iron sufficiency in an SreA-dependent manner (Schrettl et al., 2008), were only identified to be up-regulated when comparing the OE::hasA strain to OE::hasA/ΔhasD (Table 3). Specifically, the gene cluster located on chromosome III (AfuA_3g03390-03440) involved in siderophore production is up-regulated in this comparison. Genes identified in this cluster encode for the hydroxyornithine transacylase SidF, the fusarinine C NRPS SidD, the ABC transporter SitT and the fusarinine C acyltransferase SidG (Table 3).
Table 3. Genes significantly up-regulated comparing OE::hasA vs. OE::hasA/ΔhasD involved in siderophore biosynthesis and SreA-repressed.
Several of the genes negatively affected in both comparisons (OE::hasA to WT and OE::hasA to OE::hasA/ΔhasD) belong to Fe-dependent enzymes involved in primary metabolism (Figure S2; Dataset S2), supporting the hypothesis that overproduction of HAS leads these strains to experience iron-deficiency. Among those genes were eight putative oxidoreductases (AFUA_1g17430; AFUA_2g03820; AFUA_2g09355; AFUA_3g00150; AFUA_3g01070; AFUA_4g13780; AFUA_5g09720; AFUA_6g03290) (Table 4). In line with this finding we observed an up-regulation of the cytochrome C-encoding gene (cycA) in northern blot analysis (Figure 5B). Similar to what was observed for the genes belonging to secondary metabolism, the global effects of hasA overexpression were also lessened in the absence of hasD. In this case, only 101 and 194 genes were up- or down-regulated in the OE::hasA/ΔhasD strain compared to WT, with less than a third of those having differences greater than four-fold, supporting our the hypothesis that HAS itself causes the major effects observed (Dataset S1).
Table 4. Genes significantly down-regulated comparing comparing OE::hasA vs. WT and OE::hasA vs. OE::hasA/ΔhasD.
Initial characterization of the HAS biosynthetic pathway resulted in two significant findings, one being that the HAS bound iron(III) and another that HAS and the cluster transcription factor, hasA, might be involved in production of other small molecules (Yin et al., 2013). The binding with iron led to the hypothesis that HAS production might impact iron homeostasis in A. fumigatus whereas the latter observation spoke of a global influence on secondary metabolite production. To investigate these observations further, we coupled a genomic transcriptional profile analysis of hasA and hasD mutants with biochemical examination of these strains and iron homeostasis mutants, in the process uncovering a heretofore uncharacterized iron-responsive secondary metabolite network in A. fumigatus. A model of our findings is illustrated in Figure 11 where we place HAS as a component of the iron regulon that ties into global regulation of secondary metabolic clusters.
Figure 11. Iron-dependent regulation in Aspergillus fumigatus. The two iron-dependent transcription factors SreA and HapX are in the center of an iron-dependent feed-back loop controlling iron-acquisition and -consumption pathways (Haas, 2012). Iron availability induces production of secondary metabolism, with SreA and HapX acting as activator and repressor. Production of Terezine D (HAS precursor) is repressed by SreA. Overproduction of HAS caused iron starvation conditions resulting in deregulation of iron-acquisition and -consumption pathways, most likely through SreA and HapX. Overproduction of HAS has global influence on expression of secondary metabolite gene clusters.
Characterization of HAS mutants on iron deficient and iron replete medium was essential for understanding the impact of this metabolite on fungal physiology. Our first query was to determine if HAS exhibited any properties of a siderophore and thus compared growth of hasA mutants to that of ΔsidA, lacking the L-ornithine-N5-monooxygenase SidA, the first committed step required for both extracellular and intracellular siderophore biosynthesis in A. fumigatus (Schrettl et al., 2004a). Deletion of hasA had no gross impact on growth on medium lacking iron and OE::hasA in a ΔsidA background did not restore growth on iron deficient media (Figure 1). However, there were several significant aberrancies in both the OE::hasA and the OE::hasA/ΔhasG (making the penultimate iron-binding precursor astechrome) strains grown on various concentrations of iron whereas the OE::hasA/ΔhasD strain (unable to produce a iron(III) binding molecule) demonstrated a wild type growth pattern.
First among the phenotypes in the OE::hasA and OE::hasA/ΔhasG strains was the accumulation of red pigmentation in low (0.03 mM) and high (1.5 mM) iron supplementation with the red intensity proportional to the iron content (Figure 2). Pigmentation was absent in the OE::hasA/ΔhasD strain. The OE::hasA strain also showed a statistically significant submersed growth deficit in iron concentrations below 1 mM (Figure 3) which was not obvious on agar medium (Figure 1). The concomitant increase in extracellular siderophore content in these same conditions coupled with decreased intracellular siderophore content at 0 mM iron in the OE::hasA strain (Figure 4) suggested that HAS might act as a one-way sink for iron, rendering it inaccessible to the fungus (Figure 11).
This hypothesis was supported by assessing expression of key genes associated with iron starvation. In addition to siderophore production, A. fumigatus also hosts an active reductive iron assimilation system with FtrA acting as a high affinity iron permease (Schrettl et al., 2004a). Both SidA and FtrA are activated during conditions of iron starvation (Schrettl et al., 2008, 2010a) and their expression was upregulated in the OE::hasA and OE::hasA/ΔhasG mutants. On the other hand, genes encoding proteins known to consume iron, such as cytochrome C (cycA) in the respiratory pathway, are down regulated in iron starvation conditions (Oberegger et al., 2002). Here cycA mRNA was down regulated in the HAS over-production strains (Figure 5A), yet another indication that excess HAS or astechrome results in iron deficiency. Furthermore, supporting this hypothesis, the GATA-type transcription factor, SreA, that acts as repressor under iron excess conditions (Schrettl et al., 2008) was significantly down-regulated in strains over-producing HAS (Figure 5B). Global RNA-seq analysis under iron-replete conditions support these findings, as genes involved in iron acquisition that are part of the SreA-regulon are positively affected when HAS is overproduced compared to a strain unable to produce HAS or any intermediate of the pathway (Figure 11; Tables 2, 3). Interestingly, the RNA-seq data not only suggest an up-regulation of genes involved in iron acquisition but also genes involved in copper and zinc homeostasis (Table 2). A relationship between iron and zinc as well as iron and copper homeostasis was previously described in A. fumigatus, where iron starvation caused an up-regulation of zinc transporters (Yasmin et al., 2009). Up-regulation of copper transporters under iron deficient conditions can best be explained by the reductive iron-uptake system of A. fumigatus which is copper dependent (Blatzer et al., 2011b).
Additionally to the regulatory effects on iron acquisition and consumption pathways caused by HAS production, our RNA-seq data implicate that production of HAS and not just the transcription factor HasA positively and negatively influences expression of specific A. fumigatus secondary metabolite gene clusters (Figure 11; Table 1). Interestingly, while the has genes were still overexpressed in OE::hasA compared to OE::hasA/ΔhasD, the fold overexpression decreased compared to OE::hasA vs. WT suggesting that either HAS itself or HasD had some influence on has expression (Table 1). Such a situation has been noted for the gliotoxin cluster in A. fumigatus and the bikaverin gene cluster in F. fujikuroi, where cluster gene expression is influenced by formation of the final product and the respective key enzymes GliP (NRPS) and Bik1 (PKS) (Cramer et al., 2006a; Wiemann et al., 2009). The observed effect of cluster cross-talk, where production of one metabolite influences the abundance of another one, was previously reported for A. fumigatus. Specifically, deletion of either the NRPS-encoding gene responsible for fumiquinazoline production (nrps11/pesL/fqzC/fmqC) or nrps1/pes1/pesB, involved in production of an uncharacterized non-ribosomal peptide, had significant influence on production of fumigaclavines and fumitremorgins, the former depleted and the latter found to be increased (O'Hanlon et al., 2012). It is interesting to note that the same trend was observed when HAS production was increased (Table 1). Cross-talk between gliotoxin and ergothioneine has also been reported for A. fumigatus and is suggested to be caused by response to a change in redox status of the cells rather than changed precursor pools (Gallagher et al., 2012). In line with these observations, exogenously added gliotoxin resulted in differentially expressed proteins involved in redox homeostasis (Carberry et al., 2012). A relation between iron and redox homeostasis is described in pro- and eukaryotic organisms, where iron-containing proteins function as reactive oxygen species (ROS)-sensors on the one hand (Outten and Theil, 2009; Outten and Albetel, 2013) and ROS-detoxifying enzymes on the other hand (Oberegger et al., 2000, 2001; Haas, 2012). Since HAS had a major impact on iron homeostasis, it is feasible that the effects observed on secondary metabolism are a consequence of the changes in iron availability through modification of redox homeostasis.
This hypothesis, that the global impact of HAS could be in large part mediated by iron availability, was supported by both the effect of HapX and SreA loss on secondary metabolite output as well as iron feeding impacts on small molecule production (Figures 7, 9). With the exception of fumitremorgin production, all analyzed metabolites including terezine D (the first stable intermediate in the HAS pathway) were significantly increased when ambient iron concentrations were increased from 0.03 to 1 mM, suggesting that SreA might execute activating effects under iron excess conditions (Figure 8). Supporting this hypothesis, this induction of two secondary metabolites (pseurotin A and fumagillin) was lost in the ΔsreA strain (Figures 7, 8) and production was significantly reduced compared to the wild type in 1 mM iron (Figure 9). However, production of terezine D was significantly increased in the ΔsreA strain compared to the WT in 0.03 and 1 mM iron, indicating that SreA negatively affects HAS production under these conditions, similarly to iron acquisition pathways identified by Schrettl et al. (2008) (Figure 11). Deletion of sreA was also shown to induce increased iron toxicity due to uncoordinated iron uptake (Schrettl et al., 2008), possibly yet another route by which sreA loss could affect secondary metabolism.
As predicted by our model (Figure 11), deficiency of HapX alters secondary metabolite production mainly under iron deplete conditions, where it seems to function as a repressor (Figures 6, 8, 9). HapX is part of the CCAAT-Box binding complex in filamentous fungi (Kato, 2005; Hortschansky et al., 2007), and was shown to be an inducer of siderophore production in A. fumigatus (Schrettl et al., 2010a) and is directly implicated to regulate the penicillin gene cluster in A. nidulans (Litzka et al., 1996). This work presents a role for HapX in repressing expression and subsequent production of several secondary metabolite clusters.
Summarizing our combined data, we show here that production of the iron(III)-chelating metabolite HAS is strongly regulated by the iron-induced transcription factor SreA, and HAS production feeds into regulation of sreA expression, thereby eliciting an iron-depletion response in iron-replete conditions. These findings suggest overexpression of HAS and deficiency of both iron-responsive elements HapX and SreA have significant effects on secondary metabolism, thus placing iron availability as one environmental trigger of secondary metabolism (Figure 11). This impact of iron on secondary metabolism could be direct (many secondary metabolite enzymes require iron as a cofactor) or through changes in primary metabolism, possibly through alterations in amino acid pools that have previously been reported to be coordinated by iron-availability in a Hapx-dependent manner in A. fumigatus (Schrettl et al., 2008). Precedence for primary metabolism affecting secondary metabolism comes from F. fujikuroi, where mutations of specific amino acids in the glutamine synthetase had significant effects on amino acid composition and secondary metabolite production (Wagner et al., 2013). It is tempting to speculate that the effects of HAS and iron gradients on secondary metabolism are a reflection not only of iron availability but of accessible primary metabolite pools and/or redox homeostasis.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Philipp Wiemann was supported by National Institutes of Health Grant NIH R01 Al065728-01 to Nancy P. Keller. Beatrix E. Lechner and Hubertus Haas were supported by the Austrian Science Foundation Grant FWF I1346-B21. Joshua A. Baccile and Frank C. Schroeder were supported by National Institutes of Health Chemistry-Biology Interface Training Predoctoral Training Program T32GM008500.
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fmicb.2014.00530/abstract
Figure S1 | Fold enrichment of functional categories of the common gene set at least two-fold up-regulated in OE::hasA compared to WT and OE::hasA/ΔhasD. Gene sets were analyzed using the Munich Information Center for Protein Sequences Functional Catalog as described in Material and Methods.
Figure S2 | Fold enrichment of functional categories of the common gene set at least two-fold down-regulated in OE::hasA compared to WT and OE::hasA/ΔhasD. Gene sets were analyzed using the Munich Information Center for Protein Sequences Functional Catalog as described in Material and Methods.
Figure S3 | Southern analysis of ΔhasA, ΔsidA and OE::hasA/ΔsidA strains. Genomic DNA of indicated strains was subject to restriction enzyme digestion using either BamHI or HindIII. WT = parental wild type. For deletion of hasA, the WT shows the expected band of 214 bp and transformants T6 and T9 show the expected band of 6300 bp. For deletion of sidA, the WT shows the expected band of 4000 bp and transformants T1–5 (ΔsidA) and T1–4 (OE::hasA/ΔsidA) show the expected band of 2300 bp.
Table S1 | Strains used in this study.
Table S2 | Oligonucleotides used in this study.
Dataset S1 | Excel file containing expression levels based on normalized reads per kilobase per million mapped reads and log(2) values of WT, OE::hasA, and OE::hasA/ΔhasD.
Dataset S2 | Excel file containing the output of the functional enrichment analysis of the common gene set of genes at least two-fold up- and down-regulated in OE::hasA compared to WT and OE::hasA/ΔhasD.
Dataset S3 | Excel file containing areas under the curve for the secondary metabolites analyzed in this study as described in Material and Methods.
An, Z., Mei, B., Yuan, W. M., and Leong, S. A. (1997). The distal GATA sequences of the sid1 promoter of Ustilago maydis mediate iron repression of siderophore production and interact directly with Urbs1, a GATA family transcription factor. EMBO J. 16, 1742–1750. doi: 10.1093/emboj/16.7.1742
Arai, K., Sato, S., Shimizu, S., Nitta, K., and Yamamoto, Y. (1981). Metabolic products of Aspergillus terreus. VII. Astechrome: an iron-containing metabolite of the Strain IFO 6123. Chem. Pharm. Bull. 29, 1510–1517. doi: 10.1248/cpb.29.1510
Balibar, C. J., and Walsh, C. T. (2006). GliP, a multimodular nonribosomal peptide synthetase in Aspergillus fumigatus, makes the diketopiperazine scaffold of gliotoxin. Biochemistry 45, 15029–15038. doi: 10.1021/bi061845b
Berthier, E., Lim, F. Y., Deng, Q., Guo, C. J., Kontoyiannis, D. P., Wang, C. C., et al. (2013). Low-volume toolbox for the discovery of immunosuppressive fungal secondary metabolites. PLoS Pathog. 9:e1003289. doi: 10.1371/journal.ppat.1003289
Blatzer, M., Barker, B. M., Willger, S. D., Beckmann, N., Blosser, S. J., Cornish, E. J., et al. (2011a). SREBP coordinates iron and ergosterol homeostasis to mediate triazole drug and hypoxia responses in the human fungal pathogen Aspergillus fumigatus. PLoS Genet. 7:e1002374. doi: 10.1371/journal.pgen.1002374
Blatzer, M., Binder, U., and Haas, H. (2011b). The metalloreductase FreB is involved in adaptation of Aspergillus fumigatus to iron starvation. Fungal Genet. Biol. 48, 1027–1033. doi: 10.1016/j.fgb.2011.07.009
Bok, J. W., Chung, D., Balajee, S. A., Marr, K. A., Andes, D., Nielsen, K. F., et al. (2006). GliZ, a transcriptional regulator of gliotoxin biosynthesis, contributes to Aspergillus fumigatus virulence. Infect. Immun. 74, 6761–6768. doi: 10.1128/IAI.00780-06
Bullard, J. H., Purdom, E., Hansen, K. D., and Dudoit, S. (2010). Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 11:94. doi: 10.1186/1471-2105-11-94
Calvo, A. M., Bok, J., Brooks, W., and Keller, N. P. (2004). veA is required for toxin and sclerotial production in Aspergillus parasiticus. Appl. Environ. Microbiol. 70, 4733–4739. doi: 10.1128/AEM.70.8.4733-4739.2004
Carberry, S., Molloy, E., Hammel, S., O'Keeffe, G., Jones, G. W., Kavanagh, K., et al. (2012). Gliotoxin effects on fungal growth: mechanisms and exploitation. Fungal Genet. Biol. 49, 302–312. doi: 10.1016/j.fgb.2012.02.003
Cramer, R. A. J., Gamcsik, M. P., Brooking, R. M., Najvar, L. K., Kirkpatrick, W. R., Patterson, T. F., et al. (2006a). Disruption of a nonribosomal peptide synthetase in Aspergillus fumigatus eliminates gliotoxin production. Eukaryot. Cell 5, 972–980. doi: 10.1128/EC.00049-06
Cramer, R. A. J., Stajich, J. E., Yamanaka, Y., Dietrich, F. S., Steinbach, W. J., and Perfect, J. R. (2006b). Phylogenomic analysis of non-ribosomal peptide synthetases in the genus Aspergillus. Gene 383, 24–32. doi: 10.1016/j.gene.2006.07.008
da Silva Ferreira, M. E., Kress, M. R., Savoldi, M., Goldman, M. H., Hartl, A., Heinekamp, T., et al. (2006). The akuB(KU80) mutant deficient for nonhomologous end joining is a powerful tool for analyzing pathogenicity in Aspergillus fumigatus. Eukaryot. Cell 5, 207–211. doi: 10.1128/EC.5.1.207-211.2006
Davis, C., Carberry, S., Schrettl, M., Singh, I., Stephens, J. C., Barry, S. M., et al. (2011). The role of glutathione S-transferase GliG in gliotoxin biosynthesis in Aspergillus fumigatus. Chem. Biol. 18, 542–552. doi: 10.1016/j.chembiol.2010.12.022
Eisendle, M., Oberegger, H., Zadra, I., and Haas, H. (2003). The siderophore system is essential for viability of Aspergillus nidulans: functional analysis of two genes encoding l-ornithine N5-monooxygenase (sidA) and a non-ribosomal peptide synthetase (sidC). Mol. Microbiol. 49, 359–375. doi: 10.1046/j.1365-2958.2003.03586.x
Fluckinger, M., Haas, H., Merschak, P., Glasgow, B. J., and Redl, B. (2004). Human tear lipocalin exhibits antimicrobial activity by scavenging microbial siderophores. Antimicrob. Agents Chemother. 48, 3367–3372. doi: 10.1128/AAC.48.9.3367-3372.2004
Forseth, R. R., Fox, E. M., Chung, D., Howlett, B. J., Keller, N. P., and Schroeder, F. C. (2011). Identification of cryptic products of the gliotoxin gene cluster using NMR-based comparative metabolomics and a model for gliotoxin biosynthesis. J. Am. Chem. Soc. 133, 9678–9681. doi: 10.1021/ja2029987
Gallagher, L., Owens, R. A., Dolan, S. K., O'Keeffe, G., Schrettl, M., Kavanagh, K., et al. (2012). The Aspergillus fumigatus protein GliK protects against oxidative stress and is essential for gliotoxin biosynthesis. Eukaryot. Cell 11, 1226–1238. doi: 10.1128/EC.00113-12
Gardiner, D. M., Waring, P., and Howlett, B. J. (2005). The epipolythiodioxopiperazine (ETP) class of fungal toxins: distribution, mode of action, functions and biosynthesis. Microbiology 151, 1021–1032. doi: 10.1099/mic.0.27847-0
Gauthier, G. M., Sullivan, T. D., Gallardo, S. S., Brandhorst, T. T., Vanden Wymelenberg, A. J., Cuomo, C. A., et al. (2010). SREB, a GATA transcription factor that directs disparate fates in Blastomyces dermatitidis including morphogenesis and siderophore biosynthesis. PLoS Pathog. 6:e1000846. doi: 10.1371/journal.ppat.1000846
Greenshields, D. L., Liu, G., Feng, J., Selvaraj, G., and Wei, Y. (2007). The siderophore biosynthetic gene SID1, but not the ferroxidase gene FET3, is required for full Fusarium graminearum virulence. Mol. Plant Pathol. 8, 411–421. doi: 10.1111/j.1364-3703.2007.00401.x
Grundmann, A., Kuznetsova, T., Afiyatullov, S. S., and Li, S. M. (2008). FtmPT2, an N-prenyltransferase from Aspergillus fumigatus, catalyses the last step in the biosynthesis of fumitremorgin B. Chembiochem 9, 2059–2063. doi: 10.1002/cbic.200800240
Gsaller, F., Eisendle, M., Lechner, B. E., Schrettl, M., Lindner, H., Muller, D., et al. (2012). The interplay between vacuolar and siderophore-mediated iron storage in Aspergillus fumigatus. Metallomics 4, 1262–1270. doi: 10.1039/c2mt20179h
Haas, H. (2003). Molecular genetics of fungal siderophore biosynthesis and uptake: the role of siderophores in iron uptake and storage. Appl. Microbiol. Biotechnol. 62, 316–330. doi: 10.1007/s00253-003-1335-2
Haas, H., Zadra, I., Stoffler, G., and Angermayr, K. (1999). The Aspergillus nidulans GATA factor SREA is involved in regulation of siderophore biosynthesis and control of iron uptake. J. Biol. Chem. 274, 4613–4619. doi: 10.1074/jbc.274.8.4613
Heinekamp, T., Thywissen, A., Macheleidt, J., Keller, S., Valiante, V., and Brakhage, A. A. (2012). Aspergillus fumigatus melanins: interference with the host endocytosis pathway and impact on virulence. Front. Microbiol. 3:440. doi: 10.3389/fmicb.2012.00440
Hissen, A. H., Wan, A. N., Warwas, M. L., Pinto, L. J., and Moore, M. M. (2005). The Aspergillus fumigatus siderophore biosynthetic gene sidA, encoding l-ornithine N5-oxygenase, is required for virulence. Infect. Immun. 73, 5493–5503. doi: 10.1128/IAI.73.9.5493-5503.2005
Hortschansky, P., Eisendle, M., Al-Abdallah, Q., Schmidt, A. D., Bergmann, S., Thon, M., et al. (2007). Interaction of HapX with the CCAAT-binding complex–a novel mechanism of gene regulation by iron. EMBO J. 26, 3157–3168. doi: 10.1038/sj.emboj.7601752
Hwang, L. H., Mayfield, J. A., Rine, J., and Sil, A. (2008). Histoplasma requires SID1, a member of an iron-regulated siderophore gene cluster, for host colonization. PLoS Pathog. 4:e1000044. doi: 10.1371/journal.ppat.1000044
Hwang, L. H., Seth, E., Gilmore, S. A., and Sil, A. (2012). SRE1 regulates iron-dependent and -independent pathways in the fungal pathogen Histoplasma capsulatum. Eukaryot. Cell 11, 16–25. doi: 10.1128/EC.05274-11
Inglis, D. O., Binkley, J., Skrzypek, M. S., Arnaud, M. B., Cerqueira, G. C., Shah, P., et al. (2013). Comprehensive annotation of secondary metabolite biosynthetic genes and gene clusters of Aspergillus nidulans, A. fumigatus, A. niger and A. oryzae. BMC Microbiol. 13:91. doi: 10.1186/1471-2180-13-91
Itoh, T., Kushiro, T., and Fujii, I. (2012). Reconstitution of a secondary metabolite biosynthetic pathway in a heterologous fungal host. Methods Mol. Biol. 944, 175–182. doi: 10.1007/978-1-62703-122-6_12
Jahn, B., Langfelder, K., Schneider, U., Schindel, C., and Brakhage, A. A. (2002). PKSP-dependent reduction of phagolysosome fusion and intracellular kill of Aspergillus fumigatus conidia by human monocyte-derived macrophages. Cell. Microbiol. 4, 793–803. doi: 10.1046/j.1462-5822.2002.00228.x
Kato, M. (2005). An overview of the CCAAT-box binding factor in filamentous fungi: assembly, nuclear translocation, and transcriptional enhancement. Biosci. Biotechnol. Biochem. 69, 663–672. doi: 10.1271/bbb.69.663
Kato, N., Suzuki, H., Takagi, H., Asami, Y., Kakeya, H., Uramoto, M., et al. (2009). Identification of cytochrome P450s required for fumitremorgin biosynthesis in Aspergillus fumigatus. Chembiochem 10, 920–928. doi: 10.1002/cbic.200800787
Khaldi, N., Seifuddin, F. T., Turner, G., Haft, D., Nierman, W. C., Wolfe, K. H., et al. (2010). SMURF: genomic mapping of fungal secondary metabolite clusters. Fungal Genet. Biol. 47, 736–741. doi: 10.1016/j.fgb.2010.06.003
Knox, B. P., Deng, Q., Rood, M., Eickhoff, J. C., Keller, N. P., and Huttenlocher, A. (2014). Distinct innate immune phagocyte responses to Aspergillus fumigatus conidia and hyphae in zebrafish larvae. Eukaryot. Cell. 13, 1266–1277. doi: 10.1128/EC.00080-14
Leal, S. M. J., Roy, S., Vareechon, C., Carrion, S. D., Clark, H., Lopez-Berges, M. S., et al. (2013). Targeting iron acquisition blocks infection with the fungal pathogens Aspergillus fumigatus and Fusarium oxysporum. PLoS Pathog. 9:e1003436. doi: 10.1371/journal.ppat.1003436
Litzka, O., Then Bergh, K., and Brakhage, A. A. (1996). The Aspergillus nidulans penicillin-biosynthesis gene aat (penDE) is controlled by a CCAAT-containing DNA element. Eur. J. Biochem. 238, 675–682. doi: 10.1111/j.1432-1033.1996.0675w.x
Lopez-Berges, M. S., Hera, C., Sulyok, M., Schafer, K., Capilla, J., Guarro, J., et al. (2013). The velvet complex governs mycotoxin production and virulence of Fusarium oxysporum on plant and mammalian hosts. Mol. Microbiol. 87, 49–65. doi: 10.1111/mmi.12082
Maiya, S., Grundmann, A., Li, S. M., and Turner, G. (2006). The fumitremorgin gene cluster of Aspergillus fumigatus: identification of a gene encoding brevianamide F synthetase. Chembiochem 7, 1062–1069. doi: 10.1002/cbic.200600003
Medema, M. H., Blin, K., Cimermancic, P., de Jager, V., Zakrzewski, P., Fischbach, M. A., et al. (2011). antiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 39, W339–W346. doi: 10.1093/nar/gkr466
Oberegger, H., Schoeser, M., Zadra, I., Abt, B., and Haas, H. (2001). SREA is involved in regulation of siderophore biosynthesis, utilization and uptake in Aspergillus nidulans. Mol. Microbiol. 41, 1077–1089. doi: 10.1046/j.1365-2958.2001.02586.x
Oberegger, H., Schoeser, M., Zadra, I., Schrettl, M., Parson, W., and Haas, H. (2002). Regulation of freA, acoA, lysF, and cycA expression by iron availability in Aspergillus nidulans. Appl. Environ. Microbiol. 68, 5769–5772. doi: 10.1128/AEM.68.11.5769-5772.2002
Oberegger, H., Zadra, I., Schoeser, M., and Haas, H. (2000). Iron starvation leads to increased expression of Cu/Zn-superoxide dismutase in Aspergillus. FEBS Lett. 485, 113–116. doi: 10.1016/S0014-5793(00)02206-7
Oglesby-Sherrouse, A. G., Djapgne, L., Nguyen, A. T., Vasil, A. I., and Vasil, M. L. (2014). The complex interplay of iron, biofilm formation, and mucoidy affecting antimicrobial resistance of Pseudomonas aeruginosa. Pathog. Dis. 70, 307–320. doi: 10.1111/2049-632X.12132
O'Hanlon, K. A., Gallagher, L., Schrettl, M., Jochl, C., Kavanagh, K., Larsen, T. O., et al. (2012). Nonribosomal peptide synthetase genes pesL and pes1 are essential for Fumigaclavine C production in Aspergillus fumigatus. Appl. Environ. Microbiol. 78, 3166–3176. doi: 10.1128/AEM.07249-11
Oide, S., Moeder, W., Krasnoff, S., Gibson, D., Haas, H., Yoshioka, K., et al. (2006). NPS6, encoding a nonribosomal peptide synthetase involved in siderophore-mediated iron metabolism, is a conserved virulence determinant of plant pathogenic ascomycetes. Plant Cell 18, 2836–2853. doi: 10.1105/tpc.106.045633
Outten, C. E., and Albetel, A. N. (2013). Iron sensing and regulation in Saccharomyces cerevisiae: ironing out the mechanistic details. Curr. Opin. Microbiol. 16, 662–668. doi: 10.1016/j.mib.2013.07.020
Palmer, J. M., Perrin, R. M., Dagenais, T. R., and Keller, N. P. (2008). H3K9 methylation regulates growth and development in Aspergillus fumigatus. Eukaryot. Cell 7, 2052–2060. doi: 10.1128/EC.00224-08
Reeves, E. P., Reiber, K., Neville, C., Scheibner, O., Kavanagh, K., and Doyle, S. (2006). A nonribosomal peptide synthetase (Pes1) confers protection against oxidative stress in Aspergillus fumigatus. FEBS J. 273, 3038–3053. doi: 10.1111/j.1742-4658.2006.05315.x
Ruepp, A., Zollner, A., Maier, D., Albermann, K., Hani, J., Mokrejs, M., et al. (2004). The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 32, 5539–5545. doi: 10.1093/nar/gkh894
Scharf, D. H., Chankhamjon, P., Scherlach, K., Heinekamp, T., Roth, M., Brakhage, A. A., et al. (2012a). Epidithiol formation by an unprecedented twin carbon-sulfur lyase in the gliotoxin pathway. Angew. Chem. Int. Ed. Engl. 51, 10064–10068. doi: 10.1002/anie.201205041
Scharf, D. H., Heinekamp, T., Remme, N., Hortschansky, P., Brakhage, A. A., and Hertweck, C. (2012b). Biosynthesis and function of gliotoxin in Aspergillus fumigatus. Appl. Microbiol. Biotechnol. 93, 467–472. doi: 10.1007/s00253-011-3689-1
Schrettl, M., Beckmann, N., Varga, J., Heinekamp, T., Jacobsen, I. D., Jochl, C., et al. (2010a). HapX-mediated adaption to iron starvation is crucial for virulence of Aspergillus fumigatus. PLoS Pathog. 6:e1001124. doi: 10.1371/journal.ppat.1001124
Schrettl, M., Carberry, S., Kavanagh, K., Haas, H., Jones, G. W., O'Brien, J., et al. (2010b). Self-protection against gliotoxin–a component of the gliotoxin biosynthetic cluster, GliT, completely protects Aspergillus fumigatus against exogenous gliotoxin. PLoS Pathog. 6:e1000952. doi: 10.1371/journal.ppat.1000952
Schrettl, M., Bignell, E., Kragl, C., Joechl, C., Rogers, T., Arst, H. N. J., et al. (2004a). Siderophore biosynthesis but not reductive iron assimilation is essential for Aspergillus fumigatus virulence. J. Exp. Med. 200, 1213–1219. doi: 10.1084/jem.20041242
Schrettl, M., Bignell, E., Kragl, C., Sabiha, Y., Loss, O., Eisendle, M., et al. (2007). Distinct roles for intra- and extracellular siderophores during Aspergillus fumigatus infection. PLoS Pathog. 3, 1195–1207. doi: 10.1371/journal.ppat.0030128
Schrettl, M., Kim, H. S., Eisendle, M., Kragl, C., Nierman, W. C., Heinekamp, T., et al. (2008). SreA-mediated iron regulation in Aspergillus fumigatus. Mol. Microbiol. 70, 27–43. doi: 10.1111/j.1365-2958.2008.06376.x
Shimizu, K., and Keller, N. P. (2001). Genetic involvement of a cAMP-dependent protein kinase in a G protein signaling pathway regulating morphological and chemical transitions in Aspergillus nidulans. Genetics 157, 591–600.
Steffan, N., Grundmann, A., Afiyatullov, S., Ruan, H., and Li, S. M. (2009). FtmOx1, a non-heme Fe(II) and alpha-ketoglutarate-dependent dioxygenase, catalyses the endoperoxide formation of verruculogen in Aspergillus fumigatus. Org. Biomol. Chem. 7, 4082–4087. doi: 10.1039/b908392h
Sugui, J. A., Kim, H. S., Zarember, K. A., Chang, Y. C., Gallin, J. I., Nierman, W. C., et al. (2008). Genes differentially expressed in conidia and hyphae of Aspergillus fumigatus upon exposure to human neutrophils. PLoS ONE 3:e2655. doi: 10.1371/journal.pone.0002655
Szewczyk, E., Nayak, T., Oakley, C. E., Edgerton, H., Xiong, Y., Taheri-Talesh, N., et al. (2006). Fusion PCR and gene targeting in Aspergillus nidulans. Nat. Protoc. 1, 3111–3120. doi: 10.1038/nprot.2006.405
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
Unsold, I. A., and Li, S. M. (2006). Reverse prenyltransferase in the biosynthesis of fumigaclavine C in Aspergillus fumigatus: gene expression, purification, and characterization of fumigaclavine C synthase FGAPT1. Chembiochem 7, 158–164. doi: 10.1002/cbic.200500318
Voisard, C., Wang, J., McEvoy, J. L., Xu, P., and Leong, S. A. (1993). urbs1, a gene regulating siderophore biosynthesis in Ustilago maydis, encodes a protein similar to the erythroid transcription factor GATA-1. Mol. Cell. Biol. 13, 7091–7100.
Wagner, D., Wiemann, P., Huss, K., Brandt, U., Fleissner, A., and Tudzynski, B. (2013). A sensing role of the glutamine synthetase in the nitrogen regulation network in Fusarium fujikuroi. PLoS ONE 8:e80740. doi: 10.1371/journal.pone.0080740
Wallner, A., Blatzer, M., Schrettl, M., Sarg, B., Lindner, H., and Haas, H. (2009). Ferricrocin, a siderophore involved in intra- and transcellular iron distribution in Aspergillus fumigatus. Appl. Environ. Microbiol. 75, 4194–4196. doi: 10.1128/AEM.00479-09
Wallwey, C., Heddergott, C., Xie, X., Brakhage, A. A., and Li, S. M. (2012). Genome mining reveals the presence of a conserved gene cluster for the biosynthesis of ergot alkaloid precursors in the fungal family Arthrodermataceae. Microbiology 158, 1634–1644. doi: 10.1099/mic.0.056796-0
Wallwey, C., Matuschek, M., Xie, X. L., and Li, S. M. (2010). Ergot alkaloid biosynthesis in Aspergillus fumigatus: conversion of chanoclavine-I aldehyde to festuclavine by the festuclavine synthase FgaFS in the presence of the old yellow enzyme FgaOx3. Org. Biomol. Chem. 8, 3500–3508. doi: 10.1039/c003823g
Wiemann, P., Albermann, S., Niehaus, E. M., Studt, L., von Bargen, K. W., Brock, N. L., et al. (2012). The Sfp-type 4'-phosphopantetheinyl transferase Ppt1 of Fusarium fujikuroi controls development, secondary metabolism and pathogenicity. PLoS ONE 7:e37519. doi: 10.1371/journal.pone.0037519
Wiemann, P., Guo, C. J., Palmer, J. M., Sekonyela, R., Wang, C. C. C., and Keller, N. P. (2013). Prototype of an intertwined secondary metabolite supercluster. Proc. Natl. Acad. Sci. U.S.A. 110, 17065–17070. doi: 10.1073/pnas.1313258110
Wiemann, P., Willmann, A., Straeten, M., Kleigrewe, K., Beyer, M., Humpf, H. U., et al. (2009). Biosynthesis of the red pigment bikaverin in Fusarium fujikuroi: genes, their function and regulation. Mol. Microbiol. 72, 931–946. doi: 10.1111/j.1365-2958.2009.06695.x
Xue, T., Nguyen, C. K., Romans, A., Kontoyiannis, D. P., and May, G. S. (2004). Isogenic auxotrophic mutant strains in the Aspergillus fumigatus genome reference strain AF293. Arch. Microbiol. 182, 346–353. doi: 10.1007/s00203-004-0707-z
Yasmin, S., Abt, B., Schrettl, M., Moussa, T. A., Werner, E. R., and Haas, H. (2009). The interplay between iron and zinc metabolism in Aspergillus fumigatus. Fungal Genet. Biol. 46, 707–713. doi: 10.1016/j.fgb.2009.05.003
Yin, W. B., Baccile, J. A., Bok, J. W., Chen, Y., Keller, N. P., and Schroeder, F. C. (2013). A nonribosomal peptide synthetase-derived iron(III) complex from the pathogenic fungus Aspergillus fumigatus. J. Am. Chem. Soc. 135, 2064–2067. doi: 10.1021/ja311145n
Zhou, L. W., Haas, H., and Marzluf, G. A. (1998). Isolation and characterization of a new gene, sre, which encodes a GATA-type regulatory protein that controls iron transport in Neurospora crassa. Mol. Gen. Genet. 259, 532–540. doi: 10.1007/s004380050845
Keywords: Aspergillus fumigatus, secondary metabolism, iron, HapX, SreA, hexadehydroastechrome, gene regulation
Citation: Wiemann P, Lechner BE, Baccile JA, Velk TA, Yin W-B, Bok JW, Pakala S, Losada L, Nierman WC, Schroeder FC, Haas H and Keller NP (2014) Perturbations in small molecule synthesis uncovers an iron-responsive secondary metabolite network in Aspergillus fumigatus. Front. Microbiol. 5:530. doi: 10.3389/fmicb.2014.00530
Received: 11 July 2014; Accepted: 23 September 2014;
Published online: 24 October 2014.
Edited by:Biswarup Mukhopadhyay, Virginia Tech, USA
Reviewed by:Daniel Panaccione, West Virginia University, USA
Anita Sil, University of California, San Francisco, USA
Copyright © 2014 Wiemann, Lechner, Baccile, Velk, Yin, Bok, Pakala, Losada, Nierman, Schroeder, Haas and Keller. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nancy P. Keller, Department of Medical Microbiology and Immunology, University of Wisconsin-Madison, 1550 Linden Drive, Madison, WI 53706, USA e-mail: firstname.lastname@example.org
†Present address: Wen-Bing Yin, State Key Laboratory of Mycology, The Institute of Microbiology, Chinese Academy of Sciences, Beijing, China