Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Unravelling the molecular basis for light modulated cellulase gene expression - the role of photoreceptors in Neurospora crassa

Monika Schmoll12*, Chaoguang Tian13, Jianping Sun1, Doris Tisch2 and N Louise Glass1

Author affiliations

1 Plant and Microbial Biology Department, University of California, Berkeley, CA 94720, USA

2 Research Area Gene Technology and Applied Biochemistry, Vienna University of Technology, Getreidemarkt 9, 1060 Wien, Austria

3 Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin 300308, China

For all author emails, please log on.

Citation and License

BMC Genomics 2012, 13:127  doi:10.1186/1471-2164-13-127

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/13/127


Received:21 November 2011
Accepted:31 March 2012
Published:31 March 2012

© 2012 Schmoll et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Light represents an important environmental cue, which exerts considerable influence on the metabolism of fungi. Studies with the biotechnological fungal workhorse Trichoderma reesei (Hypocrea jecorina) have revealed an interconnection between transcriptional regulation of cellulolytic enzymes and the light response. Neurospora crassa has been used as a model organism to study light and circadian rhythm biology. We therefore investigated whether light also regulates transcriptional regulation of cellulolytic enzymes in N. crassa.

Results

We show that the N. crassa photoreceptor genes wc-1, wc-2 and vvd are involved in regulation of cellulase gene expression, indicating that this phenomenon is conserved among filamentous fungi. The negative effect of VVD on production of cellulolytic enzymes is thereby accomplished by its role in photoadaptation and hence its function in White collar complex (WCC) formation. In contrast, the induction of vvd expression by the WCC does not seem to be crucial in this process. Additionally, we found that WC-1 and WC-2 not only act as a complex, but also have individual functions upon growth on cellulose.

Conclusions

Genome wide transcriptome analysis of photoreceptor mutants and evaluation of results by analysis of mutant strains identified several candidate genes likely to play a role in light modulated cellulase gene expression. Genes with functions in amino acid metabolism, glycogen metabolism, energy supply and protein folding are enriched among genes with decreased expression levels in the wc-1 and wc-2 mutants. The ability to properly respond to amino acid starvation, i. e. up-regulation of the cross pathway control protein cpc-1, was found to be beneficial for cellulase gene expression. Our results further suggest a contribution of oxidative depolymerization of cellulose to plant cell wall degradation in N. crassa.

Background

Light is one of the most important environmental cues for almost all living organisms. The daily rhythms of light and darkness and concomitant alterations in temperature, humidity and activity require metabolic adaptation for optimal use of cellular resources. Transmission of the light signal is accomplished by photoreceptor proteins, which are regulated and modified in response to light [1]. The filamentous fungus Neurospora crassa is one of the best studied organisms in this respect and has become a model system for understanding the light response and circadian rhythms [2-5]. Photoreceptors and their cognate signalling cascades have numerous regulatory targets and impact almost every aspect of physiology in N. crassa and other filamentous fungi [6-9]. The N. crassa photoreceptors White Collar 1 and 2 (WC-1 and WC-2) are transcription factors of the fungal GATA zinc finger family [5,10,11]. A complex of these proteins - the White Collar Complex (WCC) - binds to consensus GATA and LRE (light response) elements within the promoters of light regulated genes [7,12-14]. Both photoreceptors contain PAS domains, which are required for homo- and heterodimerization [10,15,16]. Strains containing loss-of-function mutations in wc-1 or wc-2 are largely blind and circadian rhythmicity is perturbed, although a minor residual response to light has been observed [8]. Genes under the control of the WCC can be circadian only, light responsive only, circadian and light responsive or neither light responsive nor circadian [17]. WC-1 and WC-2 also have individual functions besides acting as a component of the WCC [1,18,19]. VVD, the third photoreceptor of N. crassa, is responsible for adaptation to light [20-22] and interacts with the WCC to alter light and clock responses [23]. VVD also contains a PAS domain and acts as a negative regulator following the initiation of light response; Δvvd mutants show an enhanced and prolonged response to light [8,20]. Recently, VVD was found to serve as a molecular memory of the brightness of the preceding day and to discriminate between high and low light, which is especially important for the correct reaction to moonlight [24].

Recently, it has been shown that light modulates transcription of cellulase genes in the industrially important species Trichoderma reesei (anamorph of Hypocrea jecorina) [25].

Subsequent studies aimed at the elucidation of the interconnection of light signalling and cellulase gene expression in T. reesei revealed that orthologs of wc-1 and wc-2 (brl1 and brl2, respectively) and components of heterotrimeric G-protein signalling are involved in light-modulated regulation of cellulase gene expression [25-29]. Light also plays a role in the sulphur requirement for growth on cellulose in T. reesei [30].

Transcriptome and secretome analysis highlighted the value of N. crassa as a model system for understanding the production of plant cell wall degrading enzymes [31]. However, the effect of light on enzyme production has not previously been evaluated in this fungus. In this study, we determined whether light regulation of transcription and production of plant cell wall degrading enzymes also occurs in N. crassa and how this regulation is accomplished. We show that strains containing deletion in genes encoding the photoreceptors WC-1, WC-2 and VVD influence cellulase gene expression. From transcriptional profiling data, we evaluate the phenotype of a number of mutants in genes that showed differential regulation in wc-1, wc-2 and vvd mutants as compared to wild type, revealing that these photoreceptors and the light response are involved in a complex adjustment of physiological processes and composition of the secreted enzyme mixture.

Results

Analysis of cellulolytic activity and growth of Δwc-1, Δwc-2 and Δvvd mutants on cellulose

Measurements of endoglucanase activity revealed that enzyme production is apparent after 16 hrs of growth when N. crassa is cultivated in liquid culture at 25°C with 2% (w/v) crystalline cellulose (Avicel) in light (Figure 1). The Δwc-2 mutant showed a significantly higher level of endoglucanase activity at 28 and 40 hrs (1.15 ± 0.25 U/ml and 1.42 ± 0.08 U/ml) as compared to wild type (WT) (0.34 ± 0.08 U/ml and 0.73 ± 0.16 U/ml). However, after 5 days of cultivation on Avicel, the Δvvd mutant showed significantly higher endoglucanase activity than WT (11.74 ± 0.73 U/ml versus 5.29 ± 0.36 U/ml), while the Δwc-1 and Δwc-2 mutants showed the lowest activity among all strains (~3 U/ml) (Figure 1). Significant differences in the amount of secreted protein among the mutant strains versus wild type were not apparent (Figure 2A, Additional file 1: Figure S1), nor were major defects in hyphal extension rates (Additional file 1: Figure S2). Whereas the Δwc-1, Δwc-2 and Δvvd strains initially produced more biomass than the wild type (28 hrs), this difference was compensated at the later time point (40 hrs) (Figure 2B). Analysis of germination of Δwc-1, Δwc-2 and Δvvd upon growth on cellulose did not reveal a statistically significant difference as compared to wild type germination rates or frequency (data not shown). The specific cellulase activity produced by the mutant strains was comparable to the WT during active growth (28 hrs) and increased with continued cultivation (40 hrs) (Figure 2C).

thumbnailFigure 1. Secretion of cellulolytic enzymes by wild type, Δwc-1, Δwc-2 and Δvvd mutants. Azo- CMC cellulase (endoglucanase) activity in culture filtrates of N. crassa wild type (WT) and mutant strains upon growth in Birds-medium with 2% (w/v) Avicel cellulose as sole carbon source at 16, 28 and 40 hrs versus 5 days. The flasks show examples for degradation of Avicel in wild type (left) and Δvvd mutant (right) strains after 40 hours of growth. At this timepoint the cellulosic substrate was almost or completely degraded by the strains used in this study. Error bars show standard deviation of at least three biological replicates.

thumbnailFigure 2. Growth characteristics of wild type, Δvvd, Δwc-1 and Δwc-2 mutants. (A) extracellular protein production. The wild type strain secreted 42.7 ± 5.1 mg/L of protein after 28 hours and 51.7 ± 2.4 mg/L of protein after 40 hours of cultivation on cellulose. (B) NaOH soluble protein of mycelia reflecting biomass accumulation of WT and mutant strains after 28 or 40 hours of growth on cellulose, relative to WT at the respective time point. (C) Specific Azo- CMC cellulase activity (cellulase activity per biomass) after 28 or 40 hours of growth on cellulose relative to wild type at the respective time point. Error bars show standard deviation of at least two biological replicates.

Additional file 1. Supplementary information. Figure S1 - Profile of total secreted proteins in wild-type and photoreceptor mutants, Figure S2 - Analysis of hyphal extension rates, Figure S3 - Experimental design for microarray analysis of photoreceptor strains, Figure S4 - Transcript abundance of major cellulase genes upon cultivation on cellulose, Figure S5 - Transcript abundance of selected cellulse and hemicellulase genes upon growth on cellulose, Table S1 - Genes with a putative function related to fenton chemistry and in oxidative depolymerization of cellulose.

Format: PDF Size: 649KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Genome wide analysis of transcriptional regulation by photoreceptors on cellulose

The Δwc-1, Δwc-2 and Δvvd mutants showed alterations in cellulase activity (above), suggesting that WC-1, WC-2 and VVD might affect transcription of plant cell wall degrading enzymes. To correlate transcriptional patterns of cellulase gene expression with cellulase activity, we performed qRT-PCR of cellulolytic genes in wild type and mutant strains cultivated under the identical conditions as used for analysis of cellulase production. Taking in account both time points (28 and 40 hours), a significantly lower expression level of the N. crassa cbh1/cel7a homologue (NCU07340; cbh-1) and the cbh2/cel6a homologue (NCU09680; cbh-2) was observed in the Δwc-1 and Δwc-2 mutants (p-values < 0.0005; Figure 3), while expression levels of cbh-1 and cbh-2 in Δvvd were comparable to the WT strain (p-values 0.934 and 0.284, respectively). To further explore how mutations in the photoreceptor genes (wc-1, wc-2 and vvd) globally affect transcriptional regulation of lignocellulolytic genes in N. crassa, we performed expression profiling experiments using full genome oligonucleotide microarrays constructed for N. crassa [32,33]. Transcriptional profiles were assessed in the mutants and wild type at two time points, when active cellulase biosynthesis was occurring (28 and 40 hrs). We chose a closed circuit design for our microarray analysis (Additional file 1 Figure S3) as it allows for robust statistical analyses [34]. Statistical support was detected for transcript levels for a total of 6154 genes (out of ~10,000) (Additional file 2: Dataset 1).

