- Open Access
Identification and characterization of the glucose dual-affinity transport system in Neurospora crassa: pleiotropic roles in nutrient transport, signaling, and carbon catabolite repression
Biotechnology for Biofuelsvolume 10, Article number: 17 (2017)
The glucose dual-affinity transport system (low- and high-affinity) is a conserved strategy used by microorganisms to cope with natural fluctuations in nutrient availability in the environment. The glucose-sensing and uptake processes are believed to be tightly associated with cellulase expression regulation in cellulolytic fungi. However, both the identities and functions of the major molecular components of this evolutionarily conserved system in filamentous fungi remain elusive. Here, we systematically identified and characterized the components of the glucose dual-affinity transport system in the model fungus Neurospora crassa.
Using RNA sequencing coupled with functional transport analyses, we assigned GLT-1 (K m = 18.42 ± 3.38 mM) and HGT-1/-2 (K m = 16.13 ± 0.95 and 98.97 ± 22.02 µM) to the low- and high-affinity glucose transport systems, respectively. The high-affinity transporters hgt-1/-2 complemented a moderate growth defect under high glucose when glt-1 was deleted. Simultaneous deletion of hgt-1/-2 led to extensive derepression of genes for plant cell wall deconstruction in cells grown on cellulose. The suppression by HGT-1/-2 was connected to both carbon catabolite repression (CCR) and the cyclic adenosine monophosphate-protein kinase A pathway. Alteration of a residue conserved across taxa in hexose transporters resulted in a loss of glucose-transporting function, whereas CCR signal transduction was retained, indicating dual functions for HGT-1/-2 as “transceptors.”
In this study, GLT-1 and HGT-1/-2 were identified as the key components of the glucose dual-affinity transport system, which plays diverse roles in glucose transport and carbon metabolism. Given the wide conservation of the glucose dual-affinity transport system across fungal species, the identification of its components and their pleiotropic roles in this study shed important new light on the molecular basis of nutrient transport, signaling, and plant cell wall degradation in fungi.
Organisms dynamically interact with their biotic and abiotic environments by adjusting cell transcriptional and post-transcriptional programs in response to external nutrient (like sugar) changes and stress stimuli. In theory, nutrient transport, sensing, and downstream signal transduction are critical for the coordination of extracellular nutrient availability with internal metabolism, development, and cell survival [1–3]. Glucose acts as a primary carbon source as well as a pivotal signal that triggers a cellular regulatory network influencing sugar transporter expression, carbon catabolism, and biomass accumulation. In particular, glucose sensing and uptake are critical for cellulase expression regulation via carbon catabolite repression (CCR) [4, 5]; therefore, elucidating the molecular basis of glucose transport is critical for plant biomass deconstruction and bio-based chemical and fuel production. The glucose transport system of Saccharomyces cerevisiae is the best documented system among fungal species, and consists of 20 different hexose carriers belonging to the major facilitator superfamily (MFS) [3, 6, 7]. S. cerevisiae possesses a dual-affinity transport system for glucose uptake that is coordinated with these hexose transporters composed of a high-affinity system (K m = 1–2 mM) and a low-affinity system (K m = 15–20 mM) . In the filamentous fungi Neurospora crassa, glucose uptake also behaves in a dual-affinity fashion, with K m values of 10–50 µM (system II) and 8–20 mM (system I) [8–10]. However, the genes of this dual system have still not been uncovered and their functions thus remain to be dissected.
Among the 20 different hexose carries in S. cerevisiae, two transporter-like glucose receptors Rgt2p and Snf3p, which sense low and high levels of external glucose, respectively, mediate a glucose signaling pathway [11, 12]. The recently identified high-affinity hexose transporter Hxt1, the homolog of Rgt2p/Snf3p in Ustilago maydis, also functions in glucose signaling ; this indicates the potential of hexose transporters to act as “transceptors” in fungi. Transceptors of other nutrients, such as nitrate, amino acids, phosphate, and cellobiose, have been found in eukaryotes [14–17], but glucose transceptors, to our knowledge, have not yet been definitively characterized, although Hxt1 in Ustilago maydis and GLUT2 in humans might possess this dual function [13, 18].
Another branch of glucose-sensing and signal modulation in fungi is mediated by the cAMP-protein kinase A (PKA) signaling pathway. In S. cerevisiae, binding of glucose to Gpr1p activates the downstream Gα protein Gpa2p, leading to stimulation of the adenylate cyclase Cyr1p, which increases cAMP levels, thereby affecting PKA activity . The PKA pathway regulates a wide range of processes in fungi, including metabolism, cell growth, circadian rhythms, germination, and conidiation (for details, see reviews [19, 20]).
In glucose signaling in N. crassa, the ortholog of yeast Gpr1p is GPR-4 (NCU06312), which has been identified as a carbon source receptor . Ligand binding to GPR-4 can stimulate the downstream Gα protein GNA-1 (NCU06493), leading to an increase in the level of cAMP produced by the activated adenylate cyclase CR-1 (NCU08377) . External glucose sensing is associated with a Rgt2p/Snf3p ortholog RCO-3 (NCU02582)-mediated pathway, in which RCO-3 appears to act as a non-transporting glucose sensor . Mutation of rco-3 leads to complete dysfunction of the low-affinity transport system and partial impairment of the high-affinity system .
In this study, we characterized the glucose dual-affinity transport system in the model fungus N. crassa. One low-affinity glucose transporter, GLT-1 (NCU01633), and two high-affinity transporters, HGT-1/-2 (NCU10021 and NCU04963), were identified as the major components of systems I and II, respectively. Simultaneous deletion of hgt-1/-2 (strain Δ2hgt) resulted in a notable increase in cellulolytic enzyme production when N. crassa was grown on cellulose. A group of carbohydrate-active enzymes (CAZys), nearly all glycolytic enzymes and asexual sporulation genes, were differentially expressed in the Δ2hgt mutant according to RNA-Seq. Based on analyses of an array of mutants and point mutations within HGT-1/-2, we hypothesized that the glucose dual-affinity transport system comprising GLT-1 and HGT-1/-2 in N. crassa was involved in glucose transportation, sensing, and downstream signaling cascades. This is the first time the glucose dual-affinity transport system of N. crassa has been systematically elucidated at the gene level. Our findings significantly improve our understanding of glucose uptake and signaling in fungi, and shed new light on plant cell wall deconstruction for cellulosic biorefinery by filamentous fungi.
Genome-wide analysis of the transcriptional responses of fungal cells to a glucose gradient
Transmembrane proteins are tightly regulated to coordinate environmental nutrient changes with intracellular metabolism and cell proliferation. Previous studies have suggested that the predicted low-affinity glucose transport system (system I) in N. crassa [8, 23, 24] is induced at high glucose levels, while glucose-limited or carbon-deprived conditions trigger de novo protein synthesis of the predicted high-affinity system (system II), which is fully activated within 1–2 h [9, 23]. These findings indicate that de novo mRNA synthesis of system II components occurs when mycelia are exposed to carbon-limited conditions. To obtain a broad view of the mode of expression and to define the functional molecular elements of the dual-affinity transport system, we conducted high-throughput sequencing (RNA-Seq) of wild-type (FGSC2489) mycelia exposed to a gradient of glucose (0, 0.05, 0.5, 2.0, 10.0%) for 1 or 2 h. To our knowledge, this is the first time transcriptional profiling analysis of a fungal response to a glucose gradient has been conducted. Sample-to-sample clustering  demonstrated the biological replicates were reliable, as evidenced by a high Spearman’s rho (>0.98, P value < 0.001) for all tested samples (Additional file 1: Figure S1a). Furthermore, the response observed under no-carbon conditions in this study was in good accordance with the published data  (Spearman’s rho = 0.96, P value < 0.001; Additional file 1: Figure S1B). Differential gene expression analysis (|GFOLD| > 1 and P value < 1 × 10−4, “Methods” section; genome-wide expression levels and differential expression analysis are described in Additional file 2: Table S1) revealed that the transcriptomic responses to 0.5% (27.8 mM), 2.0%, and 10.0% glucose were essentially similar (Table 1), indicating that a glucose threshold of 0.5% is adequate to support vegetative growth of N. crassa (note: all cited glucose percentages are w/v unless otherwise indicated). Although 96 genes with altered expression levels were detected on 10.0% glucose compared with 2.0%, no specific functional categories relevant to this mild osmotic stress condition (with reference to typical osmotic stress treatments of 3.0–8.0% NaCl in S. cerevisiae ) were enriched according to FunCat analysis  (Additional file 3: Table S2). In contrast, several categories associated with 66 genes robustly induced by 0.05% glucose for 1 h (compared with 2.0% glucose) were overrepresented (Additional file 3: Table S2). Among these categories were several reflecting cellular efforts to maintain homeostasis, including carbohydrate and fatty acid metabolism, glycolysis/gluconeogenesis, pentose utilization, sensing of external changes, and especially nutrient and ion transport. In the last category, 15 genes were annotated as transmembrane transporters based on TransportDB  (Additional file 3: Table S2). Additionally, 66 genes were predominant among the genes commonly upregulated in both 0.05% glucose for 2 h and carbon-free conditions (Fig. 1a; Additional file 3: Table S2), suggesting they belong to pioneer blocks in the response of the cellular architecture to external glucose depletion. This inference was especially supported by the strikingly decreased expression of the CCR regulator cre-1 and its corepressor rcm-1 (Additional file 2: Table S1), whereas the ortholog of Snf1, prk-10 (NCU04566), was robustly upregulated under severe glucose-depleted conditions (Additional file 2: Table S1). In A. nidulans, activated SnfA can phosphorylate CreA, leading to CreA dissociation from the CreA-Ssn6-RcoA complex and consequent attenuation of CCR . In addition to the attenuation of CCR, the cAMP-PKA glucose induction pathway seemed to have been downregulated, as both parallel cAMP synthetic processes were transcriptionally repressed. In particular, gna-1 (NCU06493) and two Ras proto-oncogenes, ras-1/-2 (NCU08823 and NCU06111), were 3- to tenfold downregulated under glucose-deprived or carbon-free conditions (Additional file 2: Table S1). In addition to cAMP-PKA inactivation, the transcriptomic downregulation also included reduced expression of DNA, RNA, and ribosomal biosynthetic genes, indicating that the fungal cells were struggling with carbon starvation (Fig. 1b; Additional file 3: Table S2). Two key autophagy genes, NCU00188 and NCU01545, encoding orthologs of ATG1/AtgA and ATG8/AtgH in S. cerevisiae and A. nidulans, respectively , showed increased expression levels under the no-carbon condition in comparison with 2.0% glucose treatment (Additional file 2: Table S1), whereas vib-1 (NCU03725), an ortholog of the crucial autolysis regulator XprG in A. nidulans , was dramatically downregulated (Additional file 2: Table S1). These data suggest complex coordination between autophagy and autolysis in N. crassa to cope with carbon starvation as these two processes might be activated independently .
Identification and characterization of the glucose dual-affinity transport system revealed it is coordinated with glucose homeostasis
One of the most sensitive responses during external glucose depletion was related to nutrient assimilation (Additional file 3: Table S2). Previous work has shown that 39 putative sugar transporters are present in the genome of N. crassa [33, 34]. Among these transporters, 26 showed robust expression levels (i.e., reads per kilobase per million mapped reads [RPKM] > 20) under at least one tested treatment condition (Fig. 2a). Clustering analysis grouped eight identified MFS transporters into cluster 1 (C1), comprising one glucose sensor, rco-3 ; two cellodextrin transporters, cdt-1/-2 ; three pentose transporters, lat-1 (also a weak glucose transporter), xat-1, and xyt-1 [33, 36]; the galacturonic acid transporter gat-1 ; and a cellobionic acid transporter, cbt-1 . C1 showed very little expression under 0.5% (or higher) glucose conditions; 17 out of 21 genes in this cluster displayed a weak response to external glucose changes, except cdt-1, xyt-1, NCU05897 (a putative fucose transporter), and NCU09287 (Fig. 2a). In contrast, C2 expression was more sensitive to glucose fluctuations. Included in C2 were the glucose transporter glt-1 , the xylose-specific transporter NCU00821 , the putative sucrose transporter NCU00450, the high-affinity glucose transporter hgt-1 , and NCU04963 (designated as high-affinity glucose transporter hgt-2; see the next paragraph). The glucose transporter glt-1 was the only transporter that was constantly expressed on high glucose (Fig. 2), suggesting a predominant role in the low-affinity transport system. In contrast, hgt-1/-2 were strongly upregulated under low-glucose or glucose-deprived conditions vs. 2.0% glucose. Furthermore, the mRNA expression levels of these two genes accounted for nearly 80% of the total expression of all sugar transporters under the no-carbon condition (Fig. 2b, c). This observed dynamic expression pattern demonstrated that hgt-1/-2 were sensitively triggered by glucose depletion (Fig. 2c). Intriguingly, glt-1 was predominantly expressed at high levels of glucose, but transiently strong upregulation was also observed under a moderately low-glucose concentration (0.05%, 1 h) (Fig. 2c), which suggests that its expression might be regulated by glucose in a dosage-controlled, bidirectional way.
Because the glt-1 mRNA was predominant on high glucose, while hgt-1/-2 mRNAs dominated glucose-limited conditions (Fig. 2), these three transporters were selected as potential candidate genes for systems I and II, respectively. GLT-1 has been characterized as a glucose transporter that can efficiently complement the growth of the S. cerevisiae null-hexose-transporter strain EBY.VW4000 . As anticipated, HGT-1/-2 functioned as glucose transporters when heterologously expressed in EBY.VW4000 (Fig. 3a, b). Uptake analysis with radiolabeled glucose in S. cerevisiae revealed a high K m value for GLT-1 (18.42 ± 3.38 mM), while the K m values for HGT-1 and HGT-2 were three orders of magnitude lower (16.13 ± 0.95 and 98.97 ± 22.02 µM, respectively) (Fig. 3c). This result was in line with a previous observation for HGT-1 . The high K m value observed for GLT-1 was comparable to the apparent affinity of system I (K m = 8–20 mM), while the values of HGT-1/-2 were consistent with that of system II (10–50 µM). Summarizing the above data (Figs. 2, 3), GLT-1 was assigned to the low-affinity system (system I) and HGT-1/-2 to the high-affinity system (system II).
To dissect the physiological functions of these major components of the glucose dual-affinity transport system, we generated double- and triple-gene deletion mutants via sexual crosses (three progenies were selected for the double mutant and two for the triple mutant in subsequent analysis; “Methods” section). Double deletion of hgt-1/-2 resulting in the mutant Δhgt-1;Δhgt-2 (designated Δ2hgt hereafter in the text) led to a significant deficiency in the uptake rate on a low concentration of glucose (1.1 mM; Fig. 4a). Further deletion of glt-1 did not alter the transport capacity of the Δ2hgt mutant (Fig. 4a). This result was in agreement with the observation that the Δglt-1 strain had similar uptake activity to the WT (Fig. 4a). This similarity suggested that transport by GLT-1 (system I) was dispensable for glucose uptake under glucose-limited conditions, and that the high-affinity system (system II) was sufficient. When the putative sucrose transporter NCU00450 or NCU05897 (also strongly induced in low glucose; Fig. 2a) was also knocked out in the Δ2hgt background, no other significant difference was observed (Additional file 4: Figure S2). This result indicated that HGT-1/-2 were the major components of system II. In addition, the fact that the Δ2hgt mutant still exhibited detectable glucose uptake indicated that some other sugar transporters might have been responsible for limited glucose uptake when the two hgt genes were deleted (Fig. 4a). For example, LAT-1, which was slightly upregulated under no-carbon conditions (Fig. 2a), has been shown to have the ability to transport glucose in addition to its preferred substrate arabinose .
Intriguingly, deletion of glt-1 did not lead to any deficiencies when N. crassa mycelia were grown on high levels of glucose (Fig. 4b, c), whereas the triple mutant Δ2hgt;Δglt-1 (Fig. 4b), in which all the major players of systems I and II were deleted, showed growth lag. This result indicated that system II could complement reduced activity of system I (Fig. 4c). To verify this conclusion, we performed a quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) assay to monitor gene expression in the Δglt-1 strain treated with high glucose. As expected, hgt-1/-2 were dramatically upregulated (approximately 60-fold) in the presence of adequate glucose (0.5%; Fig. 4d). Some other putative sugar transporters, such as NCU06358 and xyt-1 (Fig. 4d), were synchronously induced, indicating they have redundant roles in glucose uptake. In contrast, no significant increase in glt-1 expression was found in the Δ2hgt strain at low levels of glucose (Fig. 4e).
Negative influence of HGT-1/-2 on derepression of genes for plant cell wall deconstruction
Transcriptional expression data from previous reports [26, 33, 42, 43] and this study revealed a distinct utilization preference for glucose transporter systems in N. crassa: system I (glt-1) is activated when an adequate amount of the preferred carbon source is present, whereas system II (hgt-1/-2) is expressed under low-glucose, lignocellulosic, or carbon-free conditions (Additional file 5: Figure S3). Previous studies have demonstrated that system II can be derepressed by carbon-limited conditions [9, 10, 23], which was verified by our transcriptome profiling data (Fig. 2). Considering that a synergistic carbon starvation response can be induced by Avicel or plant biomass [26, 43], the high expression levels of hgt-1/-2 observed in this study likely resulted from the lifting of CCR when mycelia were exposed to plant cell walls.
Given that hgt-1/-2 are strongly derepressed under cellulosic conditions, we wondered whether these two genes play roles in cellulose utilization in N. crassa. To answer this question, we inoculated a series of hgt-deleted mutants and the WT into liquid medium with 2.0% (w/v) Avicel as the sole carbon source. Surprisingly, the Δ2hgt strain showed nearly twofold higher protein production than the WT (Fig. 5a), while only slightly improved secretion rates were found in the single hgt-deleted mutants (Additional file 6: Figure S4). The triple mutant Δ2hgt;Δglt-1 displayed hyper-production similar to that of Δ2hgt, consistent with the glucose uptake observed in Fig. 4a, with similar biomass accumulation to the WT (Additional file 7: Figure S5). To verify that the two high-affinity glucose transporters (HGT-1/-2) but not the low-affinity one (GLT-1) affected lignocellulase expression, the three transporters driven by the promoter of ccg-1, which is highly expressed under cellulose- or carbon starvation conditions [26, 44] (Additional file 2: Table S1), were individually introduced into the Δ2hgt mutant (Fig. 5b). Both high-affinity transporters, HGT-1/-2, restored the lignocellulase production and activity of Δ2hgt to WT levels, but the low-affinity transporter GLT-1 did not (Fig. 5a, c; Additional file 8: Table S3). These data demonstrate that the derepressible high-affinity glucose transport system negatively affects cellulase expression in N. crassa.
Transcriptional analysis by qRT-PCR revealed that expression levels of the cellulase regulator clr-2 (NCU08042) and the tested carbohydrate-active enzymes (CAZy genes) were considerably higher in Δ2hgt and Δ2hgt;Δglt-1 compared with the WT (Fig. 6a, b). The largest differences from the WT (Fig. 6b) were observed 3 days after conidial inoculation, at around the time protein production started to deviate upwards in the Δ2hgt mutant (Fig. 5a).
To further elucidate the genome-wide expression changes in Δ2hgt, time-course comparative transcriptome profiling was performed via RNA-Seq. Correlation analysis showed reliable biological reproducibility as indicated by sample-to-sample clustering and a high Spearman’s rho (>0.96, P value < 0.001) (Additional file 9: Figure S6; gene expression levels and differential expression analysis are shown in Additional file 10: Table S4). Overall, 629 genes were robustly upregulated in the system II-deficient mutant compared with the WT at the three tested time points (Fig. 6c; Table 2). Hierarchical clustering analysis of the 629 genes revealed three clusters (Fig. 6c). The first cluster (C1) contained 237 genes whose expression levels decreased in both the WT and the Δ2hgt mutant with increasing culture time, but were generally higher in Δ2hgt compared with the WT at each tested time point (Fig. 6c). Most upregulated CAZy genes from the 629 gene set were overrepresented in C1 (57 out of 70, P = 6.65 × 10−5 according to the one-tailed Fisher’s exact test; Additional file 11: Table S5). Also included in this group was clr-2 (Additional file 11: Table S5). Although clr-1 did not meet the strict differential analysis criteria (“Methods” section) with respect to the hgt null mutant (GFOLD = 0.43), it displayed a significantly elevated expression level at 72 h (P = 6.88 × 10−19 and raw log2[fold change] = 1.42; Additional file 10: Table S4). In N. crassa, the regulators clr-1/-2 positively control the transcription of a number of lignocellulases, among which 126 genes are defined as belonging to the “Avicel regulon” . Fifty-seven such genes were found in the 629 gene set, 53 of which were enriched in C1. This finding suggests that the upregulated genes in C1 are at least partially regulated by clr-1/-2. Functional category enrichment analysis also supported this idea (Additional file 11: Table S5). The group C1 also contained 11 MFS transporters, including glt-1. This sugar transporter and another one in C2 (xyt-1) were synergistically transcribed to high levels in the Δ2hgt mutant (Additional file 10: Table S4).
In addition to many derepressed lignocellulases, the pentose phosphate pathway was overrepresented in C1 (P = 4.03 × 10−4; Additional file 11: Table S5). This observation supports the notion that CCR was lifted when hgt-1/-2 were deleted. In fact, the expression level of inv (NCU04265), an invertase that has been widely demonstrated to be a marker for CCR in N. crassa and yeast species [45–47], was consistently high in Δ2hgt (Additional file 10: Table S4). Furthermore, cre-1, a well-characterized regulator involved in CCR that globally suppresses cellulase expression in N. crassa , was downregulated in Δ2hgt (Fig. 6a). To confirm that CCR was lifted in the Δ2hgt mutant, the glucose analog 2-deoxy-glucose (2-DG, which cannot be catalyzed during glycolysis and is a drug often used for glucose repression analysis in filamentous fungi) was added to the Vogel’s growth medium containing Avicel as the sole carbon source. Supporting the above hypothesis, the Δ2hgt strain showed less sensitivity to 2-DG relative to the WT (Fig. 6d, e). Taken together, these results indicate that HGT-1/-2 contribute to the activation of CCR in the WT at a certain time point during the cellulolytic utilization process, including activation of cre-1 expression, which triggers general repression affecting clr-1/-2 and thus suppressing lignocellulase gene expression.
Lignocellulase gene repression by HGT-1/-2: functionally coordinated but not confined to CRE-1
As shown above, exposure to carbon-limited environments such as plant cell walls causes N. crassa to use the high-affinity transporter system to uptake cellulose-derived glucose. The observation that cre-1 expression was reduced during cellulose utilization (Fig. 6a) therefore suggests that the suppression of lignocellulolytic degradation enzymes mediated by HGT-1/-2 might be downstream-regulated by CRE-1. To test this idea, cre-1 was deleted or constantly expressed in the Δ2hgt mutant background. The rate of secreted protein production by the Δ2hgt;Δcre-1 and Δ2hgt;Δglt-1;Δcre-1 mutants on Avicel was dramatically increased compared with the Δ2hgt strain (Fig. 7a). This difference was even more obvious when the Avicel substrate was replaced with cellobiose, on which deletion of cre-1 resulted in protein production starting at least 3 days earlier and ultimately twofold greater accumulation than the background strain Δ2hgt (Additional file 12: Figure S7). Examination of submerged growth phenotypes and protein gels confirmed this observation (Fig. 7c, d). In accordance with these results, the consumption of Avicel by Δ2hgt;Δcre-1 was also faster than that of Δ2hgt and the WT (Fig. 7b). Previous work has demonstrated that the single cre-1 deletion mutant utilizes the Avicel substrate 2–3 days faster than the WT , implying an independent role for cre-1 in the Δ2hgt mutant during cellulose utilization. In the Δ2hgt strain (Additional file 13: Figure S8), however, constant activation of cre-1 expression via the ccg-1 promoter or the ribosomal protein 27 promoter from Magnaporthe grisea  suppressed extracellular protein production to a level comparable to that of the WT (Fig. 7a, d). This result implies that regulation of cellulase expression via CRE-1 is coordinated with HGT-modulated repression in N. crassa. All these data, combined with the observation of reduced cre-1 expression in the Δ2hgt strain (Fig. 7a), suggest that repression by HGT-1/-2 involves cross-talk with the pivotal CCR regulator CRE-1 and potentially some other unknown factor(s).
To verify our hypothesis, comparative genome-wide expression profiling was performed via RNA-Seq between the Δ2hgt and Δ2hgt;Δcre-1 mutants (Additional file 9: Figure S6; Additional file 10: Table S4). The expression levels of selected cellulases as well as total CAZome expression were greatly elevated in the Δ2hgt;Δcre-1 strain compared with the WT (Fig. 8a, b). The upregulated genes included those encoding typical “CRE-1 regulon” proteins—an α-amylase A (NCU09805), a starch-binding enzyme (NCU08746), gh6-3 (NCU07190), the invertase inv, a glucoamylase (NCU01517), and a galactosidase (NCU00972) (Additional file 11: Table S5), which were all identified in N. crassa in a previous study  and implied the existence of CRE-1-mediated cross-talk in the Δ2hgt mutant as mentioned above. In contrast, 115 out of 128 upregulated genes in the Δ2hgt strain at the 40-h time point (Table 2) were not typical CRE-1 regulon members . Furthermore, combined deletion of hgt-1/-2 in the Δcre-1 background (strain Δ2hgt; Δcre-1) resulted in upregulation of 317 genes, of which 225 overlapped with the 629 gene set identified in Δ2hgt (Fig. 8c; Table 3). CAZy proteins were also overrepresented in this overlapping gene set (51 out of a total of 70 in the 629 gene set) (Additional file 11: Table S5; Fig. 8d). This observation cannot be completely explained by the action of the CCR regulator CRE-1, and thus further suggests that one or more additional factors might be involved in modulating the HGT-1/-2 suppression signal that leads to the repression of cellulase expression.
Significantly increased intracellular levels of cAMP in Δ2hgt under cellulosic conditions give rise to a negative regulation circuit in the cAMP-PKA pathway
As shown above, one or more unknown factors seem to participate in the connection between HGT-modulated suppression and CAZy expression. The cAMP-PKA signaling pathway, which functions in the sensing of external or internal nutrient changes, is highly conserved across living kingdoms [18, 19]. In N. crassa, misactivation of PKA triggers a number of developmental events, such as apolar growth, aerial hyphal formation, vegetative growth, and improved protein production upon cellulose; in contrast, inactivation leads to premature conidiation, increased thermotolerance, and nearly abolished circadian rhythms [49–51]. Intriguingly, 11 genes involved in glycolysis, two in gluconeogenesis, and two fermentative genes were robustly derepressed in the Δ2hgt mutant, which implies that vegetative growth was facilitated under cellulosic conditions when hgt-1/-2 were deleted (Additional file 11: Table S5). In contrast, a large number of sporulation-related genes were inactivated, including the critical regulator of minor chain formation fl (NCU08726) , the rodlet layer protein eas (NCU08457) , the conidiation-specific gene con-10 (NCU07325) , the all development-altered regulator ada-6 (NCU04866) , the Aspergillus flbC homolog NCU03184 , and the Aspergillus flbD ortholog rca-1 (NCU01312) (Additional file 11: Table S5). The role of rca-1 in conidiation is uncertain, but it was recently revealed to participate in lignocellulolytic enzyme production in N. crassa . All these observations suggest that the cAMP-PKA pathway was activated by the double deletion of hgt-1/-2. To test this hypothesis, intracellular cAMP levels were measured in the WT and the Δ2hgt strain. As anticipated, cAMP levels during the cellulolytic utilization phase were significantly higher in the hgt-deleted mutant relative to the WT (Fig. 9). Additionally, pkac-2 (NCU00682), the paralog of PKAC-1 in N. crassa , was upregulated in the hgt-deleted mutant compared with the WT at all three tested time points from 40 to 72 h (Additional file 11: Table S5), further supporting the involvement of the cAMP-PKA pathway in HGT-modulated suppression.
Replacing a conserved arginine with lysine in HGT-1/-2 leads to dysfunction of glucose transport but no change in CCR signal transduction
Various studies have found that replacing a conserved arginine with a lysine residue in yeast species and the fungal pathogen U. maydis leads to constitutive glucose signaling by hexose sensors, including Snf3p/Rgt2p, Hgt4, Hxs1, and UmHxt1 [11, 13, 58, 59]. To test whether HGT-1/-2 also showed a similar conserved function when this residue was altered, we generated the constructs HGT-1(R172K), HGT-2(R155K), and GLT-1(R167K) (Fig. 10a). Consistent with previous studies , these analogs no longer supported the growth of EBY.VW4000 when glucose was the sole carbon source (Fig. 10b). When the mutant proteins were separately introduced into the Δ2hgt strain, the extracellular protein levels were similar to the background strain on Avicel, seemingly implying complete loss of function via this point mutation (Fig. 10c). With the addition of 2-DG, however, the strains harboring the analogs of the two high-affinity glucose transporters (HGT-1/-2) were much more sensitive to 2-DG than the background strain (Δ2hgt) and the strain with mutated GLT-1. The latter two mutants, Δ2hgt and Δ2hgt::GLT-1(R167K), grew much better than the WT, Δ2hgt::HGT-1(R172K), and Δ2hgt::HGT-2(R155K) (Fig. 10f), indicating that the CCR signal was transmitted by the HGT-1/-2 analogs in the Δ2hgt::HGT-1(R172K) and Δ2hgt::HGT-2(R155K) strains. The strength of CCR in the mutants harboring HGT-1/-2 analogs was reduced compared with the WT (Fig. 10d, e), which may have been because a single hgt sequence could not fully complement the defect of the double deletion strain Δ2hgt (Fig. 5a). Taken together, these data supported the hypothesis that HGT-1/-2 possess an additional function as “glucose transceptors” in N. crassa.
Conservation of the dual-affinity transport components from saprophytes to parasitic fungi
To investigate the conservation of the glucose transport dual-affinity system in fungi, homologs of systems I and II proteins were identified using GLT-1 and HGT-1/-2 protein sequences as queries in searches against the genomes of saprophytes (A. nidulans, A. niger, A. oryzae, Myceliophthora thermophila, and T. reesei), phytopathogens (Botrytis cinerea, Chaetomium globosum, Colletotrichum graminicola, Fusarium graminearum, Magnaporthe oryzae, and U. maydis), and animal fungal pathogens (A. fumigatus and Talaromyces marneffei) (“Methods” section). A phylogenetic analysis was conducted that included hexose transporters from S. cerevisiae (Fig. 11). Except for U. maydis, most tested fungal species possessed a homolog of GLT-1 and one or two HGTs. For example, several functionally identified low- and high-affinity glucose transporters (MstE, MstC/HxtE, and HxtC in A. nidulans [60–62]; MstE, MstA, and MstH in A. niger [63, 64]; and CgHXT5, CgHXT3, and CgHXT1 in C. graminicola ) were clustered into the respective GLT-1, HGT-1, and HGT-2 groups (Fig. 11), suggesting robust conservation of this dual-affinity transport system in the kingdom Fungi. Although S. cerevisiae is widely accepted to possess a glucose dual-affinity transport system (K m = ~20 mM for the low-affinity system and ~1 mM for the high-affinity system), the affinity difference between the two systems in yeast is not so extreme as in filamentous fungi such as N. crassa, reflecting their different evolutionary trajectories shaped by different growth niches in nature.
Dual-affinity systems for nutrient transport in filamentous fungi: revision of an old story
In individual organisms, uptake affinities usually vary widely among different nutrients [8, 10, 14, 66–69]. In filamentous fungi such as A. niger and N. crassa, glucose uptake at the cellular level behaves in a dual-transport fashion [8–10, 70]. Several physiological characteristics of the glucose dual-transporter system in N. crassa are worth noting: (1) the low- and high-affinity systems (systems I and II) cover a 1000-fold difference in affinity making this system a good representative model of glucose uptake; (2) system I is a glucose diffusion system, while system II is an active, H+-co-transport mechanism; and (3) system I is constantly expressed at high glucose levels, whereas system II is subject to repression by glucose and can be de novo synthesized under low- or no-carbon conditions [8–10, 23, 24, 71]. The dual-affinity glucose transport system in N. crassa has been known about for nearly 50 years since the 1970s but has not been investigated systematically at the gene level until the present study, although some characteristics of hgt-1 and GLT-1 have been identified before [33, 40]. Here, three genes were assigned to this dual-system: GLT-1 forms the low-affinity transport system that takes up glucose under adequate glucose conditions, whereas HGT-1 and HGT-2 comprise the high-affinity system that imports sugar when the external glucose concentration is very low (Fig. 12). Whether this dual-affinity transport system, which is highly conserved across fungal species (Fig. 11), increases the adaptability of fungal systems in nature needs to be further investigated. Moreover, because both sensing and import are pivotal for microorganism growth , the finding that HGT-1/-2 perform glucose signaling in addition to their transport functions suggests a novel role for HGT-1/-2 as sensor-like transporters (the so-called “transceptors”) in the dual-affinity transport system, thereby promoting both cell growth and interaction with the environment. Together, the findings here, to our knowledge, represent the first genome-wide molecular characterization of the dual-affinity glucose transport system in filamentous fungi.
Low- vs. high-affinity systems: a trade-off between transport capacity and gene dosage
Despite the apparent exclusiveness of high-affinity and high-capacity, high-affinity transporters, such as hgt-1/-2, are always highly expressed in derepressed conditions (Fig. 2 and Additional file 5: Figure S3). Even though the low-affinity system I gene glt-1 displayed considerably higher expression than other sugar transporters in the presence of high glucose (Fig. 2), its mRNA level was not at the same order of magnitude as those of hgt-1/-2 under low-glucose conditions (e.g., RPKM = 210 for glt-1 vs. 5,700 for hgt-1; Fig. 2b). In contrast, the V max value of GLT-1 was considerably higher than those of HGT-1/-2 (30.75 ± 1.34 mmol h−1 g−1 DCW for GLT-1 vs. 26.42 ± 0.38 and 78.11 ± 4.46 µmol h−1 g−1 DCW for HGT-1 and -2). This suggests that the high-affinity transport system, consisting of low-capacity transporters, compensates for this inherent kinetic handicap by increasing the number of individual proteins to keep pace with the low-affinity system in uptake velocity. The maximum glucose uptake rate of the derepressible system II was demonstrated to be similar to that of the glucose-inducible system I . The highly expressed high-affinity transport system is therefore able to maintain nutrient homeostasis while external nutrients are being depleted, thereby providing sufficient preparation time for starvation and consequent cell recovery when nutrients are replenished .
Glucose transport by the derepressible glucose transporters HGT-1/-2 suppresses lignocellulase gene expression and vegetative growth under carbon-limited environments, thereby favoring fungal sporulation for survival
In a number of filamentous fungi, deletion of cre-1/creA/cre1 promotes the utilization of alternative carbon sources, particularly through derepression of cellulolytic enzymes [4, 73–75]. Although deletion of cre-1 in the Δ2hgt strain resulted in faster growth on Avicel (Fig. 7), simply concluding that the hyper-production of cellulases in Δ2hgt is related to the regulator CRE-1 is not warranted: single deletion of cre-1 only resulted in approximately 30% higher protein production than in the WT , far less than the 180–210% increases are seen in the Δ2hgt and Δ2hgt;Δcre-1 mutants (Fig. 7a). HGT-1/-2 appear to connect extracellular glucose signaling with the internal CCR signal, because the Δ2hgt mutant displayed moderate resistance to 2-DG inhibition (Fig. 6d, e). This HGT-mediated suppression is partly associated with CRE-1-mediated CCR, as downregulation of cre-1 in Δ2hgt occurred in the late phase of cellulose utilization (Fig. 6a; Additional file 10: Table S4). Despite these findings, exactly how HGT-1/-2 function with CRE-1 (or other unknown factors) to mediate CCR in N. crassa remains to be elucidated. Intriguingly, the transcriptional regulator vib-1, which was recently found to coordinate glucose signaling and CCR during plant cell wall degradation , was not strongly upregulated when hgt-1/-2 were deleted (Additional file 10: Table S4). Future investigations of the regulatory network consisting of VIB-1-mediated CCR and the HGT-mediated glucose response should help to decipher fungal glucose signaling and CCR regulation.
In S. cerevisiae, extracellular glucose can be sensed by two membrane transporter-like sensors, Rgt2p and Snf3p . These two sensors possess a long cytoplasmic C-terminus harboring one or two short conserved sequence blocks that are essential for glucose signal transduction . The identical homolog of Rgt2p/Snf3p in N. crassa is RCO-3 (Fig. 11), which also contains an extremely long C-terminus in the cytoplasm (Fig. 10a). Mutation of rco-3 in N. crassa leads to severe cell defects in glucose uptake as well as decreased growth in the presence of high glucose levels, but the ability to sense low levels of glucose is probably still retained . These results were corroborated in our study: qRT-PCR analysis revealed that glt-1 was dramatically downregulated (~100-fold) in Δrco-3 compared with the WT, resulting in synergistically elevated expression (ten- to 30-fold) of hgt-1/-2 and NCU06358 (Additional file 14: Figure S9). Because the derepression of high-affinity glucose transporters hgt-1/-2 was not affected in the Δrco-3 mutant, RCO-3 might act as a low-affinity glucose sensor , whereas HGT-1/-2 are high-affinity ones.
Intriguingly, glycolysis and fermentation were stimulated in Δ2hgt (Additional file 11: Table S5). This phenomenon is quite similar to the “Warburg effect” in cancer cells, where glycolysis is improperly activated by several important kinases, such as PI3K/AKT and tyrosine kinase . Similarly, the cAMP-PKA signaling pathway cannot be excluded as a possible explanation for these HGT-dependent functions (Fig. 12): (1) the internal cAMP level was significantly higher in the Δ2hgt mutant during cellulose utilization (Fig. 9); (2) pkac-2, one of only two PKA catalytic subunits reported to date in N. crassa , was strongly and constantly upregulated (Additional file 11: Table S5); and (3) asexual conidiation-related genes were extensively inactivated in Δ2hgt (Additional file 11: Table S5). In N. crassa, inactivation of PKA leads to strong derepression of submerged-culture conidiation as well as aerial sporulation [57, 77]. Future investigations of the connection between HGT and the cAMP-PKA signaling pathway as well as the specific roles of HGT in glycolysis, plant cell wall deconstruction, and sporulation will help us to understand the architecture of nutrient signaling regulation in filamentous fungi (Fig. 12).
Sporulation and quiescence are conserved strategies for survival when organisms encounter nutrient depletion or other chemical/physical stresses [78, 79]. In this regard, signaling to suppress vegetative growth and activate sporulation is a preferred mechanism for fungi to overcome severe growth environments. As revealed by this study, HGTs are strongly derepressed under carbon-limited or lignocellulolytic conditions, which lead to the inactivation of glycolysis but activation of asexual sporulation (Fig. 12). HGT-mediated repression of cellulosic growth therefore favors fungal conidiation, which might be beneficial for species survival under carbon-limited environments. Conversely, upon encountering plentiful glucose, it should be noted that N. crassa can inhibit system II  and synergistically elevate the expression of the low-affinity transporter glt-1 for cell proliferation (Fig. 12).
In this study, HGT-1/-2 and GLT-1 were identified as major derepressible and glucose-inducible components, respectively, of the high- and low-affinity glucose transport systems in N. crassa. Growth defects due to loss of glt-1 were restored by upregulation of the high-affinity transporters hgt-1/-2. In addition to their glucose-transporting functions under carbon-limited or cellulolytic conditions, HGT-1/-2 also mediate glucose signaling to connect extracellular nutrient availability with internal catabolite repression and metabolism, and thus may act as glucose transceptors. Simultaneous deletion of hgt-1 and hgt-2 leads to comprehensive derepression of a large group of genes including those encoding glycolysis enzymes and plant cell wall-degrading enzymes, which are associated with CRE-1-related CCR, transcriptional regulation by CLR-1/-2, and the cAMP-PKA signaling pathway. Given the wide conservation of these dual-affinity transport components across the fungal kingdom, future investigations of GLT-1 and HGT-1/-2, and studies of HGT-mediated signaling in other fungal systems will shed new light on long-standing questions about the physiological roles and evolutionary traits of the dual-affinity transport system, and thus inform studies on fungal glucose signaling, pathogenicity, and plant cell wall deconstruction for biorefinery.
Strains, media, and culture conditions
The N. crassa WT strain (FGSC2489) and mutants Δglt-1 (FGSC13161), Δhgt-1 (FGSC22819), Δhgt-2 (FGSC18807), ΔNCU05897 (FGSC13717), ΔNCU00450 (FGSC15906), and Δrco-3 (FGSC17928) were obtained from the Fungal Genetics Stock Center (FGSC) . The mutant Δcre-1 was a gift from the laboratory of Professor N. Louise Glass of the University of California, Berkeley. Multiple deletion mutants were generated via sexual crosses of the above-mentioned mutants (FGSC protocol, http://www.fgsc.net/Neurospora/NeurosporaProtocolGuide.htm). Artificial misexpression strains were constructed by transforming the Δ2hgt;his-3 − strain with the linearized plasmid pMF272 harboring the clock-controlled gene 1 (ccg-1) or M. grisea ribosomal protein 27 (MgRP27) [44, 48] promoters, various gene-coding sequences or point-mutated analogs, and flanking regions from the his-3 gene sequence . Transformants with histidine prototrophy were selected for further purification via microconidial separation (FGSC protocol, http://www.fgsc.net/Neurospora/NeurosporaProtocolGuide.htm). At least three purified biological transformants per selected gene were used for downstream experiments. N. crassa strains were pre-grown on slants containing 3 mL Vogel’s minimal medium  with 2.0% (w/v) sucrose as a sole carbon source for 1–2 days in darkness at 28 °C, followed by constant light for 6–8 days at room temperature to obtain mature conidia. Unless indicated, N. crassa conidia were inoculated into liquid Vogel’s minimal medium with various carbon sources at 106 conidial mL−1 with constant light shaking at 200 rpm, 25 °C. For CCR sensitivity analysis, 2-DG (Sigma-Aldrich, St. Louis, MO, USA) was added to 100 mL liquid Vogel’s minimal medium containing 2.0% Avicel PH-101 (Sigma-Aldrich) to a final concentration of 0.1%.
For media shift experiments, N. crassa cultures were pre-grown on 100 mL Vogel’s salts supplemented with 2.0% sucrose for 16 h, with the generated mycelia then filtered through six layers of gauze and immediately washed with sterilized water at least five times. The mycelia from each flask pre-culture were then transferred to 100 mL Vogel’s medium containing 0–10% glucose as the sole carbon source, followed by 60-min growth according to previous studies . The mycelia were subsequently sampled via filtration and flash frozen in liquid nitrogen for total RNA extraction. For low-glucose uptake assays, the washed mycelia were transferred to 100 mL carbon-free Vogel’s medium with shaking for 90 min . The carbon-starved mycelia were then centrifuged at 3220×g for 5 min, and approximately 4 mL wet mycelia were transferred into new medium containing 100 mL Vogel’s salts plus 0.02% glucose with shaking at 25 °C and 200 rpm. Culture supernatants (500 µL) were taken at indicated time points (0, 5, 10, 20, and 30 min) and immediately filtered through 0.22-µm membranes. Each mycelial culture used for the glucose uptake assay was harvested, dried, and weighed. Glucose concentrations were measured with a Megazyme d-Glucose Assay kit according to the manufacturer’s instructions (Megazyme, Wicklow, Ireland). Equivalent glucose uptake was defined as the amount of glucose consumed per gram of dried mycelia.
Saccharomyces cerevisiae EBY.VW4000 , a gift from Professor Eckhard Boles, was grown in YPM medium (1.0% yeast extract, 2.0% peptone, 2.0% maltose, and optionally 2.0% agar) for subsequent use. Recombinant strains were cultured in synthetic complete medium (SC) supplemented with drop-out amino acids lacking uracil (Ura−) and 2.0% glucose or maltose as the sole carbon source (with or without 2.0% agar). For the growth complementation assay, recombinant strains were incubated in liquid SC (Ura−) medium with maltose. Cells were harvested at an optical density at 600 nm (OD600) of 1.0–2.0, washed twice with distilled water, and resuspended in distilled water to an OD600 of 0.25–0.30. Serially diluted cells (×1, ×10, and ×100; 2 µL for each diluted clone) were plated on solid SC (Ura−) medium containing either maltose or glucose. Growth on maltose or glucose medium at 3 or 6 days, respectively, was photographed with a Canon camera (Canon EOS207, Japan).
The open reading frames (ORFs) of cre-1, hgt-1/-2, and glt-1 were amplified from the cDNAs of N. crassa. For misexpression of these genes in the Δ2hgt mutant, their ORFs were inserted into the multiple cloning sites of the shuttle vector pMF272  driven by either the ccg-1 or the MgRP27 promoter and tagged with enhanced green fluorescence protein (eGFP). The analogs hgt-1(R172K), hgt-2(R155K), and glt-1(R167K) were generated by site-directed mutagenesis using high-fidelity PCR polymerase (Thermo Fisher Scientific, Waltham, MA, USA). For heterologous expression in S. cerevisiae, these transporter ORFs were recombined with the yeast shuttle vector pRS426  fused with the native phosphoglycerate kinase-1 (pgk-1) promoter and the eGFP reporter gene. Cloning primers are listed in Additional file 15: Table S6. Restriction enzymes were purchased from Thermo Fisher Scientific. All recombinant plasmids were amplified in Escherichia coli strain DH5α and sequenced for gene authenticity.
RNA extraction, sequencing, and data analysis
Mycelial sampling and total RNA extraction were performed as described in a previous study . Briefly, cultured mycelia were harvested by filtration and immediately frozen in liquid nitrogen. Total RNA was isolated with Trizol reagent (Invitrogen, Carlsbad, CA, USA) and further purified using DNase I (RNeasy Mini kit, Qiagen, Hilden, Germany). The RNA concentration and OD260/OD280 were measured with a Nanodrop 2000c (Thermo Scientific), and RNA integrity was checked by agarose gel electrophoresis and using an Agilent 2100 (Agilent Technologies, Santa Clara, CA, USA).
Qualified RNA with an OD260/OD280 >1.8 and RIN (RAN Integrity Number) >7.0 was prepared according to Shenzhen BGI (Shenzhen, China) and Novogene (Tianjin, China) standard protocols, and sequenced on an Illumina HiSeq 2000/2500 platform (San Diego, CA, USA). The biological reproducibility of the RNA-Seq data was demonstrated to be high through sample-to-sample clustering  and Spearman correlation analysis (Additional file 1: Figure S1, Additional file 9: Figure S6). Read mapping and counting using TopHat2 (version 2.0.12)  and HTSeq (version 0.6.0)  was performed as previously described . The reads mapped to each transcript were used to calculate normalized transcript abundance (as RPKM ) and to perform differential gene expression analysis in GFOLD (version 1.1.0)  and DESeq2 (version 1.2.10) . Genes with a DESeq2 P value <1 × 10−4 were considered to be statistically differentially expressed at a robust significance level. Because genes with |GFOLD| > 1 are empirically more likely to be of biological importance , a combined criterion of |GFOLD| > 1 and P < 1 × 10−4 was applied for genome-wide differential gene analysis. DESeq2  was used to estimate expression fold-changes, which represent apparent expression changes that may also be useful for assessment of biological expression variations. RNA-Seq raw data are available at the Gene Expression Omnibus under accession number GSE78952.
Hierarchical clustering analysis was performed using Cluster 3.0  or the pheatmap package in R (version 1.0.2) (http://www.r-projectorg). To generate a clustering heatmap, the RPKM values of each gene were log-transformed and calculated by the complete linkage method with Euclidean distance as the similarity metric. Selected differentially expressed genes were submitted to the MIPS Functional Category Database , and significantly enriched pathways were estimated.
Quantitative real-time reverse-transcription polymerase chain reaction was performed using SYBR Green Realtime PCR master mix (Toyobo, Osaka, Japan) according to the manufacturer’s instructions on a CFX96 real-time PCR detection system (Bio-Rad, Hercules, USA). Each reaction was conducted in duplicate or triplicate. The actin gene (NCU04173) was used as an endogenous control for all experiments. All primers used in this study are listed in Additional file 15: Table S6. The relative expression level of each gene was calculated using the 2−ΔΔCt method .
Samples submerged in Vogel’s medium with 2.0% Avicel and cultured for the indicated time periods were harvested by filtration and immediately frozen in liquid nitrogen. Mycelia were finely ground in liquid nitrogen, transferred to 1.0% hydrochloric acid, briefly vortexed, and frozen at −80 °C until use. Before cAMP measurements, the samples were thawed at 4 °C and centrifuged at 12,000×g for 15 min at 4 °C. The supernatant was used for the cAMP assay following the Applied Biosystems (Waltham, MA, USA) protocol. The protein concentration of each supernatant was quantified by the Bradford method (Bio-Rad). The protein in the pellet was solubilized in 0.5% sodium dodecyl sulfate plus 0.1 M sodium hydroxide, vortexed for 30 s, and also measured using the Bradford method.
Secreted enzyme and dried mycelial weight assays
For secreted protein assays, 800 µL of each culture supernatant was collected during the growth period (2–7 days with Avicel and 1–6 days with cellobiose), centrifuged at 15,294×g for 8 min to remove mycelia, and stored at 4 °C for analysis within 24 h or at −20 °C for sodium dodecyl sulfate polyacrylamide gel electrophoresis (Novex NuPAGE Pre-cast Protein Gels; Thermo Fisher Scientific).
The total secreted protein content was determined using the Bradford method (Bio-Rad) with bovine serum albumin as a standard. Carboxymethyl cellulase and xylanase activities were measured using an Azo-CMC/xylan kit (Megazyme) according to the manufacturer’s instructions. Exoglucanase activity was assayed as previously described using p-nitrophenyl-d-cellobioside (Sigma-Aldrich) as the substrate .
Mycelia grown on sucrose, glucose, and cellobiose for designated times were harvested, dried, and weighed. Biomass dry weights of the Avicel cultures were measured according to a previous study  with a slight modification. In brief, 5 mL of thoroughly mixed culture broth was centrifuged at 3220×g for 5 min. After discarding the supernatant, 3 mL 80% (v/v) acetic acid:concentrated nitric acid (10:1, v/v) reagent was added and the mixture was boiled in water for 1 h to solubilize the fungal biomass. This procedure was repeated with a fresh 3-mL aliquot of acetic acid:nitric acid reagent. The reaction mixture (residual Avicel) was then centrifuged, dried, and weighed. Mycelial dry weight was defined as the dry weight of the original 5-mL culture minus that of the reaction mixture.
Microscopy and imaging
Visualization of eGFP-tagged glucose transporters in S. cerevisiae and N. crassa was conducted using a 100 × 1.4 NA oil immersion objective on a Leica TCS SP5 II laser scanning confocal microscope (Leica, Wetzlar, Germany). For confocal microscopy of yeast, single clones of recombinant EBY.VW4000 strains harboring various glucose transporters were inoculated into SC(Ura−) liquid medium supplemented with 2.0% maltose and grown overnight to an OD600 of 1.0–2.0. The cell cultures were centrifuged, washed twice with sterile water, and resuspended in phosphate-buffered saline solution. For microscopy of N. crassa, recombinant N. crassa was pre-grown on Vogel’s medium with 2.0% sucrose for 16 h and then briefly washed with sterilized water at least five times and resuspended in Vogel’s medium with either 2.0% Avicel, glucose, or no added carbon for 4 h with shaking at 200 rpm and 25 °C. Before confocal scanning, mycelia were treated with 1 µg mL−1 4′,6-diamidino-2-phenylindole for 15 min as needed. Images were processed using the Leica Microsystems LAS AF-TCS MP (version 2.4.1) and ImageJ (version 1.47) software.
Radiolabeled glucose transport
A glucose transport assay in yeast was performed according to a previously published protocol .
Amino acid sequences of GLT-1 and HGT-1/-2 (http://www.broadinstitute.org) were used as queries in BLASTp searches against the protein sequences in the genomes of A. nidulans, A. niger, A. oryzae, A. fumigatus, M. thermophile, T. reesei, B. cinerea, C. globosum, C. graminicola, F. graminearum, M. oryzae, U. maydis, and T. marneffei in the National Center for Biotechnology Information (NCBI) and DOE Joint Genome Institute (JGI) databases. The best hits with identity >50%, E-value <1 × 10−10, and coverage >60% were selected as the closest homologs of GLT-1 and HGT-1/-2, with the exception of two GLT-1 homologs (Um XP_758184.1: 43% identity; Mo XP_003719105.1: 40% identity). Among these highly conserved homologs, under-documented low- and high-affinity glucose transporters served as positive controls. Hexose transporters of S. cerevisiae and H. sapiens, and reported pure sensors in S. cerevisiae and N. crassa were included as outgroups. The protein sequences of these MFS transmembrane proteins were aligned with ClustalW  using the following parameters: protein weight matrix, Gonnet; gap opening penalty, 10; gap extension, 0.2; and gap distance, 5. A phylogenetic tree was constructed under the JTT amino acid substitution model by maximum likelihood in MEGA 5.2 , with 1000 bootstrap replicates.
All figures were plotted on the R program platform (http://www.r-project.org/) or Excel 2013. Unless otherwise indicated, values represent the means of at least three replicates; error bars show standard deviations.
cyclic adenosine monophosphate
carbon catabolite repression
enhanced green fluorescence protein
Fungal Genetics Stock Center
gene expression omnibus
major facilitator superfamily
protein kinase A
quantitative real-time reverse transcription polymerase chain reaction
reads per kilobase per million mapped reads
sodium dodecyl sulfate polyacrylamide gel electrophoresis
Jones AM, Xuan Y, Xu M, Wang RS, Ho CH, Lalonde S, You CH, Sardi MI, Parsa SA, Smith-Valle E, et al. Border control–a membrane-linked interactome of Arabidopsis. Science. 2014;344:711–6.
Chantranupong L, Wolfson RL, Sabatini DM. Nutrient-sensing mechanisms across evolution. Cell. 2015;161:67–83.
Horak J. Regulations of sugar transporters: insights from yeast. Curr Genet. 2013;59:1–31.
Sun J, Glass NL. Identification of the CRE-1 cellulolytic regulon in Neurospora crassa. PLoS ONE. 2011;6:e25654.
Xiong Y, Sun J, Glass NL. VIB1, a link between glucose signaling and carbon catabolite repression, is essential for plant cell wall degradation by Neurospora crassa. PLoS Genet. 2014;10:e1004500.
Reddy VS, Shlykov MA, Castillo R, Sun EI, Saier MH. The major facilitator superfamily (MFS) revisited. FEBS J. 2012;279:2022–35.
Leandro MJ, Fonseca C, Gonçalves P. Hexose and pentose transport in ascomycetous yeasts: an overview. FEMS Yeast Res. 2009;9:511–25.
Scarborough GA. Sugar transport in Neurospora crassa: II. A second glucose transport system. J Biol Chem. 1970;245:3985–7.
Neville MM, Subkind SR, Roseman S. A derepressible active transport system for glucose in Neurospora crassa. J Biol Chem. 1971;246:1294–301.
Schneider RP, Wiley WR. Kinetic characteristics of the two glucose transport systems in Neurospora crassa. J Bacteriol. 1971;106:479–86.
Ozcan S, Dover J, Rosenwald AG, Wölfl S, Johnston M. Two glucose transporters in Saccharomyces cerevisiae are glucose sensors that generate a signal for induction of gene expression. Proc Natl Acad Sci USA. 1996;93:12428–32.
Ozcan S, Dover J, Johnston M. Glucose sensing and signaling by two glucose receptors in the yeast Saccharomyces cerevisiae. EMBO J. 1998;17:2566–73.
Schuler D, Wahl R, Wippel K, Vranes M, Münsterkötter M, Sauer N, Kämper J. Hxt1, a monosaccharide transporter and sensor required for virulence of the maize pathogen Ustilago maydis. New Phytol. 2015;206:1086–100.
Ho C-H, Lin S-H, Hu H-C, Tsay Y-F. CHL1 functions as a nitrate sensor in plants. Cell. 2009;138:1184–94.
Van Zeebroeck G, Bonini BM, Versele M, Thevelein JM. Transport and signaling via the amino acid binding site of the yeast Gap1 amino acid transceptor. Nat Chem Biol. 2009;5:45–52.
Popova Y, Thayumanavan P, Lonati E, Agrochao M, Thevelein JM. Transport and signaling through the phosphate-binding site of the yeast Pho84 phosphate transceptor. Proc Natl Acad Sci USA. 2010;107:2890–5.
Znameroski EA, Li X, Tsai JC, Galazka JM, Glass NL, Cate JHD. Evidence for transceptor function of cellodextrin transporters in Neurospora crassa. J Biol Chem. 2014;289:2610–9.
Conrad M, Schothorst J, Kankipati HN, Van Zeebroeck G, Rubio-Texeira M, Thevelein JM. Nutrient sensing and signaling in the yeast Saccharomyces cerevisiae. FEMS Microbiol Rev. 2014;38:254–99.
Broach JR. Nutritional control of growth and development in yeast. Genetics. 2012;192:73–105.
Park H-S, Yu J-H. Genetic control of asexual sporulation in filamentous fungi. Curr Opin Microbiol. 2012;15:669–77.
Li L, Borkovich KA. GPR-4 is a predicted G-protein-coupled receptor required for carbon source-dependent asexual growth and development in Neurospora crassa. Eukaryot Cell. 2006;5:1287–300.
Madi L, McBride SA, Bailey LA, Ebbole DJ. rco-3, a gene involved in glucose transport and conidiation in Neurospora crassa. Genetics. 1997;146:499–508.
Schneider RP, Wiley WR. Regulation of sugar transport in Neurospora crassa. J Bacteriol. 1971;106:487–92.
Scarborough GA. Sugar transport in Neurospora crassa. J Biol Chem. 1970;245:1694–8.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Coradetti ST, Craig JP, Xiong Y, Shock T, Tian C, Glass NL. Conserved and essential transcription factors for cellulase gene expression in ascomycete fungi. Proc Natl Acad Sci USA. 2012;109:7397–402.
Rep M, Albertyn J, Thevelein JM, Prior BA, Hohmann S. Different signalling pathways contribute to the control of GPD1 gene expression by osmotic stress in Saccharomyces cerevisiae. Microbiology. 1999;145:715–27.
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–45.
Ren Q, Chen K, Paulsen IT. TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res. 2007;35:D274–9.
Brown NA, Ries LNA, Goldman GH. How nutritional status signalling coordinates metabolism and lignocellulolytic enzyme secretion. Fungal Genet Biol. 2014;72:48–63.
Sharon A, Finkelstein A, Shlezinger N, Hatam I. Fungal apoptosis: function, genes and gene function. FEMS Microbiol Rev. 2009;33:833–54.
Katz ME, Buckland R, Hunter CC, Todd RB. Distinct roles for the p53-like transcription factor XprG and autophagy genes in the response to starvation. Fungal Genet Biol. 2015;83:10–8.
Li JG, Lin LC, Li HY, Tian CG, Ma YH. Transcriptional comparison of the filamentous fungus Neurospora crassa growing on three major monosaccharides d-glucose, d-xylose and l-arabinose. Biotechnol Biofuels. 2014;7:31.
Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma L-J, Smirnov S, Purcell S, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422:859–68.
Galazka JM, Tian CG, Beeson WT, Martinez B, Glass NL, Cate JHD. Cellodextrin transport in yeast for improved biofuel production. Science. 2010;330:84–6.
Li JG, Xu J, Cai PL, Wang B, Ma YH, Benz JP, Tian CG. Functional analysis of two l-arabinose transporters from filamentous fungi reveals promising characteristics for improved pentose utilization in Saccharomyces cerevisiae. Appl Environ Microbiol. 2015;81:4062–70.
Benz JP, Protzko R, Andrich J, Bauer S, Dueber J, Somerville C. Identification and characterization of a galacturonic acid transporter from Neurospora crassa and its application for Saccharomyces cerevisiae fermentation processes. Biotechnol Biofuels. 2014;7:20.
Xiong Y, Coradetti ST, Li X, Gritsenko MA, Clauss T, Petyuk V, Camp D, Smith R, Cate JHD, Yang F, Glass NL. The proteome and phosphoproteome of Neurospora crassa in response to cellulose, sucrose and carbon starvation. Fungal Genet Biol. 2014;71:21–33.
Du J, Li S, Zhao H. Discovery and characterization of novel d-xylose-specific transporters from Neurospora crassa and Pichia stipitis. Mol BioSyst. 2010;6:2150–6.
Xie X, Wilkinson HH, Correa A, Lewis ZA, Bell-Pedersen D, Ebbole DJ. Transcriptional response to glucose starvation and functional analysis of a glucose transporter of Neurospora crassa. Fungal Genet Biol. 2004;41:1104–19.
Wieczorke R, Krampe S, Weierstall T, Freidel K, Hollenberg CP, Boles E. Concurrent knock-out of at least 20 transporter genes is required to block uptake of hexoses in Saccharomyces cerevisiae. FEBS Lett. 1999;464:123–8.
Cai PL, Gu RM, Wang B, Li JG, Wan L, Tian CG, Ma YH. Evidence of a critical role for cellodextrin transporter 2 (CDT-2) in both cellulose and hemicellulose degradation and utilization in Neurospora crassa. PLoS ONE. 2014;9:e89330.
Wang B, Cai PL, Sun WL, Li JG, Tian CG, Ma YH. A transcriptomic analysis of Neurospora crassa using five major crop residues and the novel role of the sporulation regulator rca-1 in lignocellulase production. Biotechnol Biofuels. 2015;8:21.
McNally M, Free S. Isolation and characterization of a Neurospora glucose-repressible gene. Curr Genet. 1988;14:545–51.
Allen KE, McNally MT, Lowendorf HS, Slayman CW, Free SJ. Deoxyglucose-resistant mutants of Neurospora crassa: isolation, mapping, and biochemical characterization. J Bacteriol. 1989;171:53–8.
Ziv C, Gorovits R, Yarden O. Carbon source affects PKA-dependent polarity of Neurospora crassa in a CRE-1-dependent and independent manner. Fungal Genet Biol. 2008;45:103–16.
Trumbly RJ. Glucose repression in the yeast Saccharomyces cerevisiae. Mol Microbiol. 1992;6:15–21.
Bourett TM, Sweigard JA, Czymmek KJ, Carroll A, Howard RJ. Reef coral fluorescent proteins for visualizing fungal pathogens. Fungal Genet Biol. 2002;37:211–20.
Terenzi HF, Flawia MM, Torres HN. A Neurospora crassa morphological mutant showing reduced adenylate cyclase activity. Biochem Biophys Res Commun. 1974;58:990–6.
Bruno KS, Aramayo R, Minke PF, Metzenberg RL, Plamann M. Loss of growth polarity and mislocalization of septa in a Neurospora mutant altered in the regulatory subunit of cAMP-dependent protein kinase. EMBO J. 1996;15:5772–82.
Huang G, Chen S, Li S, Cha J, Long C, Li L, He Q, Liu Y. Protein kinase A and casein kinases mediate sequential phosphorylation events in the circadian negative feedback loop. Genes Dev. 2007;21:3283–95.
Bailey-Shrode L, Ebbole DJ. The fluffy gene of Neurospora crassa is necessary and sufficient to induce conidiophore development. Genetics. 2004;166:1741–9.
Beever RE, Dempsey GP. Function of rodlets on the surface of fungal spores. Nature. 1978;272:608–10.
Berlin V, Yanofsky C. Isolation and characterization of genes differentially expressed during conidiation of Neurospora crassa. Mol Cell Biol. 1985;5:849–55.
Colot HV. A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors. Proc Natl Acad Sci USA. 2006;103:10352–7.
Kwon NJ, Garzia A, Espeso EA, Ugalde U, Yu JH. FlbC is a putative nuclear C2H2 transcription factor regulating development in Aspergillus nidulans. Mol Microbiol. 2010;77:1203–19.
Banno S, Ochiai N, Noguchi R, Kimura M, Yamaguchi I, Kanzaki S-I, Murayama T, Fujimura M. A catalytic subunit of cyclic AMP-dependent protein kinase, PKAC-1, regulates asexual differentiation in Neurospora crassa. Genes Genet Syst. 2005;80:25–34.
Brown V, Sexton JA, Johnston M. A glucose sensor in Candida albicans. Eukaryot Cell. 2006;5:1726–37.
Stasyk OG, Maidan MM, Stasyk OV, Van Dijck P, Thevelein JM, Sibirny AA. Identification of hexose transporter-like sensor HXS1 and functional hexose transporter HXT1 in the methylotrophic yeast Hansenula polymorpha. Eukaryot Cell. 2008;7:735–46.
Forment JV, Flipphi M, Ramón D, Ventura L, MacCabe AP. Identification of the mstE gene encoding a glucose-inducible, low affinity glucose transporter in Aspergillus nidulans. J Biol Chem. 2006;281:8339–46.
Forment JV, Flipphi M, Ventura L, González R, Ramón D, MacCabe AP. High-affinity glucose transport in Aspergillus nidulans is mediated by the products of two related but differentially expressed genes. PLoS ONE. 2014;9:e94662.
dos Reis TF, Menino JF, Bom VLP, Brown NA, Colabardini AC, Savoldi M, Goldman MHS, Rodrigues F, Goldman GH. Identification of glucose transporters in Aspergillus nidulans. PLoS ONE. 2013;8:e81412.
Vankuyk PA, Diderich JA, MacCabe AP, Hererro O, Ruijter GJG, Visser J. Aspergillus niger mstA encodes a high-affinity sugar/H+ symporter which is regulated in response to extracellular pH. Biochem J. 2004;379:375–83.
Sloothaak J, Odoni DI, de Graaff LH, Dos Santos VAM, Schaap PJ, Tamayo-Ramos JA. Aspergillus niger membrane-associated proteome analysis for the identification of glucose transporters. Biotechnol Biofuels. 2015;8:150.
Lingner U, Münch S, Deising HB, Sauer N. Hexose transporters of a hemibiotrophic plant pathogen: FUNCTIONAL VARIATIONS AND REGULATORY DIFFERENCES AT DIFFERENT STAGES OF INFECTION. J Biol Chem. 2011;286:20913–22.
Gjedde A. High- and low-affinity transport of d-glucose from blood to brain. J Neurochem. 1981;36:1463–71.
Rosenberg LE, Albrecht I, Segal S. Lysine transport in human kidney: evidence for two systems. Science. 1967;155:1426–8.
Levy S, Kafri M, Carmi M, Barkai N. The competitive advantage of a dual-transporter system. Science. 2011;334:1408–12.
Rodríguez-Navarro A, Ramos J. Dual system for potassium transport in Saccharomyces cerevisiae. J Bacteriol. 1984;159:940–5.
Torres NV, Riol-Cimas JM, Wolschek M, Kubicek CP. Glucose transport by Aspergillus niger: the low-affinity carrier is only formed during growth on high glucose concentrations. Appl Microbiol Biotechnol. 1996;44:790–4.
Slayman CL, Slayman CW. Depolarization of the plasma membrane of Neurospora during active transport of glucose: evidence for a proton-dependent cotransport system. Proc Natl Acad Sci USA. 1974;71:1935–9.
Youk H, van Oudenaarden A. Growth landscape formed by perception and import of glucose in yeast. Nature. 2009;462:875–9.
Li Z, Yao G, Wu R, Gao L, Kan Q, Liu M, Yang P, Liu G, Qin Y, Song X, et al. Synergistic and dose-controlled regulation of cellulase gene expression in Penicillium oxalicum. PLoS Genet. 2015;11:e1005509.
Nakari-Setala T, Paloheimo M, Kallio J, Vehmaanpera J, Penttila M, Saloheimo M. Genetic modification of carbon catabolite repression in Trichoderma reesei for improved protein production. Appl Environ Microbiol. 2009;75:4853–60.
Orejas M, MacCabe AP, Perez Gonzalez JA, Kumar S, Ramon D. Carbon catabolite repression of the Aspergillus nidulans xlnA gene. Mol Microbiol. 1999;31:177–84.
Vander Heiden MG, Cantley LC, Thompson CB. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009;324:1029–33.
Kays AM, Rowley PS, Baasiri RA, Borkovich KA. Regulation of conidiation and adenylyl cyclase levels by the Galpha protein GNA-3 in Neurospora crassa. Mol Cell Biol. 2000;20:7693–705.
Gray JV, Petsko GA, Johnston GC, Ringe D, Singer RA, Werner-Washburne M. “Sleeping beauty”: quiescence in Saccharomyces cerevisiae. Microbiol Mol Biol Rev. 2004;68:187–206.
Adams TH, Wieser JK, Yu JH. Asexual sporulation in Aspergillus nidulans. Microbiol Mol Biol Rev. 1998;62:35–54.
McCluskey K. The fungal genetics stock center: from molds to molecules. Adv Appl Microbiol. 2003;52:245–62 (Academic Press).
Margolin B, Freitag M, Selker E. Improved plasmids for gene targeting at the his-3 locus of Neurospora crassa by electroporation. Fungal Genet Newslett. 1997;44:34–6.
Vogel HJ. A convenient growth medium for Neurospora. Microb Genet Bull. 1956;13:42–6.
Freitag M, Hickey PC, Raju NB, Selker EU, Read ND. GFP as a tool to analyze the organization, dynamics and function of nuclei and microtubules in Neurospora crassa. Fungal Genet Biol. 2004;41:897–910.
Christianson TW, Sikorski RS, Dante M, Shero JH, Hieter P. Multifunctional yeast high-copy-number shuttle vectors. Gene. 1992;110:119–22.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.
Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2014;31:166–9.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Feng J, Meyer CA, Wang Q, Liu JS, Liu XS, Zhang Y. GFOLD: a generalized fold change for ranking differentially expressed genes from RNA-seq data. Bioinformatics. 2012;28:2782–8.
de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20:1453–4.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.
Thompson J, Higgins D, Gibson T. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
BW and CT conceived and designed the project. BW, JL, and JG performed experiments. BW, JG, and CT analyzed the data. BW and CT wrote the paper. JL, JG, PC, and XH participated in critical discussion and draft revision. All authors read and approved the final manuscript.
We thank Professor N. Louise Glass and Professor J. Philipp Benz for manuscript critical reading and Professor Eckhard Boles for generously supplying the hexose-transporter-null yeast strain EBY.VW4000. We also thank Ms. Lixian Wang, Mr. Yong Chen, and Mr. Jinxiao Ji for technical assistance. This work was supported by the National Natural Science Foundation of China (31670042, 31471186) and the key project of Chinese Academy of Sciences (ZDRW-ZS-2016-3).
The Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences and the authors have made a patent application for the glucose transporters GLT-1 and HGT-1/-2 for potential applications in biotechnology. There are no non-financial competing interests for any of the authors.
Availability of data and materials
All supporting data are available, including raw RNA-Seq data, which have been deposited into the Gene Expression Omnibus under accession number GSE78952.