Comprehensive Analysis of Aspergillus nidulans PKA Phosphorylome Identifies a Novel Mode of CreA Regulation.

In filamentous fungi, an important kinase responsible for adaptation to changes in available nutrients is cyclic AMP (cAMP)-dependent protein kinase (protein kinase A [PKA]). This kinase has been well characterized at a molecular level, but its systemic action and direct/indirect targets are generally not well understood in filamentous fungi. In this work, we used a pkaA deletion strain (ΔpkaA) to identify Aspergillus nidulans proteins for which phosphorylation is dependent (either directly or indirectly) on PKA. A combination of phosphoproteomic and transcriptomic analyses revealed both direct and indirect targets of PKA and provided a global perspective on its function. One of these targets was the transcription factor CreA, the main repressor responsible for carbon catabolite repression (CCR). In the ΔpkaA strain, we identified a previously unreported phosphosite in CreA, S319, which (based on motif analysis) appears to be a direct target of Stk22 kinase (AN5728). Upon replacement of CreA S319 with an alanine (i.e., phosphonull mutant), the dynamics of CreA import to the nucleus are affected. Collectively, this work provides a global overview of PKA function while also providing novel insight regarding significance of a specific PKA-mediated phosphorylation event.IMPORTANCE The cyclic AMP (cAMP)-dependent protein kinase A (PKA) signaling pathway is well conserved across eukaryotes, and previous work has shown that it plays an important role in regulating development, growth, and virulence in a number of fungi. PKA is activated in response to extracellular nutrients and acts to regulate metabolism and growth. While a number of components in the PKA pathway have been defined in filamentous fungi, current understanding does not provide a global perspective on PKA function. Thus, this work is significant in that it comprehensively identifies proteins and functional pathways regulated by PKA in a model filamentous fungus. This information enhances our understanding of PKA action and may provide information on how to manipulate it for specific purposes.

protein kinase A (PKA) pathway. In fungi, PKA signaling mediates a wide range of processes, including nutrient sensing (3), stress responses (4), regulation of metabolism (5), as well as development and pathogenicity (6,7). Although the molecular basis for PKA action is well understood, with more than a hundred different crystal structures published (8), considerably less is known about its systemic action and impact on downstream effectors (9,10). This is particularly true for filamentous fungi.
In Aspergillus nidulans, CCR is shown to be mediated mainly by the transcription factor CreA, a Cys2-His2 type DNA-binding zinc finger protein (27,28). CreA is responsible for repressing the transcription of numerous genes encoding lignocellulolytic enzymes when glucose is present, thus preventing the utilization of alternative carbon sources (29). When the cell senses glucose, CreA is imported into the nucleus, binds to DNA, and represses transcription of several genes. However, the mechanisms that direct the localization of CreA, and its affinity for binding to DNA remain somewhat unclear in filamentous fungi (12). Previous work has shown that Trichoderma reesei CreA homologue Cre1 is phosphorylated by casein kinase II at Ser 241, which is essential for DNA binding (30). However, it has been suggested that Cre1 nuclear import is regulated by a different mechanism, perhaps via interactions with other proteins or protein turnover (30). In S. cerevisiae, nuclear localization of the CreA functional homologue Mig1p is regulated through phosphorylation by the kinase Snf1p (31). In the presence of glucose, Snf1p is inactivated, resulting in dephosphorylation of Mig1 and its retention in the nucleus (32). In A. nidulans, there is a lack of evidence confirming direct phosphorylation of CreA. However, some studies suggest that kinases in A. nidulans are responsible for CreA cellular localization because the deletion of snfA or schA prevents CreA from leaving the nucleus in derepressing conditions (33). While CreA contains sites that are predicted to be phosphorylated by PKA, their functional relevance has not yet been investigated (12).
Here we utilize a comprehensive phosphoproteomic screen to identify both direct and indirect targets of PKA. One of these targets is CreA, which is also found to be regulated at a transcriptional level. We identify a unique phosphorylation site on this protein which has not been previously described. Upon further functional investigation, we find that this site is required for the proper dynamics of CreA nuclear import during CCR.