thumbnailFigure 3. Transcript abundance of NCU07340 (encoding CBH-1) and NCU09680 (encoding CBH-2/GH6-2) upon cultivation on cellulose. qRT-PCR data were normalized to transcript levels in wild type after 28 hours and are given for WT after 28 (WT 28) or 40 (WT 40) hours, Δvvd (VVD28 or VVD40), Δwc-1 (WC1 28 or WC1 40) and Δwc-2 (WC2 28 or WC2 40). l6e (NCU02702), a gene encoding a ribosomal protein, was used as control. Analysis was done in triplicates and error bars show standard deviation of these triplicates.

Additional file 2. Dataset 1. Microarray data of N. crassa wild-type and photoreceptor mutants upon growth on cellulose.

Format: XLS Size: 11.8MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Analysis of transcriptional patterns of WT and the photoreceptor mutants (Δwc-1, Δwc-2 and Δvvd) revealed 1718 genes to be significantly differentially transcribed more than 1.5 fold in at least one mutant strain as compared to WT (619 genes up-regulated, 1098 down-regulated). A comparison of genes that showed increased or decreased relative expression levels in the Δwc-1, Δwc-2 and Δvvd mutants revealed that a larger proportion of genes overlapped between the Δwc-1 and Δwc-2 mutants, especially in the gene set that showed decreased expression levels in the mutants (Figure 4). These data are consistent with the role of the WCC functioning in transcriptional activation. For example, WCC directly regulates photoinduction of al-1 (NCU00552, a phytoene desaturase; [35]) and al-2 (NCU00585, a phytoene synthase; [36]), which are components of the carotenoid biosynthetic pathway; wc-1 and wc-2 mutants show reduced levels of al-1 and al-2 [5,10,11], while vvd mutants show increased expression levels [20,37]. Consistent with these data, we observed decreased expression levels for al-1 and al-2 in Δwc-1 and Δwc-2 mutants and increased expression levels in Δvvd (Additional file 2: Dataset 1; Additional file 3: Dataset 2). In addition, transcription data of cbh-1 (NCU07340) and cbh-2 (NCU09680) are in accordance with the data revealed by qRT-PCR, which are more precise than microarray data (Additional file 1: Figure S4).

thumbnailFigure 4. Comparison of gene expression patterns among the three photoreceptor mutants. (A) Venn diagram showing the overlap among genes that showed statistically significantly increase in expression levels in the respective mutant strain(s) compared to WT. (B) Venn diagram showing the overlap among genes that showed statistically significantly decrease in expression levels in the respective mutant strain(s) compared to WT.

Additional file 3. Dataset 2. Functional category analysis of genes up- or down-regulated in photoreceptor mutants.

Format: XLS Size: 387KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Gene clustering of microarray expression profiles and functional category enrichment

Hierarchical clustering showed that the 1718 significantly regulated genes fell into 7 clusters (Figure 5; Additional file 4: Dataset 3). Cluster 1 (461 genes) contained genes that showed reduced expression levels in both Δwc-1 and Δwc-2 at 40 hrs relative to both WT and Δvvd. Functional category analysis showed an enrichment for metabolic genes in this cluster, especially of those involved in amino acid metabolism (P-value 6.24e-22), but also genes related to C-compound and carbohydrate metabolism (P-value 5.91e-12), including two cellulases of glycoside hydrolase (GH) family 7 (NCU04854 and cbh-1, NCU07340), three GH family 61 members (NCU02344, NCU02916 and NCU05969), three hemicellulases (GH family 10, NCU08189; GH family 43, NCU09652 and GH family 51, NCU02343), two β-glucosidases (NCU00130 and NCU07484) and a D-xylose reductase (NCU08384). This cluster also included early light response genes (ELRG) identified in expression profiling experiments performed under sucrose conditions [8], such as genes involved in the biosynthesis of photoprotective pigments (lipid, fatty acid and isoprenoid metabolism), vitamins, cofactors and prosthetic groups, as well as genes involved in secondary metabolism, osmotic and salt stress response. Genes responsible for energy supply, including genes encoding enzymes important in glycolysis and the tricarboxylic acid cycle, were particularly enriched in this group of genes with reduced expression levels in both Δwc-1 and Δwc-2 mutants, as were genes with functions in protein synthesis, transport and stress response.

thumbnailFigure 5. Identification of genes with similar transcription profiles in the strains used for this study by hierarchical clustering of microarray data. (A) A total of 1718 genes were hierarchically clustered based on their expression profiles across all strains and time points and seven clusters were identified. Normalized values are plotted as determined by HCE3.5 (B) Average expression profiles (schematic representation) of genes assigned to each cluster.

Additional file 4. Dataset 3. Functional category analysis of genes assigned to coregulated clusters.

Format: XLS Size: 325KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The second cluster (cluster 2; 417 genes) contained genes that showed increased expression levels in Δvvd, but with decreased expression levels in the white-collar mutants (Figure 5). This cluster showed significant enrichment of genes involved in C-compound and carbohydrate metabolism (P value = 3.54e-14). Genes in this category included the GH6 family cellulases (NCU09680 and NCU07190), a GH7 family cellulase (NCU05057), 6 GH61 family cellulases (NCU00836, NCU02240, NCU03328, NCU07760, NCU07898, NCU08760) as well as four putative acetyl xylan esterases (NCU00710, NCU04870, NCU05159 and NCU09664) and one cellobiose dehydrogenase (NCU00206). Additionally, a slight enrichment of genes involved in oxygen and radical detoxification (P value = 5.70e-4), such as the catalase NCU00355, was observed.

Genes assigned to cluster 3 (415 genes) showed a large increased expression level specifically in the Δwc-1 mutant at 40 hrs (Figure 5). This cluster was predominantly enriched in genes involved in phosphate metabolism (P value = 9.68e-04), energy metabolism (P value = 7.49e-05), ATP binding (P value = 7.76e-06) and stress response (P value = 1.64e-04). Genes assigned to cluster 4 (320 genes) showed an increase in expression level at 40 hrs, particularly in the Δwc-2 mutant (Figure 5). The functional category of polysaccharide metabolism was significantly enriched in this cluster, albeit to a lower extent than clusters 1 and 2 (P value = 1.52e-4). Interestingly, this group comprised three GH 2 family genes (NCU06781, NCU07253 and NCU08909), which encode a β-galactosidase, a β-mannosidase and a β-glucuronidase, respectively. Genes included in clusters 5 and 6 (32 and 21 genes respectively) did not show an enrichment for any functional category, while genes in cluster 7 (51 genes) showed a strong decrease in transcript abundance in the Δvvd mutant at 40 hrs and were enriched in the functions of respiration (P value = 3.55e-5) and transport of substrates (P value 6.62e-4).

Specific regulation in individual photoreceptor mutants

To obtain detailed information on the specific function of the individual photoreceptors under cellulolytic conditions, we analyzed and compared expression profiles for each mutant individually to WT (Table 1 Additional file 3: Dataset 2). In the Δwc-1 mutant, six cellulase genes (NCU02344, NCU04854, NCU05969, NCU7340, NCU7760, NCU9680), four hemicellulase genes (NCU02343, NCU05924, NCU08189, NCU09652) as well as three putative β-glucosidase genes of GH families 1 and 3 (NCU00130, NCU04952, and NCU07487) were significantly down-regulated (P values < 0.005; Figure 6A, B). Additionally, significantly reduced expression levels for genes involved in amino acid metabolism (P value = 8.13e-13), translation (P value = 8.17e-04), protein folding and stabilization (P value = 1.19e-06) as well as energy supply (P value = 5.36e-21) (Additional file 3: Dataset 2) were observed. Among genes with increased expression level in the Δwc-1 mutant, only an enrichment for C-3 compound metabolism was observed (P value = 3.84e-04).

Table 1. Transcript levels of cellulase, hemicellulase and beta-glucosidase genes in photoreceptor mutants

thumbnailFigure 6. Transcript abundance of cellulase and hemicellulase genes in wild type and mutant strains upon growth on cellulose. (A) Genes encoding cellulases for which their contribution to cellulase activity has been evaluated [31], and this study. (B) Genes encoding hemicellulases for which their contribution to cellulase activity has been evaluated [31]. The suffix (P) indicates positive effect on efficiency of the secreted cellulase mixture, (N) indicates a negative effect on the secreted cellulase mixture. (C) Hierarchical clustering of gene regulation patterns in wild type and mutant strains as analyzed using HCE3.5 with the average linkage (UPGMA) method and the Pearson correlation coefficient as distance measure.

The decreased expression levels of genes involved in plant cell wall degradation observed in the Δwc-1 mutant were even more pronounced in the Δwc-2 mutant (Table 1), with additional genes encoding cellulolytic enzymes showing decreased expression levels, including the predicted cellulases NCU00836, NCU01050, NCU02240, NCU02916, NCU03328, NCU05057 and NCU07190 as well as the hemicellulases NCU05955 and NCU07351 (Figure 6A, B; Additional file 1: Figure S5), for a total of 15 cellulase genes (out of 23 predicted;) and 8 hemicellulase genes (out of 19 predicted) [31]. From the cellulases known to be secreted into the medium [31] and hence most likely to contribute to the efficiency of the cellulolytic enzyme mixture of the mutants, two were significantly down-regulated in both Δwc-1 and Δwc-2 (NCU07340 and NCU09680), which correspond to cbh-1 and cbh-2, respectively (Figure 3). Deletion of either cbh-1 or cbh-2 affects cellulolytic activity [31]. Of the four secreted hemicellulases (NCU05924, NCU05955, NCU07326 and NCU08189 [31]), two (NCU05924, NCU08189 encoding an endo-1,4-beta xylanase) were found to be significantly down-regulated in both white collar mutants, while NCU07326 (encoding a glycoside hydrolase family 43) and NCU05955 were down-regulated only in Δwc-2. In Δvvd, NCU05955 and NCU07326 were up-regulated. However, up to ~72 hr of growth, the spectrum and amount of protein secreted into the culture filtrate did not show major defects in the Δwc-2 mutant, even though gene expression profiles for many secreted hydrolases were reduced (Figure 2A, Additional file 1: Figure S1).

