RNA-seq reveals the pan-transcriptomic impact of attenuating the gliotoxin self-protection mechanism in Aspergillus fumigatus
BMC Genomics volume 15, Article number: 894 (2014)
Aspergillus fumigatus produces a number of secondary metabolites, one of which, gliotoxin, has been shown to exhibit anti-fungal activity. Thus, A. fumigatus must be able to protect itself against gliotoxin. Indeed one of the genes in the gliotoxin biosynthetic gene cluster in A. fumigatus, gliT, is required for self-protection against the toxin- however the global self-protection mechanism deployed is unclear. RNA-seq was employed to identify genes differentially regulated upon exposure to gliotoxin in A. fumigatus wild-type and A. fumigatus ∆gliT, a strain that is hypersensitive to gliotoxin.
Deletion of A. fumigatus gliT resulted in altered expression of 208 genes (log2 fold change of 1.5) when compared to A. fumigatus wild-type, of which 175 genes were up-regulated and 33 genes were down-regulated. Expression of 164 genes was differentially regulated (log2 fold change of 1.5) in A. fumigatus wild-type when exposed to gliotoxin, consisting of 101 genes with up-regulated expression and 63 genes with down-regulated expression. Interestingly, a much larger number of genes, 1700, were found to be differentially regulated (log2 fold change of 1.5) in A. fumigatus ∆gliT when challenged with gliotoxin. These consisted of 508 genes with up-regulated expression, and 1192 genes with down-regulated expression. Functional Catalogue (FunCat) classification of differentially regulated genes revealed an enrichment of genes involved in both primary metabolic functions and secondary metabolism. Specifically, genes involved in gliotoxin biosynthesis, helvolic acid biosynthesis, siderophore-iron transport genes and also nitrogen metabolism genes and ribosome biogenesis genes underwent altered expression. It was confirmed that gliotoxin biosynthesis is induced upon exposure to exogenous gliotoxin, production of unrelated secondary metabolites is attenuated in A. fumigatus ∆gliT, while quantitative proteomic analysis confirmed disrupted translation in A. fumigatus ∆gliT challenged with exogenous gliotoxin.
This study presents the first global investigation of the transcriptional response to exogenous gliotoxin in A. fumigatus wild-type and the hyper-sensitive strain, ∆gliT. Our data highlight the global and extensive affects of exogenous gliotoxin on a sensitive strain devoid of a self-protection mechanism and infer that GliT functionality is required for the optimal biosynthesis of selected secondary metabolites in A. fumigatus.
Gliotoxin, a non-ribosomally synthesised peptide produced by Aspergillus fumigatus and related fungi, is redox-active, depletes cellular glutathione (GSH), causes DNA damage and protein modification, and is consequently capable of inhibiting mammalian, fungal and bacterial cell growth [1–5]. Conversely, it has also been demonstrated that gliotoxin presence protects against H2O2-induced oxidative stress in A. fumigatus and that gliotoxin can substitute for peroxiredoxin in mammalian cells to protect against similar oxidative stresses [6, 7]. In A. fumigatus, gliotoxin biosynthesis is encoded by a 13-gene cluster, gli, and it has been demonstrated that gliotoxin effects induction of many genes within the gli cluster in a positive feedback manner [2, 8, 9]. Thus, gliotoxin presence induces gli cluster activation via gliZ, a Zn2Cys6 binuclear transcription factor, and deletion of gliZ abolishes gliotoxin biosynthesis . Moreover, Forseth et al.  revealed that an additional nine metabolites, dependent on gliZ presence, are produced consequent to gliotoxin biosynthetic pathway functionality in A. fumigatus.
Gliotoxin exposure has been shown by qRT-PCR to either activate or induce increased expression of all genes in the gli cluster, as especially observed for A. fumigatus ∆gliP, deficient in the non-ribosomal peptide synthetase which mediates cyclo-L-Phe-L-Ser formation [9, 12]. Others have shown induction of gliG (a glutathione S-transferase), gliA (an MFS transporter) and gliT upon exposure of A. fumigatus wild-type to gliotoxin, by Northern analysis . However, definitive evidence of concomitant increased de novo gliotoxin production has not been forthcoming. Relatedly, it has been shown that transformation with A. fumigatus gliA confers resistance against exogenous gliotoxin upon Leptosphaeria maculans  while deletion of gliA in A. fumigatus renders it less resistant to exogenous gliotoxin , and Schrettl et al.  were the first to demonstrate increased GliT abundance by 2D-PAGE/MALDI-ToF analysis in A. fumigatus upon exposure to exogenous gliotoxin. However, within A. fumigatus, gliotoxin biosynthesis must be controlled to avoid manifestation of the deleterious affects of this reactive metabolite.
Since the original observations that self-protection against gliotoxin was largely mediated by the enzyme GliT, a gliotoxin oxidoreductase [2, 15], it has subsequently been found that other organisms contain similar enzymes which facilitate self-protection against related epipolythiodioxopiperazines [16, 17]. Indeed, in the bacterium, Streptomyces clavuligerus, it has been demonstrated that an oxidoreductase, HlmI, confers self-protection against the disulfide-bridge-containing, non-ribosomal peptide, holomycin . Interestingly, an RNA methyltransferase, Hom12 in Yersinia ruckeri, also enables self-protection against holomycin, since Hom12 deletion results in acquisition of a holomycin-sensitive phenotype . This clearly infers that self-protection against redox-active non-ribosomal peptides is a multi-faceted process, yet few studies have attempted molecular dissection of the process. Carberry et al.  revealed significantly elevated GSH levels in A. fumigatus ΔgliT, and speculated about exacerbation of gliotoxin toxicity, resulting from formation of the dithiol form of gliotoxin, consequent to this apparent dysregulation in the level of an important cellular reductant. Indeed, the apparent resistance of Saccharomyces cerevisiae Δgsh1, which exhibits significantly attenuated GSH levels, to exogenous gliotoxin supported this hypothesis. Interestingly, these authors also observed that S. cerevisiae Δsod1 and Δyap1 were hypersensitive to exogenous gliotoxin, suggesting that a deficient oxidative stress response sensitizes this organism to gliotoxin. Coleman et al.  further revealed that both Candida albicans and Cryptococcus neoformans were sensitive to gliotoxin exposure, however apart from an elegant demonstration of membrane damage consequent to gliotoxin exposure, no mechanistic basis of the anti-fungal effect of gliotoxin was forthcoming.
It is somewhat surprising that the concept of self-protection against gliotoxin, in fungi capable of its biosynthesis, has received scant attention since the discovery of gliotoxin in 1936 - given the reactive nature of this disulfide-containing metabolite. However the availability of powerful new technologies such as RNA-seq [18–21], provides us with a tool to address this information deficient. Consequently, we present here the first exploration of the global transcriptomic response of both A. fumigatus wild-type and ∆gliT to exogenous gliotoxin, which illuminates the important role played by gliT in mediating control of the cellular systems in the presence of this reactive metabolite.
Deletion of gliT results in altered expression of over 200 genes involved in many functions in A. fumigatus
As deletion of A. fumigatus gliT renders the strain sensitive to exogenous gliotoxin [2, 15], to achieve a better understanding of the self-protection against gliotoxin provided by A. fumigatus gliT, high throughput RNA sequencing analysis was carried out. An average of 9312 transcripts were expressed in A. fumigatus wild-type and ∆gliT (available from the European Nucleotide Archive under accession ERP001382), which is in accordance with other RNA-seq investigations of the A. fumigatus transcriptome [18, 22]. A comparison of A. fumigatus wild-type and ∆gliT revealed that the deletion of gliT resulted in the significant (p < 0.05) dysregulated expression of 208 genes, consisting of 175 up-regulated genes while 33 genes were significantly down-regulated (p < 0.05) (Figure 1A, Additional file 1). The differentially regulated genes were classified according to the Functional Catalogue (FunCat)  and KEGG pathways  to identify the functions of these genes. Following deletion of gliT, 13 2nd level, 19 3rd level and 12 4th level FunCat categories were over-represented in the up-regulated gene set, compared to 3 2nd level, 10 3rd level and 3 4th level FunCat categories enriched in the down-regulated gene set (Additional file 2: Tables S1 and S2). Categories over-represented in the up-regulated gene set include; secondary metabolism, transport and detoxification, while genes involved in C-compound and carbohydrate metabolism were enriched in the down-regulated gene set (Figure 2). Additionally, 11 KEGG pathways were over-represented in the up-regulated gene set, while genes involved in 13 KEGG pathways were enriched in the down-regulated gene set (Additional file 2: Tables S3 and S4). Carbohydrate and lipid metabolism were over-represented in the up-regulated gene set, while amino acid metabolism was enriched in the down-regulated gene set. Genes involved in secondary metabolite biosynthesis were enriched in both the up- and down-regulated gene sets.
Gliotoxin exposure alters the regulation of hundreds of genes in both A. fumigatus wild-type and ∆gliT
Following A. fumigatus wild-type exposure to exogenous gliotoxin, 164 genes were significantly differentially regulated (p <0.05) compared to control conditions (Figure 1B). Of these 164 genes, 101 genes were up-regulated when A. fumigatus wild-type was exposed to gliotoxin, while 63 genes were down-regulated (Additional file 1). In contrast, when A. fumigatus ∆gliT was exposed to exogenous gliotoxin, 1,700 genes were significantly differentially regulated (p <0.05) (Figure 1C), whereby expression of 508 genes was up-regulated in response to gliotoxin exposure and 1192 genes down-regulated (Additional file 1). Given the large transcriptomic remodelling observed in A. fumigatus ∆gliT following exogenous gliotoxin exposure, cell viability was assessed at 85%. Although a significant decrease (p =0.0019) in the cell viability of A. fumigatus ∆gliT exposed to exogenous gliotoxin, compared with the methanol control (97%) (Additional file 2: Figure S1), gross cell death was not observed.
Functional characterisation of differentially regulated genes in A. fumigatus wild-type and ∆gliTfollowing exogenous gliotoxin exposure
In classifying the differentially regulated genes in A. fumigatus wild-type upon exposure to exogenous gliotoxin, 9 2nd level, 17 3rd level and 11 4th level FunCat categories, respectively, were over-represented in the up-regulated gene set compared to 12 2nd level, 17 3rd level and 8 4th level categories in the down-regulated gene set (Additional file 2: Tables S5 and S6). Comparatively, in A. fumigatus ∆gliT exposed to exogenous gliotoxin, a greater number of FunCat categories were over-represented. Here, 24 2nd level, 23 3rd level and 12 4th level FunCat categories, respectively, were over-represented in the up-regulated gene set compared to 42 2nd level, 84 3rd level and 39 4th level categories in the down-regulated gene set (Additional file 2: Tables S7 and S8). A. fumigatus wild-type exposure to exogenous gliotoxin results in an over-representation of up-regulated expression of genes involved in secondary metabolism, transport (particularly siderophore-iron transport), detoxification processes and homeostasis processes as well as others (Figure 3). Of the down-regulated genes, some of the categories which were enriched included amino acid metabolism, carbohydrate metabolism and complex cofactor/cosubstrate/vitamin binding (Figure 3). Loss of gliT in combination with exogenous gliotoxin exposure results in the dysregulation of a large number of processes within the cell. It resulted in an over-representation of up-regulated genes involved in stress response, ribosome biogenesis and translation, and of down-regulated genes involved in metabolism of cysteine, nitrogen, sulphur and selenium metabolism, RNA synthesis, transport (including siderophore-iron transport), homeostasis and cellular import (Figure 4).
In combination with the identified FunCat categories, 2 KEGG pathways were over-represented in the up-regulated gene set, and 13 KEGG pathways in the down-regulated gene set of A. fumigatus wild-type upon exogenous gliotoxin exposure were over-represented (Additional file 2: Tables S9 and S10). Pathways involved in xenobiotic biodegradation and metabolism were enriched in the up-regulated gene set, while in the down-regulated gene set, pathways involved in the biosynthesis of secondary metabolites, glycolysis/gluconeogenesis and amino acid metabolism are over-represented. In A. fumigatus ∆gliT exposed to exogenous gliotoxin, 18 KEGG pathways were over-represented in the up-regulated gene set compared to 38 KEGG pathways in the down-regulated gene set (Additional file 2: Tables S11 and S12). Among the KEGG pathways over-represented in the up-regulated gene set were pathways involved in carbohydrate metabolism, translation and selenoamino acid metabolism. In the down-regulated gene set, some of the pathways over-represented included; biosynthesis of secondary metabolites, starch and sucrose metabolism, and amino acid metabolism.
Exogenous gliotoxin causes dysregulation of gliotoxin biosynthesis cluster gene expression in both A. fumigatus wild-type and ∆gliT
Close inspection of the 13-gene gliotoxin biosynthesis cluster [2, 3], revealed that exogenous gliotoxin caused the dysregulated expression of a number of gliotoxin biosynthetic genes in both A. fumigatus wild-type and ∆gliT (Table 1). Upon exposure to exogenous gliotoxin, five out of the 13 genes in the cluster were significantly up-regulated in A. fumigatus wild-type. Expression of A. fumigatus gliZ was up-regulated log2 3.2-fold, gliP log2 6.2 fold, while gliA and gliF were up-regulated log2 10- and log2 6.7-fold, respectively. Expression of gliT, which confers protection against exogenous gliotoxin [2, 15], was up-regulated log2 9.2-fold. Increased expression of the remaining genes in the gliotoxin biosynthetic cluster was also observed in A. fumigatus wild-type upon exogenous gliotoxin exposure, however altered expression was not significant.
Upon A. fumigatus ∆gliT exposure to exogenous gliotoxin, dysregulated expression of gliotoxin biosynthetic genes was also observed (Table 1). As in A. fumigatus wild-type, gliP, gliA and gliF expression was also up-regulated in ∆gliT upon exogenous gliotoxin addition. A. fumigatus gliP was up-regulated log2 3.4-fold, gliA was up-regulated log2 8.5-fold and A. fumigatus gliF was up-regulated log2 4.7-fold. Additionally, A. fumigatus gliM was up-regulated log2 6.9-fold in A. fumigatus ∆gliT exposed to exogenous gliotoxin, which was not observed in A. fumigatus wild-type upon gliotoxin exposure. Of the remaining gliotoxin biosynthetic genes, while not significant, with the exception of A. fumigatus gliZ and gliI which did not exhibit altered expression, and gliJ which appeared to undergo down-regulated expression, increased expression of the remaining genes was observed in A. fumigatus ∆gliT upon exposure to exogenous gliotoxin. Quantitative real-time PCR (qRT-PCR) analysis of the expression of A. fumigatus gliZ confirmed the up-regulation in A. fumigatus wild-type exposed to exogenous gliotoxin compared to the control, however it also showed increased gliZ expression in ∆gliT in response to gliotoxin (Figure 5), which was not observed in the RNA-seq analysis, possibly due to the sensitivities of the different techniques. qRT-PCR analysis of the expression of A. fumigatus gliA in both A. fumigatus wild-type and A. fumigatus ∆gliT in both the absence and presence of exogenous gliotoxin confirmed the up-regulation in response to exogenous gliotoxin (Figure 5).
The impact of gliotoxin on the expression of the gliotoxin biosynthetic genes is in accordance with that observed in other studies using different techniques [2, 9]. Moreover, feeding experiments with [13C]-phenylalanine herein confirm, for the first time, that gliotoxin biosynthesis is actually induced by exogenous gliotoxin. A significant increase (p =0.0295) was observed in the amount of [13C]-gliotoxin detectable in the wild-type A. fumigatus Af293 culture supernatants following exogenous gliotoxin addition compared to the methanol control (Figure 6).
Exogenous gliotoxin results in down-regulation of secondary metabolite gene cluster expression in A. fumigatus ∆gliT
FunCat classification of the altered gene expression in A. fumigatus wild-type and ∆gliT, respectively, revealed an enrichment of genes involved in secondary metabolism upon exogenous gliotoxin exposure. In A. fumigatus wild-type, expression of 16 genes which was up-regulated, and that of 3 genes down-regulated in response to exogenous gliotoxin, were classified by 2nd level of FunCat as being involved in secondary metabolism (Additional file 2: Tables S5 and S6). A larger number of genes were classified in this category in A. fumigatus ∆gliT in response to exogenous gliotoxin presence, where 22 up-regulated genes and 103 down-regulated genes were observed (Additional file 2: Tables S7 and S8). A closer inspection of some of the secondary metabolite genes clusters revealed significant alterations in the expression of genes in the helvolic biosynthesis cluster  and the “supercluster” on chromosome 8 which encodes the biosynthetic pathways of a number of secondary metabolites, including fumitremorgin B, fumagillin and pseurotin A [26–29] (Table 2).
Helvolic acid, a triterpene, is encoded by a 9-gene cluster on chromosome 4 . In A. fumigatus wild-type, exogenous gliotoxin results in the significant up-regulation of 2 genes, A. fumigatus osc3 (log2 3.28-fold) and sdrI (log2 3.22-fold) (Table 2). Conversely, in A. fumigatus ∆gliT exposed to exogenous gliotoxin, osc3 was significantly down-regulated (log2 6.3-fold) along with a predicted O-methyltransferase, AFUA_4G14580, which is significantly down-regulated (log2 2.7-fold) (Table 2). In addition to this, expression of cyp5081B1 is completely abrogated in A. fumigatus ∆gliT exposed to exogenous gliotoxin (Table 2). qRT-PCR analysis of A. fumigatus osc3 confirmed the observed down-regulation in A. fumigatus ∆gliT and up-regulation of A. fumigatus osc3 in A. fumigatus wild-type exposed to exogenous gliotoxin (Figure 5).
Of the 69 genes in the “supercluster” on chromosome 8, expression of two genes is significantly down-regulated in A. fumigatus wild-type in response to exogenous gliotoxin (Table 2). However, in A. fumigatus ∆gliT, when exposed to exogenous gliotoxin, expression of 26 genes from the “supercluster” was down-regulated (Table 2). Closer inspection of the fumitremorgin B biosynthetic genes revealed that expression of the non-ribosomal peptide synthetase (NRPS), A. fumigatus ftmA , was significantly down-regulated (log2 2.78-fold) in A. fumigatus wild-type in exogenous gliotoxin presence (Table 2). While in A. fumigatus ∆gliT, ftmA (log2 6.65-fold), ftmC (log2 5.68-fold), ftmD (log2 5.87-fold), ftmG (log2 6.52-fold) and ftmI (log2 1.45-fold) expression was significantly down-regulated, and ftmE expression was completely inhibited in response to exogenous gliotoxin (Table 2). qRT-PCR analysis of A. fumigatus ftmA confirmed decreased expression in both wild-type and ∆gliT following exposure to exogenous gliotoxin (Figure 5). Determination of the levels of the fumitremorgins and associated compounds in A. fumigatus wild-type and ∆gliT cultured in secondary metabolite-inducing conditions (96 h in Czpaeks-Dox Broth) revealed significant alterations in the levels of a number of cognate metabolites. Specifically, brevianamide F levels were significantly increased in A. fumigatus ∆gliT compared to wild-type (p =0.0243), while the levels of both tryprostatin A and tryprostatin B were significantly decreased (p =0.008 and p =0.0453) in A. fumigatus ∆gliT compared to wild-type (Figure 7). There was no significant difference determined in the level of fumitremorgin C between A. fumigatus wild-type and ∆gliT (Figure 7).
There was no significant dysregulation of the fumagillin biosynthetic genes in A. fumigatus wild-type in exogenous gliotoxin presence. However, when A. fumigatus ∆gliT was exposed to exogenous gliotoxin, expression of 12 fumagillin biosynthetic genes in the cluster was significantly down-regulated (p <0.05) (Table 2). Some of the down-regulated genes include a polyketide synthase (PKS) (A. fumigatus fma-PKS)  and a putative PKS (AFUA_8G00490)  which were down-regulated log2 7.20- and log2 5.65-fold, respectively, A. fumigatus metAP/fpaII expression was down-regulated log2 5.20-fold and fma-TC/fmaA was down-regulated log2 8.37-fold. The decreased expression of A. fumigatus fma-PKS in A. fumigatus ∆gliT exposed to exogenous gliotoxin was confirmed by qRT-PCR analysis (Figure 5). As was the case for fumitremorgin B and fumagillin, none of the pseurotin A biosynthetic genes were differentially regulated in A. fumigatus wild-type in exogenous gliotoxin presence (Table 2). However, in A. fumigatus ∆gliT exposed to exogenous gliotoxin, with the exception of AFUA_8G00570, expression of all of the pseurotin A biosynthetic genes was significantly down-regulated (Table 2). Expression of AFUA_8G00530 and A. fumigatus psoA/nrps14, the PKS-NRPS hybrid , were down-regulated log2 7.45- and log2 5.79-fold respectively, while A. fumigatus psoC, psoD and elfB  expression was down-regulated log2 7.41-, log2 6.78- and log2 8.21-fold respectively. qRT-PCR confirmed the decreased expression of A. fumigatus psoA/nrps14 in A. fumigatus ∆gliT upon exogenous gliotoxin exposure (Figure 5). Relevantly, under secondary metabolite inducing growth conditions, production of fumagillin and pseurotin A was significantly reduced (p =0.0471 and p =0.0297, respectively) in A. fumigatus ∆gliT compared to wild-type (Figure 7).
The methyltransferase, laeA, is a global regulator of secondary metabolism in A. fumigatus, which fully or partially regulates expression of multiple secondary metabolite gene clusters, including those encoding gliotoxin, helvolic acid, pseurotin A, fumagillin and fumitremorgin biosynthesis [26, 27, 32, 33]. Although dysregulation of secondary metabolite gene expression in A. fumigatus ∆gliT compared to A. fumigatus was observed after exposure to gliotoxin, laeA expression was not differentially regulated in A. fumigatus wild-type in exogenous gliotoxin presence, but was significantly down-regulated (p =0.015; log21.31-fold) in A. fumigatus ∆gliT under identical conditions. Although significantly down-regulated, it was outside the cut-off of log21.5-fold change. qRT-PCR analysis for A. fumigatus laeA did not show any differential regulation in either A. fumigatus wild-type or A. fumigatus ∆gliT when challenged with exogenous gliotoxin (Figure 5). This suggests that laeA expression is not solely responsible for the altered expression of secondary metabolite genes, particularly in A. fumigatus ∆gliT when it is exposed to exogenous gliotoxin. Expression of A. fumigatus gliT in ∆laeA  exposed to exogenous gliotoxin for 3 h was analysed. A. fumigatus gliT expression increased following exogenous gliotoxin exposure in ∆laeA (Figure 8) indicating that loss of laeA does not affect gliT expression in the presence of exogenous gliotoxin. Recently, the global regulatory velvet gene, A. fumigatus veA, has been implicated in secondary metabolite biosynthesis regulation, particularly with respect to gliotoxin  and fumagillin . However, from the data presented here, A. fumigatus veA is not differentially regulated in either A. fumigatus wild-type or A. fumigatus ∆gliT in exogenous gliotoxin presence. A. fumigatus gliT expression in A. fumigatus ∆veA  was assessed following exposure to exogenous gliotoxin and was increased in ∆veA in the presence of exogenous gliotoxin (Figure 8). A. fumigatus gliT expression was higher in the methanol control (Figure 8) compared to the other mutants under the same conditions. Expression of A. fumigatus gliT in ∆veA grown in media only was lower than the methanol control (data not shown), suggesting that methanol induced increased gliT expression in A. fumigatus ∆veA. The increased gliT expression in A. fumigatus ∆veA following methanol addition was only observed in this strain, and was not observed in other deletion strains generated in this background (e.g. ∆laeA).
The expression of A. fumigatus gliT was also investigated in deletion mutants with abrogated fumitremorgin B, fumagillin and pseurotin A production namely, A. fumigatus ∆fapR, ∆fmaA and ∆psoA, respectively . As was the case in A. fumigatus laeA, expression of A. fumigatus gliT increased in these mutants following exogenous gliotoxin exposure (Figure 8) indicating that loss of production of these metabolites does not affect gliT expression.
Siderophore-iron transport and siderophore biosynthesis is affected by exogenous gliotoxin and to a greater extent in A. fumigatus ∆gliT
Functional classification of the differentially regulated gene set revealed an enrichment of genes involved in siderophore-iron transport in both A. fumigatus wild-type and A. fumigatus ∆gliT when exposed to exogenous gliotoxin. In A. fumigatus wild-type exposed to exogenous gliotoxin, expression of six siderophore-iron transport genes was up-regulated, while in A. fumigatus ∆gliT, expression of two siderophore-iron transport genes was up-regulated and 12 genes are down-regulated (Table 3). A. fumigatus sidF was de novo expressed in A. fumigatus wild-type upon exposure to exogenous gliotoxin. Expression of A. fumigatus mirD was up-regulated in both A. fumigatus wild-type (log2 2.57-fold) and A. fumigatus ∆gliT (log2 3.07-fold) in exogenous gliotoxin presence. Interestingly, three genes that were significantly up-regulated in A. fumigatus wild-type exposed to exogenous gliotoxin were significantly down-regulated in A. fumigatus ∆gliT, when challenged with exogenous gliotoxin. A. fumigatus fre7 was significantly up-regulated log2 2.80-fold in A. fumigatus wild-type, but significantly down-regulated log2 2.77-fold in A. fumigatus ∆gliT. AFUA_3G13670 and AFUA_3G13700, significantly up-regulated log2 7.03- and log2 7.62-fold respectively, in A. fumigatus wild-type in exogenous gliotoxin presence, were significantly down-regulated upon exogenous gliotoxin exposure log2 7.06- and log2 7.82-fold in A. fumigatus ∆gliT. Expression of AFUA_3G13700 in both A. fumigatus wild-type and ∆gliT exposed to exogenous gliotoxin, was confirmed by qRT-PCR (Figure 5).
The observed differential regulation of genes involved in siderophore-iron transport prompted us to review whether expression of siderophore biosynthetic genes  was affected by exogenous gliotoxin in A. fumigatus wild-type and A. fumigatus ∆gliT. Indeed we have already noted that sidF expression, a transacylase that transfers anhydromevalonyl to hydroxyornithine during extracellular siderophore biosynthesis , was activated in A. fumigatus wild-type following exogenous gliotoxin exposure. In A. fumigatus ∆gliT, expression of two genes involved in siderophore biosynthesis was differentially regulated in response to exogenous gliotoxin (Table 3). A. fumigatus sidH, involved in providing the anhydromevalonyl-CoA moiety for extracellular siderophore biosynthesis , was de novo expressed while A. fumigatus sidA, which is required for the first step of siderophore biosynthesis , was significantly up-regulated log2 2.16-fold (p =0.0005). qRT-PCR analysis of the expression of A. fumigatus sidH confirmed the increased expression in A. fumigatus ∆gliT following exogenous gliotoxin exposure (Figure 5).
Nitrogen metabolism is down-regulated in A. fumigatus ∆gliTin response to exogenous gliotoxin
In A. fumigatus ∆gliT exposed to exogenous gliotoxin, expression of 55 genes involved in nitrogen metabolism was down-regulated, while expression of 4 genes was up-regulated (Table 4). In contrast, 3 genes were up-regulated in A. fumigatus wild-type and 8 genes were down-regulated in response to exogenous gliotoxin (Table 4). In A. fumigatus ∆gliT, expression of two genes, cyp5081B1, a putative cytochrome P450 monooxygenase already mentioned in helvolic acid biosynthesis , and AFUA_6G00412, predicted to have amino acid transmembrane transporter activity  was completely abrogated in response to exogenous gliotoxin. Expression of A. fumigatus optB and cpsI, which are induced when BSA is the sole nitrogen source , was down-regulated log2 5.399- and log2 3.98-fold, respectively, in A. fumigatus ∆gliT exposed to exogenous gliotoxin. The decreased expression of A. fumigatus optB in A. fumigatus ∆gliT exposed to exogenous gliotoxin was confirmed by qRT-PCR (Figure 5). When utilising complex nitrogen sources, A. fumigatus secretes proteases, and increases expression of amidase, aminotransferase, and amino acid and peptide transporter genes [40, 41]. These secreted proteases are regulated by a conserved regulatory factor, prtT [40, 41], which is decreased log2 3.66-fold in expression in ∆gliT but unchanged in wild-type following exogenous gliotoxin exposure. Indeed, further examination revealed decreased expression of 14 proteases, 3 amidases and 7 amino acid and peptide transporters in ∆gliT exposed to exogenous gliotoxin, while expression of one amidase, AFUA_5G09140, was increased log2 3.37-fold (Table 4). Expression of both nitrate transporters, crnA and nrtB, which are predicted gene pairs , was significantly down-regulated (log2 2.689- and log2 4.504-fold) in A. fumigatus ∆gliT upon exogenous gliotoxin exposure. Expression of a number of genes with predicted involvement in the oxidation-reduction process was down-regulated in A. fumigatus ∆gliT only following exogenous gliotoxin exposure. Among these genes were nitronate monooxygenases (AFUA_2G17430 and AFUA_4G07940 down-regulated log2 4.444- and log2 2.304-fold) and copper ion binding domain genes (AFUA_1G13440, AFUA_5G01470 and AFUA_5G07360; down-regulated log21.692-, log2 4.125- and log2 3.948-fold, respectively). In A. fumigatus wild-type exposed to exogenous gliotoxin, two predicted carbon-nitrogen ligases, AFUA_4G04160 and AFUA_4G04170, were down-regulated log2 1.926- and log2 2.622-fold, respectively. These carbon-nitrogen ligases are not differentially regulated in A. fumigatus ∆gliT upon exogenous gliotoxin exposure, suggesting that loss of gliT may hinder the down-regulation of these genes in response to gliotoxin.
Exogenous gliotoxin exposure alters gene expression involved in ribosome biogenesis and translation in A. fumigatus ∆gliTwhile transcription related genes are down-regulated
Dysregulation of ribosome biogenesis and translation was observed in A. fumigatus ∆gliT following exogenous gliotoxin exposure, whereby expression of 20 genes was up-regulated and two genes down-regulated (Table 5). Of the 20 genes up-regulated in A. fumigatus ∆gliT in response to exogenous gliotoxin, 7 genes encode 60S ribosomal proteins and two encode 40S ribosomal proteins, all of which are increased log21.5-1.7-fold in expression (Table 5). The two genes down-regulated in response to exogenous gliotoxin in A. fumigatus ∆gliT have already been mentioned as secondary metabolite biosynthetic genes. A. fumigatus metAP in the fumagillin biosynthetic gene cluster  was down-regulated log2 4.340-fold and A. fumigatus elfB in the pseurotin A biosynthetic gene cluster was down-regulated log2 7.265-fold [27, 31]. In contrast to A. fumigatus ∆gliT, only two genes involved in ribosome biogenesis and translation were up-regulated in A. fumigatus wild-type exposed to exogenous gliotoxin (Table 5). AFUA_5G07340 was up-regulated log2 1.909-fold while AFUA_2G12150 was up-regulated log2 2.417-fold. In addition to the dysregulation of ribosome biogenesis and translation, 44 genes involved in transcription were also down-regulated in A. fumigatus ∆gliT in response to exogenous gliotoxin, compared to 3 down-regulated genes in wild-type exposed to exogenous gliotoxin (Table 6). A proteomic investigation of A. fumigatus ∆gliTATCC26933 reflected this alteration in translation whereby 6 proteins were uniquely expressed or increased in abundance and 15 proteins were absent or decreased in abundance when exposed to exogenous gliotoxin (Tables 7 and 8). Loss of GliT negatively impacts expression of genes involved in transcription processes and results in significantly decreased transcription of genes in the presence of exogenous gliotoxin, culminating in disruption of translation.
To dissect the role played by the gliotoxin oxidoreductase, gliT, in self-protection against gliotoxin, a comparative transcriptomic analysis of the impact of exogenous gliotoxin on A. fumigatus wild-type and ∆gliT was carried out via RNA-seq analysis. To reduce the affects of endogenous gliotoxin, a low gliotoxin-producing background strain, A. fumigatus ATCC46645  was employed, along with Sabouraud-Dextrose media which is non-permissive for gliotoxin biosynthesis . Exogenous gliotoxin exposure resulted in changed expression of 164 genes in A. fumigatus wild-type. However, altered expression of over 1700 genes was observed in A. fumigatus ∆gliT. Closer inspection revealed alterations in expression of clusters encoding secondary metabolites, particularly gliotoxin, helvolic acid, fumitremorgin, fumagillin and pseurotin A biosynthesis genes, siderophore biosynthesis genes, ribosome biogenesis genes and genes involved in translation and nitrogen metabolism.
Exogenous gliotoxin induces the expression of gliotoxin biosynthetic genes in both A. fumigatus wild-type and ∆gliT. Expression of gliZ was increased in expression in wild-type, along with gliP, the bimodular non-ribosomal peptide synthetase that produces the cyclo-L-Phe-L-Ser diketopiperazine intermediate [9, 10, 43]. Up-regulation of gliZ and gliP expression in particular, along with gliF suggests that exogenous gliotoxin induces de novo gliotoxin production in wild-type. Indeed this has been shown to be the case in A. fumigatus Af293 where addition of exogenous gliotoxin under gliotoxin-inducing culture conditions resulted in detection of de novo gliotoxin production determined by the use of [13C]-Phenylalanine and LC-MS analysis. Recently, cyclo-L-Phe-L-Ser was determined to be a major component of the metabolome and was detected in lung tissue of infected mice [11, 44]. The production of this compound may be increased consequent to exogenous gliotoxin exposure in wild-type due to up-regulation of gliP expression. It has been demonstrated that A. fumigatus gliT expression is not solely under the control of gliZ , and so up-regulation of gliZ expression is not absolutely required to induce gliT expression for the purpose of self-protection. While gliZ expression is not significantly altered in ∆gliT following exogenous gliotoxin exposure in the RNA-seq analysis, qRT-PCR demonstrated increased expression under these conditions. Expression of a number of gli genes was up-regulated in ∆gliT following exogenous gliotoxin exposure, including gliM, which was not altered in wild-type. A. fumigatus gliM is a predicted O-methyltransferase  and methylation has been proposed as a method of self-protection against disulfide bridge-containing metabolites. Holomycin contains a disulfide bridge similar to gliotoxin [16, 17, 45] and S-methylation has been proposed as an alternative method of self-protection against this antibiotic and its biosynthetic intermediates in S. clavuligerus upon deletion of the dithiol oxidase, HlmI, which is functionally homologous to A. fumigatus gliT [17, 45]. In Y. ruckeri, an RNA methyltransferase, Hom12, methylates RNA in a proposed self-protection mechanism against the cytotoxic effects of holomycin . Interestingly, in A. fumigatus ∆gliT following exogenous gliotoxin exposure, a tRNA methyltransferase, AFUA_4G12280, is up-regulated log2 2.56-fold, while its expression is unaltered in wild-type in the presence of exogenous gliotoxin. The increased expression of these methyltransferases, among others, in A. fumigatus ∆gliT following exposure to exogenous gliotoxin may suggest possible alternative functionalities, or self-protection mechanisms against gliotoxin in the absence of GliT. Interestingly, a newly-identified methyltransferase, gliotoxin methyltransferase A (GtmA), has been demonstrated to play a role in attenuating gliotoxin biosynthesis . GtmA (AFUA_2G11120) expression is significantly up-regulated (p <0.00005) by gliotoxin exposure (Additional file 1).
In addition to up-regulation of the gli cluster, dysregulation of other secondary metabolite gene clusters was observed upon exogenous gliotoxin addition, particularly to A. fumigatus ∆gliT. Up-regulated expression of two genes; osc3 and sdrI, from the gene cluster encoding the fusidane antibiotic helvolic acid was observed upon exogenous gliotoxin addition to A. fumigatus wild-type . However, exogenous gliotoxin addition to A. fumigatus ∆gliT resulted in down-regulated expression of two genes and abrogation of the expression of cyp5081D1, suggesting that gliT deletion, along with exogenous gliotoxin exposure negatively regulates the helvolic acid gene cluster. In an A. fumigatus gliZ complemented strain, helvolic acid production was detectable at 37°C, whereas no helvolic acid was detectable in wild-type or ∆gliZ strains at an identical growth temperature , suggesting that gliZ may be involved in regulating helvolic acid production. Indeed in the present study in wild-type, gliZ expression was up-regulated upon gliotoxin exposure and the helvolic acid biosynthetic genes, including osc3 which encodes the protein that catalyses the first biosynthetic step of helvolic acid biosynthesis , are also up-regulated, whereas in ∆gliT, the helvolic acid biosynthesis genes are down-regulated. However, we cannot unambiguously assign gliZ functionality to hevolic acid production because consistently discrepant gliZ expression data was obtained by RNA-seq and qRT-PCR (Figure 5).
The combined loss of gliT and exposure to exogenous gliotoxin resulted in attenuated expression of the biosynthetic genes for the secondary metabolites: fumitremorgin B, fumagillin, and pseurotin A. Secondary metabolite production in fungi is generally dependant on the presence of a key backbone enzyme, namely a PKS or a NRPS . In A. fumigatus ∆gliT exposed to exogenous gliotoxin, the NRPS ftmA in the fumitremorgin B cluster, the PKS, fma-PKS in the fumagillin cluster and the PKS-NRPS hybrid psoA/nrps14 in the pseurotin A cluster were all significantly down-regulated suggesting production of these secondary metabolites is also decreased. A. fumigatus ftmA encodes the enzyme required for the first biosynthetic step of fumitremorgin synthesis, the synthesis of the diketopiperazine, brevianamide F , and is significantly down-regulated log2 6.640-fold in A. fumigatus ∆gliT following exogenous gliotoxin addition. Expression of the cytochrome P450 A. fumigatus ftmE, which was absent in ∆gliT exposed to exogenous gliotoxin, encodes the enzyme responsible for formation of fumitremorgin C , while another cytochrome P450, ftmG, which encodes the enzyme that dihydroxylates the fumitremorgin B intermediate was significantly down-regulated log2 6.519-fold . We hypothesise that fumitremorgin B biosynthesis is decreased, if not abrogated, in A. fumigatus ∆gliT when exposed to exogenous gliotoxin given the decreased expression of genes essential for synthesis. While determination of fumitremorgin B levels was not successful, measurement of fumitremorgin C and related compounds  was carried out in A. fumigatus wild-type and ∆gliT cultured under secondary metabolite-inducing conditions. Significant reductions in the levels of tryprostatin A and tryprostatin B were observed in A. fumigatus ∆gliT compared to wild-type, while there was no significant difference in levels of fumitremorgin C between the two strains. Therefore, GliT loss disrupts the production of brevianamide F, tryprostatin A and tryprostatin B, from the fumitremorgin biosynthetic pathway in A. fumigatus. The PKS, fma-PKS, down-regulated log2 7.287-fold in A. fumigatus ∆gliT, is essential for fumagillin biosynthesis . Additionally, the C6 type transcription factor fapR/fumR, which controls the expression of the other fumagillin cluster genes , was down-regulated log2 5.826-fold in A. fumigatus ∆gliT following exogenous gliotoxin addition. Decreased expression of these two genes, in addition to others in the biosynthetic gene cluster, suggests down-regulation of fumagillin biosynthesis as a direct consequence of A. fumigatus gliT loss combined with exogenous gliotoxin stress. Indeed supporting the hypothesis of A. fumigatus gliT involvement in facilitating fumagillin biosynthesis, under secondary metabolite-inducing conditions we observed a significant decrease in the production of fumagillin in A. fumigatus ∆gliT compared to wild-type. It has been determined that the biosynthetic genes for fumagillin and pseurotin A are physically intertwined  and these authors also revealed that fapR, which controls expression of the fumagillin biosynthesis genes, also controls the expression of pseurotin A biosynthesis genes . Indeed, with the exception of one gene, expression of the pseurotin A biosynthetic cluster was significantly down-regulated in A. fumigatus ∆gliT upon exogenous gliotoxin addition. psoA/nrps14 is essential for pseurotin A biosynthesis, as demonstrated by Maiya et al. , while over-expression of this gene increased pseurotin A accumulation. In A. fumigatus ∆gliT exposed to exogenous gliotoxin, psoA/nrps14 expression was down-regulated log2 5.826-fold, while in the proteomic investigation, PsoA/nrps14 was decreased log2 3.032-fold in abundance. In addition to PsoA/nrps14, two other proteins required for pseurotin A biosynthesis  were decreased in abundance upon exogenous gliotoxin exposure in A. fumigatus ∆gliT. PsoF (AFUA_8G00440) and PsoC (AFUA_8G00550) were decreased log2 2.461- and log2 2.016-fold providing further support to the hypothesis that pseurotin A biosynthesis is down-regulated in the absence of A. fumigatus gliT when challenged with exogenous gliotoxin. As was the case with fumagillin, pseurotin A levels were significantly decreased in A. fumigatus ∆gliT compared to wild-type when cultured under secondary metabolite-inducing conditions again suggesting possible GliT involvement in enabling pseurotin A biosynthesis.
The biosynthetic gene clusters of fumitremorgin, fumagillin and pseurotin are under the control of the global regulator, A. fumigatus laeA [26, 27]. A. fumigatus laeA is a methyltransferase and regulates chromatin remodelling through this methyltransferase activity [26, 32, 49]. In A. fumigatus ∆gliT treated with exogenous gliotoxin, laeA expression was not significantly altered, suggesting that it is not responsible for the down-regulation of these gene clusters. In A. nidulans, LaeA forms a trimeric complex with two members of the velvet protein family; VeA and VelB, and this complex up-regulates secondary metabolism and sexual development . The complex was subsequently identified and characterised in A. fumigatus . A. fumigatus veA has been demonstrated to regulate fumagillin and gliotoxin production [20, 22], however veA expression was not significantly altered in A. fumigatus ∆gliT following exogenous gliotoxin production. The bZip transcriptional enhancer RsmA has been shown to positively regulate gliotoxin biosynthesis , however the expression of A. fumigatus rsmA was not altered significantly in A. fumigatus ∆gliT, in the absence or presence of exogenous gliotoxin. Overall, our findings lead us to postulate that the gliotoxin biosynthetic/self-protection capacity (i.e., GliT functionality) is necessary for optimal biosynthesis of selected secondary metabolites in A. fumigatus.
Altered expression of siderophore-iron transport and siderophore biosynthetic genes in both A. fumigatus wild-type and ∆gliT following exogenous gliotoxin exposure suggests a disruption in iron homeostasis or iron sensing. Iron is an essential nutrient that is required for many cellular processes, including as a cofactor for numerous enzymes . In both wild-type and ∆gliT, siderophore biosynthesis gene expression was up-regulated in response to exogenous gliotoxin indicating that consequent to exogenous gliotoxin exposure, there is an increased requirement for iron despite sufficient iron availability in the culture media, or else a deficit in iron-sensing. A. fumigatus sidF expression was activated in wild-type, while in ∆gliT, sidH expression was activated and that of sidA was up-regulated . Interestingly, up-regulation of sidA expression has also been observed in ∆metR, an A. fumigatus mutant deficient in the transcription factor that regulates sulfur assimilation, under iron sufficient but sulfur deficient conditions . Regulatory cross-talk between secondary metabolism and iron requirement has been reported whereby in a laeA mutant deficient in gliotoxin production, decreased expression of a number of the siderophore biosynthesis genes was observed under high iron conditions . These authors concluded that laeA was also involved in regulating expression of the siderophore biosynthetic genes and in particular, sidD. As discussed, we have observed that the combined effect of exogenous gliotoxin exposure and gliT deletion has a significant impact on secondary metabolism, despite laeA expression being unaffected, and therefore it is interesting that siderophore biosynthesis is also affected. In A. fumigatus wild-type, expression of six siderophore-iron transport genes was up-regulated in response to exogenous gliotoxin, in contrast to up-regulation of two siderophore-iron transport genes in ∆gliT. Interestingly, in ∆gliT, expression of 12 siderophore-iron transport genes was down-regulated. Taken together, this suggests that deletion of gliT in combination with exogenous gliotoxin exposure results in a disruption of, or decrease in, siderophore-iron transport. Amich et al.  noted up-regulation of siderophore-iron transport genes in ∆metR under sulfur-limited but iron replete conditions. Despite observing increased expression of both siderophore biosynthetic genes and siderophore-iron transport genes in ∆metR, suggesting iron starvation, the authors noted increased levels of ferricrocin, the intracellular siderophore that is used for transport and storage [38, 51]. Iron is utilised in many processes one of which is iron-sulfur cluster biosynthesis. Iron-sulfur clusters are inorganic cofactors involved in cellular processes including enzyme activity regulation, mitochondrial respiration, ribosome biosynthesis and cofactor biosynthesis . Translocation of iron-sulfur clusters requires glutathione, and depletion has been shown to induce an iron starvation-like response in Saccharomyces cerevisiae [53–55]. It has previously been shown that exogenous gliotoxin exposure to both A. fumigatus wild-type and ∆gliT results in decreased GSH levels  and so it is conceivable that this decrease in cellular GSH could impact on iron homeostasis or iron sensing.
A. fumigatus encodes two transcription factors, hapX and sreA, that maintain iron homeostasis whereby hapX represses sreA expression and subsequently iron-consuming pathways, and activates siderophore biosynthesis during iron-starvation, while sreA represses hapX during iron-sufficient conditions in a negative feedback loop [35, 56]. Interestingly, induction of gliT has been shown before when cultures were shifted from iron-limited to replete conditions in both wild-type and a sreA deletion strain . Although expression of hapX and sreA is unchanged in both wild-type and ∆gliT upon exogenous gliotoxin exposure, there is significant interplay between gliT, sulfur and iron as demonstrated by the altered expression of siderophore biosynthesis and siderophore-iron transport genes and the decreased GSH levels  in A. fumigatus ∆gliT.
FunCat analysis identified significant enrichment of genes involved in nitrogen metabolism in the down-regulated gene set in ∆gliT exposed to exogenous gliotoxin. Fungi can utilise various sources of nitrogen, from easily assimilated sources (e.g., ammonium and glutamate), to more complex secondary sources including amino acids and proteins [40, 58]. Following exposure to exogenous gliotoxin, the expression of prtT, the conserved regulator of secreted proteases, was decreased and consequentially, a number of proteases previously identified to be under its control were also decreased in expression, including alp1 and pep1 . Additionally, a number of other proteases underwent decreased expression in ∆gliT when exposed to exogenous gliotoxin. These proteases, some of which include dppV, dppVI, cpdS and gprA, have been shown to be induced when BSA is the sole nitrogen source . In order to conserve energy, fungi will preferentially utilise nitrogen sources that are easily assimilated over complex nitrogen sources . Sabouraud-Dextrose medium contains a pancreatic digest of casein and a peptic digest of animal tissue. As the media already contained digested proteins, we conclude that expression of these secreted proteases is significantly down-regulated in ∆gliT exposed to exogenous gliotoxin in order to conserve energy. Interestingly, a number of genes involved in amino acid and peptide transport are also decreased in expression with the exception of one amidase, AFUA_5G09140, the expression of which is up-regulated in ∆gliT following exogenous gliotoxin exposure. This was surprising as, despite the availability of amino acids and peptides in the media, the expression of genes encoding enzymes required for uptake of these amino acids and peptides are down-regulated. It is possible that sufficient enzymes are present for adequate uptake of nitrogen sources to compensate for the decreased expression of these genes, or that the down-regulation of these genes is a consequence of the down-regulation of prtT and the other secreted proteases suggesting that there is a decreased nitrogen requirement in ∆gliT, when exposed to exogenous gliotoxin.
Loss of A. fumigatus gliT had a significant impact on the transcriptome when challenged with exogenous gliotoxin whereby expression of 1,700 genes was altered, of which 1192 were down-regulated, that was not observed in wild-type where 164 genes had altered expression. This suggests that consequent to GliT absence, transcription is suppressed in the presence of exogenous gliotoxin. This is further supported by the decreased expression of 44 genes involved in transcription processes, many of which are transcription factors. Further to the effects of gliT deletion on transcription, translation is also disrupted in A. fumigatus ∆gliT following exposure to exogenous gliotoxin. Altered expression of genes required for ribosome biogenesis and translation processes was observed in ∆gliT. To support the hypothesis that translation is disrupted, the LFQ proteomic investigation of A. fumigatus ∆gliT did not reflect the large transcriptome changes, as only 6 proteins with increased abundance or uniquely present, and 15 proteins absent or with decreased abundance were identified. It is important to note that despite the low number of proteins altered in abundance, there is agreement between the RNA-seq analysis and the proteomic analysis of A. fumigatus ∆gliT exposed to exogenous gliotoxin, whereby the genes encoding one of the proteins increased in abundance and 8 of the proteins decreased in abundance were up-regulated and down-regulated, respectively.
We present the first global investigation of the transcriptional response to exogenous gliotoxin in A. fumigatus wild-type and ∆gliT employing RNA-seq analysis. While exogenous gliotoxin elicits some transcriptome remodelling in wild-type, in ∆gliT, the transcriptional response is 10-fold that of wild-type, with approximately 70% of these altered genes decreased in expression. We found that the combined loss of gliT and exogenous gliotoxin exposure results in decreased expression of a number of secondary metabolite genes from the biosynthetic clusters of helvolic acid, fumitremorgin B, pseurotin A and fumagillin despite the unchanged expression of laeA, the global regulator of secondary metabolism. Thus, GliT functionality may extend to enhancing the biosynthesis of selected secondary metabolites in A. fumigatus. In addition to this, the decreased expression of many transcription factor genes, along with genes involved in siderophore-iron transport and siderophore biosynthesis and nitrogen metabolism indicates that exogenous gliotoxin induces a starvation-like response despite the use of a rich media. Furthermore, the combined RNA-seq and proteomic analysis revealed deletion of gliT abrogated transcription and disrupts translation when exposed to exogenous gliotoxin. Overall, this study provides a detailed overview of the response to exogenous gliotoxin in resistant and sensitive A. fumigatus strains, enhances our understanding of the manner in which gliotoxin exerts its affects as a toxin, and provides a unique glimpse into cross-talk between apparently unrelated secondary metabolite gene clusters.
Aspergillus fumigatusstrain information and growth conditions
Conidia were maintained on malt extract agar plates. A. fumigatus wild-type and mutant strains (1 × 106 cfu/ml) were cultured in Sabouraud-Dextrose media at 37°C with shaking at 200 rpm for 21 h in duplicate or triplicate. Gliotoxin (5 μg/ml final) or MeOH (solvent control) was added and the cultures were incubated for a further 3 h. The mycelia were harvested through miracloth, washed with DEPC water and snap frozen in liquid N2.
RNA extraction and mRNA isolation
RNA was isolated from mycelia, ground to a fine powder in liquid N2, using the RNeasy™ Plant Mini Kit (Qiagen), according to the manufacturer’s instructions. RNA integrity was analysed using an Agilent 2100 Bioanalyzer™ and an Agilent RNA 6000 Nano Kit following the manufacturer’s recommendation.
Library preparation and sequencing
A library was independently prepared for each biological replicate. Two protocols (TruSeq and Illumina mRNA-seq kit) were used for preparing the Illumina transcriptome libraries. For both protocols, polyadenylated mRNA was purified from total RNA using oligo-dT dynabead selection followed by metal ion hydrolysis fragmentation with an RNA fragmentation solution supplied in kits. First strand synthesis, primed using random oligonucleotides, was followed by second strand synthesis with RNaseH and DNA polI to produce double-stranded cDNA using the Illumina mRNA Seq kit or the TruSeq Illumina kit. Template DNA fragments were end-repaired with T4 and Klenow DNA polymerases and blunt-ended with T4 polynucleotide kinase. A single 3′ adenosine was added to the repaired ends using Klenow exo- and dATP to reduce template concatemerization and adapter dimer formation, and to increase the efficiency of adapter ligation. Adapters (containing primer sites for sequencing, and index sequences when using the TruSeq protocol) were then ligated. Libraries made with the TruSeq protocol were amplified by PCR using KAPA HiFi Polymerase (to enrich for properly ligated template strands, to generate enough DNA, and to add primers for flowcell surface annealing). AMPure SPRI beads were used to purify amplified templates before pooling based on quantification using an Agilent Bioanalyser chip. Pooled TruSeq libraries were then pooled and size selected using the Caliper. After adaptor ligation, individual libraries made with the Illumina mRNA-seq kit were size selected using the CaliperLabChip before PCR amplification followed by AMPure SPRI bead clean up and removal of adaptors with a second Caliper run. KAPA Illumina SYBR Fast qPCR kit was used to quantify the Illumina mRNA-seq libraries before pooling. No qPCR was necessary with the TruSeq libraries and instead a final Agilent Bioanalyser chip analysis was run to confirm the dilution of the final pool. The libraries were sequenced on the Illumina HiSeq platform with a read length of 75 bp paired-ended according to manufacturer’s instructions.
The RNA-seq paired-end reads for each biological replicate were aligned independently using Tophat v2.0.4  to the A. fumigatus Af293 (CADRE 3a) reference genome sequence with default parameters. The numbers of fragments mapped per replicate is given in Additional file 2: Table S13. Differential gene expression analysis was carried out for each sample independently using Cufflinks (cuffdiff) v2.0.2  with default parameters against the gene set in Ensembl Genomes release 14 (CADRE genebuild 3a). The sequencing data has been submitted to the European Nucleotide Archive (ENA) under accession ERP001382 (https://www.ebi.ac.uk/ena/data/view/ERP001382&display=html).
The differentially regulated genes were analysed using FungiFun  to establish association with any functions or pathways in comparison with the non-differentially regulated genes. Gene enrichment analysis was carried out on the up-regulated and down-regulated genes, respectively, for FunCat  categories and the A. fumigatus annotated KEGG pathways . A p value cut-off of 0.05 was used.
Real time PCR (qRT-PCR)
RNA samples were DNase treated using a DNase kit supplied by Sigma-Aldrich. cDNA synthesis was performed using qScript™ cDNASuperMix (Quanta Biosciences) following the kit instructions. Primers used in this study are listed in Additional file 2: Table S14. The constitutively expressed gene, A. fumigatus calmodulin (calm)  was used as a reference gene. qRT-PCR was performed on the LightCycler®2 480 Real-Time PCR System using the LightCycler®1 Sybr Green 1 Master Mix (Roche) as described previously . qRT-PCR reactions for each gene were analysed in triplicate and were carried out for 40 cycles.
Feeding experiments with [13C]-phenylalanine
A. fumigatus Af293 was cultured in Czapek-Dox Broth at 37°C with shaking at 200 rpm in duplicate for 24 h before addition of gliotoxin (5 μg/ml final) or methanol (solvent control) and [13C]-phenylalanine (10 μg/ml final) or water (solvent control). Cultures were incubated again at 37°C with shaking at 200 rpm. Culture supernatant was removed after 48 h and again after 72 h. Supernatants were chloroform extracted, dried and resolubilised in methanol prior to LC-MS analysis as described previously .
Whole proteome analysis
A. fumigatus ∆gliTATCC26933 was cultured in Sabouraud-Dextrose media for 21 h followed by gliotoxin (5 μg/ml final) or methanol addition for 3 h (n =4 biological replicates for all specimens). Mycelial lysates were prepared in lysis buffer (100 mM Tris–HCl, 50 mM NaCl, 20 mM EDTA, 10% (v/v) Glycerol, 1 mM PMSF, 1 μg/ml pepstatin A, pH 7.5) with grinding, sonication and clarified using centrifugation. The resultant protein lysates were precipitated using TCA/acetone and resuspended in 8 M Urea. After DTT reduction and IAA-mediated alkylation , sequencing grade trypsin combined with ProteaseMax surfactant was added. Digested samples were desalted prior to analysis using C18 spin columns (Thermo Scientific). All peptide mixtures were analysed via a Thermo Scientific Q-Exactive mass spectrometer coupled to a Dionex RSLCnano. LC gradients ran from 14–35%B (A: 0.1% (v/v) formic acid, B: 80% (v/v) acetonitrile, 0.1% (v/v) formic acid) over 2 h, and data was collected using a Top15 method for MS/MS scans. Comparative proteome abundance and data analysis was performed using MaxQuant software (Version 220.127.116.11) , with Andromeda used for database searching and Perseus (Version 18.104.22.168) used to organise the data. Carbamidomethylation of cysteines was set as a fixed modification, while oxidation of methionines and acetylation of N-termini were set as variable modifications. The maximum peptide/protein false discovery rates (FDR) were set to 1% based on comparison to a reverse database. The Label-Free Quantification (LFQ) algorithm was used to generate normalised spectral intensities and infer relative protein abundance. Proteins that matched to a contaminants database or the reverse database were removed and proteins were only retained in final analysis if detected in at least three replicates from at least one sample. Quantitative analysis was performed using a t-test to compare pairs of samples, and proteins with significant change in abundance (p value <0.05; fold change ≥2) were included in the quantitative results. Qualitative analysis was also performed, to detect proteins that were found in at least 3 replicates of a particular sample, but undetectable in the comparitor sample.
Secondary metabolite analysis
A. fumigatus wild-type and ∆gliT46645 were cultured in Czapeks-Dox broth for 96 h at 37°C at 200 rpm in duplicate. Culture supernatants were organically extracted using an equal volume of chloroform and the extracts were dried by rota-evaporation and resuspended in methanol. Extracts were diluted 1/10 in 0.01% formic acid prior to analysis by LC-MS as described previously .
Cell viability analysis
Following culturing and treatment as described above, mycelia were removed from the culture supernatant and washed with water. They were resuspended in 10 mM HEPES before aliquots were plated on Sabouraud dextrose agar and incubated at 37°C for 24 h. Viability was recorded.
Quantitative real-time PCR
Non-ribosomal peptide synthetase
Bernardo PH, Brasch N, Chai CLL, Waring P: A novel redox mechanism for the glutathione-dependent reversible uptake of a fungal toxin in cells. J Biol Chem. 2003, 278: 46549-46555.
Schrettl M, Carberry S, Kavanagh K, Haas H, Jones GW, Brien O, Nolan A, Stephens J, Fenelon O, Doyle S: Self-protection against Gliotoxin — a component of the Gliotoxin Biosynthetic Cluster, GliT, completely Protects Aspergillus fumigatus against Exogenous Gliotoxin. PLoS Pathog. 2010, 6: e1000952-doi:10.1371/journal.ppat.1000952
Gardiner DM, Howlett BJ: Bioinformatic and expression analysis of the putative gliotoxin biosynthetic gene cluster of Aspergillus fumigatus. FEMS Microbiol Lett. 2005, 248: 241-248.
Coleman JJ, Ghosh S, Okoli I, Mylonakis E: Antifungal activity of microbial secondary metabolites. PLoS One. 2011, 6: e25321-
Carberry S, Molloy E, Hammel S, O’Keeffe G, Jones GW, Kavanagh K, Doyle S: Gliotoxin effects on fungal growth: mechanisms and exploitation. Fungal Genet Biol. 2012, 49: 302-312.
Gallagher L, Owens RA, Dolan SK, O’Keeffe G, Schrettl M, Kavanagh K, Jones GW, Doyle S: The Aspergillus fumigatus protein GliK protects against oxidative stress and is essential for gliotoxin biosynthesis. Eukaryot Cell. 2012, 11: 1226-1238.
Shim H, Sup J, Kim J, Won S, Kwon HJ: Discovery of gliotoxin as a new small molecule targeting thioredoxin redox system. Biochem Biophys Res Commun. 2007, 359: 523-528.
Gardiner DM, Cozijnsen AJ, Wilson LM, Pedras MSC, Howlett BJ: The sirodesmin biosynthetic gene cluster of the plant pathogenic fungus Leptosphaeria maculans. Mol Microbiol. 2004, 53: 1307-1318.
Cramer RA, Gamcsik MP, Brooking RM, Najvar LK, Kirkpatrick WR, Patterson TF, Balibar CJ, Graybill JR, Perfect JR, Abraham SN, Steinbach WJ: Disruption of a nonribosomal peptide synthetase in Aspergillus fumigatus eliminates gliotoxin production. Eukaryot Cell. 2006, 5: 972-980.
Bok JW, Chung D, Balajee SA, Marr KA, Andes D, Nielsen KF, Frisvad JC, Kirby KA, Keller NP: GliZ, a transcriptional regulator of Gliotoxin Biosynthesis, contributes to Aspergillus fumigatus Virulence. Infect Immun. 2006, 74: 6761-6768.
Forseth RR, Fox EM, Chung D, Howlett BJ, Keller NP, Schroeder FC: Identification of cryptic products of the gliotoxin gene cluster using NMR-based comparative metabolomics and a model for gliotoxin biosynthesis. J Am Chem Soc. 2011, 133: 9678-9681.
Balibar CJ, Walsh CT: GliP, a multimodular nonribosomal peptide synthetase in Aspergillus fumigatus, makes the diketopiperazine scaffold of gliotoxin. Biochemistry. 2006, 45: 15029-15038.
Gardiner DM, Jarvis RS, Howlett BJ: The ABC transporter gene in the sirodesmin biosynthetic gene cluster of Leptosphaeria maculans is not essential for sirodesmin production but facilitates self-protection. Fungal Genet Biol. 2005, 42: 257-263.
Wang D, Toyotome T, Muraosa Y, Watanabe A, Wuren T, Bunsupa S, Aoyagi K, Yamazaki M, Takino M, Kamei K: GliA in Aspergillus fumigatus is required for its tolerance to gliotoxin and affects the amount of extracellular and intracellular gliotoxin. Med Mycol. 2014, 52: 506-518.
Scharf D, Remme N: Transannular disulfide formation in gliotoxin biosynthesis and its role in self-resistance of the human pathogen Aspergillus fumigatus. J Am Chem Soc. 2010, 132: 10136-10141.
Qin Z, Baker AT, Raab A, Huang S, Wang T, Yu Y, Jaspars M, Secombes CJ, Deng H: The fish pathogen Yersinia ruckeri produces holomycin and uses an RNA methyltransferase for self-resistance. J Biol Chem. 2013, 288: 14688-14697.
Li B, Walsh C: Streptomyces clavuligerus HlmI is an intramolecular disulfide-forming dithiol oxidase in holomycin biosynthesis. Biochemistry. 2011, 50: 4615-4622.
Gibbons JG, Beauvais A, Beau R, McGary KL, Latgé J-P, Rokas A: Global transcriptome changes underlying colony growth in the opportunistic human pathogen Aspergillus fumigatus. Eukaryot Cell. 2012, 11: 68-78.
Yu J, Fedorova ND, Montalbano BG, Bhatnagar D, Cleveland TE, Bennett JW, Nierman WC: Tight control of mycotoxin biosynthesis gene expression in Aspergillus flavus by temperature as revealed by RNA-Seq. FEMS Microbiol Lett. 2011, 322: 145-149.
Dhingra S, Andes D, Calvo AM: VeA regulates conidiation, gliotoxin production, and protease activity in the opportunistic human pathogen Aspergillus fumigatus. Eukaryot Cell. 2012, 11: 1531-1543.
Amich J, Schafferer L, Haas H, Krappmann S: Regulation of sulphur assimilation is essential for virulence and affects iron homeostasis of the human-pathogenic mould Aspergillus fumigatus. PLoS Pathog. 2013, 9: e1003573-
Dhingra S, Lind AL, Lin H-C, Tang Y, Rokas A, Calvo AM: The fumagillin gene cluster, an example of hundreds of genes under veA control in Aspergillus fumigatus. PLoS One. 2013, 8: e77147-
Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, Tetko I, Güldener U, Mannhaupt G, Münsterkötter M, Mewes HW: The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004, 32: 5539-5545.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30.
Mitsuguchi H, Seshime Y, Fujii I, Shibuya M, Ebizuka Y, Kushiro T: Biosynthesis of steroidal antibiotic fusidanes: functional analysis of oxidosqualene cyclase and subsequent tailoring enzymes from Aspergillus fumigatus. J Am Chem Soc. 2009, 131: 6402-6411.
Perrin RM, Fedorova ND, Bok JW, Cramer RA, Wortman JR, Kim HS, Nierman WC, Keller NP: Transcriptional regulation of chemical diversity in Aspergillus fumigatus by LaeA. PLoS Pathog. 2007, 3: e50-
Wiemann P, Guo C-J, Palmer JM, Sekonyela R, Wang CCC, Keller NP: Prototype of an intertwined secondary-metabolite supercluster. Proc Natl Acad Sci U S A. 2013, 110: 17065-17070.
Maiya S, Grundmann A, Li X, Li S-M, Turner G: Identification of a hybrid PKS/NRPS required for pseurotin A biosynthesis in the human pathogen Aspergillus fumigatus. Chembiochem. 2007, 8: 1736-1743.
Maiya S, Grundmann A, Li S-M, Turner G: The fumitremorgin gene cluster of Aspergillus fumigatus: identification of a gene encoding brevianamide F synthetase. Chembiochem. 2006, 7: 1062-1069.
Lin H-C, Chooi Y-H, Dhingra S, Xu W, Calvo AM, Tang Y: The fumagillin biosynthetic gene cluster in Aspergillus fumigatus encodes a cryptic terpene cyclase involved in the formation of β-trans-bergamotene. J Am Chem Soc. 2013, 135: 4616-4619.
Carberry S, Neville CM, Kavanagh KA, Doyle S: Analysis of major intracellular proteins of Aspergillus fumigatus by MALDI mass spectrometry: identification and characterisation of an elongation factor 1B protein with glutathione transferase activity. Biochem Biophys Res Commun. 2006, 341: 1096-1104.
Bok JW, Keller NP: LaeA, a Regulator of Secondary Metabolism in Aspergillus spp. Eukaryot Cell. 2004, 3: 527-535.
Lodeiro S, Xiong Q, Wilson WK, Ivanova Y, Smith ML, May GS, Matsuda SPT: Protostadienol biosynthesis and metabolism in the pathogenic fungus Aspergillus fumigatus. Org Lett. 2009, 11: 1241-1244.
Park H-S, Bayram O, Braus GH, Kim SC, Yu J-H: Characterization of the velvet regulators in Aspergillus fumigatus. Mol Microbiol. 2012, 86: 937-953.
Haas H: Iron - a key nexus in the Virulence of Aspergillus fumigatus. Front Microbiol. 2012, 3: 28-
Schrettl M, Bignell E, Kragl C, Sabiha Y: Distinct roles for intra-and extracellular siderophores during Aspergillus fumigatus infection. PLoS Pathog. 2007, 3: 1195-1207.
Yasmin S, Alcazar-Fuoli L, Gründlinger M, Puempel T, Cairns T, Blatzer M, Lopez JF, Grimalt JO, Bignell E, Haas H: Mevalonate governs interdependency of ergosterol and siderophore biosyntheses in the fungal pathogen Aspergillus fumigatus. Proc Natl Acad Sci U S A. 2012, 109: E497-E504.
Eisendle M, Oberegger H, Zadra I, Haas H: The siderophore system is essential for viability of Aspergillus nidulans: functional analysis of two genes encoding l-ornithine N 5-monooxygenase (sidA) and a non-ribosomal peptide synthetase (sidC). Mol Microbiol. 2003, 49: 359-375.
Arnaud MB, Chibucos MC, Costanzo MC, Crabtree J, Inglis DO, Lotia A, Orvis J, Shah P, Skrzypek MS, Binkley G, Miyasato SR, Wortman JR, Sherlock G: The Aspergillus Genome Database, a curated comparative genomics resource for gene, protein and sequence information for the Aspergillus research community. Nucleic Acids Res. 2010, 38 (Database issue): D420-D427.
Hartmann T, Cairns TC, Olbermann P, Morschhäuser J, Bignell EM, Krappmann S: Oligopeptide transport and regulation of extracellular proteolysis are required for growth of Aspergillus fumigatus on complex substrates but not for virulence. Mol Microbiol. 2011, 82: 917-935.
Bergmann A, Hartmann T, Cairns T, Bignell EM, Krappmann S: A regulator of Aspergillus fumigatus extracellular proteolytic activity is dispensable for virulence. Infect Immun. 2009, 77: 4041-4050.
Nishida H: Conservation of nucleosome positions in duplicated and orthologous gene pairs. ScientificWorldJournal. 2012, 2012: 298174-
Kupfahl C, Heinekamp T: The gliP gene of Aspergillus fumigatus results in loss of gliotoxin production but has no effect on virulence of the fungus in a low‒dose mouse infection model. Mol Microbiol. 2006, 62: 292-302.
Sekonyela R, Palmer JM, Bok J-W, Jain S, Berthier E, Forseth R, Schroeder F, Keller NP: RsmA regulates Aspergillus fumigatus gliotoxin cluster metabolites including cyclo(L-Phe-L-Ser), a potential new diagnostic marker for invasive aspergillosis. PLoS One. 2013, 8: e62591-
Li B, Forseth RR, Bowers AA, Schroeder FC, Walsh CT: A backup plan for self-protection: S-methylation of holomycin biosynthetic intermediates in Streptomyces clavuligerus. Chembiochem. 2012, 13: 2521-2526.
Dolan SK, Owens RA, O'Keeffe G, Hammel S, Fitzpatrick DA, Jones GW, Doyle S: Regulation of nonribosomal peptide synthesis: bis -thiomethylation attenuates gliotoxin biosynthesis in Aspergillus fumigatus. Chem Biol. 2014, 21: 999-1012.
Sanchez JF, Somoza AD, Keller NP, Wang CCC: Advances in Aspergillus secondary metabolite research in the post-genomic era. Nat Prod Rep. 2012, 29: 351-371.
Kato N, Suzuki H, Takagi H, Asami Y, Kakeya H, Uramoto M, Usui T, Takahashi S, Sugimoto Y, Osada H: Identification of cytochrome P450s required for fumitremorgin biosynthesis in Aspergillus fumigatus. Chembiochem. 2009, 10: 920-928.
Bok J, Balajee S, Marr K: LaeA, a regulator of morphogenetic fungal virulence factors. Eukaryot Cell. 2005, 4: 1574-1582.
Bayram O, Krappmann S, Ni M, Bok JW, Helmstaedt K, Valerius O, Braus-Stromeyer S, Kwon N-J, Keller NP, Yu J-H, Braus GH: VelB/VeA/LaeA complex coordinates light signal with fungal development and secondary metabolism. Science. 2008, 320: 1504-1506.
Wallner A, Blatzer M, Schrettl M, Sarg B, Lindner H, Haas H: Ferricrocin, a siderophore involved in intra- and transcellular iron distribution in Aspergillus fumigatus. Appl Environ Microbiol. 2009, 75: 4194-4196.
Xu XM, Møller SG: Iron-sulfur clusters: biogenesis, molecular mechanisms, and their functional significance. Antioxid Redox Signal. 2011, 15: 271-307.
Rutherford JC, Ojeda L, Balk J, Mühlenhoff U, Lill R, Winge DR: Activation of the iron regulon by the yeast Aft1/Aft2 transcription factors depends on mitochondrial but not cytosolic iron-sulfur protein biogenesis. J Biol Chem. 2005, 280: 10135-10140.
Sipos K, Lange H, Fekete Z, Ullmann P, Lill R, Kispal G: Maturation of cytosolic iron-sulfur proteins requires glutathione. J Biol Chem. 2002, 277: 26944-26949.
Kumar C, Igbaria A, D’Autreaux B, Planson A-G, Junot C, Godat E, Bachhawat AK, Delaunay-Moisan A, Toledano MB: Glutathione revisited: a vital function in iron metabolism and ancillary role in thiol-redox control. EMBO J. 2011, 30: 2044-2056.
Schrettl M, Beckmann N, Varga J, Heinekamp T, Jacobsen ID, Jöchl C, Moussa TA, Wang S, Gsaller F, Blatzer M, Werner ER, Niermann WC, Brakhage AA, Haas H: HapX-mediated adaption to iron starvation is crucial for virulence of Aspergillus fumigatus. PLoS Pathog. 2010, 6: e1001124-
Schrettl M, Kim HS, Eisendle M, Kragl C, Nierman WC, Heinekamp T, Werner ER, Jacobsen I, Illmer P, Yi H, Brakhage AA, Haas H: SreA-mediated iron regulation in Aspergillus fumigatus. Mol Microbiol. 2008, 70: 27-43.
Lee IR, Morrow CA, Fraser JA: Nitrogen regulation of virulence in clinically prevalent fungal pathogens. FEMS Microbiol Lett. 2013, 345: 77-84.
Wong KH, Hynes MJ, Davis MA: Recent advances in nitrogen regulation: a comparison between Saccharomyces cerevisiae and filamentous fungi. Eukaryot Cell. 2008, 7: 917-925.
Priebe S, Linde J, Albrecht D, Guthke R, Brakhage AA: FungiFun: A web-based application for functional categorization of fungal genes and proteins. Fungal Genet Biol. 2011, 48: 353-358.
Burns C, Geraghty R, Neville C, Murphy A, Kavanagh K, Doyle S: Identification, cloning, and functional expression of three glutathione transferase genes from Aspergillus fumigatus. Fungal Genet Biol. 2005, 42: 319-327.
O’Hanlon KA, Cairns T, Stack D, Schrettl M, Bignell EM, Kavanagh K, Miggin SM, O’Keeffe G, Larsen TO, Doyle S: Targeted disruption of nonribosomal peptide synthetase pes3 augments the virulence of Aspergillus fumigatus. Infect Immun. 2011, 79: 3978-3992.
Collins C, Keane TM, Turner DJ, O’Keeffe G, Fitzpatrick DA, Doyle S: Genomic and proteomic dissection of the ubiquitous plant pathogen, Armillaria mellea: toward a new infection model system. J Proteome Res. 2013, 12: 2552-2570.
Cox J, Mann M: MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008, 26: 1367-1372.
This work was supported by a Science Foundation Ireland Award to SD (SFI/11/PI/1188). SH was funded by the Irish Research Council. RAO was supported by 3U Biomedical Research (DCU-NUI Maynooth-RCSI). qRT-PCR instrumentation was funded by Science Foundation Ireland (SFI/07/RFP/GEN/F571/EC07). LFQ mass spectrometry facilities were funded by Science Foundation Ireland (12/RI/2346(3)). Agilent 6340 Ion Trap mass spectrometer was funded by the Higher Education Authority. We thank Prof. Nancy Keller (U. Wisconsin-Madison) for provision of relevant mutant strains.
The authors declare they have no competing interests.
GOK, SH, RAO, GWJ, and SD conceived and designed the experiments. SH and GOK extracted the RNA. TMK carried out library prep, sequencing and data mapping. Data analysis was performed by GOK, RAO, SH, TMK, DAF and SD. GOK and SH performed qPCR analysis. GOK carried out gliotoxin and secondary metabolite measurements and viability analysis. RAO performed the proteomics experimentation. GOK, GWJ and SD wrote the paper with contributions from RAO and TMK. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Gene sequencing information. Combined FPKM values for all replicates, Log2 fold changes and statistical significance data for all genes. Processed RNA-seq data. Relative expression of genes in A. fumigatus +/− gliotoxin; Relative expression of genes in A. fumigatus ∆gliT +/− gliotoxin. (ZIP 386 KB)
Additional file 2: Additional study information. Additional tables containing the functional categorisation of differentially regulated genes, Illumina RNA-seq summary statistics and primers used in the study. Additional Figure S1 is included. FunCat and KEGG categories of genes altered in expression; Summary Statistics for Illumina RNA-seq; Primers used for qRT-PCR; cell viability data. (DOCX 116 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
O’Keeffe, G., Hammel, S., Owens, R.A. et al. RNA-seq reveals the pan-transcriptomic impact of attenuating the gliotoxin self-protection mechanism in Aspergillus fumigatus. BMC Genomics 15, 894 (2014). https://doi.org/10.1186/1471-2164-15-894