RESULTS AND DISCUSSION
To better understand the systemic action of PKA in A. nidulans, we carried out phosphoproteomic and transcriptomic differential analyses of two A. nidulans strains: a pkaA deletion strain (ΔpkaA) and its isogenic parent (34).
Phosphoproteomic analysis. To conduct phosphoproteomic analysis, we grew fungi in glucose to induce the PKA pathway. Our analysis resulted in the identification of 847 differentially phosphorylated peptides in the ΔpkaA strain (see Table S1 in the supplemental material). For comparison, Franck et al. (35) found 485 peptides differentially phosphorylated in Magnaporthe oryzae in a comparison between the wild type (WT) and a cpkA (PKA) deletion strain. Among the 847 peptides identified here, only 5 phosphopeptides (corresponding to 5 proteins) were found exclusively in the ΔpkaA strain. These five phosphopeptides included AN8281, AN10591, AN8035.2, and AN7474 (uncharacterized proteins) and AN1948 (SPA10), which is a septal pore-associated protein (36). In contrast, 842 identified phosphopeptides (from 554 proteins) were found exclusively in the control strain.
Because the presence or absence of PKA is expected to affect the transcription of a large number of genes, we used RNA-seq analysis to select for proteins within our phosphoproteomic data set that were differentially phosphorylated but transcribed at similar levels. This approach increases the probability that a differentially phosphorylated peptide was not an artifact of different protein expression levels. When only similarly transcribed genes (0.5 Ͻ ΔpkaA/WT Ͻ 2.0) are considered, we find 337 unique phosphosites from 160 proteins that are differentially phosphorylated (Fig. 1A).
Enriched motifs among the set of PKA targets identified by phosphoproteomics were identified using motif X (Fig. 1B)  and ΔpkaA strains that were equally transcribed. (B) Ten most prevalent phosphorylation motifs (motif X) from differentially phosphorylated peptides. (C) Gene Ontology (GO) categories for proteins identified for biological function.
Novel Mode of CreA Regulation ® by a mitogen-activated protein kinase (MAPK) and ERK (40). We note that [SxxxSP] (motif 5), [RxxSxxS] (motif 6), [PS] (motif 9), and [SxP] (motif 10) appear to be novel as they have not been previously identified in the literature (37). We also note that the kinases predicted to phosphorylate the known motifs all play significant roles in the PKA pathway, acting as important mediators between active PKA and the cellular response to this activation.
All similarly transcribed and differentially phosphorylated proteins were categorized according to Gene Ontology (GO) for Biological Process (Fig. 1C). Consistent with previous findings, PKA impacts a broad range of proteins that are involved in a number of different biological processes. We find that PKA has a strong impact on a number of other kinases, as 28 were found to be differentially phosphorylated (Table S2). Networ-KIN 3.0 was used to predict which kinase phosphorylates each phosphosite. Of these 28 kinases, 4 were predicted to be phosphorylated directly by PKA, including a hexokinase (AN7459), a Pom1 kinase homologue pomA (AN7678), the regulatory subunit of PKA (AN4987), and a choline kinase (AN0929). Of these, the first three (AN7459, AN7678, and AN4987) have phosphorylation motifs matching motif 4 in Fig. 1B. However, AN0929 was phosphorylated at a different motif (RTVS*).
To identify additional direct PKA substrates, we searched for all differentially phosphorylated peptides that contained the classical PKA phosphorylation motif R/K-R/K-x-S and found 71 that matched it. GO categories for these phosphopeptides show that they are from proteins primarily related to cellular process (Fig. 2). Among these proteins, only one is a known PKA target, StuA. This implies that the remaining 70 are potentially new PKA substrates, none of which have been previously identified in the literature (Table S3). We assume the phosphopeptides with other motifs were phosphorylated indirectly, by another kinase, that itself may be a target of PKA.
One of the 70 proteins was a hexokinase, AN7459. A. nidulans contains two confirmed hexokinase-encoding genes (41,42); AN7459 (hxkA), and AN2638 (hxkB). The AN4255 (hxkC) and AN2180 (xprF) genes encode atypical noncatalytic hexokinases that regulate protease production (40,43). In our results, AN7459 was found to be equally transcribed and phosphorylated only in the WT, and AN4255 was differentially transcribed (three times less in the PKA Ϫ strain). AN2638 and AN2180 were equally  transcribed and not found to be phosphorylated in this work. de Assis et al. (29) reported a reduction in hexokinase activity in a pkaA knockout strain of A. nidulans compared to the WT, showing that PKA controls the activity of hexokinase. However, the mechanism by which hexokinase control occurs was not described. Therefore, according to the results obtained in this work, hexokinase activity could be controlled transcriptionally or by phosphorylation/dephosphorylation of HxkA or possibly HxkC.
Among the proteins that were equally transcribed, 19 transcription factors (TFs) were observed to be differentially phosphorylated ( Table 1). One of the TFs identified as being phosphorylated in the presence of PKA is PacC. PacC is a zinc finger transcription factor, and in A. nidulans, it mediates regulation of both acid-and alkalineexpressed genes by ambient pH. In Cryptococcus neoformans, Rim101/PacC is responsible for adaptive response to changes in pH and has been shown to share downstream targets with PKA (44). Previous work and our data strongly suggest that PacC is an indirect target of PKA. Among the 19 transcription factors identified here, two (AreA [AN8667] amd AreB [AN6221]) are involved with regulation of nitrogen metabolism. Since PKA is known to affect different aspects of cellular growth and to be involved in the regulation of nitrogen metabolism (45), our results imply that these TFs are downstream effectors of the PKA pathway and further explain the broad action of PKA activation.
To validate our phosphoproteomic data and provide more information regarding substrates that may be direct targets of PKA, we used anti-GFP antibodies to carry out an immunoprecipitation of PKA. In this experiment, the proteins directly interacting with PKA were coprecipitated and identified via MS/MS. Among a total of 217 proteins identified (Table S4), 32 proteins were common in both phosphoproteomic and immunoprecipitation experiments (Table S5). While coimmunoprecipitation is often used to study protein-protein interactions, it frequently presents a limitation with contaminant proteins being coprecipitated as well (46), but it does provide insight regarding potentially direct PKA substrates.
To investigate PKA targets in Candida albicans, Cao et al. mutated the two PKA catalytic subunits and used differential phosphoproteomic analysis to identify proteins in which phosphorylation was affected by the mutations (43). Identified proteins were subjected to in silico analysis to determine potentially direct substrates of PKA. To do this, they analyzed the phosphorylated motif (R/K-R/K-x-S/T) and found 148 proteins that are likely direct PKA substrates. Among the proteins listed as potential PKA targets, there were five orthologs also identified in this study. The proteins identified in common are a kinase (AN10515/PRK1), a protein related to biosynthesis of pyrimidine (AN0565/pyrABCN), an ATP-dependent RNA helicase subunit 2 (AN8722/SUB2), a protein repressed during the mating process (AN0446/PSP1), and a protein required for cell wall chitin distribution, morphology, and hyphal growth (AN0979/BNI4) (43). Transcriptomic analysis. We used an RNA-seq approach to assess differential gene expression in the pkaA deletion strain (ΔpkaA) and its isogenic parent. Because PKA is involved in many processes, deletion of pkaA affected gene expression very broadly. RNA-seq results showed that 125 TFs were differentially transcribed. A total of 44 TFs were downregulated and 81 were upregulated in the absence of PKA (Table S6). Of the TFs that were downregulated, CreA, an important carbon catabolite repressor in A. nidulans (14), had a 3.6-fold-lower expression level in the knockout strain. Since CreA represses transcription of cellulases and xylanases, we observed that the levels of expression of several genes related to carbohydrate metabolism change significantly [Ϫ1.0 Ͻ log 2 (ΔpkaA/wt) Ͻ 1.0] in the PKA deletion strain compared to the WT (Fig. 3). Our data also imply that PKA regulation occurs by both differential phosphorylation of TFs and altered TF transcription levels.
CreA regulation. Although our results showed numerous proteins that were equally transcribed and differentially phosphorylated, most of the proteins were uncharacterized or of unknown function. Therefore, we selected CreA, a transcription factor that has been well characterized in terms of the genes it regulates, the conditions in which it is active, and its functional localization, to study the effects of the phosphosites identified in our results. In the ΔpkaA strain, CreA was shown to be regulated both transcriptionally (i.e., 3.6-fold-lower expression in the ΔpkaA strain) and posttranslationally (S319 phosphorylated in the wild type and not in the ΔpkaA mutant). We note that in a previous A. nidulans study, we created a TAP-tagged version of CreA (47). We then grew this CreA:TAP strain in xylan, transferred it to glucose for various lengths of time (i.e., 5, 10, 15, and 30 min), and then immunoprecipitated CreA. We saw no physical interaction between CreA and the catalytic subunit of PKA (47). This evidence and the fact that we did not observe CreA when we immunoprecipitated PkaA:GFP (described above) strongly imply that CreA is not physically interacting with PkaA nor is it interacting in a very transient and/or weak manner. Further, a motif-based analysis (NetworKIN 3.0) predicts S319 to be a direct substrate of Stk22 kinase (AN5728). These To test the hypothesis that phosphorylation of S319 plays a critical role in the ability of CreA to mediate CCR, we mutated CreA residue S319 (Fig. 4A), identified in our phosphoproteomic analysis, to an alanine (S319A) and fused this construct to GFP (CreA S319A ::GFP). Subsequently, we investigated the expression of cellulase activity (Fig. 4B) and CreA::GFP nuclear localization ( Fig. 4C and D; see also Fig. S1 in the supplemental material) in this strain.
In A. nidulans, CreA consists of two C2H2-type zinc fingers (necessary for DNA binding), an alanine-rich region, an acidic region, and a conserved region, followed by a region that is required for the repression mediated by CreA (Fig. 4A) (28,48). Only the repressing region has been functionally characterized, and few studies have investigated CreA posttranslational regulation (28). As the main repressor in CCR, CreA inhibits the expression of numerous hydrolytic enzymes. In T. reesei, dephosphorylation of the CreA homologue Cre1 at residue S241 is necessary for binding to target DNA and repression of transcription (30). In another study, Alam et al. (49) have identified two peptides that were phosphorylated in CreA when A. nidulans was cultivated under repressing conditions and not phosphorylated under derepressing conditions. Among the two phosphopeptides, they found the sites S289, S312, and S319 to be phosphor- Novel Mode of CreA Regulation ® ylated only in the presence of glucose. To characterize the mutant CreA S319A strain phenotype, we assayed for cellulase activity (Fig. 4B), xylanase activity (Fig. S2B), radial growth in carboxymethyl cellulose, allyl alcohol, and 2-deoxy-D-glucose (Fig. S2). In all cases, we find no phenotypic difference between the wild-type CreA and CreA S319A strains. In accordance with a previous study (29), the pkaA deletion strain showed a significant increase in cellulase activity which was expected since transcription of several glucanases was higher in the deletion strain (Fig. 3). However, there was no difference in activity in the CreA S319A strain (Fig. 4B). CreA is known to display nuclear localization in the presence of glucose (33), a repressive condition. However, upon derepression, CreA is dispersed throughout the cytoplasm (33). We found that localization of CreA S319A was unaffected during growth in glucose (Fig. 4C). However, when the fungus was cultivated in derepressing conditions (i.e., ethanol) for 16 h and then switched to glucose, CreA S319A ::GFP did not migrate to the nucleus at the same rate as the CreA::GFP in the control (WT) strain ( Fig. 4D and Fig. S1). We used a simple, saturation, kinetic model to estimate saturation constants for both CreA::GFP and CreA S319A ::GFP (equation 1).
where Nuclear GFP (%) is the fraction of nuclei containing GFP (i.e., CreA), Y max is assumed to be 100%, and K m is the saturation constant. Values for K m for the wild-type control (CreA::GFP) and CreA phosphonull mutant (CreA S319A ::GFP) are 17.8 and 4.2 s, respectively, implying that the phosphonull CreA takes approximately 4 times longer to enter the nucleus, and thus, the CreA S319 phosphosite appears important for efficient CreA nuclear import. However, the mutation does not completely prevent CreA from travelling to the nucleus. Previous work has shown that the deletion of the conserved, repressor, and alanine-rich subdomains affect CreA localization (28). The phosphosite investigated here, S319, is located very close to the conserved domain and is also proximal to the repressing domain. In one study, the deletion of the conserved subdomain abolished growth on cellulose, which was not observed for CreA S319A (28). In addition, growth in glucose or cellulose was not affected by deletion of the repressor domain (28). Although the mechanism of CreA nuclear import remains incompletely understood, our results suggest that S319 phosphorylation regulates the dynamics of CreA nuclear transport.
To determine whether expression levels of GFP-tagged, phosphonull CreA S319A were similar to those of GFP-tagged CreA we grew fungi in minimal medium with glucose, isolated cellular proteins, and performed Western blotting (Fig. S3). We find a significant increase in expression of GFP-tagged phosphonull CreA S319A but similar stabilities of CreA in both strains. While it is not clear why there was an increase in expression, both strains showed similar cellulase activity (Fig. 4B), implying that this increase did not impact cellular phenotype.
Conclusion. In this work, we were able to comprehensively identify both direct and indirect targets of PKA phosphorylation in A. nidulans, as well as genes whose expression depends on the presence of PKA. In the PKA deletion strain, transcription factor CreA was found to be down phosphorylated at S319. Recently, we have immunoprecipitated A. nidulans CreA:TAP under repressing and derepressing conditions using glucose or xylan, respectively, as the single carbon source (47). In both conditions, we have not seen PkaA interacting physically with CreA (47). Together with the results shown here, these data strongly indicate that PkaA is not physically interacting with CreA or that this interaction is very transient and/or weak. Thus, it is quite likely that CreA S319 is phosphorylated by another kinase regulated by PkaA and thus is an indirect target of PKA. Functional analysis of CreA S319 revealed its importance for timely nuclear import. Besides providing a broad overview of the PKA regulatory network in a model filamentous fungus, our results will also facilitate detailed functional investigation of PKA phosphorylation events and their role(s) in growth and development.
Extraction and digestion of intracellular proteins for phosphoproteomics. Mycelia from FGSC A1405 (WT) and FGSC A1358 (PKA deletion) were harvested after 20 h of growth (during logarithmic growth phase; growth curves shown in Fig. S4 in the supplemental material) by passing the culture through a Buchner funnel and immediately freezing mycelia in liquid nitrogen. Frozen mycelia were ground into powder by using a mortar and pestle, and RIPA buffer (Thermo Fisher Scientific) was used to extract protein with a 15-min incubation at 4°C and constant agitation. The mixture was centrifuged at 10,000 ϫ g for 20 min, and the supernatant was collected and processed for protein concentration by BCA protein assay (Thermo-Fisher Scientific Inc., Rockford, IL). Five hundred micrograms of protein was purified by 3-kDa-molecularweight cutoff ultrafiltration (Centriprep YM-3; Millipore, Billerica, MA) followed by trichloroacetic acid (TCA) precipitation. TCA was added to protein extract at 10% final concentration (wt/vol, made fresh), vortexed for 15 s, and placed on ice for a minimum of 20 min. The samples were centrifuged at 14,000 ϫ g and 4°C for 15 min, and the supernatant was removed. The pellet was washed with cold acetone and then centrifuged at 14,000 ϫ g and 4°C for 10 min. The acetone was discarded, and the samples were air dried. Protein (500 g) was reduced with Tris(2-carboxyethyl)phosphine, alkylated with iodoacetamide, and digested with sequencing-grade trypsin (Promega, Madison, WI). Peptides were lyophilized, and phosphopeptides were enriched using Pierce TiO2 Phosphopeptide Enrichment and Clean-up kit (Thermo Fisher Scientific) according to the manufacturer's instructions. Samples were then lyophilized. . The gradient elution started at 5% mobile phase B, increased to 35% at 60 min and then 80% at 65 min, and was held at 80% for 5 min before a 25-min reequilibration at 5%. The eluting peptides were interrogated with an Orbitrap Fusion mass spectrometer running a data-dependent LC-MS/MS method with the Top Speed decisions selection. FTMS1 spectra were collected using the following parameters: scan range of 350 to 1,800 m/z, resolving power of 120K, AGC target 4E5, and maximum injection time of 50 ms. ITMS2 spectra were collected using the following parameters: rapid scan rate, CID NCE 35, 1.6 m/z isolation window, AGC target 1E4, and maximum injection time of 50 ms. MS2 precursors were selected for a 3-s cycle. Precursors with an assigned monoisotopic m/z and a charge state of 2 to 7 were interrogated. Precursors were filtered using a 60-s dynamic exclusion window.
Data analysis for phosphoproteomics. Raw mass spectrometric data were analyzed in the MaxQuant environment (50) v.1.5.0.0 and employed Andromeda for database search (51). The MS/MS spectra were matched against the Aspergillus nidulans Uniprot FASTA database. Enzyme specificity was set to trypsin, and the search included cysteine carbamidomethylation as a fixed modification and N-acetylation of protein, oxidation of methionine, and/or phosphorylation of Ser, Thr, and Tyr residues (STY) as variable modifications. Up to two missed cleavages were allowed for protease digestion, and peptides had to be fully tryptic. The maximum peptide and site false-discovery rates were specified as 0.01. We used Student test for testing differences between the means of two conditions. Two-tailed hypothesis test at the 0.05 level of significance was applied.
RNA-seq analysis. Samples for transcriptomics were obtained from the same cultures used for phosphoproteomic samples. Harvested cells were immediately frozen, crushed, and prepped for RNA-seq analysis. For RNA-seq analysis, sequencing libraries were constructed from RNA samples using TruSeq RNA Sample Preparation kit v2 (Illumina, Inc., San Diego, CA), and 100-bp reads were generated on the Illumina Genome Analyzer IIx platform (Illumina, Inc., San Diego, CA). RNA-seq output data were checked for base quality and aligned to the reference genome of Aspergillus nidulans (www.aspgd.org/) using bowtie2, and the read counts were normalized using HT-Seq. Differential expression analysis was done using DESeq software (52). Genes were called differentially expressed based on log 2 fold change, and significance is estimated for P values less than 0.01.
GFP tag protein purification. The GFP-tagged protein complex was purified from crude protein extracts incubating 4 g ground mycelia in 4 ml B250 buffer described above. The total protein lysate was incubated with 40 l GFP-Trap magnetic beads (Chromotek) for 4 h at 4°C. Beads were harvested using magnetics hack and washed three times using B250 buffer without DTT, being added only in the last wash step. The precipitated complex was digested using sequencing-grade modified trypsin (catalog no. V5117; Promega), and the desalinization step was made using Zip-Tip (catalog no. ZTC18S096; Millipore) according to the manufacturer's instructions.
LC-MS/MS protein identification for pka::GFP pulldown. Digested peptides were separated using reversed-phase liquid chromatography with an RSLCnano Ultimate 3000 system (Thermo Scientific) followed by mass identification with an Orbitrap Velos Pro mass spectrometer (Thermo Scientific). Chromatographically separated peptides were on-line ionized by nanoelectrospray ionization (nESI) using the Nanospray Flex Ion Source (Thermo Scientific) at 2.4 kV and continuously transferred into the mass spectrometer. Full scans within m/z of 300 to 1,850 were recorded by the Orbitrap-FT analyzer at a resolution of 30.000 (using m/z 445.120025 as lock mass) with parallel data-dependent top 10 MS2fragmentation in the LTQ Velos Pro linear ion trap. LC-MS method programming and data acquisition were performed with the software XCalibur 2.2 (Thermo Scientific) and method/raw data validation with the program RawMeat 2.1 (Vast Scientific).
MS/MS2 data processing for protein analysis and identification was done with either MaxQuant quantitative proteomic software in conjunction with Perseus software for statistical analysis or the Proteome Discoverer 1.3 (PD) (Thermo Scientific) and the Discoverer Daemon 1.3 (Thermo Scientific) software using the Sequest (and/or Mascot) peptide analysis algorithm(s) and organism-specific taxondefined protein databases extended by the most common contaminants.
The interactions identified were submitted to filtering subtracting unspecific interactions from AGB551 as the genetic background strain grown under the same experimental conditions. MS interactions from duplicates for each time point with at least two unique peptides identified per protein were applied during the analysis; PkaA was identified in all pulldown assays as an internal control for tagged protein precipitation.
GFP-Trap immunoprecipitation. Total protein extracts from mycelia were obtained by extraction from ground mycelia with B250 extraction buffer (250 mM NaCl, 100 mM Tris-HCl [pH 7.5], 10% glycerol, 1 mM EDTA, and 0.1% NP-40) supplemented with 1.5 ml/liter of 1 M DTT, 1 pill/10 ml of the Completemini Protease Inhibitor Cocktail EDTA-free (Roche), 3 ml/liter of 0.5 M benzamidine, 100 l/10 ml of phosphatase inhibitors P0044 (Sigma), and 10 ml/liter of 100 mM PMSF. GFP-Trap beads (40 l) were equilibrated with 500 l of extraction buffer and then centrifuged at 3,000 rpm at 4°C for 1 min; all the next wash steps were done under the same conditions. Five milligrams of total protein extract was incubated with GFP-Trap beads for 4 h at 4°C in a rotator shaker; after the incubation, the beads were washed two times using extraction buffer without DTT and one more additional wash with DTT. The supernatant was removed, and the GFP-Trap beads were incubated with 40 l of 1ϫ SDS sample buffer and boiled at 95°C for 5 min. The proteins (total) released from the beads were loaded into SDSpolyacrylamide gel (12%) for Western blotting as previously described (47).
Fluorescence microscopy for nuclear localization. For microscopy experiments, spores were grown on coverslips for image analysis. Briefly, coverslips bearing freshly harvested spores were inoculated in a small petri dish filled with 3 ml minimal medium supplemented with 1% Avicel or 1% glucose for 16 h at 22°C. When mycelia were grown in minimal medium supplemented with 1% Avicel or 1% ethanol, 1% glucose was added after 16 h. Images were taken after 2.5, 5, 10, 15, 30, 60, 80, and 100 min. Mycelia mounted on coverslips were washed with phosphate-buffered saline (PBS) (140 mM NaCl, 2 mM KCl, 10 mM NaHPO 4 , 1.8 mM KH 2 PO 4 [pH 7.4]). The mycelia were then stained with 100 ng/ml Hoechst 33258 (Molecular Probes) for 2 min. The mycelia were washed again in water and examined using a Zeiss epifluorescence microscope with excitation wavelengths of 359 and 498 nm and emission wavelengths of 461 and 516 nm for Hoechst and GFP, respectively. Phase-contrast bright-field and fluorescent images were captured with an AxioCam camera (Carl Zeiss) and processed using AxioVision software version 3.1. Representative images (Fig. S1) show nuclear localization. Brightness and contrast have been altered identically in all images to better show CreA::GFP localization. To assess the rate of nuclear localization, data were subjected to nonlinear regression (equation 1) using PolyMath software (CACHE Corp.) to determine values of K m .
Cellulolytic activity. Enzymatic activity was measured by colorimetric assay using Avicel as a substrate. The reducing sugars were determined by the Miller procedure using glucose as the control (53). The reaction mixture (0.05 ml substrate [1%, wt/vol] in 50 mM sodium acetate buffer [pH 5.0] and 0.05 ml enzyme solution) was incubated at 50°C for 30 min. The reaction was stopped by adding 0.1 ml of dinitrosalicylic acid (DNS) and immediately boiling for 5 min. Quantification of the reducing sugars released as a result of enzyme activity was estimated by A 540 measurements, where one unit of enzymatic activity was defined as the amount of enzyme that produced 1 mol min Ϫ1 of reducing sugars.
Site-directed mutagenesis of CreA. The substitution of CreA S319 with alanine was introduced by site-directed mutagenesis, using oligonucleotides containing the mutation. The DNA cassette was constructed containing 1 kb of the upstream creA region, creA mutated, pyrG, and 1 kb of the downstream creA region. The plasmid pCreA/GFP containing the whole cassette was built using plasmid pUC19 as a backbone via Gibson assembly (54). The mutated CreA (S319A) was constructed by inverse PCR (55) of the pCreA/GFP plasmid, generating the construct pCreAS319A/GFP containing the mutated CreA gene. After fully sequencing pCreAS319A/GFP, the CreAS319A/GFP/pyrG cassette was PCR ampli-fied and used to transform A. nidulans TNO2A3. Genomic DNA extraction was conducted for six transformants and submitted for sequencing to verify the absence of any additional base changes introduced into creA as a result of transformation.
Data availability. The RNA-seq data from this publication have been deposited into the GEO database (https://www.ncbi.nlm.nih.gov/geo/) and assigned the identifier or accession no. GSE116579.

SUPPLEMENTAL MATERIAL
Supplemental material for this article may be found at https://doi.org/10.1128/mBio .02825-18.   Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.
We declare that we have no conflict of interest.