In addition to an enrichment for genes within the functional category C-compound and carbohydrate metabolism that showed decreased expression levels in Δwc-2 (P-value 3.67e-25), expression levels for genes encoding proteins involved in protein synthesis (P value = 2.45e-35), energy supply (P value = 1.31e-19) and amino acid metabolism (P value = 6.44e-14) were decreased. Genes involved in secretion were not affected in expression level in the white collar mutants, in agreement with earlier data on cellulase expression in T. reesei, for which also no enhancement of secretion was found under cellulase inducing conditions or in high performance mutants [38].

Among the genes that showed increased expression levels in the Δwc-2 mutant, genes with functionalities of lysosomal and vacuolar protein degradation (P-value 3.03e-4) as well as polysaccharide degradation (P-value 4.95e-5) were significantly enriched. In addition, a few genes related to polysaccharide metabolism showed increased expression levels, for example GH family 2 genes NCU06871, NCU07253 and NCU08909, a glycoside transferase NCU08226 and a pectate lyase NCU09791; no predicted cellulase or hemicellulase genes were identified.

In the Δvvd mutant, the gene set that showed 1.5 fold significant increase in expression level included those responsible for C-compound and carbohydrate degradation, such as four cellulases (NCU00762, NCU02240, NCU02344 and the most highly produced endoglucanase in N. crassa - NCU07190) and three hemicellulases (NCU02855, NCU05924 and NCU07326) (Table 1). Of these genes, NCU00762 [39] encodes a secreted endoglucanase (gh5-1) [31,40]; deletion of NCU00762 significantly reduced endoglucanase activity in N. crassa [31].

Among the down-regulated genes in the Δvvd mutant, enrichment for functions in energy supply (P value 9.65e-05) and phosphate metabolism (P value 5.81e-04) was observed. With respect to genes within C-compound and carbohydrate metabolism, we also found a decrease in transcription of the genes involved in transport (for example of carboxylic acid, glycerol or hexoses; P value = 5.78e-03) and metabolism of energy reserves (for example glycogen synthase; P value = 1.31e-03). Thus, the expression patterns observed in the Δwc-1 and Δwc-2 mutants are more similar to each other, especially in the reduced expression levels of genes within the C-compound and carbohydrate metabolism functional category, while an opposite pattern was observed in the Δvvd mutant. This conclusion is also reflected upon hierarchical clustering of expression patterns (Figure 6C). Additionally this analysis showed the increasing importance of WC-1 and WC-2 after longer cultivation times on cellulose. However, our transcriptional data suggests that the function of WC-1 and WC-2 is not exerted solely as the WCC towards the transcriptional response to growth on cellulose, but that WC-2, especially, has a more important effect on expression patterns and physiology than WC-1.

Identification of genes regulated by WC-1 or WC-2 versus the WCC

Because of the expression differences in the Δwc-1 or Δwc-2 mutants, we evaluated whether genes sets were regulated (directly or indirectly) by the WCC-complex (similar regulation in both mutants) or regulated (directly or indirectly) by WC-1 or WC-2 individually, but not by the WCC (Additional file 5 Dataset 4). Direct targets of WCC have previously been identified by chromatin-immunoprecipitation under sucrose grown conditions [7]. We identified a 27-gene overlap (in the total dataset of 290 genes showing differences in expression level in both Δwc-1 and Δwc-2 mutants versus wild type) to be direct targets of WCC (Additional file 5: Dataset 4). Besides several metabolic genes (a chitinase NCU07035, a GH family 2 protein NCU08909 and a pectate lyase NCU09791), these overlapping targets include vvd [20,41], the carbon catabolite repressor gene cre-1 (NCU08807), two genes involved in conidiation (csp-1 and con-8) and a hydrophobin gene (ccg-2). Cellulases or hemicellulases were not identified among the direct targets of WCC upon growth on sucrose [7].

Additional file 5. Dataset 4. Data on specific regulation by WC-1 and/or WC-2.

Format: XLS Size: 322KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Our data suggest that WC-1 and WC-2 may also have independent functions. We consequently analyzed which functionalities were affected individually in expression profiles from the wc-1 or wc-2 photoreceptor mutants but not by the WCC. Within the WCC regulated gene set, 221 genes showed decreased expression levels, while 69 genes showed increased expression levels (Figure 4; Additional file 5: Dataset 4). In the Δwc-1 mutant, 105 genes showed decreased expression levels and 216 genes showed increased expression levels. In the Δwc-2 mutants, 607 genes showed reduced expression level, while 203 genes showed increased expression levels (Figure 4; Additional file 5: Dataset 4). Metabolic functions including C-compound and carbohydrate metabolism (Δwc-1 P-value = 2.17e-08; Δwc-2 P value = 5.00e-09) were strongly down-regulated in these strains individually as were genes contributing to energy supply (Δwc-1 P-value = 1.89e-07; Δwc-2 P value = 7.54e-08). These functionalities were also enriched in putative WCC-targets, but the individual regulation of components of these pathways indicates additional and distinct roles of WC-1 and WC-2 in these processes.

In the Δwc-2 mutant, the unfolded protein response pathway was significantly down-regulated (P value = 1.66e-04). The most striking difference, however, was the considerable enrichment of the functionalities of ribosome biogenesis, ribosomal proteins and translation (P-values < 1e-33) among genes that showed decreased expression levels in the Δwc-2 mutant only. In the gene set that showed increased expression levels in the Δwc-2 mutant, nine genes involved in polysaccharide metabolism (which are not increased in Δwc-1) were identified, including a gene encoding a predicted glycosyl transferase family 2 enzyme (NCU08226) and a putative GH16 enzyme (NCU05789) (Table 2). Both genes are members of cluster 5 (Figure 5). We tested specific cellulase activity in strains carrying deletions of NCU08226 or NCU05789 (Figure 7). However, both mutants were not significantly different than WT, although lack of NCU08226 caused significantly increased biomass production.

Table 2. Genes selected for analysis with respect to cellulase expression and deletion strains

thumbnailFigure 7. Analysis of selected mutant strains lacking genes with expression profiles or functions of interest for cellulase gene expression. Strains as described in Table 2 are presented according to the clusters assigned to the respective gene. Designation of strains is valid for both columns. Cultivation of strains was performed on Avicel cellulose (2% w/v) for 40 hrs. Specific cellulase activity (cellulase per biomass) and biomass content were determined and normalized to WT. Every measurement was done with at least two biological and two technical replicates. Eight biological replicates were considered for wild type. Error bars show standard deviation of the measurements considering all biological replicates.

Impact of WC-1, WC-2 and VVD on known cellulase regulating genes

The carbon catabolite repression factor cre-1 [42,43] is subject to light-regulation by the WCC [7,8], as is the cross pathway control regulator cpc-1 [32], which plays an important role in regulation of amino acid metabolism. In T. reesei, cellulase regulation is affected by homologues of the HAP-complex [44], the transcriptional activator xyr1 [45] and repressor ace1 [46]. While neither the HAP-complex homologues nor the xyr1 homologue showed altered transcriptional patterns in the N. crassa photoreceptor mutants as compared to wild type, an ace1 homolog (NCU09333) showed increased expression levels in the Δwc-2 mutant after 40 hrs. Similarly, transcriptional levels of cpc-1 and cre-1 increased in both white-collar mutants after 40 hrs (Additional file 2: Dataset 1).

Genomic and functional analyses of the high cellulase producing T. reesei mutant strain RutC30 showed that the lack of functional CRE1 is an important factor in increased cellulase activity [47-49]. In N. crassa, a deletion of cre-1 (NCU08807 [43]) resulted in an increase of 63 ± 14% in biomass-specific cellulase activity compared to WT (Figure 7). Of the 72 genes altered or missing in RutC30 [47-49], 12 were detected in our analysis, of which 8 showed at least a 1.5-fold difference in regulation in at least in one mutant strain. Among them were two genes belonging to cluster 3 (increased expression level in Δwc-1) encoding a MFS multidrug transporter, NCU00306 and a transcription factor, NCU00289 and one gene in cluster 4 (increased expression level in Δwc-2) encoding a predicted nitrilase, NCU01838. The remaining 5 genes (NCU02771, uroporphyrinogen decarboxylase; NCU01332, a vacuolar membrane ATPase; NCU06294, hypothetical protein; NCU06961, an exopolygalacturonase belonging to glycosyl hydrolase family 28 and NCU07129, a predicted amino acid permease) showed increased expression levels in the Δvvd and decreased expression levels in the white-collar mutants (cluster 2).

Interestingly, we found that transcription of NCU05137, encoding a highly conserved non-anchored cell wall protein (NCW-1), the deletion of which causes significantly enhanced secretion of cellulolytic enzymes [31,50], was co-regulated with the known carbon catabolite repression protein encoding cre-1 (NCU08807) and the cross-pathway control protein encoding cpc-1 (NCU04050), which impacts the expression of a large number of target genes [32]. These observations suggest that NCU05137 encodes an auxiliary secreted/cell wall protein important for modulation of extracellular cellulose degrading capacity.

Genes likely to contribute to altered cellulose degradation in photoreceptor mutants

Previously, 14 mutants in genes encoding proteins secreted under cellulolytic conditions were evaluated for cellulase activity [31]; supernatants from most mutants showed near WT cellulolytic activity, with the exception of Δcbh-1 (NCU07340; cellobiohydrolase), Δgh6-2 (NCU09680; cellobiohydrolase), Δgh5-1 (NCU00762; endoglucanase), Δgh3-4 (NCU04952; β-glucosidase) and ncw-1 (NCU05137; cell wall associated protein of unknown function [50]). We selected mutants in three additional genes encoding GH family 61 enzymes (NCU02240, NCU07898 and NCU02344) and one GH family 11 hemicellulase (NCU02855) belonging to cluster 2 (up in Δvvd) for further analysis; no significant alteration in specific cellulase gene expression was observed for mutants in any of the three GH 61 family cellulases (Figure 7). However, deletion of NCU02855 (endo-1,2-beta xylanase A) resulted in considerably increased growth on cellulose (Figure 7) and cellulase secretion, albeit specific cellulase activity was not increased. Among other genes assigned to cluster 2, we also identified acetyl xylan esterases, which have been reported to enhance the performance of T. reesei CBH1/Cel7a [51]. However, a strain carrying a deletion of NCU05159, which encodes an acetyl xylan esterase, did not significantly affect specific cellulase activity, while a mutant in NCU10045 (encoding pectin esterase), showed only a slight increase in specific cellulase activity (Figure 7).

