Scholarly article on topic 'The impact of a single-nucleotide mutation of bgl2 on cellulase induction in a Trichoderma reesei mutant'

The impact of a single-nucleotide mutation of bgl2 on cellulase induction in a Trichoderma reesei mutant Academic research paper on "Biological sciences"

0
0
Share paper
Academic journal
Biotechnol Biofuels
OECD Field of science
Keywords
{""}

Academic research paper on topic "The impact of a single-nucleotide mutation of bgl2 on cellulase induction in a Trichoderma reesei mutant"

Shida et al. Biotechnol Biofuels (2015) 8:230 DOI 10.1186/s13068-015-0420-y

Biotechnology for Biofuels

RESEARCH Open Access

The impact of a single-nucleotide mutation of bgl2 on cellulase induction in a Trichoderma reesei mutant

Yosuke Shida1, Kaori Yamaguchi1, Mikiko Nitta1,2, Ayana Nakamura1, Machiko Takahashi1,

Shun-ichi Kidokoro1, Kazuki Mori3, Kosuke Tashiro3, Satoru Kuhara3, Tomohiko Matsuzawa4, Katsuro Yaoi4,

Yasumitsu Sakamoto5, Nobutada Tanaka6, Yasushi Morikawa1 and Wataru Ogasawara1*

Abstract

Background: The filamentous fungus Trichoderma reesei (anamorph of Hypocrea jecorina) produces increased cellulase expression when grown on cellulose or its derivatives as a sole carbon source. It has been believed that P-glucosidases of T. reesei not only metabolize cellobiose but also contribute in the production of inducers of cellulase gene expression by their transglycosylation activity. The cellulase hyper-producing mutant PC-3-7 developed in Japan has enhanced cellulase production ability when cellobiose is used as the inducer. The comparative genomics analysis of PC-3-7 and its parent revealed a single-nucleotide mutation within the bgl2 gene encoding intracellular P-glucosidase II (BGLII/Cella), giving rise to an amino acid substitution in PC-3-7, which could potentially account for the enhanced cellulase expression when these strains are cultivated on cellulose and cellobiose.

Results: To analyze the effects of the BGLII mutation in cellulase induction, we constructed both a bgl2 revertant and a disruptant. Enzymatic analysis of the transformant lysates showed that the strain expressing mutant BGLII exhibited weakened cellobiose hydrolytic activity, but produced some transglycosylation products, suggesting that the SNP in bgl2 strongly diminished cellobiase activity, but did not result in complete loss of function of BGLII. The analysis of the recombinant BGLII revealed that transglycosylation products might be oligosaccharides, composed probably of glucose linked P-1,4, P-1,3, or a mixture of both. PC-3-7 revertants of bgl2 exhibited reduced expression and induc-ibility of cellulase during growth on cellulose and cellobiose substrates. Furthermore, the effect of this bgl2 mutation was reproduced in the common strain QM9414 in which the transformants showed cellulase production comparable to that of PC-3-7.

Conclusion: We conclude that BGLII plays an important role in cellulase induction in T. reesei and that the bgl2 mutation in PC-3-7 brought about enhanced cellulase expression on cellobiose. The results of the investigation using PC-3-7 suggested that other mutation(s) in PC-3-7 could also contribute to cellulase induction. Further investigation is essential to unravel the mechanism responsible for cellulase induction in T. reesei.

Keywords: Trichoderma reesei, Hypocrea jecorina, P-Glucosidase, Cellobiose, Cellulase induction, Gene regulation

Background

Plant cell walls comprise the most abundant biomasses on the earth and consist predominantly of cellulose. In nature, this biomass is utilized by a variety of cellulolytic

'Correspondence: owataru@vos.nagaokaut.ac.jp

1 Department of Bioengineering, Nagaoka University of Technology,

1603-1 Kamitomioka, Nagaoka, Niigata 940-2188, Japan

Full list of author information is available at the end of the article

organisms, which thereby play a key role in carbon recycling into the ecosystem. Among them, filamentous fungi are considered extremely efficient degraders of plant biomass and express a large amount of cellulases, which comprise three types of enzymes: cellobiohydrolases (EC 3.2.1.91) and endoglucanases (EC 3.2.1.4), which act syn-ergistically to degrade cellulose to cello-oligosaccharides (mainly cellobiose), and p-glucosidases (EC3.2.1.21),

