- Open Access
Disruption of gul-1 decreased the culture viscosity and improved protein secretion in the filamentous fungus Neurospora crassa
- Liangcai Lin†1,
- Zhiyong Sun†1,
- Jingen Li1,
- Yong Chen1,
- Qian Liu1,
- Wenliang Sun1 and
- Chaoguang Tian1Email author
© The Author(s) 2018
- Received: 6 March 2018
- Accepted: 12 June 2018
- Published: 16 June 2018
The cellulolytic fungus Neurospora crassa is considered a potential host for enzyme and bioethanol production. However, large scale applications are hindered by its filamentous growth. Although previous investigations have shown that mycelial morphology in submerged culture can be controlled by altering physical factors, there is little knowledge available about the potential for morphology control by genetic modification.
In this study, we screened morphological mutants in the filamentous fungus N. crassa. Of the 90 morphological mutants screened, 14 mutants exhibited considerably higher viscosity compared with that of the wild type strain, and only two mutants showed low-viscosity morphologies in submerged culture. We observed that disruption of gul-1 (NCU01197), which encodes an mRNA binding protein involved in cell wall remodeling, caused pellet formation as the fermentation progressed, and resulted in the most significant decrease in viscosity of culture broth. Moreover, over-expression of gul-1 caused dramatically increased viscosity, suggesting that the gul-1 had an important function in mycelial morphology during submerged cultivation. Transcriptional profiling showed that expression of genes encoding eight GPI-anchored cell wall proteins was lowered in Δgul-1 while expression of genes associated with two non-anchored cell wall proteins was elevated. Meanwhile, the expression levels of two hydrophobin genes were also significantly altered. These results suggested that GUL-1 affected the transcription of cell wall-related genes, thereby influencing cell wall structure and mycelial morphology. Additionally, the deletion of gul-1 caused increased protein secretion, probably due to a defect in cell wall integrity, suggesting this as an alternative strategy of strain improvement for enzyme production. To confirm practical applications, deletion of gul-1 in the hyper-cellulase producing strain (∆ncw-1∆Ncap3m) significantly reduced the viscosity of culture broth.
Using the model filamentous fungus N. crassa, genes that affect mycelial morphology in submerged culture were explored through systematic screening of morphological mutants. Disrupting several candidate genes altered viscosities in submerged culture. This work provides an example for controlling fungal morphology in submerged fermentation by genetic engineering, and will be beneficial for industrial fungal strain improvement.
- Neurospora crassa
- Mycelial morphology
- Protein secretion
Filamentous fungi are widely used for industrial scale production of antibiotics, proteins and many other useful chemicals [1–3]. Unlike bacteria and yeast, filamentous fungi are morphologically complex microorganisms, which often leads to lots of process engineering problems during the fermentation process. When grown in submerged culture, fungi display different morphological forms, such as freely dispersed mycelia, clumps and pellets . In general, the clump form predominates in fungal fermentations, resulting in high viscosity and a fermentation broth that shows several non-Newtonian characteristics. Although these disadvantages may not significant impact preliminary Erlenmeyer flask batch fermentations, when scaled up there are negative impacts on nutrient consumption and oxygen uptake, and a decrease in productivity. The simplest strategy to overcome this problem is to increase agitation speed. However, high impeller speed increases power consumption, and produces high shear stress that often damages fungal mycelia reducing product yield . In contrast, the pellet growth form leads to Newtonian rheological behavior and low viscosities. As expected, lower power inputs are required for achieving sufficient agitation and mass transfer. Thus, the pellet morphology is preferred for many industrial processes, and it would greatly facilitate cost-effective production if one could design a way of controlling mycelial morphology during the fermentation process.
Many studies have reported on the effects of cultivation conditions and physical factors on fungal morphology, including media constituents, inoculum size, pH, temperature, agitation systems and the type of fermenter [5–8]. Since the genetics of morphogenesis is very complicated in filamentous fungi, there is little knowledge available about morphology control through genetic modification. In recent years, chemical and UV mutagenesis strategies have been widely used to create mutants that show low-viscosity morphology in submerged culture . However, random mutagenesis is quite time-consuming, and it is also difficult to determine which gene is responsible for a desired phenotype.
In the past century, the filamentous fungus Neurospora crassa has been used as a model organism for genetic and molecular studies , leading to a range of genetic techniques and tools . Compared with other filamentous fungi, only N. crassa possesses a near full genome deletion strain collection. Recently, N. crassa was deployed as a model to unravel mechanisms of lignocellulase expression and regulation due to its capacity to secrete lots of enzymes involved in lignocellulose utilization . In addition, N. crassa has been reported to be a potential alternative candidate for heterologous protein expression and bioethanol production . Several approaches have been employed to improve cellulase production, including regulation of the cellulase induction pathway, regulation of the unfolded protein response pathway, and by enhancing protein secretion [14–16]. Furthermore, N. crassa can be easily sexually crossed to generate multiple-gene mutants, which facilitates construction of hyper-cellulase producers. However, to our knowledge, N. crassa has not been used in industrial scale processes. One of the most limiting factors is its mycelial morphology in submerged culture. The N. crassa wild type strain exhibits a clump type of morphology, which results in high viscosity and mass transfer limitations. Although it is reported that N. crassa can grow as dispersed pellets following the addition of Triton N-101 or anionic polymer carboxypolymethylene [17, 18], investigations into the molecular controls of morphology could have great significance in industrial application.
Therefore, we screened 90 N. crassa morphological mutants to identify putative genes affecting mycelial morphology during submerged cultivation. Several candidate genes affected the viscosity of media when disrupted in N. crassa. Intriguingly, we found that disruption of gul-1 (NCU01197), the homolog of Saccharomyces cerevisiae ssd1, led to a significant change in viscosity. In S. cerevisiae, SSD1 plays critical roles in maintaining cell wall integrity , and has been found to specifically associate with the mRNAs of proteins that are involved in cell wall remodeling . SSD1 is also negatively regulated by the NDR family kinase Cbk1 through phosphorylation . Deletion of ssd1 alters the composition and cell wall architecture of the yeast cell surface . Similar to the functions described for SSD1, GUL-1 is an important component of the NDR kinase COT-1 pathway, which is involved in the regulation of polarized growth in N. crassa [23, 24]. The NDR family kinase COT-1, which is a homolog of the S. cerevisiae Cbk1, is required for normal hyphal elongation . The temperature-sensitive cot-1 mutant exhibits cessation of tip extension and hyperbranching after being shifted to restrictive temperature. Disruption of gul-1 can partially suppress the growth defects of cot-1. Thus, gul-1 has been initially identified as a dominant modifier of cot-1 . Further investigation demonstrated that the Ste20 kinase POD-6 is involve in polar tip extension and also acts in the NDR kinase COT-1 pathway. The growth defects in pod-6 and cot-1 can be suppressed by the common extragenic suppressor, gul-1, suggesting that gul-1 may function as a downstream mediator in this pathway . Recently, Herold and Yarden provided strong evidence to support this hypothesis. Their results illuminated that the NDR family kinase COT-1 affects hyphal growth by influencing the transcript levels of cell wall remodeling genes, which is mediated by gul-1 . Despite the previous studies described above, the mechanism of how gul-1 affects mycelial morphology in submerged culture remains unclear. In this study, we demonstrated that gul-1 alters cell wall structure by influencing the expression of cell wall protein genes, which leading to changes in mycelial morphology. In addition, we also illustrated that morphology engineering enables enhanced protein secretion. Our findings enhance our understanding of the genetics of morphogenesis in filamentous fungi and provided a novel strategy for morphological engineering of filamentous fungi by genetic modification.
A low-viscosity mutant from a morphological mutant library
List of morphological mutants in Neurospora crassa
Increased vs WT (%)
267.5 ± 3.5
Thyroid hormone receptor interactor 12
324.5 ± 34.6
340.5 ± 12.0
430.0 ± 127.3
418.5 ± 9.2
Glycosylhydrolase family 2–2
466.0 ± 36.8
407.5 ± 88.5
599.0 ± 58.0
Secretory phospholipase A2
610.0 ± 56.6
458.0 ± 116.0
Fatty acid oxygenase
630.0 ± 42.4
495.0 ± 41.4
700.0 ± 28.3
42.0 ± 2.8
Heat shock protein 30
261.0 ± 26.9
324.5 ± 54.4
345.0 ± 83.4
366.5 ± 75.7
380.0 ± 127.3
Serine/threonine protein kinase-4
382.5 ± 95.5
401.0 ± 79.2
406.0 ± 19.8
473.5 ± 79.9
High affinity nickel transporter nic1
498.0 ± 46.7
Two-component system response regulator
615.0 ± 62.2
644.5 ± 21.9
Regulatory particle, non-ATPase-like-2
322.5 ± 64.3
Regulatory particle, ATPase-like-3
370.5 ± 43.1
Ribosomal protein L14
375.0 ± 77.8
rRNA 2′-O-methyltransferase fibrillarin
440.0 ± 67.9
453.5 ± 50.2
Proteasome catalytic alpha-2
490.0 ± 155.6
Pre-mRNA-splicing factor rse-1
502.5 ± 51.6
C2H2 conidiation transcription factor FlbC
543.0 ± 28.3
456.0 ± 33.9
Proteasome catalytic alpha-5
580.5 ± 26.2
Proteasome catalytic alpha-4
607.5 ± 41.7
Proteasome catalytic beta-4
651.5 ± 68.6
720.0 ± 14.1
351.0 ± 38.2
371.5 ± 64.3
448.5 ± 78.5
60S ribosome subunit biogenesis protein NIP7
484.5 ± 17.7
485.0 ± 77.8
Proteasome catalytic beta-6
611.0 ± 55.2
633.0 ± 53.7
267.5 ± 17.7
310.5 ± 0.7
Topogenesis of outer membrane beta barrel protein 37
336.5 ± 19.1
332.0 ± 72.1
357.0 ± 42.4
Autoinducer 2 sensor kinase/phosphatase luxQ
376.0 ± 1.4
Mitogen activated protein kinase-1
385.5 ± 33.2
387.0 ± 83.4
Cullin binding protein CanA
397.0 ± 49.5
Glycosylhydrolase family 2–3
436.5 ± 19.1
636.0 ± 58.0
462.9 ± 18.3
306.3 ± 33.6
331.1 ± 26.7
Ubiquitin conjugating enzyme
259.0 ± 29.8
Regulator of conidiation in Aspergillus-1
336.1 ± 33.8
372.0 ± 31.1
283.8 ± 22.9
262.7 ± 24.5
Nuclear import and export protein Msn5
306.3 ± 33.6
334.9 ± 21.4
269.0 ± 43.9
Anchored cell wall protein-10
277.7 ± 45.7
305.0 ± 21.2
NIMA-interacting protein TinC
304.9 ± 21.0
Transport of copper-2
347.3 ± 18.0
364.6 ± 20.6
Ubiquitin-conjugating enzyme E
291.3 ± 12.4
473.0 ± 38.1
317.4 ± 17.8
156.9 ± 32.6
Extracellular developmental signal biosynthesis protein FluG
331.1 ± 26.7
Ubiquitin-conjugating enzyme E
331.0 ± 15.6
Ubiquitin conjugating enzyme Ubc14
364.6 ± 20.6
Abundant perithecial protein
339.9 ± 28.5
373.4 ± 23.4
329.7 ± 27.9
ABC transporter CDR4
353.3 ± 47.1
Serine/threonine protein kinase-12
323.7 ± 37.2
416.9 ± 9.8
b-ZIP transcription factor IDI4
331.1 ± 26.7
361.0 ± 26.9
Vegetative incompatibility blocked-1
415.7 ± 20.2
C2H2 finger domain-containing protein FlbC
356.0 ± 19.8
GUL-1 is involved in mycelial morphology in N. crassa
In order to determine the extent of the morphological changes conferred by gul-1 disruption, we conducted a microscopic analysis of Δgul-1 grown under permissive temperature. When grown on agar plates at 25 °C, there were significant differences in hyphal elongation and branching frequency between Δgul-1 and WT. The Δgul-1 exhibited slow growth rate and hyperbranching (Additional file 1: Figure S1), which is consistent with previous observations . These results indicated that disruption of gul-1 substantially impact hyphal development in N. crassa.
It is well established that mycelial biomass significantly affects the apparent viscosity of fermentation broth. Thus, the biomass of Δgul-1 and WT were determined as described . The results revealed that the dry weight of Δgul-1 was similar to that of the WT (Additional file 3: Figure S3). Thus, the change in rheological properties was due to the morphological change in mycelia, and not to a reduction in biomass.
In addition, fungal morphology in submerged cultures is affected by the type of carbon source . The Δgul-1 grew as pellets, and exhibited dramatically reduced viscosity (> 50%) compared with that of WT when grown on glucose, xylose and fructose (Additional file 4: Figure S4) indicated that growth form is independent of the type of carbon source. This suggests that gul-1 might be considered as a universal target for morphology control under different cultivation conditions.
To the best of our knowledge, this is the first report indicating that N. crassa can change morphological forms in submerged culture (from clump type to pellet type) through genetic modification.
Expression level of gul-1 affected culture viscosity
Subcellular localization of the GUL-1 protein in N. crassa
Sensitivity to cell wall inhibitors and oxidative agents
Oxidative stress was provided by hydrogen peroxide (H2O2), diamide or methyl-viologen. These oxidants elicit oxidative damage to fungal cells in different ways. H2O2 is a common oxidative agent able to react with many biological molecules including proteins and DNA . Diamide triggers oxidative stress by oxidizing intracellular glutathione . Methyl-viologen is a redox cycling reagent that produces superoxide anions under aerobic condition . All three oxidative agents partially inhibited the elongation of hyphae in both Δgul-1 and WT strains; only diamide was more effective on Δgul-1 than the WT (Fig. 5). The growth of Δgul-1 on MM agar containing 9 μg/mL diamide was reduced by 38.6 ± 0.4% as compared to 24.8 ± 0.5% for the WT. These results indicate that gul-1 has a role in glutathione-mediated antioxidant processes, but not in general oxidative stress protection of N. crassa.
Δgul-1 grown on MM medium supplemented with NaCl exhibited no greater sensitivity to osmotic stress than the wild type strain (Fig. 5), indicating that gul-1 is not required for osmotic stress tolerance.
Deletion of gul-1 improved protein secretion
Comparative analysis of transcriptomes of WT and the Δgul-1 strain
In addition, the expression levels of major cellulase genes were not affected by the gul-1 deletion, except that the extracellular β-glucosidase gene (NCU04952) had significantly higher expression level in Δgul-1 (Fig. 7b and c). These results are consistent with the enzyme activity assay shown in Fig. 6.
Genetic engineering of N. crassa for hyper-cellulase production and lower viscosity
In fungal fermentations, high apparent viscosity and non-Newtonian behavior of broths are critical problems in industrial application. Thus, controlling the morphology of filamentous fungus in submerged fermentations is a major challenge. Many strategies have been applied to solve these problems, principally by optimization of fermentation process parameters [5–8]. However, morphology control by altering chemical and physical factors is not sufficient. Moreover, these optimized parameters are not always suited for achieving maximum productivity. Recently, several attempts have been employed to control mycelial morphology by genetic modification. As previously reported, the motor proteins played critical roles in polarized growth in Aspergillus spp. Disruption of kinesin motor KipA in Aspergillus glaucus caused pellet growth in submerged cultivation and the aspergiolide A production was 82% higher than that of wild type strain . Furthermore, disturbance of cell wall synthesis had significant influences on mycelial morphology . In Aspergillus orzyae, disruption of chsB had no effect on enzyme production, but the viscosity was lower than that of parental strain . In Penicillium chrysogenum, down-regulation of chitin synthase gene resulted in pellet growth in submerged culture, and also caused enhanced penicillin production . Thus, rationally shaping mycelial morphology through genetic modification seems to be an attractive strategy. However, large-scale generation and screening of morphological mutants is difficult to achieve. Fortunately, N. crassa possesses a near full genome deletion strain collection . In this work, we systematically screened the 90 N. crassa morphological mutants to identify the critical genes that could influence mycelial morphology during submerged cultivation. The viscosity of fermentation broth was used to reflect mycelial morphology change. Fourteen mutants showed increased viscosity compared to the WT strain. Among these mutants, we found several genes involved in the proteasome pathway, e.g., NCU01368, NCU05295, NCU09366 and NCU06440. All these four genes encode different subunits of the proteasome, implying that interruption of the protein degradation pathway might influence mycelial morphology. Moreover, NCU00824 encodes a histone deacetylase which possibly contributes to chromatin structure modification. Previous work in the yeast C. albicans has demonstrated that histone deacetylase interacting with transcription factors plays a critical role in morphological transition . The high-viscosity phenotype induced by a defect in NCU00824 implies a close link between chromatin structure modification and morphology control in filamentous fungi.
Of most interest, we found that disruption of gul-1 (NCU01197) in N. crassa led to favorable morphological changes in submerged culture. The Δgul-1 mutant grew as pellets and showed a low viscosity phenotype. Furthermore, scanning electron microscopy observations revealed that the cell wall structure was altered in the gul-1 mutant. Our results with Congo Red and Calcofluor White also illustrated that gul-1 is required for cell wall integrity. Based on these results, we hypothesized that pellet formation is closely associated with changes in cell wall structure. In filamentous fungi, β-(1,3)-glucan and chitin are two major structural components of the fungal cell wall. However, according to our results, the transcript levels of genes encoding chitin and glucan synthases (chs1-7 and gls-1), which are involved in cell wall synthesis, were not changed by the gul-1 deletion (Additional file 6: Table S2). Moreover, the expression levels of chit-1 (NCU02184) and gh18-5 (NCU04554), which are involved in cell wall maintenance and remodeling , were also not altered in Δgul-1 (Additional file 6: Table S2). It suggested that these genes are probably not involved in conferring the gul-1 phenotype. On the other hand, cell wall proteins are the important components of the fungal cell wall, and also play important roles in cell wall biosynthesis and interactions with external environment . Intriguingly, the transcript abundance of several cell wall protein genes was found to be significantly altered in Δgul-1. Among these genes, the GPI-anchored cell wall β-1,3-endoglucanase gene bgt-2 (NCU09175), which is seemingly involved in cell wall remodeling, was down-regulated 2.1-fold in Δgul-1. A recent study also showed the cell wall architecture is altered in the Δbgt-2 mutant . Likewise, the transcript abundances of GPI-anchored glycoside hydrolases NCU05974 and NCU06781, the homologs of S. cerevisiae crh1 and Aspergillus fumigatus gel2, respectively, were also lower in Δgul-1. It has been reported that both S. cerevisiae crh1 and A. fumigatus gel2 are required for normal cell wall structure. Loss of S. cerevisiae crh1 causes hypersensitivity to Congo Red and Calcofluor White, suggesting defects in cell wall structure . In A. fumigatus, the gel2 disruption strain exhibits slower growth and an altered cell wall composition . These implicated that down-regulation of both NCU05974 and NCU06781 may results in an abnormal cell wall structure, leading to influencing mycelial morphology during submerged cultivation. Besides that, other GPI-anchored cell wall proteins, which are of unknown function, are annotated as ACW (GPI-anchored cell wall protein) . They may serve as important structural components of the wall, and play critical roles in maintaining cell wall stability and function . For example, AfuEcm33, a homolog of the N. crassa GPI-anchored cell wall protein-1 (acw-1), has been shown to be involved in maintaining cell wall integrity and virulence in A. fumigatus . It is worth noting that disruption of AfuEcm33 results in increased cell–cell adhesion. The disruption strain forms large ‘star-shaped’ clusters of germinating conidia during submerged cultivation . In Candida albicans, a similar aggregate phenotype has been also observed in the Caecm33 deletion strain . In this study, the GPI-anchored cell wall protein ACW-1 encoded by NCU08936 was found to be down-regulated in Δgul-1, suggesting that ACW-1 might, at least in part, contribute to pellet formation by increased cell–cell adhesion. Taken together, our results indicated that gul-1 affects mycelial morphology during submerged cultivation by influencing the expression of cell wall-related protein. However, a detailed mechanistic description of how the gul-1 influences these genes expression requires further investigation.
According to our results, the deletion of gul-1 led to a hypersecretion phenotype in N. crassa. When grown on 2% (w/v) Avicel, the extracellular protein secreted by the Δgul-1 mutant was about 25% higher than that in the WT, while only β-glucosidase activity was found to be improved. These results were also confirmed by our transcriptomic analysis (Fig. 7b and c). However, the increase in activity was disproportionate to the increase in total protein, suggesting that the effect of gul-1 on protein secretion is not due to increased transcription levels. Thus, we hypothesized that the defects in cell wall integrity benefit protein secretion. According to the ‘bulk flow’ hypothesis, growing hyphal tips are filled with vesicles that contain extracellular proteins, and these proteins get carried to the outer region of the wall for secretion . However, some extracellular proteins might become trapped in the rigid cell wall. These ‘trapped’ enzymes could be released into an external medium through a weakening of cell wall integrity , which may explain the Δgul-1 mutant having a higher protein secretion capacity. Similar results could be observed in the N. crassa slime mutant which has serious defects in cell wall synthesis .
Our previous study demonstrated that the disruption of ncw-1 led to enhanced cellulase production in N. crassa, but unfortunately the viscosity of ∆ncw-1 was increased two fold . A similar phenotype was also observed in the ∆ncw-1∆Ncap3m strain. Having identified the impact of gul-1 disruption on mycelial morphology during submerged cultivation, we applied this approach to the hyper-cellulase producing strain ∆ncw-1∆Ncap3m. The triple deletion strain (∆ncw-1∆Ncap3m∆gul-1) grew in pellet form and exhibited a low viscosity phenotype, indicating that gul-1 is a promising target gene for morphological improvement of hyper-cellulase-producing strains. In addition, phylogenetic analyses demonstrated that GUL-1 is highly conserved in other industrial filamentous fungi, including Aspergillus niger (XP_001392470, identity 66%), Myceliophthora thermophila (XP_003661716, 75% identity) and Trichoderma reesei (XP_006964501, 78% identity) (Additional file 8: Figure S5). Thus, manipulation of gul-1 homologs could be a highly useful strategy for morphology engineering of industrial filamentous fungi.
In this work, we systematically screened morphological mutants of N. crassa and identified several potential genes that affect mycelial morphology in submerged cultivation. Deletion of one of these genes, gul-1, resulted in the pellet form and enhanced protein secretion in N. crassa. Transcriptional profiling revealed that a deletion of gul-1 led to an alteration in transcription of genes encoding cell wall proteins, which may be the reason for changes in mycelial morphology during submerged cultivation. Furthermore, a significant reduction in viscosity of broth culture was obtained when gul-1 was deleted from the hyper cellulolytic strain ∆ncw-1∆Ncap3m. These results indicated that gul-1 might be a potential target to control mycelial morphology in submerged cultivation. This work provides a novel strategy for morphological engineering of filamentous fungi through genetic modification, especially for industrially relevant species.
Strains, growth conditions and microscopy
Neurospora crassa strains were obtained from Fungal Genetics Stock Centre (FGSC), including wild type (WT, FGSC 2489, A), a his-3 mutant strain (FGSC 9716, his-3, a) and the morphological mutants. The his-3; Δgul-1, a strain was obtained by crossing FGSC 9716 with FGSC 11288 (NCU01197, A). Multiple deletion strains were made by performing sequential crosses. All of the constructed strains were verified by PCR.
To obtain conidia, N. crassa was grown on slant tubes containing Vogel’s minimal media (1 × Vogel’s salts, 2% w/v sucrose and 1.5% w/v agar) for 10 days at 28 °C. For flask culture, 10-day-old conidia of N. crassa strains were suspended in sterile water and inoculated using a final concentration of 1 × 105 per mL into a 250-mL Erlenmeyer flask containing 100 mL Avicel medium (1 × Vogel’s salts, 2% w/v Avicel, 0.75% w/v yeast extract and 0.2% v/v Tween 80), then cultured for 7 days at 25 °C, with shaking at 200 rpm under constant light. The observations were carried out on 7 consecutive days at 24 h intervals.
Microscopic observations were carried out using an OLYMPUS BX51 microscope equipped with a QImagingRetiga 2000R camera (QImaging, Surrey, Canada), and analyzed with the Image-Pro Express 6.3 software. The images were also acquired by an Olympus SZX-7 stereomicroscopy with a digital camera attached.
To compare the differences in cell surface structure between the WT and Δgul-1, fungal conidia were incubated in Avicel medium at 25 °C for 5 days. The harvested mycelia were fixed in 0.1 M phosphate buffer (pH 7.2) containing a final concentration of 2.5% (w/v) glutaraldehyde. The samples were rinsed in thrice-distilled PBS, and dehydrated with a series of ethanol solutions (30–100%, v/v). After which the samples were dried by critical point drying (Leica CPD300) and coated with platinum (Hitachi E-1045). Observations were conducted using a scanning electron microscope (Hitachi SU8010) operating at 1 kV.
Plasmid construction and transformation
Genomic DNA from the WT strain was extracted as described previously . To complement Δgul-1, the open reading frame (ORF) of gul-1 was cloned by PCR with primers ORF-F and ORF-R, and the fragment was inserted into the XbaI and PacI sites of pMF272 to form pMF272-Pc-gul-1-gfp. The promoter of gul-1 was amplified with primers PF and PR, and then cloned into the NotI and XbaI sites of pMF272-Pc-gul-1-gfp to form pMF272-Pn-gul-1-gfp. The primers used in this study are listed in Additional file 9: Table S4. All PCR products in this study were cloned into the pGEM-T vector (Promega) and sequenced for confirmation.
Plasmid DNA (10 μg) was transformed into a (his-3; Δgul-1, a) strain as described previously . The constructs were targeted to the his-3 locus by homologous recombination, and correct integration was confirmed by GFP fluorescence and PCR. To recover homokaryotic strains, His+GFP+ transformants were crossed with FGSC 11288 (NCU01197, A). Progeny were selected for histidine prototrophy and GFP fluorescence. The resulting complemented strains were termed Pc-gul-1 and Pn-gul-1.
Enzyme activity assay
For protein assays, 5 mL of culture supernatants were collected at each time point, centrifuged at 12,000 rpm for 10 min to remove mycelia and stored at 4 °C for further analysis. Total extracellular proteins were quantified using a Bio-Rad Protein Assay kit based on absorbance at 595 nm, with bovine serum albumin as the standard. The endoglucanase activity was determined using the azo-CMC assay kit according to manufacturer’s instruction (Megazyme). We used p-nitrophenyl-β-d-glucopyranoside (pNPG) for the measurement of β-glucosidase activity as described previously . Briefly, a 250 μL culture supernatant that was diluted with 50 mM sodium citrate (pH 4.8) was added to 250 μL of 1 mg/mL pNPG, and then incubated immediately for 10 min at 50 °C. The reaction mixture was terminated by adding 500 μL of 1 M Na2CO3 and determined at an absorbance of 420 nm. Inactive enzyme, which was boiled at 100 °C for 10 min, was used as a control. pNP was used to generate a standard curve. One unit (U) of β-glucosidase activity was defined as the amount of pNP released from the substrate per minute using 1 mL enzyme under the standard assay conditions.
Fungal biomass assay and rheological measurement
The dry weight of fungal biomass was indirectly measured as previously described . Briefly, a 5-mL culture broth was centrifuged at 4000 rpm for 10 min, the supernatant was discarded, and the pellet was dried and weighed. The dried residue was resuspended in 3 mL of acetic nitric reagent (80:20, v/v), and bathed with boiling water for 2 h to remove mycelia. The residual Avicel was washed thoroughly, then dried, and reweighed. Mycelial dry weight was defined as the dry weight of the original 5-mL culture minus that of residual Avicel.
Rheological properties of fermentation broth were measured with a Brookfield digital HBDV-III Ultra rheometer according to the manufacturer’s instructions. Viscosity was measured as the torque exerted by the mycelial sample on the rheometer impeller at 50 rpm. Torque was recorded after the reading had stabilized (45–60 s).
Oxidative stress, high osmotic stress and cell wall integrity test
For plate assays, 5 μL aliquots of conidial suspensions (107 mL−1) were applied onto Vogel’s solid medium (MM: 1 × Vogel’s salts, 2% w/v sucrose and 1.5% w/v agar) supplemented with H2O2 (10, 20 mM), diamide (6, 9 μg/mL), methyl-viologen (4.5, 9 μg/mL), NaCl (0.5, 1.0 M), Congo Red (1, 2 mg/mL) or Calcofluor White (200, 400 μg/mL). Plates were incubated at 28 °C for 18 h. All colonies were cross-measured to determine their growth indices under each stress.
RNA extraction, sequencing and data analysis
For gene expression analysis, conidia of ∆gul-1 and WT strains were separately inoculated into Avicel medium (1 × Vogel’s salt, 2% w/v crystalline cellulose, 0.75% w/v yeast extract and 0.2% v/v Tween 80) and batch cultured for 3 days at 25 °C with constant light. Mycelia were harvested by vacuum filtration and frozen immediately in liquid nitrogen. Total RNA from frozen samples was isolated using Trizol reagent (Invitrogen) in accordance with the manufacturer’s protocol, and further treated with DNase I (RNeasy Mini Kit, Qiagen). The RNA integrity was checked by agarose gel electrophoresis. RNA-seq was performed on the Illumina HiSeqTM 2000 platform of Beijing Genomics Institute (BGI, Shenzhen, China). Filtered clean reads were mapped against predicted transcripts from the N. crassa OR74R genome (version 12)  using TopHat (version 2.0.8b) . The alignment results were stored in SAM format files for subsequent analysis. The counts of reads that uniquely mapped to only one gene were calculated for each gene by Htseq-count (http://www-huber.embl.de/users/anders/HTSeq) using SAM files and genome annotation as input. Abundance for each transcript was calculated using the reads per kilobase per million (RPKM) . Genes with altered expression were performed by using R package NOISeq (version 2.6.0)  (Q value ≥ 0.95 used as threshold, which approximately corresponds to a |log2 ratio| ≥ 1). To discover significantly upregulated and downregulated genes between WT and Δgul-1, only the genes with relative high transcriptional abundance (RPKM values above 20 in at least one strain) went into further analysis . Expression differences detected by RNA-seq were also verified by RT-PCR. Profiling data are listed in Additional file 6: Table S2, and the sequence data produced in this study can be accessed (GEO: GSE113321).
Quantitative real-time PCR
Quantitative real-time PCR (qRT-PCR) was performed using a CFX96 real-time PCR detection system (Bio-Rad) with reagents from TOYOBO (One-step qPCR Kit). All assays were performed in triplicate with actin (NCU04173) as the endogenous standard, according to the manufacturer’s instruction. All primers used in this study were listed in Additional file 9: Table S4.
Putative orthologs of gul-1 in selected fungal species were identified as best reciprocal BLAST hits from the National Center for Biotechnology Institute (NCBI) protein database. The phylogenetic analysis was carried out using the Neighbor-Joining in MEGA4 software. Bootstrap values were adjacent to each internal node, representing the percentage of 1000 bootstrap replicates.
Unless otherwise noted, all experiments were performed in triplicate and statistical tests for significance were determined by two-tailed Student’s t test.
LL and CT conceived and designed the project. LL, ZS, YC, QL and WS performed the experiments. LL, JL and CT analyzed the data. LL and CT wrote the manuscript. All authors read and approved the final manuscript.
We thank Dr. Raymond J St. Leger from University of Maryland, College Park for revising the manuscript. We also thank Dr. Timothy C. Cairns for the critical discussions and Huanhuan Zhai for the technical assistance.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this manuscript (and its Additional files).
Consent for publication
Ethics approval and consent to participate
This work was supported by the National Natural Science Foundation of China (31501007 and 31471186).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Nevalainen KM, Te’o VS, Bergquist PL. Heterologous protein expression in filamentous fungi. Trends Biotechnol. 2005;23:468–74.View ArticleGoogle Scholar
- Meyer V. Genetic engineering of filamentous fungi—progress, obstacles and future trends. Biotechnol Adv. 2008;26:177–85.View ArticleGoogle Scholar
- Ward OP. Production of recombinant proteins by filamentous fungi. Biotechnol Adv. 2012;30:1119–39.View ArticleGoogle Scholar
- Li ZJ, Shukla V, Fordyce AP, Pedersen AG, Wenger KS, Marten MR. Fungal morphology and fragmentation behavior in a fed-batch Aspergillus oryzae fermentation at the production scale. Biotechnol Bioeng. 2000;70:300–12.View ArticleGoogle Scholar
- Amanullah A, Christensen LH, Hansen K, Nienow AW, Thomas CR. Dependence of morphology on agitation intensity in fed-batch cultures of Aspergillus oryzae and its implications for recombinant protein production. Biotechnol Bioeng. 2002;77:815–26.View ArticleGoogle Scholar
- Chao Y, Singh D, Yu L, Li Z, Chi Z, Chen S. Secretome characteristics of pelletized Trichoderma reesei and cellulase production. World J Microbiol Biotechnol. 2012;28:2635–41.View ArticleGoogle Scholar
- Nair RB, Lennartsson PR, Taherzadeh MJ. Mycelial pellet formation by edible ascomycete filamentous fungi, Neurospora intermedia. AMB Express. 2016;6:31.View ArticlePubMedPubMed CentralGoogle Scholar
- Abd-Aziz S, Fernandez CC, Salleh MM, Illias RM, Hassan MA. Effect of agitation and aeration rates on chitinase production using Trichoderma virens UKM1 in 2-l stirred tank reactor. Appl Biochem Biotechnol. 2008;150:193–204.View ArticleGoogle Scholar
- Bocking SP, Wiebe MG, Robson GD, Hansen K, Christiansen LH, Trinci AP. Effect of branch frequency in Aspergillus oryzae on protein secretion and culture viscosity. Biotechnol Bioeng. 1999;65:638–48.View ArticleGoogle Scholar
- Roche CM, Loros JJ, McCluskey K, Glass NL. Neurospora crassa: looking back and looking forward at a model microbe. Am J Bot. 2014;101:2022–35.View ArticleGoogle Scholar
- Dunlap JC, Borkovich KA, Henn MR, Turner GE, Sachs MS, Glass NL, McCluskey K, Plamann M, Galagan JE, Birren BW, et al. Enabling a community to dissect an organism: overview of the Neurospora functional genomics project. Adv Genet. 2007;57:49–96.PubMedPubMed CentralGoogle Scholar
- Tian C, Beeson WT, Iavarone AT, Sun J, Marletta MA, Cate JH, Glass NL. Systems analysis of plant cell wall degradation by the model filamentous fungus Neurospora crassa. Proc Natl Acad Sci USA. 2009;106:22157–62.View ArticleGoogle Scholar
- Dogaris I, Mamma D, Kekos D. Biotechnological production of ethanol from renewable resources by Neurospora crassa: an alternative to conventional yeast fermentations? Appl Microbiol Biotechnol. 2013;97:1457–73.View ArticleGoogle Scholar
- Lin L, Chen Y, Li J, Wang S, Sun W, Tian C. Disruption of non-anchored cell wall protein NCW-1 promotes cellulase production by increasing cellobiose uptake in Neurospora crassa. Biotechnol Lett. 2016;39:545–51.View ArticleGoogle Scholar
- Fan F, Ma G, Li J, Liu Q, Benz JP, Tian C, Ma Y. Genome-wide analysis of the endoplasmic reticulum stress response during lignocellulase production in Neurospora crassa. Biotechnol Biofuels. 2015;8:66.View ArticlePubMedPubMed CentralGoogle Scholar
- Pei X, Fan F, Lin L, Chen Y, Sun W, Zhang S, Tian C. Involvement of the adaptor protein 3 complex in lignocellulase secretion in Neurospora crassa revealed by comparative genomic screening. Biotechnol Biofuels. 2015;8:124.View ArticlePubMedPubMed CentralGoogle Scholar
- Rajan JS, Virkar PD. Induced pelletized growth of Neurospora crassa for tyrosinase biosynthesis in airlift fermenters. Biotechnol Bioeng. 1987;29:770–2.View ArticleGoogle Scholar
- Wen SuW, Jun He B. Secreted enzyme production by fungal pellets in a perfusion bioreactor. J Biotechnol. 1997;54:43–52.View ArticleGoogle Scholar
- Kaeberlein M, Guarente L. Saccharomyces cerevisiae MPT5 and SSD1 function in parallel pathways to promote cell wall integrity. Genetics. 2002;160:83–95.PubMedPubMed CentralGoogle Scholar
- Hogan DJ, Riordan DP, Gerber AP, Herschlag D, Brown PO. Diverse RNA-binding proteins interact with functionally related sets of RNAs, suggesting an extensive regulatory system. PLoS Biol. 2008;6:e255.View ArticlePubMedPubMed CentralGoogle Scholar
- Jansen JM, Wanless AG, Seidel CW, Weiss EL. Cbk1 regulation of the RNA-binding protein Ssd1 integrates cell fate with translational control. Curr Biol. 2009;19:2114–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Wheeler RT, Kupiec M, Magnelli P, Abeijon C, Fink GR. A Saccharomyces cerevisiae mutant with increased virulence. Proc Natl Acad Sci USA. 2003;100:2766–70.View ArticleGoogle Scholar
- Seiler S, Vogt N, Ziv C, Gorovits R, Yarden O. The STE20/germinal center kinase POD6 interacts with the NDR kinase COT1 and is involved in polar tip extension in Neurospora crassa. Mol Biol Cell. 2006;17:4080–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Herold I, Yarden O. Regulation of Neurospora crassa cell wall remodeling via the cot-1 pathway is mediated by gul-1. Curr Genet. 2017;63:145–59.View ArticleGoogle Scholar
- Yarden O, Plamann M, Ebbole DJ, Yanofsky C. cot-1, a gene required for hyphal elongation in Neurospora crassa, encodes a protein kinase. EMBO J. 1992;11:2159–66.View ArticleGoogle Scholar
- Terenzi HF, Reissig JL. Modifiers of the cot gene in Neurospora: the gulliver mutants. Genetics. 1967;56:321–9.PubMedPubMed CentralGoogle Scholar
- 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:10352–7.View ArticleGoogle Scholar
- Watters MK, Boersma M, Johnson M, Reyes C, Westrick E, Lindamood E. A screen for Neurospora knockout mutants displaying growth rate dependent branch density. Fungal Biol. 2011;115:296–301.View ArticleGoogle Scholar
- Papagianni M. Fungal morphology and metabolite production in submerged mycelial processes. Biotechnol Adv. 2004;22:189–259.View ArticleGoogle Scholar
- Tanaka S, Yamada K, Yabumoto K, Fujii S, Huser A, Tsuji G, Koga H, Dohi K, Mori M, Shiraishi T, et al. Saccharomyces cerevisiae SSD1 orthologues are essential for host infection by the ascomycete plant pathogens Colletotrichum lagenarium and Magnaporthe grisea. Mol Microbiol. 2007;64:1332–49.View ArticleGoogle Scholar
- Kurischko C, Kuravi VK, Herbert CJ, Luca FC. Nucleocytoplasmic shuttling of Ssd1 defines the destiny of its bound mRNAs. Mol Microbiol. 2011;81:831–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Ba ANN, Pogoutse A, Provart N, Moses AM. NLStradamus: a simple Hidden Markov Model for nuclear localization signal prediction. BMC Bioinform. 2009;10:202.View ArticleGoogle Scholar
- Kopecka M, Gabriel M. The influence of congo red on the cell wall and (1→3)-beta-d-glucan microfibril biogenesis in Saccharomyces cerevisiae. Arch Microbiol. 1992;158:115–26.View ArticleGoogle Scholar
- Ram AF, Klis FM. Identification of fungal cell wall mutants using susceptibility assays based on Calcofluor white and Congo red. Nat Protoc. 2006;1:2253–6.View ArticleGoogle Scholar
- Collinson LP, Dawes IW. Inducibility of the response of yeast cells to peroxide stress. J Gen Microbiol. 1992;138:329–35.View ArticleGoogle Scholar
- Wemmie JA, Steggerda SM, Moye-Rowley WS. The Saccharomyces cerevisiae AP-1 protein discriminates between oxidative stress elicited by the oxidants H2O2 and diamide. J Biol Chem. 1997;272:7908–14.View ArticleGoogle Scholar
- Srijaruskul K, Charoenlap N, Namchaiw P, Chattrakarn S, Giengkam S, Mongkolsuk S, Vattanaviboon P. Regulation by SoxR of mfsA, which encodes a major facilitator protein involved in paraquat resistance in Stenotrophomonas maltophilia. PLoS ONE. 2015;10:e0123699.View ArticlePubMedPubMed CentralGoogle Scholar
- Allgaier S, Taylor RD, Brudnaya Y, Jacobson DJ, Cambareri E, Stuart WD. Vaccine production in Neurospora crassa. Biologicals. 2009;37:128–32.View ArticleGoogle Scholar
- Zhang J, Zhang J. The filamentous fungal pellet and forces driving its formation. Crit Rev Biotechnol. 2016;36:1066–77.View ArticleGoogle Scholar
- 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:768–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Pardo M, Monteoliva L, Vazquez P, Martinez R, Molero G, Nombela C, Gil C. PST1 and ECM33 encode two yeast cell surface GPI proteins important for cell wall integrity. Microbiology. 2004;150:4157–70.View ArticleGoogle Scholar
- Dynesen J, Nielsen J. Surface hydrophobicity of Aspergillus nidulans conidiospores and its role in pellet formation. Biotechnol Prog. 2003;19:1049–52.View ArticleGoogle Scholar
- Znameroski EA, Glass NL. Using a model filamentous fungus to unravel mechanisms of lignocellulose deconstruction. Biotechnol Biofuels. 2013;6:6.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- Phillips CM, Iavarone AT, Marletta MA. Quantitative proteomic approach for cellulose degradation by Neurospora crassa. J Proteome Res. 2011;10:4177–85.View ArticleGoogle Scholar
- Cai M, Zhang Y, Hu W, Shen W, Yu Z, Zhou W, Jiang T, Zhou X, Zhang Y. Genetically shaping morphology of the filamentous fungus Aspergillus glaucus for production of antitumor polyketide aspergiolide A. Microb Cell Fact. 2014;13:73.View ArticlePubMedPubMed CentralGoogle Scholar
- Bowman SM, Free SJ. The structure and synthesis of the fungal cell wall. BioEssays. 2006;28:799–808.View ArticleGoogle Scholar
- Muller C, Hansen K, Szabo P, Nielsen J. Effect of deletion of chitin synthase genes on mycelial morphology and culture viscosity in Aspergillus oryzae. Biotechnol Bioeng. 2003;81:525–34.View ArticleGoogle Scholar
- Liu H, Zheng Z, Wang P, Gong G, Wang L, Zhao G. Morphological changes induced by class III chitin synthase gene silencing could enhance penicillin production of Penicillium chrysogenum. Appl Microbiol Biotechnol. 2013;97:3363–72.View ArticleGoogle Scholar
- Kim J, Lee JE, Lee JS. Histone deacetylase-mediated morphological transition in Candida albicans. J Microbiol. 2015;53:805–11.View ArticleGoogle Scholar
- Kovacs Z, Szarka M, Kovacs S, Boczonadi I, Emri T, Abe K, Pocsi I, Pusztahelyi T. Effect of cell wall integrity stress and RlmA transcription factor on asexual development and autolysis in Aspergillus nidulans. Fungal Genet Biol. 2013;54:1–14.View ArticleGoogle Scholar
- Martinez-Nunez L, Riquelme M. Role of BGT-1 and BGT-2, two predicted GPI-anchored glycoside hydrolases/glycosyltransferases, in cell wall remodeling in Neurospora crassa. Fungal Genet Biol. 2015;85:58–70.View ArticleGoogle Scholar
- Cabib E, Blanco N, Grau C, Rodriguez-Pena JM, Arroyo J. Crh1p and Crh2p are required for the cross-linking of chitin to beta(1-6)glucan in the Saccharomyces cerevisiae cell wall. Mol Microbiol. 2007;63:921–35.View ArticleGoogle Scholar
- Mouyna I, Morelle W, Vai M, Monod M, Lechenne B, Fontaine T, Beauvais A, Sarfati J, Prevost MC, Henry C, Latge JP. Deletion of GEL2 encoding for a beta(1-3)glucanosyltransferase affects morphogenesis and virulence in Aspergillus fumigatus. Mol Microbiol. 2005;56:1675–88.View ArticleGoogle Scholar
- Bowman SM, Piwowar A, Al Dabbous M, Vierula J, Free SJ. Mutational analysis of the glycosylphosphatidylinositol (GPI) anchor pathway demonstrates that GPI-anchored proteins are required for cell wall biogenesis and normal hyphal growth in Neurospora crassa. Eukaryot Cell. 2006;5:587–600.View ArticlePubMedPubMed CentralGoogle Scholar
- Romano J, Nimrod G, Ben-Tal N, Shadkchan Y, Baruch K, Sharon H, Osherov N. Disruption of the Aspergillus fumigatus ECM33 homologue results in rapid conidial germination, antifungal resistance and hypervirulence. Microbiology. 2006;152:1919–28.View ArticleGoogle Scholar
- Martinez-Lopez R, Monteoliva L, Diez-Orejas R, Nombela C, Gil C. The GPI-anchored protein CaEcm33p is required for cell wall integrity, morphogenesis and virulence in Candida albicans. Microbiology. 2004;150:3341–54.View ArticleGoogle Scholar
- Wessels JGH. Tansley Review No. 45 Wall growth, protein excretion and morphogenesis in fungi. New Phytol. 1993;123:397–413.View ArticleGoogle Scholar
- Peberdy JF. Protein secretion in filamentous fungi—trying to understand a highly productive black box. Trends Biotechnol. 1994;12:50–7.View ArticleGoogle Scholar
- Bigger CH, White MR, Braymer H. Ultrastructure and invertase secretion of the slime mutant of Neurospora crassa. J Gen Microbiol. 1972;71:159–66.View ArticleGoogle Scholar
- Lin L, Fang W, Liao X, Wang F, Wei D, St Leger RJ. The MrCYP52 cytochrome P450 monoxygenase gene of Metarhizium robertsii is important for utilizing insect epicuticular hydrocarbons. PLoS ONE. 2011;6:e28984.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang B, Cai P, Sun W, Li J, Tian C, Ma Y. 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422:859–68.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.View ArticleGoogle Scholar
- Tarazona S, Garcia-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: a matter of depth. Genome Res. 2011;21:2213–23.View ArticlePubMedPubMed CentralGoogle Scholar
- McIntyre LM, Lopiano KK, Morse AM, Amin V, Oberg AL, Young LJ, Nuzhdin SV. RNA-seq: technical variability and sampling. BMC Genomics. 2011;12:293.View ArticlePubMedPubMed CentralGoogle Scholar