Genes encoding hydrophobins have been reported to be regulated by photoreceptors in T. atroviride [52] and by light in T. reesei [28,53]. Although hydrophobins are not known to have a major influence in cellulase efficiency [54], deletion of the strongly regulated NCU08192 (related to the M. grisea hydrophobin MHP1 [55], cluster 4), that showed increased expression levels in both white collar mutants, unexpectedly had significantly increased specific cellulase activity (Figure 7). A strain containing a mutation in another hydrophobin-like protein belonging to cluster 2 (ccg-14/NCU07787), encoding a ceratoplatanin-like small protein [56,57] displayed cellulolytic activity indistinguishable from WT (Figure 7). An additional gene in cluster 2, NCU02663, encoding an L-lysine 2,3 aminomutase, showed an increase in relative expression of ~50-fold in the Δvvd mutant compared to the white-collar mutants. However, lack of NCU02663 (FGSC16318), did not result in a significant alteration in cellulase activity (Figure 7).

From cluster 3 (genes with increased expression level in Δwc-1), two putative heatshock proteins (NCU00104 and NCU07232) and a gene encoding an uncharacterized cell wall bound protein [50] (NCU00716), were selected for further analysis. Strains containing deletions of NCU00104 and NCU07232 displayed WT cellulase activity. However, deletion of the hypothetical protein encoding NCU00716 caused a significant increase in specific cellulase activity (Figure 7). Since this gene was clearly down-regulated in Δvvd and Δwc-2, a contribution to enhanced cellulase efficiency in these strains is possible.

Only three genes were found to have significantly reduced expression levels in Δwc-1, Δwc-2 and Δvvd, albeit in part less than 1.5 fold. One strongly down regulated gene encoded the alpha type peptide pheromone precursor PPG-1 (ccg-4; NCU02500, cluster 4); a Δccg-4 mutant showed a slight decrease in cellulase activity (Figure 7). Two additional genes were predicted to be involved in glycogen metabolism: NCU02797, encoding a predicted UTP glucose-1-phosphate uridylyltransferase and NCU06687, encoding glycogen synthase. Since a link between glycogen metabolism and cellulase production has been proposed in T. reesei [27,58], we evaluated the cellulolytic phenotype of a ΔNCU06687 mutant (cluster 1); a significant reduction in both biomass and cellulase activity was observed (Figure 7). To further test for the relevance of glycogen metabolism in cellulase production, we also analyzed the cellulolytic activity of a strain lacking the gene for the catabolic enzyme glycogen phosphorylase gph-1 (NCU07027) (cluster 4). Similar to the glycogen synthase mutant (ΔNCU06687), specific cellulase activity was significantly reduced in the ΔNCU07027 mutant (Figure 7).

Co-regulation of genes with cbh-1 (NCU07340)

The glycoside hydrolase family 7 cellulase CBH-1 (NCU07340) is a homologue of T. reesei Cel7a/CBH1, which is the major cellulase in this fungus [59]. Deletion of cbh-1 causes the most severe effect among cellulase genes in terms of both reduced growth on cellulose and decreased cellulase production in N. crassa [31]. We therefore assessed co-regulation of genes with cbh-1 in WT versus the photoreceptor mutants to interconnect regulation of cellulose degradation and the light response. Hierarchical clustering of expression data revealed 70 genes that showed co-regulation with cbh-1 (Additional file 4: Dataset 3). Among these genes we found one additional cellulase besides cbh-1 (NCU02344), one hemicellulase (NCU09652), one beta-glucosidase (NCU07487), a mannosyltransferase (NCU07338) and the aldose epimerase GAL10 (NCU09705). Interestingly, besides several amino acid biosynthetic genes in this cluster (NCU01300, NCU01666, NCU06724, NCU08409 and NCU09320), we also found a sulphate adenylyltransferase (NCU01985). This gene is required for sulphate metabolism, which was shown to be important for cellulase gene expression in T. reesei [30]. This finding supports the relevance of sulphur metabolism for cellulase gene expression to be conserved also in N. crassa.

Functional category analyses of these 70 genes showed significant enrichment in C-compound and carbohydrate metabolism (P value = 9.61e-04), as expected, but also amino acid metabolism (P value = 9.76e-04), which may play a role in cellulose degradation due to its specific regulation in the photoreceptor mutants. Additionally, we found genes responsible for protein synthesis (P value = 4.91e-04) and translation to be enriched among these 70 genes, which suggests posttranslational regulation and/or modification of cellulose degrading enzymes. The enrichment of mitochondrial functions (P value = 2.12e-04) is in accordance with earlier findings on the importance of the physiological state of the mitochondria for cellulase gene expression [60].

Within this 70-gene cluster, we evaluated cellulose activity in two additional deletion strains, one lacking NCU00829, a ferric reductase and a ΔNCU06724 (glutamine synthetase) mutant. Both strains showed a significant decrease in specific cellulase activity (Figure 7). Interestingly, NCU00829, a putative ferric reductase, also showed an increase in relative expression level in the Δvvd mutant (Additional file 2: Dataset 1), suggesting a contribution of oxidative depolymerization of cellulose in N. crassa, which may be important for the efficient cellulose degradation.

Involvement of oxidative depolymerization and Fenton chemistry in degradation of cellulose by N. Crassa

The transcriptional pattern of cellulose degrading enzymes in the white collar mutants as well as the positive effect of the ferric reductase encoded by NCU00829 on cellulase activity led us to hypothesize that fenton chemistry, in which an oxidative process characterized by the formation of highly reactive hydroxyl radicals by reaction of Fe(II) and H2O2, might contribute to cellulose depolymerization, as has been suggested for the basidiomycete Postia placenta [61]. We therefore analyzed the transcriptional patterns of genes putatively involved in this process such as ferric reductases, multicopper oxidases and GMC oxidoreductases (Additional file 1: Table S1) and selected genes encoding different functionalities (ferric reductase, catalase, quinone oxidoreductase, multicopper oxidase, copper superoxide dismutase) putatively involved in oxidative depolymerization of cellulose or its regulation [61] for further analysis. The selected genes were significantly regulated more than 1.5 fold in at least one of the photoreceptor mutant strains and clustered with genes up-regulated in Δvvd (cluster 2) and included NCU00829 (ferric reductase; see above), NCU00355 (catalase-3), NCU02948 (ncw-4, encoding the type IV quinone oxidoreductase), NCU03013 (predicted copper superoxide dismutase anchored cell wall protein, acw-10) and NCU05113 (multicopper oxidase). Proteins encoded by NCU00355, NCU02948 and NCU03013 were reported to be associated with the cell wall in N. crassa [50]. Specific cellulase activity decreased significantly in strains containing deletions of NCU00355 and NCU03013, while a strain lacking NCU02948 (cluster 2) did not differ significantly from wild type (Figure 7). A strain containing a deletion of NCU05113 did not show a statistically significantly altered specific cellulase activity, but showed a considerably increased biomass formation upon growth on cellulose (Figure 7). Since deletion of three of the five genes led to a decrease in specific cellulase activity, fenton chemistry may take part in the attack on cellulose in N. crassa.

Recent reports showed contributions of a cellobiose dehydrogenase and a GH61 protein to oxidative cleavage of cellulose microfibrils [62-65]. We identified a cellobiose dehydrogenase (NCU00206) belonging to cluster 2 [65,66]. Unfortunately, only a heterokaryotic deletion strain (FGSC16368) was available. To evaluate whether other probable/predicted cellobiose dehydrogenases (CBDH) affected cellulase activity, we analyzed a homokaryotic strain lacking a probable CBDH (NCU05923) and three homokaryotic mutants with carrying deletion of a gene related to CBDH proteins (NCU01873, NCU05595 and NCU08432). However, none of these mutants affected cellulase activity, although the strain carrying a deletion of NCU05923 caused significantly decreased biomass production on cellulose (Figure 7) suggesting a function beneficial for growth on cellulose.

Assessment of the role of amino acid metabolism

Transcriptional profiling analyses suggested that genes involved in amino acid metabolism may be regulated by WC-1 and WC-2 during growth on cellulose. This regulation has already been reported for growth on glucose [8] and may thus represent a general phenomenon correlated to the light response. A study with the T. reesei mutant RutC30 showed that decreased ability to utilize amino acids correlates with increased cellulase gene expression [48]. Additionally, an interrelationship between methionine- and sulphur signalling as well as a clear difference of the respective effect in light and darkness has been shown for T. reesei [30]. Since the response to amino acid starvation is mainly regulated by CPC-1 in N. crassa [32,67], which directly or indirectly targets more than 400 genes [32], we hypothesized that carbohydrate metabolism on cellulose might also be a target of CPC-1. Consequently, we compared the CPC-1 targets [32] with the genes regulated by both WC-1 and WC-2 (290 total). Indeed, 48 genes were found to be targets of both CPC-1 and WCC. Functional analysis of these 48 genes revealed significant enrichment in functionalities in both amino acid metabolism (P value = 1.99e-06) and C-compound and carbohydrate metabolism (P value = 7.10e-04) (Additional file 6: Dataset 5). These data support the hypothesis that the light signalling pathway is interconnected with amino acid metabolism with one of the output targets being cellulase gene expression.

Additional file 6. Dataset 5. Data on CPC-1 target genes regulated by WC-1 and WC-2.

Format: XLS Size: 56KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