O© 2015 Shida et al. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License Centfcll (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.

which hydrolyze cellobiose into glucose. In addition, the recently discovered class of enzymes known as lytic polysaccharide monooxygenases (LPMOs) has been recognized as an effective auxiliary enzyme for cellulose degradation [1, 2]. Trichoderma reesei (an anamorph of Hypocrea jecorina) is a potent cellulolytic fungus that expresses a complete set of cellulases to degrade cellulose to glucose and is one of the best-studied cellulolytic fungi [3]. In addition, its remarkable ability to produce proteins has made this fungus an important commercial source of cellulases.

Generally, T. reesei expresses cellulase when cellulose is available as the sole carbon source. In addition to cellulose, some soluble disaccharides (cellobiose, a-sophorose, lactose, etc.) also trigger cellulase expression by this organism (reviewed by Bisaria and Mishra, and Kubicek et al. [4, 5]). Among soluble inducers, a-sophorose, a glucosyl-p-1,2-glucoside, shows the highest cellulase-inducing ability and is thought to be the transglycosylation product of p-glucosidase [6]. Under cellulase-inducing conditions, the genes encoding cel-lulases are transcribed coordinately. This phenomenon suggests that a common regulatory machinery controls cellulase expression [7, 8]. To date, several transcription regulators that control cellulase gene expression have been isolated from T. reesei. These include Xyrl, a key activator of cellulase and hemicellulase expression [9, 10], the additional activator ACEII [11] and ACEIII [12], the repressor ACEI [13, 14], and the p-glucosidase activator BglR [15]. In addition, the wide domain carbon catabolite repressor Cre1, which represses cellulase gene expression during growth on glucose, has also been reported [16-18].

Because of the industrial usefulness of T. reesei, cel-lulase hyper-secreting mutants have been isolated by systematic screening strategies involving mutagenesis mediated by UV irradiation or chemical mutagens [19]. In Japan, a T. reesei mutant lineage has been developed with the support of a national project. PC-3-7, a ninth generation mutant of QM6a, exhibits enhanced protein production when grown on a broad carbon source and exhibits particularly elevated cellulase induction by L-sorbose [20, 21]. Recently, it was demonstrated that PC-3-7 induces higher cellulase production on cellobi-ose, which had been considered a poor cellulase inducer for T. reesei [15].

Although factors such as the protein production ability, response to different carbon sources, and the mechanism of cellulase regulation have been analyzed in detail in PC-3-7 [8, 15, 20-30], the genetic and genomic mutation(s) underlying these phenotypes have not been determined as yet. To answer these questions, we carried out a comparative genomic analysis of the Japanese

T. reesei mutant lineage using high-throughput genomic DNA sequencing technologies. We have identified 154 single-nucleotide polymorphisms (SNPs) between QM6a and PC-3-7 [30]. One SNP among them was found in crel, encoding the carbon catabolite repressor, Cre1. Comparison of the PC-3-7 genome sequence with that of its direct parent strain, KDG-12, identified 19 genes with SNPs in the promoter, coding (leading to amino acid substitutions), or terminator regions. One of these genes encodes BglR, which is a novel fungal type Zn(II)2-Cys6 transcription regulator [15]. The effect of these SNPs in crel and bglr on cellulase production was analyzed by a combination of gene complementation and disruption techniques [30]. Electrophoresis mobility shift assay revealed that the mutation of Cre1 affected its DNA binding property. Furthermore, from the comparison between PC-3-7 and Cre1 disruptant of PC-3-7 it was assumed that Cre1 mutation gave rise to partial release from glucose repression of cel-lulase gene expression. The BglR mutation brought about a partial loss of its function, which correlated with the reduced expression of p-glucosidase genes. Consequently, weakened BGL activity leads to decreased production of glucose via cellobiose hydrolysis, resulting in an indirect release from carbon catabolite repression of cellulase expression by glucose. However, it is difficult to explain why PC-3-7 expressed such an extremely high amount of protein in response to only these two mutations.

One gene among the nine genes containing a SNP leading to an amino acid substitution in PC-3-7 is the bgl2 (cella) gene encoding the intracellular p-glucosidase, BGLII (Cel1a) [15]. BGLII belongs to the glycoside hydrolase (GH) family 1 (Carbohydrate Active Enzyme database; [31]). The bgl2 gene has been cloned and the enzymatic property of recombinant BGLII has been characterized and it has been demonstrated that BGLII was localized in T. reesei cells by an over-expression experiment [32]. T. reesei also produces Cel3B, another intracellular GH1 p-glucosidase. Most recently, it was reported that BGLII and Cel3B play an important role in cellulase induction on lactose [33]. However, the physiological role of BGLII remains unclear with regards to the cellulase induction on cellulose. In T. reesei genome, there are nine genes encoding GH3 p-glucosidase [34]. It has been suggested that extracellular BGLs belonging to the GH family 3 not only hydrolyze cellobiose to glucose, but also coordinate cellulase induction to convert cellobiose into cellulase inducers such as a-sophorose by transglycosylation [35, 36]. As mentioned above, PC-3-7 has enhanced cellulase-producing ability on cellobiose and there is a SNP in bgl2. These facts led us to speculate that BGLII plays an important role in cellulase induction during growth on cellobiose. In the present study, we investigated the effects of bgl2 complementation and

disruption in the mutant strain PC-3-7 and the standard strain QM9414 to analyze the impact of the mutation in bgl2 on cellulase induction.

Results

Single-nucleotide point mutation in bgl2 and the resulting amino acid substitution

Comparative genomic analysis of the genes encoding BGLII (bgl2) in KDG-12 and PC-3-7 identified a single-nucleotide difference at position 1298, where guanine was substituted by thymine. This mutation was confirmed by the sequencing of a DNA fragment comprising bgl2 amplified by PCR using PC-3-7 genomic DNA as a template. This mutation resulted in the amino acid substitution of V409F. Figure 1 shows a multiple sequence alignment of BGLII and p-glucosidase genes from the GH family 1 of Trichoderma species. The amino acid corresponding to V409 of T. reesei BGLII is completely conserved and the amino acid sequences flanking it are also highly conserved. Based on this, we surmised that the mutation of V409F could have resulted in a change in the activity of BGLII.

Complementation and disruption of bgl2 in PC-3-7

In order to investigate the effect of mutation of BGLII on cellulase induction and the activity of BGLII in T. reesei, we constructed two types of transformants. One was a complemented strain in which mutated bgl2 was replaced by parental bgl2 from the QM9414 strain and the second

was a bgl2 disruptant. To transform PC-3-7, a bgl2 recovery cassette and a bgl2 deletion cassette were constructed and introduced into PC-3-7 with a pyr4 and tku70 negative background. The resulting bgl2 revertant and bgl2 disruptant were designated as PC-Wbgl2 and PC-Abgl2, respectively.

The strains PC-3-7, PC-Wbgl2, and PC-Abgl2 were cultivated on medium containing Avicel for 3 days and the cellobiase activity was measured in cell-free extracts prepared from each strain. PC-Wbgl2 showed cellobi-ase activity that was four times that of PC-3-7, whereas PC-Abgl2 exhibited lower activity than PC-3-7 (Fig. 2a). These results indicate that BGLIIV409F did not lead to a complete loss of p-glucosidase activity. To determine whether PC-Wbgl2 retained transglycosylation activity, cell-free extracts prepared from each strain were incubated under high cellobiose concentrations. Consequently, formation of transglycosylation products was observed in the presence of the PC-Wbgl2 cellfree extract. PC-3-7 extracts showed reduced activity, but some transglycosylation products were observed, whereas no such compounds were detected in the presence of PC-Abgl2 extracts (Fig. 2b).

Analysis of transglycosylation products by BGLII

The results from PC-3-7 and its transformants revealed that the formation of transglycosylation products from cellobiose was likely due to BGLII or BGLIIV409F.

T. reesei (BAA74959) 387- d jjvk yneyiramvt ril veld

T. viride (AAQ21348) 383- d jJvK yneyiramvt U% veld

T. longibrachiatum (ACZ34300) 386- d jflvk ydeyiramvt ril veld

T. atroviride (EHK41167) 386- e jnvk yndyiramvt ril vtld

T. virens (EHK16702) 386- e |Svk yseyiramft U% atld

T. harzianum (AFD01233) 380- d jsvn yneyiramft 111 mld

T. harzianum (ABN50090) 386- d iSvN yneyiramft ul MID

T. reesei (AAP57758) 403- P |£irq fdshldsisk itqd

yfaw yfaw yfaw yfaw yfaw yfaw

adhy-428 adsm-428 ads3-428

¿1 ™ ¡¡ADgj-428

3SSD®-445

Fig. 1 Alignment of protein sequences belonging to glycoside hydrolase family 1 from Trichoderma. Alignment was determined by the MUSCLE multiple sequence alignment program [60]. Conserved amino acids are indicated by a black background.The T. reesei BGLII (V409F) point mutation is indicated by the box at the top. The secondary structures of BGLII are also indicated at the top. The glycone and aglycone binding sites are indicated by open circles and filled diamonds, respectively [37]. The protein sequences obtained from the NCBI database are: T. reesei BGLII (Accession number: BAA74959), T. viride beta-glucosidase 2 (Accession number: AAQ21384), T. longibrachiatum beta-glucosidase II (Accession number: ACZ34300), T. atroviride glycoside hydrolase family 1 protein (Accession number: EHK41167), T. virens glycoside hydrolase family 1 protein (Accession number: EHK16702), T. harzianum 1,4-beta-glucosidase (Accession number: AFD01233), T. harzianum beta-1,4-glucosidase (Accession number: ABN50090), and T. reesei Cel1b (Accession number AAP57758)

0.8 ■

0.6 ■

a> (0 я

ï 0.4 ■ ®

0.2 ■

<p t / О

га с

га "¡5

250000

200000 150000 Л 100000 50000

250000 200000 150000 100000 50000 0

35000 30000 25000 20000 15000 10000 5000 0

PC-3-7

PC-Wbgl2

30 40 50 retention time (min)

Fig. 2 Effect of bgl2 complementation and disruption in T. reesei PC-3-7. a Specific activity of the intracellular cellobiase from transformants PC-3-7, PC-Wbgl2, and PC-Abgl2. Cellobiase activity represents the mean of triplicate experiments. Error bars indicate standard deviation. b HPLC analysis of transglycosylation products using cell-free extracts from transformants. In the standard chromatogram, G3, G2, and G1 represent cellotriose, cellobiose, and glucose peaks, respectively. Sp represents a-sophorose. Putative transglycosylation products are indicated by arrows

However, there is also the possibility of contamination from other enzymes because we used a cell-free extract for the analysis of enzyme activity. To rule out this possibility and to measure the transglycosylation activity of BGLII exclusively, wild-type BGLII was expressed in E. coli. cDNA of bgl2 was cloned into the expression vector of E. coli, the His-tagged recombinant wild-type BGLII, rBGLIIwt, was expressed by the induction strategy, and rBGLIIwt was purified by affinity chromatography. For the HPLC analysis of transglycosylation products of rBGLIIwt, not only glucose, cellobiose, cellotriose, and

a-sophorose, but also cellotetraose, laminaribiose, and gentiobiose were used as the standard substance. HPLC data (Fig. 3a) showed that the pattern of transglyco-sylation products was almost the same as that from the cell-free extract of PC-3-7. In addition, the peak around 50 min of retention time was consistent with the peak of laminaribiose. To investigate the degree of polymerization of transglycosylation products, the sample was applied to a ligand exchange and size exclusion column. The analysis revealed that some kinds of oligosaccharides were present in the reactant and the larger molecules

retention time (min)

Fig. 3 Chromatogram of transglycosylation products by rBGLIIwt. HPLC analysis of transglycosylation products using purified recombinant wildtype BGLII in E. coli, a The chromatogram generated by the Prominence HPLC system. In the standard chromatogram, G4, G3, G2, and G1 represent cellotetraose, cellotriose, cellobiose, and glucose peaks, respectively. Other ß-disaccharides, a-sophorose, laminaribiose, and gentiobiose are represented by Sp, Lm, and Ge, respectively. b The chromatogram of size exclusion chromatography. Cello-oligosaccharides were used as the standard substance. Putative transglycosylation products are indicated by arrows. Asterisk refers to the peak of a buffer component in the reactant

seemed to be at least tetrasaccharides (Fig. 3b). In order to analyze the transglycosylation products further, the samples were subjected to TLC together with cello-oligo-saccharides, laminari-oligosaccharides, and a-sophorose as the standard material. Figure 4 shows the results of TLC in which cellobiose as the substrate was reduced and glucose as the hydrolysis product was increased according to the increase in the amount of protein in the reaction. In addition, spots having the same migration patterns as that of laminaribiose, laminaritetraose, and cellotetraose were observed. However, few spots that were not consistent with those of the standard substances were also noted.

Cellulase production by PC-3-7, PC-Wbgl2, and PC-Abgl2

From the enzymatic analysis, it was evident that the SNP of bgl2 reduced the cellobiase activity of BGLIIV409F, which appeared to be responsible for three quarters of intracellular cellobiose degrading activity. PC-3-7 shows enhanced cellulase production on cellobiose as the sole carbon source. When T. reesei is cultivated on cellulose, it has been inferred that cellobiose is supplied continuously as cellulose is degraded by cellulase. Therefore, we attempted to determine whether mutation of bgl2 affects cellulase production on cellulose or cellobiose. PC-3-7, PC-Wbgl2, and PC-Abgl2 were grown on media containing Avicel or cellobiose and a time-course of CMCase activities was analyzed. In the case of Avicel culture, the

L2 -G2 -

L4 -G4

G L2 L3 L4 L5 L6 L7

G4 G3 G2 G Sp

rBGLII G6 G5

(jjg in reaction)

Fig. 4 TLC Chromatogram of transglycosylation products by rBGLII. Serial amounts (5-0.075 |jg by two-fold dilution) of rBGLII in 50 |jL of 20 % cellobiose were subjected to TLC. The reaction mixtures were developed by silica gel plates in n-BuOH-AcOH-H2O (3:1:1, by volume) and stained with 0.5 % thymol in EtOH/H2SO4. G and Sp represent Glucose and a-sophorose, respectively. L2-7 and G2-6 represent laminari-oligosaccharides and cello-oligosaccharides. Arrows represent peaks consistent with the standard substance. Open arrow heads represent unknown transglycosylation products

rate of cellulase production by PC-3-7 was faster than that observed in the other two strains at the early stage of cultivation (Fig. 5a). Nevertheless, the other two strains caught up and a comparable final level of production was achieved by all three strains. This comparable cellulase activity could also be observed in the protein patterns seen on SDS-PAGE (Additional file 1: Figure S1). When transformants were grown on cellobiose, however, differences between the three strains began to emerge. PC-3-7 showed a faster production rate at the early cultivation stage similar to what we observed on the Avicel medium, but cellulase production rate of PC-Wbgl2 was significantly slower and the final activity was lower than PC-3-7, whereas PC-Abgl2 showed slower production rate initially, but the final activity was higher than that of PC-3-7 (Fig. 5b). SDS-PAGE analysis also showed that cellulase production of PC-Wbgl2 was lower than that of PC-3-7 or PC-Abgl2 (Additional file 1: Figure S1). These data suggest not only that SNP of bgl2 and the resulting BGLIIV409F mutation affected cellulase production but that another factor(s) responsible for cellobiose signaling may also exist in PC-3-7.

To investigate the transcriptional response of cellulase gene toward several carbon sources without taking cell growth into account, induction experiments were carried out by using the mycelia transfer method. Each transformant was pregrown on glucose and myce-lia were transferred to an induction solution containing a-sophorose, cellobiose, or Avicel as an inducer. For the target gene of expression analysis, cbhl and egll were selected as the major cellulase genes among T. reesei cel-lulase genes that were coordinately expressed with other cellulase genes in T. reesei. xyrl, the gene encoding Xyr1 that is the key activator of T. reesei cellulase genes, was also chosen as the target gene. Total cDNA derived from induced mycelia was used as a template and expression amounts of cbhl, egll, bgl2, and xyrl were determined by quantitative real-time PCR. When a-sophorose, a strong cellulase inducer for T. reesei, was used, the expression amounts of cbhl and xyrl were nearly the same between all three strains, although cbhl expression at time point of 180 min was higher in PC-Wbgl2 compared to the other two strains (Fig. 6a, d). In the case of bgl2 expression, the level of bgl2 of PC-3-7 was higher than that of PC-Wbgl2 at all time points (Fig. 6c). This suggests that bgl2 expression was partly different from cbhl and egll (Fig. 6b). As compared with a-sophorose induction, however, significant difference in gene expression pattern was observed when cellobiose or Avicel was used as an inducer. Compared to PC-3-7, the transcription of cbhl, egll, bgl2, and xyrl of PC-Wbgl2 was higher than that of PC-3-7 at the early stages of induction (after 30-60 min) on cellobiose (Fig. 7), but the pattern was completely

a Avicel

0.0 8 0

60 50 40 30 20 10 0

b cellobiose

1 2 Culture period (day)

1 2 Culture period (day)

-О-PC-3-7 -OPC-Wbg!2 -A-PC-Abg^

1 2 Culture period (day)

Fig. 5 Cellulase productivity of PC-3-7, PC-Wbgl2, PC-Abgl2. Cellulase activity of transformants grown on Avicel (a) or cellobiose (b), measured as CMCase, FPase, and pNPCase activity. As displayed in the panel, the open circle represents PC-3-7. The open square and the open triangle represent PC-Wbgl2 and PC-Abgl2, respectively. Values represent the mean of triplicate experiments. Error bars indicate standard deviation. X-axis shows the enzyme activity per 1 mL culture medium

reversed after 180 min and PC-3-7 showed significantly higher expression than PC-Wbgl2. Similar expression pattern was observed on Avicel as well (Additional file 1: Figure S2). Contrary to cellulase production, a significant delay in transcription of all the other genes was observed in PC-Abgl2 as compared with the other two strains in cellobiose induction (Fig. 7). Furthermore, transcription was not detected from mycelia induced on Avicel (Additional file 1: Figure S2).

Effects of the bgl2 mutation on T. reesei QM9414 strain

PC-3-7 was derived from QM9414 via five stages of mutation and screening and there are at least 154 SNPs within the promoter, coding region, and terminator regions of the genome when compared to wild-type

QM6a. These SNPs and the mutations in crel and bglr might affect cellulase production and their responses to growth on distinct carbon sources in PC-3-7, making it difficult to evaluate the effect of BGLIIV409F in this strain. Therefore, we chose the parental strain QM9414 with a reduced mutation background compared to that of mutants such as PC-3-7 to analyze the effects of BGLI-IV409F. The mutated bgl2 expression cassette and the bgl2 disruption cassette were introduced into QM9414. The resulting BGLIIV409F expressing strain and bgl2 disrupted strain were called QM-Mbgl2 and QM-Abgl2, respectively. The intracellular cellobiase activities of QM9414, QM-Mbgl2, and QMAbgl2 following transglycosylation activities are shown in Additional file 1: Figure S3. As was the case with the transformants of PC-3-7, bgl2 mutation

<л ф

200 180 160 140 120 100 80 60 40 20 0

a cbhl

b egll

С bgl2

d xyrl

Induction period (min)

¡3 PC-3-7 □ PC-Wbgl2 S PC-Abgl2

Fig. 6 Gene expression profile of PC-3-7 transformants in response to a-sophorose induction. Genes analyzed were cbh1 (a), egl1 (b), bgl2 (c), and xyr1 (d). Stippled gray bar represents PC-3-7, solid gray bar represents PC-Wbgl2, and the shaded gray bar indicates PC-Abgl2. Values represent the relative expression of each gene normalized to act1 as an internal control. Values represent the means of triplicate experiments. Error bars indicate standard deviations

and disruption resulted in a significant decrease in the intracellular cellobiase activity (Additional file 1: Figure S3A). In addition, the activity of QM-Mbgl2 was also slightly higher than that of QM-Abgl2. Therefore, it was clear that intracellular cellobiase activity was affected by bgl2 mutation or disruption. In terms of transglycosylation activity, despite its low level, QM-Mbgl2 produced a transglycosylation product (Additional file 1: Figure S3B).

To investigate cellulase production in the transformants, they were cultivated on Avicel or cellobiose as the sole carbon source. When Avicel was used as the carbon source (Fig. 8a), the rate of cellulase production by QM-Abgl2 was initially slower than that observed by QM9414, but toward the end of the culture period, the rate picked up and there was virtually no difference between the two strains with regards to cellulase production. Surprisingly, cellulase production by QM-Mbgl2

was significantly greater than QM9414 and was comparable to that observed by PC-3-7. This observation was supported by the data obtained from SDS-PAGE analysis of the culture filtrate (Additional file 1: Figure S1). Generally, cellobiose is regarded as a poor inducer of cellulase production in QM9414, and we indeed observed significantly less cellulase activity when QM9414 was cultured on cellobiose. However, when cellobiose was used as the poor cellulase inducer for QM9414, QM-Mbgl2 showed significantly more robust protein production than did QM9414 (Fig. 8b). QM-Abgl2 exhibited significantly lower cellulase activity compared to the other two strains. In the case of FPase from QM9414 transformant, activity of each strain was under the detection limit. The cell biomass of these strains was comparable when cultivated on cellobiose, indicating the absence of any influence of the bgl2 mutation or disruption (data not shown).

<u 0£

a cbhl

E 5 (0

с 0 О

(fl tfl <U

.=; 0.8

0.6 0.4 0.2 0.0

С bgl2

5.0 г b egl1 4.0 3.0 2.0 1.0 0.0

2.5 2.0 1.5 1.0 0.5 0.0

30 d xyrl

60 180 30

Induction period (min)

f7:] PC-3-7 □ PC-Wbgl2 ^ PC-Abg!2

Fig. 7 Gene expression profile of PC-3-7 transformants following cellobiose induction. Genes analyzed were cbh1 (a), egl1 (b), bgl2 (c),

and xyr1 (d). Stippled gray bar represents PC-3-7, solid gray bar represents PC-Wbgl2, and the shaded gray barindicates PC-Abgl2. Values represent

the means of triplicate experiments. Error bars indicate standard deviations

Cellulase gene expression following culture on Avicel

To analyze cellulase gene expression in the QM-series during Avicel cultivation, quantitative real-time PCR was performed using total cDNA obtained from myce-lia grown on Avicel culture as a template (Fig. 9a-d). In the QM-series, expression of all genes was greater in QM-Mbgl2 than in the other strains at all the time points examined. In the QM-Abgl2 strain, the expression of cbh1 and egl1 was lower at early time points than that of QM9414, but the final level of expression was comparable to what was observed in QM9414 (Fig. 9a, b). The pattern of expression of xyr1 in each strain derived from QM9414 (Fig. 9d) was nearly identical to the pattern of cellulase gene expression. In addition, the pattern of gene expression appeared to be similar to the cellulase production in each strain (Fig. 8a). In the case of PC-3-7, the final transcription level of cellulase gene was higher than that of PC-Wbgl2, although not at all time points.

Interestingly, the expression of cbh1 and egl1 was the greatest in the PC-Abgl2 strain among the PC-series. In contrast, the expression of cbh1 and egl1 in QM-Abgl2 did not exceed that in QM-9414.

Discussion

In this study, we evaluated the effect of a single-nucleo-tide mutation in the bgl2 gene of the cellulase hyper-producing mutant PC-3-7 developed in Japan on cellulase production. By comparing amino acid sequences between GH1 proteins from Trichoderma species, we noted that V409 was strictly conserved. V409 is located within a loop between the 7th a-helix and the 8th p-sheet. Hence, V409 is close to the glycone or aglycone recognition site [37], and substitution of V409 to phenylalanine, which is a large amino acid, could lead to the alteration of substrate recognition and enzyme activity. Through this study, we concluded that the hydrolysis activity of

a Avicel

^ 40 _J

3 30 >>

.=; 20

35 30 25 20 15 10 5 0

b cellobiose

y it 0.4

о 0.3

CMCase

Culture period (day)

1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0

pNPCase

QM-9414

QM-Mbgl2

QM-Abgl2

Culture period (day)

Fig. 8 Cellulase productivity of QM9414, QM-Mbgl2, QM-Abgl2. Cellulase activity of transformants grown on Avicel (a) or cellobiose (b), measured as CMCase FPase, and pNPCase activity. Activities of QM9414 (filledsquares), QM-Mbgl2 (filledcircles), and QM-Abgl2 (filled triangles) are shown. Values represent the mean of triplicate experiments. Error bars indicate standard deviation. X-axis shows the enzyme activity per 1 mL culture medium

BGLIIV409F was markedly low and that transglycosylation activity, although weak, remained. From a structural point of view, V409 is localized in close proximity to the surface of the BGLII protein. However, structural analysis by a calculation method of solvent-accessible surface [38] showed that the solvent accessibility of V409 was only 1 % (Fig. 10a). If the loop and the following 6th a-helix between F323 and Y360 comprising a part of the catalytic pocket of BGLII were removed, the solvent accessibility of V409 will increase to 30 % (Fig. 10b). This means that V409 interacts with the 6th a-helix and the substitution V409F is expected to affect the static and/or dynamic structure of the loop constituting a part of the substrate pocket indirectly through the 6th a-helix. Therefore, V409 is a very interesting residue from the standpoint of the protein design with a physical perturbation method [39]. We observed reduced cellulase production by reverting the bgl2 mutation of PC-3-7 and a significantly enhanced cellulase production when the mutation was introduced into the bgl2 gene of strain QM9414. These

facts suggest that the single-point mutation in the gene encoding BGLIIV409F affected cellulase expression.

In T. reesei, expression of cellulase genes is induced by cellulose and it has been proposed that this phenomenon is caused by the presence of trace soluble sugars derived from cellulose by constitutively expressed cellulase activity. This theory is supported by two findings: (1) trace amounts of mRNA corresponding to two major cellulase genes (cbh1 and cbh2) are detected in unin-duced cells and (2) the addition of anti-cellulase antibody hindered cellulase formation during growth on cellulose [40, 41]. Cellobiose is the primary end product generated from cellulose degradation by cellulases. It has been shown that cellulase production is induced by cel-lobiose in fungi, including T. reesei [42]. In some cases, cellulase formation is not observed during growth on cel-lobiose due to the hydrolysis of cellobiose into glucose by extracellular ß-glucosidase. Glucose inhibits cellobiose uptake [43] and could give rise to carbon catabolite repression of cellulase expression. Some reports describe

le 100

ГО 0

(Л «

S 4.0 JS

* 3.0 2.0 1.0 0.0

c bgl2

45 40 35 30 25 20 15 10 5 0

b egll

d xyrl

fflî^

Culture period (hr)

□ PC-3-7 □ PC-WTbgl2 S PC-Abgl2

□ QM9414 □ QM-Mbgl2 S QM-Abgl2

Fig. 9 Gene expression in PC-3-7 and QM9414 transformants grown on Avicel. Genes analyzed were cbh1 (a), egl1 (b), bgl2 (c), and xyr1 (d). PC-3-7 series are expressed in gray and the QM-series are represented in white. The relative expression of each gene is normalized that of act1 as an internal control. Values represent the mean of triplicate experiments. Error bars indicate standard deviations

that the addition of a ß-glucosidase inhibitor to the culture, knock-out of the major extracellular ß-glucosidase, or lowering the pH of the culture medium to reduce ß-glucosidase activity resulted in an enhanced cellulase production [36, 44, 45]. Since the Km of ß-glucosidase toward cellobiose is 20-40-fold higher than that of ß-glucoside permease [43], cellobiose is released from cellulose in concentrations that favor uptake rather than hydrolysis. Therefore, it is possible that cellobiose is the in vivo inducer in T. reesei.

The most potent inducer of cellulase gene expression is a-sophorose. It has been proposed that a-sophorose is converted from cellobiose by transglycosylation [6]. Disruption of the gene encoding the major extracellular ß-glucosidase, BGLI, led to a delay in the production of cellulase on cellulose but not on a-sophorose.

Previously, it has been shown that a BGLI overproducing strain showed higher cellulase production than the parental strain under non-saturating concentrations of a-sophorose. However, nojirimycin, a ß-glucosidase inhibitor, inhibited cellulase production on a-sophorose [35]. Therefore, a-sophorose might not be the final inducer for cellulase induction and it is possible that some other ß-glucosidase(s) are involved in true inducer formation. Recently, a comprehensive analysis of cellu-lase gene expression on sophorose and cellulose has been reported [46]. In this report, the authors proposed that both cellobiose and sophorose act as co-inducers of cellulase formation. Furthermore, additional compound with low molecular weight sugar, L-sorbose, promotes cellulase gene expression [8]. This suggests that several signal transduction pathways responsible for each of these

Fig. 10 The active site and the mutation point in the three-dimensional structure of Trichoderma reesei BGLII (PDBID: 3AHY). The catalytic acid/base, Glu165, the catalytic nucleophile, Glu367, and the mutation residue in this study, Val409 are shown as CPK model. All other residues are displayed as ribbon model in a and residues from Phe323 to Tyr360 are omitted in b

inducers might control cellulase expression in parallel. Although bgl2 expression is responsible for a-sophorose [32], it is partially released from the control of Xyrl, the xylanase and cellulase gene activator [26]. This means bgl2 expression is controlled by some other factor on a-sophorose. Therefore, it is possible that BGLII is located upstream of the cellobiose signaling cascade in T. reesei. Recently, however, a report focusing on the physiological role of BGLII described that BGLII participated in the rapid induction of cellulase genes by cellulose and cellobiose [47]. In our study, a significant difference was observed in cellulase gene expression between PC-3-7 and PC-Wbgl2 strains following cellobiose induction, but not following a-sophorose induction. Cellobiase activity in cell-free extracts of the PC-3-7 strain was greatly reduced, but the transglycosylation activity still remained. An analysis of BGLII expressed in a heterologous host showed previously that it could produce a-sophorose from cellobiose or glucose [32]. In this study, we analyzed the transglycosylation products of BGLII using a recombinant enzyme. From the data of TLC (Fig. 4), the major products appeared to be Laminaribiose, laminaritetraose, and cellotetraose and not a-sophorose. In addition to these, unknown substances composed probably of three to five glucose units were also detected. These facts imply that BGLII can produce not only disaccharides but also oligosaccharides composed probably of glucose linked (-1,4, (3-1,3, or a mixture of both. Further investigation is necessary to clarify the identity and characteristics of these larger molecular substances.

The cellobiose that is taken up accumulates in the fungal cell and a fraction of it is converted into a-sophorose or other substances that act as cellulase inducers. Furthermore, since glucose production from cellobiose was lowered, carbon catabolite repression by glucose might be hindered. Partial release of cellulase gene expression from glucose repression by mutation of the repressor Crel [30] might help to enhance cellobiose induction. This is the first report to show that the disruption of bgl2 in T. reesei mutant gives rise to elevated cellulase production on cellobiose. This phenomenon suggests that the cellobiose response of PC-3-7 is affected not only by the bgl2 mutation but also by another mutated gene(s), because there are at least 45 mutations that caused a change in amino acid sequence in the PC-3-7 genome.

PC-3-7, a derivative of QM9414, has many SNPs in its genome compared with QM6a [30]. Therefore, it was necessary to take mutations in PC-3-7 into account for analyzing the effect of bgl2 mutation. Mutations of CRE1 that are involved in carbon catabolite repression and/or BglR related to (-glucosidase expression especially affect the cellulase gene induction in PC-3-7. In our study, we selected the common strain, QM9414, as a recipient of bgl2 gene disruption or mutation because it was closer to the wild type and many mutants have been developed from this strain around the world [19]. Our results show that cellobiose is a weak inducer of cellulase expression for QM9414. Mutation of bgl2 (QM-Mbgl2) led to two-fold higher cellulase production than QM9414 despite the low level of cellulase activity.

QM-Mbgl2 produced comparable levels of cellulase as that of PC-3-7. It was surprising that the bgl2 mutation would bring higher cellulase productivity to QM9414. As mentioned above, cellobiose could be an in vivo inducer in T. reesei. As a result, a significantly low activity of intracellular ß-glucosidase activity of QM-Mbgl2 might have led to the accumulation of cellobiose in the fungal cell. In addition, the transglycosylation activity that was remaining might have caused the formation of induc-ers such as a-sophorose from accumulated cellobiose by transglycosylation. Previous report showed that bgl2 disruption in QM9414 (QM-Abgl2) results in a delay in cellulase production [47]. In our experiments, a delay in cellulase production was also observed in QM-Abgl2 and PC-Abgl2. However, the final cellulase production level of QM-Abgl2 was comparable to that of QM9414, whereas the final production level of PC-Abgl2 was equal to PC-3-7, which possessed the bgl2 mutation, and higher than bgl2 revertant of PC-3-7 (PC-Wbgl2). This result could be explained by the glucose produced by intracellular ß-glucosidase but not BGLII. PC-3-7 was partially relieved of carbon catabolite repression but this machinery works normally in QM9414. In addition, it is also possible that the recently isolated transcription regulator BglR also aids in ß-glucosidase expression, resulting in the promotion of glucose production in the cell.

In T. reesei, the key regulator of cellulase and xylanase genes is Xyr1. Previously, Derntl et al. reported a correlation between the levels of the xyr1 and cbh1/cbh2 transcripts where elevated levels of xyr1 seem to correlate directly with the up-regulation of cbh1 and cbh2 transcription [48]. In this study, the expression of the cellulase genes, cbh1 and cbh2, seemed to reflect xyr1 expression in cells cultivated on cellulose. However, when Avicel induction was performed on mycelia of PC-Abgl2 previously grown on glucose (Additional file 1: Figure S2), the expression of cellulase genes and xyr1 was undetectable. In the cultivation experiment, conidiospores were inoculated directly into the medium. As T. reesei conidia contain cellulases and ß-glucosidase [49, 50], it has been proposed that conidial-bound enzymes may carry out an initial degradation of cellulose. However, in the early induction stage of the transfer experiment, trace amounts of cellobiose might be released from cellulose by constitutively expressed cellulase (mentioned above). Thus, it is thought that cellobiose is taken up by the cell and converted to an inducer or hydrolyzed to glucose. When bgl2 was deleted in PC-Abgl2, the absence of an inducer might be the cause of the delay in PC-Abgl2 cellulase induction. Furthermore, it is possible that cellobi-ose itself served as the inducer. Recently, a mutation in the Xyr1 regulatory domain was reported in a T. reesei mutant used in industry, which had a glucose-blind phenotype [48]. In this strain, the basal level of expression of hemicellulase

and cellulase genes was elevated and the authors suggested a direct effect of glucose upon Xyrl. A functional analysis of XlnR in Aspergillus niger revealed the presence of a glucose inhibitory domain in the C-terminal region [51]. In addition, it has also been reported that xylose triggers the reversible phosphorylation of XlnR of Aspergillus oryzae [52]. If cellulase inducers (cellobiose, a-sophorose, etc.) directly interact with Xyrl, intracellular accumulation of the inducer caused by the bgl2 mutation might represent an artificially evolved ability for a T. reesei mutant in which cellulase production via Xyrl is elevated.

Conclusion

In this study, we found that the T. reesei mutant PC-3-7 developed in Japan possessed a mutated bgl2 gene. Complementation of the mutated bgl2 gene proved that this mutation was one of the means by which PC-3-7 cells acquired high cellulase expression when cultivated on cellulose and cellobiose. However, another mutation in PC-3-7 may also affect the role of cellobiose in cellulase production (Fig. 11). Furthermore, a single-nucleotide change of bgl2 in the standard strain QM9414 brought about a positive impact upon cellulase production. Our results show not only a path to engineer T. reesei for industrial use, but also that mutants developed by our great predecessors have the keys to solve the regulatory puzzle of T. reesei cellulase induction.

Methods

Fungal strains and culture conditions

T. reesei strains QM9414 (ATCC26921) and PC-3-7 (ATCC66579) used in this study were obtained from Kyowa Hakko Bio Co., Ltd. (Tokyo, Japan). T. reesei strain PC-3-7 is a cellulase hyper-producing mutant with enhanced ability to respond to cellulase induction on L-sorbose [21]. Strains were grown on Difco™ Potato Dextrose Agar (PDA; BD, NJ, USA) plates and conidi-ospores were stored in 0.9 % NaCl solution containing 10 % glycerol until use. For cellulase enzyme production, 1 x 107 conidiospores were inoculated in growth medium [20] containing Avicel cellulose (Funacel; Funakoshi ^ Tokyo, Japan) or cellobiose (Sigma-Aldrich, MO, USA) as the sole carbon source and were grown for appropriate durations at 28 °C on a rotary shaker (220 rpm). To carry out induction experiments using the resting cell of T. reesei, mycelia grown on growth medium containing glucose for 48 h were collected by filtration using Miracloth (Merck Millipore, MA, USA) and washed twice with induction medium [23] containing no carbon source. The washed mycelia were transferred to induction medium containing an appropriate carbon source and incubated for appropriate periods at 28 °C with shaking (220 rpm). After incubation, mycelia were

О : glucose GO: cellobiose ^ : sophorose ^p: laminaribiose ООО--: unknown oligosaccharides

Fig. 11 The schematic representation of the model of enhanced cellulase production responsive to cellobiose in T. reesei PC-3-7. Extracellular cel-lobioses are hydrolyzed to glucoses or transglycosylated by extracellular or cell wall-associated BGLs. Glucose causes carbon catabolite repression (CCR) of cellulase expression. However, Crel controlling CCR has mutated, and therefore PC-3-7 is relieved from glucose repression. Alternatively, cellobioses are incorporated in the cell and may be accumulated in the cell because BGLII, the major intracellular BGL of T. reesei has mutated and therefore possesses significantly reduced hydrolyzing activity. The mutated BGLII still has some residual transglycosylation activity and could generate possible cellulase inducers from cellobiose. In addition, other factor(s) that is mutated in the PC-3-7 gene through strain development from QM9414 may affect the role of cellobiose in cellulase induction. The combination of these factors might be the cause for the high cellulase production of PC-3-7

collected as described above and frozen in liquid N2 for mRNA extraction.

Construction of DNA fragments for T. reesei transformation

To construct the tku70 deletion cassette, three DNA fragments were amplified by PCR. A fragment of about 3.2 kbp containing amdS encoding acetoamidase from Aspergillus nidulans, together with its upstream and downstream region was amplified using pBamdS as a template. The second fragment comprised approximately 1.6 kbp of the tku70 upstream region and its template was the genomic DNA of PC-3-7. This fragment had a 15 bp overlap with pUC118 and the amdS fragment in their respective termini. The third fragment comprised approximately 1.5 kbp of the downstream region of tku70 and introduced a 15 bp overlap with the amdS fragment and pUC118 in their respective termini. These three fragments and pUC118 linearized by £coRI digestion were fused to create the plasmid pUdtku70S harboring the

tku70 deletion cassette by using In-Fusion Cloning Kit (Takara-Bio, Shiga, Japan).

To construct the pyr4 disruption cassette, an approximately 4.7 kbp DNA fragment containing 2.3 kbp of the pyr4 upstream region, the pyr4 coding region, and the pyr4 downstream region was amplified using sequence-specific primers. The amplicon was cloned into the HincII site of pUC118 to create pUpyr4F. Subsequently, an NcoI fragment within the pyr4 coding region was removed. Finally, the plasmid pUCdpyr4 was created by self-ligation of the remaining DNA fragment. To construct another type of pyr4 disruption cassette, an amdS fragment was inserted between AatII sites existing within the coding region of pyr4 of pUCdpyr4, resulting in the plasmid pUCdpyr4S.

The bgl2 disrupting cassette and the wild-type bgl2 expression cassette were constructed as follows: an approximately 5.3 kbp DNA fragment containing the bgl2 coding region, with its 1.6 kbp upstream region and 2.2

kbp downstream region, was amplified with sequence-specific primers using QM9414 genomic DNA as a template. AflII sites were introduced to both termini of the PCR product to excise the DNA fragment. This amplicon was inserted into the HincII site of pUC118 to obtain pUbgl2. The pyr4 marker fragment SalI-excised from pUpyr4F was inserted into the ApaI site within the downstream region of bgl2. The resulting plasmid, pUbgl2W contained the wild-type bgl2 construct. The pyr4 marker was inserted into XhoI sites existing within the bgl2 coding and downstream regions of pUbgl2. The resulting plasmid was pUbgl2D, which contained a disrupted bgl2. In addition to these two plasmids, a plasmid containing a mutated bgl2 expression cassette was constructed by site-directed mutagenesis PCR using pUbgl2W as a template. The resulting plasmid was pUbgl2M.

PCR primers used for plasmid construction are listed in Additional file 1: Table S1. PrimeSTAR® HS DNA polymerase (Takara-Bio, Shiga, Japan) was used for DNA amplification. All PCR products were sequenced to confirm that there were no PCR errors.

Transformation of T. reesei

Transformation of T. reesei PC-3-7 was performed using the protoplast-PEG method as described previously [53] with a modification in which Yatalase (Takara-Bio, Shiga, Japan) was used as a protoplasting enzyme instead of Novozyme 234 (Novozymes, Bagsv^rd, Denmark). tku70 disrupting cassette released from pUdtku70S by EcoRI was introduced into PC-3-7. Transformed protoplasts were plated on the minimal medium containing acetoam-ide as the sole nitrogen source. Candidates of transformants possessing the ability to grow on acetoamide were selected and a single colony was isolated on the same medium. After single colony isolation twice, one transformant was confirmed for tku70 disruption by Southern analysis and was named PC3-7Ak. The pyr4 disruption cassette was excised from pUdpyr4 by HindIII and DraI and introduced into PC-3-7Ak. Candidate transformants were screened on the minimum medium containing 5-flruoroorotic acid and uridine. Finally, selected uridine auxotrophic transformant was confirmed for pyr4 disruption by southern analysis and was named PC3-7AKP. QM9414 was transformed by thepyr4 disrupting cassette excised from pUCdpyr4S by HindIII and XbaI. Candidates of transformants possessing the ability to grow on acetoamide were selected. The pyr4 disruptant confirmed by Southern analysis was named QM9414AP and a single colony was isolated.

To create T. reesei PC-3-7 strain in which bgl2 is disrupted and wild-type bgl2 is expressed, bgl2 disruption cassette and wild-type bgl2 expression cassette were excised from pUbgl2D and pUbgl2 W, respectively, by

AflII treatment and introduced into PC-3-7AKP. Transformed protoplasts were plated on minimal medium without uridine to screen for cells with uridine autotro-phy. Candidate transformants were streaked twice on minimal medium without uridine to obtain stable transformants. Finally, the homologous integrants with single copy of bgl2 disruption or wild-type bgl2 expression cassette were selected by the combination of PCR and Southern hybridization analysis (Additional file 1: Figure S4).

To create QM9414 in which bgl2 is either disrupted or the mutant bgl2 is expressed, the bgl2 disruption cassette and mutant bgl2 expression cassette from pUbgl2D and pUbgl2M, respectively, were introduced into QM9414AP, and transformants were screened as described above.

Preparation of recombinant BGLII in E. coli

The mutated bgl2 cDNA was amplified using total cDNA derived from PC-3-7 as a template and cloned into pET22b (Merck Millipore, MA, USA) to add 6xHis-tag to the C-terminal of the expressed protein. The resulting plasmid pETbgl2F was used as a template for site-directed mutagenesis to create pETbgl2V that expresses wild-type BGLII. E. coli CodonPlus(DE3)-RP (Agilent Technologies, CA, USA) was transformed by pETbgl2V and used as the host of protein expression. E. coli harboring pETbgl2V was cultivated on the medium of Overnight Express™ Autoinduction system (Merck Millipore, MA, USA) according to the manufacturer's instructions. After cultivation, cell lysate of harvested E. coli cells was prepared using xTracter Buffer (Takara-Bio, Shiga, Japan) and recombinant wild-type BGLII (rBGLIIwt) was purified by TALON® Metal Affinity Resins (Takara-Bio, Shiga, Japan) as described in the instructions.

Enzyme assays and measurement of cell biomass

Cellulase activity was determined using 1.0 % carboxy-methylcellulose (CMC; Sigma-Aldrich, MO, USA), filter paper (No. 1, Whatman, LC, UK), and 4-nitrophenyl-ß-D-cellobioside (pNPC; Sigma-Aldrich, MO, USA) as a substrate in 50 mM sodium acetate buffer (pH 5.0) at 50 °C for 15, 60, and 10 min, respectively. The reducing sugars produced during the enzyme reaction were determined by the Somogyi-Nelson method [54, 55] for CMCase and DNS method [56] for FPase. One unit of activity was defined as the amount of enzyme that produced 1 ^mol of reducing sugars per minute in glucose equivalents. In the case of pNPCase, one unit of activity was defined as the amount of enzyme that produced 1 ^mol of 4-nitrophe-nol. The cellobiase activity was determined in 20 mM of phosphate buffer (pH 6.5) at 45 °C with the final concentration of 20 mM of cellobiose as a substrate. The amount

of glucose released by enzyme reaction was determined by Glucose C2 test Wako (Wako Pure Chemical, Osaka, Japan). One unit of cellobiase activity was defined as the amount of enzyme that produced 2 |imol of glucose per minute.

To determine cell volume, an aliquot was taken from the medium and filtered through a small piece of Mira-cloth. Thereafter, the mycelia were placed in a dryer for 12 h and the weight of mycelia was measured.

Protein concentration of the enzyme sample was measured by Bradford method [57] and bovine gamma globulin was used as a standard. Sodium dodecyl sul-fate-polyacrylamide gel electrophoresis (SDS-PAGE) was carried out using 12.5 % polyacrylamide gel slabs as described [58]. Proteins were stained with Coomassie brilliant blue R-250 or with silver stain kit (Wako Pure Chemical, Osaka, Japan). The molecular-mass markers used were the Precision Plus Dual Standard Marker kit for SDS-PAGE (Bio-Rad Laboratories, CA, USA).

Analysis of transglycosylation products

Trichoderma reesei transformants were grown on Avicel for 3 days and mycelia were disrupted using Multi Beads Shocker (Yasui Kikai, Osaka, Japan) in 50 mM of phosphate buffer (pH 6.5). Transglycosylation reaction was performed as follows. 50 |g of cell lysate were mixed in 100 |L of 20 % cellobiose in 50 mM of phosphate buffer (pH 6.5). The reaction mixture was incubated at 40 °C for 24 h. Reaction was stopped by boiling the reaction mixture in a water bath and was filtered using a membrane (0.2 |im, DISMIC®13HP, ADVANTEC, Tokyo, Japan). For recombinant BGLII, 100 |g of purified protein was mixed in 1000 |L of 20 % cellobiose in 50 mM of phosphate buffer (pH 6.5).

The transglycosylation reactants were analyzed by Prominence HPLC system (Shimadzu, Kyoto, Japan) equipped with Shim-Pack ISA-07/S2504 and Shim-Pack guard column ISA. A 10 |L of reactant was injected into the HPLC system with a Shimadzu SIL-20AC HT auto sampler maintained at a constant 4 °C. The column was kept at 65 °C with Shimadzu CTO-20AC column oven. Separation was carried out using a gradient of 0.1 M potassium borate buffer (pH8.0) and 0.4 M potassium borate buffer (pH 9.0) as the mobile phase at a constant flow rate of 0.6 mL/min, with 10 g/L of arginine and 30 g/L of boric acid as reaction reagents at a constant flow of 0.5 mL/min at 150 °C. Sugars were detected with RF-20AXS detector (Ex. at 320 nm, Em. at 430 nm). Glucose, cellobiose, cellotriose, cellotetraose, Laminaribiose,

Gentiobiose, and a-sophorose were used as standard sugars. Transglycosylation products were also analyzed by another HPLC column, Sugar 802KS (Showa Denko, Tokyo, Japan). Separation was performed using H2O as the mobile phase at a constant flow 1.0 mL/min. Sugars were detected with RID-10A detector (Shimadzu, Kyoto, Japan). Glucose and cello-oligosaccharides (di- to tetrasaccharide) were used as standard sugars. For further analysis of transglycosylation products, Thin-layer chromatography was performed as described previously [59]. Transglycosylation reaction was performed using serial amounts of recombinant BGLII (5-0.075 |g by two-fold dilutions) in 50 |L of 20 % cellobiose in 50 mM of phosphate buffer (pH 6.5). Glucose, sophorose, cello-oligosaccharides (di- to hexasaccharide), and laminari-oligosaccharides (di- to heptasaccharide) were used as standard sugars.

RNA extraction, reverse transcription, and first-strand cDNA synthesis

Total RNA was extracted from frozen mycelia by a modified hot-phenol method using TRIzol® (Thermo Fisher Scientific, MA, USA) for additional purification. The lysate was applied to an RNA Spin Mini column (Buckinghamshire, UK) to remove genomic DNA and further purified following the manufacturer's instructions. Total RNA (1 |g) was reverse-transcribed and complementary DNA was synthesized using a Transcriptor First-Strand cDNA Synthesis Kit (Roche Applied Science, Bavaria, Germany).

Quantitative real-time PCR

Quantitative real-time PCR was performed using a Light-Cycler® 480 System (Roche Applied Science, Bavaria, Germany). Amplification reactions were performed in a final volume of 20 |L using a LightCycler® 480 SYBR Green I Master kit (Roche Applied Science, Bavaria, Germany) with 0.5 | M forward primer, 0.5 | M reverse primer, and 2 |l of cDNA (100-fold dilutant of synthesized cDNA). Thermal cycling was conducted under the following conditions: 5 min at 95 °C followed by 45 cycles of 10 s at 95 °C, 10 s at 60 °C, and 10 s at 72 °C. Assays were performed in triplicate with non-amplification controls. The specificity of the PCR amplification was documented by melting curve analysis. Relative expression levels were calculated as the log2 of ACt value by subtracting the Ct value of the housekeeping gene (actl, the gene encoding actin) from that of the target gene. PCR primers used for expression analysis are listed in Additional file 1: Table S2.

Received: 17 September 2015 Accepted: 15 December 2015 Published online: 30 December 2015

Additional file

Additional file 1: Figure S1. Secreted protein by transformant of PC-3-7 and QM9414. A: Avicel cultivation; 10 pL of the supernatant from day 5 of PC-3-7 series and day 6 of QM9414 series were subjected to SDS-PAGE. Gels were stained by Coomassie Brilliant Blue. B: Cellobiose cultivation; 20 pL of supernatant from day 2.5 of PC-3-7 series and day 3 of QM9414 series were subjected to SDS-PAGE. Gel of Avicel cultivation was stained by Coomassie Brilliant Blue and that of cellobiose cultivation was silver stained. Figure S2. Gene expression profile of PC-3-7 transformants on Avicel induction. Genes analyzed were cbh1 (A), egl1 (B), bgl2 (C), and xyr1 (D). Stippled gray bar represents PC-3-7, solid gray bar represents PC-Wbgl2, and the shaded gray bar indicates PC-Abgl2. Values represent the relative expression of each gene normalized to act1 as an internal control. Values represent the means of triplicate experiments. Error bars indicate standard deviations. Figure S3. Effect of bgl2 mutation and disruption in T. reesei QM9414. A: specific activity of the intracellular cellobiase from transformants QM9414, QM-Mbgl2 and QM-Abgl2. Cellobiase activity is derived from the mean of triplicate experiments. Error bars indicate standard deviation. B: HPLC analysis of transglycosylation products using cell-free extracts from transformants. Details are as in Fig. 2. Putative transglycosylation products are indicated by arrows. Figure S4. Southern analysis of transformants for bgl2 analysis. A: the schematic representation of genomic structure of each transformants. B: results of hybridization by gene specific probe for bgl2 and pyr4. Genomic DNA of each transformants was digested by SacI. Lane M represents molecular marker. Table S1. PCR primers used for plasmid construction. Table S2. primers used for real-time quantitative PCR.

Abbreviations

BGL: ß-glucosidase; SNP: single-nucleotide polymorphism; GH family: glycoside hydrolase family; PCR: polymerase chain reaction; TLC: thin-layer chromatography; SDS-PAGE: sodium dodecyl sulfate Polyacrylamide gel electrophoresis.

Authors' contributions

Yosuke Shida conceived of this study and coordinated the manuscript draft and revision. KY, MN, AN, MT, TM executed the experimental work and data analysis. YM helped the interpretation of experimental data. KM, KT, SK helped with genome data analysis. SK, KY, Yasumitsu Sakamoto, NT helped with structural data analysis. WO participated in the design of the study and helped to revise the manuscript. All authors read and approved the final manuscript.

Author details

1 Department of Bioengineering, Nagaoka University of Technology, 1603-1 Kamitomioka, Nagaoka, Niigata 940-2188, Japan. 2 Japan Science and Technology Agency (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan. 3 Department of Genetic Resources Technology, Faculty of Agriculture, Kyushu University, 6-10-1 Hakozaki, Higashi-ku, Fukuoka 812-8581, Japan. 4 Bioproduction Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba Central 6, 1-1-1 Higashi, Tsukuba, Ibaraki 305-8566, Japan. 5 School of Pharmacy, Iwate Medical University, 2-1-1 Nishito-kuta, Yahaba, Iwate 028-3694, Japan. 6 School of Pharmacy, Showa University, 1-5-8 Hatanodai, Shinagawa-ku, Tokyo 142-8555, Japan.

Acknowledgements

The authors thank Dr. Satoshi Nakagawa from Kyowa Hakko Bio Co., Ltd. and Dr. Hirakawa Hideki from Kazusa DNA Res. Inst. for their valuable comments on comparative genome analysis, as well as Dr. Hiroaki Motoyama and Yoshi-yuki Yonetani from Kyowa Hakko Bio Co., Ltd. for their generosity in providing the T. reesei mutants. A part of this research was supported by Grants-in-Aid for Scientific Research from Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan (No. 23603002) and by a SATREPS Grant from Japan Science and Technology Agency (JST) and Japan International Cooperation Agency (JICA).

Competing interests

The authors declare that they have no competing interests.

References

1. Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen J-C, Brown K,

et al. Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: structure and function of a large, enigmatic family. Biochemistry. 2010;40(15):3305-16.

2. Morgenstern I, Powlowski J, Tsang A. Fungal cellulose degradation by oxidative enzymes: from dysfunctional GH61 family to powerful lytic polysaccharide monooxygenase family. Brief Funct Genomics. 2014;13(6):471-81.

3. Schuster A, Schmoll M. Biology and biotechnology of Trichoderma. Appl Microbiol Biotechnol. 2010;87(3):787-99.

4. Bisaria VS, Mishra S. Regulatory aspects of cellulase biosynthesis and secretion. Crit Rev Biotechnol. 1989;9(2):61-103.

5. Kubicek CP, Messner R, Gruber F, Mach RL, Kubicek-Pranz EM. The Trichoderma cellulase regulatory puzzle: from the interior life of a secretory fungus. Enzyme Microb Technol. 1993;15(2):90-9.

6. Vaheri M, Leisola M, Kauppinen V. Transglycosylation products of cellulase system of Trichoderma reesei. Biotechnol Lett. 1979;1(1):41-6.

7. Ilmen M, Saloheimo A, Onnela M-L, Penttilä M. Regulation of cellulase gene expression in the filamentous fungus Trichoderma reesei. Appl Environ Microbiol. 1997;63(4):1298-306.

8. Nogawa M, Goto M, Okada H, Morikawa Y. L-Sorbose induces cellulase gene transcription in the cellulolytic fungus Trichoderma reesei. Curr Genet. 2001;38(6):329-34.

9. Stricker AR, Grosstessner-Hain K, Würleitner E, Mach RL. Xyr1 (xylanase regulator 1) regulates both the hydrolytic enzyme system and D-xylose metabolism in Hypocreajecorina. Eukaryot Cell. 2006;5(12):2128-37.

10. Rauscher R, Würleitner E, Wacenovsky C, Aro N, Stricker AR, Zeilinger S, et al. Transcriptional regulation of xyn1, encoding xylanase I, Hypocrea jecorina. Eukaryot Cell. 2006;5(3):447-56.

11. Aro N, Saloheimo A, Ilmen M, Penttilä M. ACEII, a novel transcriptional activator involved in regulation of cellulase and xylanase genes of Trichoderma reesei. J Biol Chem. 2001;276(26):24309-14.

12. Häkkinen M, Valkonen MJ, Westerholm-Parvinen A, Aro N, Arvas M, Vitikainen M, et al. Screening of candidate regulators for cellulase and hemicellulase production in Trichoderma reesei and identification of a factor for cellulase production. Biotechnol Biofuels. 2014;7(1):14.

13. Saloheimo A, Aro N, Ilmen M, Penttilä M. Isolation of the ace1 gene encoding a Cys(2)-His(2) transcription factor involved in regulation of activity of the cellulase promoter cbh1 of Trichoderma reesei. J Biol Chem. 2000;275(8):5817-25.

14. Aro N, Ilmen M, Saloheimo A, Penttilä M. ACEI of Trichoderma reesei is a repressor of cellulase and xylanase expression. Appl Environ Microbiol. 2003;69(1):56-65.

15. Nitta M, Furukawa T, Shida Y, Mori K, Kuhara S, Morikawa Y, et al. A new Zn(II)(2)Cys(6)-type transcription factor BglR regulates ß-glucosidase expression in Trichoderma reesei. Fungal Genet Biol. 2012;49(5):388-97.

16. Strauss J, Mach RL, Zeilinger S, Hartler G, Stöffler G, Wolschek M, et al. Cre1, the carbon catabolite repressor protein from Trichoderma reesei. FEBS Lett. 1995;376(1-2):103-7.

17. Takashima S, Iikura H, Nakamura A, Masaki H, Uozumi T. Analysis of Cre1 binding sites in the Trichoderma reesei cbh1 upstream region. FEMS Microbiol Lett. 1996;145(3):361-6.

18. Ilmen M, Thrane C, Penttilä M. The glucose repressor gene cre1 of Trichoderma: isolation and expression of a full-length and a truncated mutant form. Mol Gen Genet. 1996;251(4):451-60.

19. Mäntylä A, Paloheimo M, Suominen P. Industrial mutants and recombinant strains of Trichoderma reesei. In: Harman GE, Kubicek CP, editors. Trichoderma and Gliocladium. Bristol: Taylor and Francis Inc; 1998. p. 291-304.

20. Kawamori M, Morikawa Y, Ado Y, Takasawa S. Production of cellulases from alkali-treated bagasse in Trichoderma reesei. Appl Microbiol Biotechnol. 1986;24(6):454-8.

21. Kawamori M, Morikawa Y, Takasawa S. Induction and production of cellulases by L-sorbose in Trichoderma reesei. Appl Microbiol Biotechnol. 1986;24(6):49-453.

22. Morikawa Y, Ohashi T, Mantani O, Okada H. Cellulase induction by lactose in Trichoderma reesei PC-3-7. Appl Microbiol Biotechnol. 1995;44(1-2):106-11.

23. Xu J, Nogawa M, Okada H, Morikawa Y. Regulation of xyn3 gene expression in Trichoderma reesei PC-3-7. Appl Microbiol Biotechnol. 2000;54(3):370-5.

24. Shida Y, Furukawa T, Ogasawara W, Kato M, Kobayashi T, Okada H, et al. Functional analysis of the egl3 upstream region in filamentous fungus Trichoderma reesei. Appl Microbiol Biotechnol. 2007;78(3):515-24.

25. Furukawa T, Shida Y, Kitagami N, Ota Y, Adachi M, Nakagawa S, et al. Identification of the cis-acting elements involved in regulation of xyla-nase III gene expression in Trichoderma reesei PC-3-7. Fungal Genet Biol. 2008;45(7):1094-102.

26. Furukawa T, Shida Y, Kitagami N, Mori K, Kato M, Kobayashi T, et al. Identification of specific binding sites for XYR1, a transcriptional activator of cellulolytic and xylanolytic genes in Trichoderma reesei. Fungal Genet Biol. 2009;46(8):564-74.

27. Rahman Z, Shida Y, Furukawa T, Suzuki Y, Okada H, Ogasawara W, et al. Application of Trichoderma reesei cellulase and xylanase promoters through homologous recombination for enhanced production of extracellular beta-glucosidase I. Biosci Biotechnol Biochem. 2009;73(5):1083-9.

28. Rahman Z, Shida Y, Furukawa T, Suzuki Y, Okada H, Ogasawara W, et al. Evaluation and characterization of Trichoderma reesei cellulase and xyla-nases promoters. Appl Microbiol Biotechnol. 2009;82(5):899-908.

29. Porciuncula JDO, Furukawa T, Shida Y, Mori K, Kuhara S, Morikawa Y,

et al. Identification of major facilitator transporters involved in cellulase production during lactose culture of Trichoderma reesei PC-3-7. Biosci Biotechnol Biochem. 2013;77(3):534-43.

30. Porciuncula JDO, Furukawa T, Mori K, Shida Y, Hirakawa H, Tashiro K, et al. Single nucleotide polymorphism analysis of a Trichoderma reesei hyper-cellulolytic mutant developed in Japan. Biosci Biotechnol Biochem. 2013;77(3):534-43.

31. Lombard V, Ramulu HG, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(D1):D490-5.

32. Saloheimo M, Ylösmäki E, Ward M, Penttilä M. Enzymatic properties and intracellular localization of the novel Trichoderma reesei beta-glucosidase BGLII (Cel1A). Appl Environ Microbiol. 2002;68(9):4546-53.

33. Xu J, Zhao G, Kou Y, Zhang W, Zhou Q, Chen G, Liu W. Intracellular ß-glucosidases CEL1a and CEL1b are essential for cellulase induction on lactose in Trichoderma reesei. Eukaryot Cell. 2014;13(8):1001-13.

34. Häkkinen M, Arvas M, Oja M, Aro N, Penttilä M, Saloheimo M, et al. Reannotation of the CAZy genes of Trichoderma reesei and transcription in the presence of lignocellulosic substrates. Microb Cell Fact. 2012;11:134.

35. Mach RL, Seiboth B, Myasnikov A, Gonzalez R, Strauss J, Harkki AM, et al. The bgl1 gene of Trichoderma reesei QM9414 encodes an extracellular, cellulose-inducible beta-glucosidase involved in cellulase induction by sophorose. Mol Microbiol. 1995;16(4):687-97.

36. Fowler T, Brown RD. The bgl1 gene encoding extracellular beta-glu-cosidase from Trichoderma reesei is required for rapid induction of the cellulase complex. Mol Microbiol. 1992;6(21):3225-35.

37. Lee HL, Chang CK, Jeng WY, Wang AHJ, Liang PH. Mutations in the substrate entrance region of ß-glucosidase from Trichoderma reesei improve enzyme activity and thermostability. Prot Eng Des Sel. 2012;25(11):733-40.

38. Lee B, Richards FM. The interpretation of protein structures: estimation of static accessibility. J Mol Biol. 1971;55(3):379-400.

39. Kidokoro S. Design of protein function by physical perturbation method. Adv Biophys. 1998;35:121-43.

40. El-Gogary S, Leite A, Crivellaro O, Eveleigh DE, El-Dorry H. Mechanism by which cellulose triggers cellobiohydrolase I gene expression in Trichoderma reesei. Proc Natl Acad Sci. 1989;86(16):6138-41.

41. Carle-Urioste JC, Escobar-Vera J, El-Gogary S, Henrique-Silva F, Torigoi E, Crivellaro O, et al. Cellulase induction in Trichoderma reesei by cellulose requires its own basal expression. J Biol Chem. 1997;272(15):10169-74.

42. Aro N, Pakula T, Penttilä M. Transcriptional regulation of plant cell wall degradation by filamentous fungi. FEMS Microbiol Rev. 2005;29(4):719-39.

43. Kubicek CP, Messner R, Gruber F, Mandels M, Kubicek-Pranz EM. Triggering of cellulase biosynthesis by cellulose in Trichoderma reesei. Involvement of

a constitutive, sophorose-inducible, glucose-inhibited beta-diglucoside permease. J Biol Chem. 1993;268(26):19364-8.

44. Fritscher C, Messner R. Cellobiose metabolism and cellobiohydrolase I biosynthesis by Trichoderma reesei. Exp Mycol. 1990;14(4):405-15.

45. Sternberg D, Mandels GR. Induction of cellulolytic enzymes in Trichoderma reesei by sophorose. J Bacteriol. 1979;139(3):761-9.

46. Dos Santos Castro L, Pedersoli WR, Antonieto ACC, Steindorff AS, Silva-Rocha R, Martinez-Rossi NM, et al. Comparative metabolism of cellulose, sophorose and glucose in Trichoderma reesei using high-throughput genomic and proteomic analysis. Biotechnol Biofuels. 2014;7:41.

47. Zhou Q, Xu J, Kou Y, Lv X, Zhang X, Zhao G, et al. Differential involvement of ß-glucosidases from Hypocreajecorina in rapid induction of cellulase genes by cellulose and cellobiose. Eukaryot Cell. 2012;11(11):1371-81.

48. Derntl C, Gudynaite-Savitch L, Calixte S, White T, Mach RL, Mach-Aigner AR. Mutation of the Xylanase regulator 1 causes a glucose blind hydrolase expressing phenotype in industrially used Trichoderma strains. Biotechnol Biofuels. 2013;6(1):62.

49. Kubicek CP. Involvement of a conidial endoglucanase and a plasma-membrane-bound beta-glucosidase in the induction of endoglu-canase synthesis by cellulose in Trichoderma reesei. J Gen Microbiol. 1987;133(6):1481-7.

50. Messner R, Kubicek-Pranz EM, Gsur A, Kubicek CP Cellobiohydrolase II is the main conidial-bound cellulase in Trichoderma reesei and other Trichoderma strains. Arch Microbiol. 1991;155(6):601-6.

51. Hasper AA, Trindade LM, van der Veen D, van Ooyen AJJ, de Graaff LH. Functional analysis of the transcriptional activator XlnR from Aspergillus niger. Microbiology. 2004;150(5):1367-75.

52. Noguchi Y, Tanaka H, Kanamaru K, Kato M, Kobayashi T. Xylose triggers reversible phosphorylation of XlnR, the fungal transcriptional activator of xylanolytic and cellulolytic genes in Aspergillus oryzae. Biosci Biotechnol Biochem. 2011;75(5):953-9.

53. Penttilä M, Nevalainen H, Rättö M, Salminen E, Knowles J. A versatile transformation system for the cellulolytic filamentous fungus Tricho-derma reesei. Gene. 1987;61(2):155-64.

54. Somogyi M. Notes on sugar determination. J Biol Chem. 1926;70(3):599-612.

55. Nelson N. A photometric adaptation of the Somogyi method for the determination of glucose. J Biol Chem. 1944;153(2):375-80.

56. Miller L. Use of dinitrosalicylic acid reagent for determination of reducing sugar. Anal Chem. 1959;31:426-8.

57. Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72(1-2):248-54.

58. Laemmli UK. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970;227:680-5.

59. Takahashi M, Yoshioka K, Imai T, Miyoshi Y, Nakano Y, Yoshida K, et al. Degradation and synthesis of ß-glucans by a Magnaportheoryzae endotrans-glucosylase, a member of the glycoside hydrolase 7 family. J Biol Chem. 2013;288(19):13821-30.

60. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004;32(5):1792-7.

Submit your next manuscript to BioMed Central and we will help you at every step:

• We accept pre-submission inquiries

• Our selector tool helps you to find the most relevant journal

• We provide round the clock customer support

• Convenient online submission

• Thorough peer review

• Inclusion in PubMed and all major indexing services

• Maximum visibility for your research

Submit your manuscript at www.biomedcentral.com/submit

Central