Previous Article | Next Article ![]()
Eukaryotic Cell, May 2006, p. 794-805, Vol. 5, No. 5
1535-9778/06/$08.00+0 doi:10.1128/EC.5.5.794-805.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Gillian S. Khew, and
Meryl A. Davis
Department of Genetics, University of Melbourne, Parkville 3010, Australia
Received 16 December 2005/ Accepted 20 February 2006
|
|
|---|
|
|
|---|
In S. cerevisiae, oleate utilization depends on peroxisomal beta-oxidation. A wide range of peroxisomal proteins and the enzymes necessary for complete conversion of oleate to acetyl-CoA are induced by oleate (29). Acetyl-CoA produced in the peroxisome is converted to acetyl-carnitine for export to the cytosol and the mitochondrion (18, 39, 63, 79, 80). The related transcriptional activators Oaf1 and Pip2 containing DNA binding domains of the Zn2-Cys6 binuclear class together with Adr1, a C2H2 zinc finger protein, are responsible for this induction via the cis-acting elements (ORE) and UAS1, respectively (4, 25, 34, 35, 48, 57, 58). The genes for metabolism of the resulting acetyl-CoA by the glyoxalate bypass and gluconeogenesis is controlled by the Zn2-Cys6 binuclear proteins Cat8 and Sip4 acting at carbon source-responsive elements (CSREs) in the 5' region of the relevant genes (26, 27, 54, 55, 56, 81). Some of these genes are also regulated by Adr1 (86). Growth on ethanol or acetate as sole carbon sources is also dependent on the Cat8, Sip4, and Adr1 activators as well as the Snf1 kinase (82, 86). In the presence of glucose the Mig1 repressor represses the expression of these genes (23). Growth on acetate depends on acetyl-CoA synthetase-dependent formation of acetyl-CoA in the cytosol. This is converted to acetyl-carnitine by a cytosolic acetyl-carnitine transferase and transported into the mitochondrion, where it is converted back to acetyl-CoA by a mitochondrial acetyl-carnitine transferase for metabolism by the tricarboxylic acid cycle (18, 62, 80).
In the filamentous fungus Aspergillus nidulans, mutations in the facA gene, encoding acetyl-CoA synthetase, and in facC, encoding cytosolic acetyl-carnitine transferase, result in an inability to use acetate as a sole carbon source. These genes are induced by acetate via the FacB activator, which is related to the Cat8 and Sip4 proteins of S. cerevisiae (10, 31, 61, 67, 74, 76, 77). Mutations in acuD and acuE, which encode the glyoxalate bypass enzymes isocitrate lyase and malate synthase, respectively, are also unable to use acetate, but in addition are unable to use fatty acids as sole carbon sources (2, 3, 21, 36, 61). These genes are also regulated by acetate induction via FacB, but facB mutations do not prevent growth on fatty acids, and these genes are, in addition, regulated by fatty acid induction (2, 5). Similarly, acuJ encoding an acetyl-carnitine transferase located in both mitochondria and peroxisomes is required for growth on both acetate and fatty acids because it is necessary for the shuttling of acetyl-CoA via the transporter AcuH (14, 50, 67). The acuJ gene is also subject to FacB-mediated acetate induction as well as fatty acid induction (67; M. J. Hynes, A. Andrianopoulos, S. Delimitrou, S. L. Murray, H. Sealy-Lewis, and M. A. Davis, unpublished data). This is also the case for idpA, the gene for NADP-isocitrate dehydrogenase located in peroxisomes and in mitochondria (68).
A. nidulans is able to utilize short-chain and long-chain fatty acids, and recently it has been shown that beta-oxidation pathways exist in both mitochondria and peroxisomes, as is the case for mammals and plants. There is evidence for fatty acid induction of some of the relevant enzymes in A. nidulans (49, 78) and in Neurospora crassa (28, 38). Nothing is known about the mechanism or the regulatory genes involved. It is therefore predicted that there are three classes of genes: (i) acetate specific, regulated by FacB; (ii) acetate and fatty acid induced and controlled by FacB and unknown fatty acid regulators, and (iii) fatty acid specific induced only by the fatty acid regulators. This latter class is predicted to include enzymes required for fatty acid breakdown via beta-oxidation as well as proteins involved in peroxisome biogenesis and function.
We have sought to investigate this regulatory circuit by isolating mutants unable to use fatty acids as carbon sources but still capable of growth on acetate and by cloning the corresponding genes. Among the genes identified are ones predicted to encode Zn2-Cys6 binuclear proteins. These have been shown to control the expression of genes involved in fatty acid catabolism but not to affect expression of genes involved specifically in acetate utilization. One of these genes affects induction by short-chain (C2 to C6) fatty acids and is not highly conserved in other fungi. Two other genes, one affecting short-chain induction and one affecting induction by all chain lengths, encode related proteins with highly conserved orthologs found in other euascomycetes and a single ortholog in the hemiascomycetes Candida albicans, Debaryomyces hansenii, and Yarrowia lipolytica but not in other hemiascomycete lineages. These proteins have been found to bind in vitro to sequences found to be overrepresented in the 5' region of a large number of genes predicted to be involved in fatty acid breakdown in the species which contain predicted orthologs but not in the 5' region of the corresponding genes in S. cerevisiae. We have therefore identified a highly conserved regulatory circuit controlling the use of fatty acids in ascomycete fungi.
|
|
|---|
Molecular techniques. Standard methods for DNA manipulations, RNA isolation, nucleic acid blotting, and hybridization have been described previously (12, 60, 75). The sequences of oligonucleotides used are given in Table 1.
|
View this table: [in a new window] |
TABLE 1. Primers used in this study
|
Isolation of mutants. A conidial suspension of strain MH9975 was treated with diepoxyoctane as described previously (32), and dilutions plated on glucose-minimal ammonium tartrate medium with 0.08% sodium deoxycholate were added to restrict the colony size. Colonies were velvet replicated to minimal medium containing 20 mM sodium butyrate as the sole carbon source with 10 mM ammonium chloride as the nitrogen source. Potential mutants affected in butyrate utilization were isolated and tested for their ability to use L-glutamate (50 mM) and acetate (50 mM) as carbon sources. Mutants specifically affected in butyrate utilization were isolated for further study. This yielded strains containing the scfA7 and farB7 mutations. One of the butyrate-nonutilizing mutants (like a number of others) was completely inhibited in growth on medium containing fatty acids. This mutant subsequently was shown to have a mutation in a gene (pexF) encoding an ortholog of the peroxin Pex6 (42, 72, 83). This strain was inoculated onto minimal medium containing 0.5% lactose and 0.2% Tween 80, as a source of the fatty acid oleate, and ammonium chloride. Spontaneous resistant sectors appearing after incubation for 3 to 5 days were isolated. Outcrossing of these yielded isolates separated from the original pexF mutation that were unable to use fatty acids as sole carbon sources but did not show the extreme pexF phenotype. This yielded a strain with the farA1 mutation.
Molecular cloning and sequencing of genes. Strains containing the proposed regulatory mutations were crossed to a strain containing the pyrG89 mutation to generate PyrG double mutants. These strains were used for cloning the genes by functional complementation by transforming with a genomic library in the autonomously replicating vector pRG3AMA1 (51). Transformants selected for uracil-uridine prototrophy were velvet replicated to medium containing butyrate as the sole carbon source, complementing transformants were selected and purified, and genomic DNA was prepared. Plasmid DNA was recovered by transforming this DNA into E. coli, selecting for ampicillin resistance. Alternatively, direct PCR on genomic DNA was used to recover insert DNA from complementing plasmids. Subcloning, sequencing, and comparison with the A. nidulans database (http://www.broad.mit.edu/annotation/fungi/aspergillus/) showed that scfA corresponded to the annotated gene AN1303.2, farA corresponded to AN7050.2, and farB corresponded to AN1425.2.
cDNA corresponding to the genes was synthesized from total RNA using gene-specific primer pairs LT1:RT5' and LT1:RT3', BUZ7RT5' and BUZ7RT3', and BUX75'RT and BUX7RTPCR3' for farA, farB, and scfA, respectively, using the Superscript II (Invitrogen) one-step reverse transcriptase PCR kit. Products were cloned into the plasmid pGEMTEasy (Promega Corp.) and sequenced. This led to a reannotation of intron positions in the genes (see Fig. 2 and 7, below).
![]() View larger version (33K): [in a new window] |
FIG. 2. Structures of farA and farB genes and comparisons with orthologs in other fungi. (A) The farA gene, showing intron positions and the sequences coding for the indicated domains. farA was generated by replacing sequences from +225 to +2670 (relative to the start codon) with the A. nidulans pyrG gene followed by gene replacement. (B) Structure of the farB gene, showing intron positions and the sequences coding for the indicated domains. farB was generated by replacing sequences from 227 to +2368 (relative to the start codon) with the A. nidulans pyrG gene followed by gene replacement. (C) Comparison of the Zn2-Cys6 binuclear cluster domains of the FarA and FarB proteins in filamentous fungi and the orthologs in hemiascomycetes compared to A. nidulans FarA (Yeast group). Identical residues present in at least 60% of the sequences are indicated by black boxes, whereas gray shading represents similar residues. Overall comparisons of the proteins, database accession numbers, and intron positions in the genes are presented in the Fig. S1 to S3 in the supplemental material. (D) Unrooted neighbor joining distance tree showing the relationships of the FarA, FarB, and yeast proteins to each other. Sequences were aligned using ClustalW (70) and distance matrices were calculated (19), which were used to create the tree in NJPlot (52). Species abbreviations: A.n, Aspergillus nidulans; A.f; Aspergillus fumigatus; S.n, Stagonospora nodorum; F.g, Fusarium graminearium; N.h, Nectria hematococca; M.g, Magnaporthe grisea; N.c, Neurospora crassa; C.g, Chaetomium globosum; B.c, Botrytis cineria; C.a, Candida albicans; D.h, Debaryomyces hansenii; Y.l, Yarrowia lipolytica.
|
![]() View larger version (22K): [in a new window] |
FIG. 7. Structure of the scfA gene. (A) The scfA gene, showing intron positions and the sequences coding for the indicated domains. scfA was generated by replacing sequences from 295 to +1470 (relative to the start codon) with the A. nidulans pyrG gene followed by gene replacement. (B) Comparison of the Zn2-Cys6 binuclear cluster domains of scfA with proposed orthologs in N. crassa and A. fumigatus. Identical residues present in at least 60% of the sequences are indicated by black boxes, whereas gray shading represents similar residues. Overall comparisons of the proteins and database accession numbers are presented in Fig. S5 in the supplemental material.
|
Expression of MBP fusion proteins. Fragments of farA, farB, and scfA genes were expressed as maltose binding protein (MBP) fusion proteins using the pMAL system (New England BioLabs). The cDNA clones of farA, farB, and scfA, pSM6128, pSM6123, and pSM5777, respectively, were used as templates for PCRs using primers to introduce a restriction site enabling the cloning of complete cDNAs of farA, farB, and scfA into the pMALC2 polylinker in the same reading frame as malE. The resulting plasmids, pSM6350 (encoding MBP fused to amino acids 1 to 399 of FarA), pSM6351 (encoding MBP fused to amino acids 1 to 514 of FarB), and pSM6229 (encoding MBP fused to amino acids 1 to 410 of ScfA), together with the previously described pRT2013 (encoding MBP fused to amino acids 4 to 417 of FacB [74]) were transformed into E. coli strain BL21-DE3 or Rosetta (Novagen) for expression. An overnight culture was diluted 1:100 in M9ZB medium, with the addition of 1 mM MgSO4, 0.4% glucose, and grown at 37°C with aeration to an optical density at 600 nm of 0.6. Expression of fusion plasmids was induced by the addition of isopropyl-ß-D-thiogalactopyranoside to a 1 mM final concentration for 2 h. Cells were harvested, and crude extracts were produced as described in the pilot experiment of the pMAL protein fusion and purification system manual (New England BioLabs). Protein concentrations were determined using the Bio-Rad protein assay reagent (Bio-Rad). Expression of fusion proteins of the predicted sizes was visualized by sodium dodecyl sulfate-polyacrylamide gel electrophoresis.
Electrophoretic mobility shift assays (EMSAs). Oligonucleotide probes were prepared, and binding reactions were carried out as previously described except that no poly(dI-dC) or Mg2+ was added (74). Electrophoresis of samples occurred on a nondenaturing 5% polyacrylamide gel in 0.5x Tris-borate-EDTA at 100 to 150 V. Gels were blotted onto 3MM paper, dried in a vacuum dryer, and exposed to X-ray film.
Nucleotide sequence accession numbers. The sequences determined in this study have been deposited in GenBank and assigned the accession numbers DQ386636 [GenBank] (scfA), DQ386637 [GenBank] (farA), and DQ386638 [GenBank] (farB).
|
|
|---|
|
View larger version (33K): [in a new window] |
FIG. 1. Growth of regulatory gene mutants on fatty acids as sole carbon sources. The following carbon sources were added to minimal medium with 10 mM ammonium chloride as the nitrogen source: glucose (1%); acetate (50 mM); butyrate and valerate (10 mM); hexanoate (5 mM), propionate, Tween 20, and Tween 80 (0.1%); lauric, myristic, palmitic, heptadecanoic, stearic, oleate, elaidic, and erucic acids (2.5 mM). Growth was for 2 to 3 days at 37°C.
|
The farA and farB genes were cloned by complementation by the isolation of double mutants containing the pyrG89 mutation, transforming with a genomic library in the autonomously replicating vector pRG3AMA1 (51), selecting for PyrG+ transformants, and replica plating to medium with butyrate as the sole carbon source. Colonies capable of growth were purified and DNA isolated. The corresponding autonomously replicating plasmids were recovered by transformation into E. coli. Sequencing of insert DNA in the recovered plasmids showed that farA corresponded to the annotated gene AN7050.2, while farB corresponded to AN1425.2 in the genome sequence (http://www.broad.mit.edu/annotation/fungi/aspergillus/). cDNA of each gene, made by reverse transcription-PCR, was cloned and sequenced, leading to a reannotation of intron positions in each gene (Fig. 2A and B). Each gene was deleted by replacement of sequences with the pyrG gene as shown in Fig. 2. The generated mutations segregated as single genes, and each deletion mutant was complemented by subclones of the relevant cloned genes, indicating that no additional events occurred during the isolation of the deletions.
The corresponding phenotypes with respect to fatty acid utilization of each gene deletion were identical to those of the original isolated mutants (Fig. 1B). More extensive testing of the deletion mutants showed that farB
resulted in loss of utilization of short-chain fatty acids but not inhibition of growth and did not affect growth on long-chain fatty acids. The farA deletion resulted in loss of the utilization of short-chain fatty acids but also resulted in loss of utilization and inhibition on all fatty acids tested longer than C10 (lauric acid) up to C22 (erucic acid), including the odd-chain-length heptadecanoate (C17) and the unsaturated oleic and elaidic acids (C18). The loss-of-function facB101 mutation (2, 76) did not affect growth on fatty acids apart from a slight effect on butyrate. Growth on propionate (C3) was inhibited in the mutants as well as the wild type and the facB101 mutant was resistant, consistent with previous results showing that propionyl-CoA generated from propionate is inhibitory to growth (6).
The farA and farB genes encode conserved proteins containing fungal-specific transcription factor domains. Each of the predicted proteins encoded by the two genes were found to contain both a Zn2-Cys6 binuclear cluster domain and a conserved central domain characteristic of fungal transcription factors containing this DNA binding domain (Fig. 2A, B, and C) (62, 73). Database searches revealed predicted conserved orthologs corresponding to each of these genes in filamentous euascomycetes, including a number of plant pathogens (see Fig. S1 and S2 in the supplemental material). In the case of farA the first four introns, including two introns within the sequence encoding the Zn2-Cys6 domain, were found to be conserved in all genes and to have led to extensive misannotation in the databases (see Fig. S4 in the supplemental material). The FarA and FarB proteins showed similarity to each other, and each of the species also had a single copy of each gene, with the exception of Aspergillus fumigatus, which was found to have two predicted copies of the farB gene, apparently arising due to a segmental duplication of approximately 22.6 kbp encoding seven hypothetical annotated genes from chromosome 8 to chromosome 1 (see Fig. S2 in the supplemental material). A single copy of a related gene was found in the hemiascomycete species Candida albicans, Debaryomyces hansenii, and Yarrowia lipolytica, but not in Ashbya (Eremothecium) gossypii nor in S. cerevisiae (see Fig. S3 in the supplemental material). Y. lipolytica had one of the introns conserved in other species, while the C. albicans and D. hansenii genes lacked introns (see Fig. S4 in the supplemental material). The neighbor joining tree constructed showed that the hemiascomycete proteins were most closely related to the FarA group, which were tightly clustered in the euascomycetes. The FarB group was more divergent (Fig. 2D).
In the filamentous ascomycetes, the proposed Zn2-Cys6 binuclear DNA binding domains of these proteins were found to be extraordinarily conserved, particularly within the FarA group, and to be related between groups and with the corresponding domains in the hemiascomycete proteins (Fig. 2C). Of particular note was the finding that the predicted proteins were highly similar to the cutinase transcription factors CTF
and CTFß of Nectria hematococca. These proteins have been found to bind to sequences in the 5' region of cutinase genes in this plant pathogen (43, 44). Cutin is a polymer containing hydroxy-oleate, strongly suggesting a connection with fatty acid breakdown.
Effects of the farA and farB genes on induction of an acuJ-lacZ reporter construct.
The acuJ gene has been found to encode a carnitine-acetyltransferase located to peroxisomes and mitochondria and is regulated by FacB-dependent acetate induction and induced by fatty acids (50, 67; Hynes et al., unpublished). A sequence from 82 to 349 from the acuJ gene was inserted upstream of the minimal gpd promoter driving lacZ expression (53), and this was inserted in single copy at the argB locus (Fig. 3A). For studies of induction, proline was used as a noninducing limiting carbon source control. The original farA and farB mutations as well as the deletions were found to affect induction with farB, reducing induction by short-chain fatty acids but not by oleate or Tween 20 (a source of lauric, palmitic, and myristic acids), while the farA mutations led to elevated levels of noninduced expression but reduced induction by all fatty acids tested (Fig. 3B, C, and H). This was also observed for the farA
farB
double mutant (Fig. 3G). In short time course induction experiments, butyrate and oleate induction was detectable at about 60 min (Fig. 3D and E). Acetate and propionate were weak inducers, and this induction was not affected by the facB101 mutation. In fact, stronger induction in this background as well as in a facA303 background (lacking acetyl-CoA synthetase) was observed (Fig. 3G). This suggested that loss of activation of acetate and propionate to the respective acyl-CoA might enhance the weak induction observed. Induction was observed in the presence of 1% glucose, and a loss of function mutation in the creA gene, encoding the major glucose repressor protein (15), only slightly increased the response to an inducer (Fig. 3F).
![]() View larger version (49K): [in a new window] |
FIG. 3. Effects of regulatory gene mutations on expression of an acuJ-lacZ reporter. (A) A fragment of the acuJ 5' region was inserted upstream of a minimal promoter (gpd mini) driving lacZ expression. The black bars represent CCGAGG sequences present in this sequence. (B to H) Mycelia of strains of the indicated genotypes were grown for 16 to 18 h in 1% glucose-10 mM ammonium tartrate medium and then transferred to minimal medium with the indicated carbon sources together with 10 mM ammonium chloride for 6 h or for the indicated times (D and E) before harvesting. Mycelia were extracted and assayed for ß-galactosidase as previously described (13). The specific activity is expressed in Miller units per minute per milligram of protein. The bars represent standard errors. The concentrations of carbon sources were 10 mM for proline, acetate, butyrate, valerate, and oleate, 5 mM for hexanoate, 0.1% for propionate and Tween 20, and 2.5 mM for myristic, palmitic, heptadecanoic, elaidic, stearic, and erucic acids.
|
and farB
strains, which we interpret as indicating that butyrate must be metabolized to acetyl-CoA to result in induction via FacB. Long-chain (Tween 20 and oleate) induction was not observed. Similarly, genes required for both acetate and fatty acid utilization, acuD, acuE, acuH, and acuJ, were induced by acetate and butyrate in a FacB-dependent fashion. However, these genes were also long-chain fatty acid inducible, dependent on FarA. The other genes investigated were predicted to be involved specifically in fatty acid utilization, and these were not induced by acetate but were induced by butyrate (FarB and FarA dependent) and Tween 20 and oleate (FarA dependent).
![]() View larger version (67K): [in a new window] |
FIG. 4. Analysis of the effects of regulatory gene mutations on expression of various genes. RNAs were extracted from strains of the indicated genotypes grown for 16 h in 1% glucose-10 mM ammonium tartrate medium and then transferred to medium with the indicated carbon sources for 6 h. Abbreviations: pro, 10 mM L-proline; ace, 10 mM acetate; but, 10 mM n-butyrate; ole, 10 mM oleate; Tw20, 0.1% Tween 20. The results represent three different Northern blot membranes probed multiple times with sequences corresponding to the indicated genes with their proposed functions. Probes were made by labeling either gel-purified restriction fragments or PCR products generated using primers based on genome sequences. H3 represents an EcoRI fragment of the histone H3 clone (17), and rRNA refers to the large rRNA species observed by ethidium bromide staining.
|
![]() View larger version (26K): [in a new window] |
FIG. 5. EMSA analysis of DNA binding by expressed MBP fusion proteins. (A) Summary of the key DNA binding studies of cutinase transcription factors in N. hematococca (43, 44). Binding was detected to the probe Pal2 but not Pal2M1 or Pal2M2. Methylation protection studies showed the protection of the G residues, indicated by asterisks. (B) Probes used in EMSAs in this study. Complementary pairs of oligonucleotides were designed based on sequences from the 5' regions of the genes acuJ, echA, acuE, and pexK, with Spe1/Xba1- or BamH1/BglII-compatible ends added (in lowercase letters). The database designation of these genes is given in Table S1 in the supplemental material. The conserved core CCGAGG sequence is indicated in bold, and only one strand is shown. Probe J3-4 is identical to J1-2 except for the modification of the core CCGAGG to CCcgGG. (C and D) EMSA of extracts (20 µg of protein added to binding reaction mixtures) to the indicated probes together with no-extract controls. Free probe is designated as F. Binding is represented by the decreased mobility of the probe in the presence of protein, compared to the free probe. (E) Specific competition of binding of FarB and FarA to probe J7-8. Binding was performed with MBP-FarA and MBP-FarB extracts (20 µg of protein) in the presence of increasing amounts of unlabeled competitor DNA. Competitors were the unlabeled probes J1-2 and J3-4, increasing from 10 to 100 to 500 times the concentration of labeled probe.
|
Unlabeled DNA (J3-4) containing the mutant sequence CCcgGG replacing the core sequence showed a loss of competition for binding of MBP-FarA and MBP-FarB with the J7-8 probe, although some competition for MBP-FarB was observed (Fig. 5E). Overall the results showed, in agreement with previous observations for the N. hematococca CTF proteins, that the sequence CCGAGG (complement, CCTCGG) is the core binding sequence for these conserved proteins. It was therefore predicted that this sequence would be found in the 5' regions of genes induced by fatty acids.
Analysis of 5' sequences of fungal genes for predicted binding sites. Sequences from the 5' region (1 kb upstream of the start codon) of a large number of genes from the A. nidulans genome were scanned for the CCTCGG sequence. One or more sequences were found in nearly all genes predicted to be involved in fatty acid metabolism (see Table S1 in the supplemental material). These included genes determining proteins involved in peroxisome biogenesis or function, enzymes of beta-oxidation, including those predicted to be in mitochondria or peroxisomes, metabolite transporters, and glutathione peroxidases (but not catalases) involved in dealing with oxidative stress. The motif was found to occur in either orientation, and at least one copy was most commonly found within 300 bp of the start codon.
As part of the comparative analysis of the genomes of A. nidulans, Aspergillus oryzae, and A. fumigatus (22), a three-way alignment of 1 kb upstream of orthologous genes using Mlagan was performed. The CCTCGG motif was found to commonly occur in high-scoring conserved sequences upstream of relevant genes (these are indicated in Table S1 in the supplemental material). Furthermore, analysis of the patterns of occurrence of these conserved sequences ("conpat" analysis) identified the CCTCGG sequence to be enriched in genes for lipid metabolism and peroxisome localization (see Table 3 in reference 22).
Analysis of the 5' regions of genes required for acetate utilization showed the presence of FacB sites in all genes, while the CCTCGG motif was present only in genes also regulated by fatty acids (Fig. 6). A more limited analysis of 1 kb upstream of genes from other species containing farA and farB orthologs showed that the CCTCGG motif was overrepresented in those genes that might be expected to be regulated by fatty acids but not in acetate-specific genes (Table 2; see also Table S1 in the supplemental material). In the case of the plant pathogen M. grisea, this included members of the cutinase gene family, and the motif was also found upstream of cutinase genes as well as a predicted lipase in A. nidulans. S. cerevisiae lacks farA and farB orthologs, and fatty acid induction is dependent on the Oaf1/Pip2 heterodimer (29, 34, 35, 58). Consistent with this, the CCTCGG motif was not overrepresented in the 5' region of relevant genes from S. cerevisiae (Table 2).
![]() View larger version (12K): [in a new window] |
FIG. 6. Positions of FacB binding sites and CCTCGG motifs in the 5' regions of genes involved in acetate utilization in A. nidulans. FacB binding sites are based on those previously proposed (67, 74) or by sequence inspection in the case of acuJ. Vertical lines represent 100-bp intervals upstream of the start codon for each gene.
|
|
View this table: [in a new window] |
TABLE 2. Presence of CCTCGG motif in the 5' region of relevant genes from fungal speciesa
|
strains were found to be the same as for scfA7 mutant strains and very similar to those for farB
strains (Fig. 1B). Furthermore, scfA
affected short-chain but not long-chain fatty acid induction of the acuJ-lacZ reporter (Fig. 3) as well as showing a similar pattern to the effects of farB
on transcription of the genes studied by Northern blotting (Fig. 4). Expression of ScfA as an MBP fusion protein followed by EMSA with the probes used to detect FarA and FarB binding failed to give binding (Fig. 5). Additional probes based on sequences present in the acuJ 82 to 349 region have so far also failed to give significant binding (results not shown). |
|
|---|
The FarA and FarB proteins are highly conserved in filamentous ascomycetes and have related DNA binding domains. The previously described cutinase transcription factors from N. hematococca (43, 44), which were defined by DNA binding studies using 5' sequences from cutinase genes, are clearly orthologs and, in agreement with these studies, we have shown that expressed fusion proteins of FarA and FarB bind in vitro to sequences containing the core CCTCGG sequence. Analysis of fungal genomes has shown that this sequence is present in one or more copies within 1 kb of the predicted start codon of a large number of genes predicted to encode proteins involved in fatty acid metabolism, providing strong support for these regulatory genes being involved in induction by fatty acids. Furthermore, comparative analyses of the genomes of Aspergillus species indicate a conservation of this sequence in the upstream region of orthologous genes enriched for lipid metabolism and peroxisomal functions (22). Our findings on the in vivo functions of FarA and FarB provide a striking validation of the value of the phylogenetic footprinting methods used to determine potential functional motifs. The ease of detection of this simple sequence makes it possible to predict in many ascomycete species those genes likely to be regulated by fatty acid induction via these proteins.
In addition to the glyoxalate cycle enzymes, we have found direct evidence for FarA and FarB regulating genes involved in the shuttling of acetyl-CoA between peroxisomes and mitochondria (acuJ and acuH) as well as genes encoding both mitochondrial and peroxisomal beta-oxidation enzymes. We have also found that an ortholog of S. cerevisiae PEX11, predicted to encode a peroxin required for peroxisome proliferation (42, 72, 83), is regulated by FarA and FarB, and we have shown that deletion of this gene in A. nidulans results in reduced growth on fatty acids (G. S. Khew and M. J. Hynes, unpublished data). The CCTCGG motif was found upstream of a number of predicted orthologs for key peroxisomal proteins, including Pex3, which is absolutely required for de novo peroxisome formation (30, 41), as well as Pex1, Pex6, and Pex5, which are required for protein import into the peroxisomal matrix. We have isolated strains containing mutations in these genes and found that the resulting loss of peroxisomal functions leads to mild developmental defects (e.g., reduced sporulation) and complete inhibition of growth by both long- and short-chain fatty acids (Khew and Hynes, unpublished). The farA
and farB
strains do not show these phenotypes. This indicates that the Pex phenotypes result from mislocalization of induced fatty acid metabolic activities, and these effects are not observed when induction is lost. In addition, it is possible that induction of these peroxisomal functions is not absolutely required for expression but rather serves to increase the capacity for peroxisome formation and function.
Filamentous ascomycetes are very versatile in their ability to use diverse carbon sources. The finding that the conserved CCTCGG motif is present in the upstream region of the cutinase gene families of the plant pathogens M. grisea and N. hematococca as well as the saprophyte A. nidulans indicates that regulation of these genes is via fatty acid induction. We have not extensively scanned potential lipase genes, but in one case in A. nidulans we found the 5' motif. In an expressed sequence tag analysis of the insect pathogen Metarhizium anisopliae grown under conditions of maximal cuticle degradation, lipases and phospholipases were detected and, significantly, an ortholog of the N. hematococca cutinase transcription factor was found (20). This result is now understandable, since it is highly likely that this fungal species uses orthologs of farA and farB to control enzymes involved in lipid degradation. A broad role for these transcription factors in controlling the expression of enzymes for degradation of sources of fatty acids in filamentous fungi is strongly suggested.
At this stage many questions relating to the mechanism of action of these transcription factors remain to be answered. Deletion of farA results in loss of the ability to use all fatty acids tested and loss of induction by fatty acids. However, the noninduced activity of the acuJ-lacZ reporter was found to be elevated, implying a repressing function for FarA in the absence of inducer. Deletion of farB abolishes growth on and induction by short-chain fatty acids but does not affect long-chain induction. Somewhat increased induction by long-chain fatty acids of the acuJ-lacZ reporter was also observed in the farB
background (Fig. 3). Given that FarA and FarB recognize the same DNA sequence, this may indicate competition for binding between inducer-activated and nonactivated proteins. Deletion of scfA gives a very similar phenotype to that of the farB
mutant in causing loss of growth on and induction by short-chain fatty acids but does not result in increased long-chain fatty acid induction. We have so far been unable to determine an ScfA binding site, and it is puzzling that, in stark contrast to FarA and FarB, highly conserved orthologs of this protein have not been found in other species. The possibilities that ScfA acts upstream of FarB by controlling the production of the short-chain inducing signal or that it is required for FarB expression are currently being examined. The nature of the inducing signal and the discrimination between short-chain and long-chain inducers is of considerable interest. Response of the acuJ-lacZ reporter to inducer is relatively rapid (Fig. 3), perhaps indicating that the fatty acid inducer does not require extensive metabolism to generate the signal. In S. cerevisiae it has been proposed that the activators Oaf1 and Pip2 respond directly to oleate as inducer, but direct evidence is lacking (29). Zn2-Cys6 binuclear cluster proteins have been proposed to commonly respond directly to inducing metabolites (65). It should also be noted that as beta-oxidation of long-chain fatty acids proceeds, the shortening of chain length might result in an increasing level of short-chain induction. Analysis of the A. nidulans genome reveals multiple genes potentially carrying out enzymatic steps in beta-oxidation, and the substrate and induction specificity of these may be complex.
The finding that the hemiascomycete yeasts C. albicans, D. hansenii, and Y. lipolytica contain orthologs of FarA is intriguing. Other hemiascomycetes have orthologs of the well-studied Oaf1 and Pip2 genes of S. cerevisiae. The overrepresentation of the CCTCGG motif in the 5' region of relevant genes in these species, but not in S. cerevisiae, provides strong circumstantial evidence for the functional role of these FarA orthologs. The phylogenetic analysis presented (Fig. 2D) suggests one scenario in which a FarA orthog involved in long-chain induction was present in a progenitor of the ascomycetes. In filamentous euascomycetes, this was duplicated and diverged, resulting in the FarB genes being involved in a greater metabolic capacity to respond to short-chain substrates. FarA may have been lost and replaced by Oaf1/Pip2 orthologs in the other hemiascomycetes. Peroxisome biogenesis and alkane utilization in Y. lipolytica have been the subject of much study (e.g., references 24, 71, and 72), and investigation of the regulatory role of the FarA ortholog in this species is required.
Recent studies have suggested that peroxisomal functions and lipid metabolism are important for fungal pathogenesis (33, 37, 45-47, 66, 69, 85). Our discovery in plant and animal pathogens of transcription factors that are highly likely to be involved in controlling these processes is therefore significant. Pathogens may use these regulatory proteins to respond to host fatty acid signals to upregulate functions necessary for the establishment of infection. Deletion of the genes in these species and testing pathogenicity will therefore be of considerable interest.
. .
Assistance by Khanh Nguyen, Rosemary Genovese, Sophie Delimitrou, and Oliver Evans is acknowledged.
Supplemental material for this article may be found at http://ec.asm.org/. ![]()
Present address: Department of Physiology, University of Melbourne, Parkville 3010, Australia. ![]()
|
|
|---|
strain. J. Biol. Chem. 280:34350-34357.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»