To test this hypothesis, we evaluated the cellulolytic phenotype of mutants in genes involved in regulation of amino acid metabolism. The transcription of cpc-1 (NCU04050), [32,67] was enhanced in the Δwc-1 and Δwc-2 mutants. A strain lacking cpc-1 showed significantly decreased cellulase activity and lower biomass formation on cellulose (Figure 7). Similarly, strains carrying deletions of genes identified within the CPC-1 regulon [32], NCU04482 and NCU00365 (both hypothetical proteins) caused a slight reduction in specific cellulase activity (Figure 7). As noted above, a ΔNCU06724 mutant (encoding a glutamine synthetase and co-regulated with NCU07340; Additional file 4: Dataset 3) showed significantly decreased specific cellulase activity (Figure 7), while a strain carrying a deletion of NCU03935, encoding a homoserine dehydrogenase (cluster 1) was not significantly different from WT. Other genes shown to be responsive to amino acid starvation ([32]; Fungal Gene Expression Database: http://bioinfo.townsend.yale.edu/index.jsp webcite) and regulated in the photoreceptor mutants showed a significant decrease in specific cellulase activity (NCU07027, glycogen phosphorylase and NCU00355, catalase, see above and Figure 7). An additional member of this cluster, NCU03753 (glucose repressible protein, ccg-1), showed a slight decrease in specific cellulase activity. These results suggest a link between response to amino acid starvation and modulation of cellulase gene expression in N. crassa as well as in T. reesei.

Discussion

The preferred natural habitats of both Trichoderma and Neurospora are humid tropical and subtropical forests [68,69] and hence utilization of cellulosic substrates is crucial in their life style. Using resources available for N. crassa together with extensive information about cellulose degradation in T. reesei [68,70,71], we aimed in this study to obtain a more detailed understanding of the physiology of cellulase gene expression across fungi. We show that photoreceptors not only modulate cellulase gene expression in T. reesei [25,26,72] but also in N. crassa. Moreover, the resulting effect seems to be largely similar in that transcript levels of the major cellulases tend to be decreased in strains lacking the photoreceptors, while secreted levels of cellulolytic enzymes are increased. Therefore, we conclude that light-modulated expression of cellulolytic enzymes is a conserved phenomenon, most probably providing an evolutionary advantage to filamentous fungi in their natural habitat.

Specific cellulase activity was slightly higher in all three photoreceptor mutants early in utilization of cellulose compared to WT (28 and 40 hrs), but the Δvvd mutant showed significantly more cellulase activity after 5 days of growth on Avicel. These data suggest that, over time, WC-1 and WC-2 promote cellulose utilization, while VVD negatively antagonizes it. Consistent with this hypothesis, the Δwc-1 and Δwc-2 mutants showed significant decrease in expression levels for the major predicted cellulase and hemicellulase genes. In case of Δvvd, both the enrichment of carbohydrate metabolism genes, as well as the beneficial effect of several genes on cellulose degradation, as revealed by mutant analysis, indicate that the reason for enhanced cellulose degradation by the Δvvd mutant lies at least in part in increased expression level of hydrolytic enzymes. From studies of circadian/light regulation by photoreceptors in N. crassa, it has been proposed that a hierarchical cascade of transcription factors regulated by WCC mediate light and circadian rhythm signaling networks, which in turn regulate 5-20% of the genes in the Neurospora genome [7,8,73]. The activation of light-induced gene expression by WCC is transient, so that levels of mRNA decrease in continuous light to reach a steady state level [74]. The attenuation of light-induced gene regulation by WCC is termed photoadaptation and requires VVD [20,37]. Photoadaptation occurs via light-dependent physical interaction of VVD and WC-1 LOV domains [23], resulting in inactivation of the WCC. Thus, in the Δvvd mutant, photoadaptation is not functioning, which results in increased and prolonged expression levels of genes regulated by WCC, including those encoding cellulolytic enzymes and genes involved in energy metabolism (Figure 8).

thumbnailFigure 8. Model for regulation of plant cell wall degrading capacity in N. crassa by WC-1/WC-2/WCC and VVD. VVD, WC-1 and WC-2 negatively affect efficiency of cellulose degradation during early stages of growth. Regulation of cellulolytic enzymes and energy supply is subject to photoadaptation by VVD, while carbon catabolite repression as mediated by CRE-1 is regulated by the WCC via VVD. Independent regulatory functions of VVD influencing cellulolytic efficiency include regulation of oxidative depolymerization, while for WC-1 and WC-2 these functions are most obvious in regulation of specific ribosomes. The WCC moreover specifically targets amino acid metabolism. Since the effect of the light response machinery is likely to be indirect, their output may involve additional target transcription factors.

Induction of cellulolytic genes and enzyme activity requires release from carbon catabolite repression (CCR), which is mediated in part by CRE-1 in filamentous fungi, which directly inhibits expression of genes encoding important plant cell wall degrading enzymes [43,70,75-77]. CRE-1 is also a target of the WCC [7]. In this respect, the negative effect of all three photoreceptors on transcription of cre-1 after 40 hours (Additional file 2: Dataset 1) indicates that their regulatory contribution to CCR is mediated by the WCC via VVD and not (only) subject to photoadaptation (Figure 8). Besides release from CCR, optimal expression of genes encoding plant cell wall degrading enzymes and cellulolytic activity requires the presence of inducer molecules. Our data suggests that the WCC regulates genes involved in plant cell wall degradation indirectly rather than directly, perhaps via activation of a target transcription factor (Figure 8). Consistent with this model is the observation that one of the transcription factors identified as a direct target of the WCC is required for growth of N. crassa on cellulose (Glass et al., unpublished results). In addition, regulatory output of VVD includes also WCC-independent targets, such as genes involved in oxidative depolymerization, which are likely to be important contributors to extracellular cellulase activity.

Our data do not support a model for the regulation of plant cell wall degradation by photoreceptors in which they exert their function - at least in part - via VVD, which is their regulatory target. The effects of lack of VVD considerably differ from those caused by deletion of wc-1 or wc-2. Moreover, VVD, WC-1 and WC-2 also individually have distinct targets, hence suggesting a complex regulation scheme, which can respond to different physiological cues. An interrelationship of the light response pathway with pathways sensing nutritional signals has been shown in T. reesei [9,78]. The distinct function of each of the three photoreceptors could reflect integration of multiple environmental signals with the light signal that enable optimal growth on cellulose.

Chen and coworkers [8] found genes involved in metabolic processes to be predominantly among the late light regulated genes (LLRGs), thus WCC may contribute to maintenance of appropriate energy supply in response to environmental conditions and substrate availability. Notably, energy metabolism, including glycogen metabolism, is positively regulated by all three photoreceptors and is likely to be perturbed in the mutants. The early enhanced extracellular cellulase activity and the considerably lower cellulase activity produced by Δwc-1 and Δwc-2 compared to the wild type after five days of growth (Figure 1) may partly be a reflection of a perturbation of energy metabolism in the Δwc-1 and Δwc-2 mutants. Consistent with this hypothesis, the lack of WC-1 or WC-2 also causes distinctly altered transcription patterns of numerous ribosomal genes and decreased transcription of genes involved in amino acid metabolism. An interconnection between cellulase activity and amino acid metabolism has been shown in T. reesei [48]; our data supports the conservation of this interconnection in N. crassa. The fact that the major regulator of starvation for amino acids, CPC-1, has 443 direct or indirect targets in N. crassa [32] reflects the importance of this regulator; a Δcpc-1 mutant showed reduced ability to grow on cellulose and reduced cellulolytic activity (Figure 7). The integration of a response to amino acid starvation with enhanced production of enzymes degrading extracellular substrates could be a useful strategy to acquire the missing nutrients present in plant cells besides cellulose. Interestingly, genes encoding ribosomal proteins are coordinately regulated during carbon source shifts [79-81] and this response is not mediated by CRE-1 [82]. Since photoreceptors have been shown to control sensitivity to and utilization of carbon sources [83] the significant enrichment of ribosomal gene regulation likely reflects altered or decreased sensitivity to the extracellular carbon source. The slightly enhanced cellulase activity during early growth of Δwc-1 and Δwc-2 may therefore be a consequence of the attempt to compensate for this defect by adjustment of the ribosome pool [84,85].

Our study shows that the adjustment of the expression of cellulolytic genes triggered by light and the photoreceptors is sophisticated (Table 2). We evaluated the celluloytic phenotype of a number of genes with either predicted roles in plant cell wall degradation, or which were differentially regulated in the mutants. The majority of the mutants did not display cellulolytic activity that was significantly different from WT. However, it has been shown that deletion of a single cellulase gene does not necessarily causes decrease the capacity to degrade cellulose [31] and that additional genes contribute to the efficiency of the cellulase mixture [40]. In fact, Tian and coworkers [31] showed that deletion of certain cellulase genes can result in an increase in cellulase activity, suggesting compensatory mechanisms are occurring. We predict that the complex machinery of N. crassa allows for a sophisticated fine-tuning plant cell wall degrading enzyme repertoire in response to environmental conditions, which is also reflected in the regulatory patterns revealed in the photoreceptor mutants.

Although the transcriptional regulatory effects on predicted cellulolytic genes in the Δwc-1 and Δwc-2 mutants were similar, specific effects were also observed. For example, the Δwc-2 mutant showed decreased transcription of six additional cellulolytic genes (Table 1). The importance of WC-2 for interaction of WC-1 with FRQ [86] may be one reason for the more severe effect on transcriptional regulation that we observed in the Δwc-2 mutant. Similarly, data from T. atroviride showed that over-expression of brl2 (the homologue of wc-2) exerted a stronger and partially contrary effect to over-expression of brl1 (homologue of wc-1) [87]. The fact that the (direct and indirect) targets of WCC upon growth on cellulose as revealed in this study, only show a relatively small overlap with those identified by Smith and coworkers [7] confirms that WC-1 and WC-2 are important factors for the interconnection between light response and regulation of carbon source utilization. These transcription factors can consequently be expected to have carbon source dependent functions in addition to their function in light response and circadian rhythmicity as was also suggested for brl1 and brl2 in T. atroviride [83]. The identification of direct targets of WC-1, WC-2 and the WCC under lignocellulolytic conditions is an important avenue for future studies.

Conclusions

In this study aimed at elucidation of the molecular mechanisms causing light modulated cellulase gene expression, we showed that this phenomenon is conserved between N. crassa and T. reesei, and most probably also in other ascomycete species. In summary, all three photoreceptors of N. crassa execute distinct roles in regulation of cellulase utilization and are hence crucial determinants of the interconnection between light response and carbon metabolism.

Methods

Strains and culture conditions

The wild type laboratory strain, FGSC 2489, used for the N. crassa genome project, as well as the deletion mutants FGSC11556 (Δvvd), FGSC11712 (Δwc-1) and FGSC11124 (Δwc-2) obtained from the Fungal Genetics Stock Center (FGSC) [88] were used for microarray analysis [89]. Additional deletion mutants used for screening were obtained from the FGSC. Strains were inoculated onto slants containing Vogel's minimal medium [90] at 25°C for 7 days in constant light conditions (approx. 300 lux) to suppress synchronous gene expression associated with circadian rhythms. Conidia were harvested with water and inoculated into 50 ml of Bird's medium [91] with 2% (w/v) Avicel cellulose (Sigma) as sole carbon source in 250 ml Erlenmayer flasks at a final concentration of 106 conidia/ml. Cells were grown under constant light at 25°C for 28 or 40 hours at 200 rpm. Growth in race tubes was analyzed using Bird's medium without carbon source (control) and with 1% of sucrose or carboxymethylcellulose (CMC) as carbon source. Measurements were done in triplicates and hyphal extension was measured daily.

Measurements of cellulase activity, biomass and secreted protein

For measurement of specific endo-1,4-β-glucanase (cellulase) activity in culture filtrates of wild type and mutant N. crassa strains, azo-CM-cellulose (article number S-ACMCL; Megazyme, Bray, Ireland) was used according to the manufacturer's instructions. Biomass accumulation after 28 or 40 hrs of growth in the presence of cellulose was determined as described earlier [25]. Briefly, mycelia including residual cellulose was harvested, squeezed dry with Whatman filter paper, frozen in liquid nitrogen, crushed and suspended in 5 ml of 0.1 N NaOH. After sonication for 2 minutes twice, samples were incubated for 3 hours at room temperature and centrifuged at 14000 g for 10 minutes. Protein content in the supernatant was determined using the Bio-Rad Protein Assay (Bio-Rad, Hercules, US), which reflects biomass content. At least two biological and two technical replicates were considered for calculation of results for every analysis.

RNA isolation and qRT-PCR

Strains were grown under the conditions specified with the respective experiments and fungal mycelia were harvested by filtration and frozen in liquid nitrogen. Total RNA was isolated by using TRIzol reagent (Invitrogen Life Technologies), and purified by using the RNeasy kit (QIAGEN) according to the manufacturer's protocols. cDNA preparation for qRT-PCR was done using the RevertAid H- cDNA Kit (No. K1631; Fermentas, Vilnius, Lithuania). qRT-PCR was performed using the iQ SYBR Green Supermix containing fluorescein (Bio-Rad, Hercules, USA) and the iCycler iQ5 Real-time PCR System (Bio-Rad). The l6e gene (NCU02702; primers: RT_NC_L6eF (5' CAGAAATGGTACCCTGCTGAGG 3') and RT_NC_L6eR (5' GCGGATGGTCTTGCGG 3')) was used as reference gene [29,92]. Primers used for analysis of cbh-1 (NCU07340) were RT_NC_cbh1F (5' CACCACCATCGAACAGCAC 3') and RT_NC_cbh1R (5' CAGTCTTGCCCTCACCGTAG 3') and for cbh-2 (NCU09680) RT_NC_cbh2F (5' CCCATCACCACTACTACC 3') and RT_NC_cbh2R (5' CCAGCCCTGAACACC 3'). Analysis of results was performed as described in [29] and for statistical analyses of expression ratios REST© Software was used with a confidence interval of 95% [93].

cDNA labeling for microarray analysis

For cDNA synthesis and labeling, the Pronto kit (catalog no. 40076, Corning) was used. Independent samples from two biological replicates were used. Briefly, cDNA was synthesized from a mixture containing 10 μg total RNA and oligo(dT)primer, ChipShot reverse transcriptase, and aminoallyl-deoxynucleoside triphosphate and incubated at 42°C for 2 h. The cDNA was purified by using ChipShot membrane columns. The dyes Cy3 and Cy5 (Amersham; catalog no. RPN5661) were incorporated into cDNA by adding Cy3 or Cy5 monofunctional N-hydroxy-succinimide ester dye to the cDNA solution for 1 h at 22°C. The cDNA was subsequently dried under vacuum and used for hybridization.

Hybridization and image acquisition

For this analysis, N. crassa whole genome 70-mer oligonucleotide arrays were used as described earlier [33]. A GenePix 4000B scanner (Axon Instruments, CA) was used to acquire images and GenePix Pro6 software was used to quantify hybridization signals. Low-quality spots were flagged automatically by GenePix software and additionally each slide was expected manually.

Experimental design and data analysis

We chose a closed-circuit design for microarray comparisons of WT and recombinant strains at the two different time-points (Additional file 1: Figure S3). Circuit designs for microarrays are statistically robust and improve resolution in identifying differentially regulated genes compared to designs for microarrays that use a universal reference [94]. Data analysis was performed essentially as described in [31-33]. Expression data has been submitted to GEO (http://www.ncbi.nlm.nih.gov/gds webcite; Accession Number = GSE32871). For further analysis hybridized spots with a mean fluorescence intensity for at least one of the Cy3 and Cy5 dyes that was greater than the mean background intensity plus three standard deviations were scored, if less than 0.02% of the pixels were saturated.

The relative expression data of genes was estimated by the BAGEL software (34). BAGEL explores the ratio for all samples using a Markov Chain Monte Carlo (MCMC) approach in a Bayesian framework, and infers relative expression levels and statistical significance from the parameter values it samples. BAGEL also provides gene-by-gene 95% credible intervals, values in which 95% of the samples from the chain are bounded. Comparisons of gene expression levels in this study that were identified with at least 1.5 fold change at significant differences of the 95% credible level (p-value = 0.05), in WT versus mutant strains (Δwc-1, Δwc-2 and Δvvd) at the conditions tested. Significant distance in gene expression was determined essentially as described in [31,32] using the PERL script "errorbar.pl" for expression differences as described in the text.

These inferred levels of gene expression were clustered using Hierarchical Clustering Explorer 3.5 (HCE; Human-Computer Interaction Lab, University of Maryland, College Park, USA) and applying default parameters (Euclidian Distance model).

The functional catalog FunCat created by MIPS [95] was mined to associate functional annotations with Neurospora genes http://mips.helmholtz-muenchen.de/genre/proj/ncrassa/Search/Catalogs/searchCatfirstFun.html webcite. The statistically significant overrepresentation of gene groups in functional categories relative to the whole genome was determined with the MIPS online FunCat system using hypergeometric distribution for P value calculation http://mips.gsf.de/proj/funcatDB/help_p-value.html webcite. We set the threshold for significant enrichment to P < 0.001.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MS carried out experimental work and data analysis and drafted the manuscript. CG supervised microarray analysis and participated in data analysis. JS and DT participated in cultivation and analysis of cellulase expression. DT performed qRT-PCR analyses. MS and NLG collaborated in designing the study and writing the final version of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Research of MS and DT was supported by the Austrian Academy of Sciences (APART fellowship 11212, Doc-fFORTE fellowship) and grants of the Austrian Research Fund (FWF) to MS (P21072, V152-B20). Research of CT was partially supported by China 973 Program (2011CB707403) and Chinese Academy of Sciences Grant (KSCX2-EW-G-8). Research of NLG was supported by a National Institutes of Health (NIH) Neurospora Program Project Grant GM068087.

References

  1. Corrochano LM: Fungal photoreceptors: sensory molecules for fungal development and behaviour.

    Photochem Photobiol Sci 2007, 6(7):725-736. PubMed Abstract | Publisher Full Text OpenURL

  2. Brunner M, Kaldi K: Interlocked feedback loops of the circadian clock of Neurospora crassa.

    Mol Microbiol 2008, 68(2):255-262. PubMed Abstract | Publisher Full Text OpenURL

  3. Davis RH, Perkins DD: Timeline: Neurospora: a model of model microbes.

    Nat Rev Genet 2002, 3(5):397-403. PubMed Abstract | Publisher Full Text OpenURL

  4. Dunlap JC, Loros JJ: The Neurospora circadian system.

    J Biol Rhythms 2004, 19(5):414-424. PubMed Abstract | Publisher Full Text OpenURL

  5. Dunlap JC, Loros JJ: Neurospora Photoreceptors. In Handbook of Photosensory Receptors. Edited by Briggs WR, Spudich JL. Weinheim: WILEY-VCH; 2005:371-390. OpenURL

  6. Schmoll M: Assessing the relevance of light for fungi: implications and insights into the network of signal transmission.

    Adv Appl Microbiol 2011, 76:27-78. PubMed Abstract | Publisher Full Text OpenURL

  7. Smith KM, Sancar G, Dekhang R, Sullivan CM, Li S, Tag AG, Sancar C, Bredeweg EL, Priest HD, McCormick RF, et al.: Transcription factors in light and circadian clock signaling networks revealed by genome-wide mapping of direct targets for Neurospora WHITE COLLAR COMPLEX.

    Eukaryot Cell 2010, 9(10):1549-1556. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Chen CH, Ringelberg CS, Gross RH, Dunlap JC, Loros JJ: Genome-wide analysis of light-inducible responses reveals hierarchical light signalling in Neurospora.

    EMBO J 2009, 28(8):1029-1042. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Tisch D, Schmoll M: Light regulation of metabolic pathways in fungi.

    Appl Microbiol Biotechnol 2010, 85(5):1259-1277. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Ballario P, Vittorioso P, Magrelli A, Talora C, Cabibbo A, Macino G: White collar-1, a central regulator of blue light responses in Neurospora, is a zinc finger protein.

    EMBO J 1996, 15(7):1650-1657. PubMed Abstract | PubMed Central Full Text OpenURL

  11. Linden H, Macino G: White collar 2, a partner in blue-light signal transduction, controlling expression of light-regulated genes in Neurospora crassa.

    EMBO J 1997, 16(1):98-109. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Carattoli A, Cogoni C, Morelli G, Macino G: Molecular characterization of upstream regulatory sequences controlling the photoinduced expression of the albino-3 gene of Neurospora crassa.

    Mol Microbiol 1994, 13(5):787-795. PubMed Abstract | Publisher Full Text OpenURL

  13. Froehlich AC, Liu Y, Loros JJ, Dunlap JC: White Collar-1, a circadian blue light photoreceptor, binding to the frequency promoter.

    Science 2002, 297(5582):815-819. PubMed Abstract | Publisher Full Text OpenURL

  14. He Q, Liu Y: Molecular mechanism of light responses in Neurospora: from light-induced transcription to photoadaptation.

    Genes Dev 2005, 19(23):2888-2899. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Ballario P, Talora C, Galli D, Linden H, Macino G: Roles in dimerization and blue light photoresponse of the PAS and LOV domains of Neurospora crassa white collar proteins.

    Mol Microbiol 1998, 29(3):719-729. PubMed Abstract | Publisher Full Text OpenURL

  16. Cheng P, Yang Y, Gardner KH, Liu Y: PAS domain-mediated WC-1/WC-2 interaction is essential for maintaining the steady-state level of WC-1 and the function of both proteins in circadian clock and light responses of Neurospora.

    Mol Cell Biol 2002, 22(2):517-524. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Dong W, Tang X, Yu Y, Nilsen R, Kim R, Griffith J, Arnold J, Schuttler HB: Systems biology of the clock in Neurospora crassa.

    PLoS One 2008, 3(8):e3105. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Cheng P, Yang Y, Wang L, He Q, Liu Y: White collar-1, a multifunctional Neurospora protein involved in the circadian feedback loops, light sensing, and transcription repression of wc-2.

    J Biol Chem 2003, 278(6):3801-3808. PubMed Abstract | Publisher Full Text OpenURL

  19. Crosthwaite SK, Dunlap JC, Loros JJ: Neurospora wc-1 and wc-2: transcription, photoresponses, and the origins of circadian rhythmicity.

    Science 1997, 276(5313):763-769. PubMed Abstract | Publisher Full Text OpenURL

  20. Heintzen C, Loros JJ, Dunlap JC: The PAS protein VIVID defines a clock-associated feedback loop that represses light input, modulates gating, and regulates clock resetting.

    Cell 2001, 104(3):453-464. PubMed Abstract | Publisher Full Text OpenURL

  21. Schwerdtfeger C, Linden H: Blue light adaptation and desensitization of light signal transduction in Neurospora crassa.

    Mol Microbiol 2001, 39(4):1080-1087. PubMed Abstract | Publisher Full Text OpenURL

  22. Schwerdtfeger C, Linden H: Vivid is a flavoprotein and serves as a fungal blue light photoreceptor for photoadaptation.

    EMBO J 2003, 22(18):4846-4855. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Hunt SM, Thompson S, Elvin M, Heintzen C: VIVID interacts with the WHITE COLLAR complex and FREQUENCY-interacting RNA helicase to alter light and clock responses in Neurospora.

    Proc Natl Acad Sci USA 2010, 107(38):16709-16714. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Malzahn E, Ciprianidis S, Kaldi K, Schafmeier T, Brunner M: Photoadaptation in Neurospora by competitive interaction of activating and inhibitory LOV domains.

    Cell 2010, 142(5):762-772. PubMed Abstract | Publisher Full Text OpenURL

  25. Schmoll M, Franchi L, Kubicek CP: Envoy, a PAS/LOV domain protein of Hypocrea jecorina (anamorph Trichoderma reesei), modulates cellulase gene transcription in response to light.

    Eukaryot Cell 2005, 4(12):1998-2007. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Castellanos F, Schmoll M, Martinez P, Tisch D, Kubicek CP, Herrera-Estrella A, Esquivel-Naranjo EU: Crucial factors of the light perception machinery and their impact on growth and cellulase gene transcription in Trichoderma reesei.

    Fungal Genet Biol 2010, 47(5):468-476. PubMed Abstract | Publisher Full Text OpenURL

  27. Schmoll M, Schuster A, do Nascimento Silva R, Kubicek CP: The G-alpha protein GNA3 of Hypocrea jecorina (anamorph Trichoderma reesei) regulates cellulase gene expression in the presence of light.

    Eukaryot Cell 2009, 8(3):410-420. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Seibel C, Gremel G, Silva RD, Schuster A, Kubicek CP, Schmoll M: Light-dependent roles of the G-protein subunit GNA1 of Hypocrea jecorina (anamorph Trichoderma reesei).

    BMC Biol 2009, 7(1):58. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Tisch D, Kubicek CP, Schmoll M: New insights into the mechanism of light modulated signaling by heterotrimeric G-proteins: ENVOY acts on gna1 and gna3 and adjusts cAMP levels in Trichoderma reesei (Hypocrea jecorina).

    Fungal Genet Biol 2011, 48(6):631-640. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Gremel G, Dorrer M, Schmoll M: Sulphur metabolism and cellulase gene expression are connected processes in the filamentous fungus Hypocrea jecorina (anamorph Trichoderma reesei).

    BMC Microbiol 2008, 8(1):174. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  31. Tian C, Beeson WT, Iavarone AT, Sun J, Marletta MA, Cate JHD, Glass NL: Systems analysis of plant cell wall degradation by the model filamentous fungus Neurospora crassa.

    Proc Natl Acad Sci USA 2009, 106(52):22157-22162. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Tian C, Kasuga T, Sachs MS, Glass NL: Transcriptional profiling of cross pathway control in Neurospora crassa and comparative analysis of the Gcn4 and CPC1 regulons.

    Eukaryot Cell 2007, 6(6):1018-1029. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Kasuga T, Townsend JP, Tian C, Gilbert LB, Mannhaupt G, Taylor JW, Glass NL: Long-oligomer microarray profiling in Neurospora crassa reveals the transcriptional program underlying biochemical and physiological events of conidial germination.

    Nucleic Acids Res 2005, 33(20):6469-6485. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Townsend JP, Hartl DL: Bayesian analysis of gene expression levels: statistical quantification of relative mRNA level across multiple strains or treatments.

    Genome Biol 2002, 3:71. OpenURL

  35. Schmidhauser TJ, Lauter FR, Russo VE, Yanofsky C: Cloning, sequence, and photoregulation of al-1, a carotenoid biosynthetic gene of Neurospora crassa.

    Mol Cell Biol 1990, 10(10):5064-5070. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Schmidhauser TJ, Lauter FR, Schumacher M, Zhou W, Russo VE, Yanofsky C: Characterization of al-2, the phytoene synthase gene of Neurospora crassa. Cloning, sequence analysis, and photoregulation.

    J Biol Chem 1994, 269(16):12060-12066. PubMed Abstract | Publisher Full Text OpenURL

  37. Shrode LB, Lewis ZA, White LD, Bell-Pedersen D, Ebbole DJ: vvd is required for light adaptation of conidiation-specific genes of Neurospora crassa, but not circadian conidiation.

    Fungal Genet Biol 2001, 32(3):169-181. PubMed Abstract | Publisher Full Text OpenURL

  38. Foreman PK, Brown D, Dankmeyer L, Dean R, Diener S, Dunn-Coleman NS, Goedegebuur F, Houfek TD, England GJ, Kelley AS, et al.: Transcriptional regulation of biomass-degrading enzymes in the filamentous fungus Trichoderma reesei.

    J Biol Chem 2003, 278(34):31988-31997. PubMed Abstract | Publisher Full Text OpenURL

  39. Sun J, Phillips CM, Anderson CT, Beeson WT, Marletta MA, Glass NL: Expression and characterization of the Neurospora crassa endoglucanase GH5-1.

    Protein Expr Purif 2011, 75(2):147-154. PubMed Abstract | Publisher Full Text OpenURL

  40. Phillips CM, Iavarone AT, Marletta MA: A quantitative proteomic approach for cellulose degradation by Neurospora crassa.

    J Proteome Res 2011, 10(9):4177-4185. PubMed Abstract | Publisher Full Text OpenURL

  41. Lewis ZA, Correa A, Schwerdtfeger C, Link KL, Xie X, Gomer RH, Thomas T, Ebbole DJ, Bell-Pedersen D: Overexpression of White Collar-1 (WC-1) activates circadian clock- associated genes, but is not sufficient to induce most light-regulated gene expression in Neurospora crassa.

    Mol Microbiol 2002, 45(4):917-931. PubMed Abstract | Publisher Full Text OpenURL

  42. Nakari-Setälä T, Paloheimo M, Kallio J, Vehmaanperä J, Penttilä M, Saloheimo M: Genetic modification of carbon catabolite repression in Trichoderma reesei for improved protein production.

    Appl Environ Microbiol 2009, 75(14):4853-4860. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Sun J, Glass NL: Identification of the CRE-1 cellulolytic regulon in Neurospora crassa.

    PLoS One 2011, 6(9):e25654. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Zeilinger S, Ebner A, Marosits T, Mach R, Kubicek CP: The Hypocrea jecorina HAP 2/3/5 protein complex binds to the inverted CCAAT-box (ATTGG) within the cbh2 (cellobiohydrolase II-gene) activating element.

    Mol Genet Genomics 2001, 266(1):56-63. PubMed Abstract | Publisher Full Text OpenURL

  45. Stricker AR, Grosstessner-Hain K, Wurleitner E, Mach RL: Xyr1 (Xylanase Regulator 1) regulates both the hydrolytic enzyme system and D-Xylose metabolism in Hypocrea jecorina.

    Eukaryot Cell 2006, 5(12):2128-2137. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Aro N, Ilmen M, Saloheimo A, Penttila M: ACEI of Trichoderma reesei is a repressor of cellulase and xylanase expression.

    Appl Environ Microbiol 2003, 69(1):56-65. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Le Crom S, Schackwitz W, Pennacchio L, Magnuson JK, Culley DE, Collett JR, Martin J, Druzhinina IS, Mathis H, Monot F, et al.: Tracking the roots of cellulase hyperproduction by the fungus Trichoderma reesei using massively parallel DNA sequencing.

    Proc Natl Acad Sci USA 2009, 106(38):16151-16156. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Seidl V, Gamauf C, Druzhinina IS, Seiboth B, Hartl L, Kubicek CP: The Hypocrea jecorina (Trichoderma reesei) hypercellulolytic mutant RUT C30 lacks a 85 kb (29 gene-encoding) region of the wild-type genome.

    BMC Genomics 2008, 9:327. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  49. Strauss J, Mach RL, Zeilinger S, Hartler G, Stoffler G, Wolschek M, Kubicek CP: Cre1, the carbon catabolite repressor protein from Trichoderma reesei.

    FEBS Lett 1995, 376(1-2):103-107. PubMed Abstract | Publisher Full Text OpenURL

  50. Maddi A, Bowman SM, Free SJ: Trifluoromethanesulfonic acid-based proteomic analysis of cell wall and secreted proteins of the ascomycetous fungi Neurospora crassa and Candida albicans.

    Fungal Genet Biol 2009, 46(10):768-781. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Selig MJ, Knoshaug EP, Adney WS, Himmel ME, Decker SR: Synergistic enhancement of cellobiohydrolase performance on pretreated corn stover by addition of xylanase and esterase activities.

    Bioresour Technol 2008, 99(11):4997-5005. PubMed Abstract | Publisher Full Text OpenURL

  52. Rosales-Saavedra T, Esquivel-Naranjo EU, Casas-Flores S, Martinez-Hernandez P, Ibarra-Laclette E, Cortes-Penagos C, Herrera-Estrella A: Novel light-regulated genes in Trichoderma atroviride: a dissection by cDNA microarrays.

    Microbiology 2006, 152(Pt 11):3305-3317. PubMed Abstract | Publisher Full Text OpenURL

  53. Nakari-Setala T, Aro N, Ilmen M, Munoz G, Kalkkinen N, Penttila M: Differential expression of the vegetative and spore-bound hydrophobins of Trichoderma reesei - cloning and characterization of the hfb2 gene.

    Eur J Biochem 1997, 248(2):415-423. PubMed Abstract | Publisher Full Text OpenURL

  54. Bailey MJ, Askolin S, Hörhammer N, Tenkanen M, Linder M, Penttilä M, Nakari-Setälä T: Process technological effects of deletion and amplification of hydrophobins I and II in transformants of Trichoderma reesei.

    Appl Microbiol Biotechnol 2002, 58(6):721-727. PubMed Abstract | Publisher Full Text OpenURL

  55. Kim S, Ahn IP, Rho HS, Lee YH: MHP1, a Magnaporthe grisea hydrophobin gene, is required for fungal development and plant colonization.

    Mol Microbiol 2005, 57(5):1224-1237. PubMed Abstract | Publisher Full Text OpenURL

  56. Djonovic S, Vargas WA, Kolomiets MV, Horndeski M, Wiest A, Kenerly CM: A proteinaceous elicitor Sm1 from the beneficial fungus Trichoderma virens is required for induced systemic resistence in maize.

    Plant Physiol 2007, 145(3):875-889. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Vargas WA, Djonovic S, Sukno SA, Kenerley CM: Dimerization controls the activity of fungal elicitors that trigger systemic resistance in plants.

    J Biol Chem 2008, 283(28):19804-19815. PubMed Abstract | Publisher Full Text OpenURL

  58. Sestak S, Farkas V: Metabolic regulation of endoglucanase synthesis in Trichoderma reesei: participation of cyclic AMP and glucose-6-phosphate.

    Can J Microbiol 1993, 39(3):342-347. PubMed Abstract | Publisher Full Text OpenURL

  59. Seiboth B, Hakola S, Mach RL, Suominen PL, Kubicek CP: Role of four major cellulases in triggering of cellulase gene expression by cellulose in Trichoderma reesei.

    J Bacteriol 1997, 179(17):5318-5320. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Abrahao-Neto J, Rossini CH, el Gogary S, Henrique-Silva F, Crivellaro O, el Dorry H: Mitochondrial functions mediate cellulase gene expression in Trichoderma reesei.

    Biochemistry 1995, 34(33):10456-10462. PubMed Abstract | Publisher Full Text OpenURL

  61. Martinez D, Challacombe J, Morgenstern I, Hibbett D, Schmoll M, Kubicek CP, Ferreira P, Ruiz-Duenas FJ, Martinez AT, Kersten P, et al.: Genome, transcriptome, and secretome analysis of wood decay fungus Postia placenta supports unique mechanisms of lignocellulose conversion.

    Proc Natl Acad Sci USA 2009, 106(6):1954-1959. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Langston JA, Shaghasi T, Abbate E, Xu F, Vlasenko E, Sweeney MD: Oxidoreductive cellulose depolymerization by the enzymes cellobiose dehydrogenase and glycoside hydrolase 61.

    Appl Environ Microbiol 2011, 77(19):7007-7015. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JC, Johansen KS, Krogh KB, Jorgensen CI, Tovborg M, Anthonsen A, et al.: Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components.

    Proc Natl Acad Sci USA 2011, 108(37):15079-15084. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Beeson WT, Phillips CM, Cate JH, Marletta MA: Oxidative cleavage of cellulose by fungal copper-dependent polysaccharide monooxygenases.

    J Am Chem Soc 2012, 134(2):890-892. PubMed Abstract | Publisher Full Text OpenURL

  65. Phillips CM, Beeson WT, Cate JH, Marletta MA: Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by Neurospora crassa.

    ACS Chem Biol 2011, 6(12):1399-1406. PubMed Abstract | Publisher Full Text OpenURL

  66. Zhang R, Fan Z, Kasuga T: Expression of cellobiose dehydrogenase from Neurospora crassa in Pichia pastoris and its purification and characterization.

    Protein Expr Purif 2011, 75(1):63-69. PubMed Abstract | Publisher Full Text OpenURL

  67. Paluh JL, Orbach MJ, Legerton TL, Yanofsky C: The cross-pathway control gene of Neurospora crassa, cpc-1, encodes a protein similar to GCN4 of yeast and the DNA binding domain of the oncogene v-jun-encoded protein.

    Proc Natl Acad Sci USA 1988, 85(11):3728-3732. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  68. Schuster A, Schmoll M: Biology and biotechnology of Trichoderma.

    Appl Microbiol Biotechnol 2010, 87(3):787-799. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  69. Turner BC, Perkins DD, Fairfield A: Neurospora from natural populations: a global study.

    Fungal Genet Biol 2001, 32(2):67-92. PubMed Abstract | Publisher Full Text OpenURL

  70. Kubicek CP, Mikus M, Schuster A, Schmoll M, Seiboth B: Metabolic engineering strategies for improvement of cellulase production by Hypocrea jecorina.

    Biotechnol Biofuels 2009, 2(1):19. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  71. Tisch D, Schmoll M: Novel approaches to improve cellulse biosynthesis for biofuel production - adjusting signal transduction pathways in the biotechnological workhorse Trichoderma reesei. In Biofuel production - Recent developments and prospects. Edited by dos Santos Bernardes MA. Intech: Rijeka, Croatia; 2011:199-224. OpenURL

  72. Gyalai-Korpos M, Nagy G, Mareczky Z, Schuster A, Reczey K, Schmoll M: Relevance of the light signaling machinery for cellulase expression in Trichoderma reesei (Hypocrea jecorina).

    BMC Res Notes 2010, 3:330. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  73. Schafmeier T, Diernfellner AC: Light input and processing in the circadian clock of Neurospora.

    FEBS Lett 2011, 585(10):1467-1473. PubMed Abstract | Publisher Full Text OpenURL

  74. Sommer T, Chambers JA, Eberle J, Lauter FR, Russo VE: Fast light-regulated genes of Neurospora crassa.

    Nucleic Acids Res 1989, 17(14):5713-5723. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  75. Aro N, Pakula T, Penttila M: Transcriptional regulation of plant cell wall degradation by filamentous fungi.

    FEMS Microbiol Rev 2005, 29(4):719-739. PubMed Abstract | Publisher Full Text OpenURL

  76. de Vries RP, Visser J: Aspergillus enzymes involved in degradation of plant cell wall polysaccharides.

    Microbiol Mol Biol Rev 2001, 65(4):497-522. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  77. Ruijter GJ, Visser J: Carbon repression in Aspergilli.

    FEMS Microbiol Lett 1997, 151(2):103-114. PubMed Abstract | Publisher Full Text OpenURL

  78. Schmoll M, Esquivel-Naranjo EU, Herrera-Estrella A: Trichoderma in the light of day - physiology and development.

    Fungal Genet Biol 2010, 47(11):909-916. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  79. Cujec TP, Tyler BM: Nutritional and growth control of ribosomal protein mRNA and rRNA in Neurospora crassa.

    Nucleic Acids Res 1996, 24(5):943-950. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  80. Shi YG, Tyler BM: Coordinate expression of ribosomal protein genes in Neurospora crassa and identification of conserved upstream sequences.

    Nucleic Acids Res 1991, 19(23):6511-6517. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  81. Sturani E, Costantini MG, Zippel R, Alberghina FA: Regulation of RNA synthesis in Neurospora crassa. An analysis of a shift-up.

    Exp Cell Res 1976, 99(2):245-252. PubMed Abstract | Publisher Full Text OpenURL

  82. de la Serna I, Ng D, Tyler BM: Carbon regulation of ribosomal genes in Neurospora crassa occurs by a mechanism which does not require Cre-1, the homologue of the Aspergillus carbon catabolite repressor, CreA.

    Fungal Genet Biol 1999, 26(3):253-269. PubMed Abstract | Publisher Full Text OpenURL

  83. Friedl MA, Schmoll M, Kubicek CP, Druzhinina IS: Photostimulation of Hypocrea atroviridis growth occurs due to a cross-talk of carbon metabolism, blue light receptors and response to oxidative stress.

    Microbiology 2008, 154(Pt 4):1229-1241. PubMed Abstract | Publisher Full Text OpenURL

  84. Mauro VP, Edelman GM: The ribosome filter hypothesis.

    Proc Natl Acad Sci USA 2002, 99(19):12031-12036. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  85. Mauro VP, Edelman GM: The ribosome filter redux.

    Cell Cycle 2007, 6(18):2246-2251. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  86. Denault DL, Loros JJ, Dunlap JC: WC-2 mediates WC-1-FRQ interaction within the PAS protein-linked circadian feedback loop of Neurospora.

    EMBO J 2001, 20(1-2):109-117. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  87. Esquivel-Naranjo EU, Herrera-Estrella A: Enhanced responsiveness and sensitivity to blue light by blr-2 overexpression in Trichoderma atroviride.

    Microbiology 2007, 153(Pt 11):3909-3922. PubMed Abstract | Publisher Full Text OpenURL

  88. McCluskey K: The Fungal Genetics Stock Center: from molds to molecules.

    Adv Appl Microbiol 2003, 52:245-262. PubMed Abstract OpenURL

  89. Colot HV, Park G, Turner GE, Ringelberg C, Crew CM, Litvinkova L, Weiss RL, Borkovich KA, Dunlap JC: A high-throughput gene knockout procedure for Neurospora reveals functions for multiple transcription factors.

    Proc Natl Acad Sci USA 2006, 103(27):10352-10357. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  90. Vogel HJ: A convenient growth medium for Neurospora (medium N).

    Microbial Genetics Bulletin 1956, 13:42-43. OpenURL

  91. Metzenberg RL: Bird medium: an alternative to Vogel medium.

    Fungal Genet Newsl 2004, 51:19-20. OpenURL

  92. Froehlich AC, Noh B, Vierstra RD, Loros J, Dunlap JC: Genetic and molecular analysis of phytochromes from the filamentous fungus Neurospora crassa.

    Eukaryot Cell 2005, 4(12):2140-2152. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  93. Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.

    Nucleic Acids Res 2002, 30(9):e36. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  94. Vinciotti V, Khanin R, D' Alimonte D, Liu X, Cattini N, Hotchkiss G, Bucca G, de Jesus O, Rasaiyaah J, Smith CP, et al.: An experimental evaluation of a loop versus a reference design for two-channel microarrays.

    Bioinformatics 2005, 21:492-501. PubMed Abstract | Publisher Full Text OpenURL

  95. Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, Tetko I, Guldener U, Mannhaupt G, Munsterkotter M, et al.: The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes.

    Nucleic Acids Res 2004, 32(18):5539-5545. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL