Skip to main content

Lytic polysaccharide monooxygenases and other oxidative enzymes are abundantly secreted by Aspergillus nidulans grown on different starches

Abstract

Background

Starch is the second most abundant plant-derived biomass and a major feedstock in non-food industrial applications and first generation biofuel production. In contrast to lignocellulose, detailed insight into fungal degradation of starch is currently lacking. This study explores the secretomes of Aspergillus nidulans grown on cereal starches from wheat and high-amylose (HA) maize, as well as legume starch from pea for 5 days.

Results

Aspergillus nidulans grew efficiently on cereal starches, whereas growth on pea starch was poor. The secretomes at days 3–5 were starch-type dependent as also reflected by amylolytic activity measurements. Nearly half of the 312 proteins in the secretomes were carbohydrate-active enzymes (CAZymes), mostly glycoside hydrolases (GHs) and oxidative auxiliary activities (AAs). The abundance of the GH13 α-amylase (AmyB) decreased with time, as opposed to other starch-degrading enzymes, e.g., the GH13 AmyF, GH15 glucoamylases (GlaA and GlaB), and the GH31 α-glucosidase (AgdE). Two AA13 LPMOs displayed similar secretion patterns as amylolytic hydrolases and were among the most abundant CAZymes. The starch-active AnLPMO13A that possesses a CBM20 carbohydrate-binding module dominated the starch-binding secretome fraction. A striking observation is the co-secretion of several redox-active enzymes with the starch-active AA13 LPMOs and GHs, some at high abundance. Notably nine AA9 LPMOs, six AA3 sub-family 2 (AA_2) oxidoreductases, and ten AA7 glyco-oligosaccharide oxidases were identified in the secretomes in addition to other non-CAZyme oxidoreductases.

Conclusions

The co-secretion and high abundance of AA13 LPMOs are indicative of a key role in starch granule deconstruction. The increase in AA13 LPMO abundance with culture time may reflect accumulation of a more resistant starch fraction towards the later stages of the culture. The identification of AmyR sites upstream AA13 LPMOs unveils co-regulation of LPMOs featuring in starch utilization. Differential deployment of amylolytic hydrolases and LPMOs over time suggests additional regulatory mechanisms. The abundant co-secretion of distinct AA3 and AA7 oxidoreductases merits further studies into their roles and possible interplay with LPMOs and other enzymes in the deconstruction of starchy substrates. The study reports for the first time the biological significance of LPMOs in starch degradation and the temporal interplay between these and amylolytic hydrolases.

Background

Starch is one of the most abundant renewable biopolymers in nature [1]. Annually two billion tons of starch crops are harvested worldwide, making it an attractive resource for industrial applications such as production of first generation biofuels, pharmaceuticals, textiles, detergents, paper, and food [1]. Starch consists of the two α-glucan polymers, amylose and amylopectin. Amylose is an essentially linear polymer of α-(1,4)-linked glucosyl units, while amylopectin which constitutes 70–80 % of starch granules, is a branched macromolecule having α-(1,4)-glucan chains branched with approximately 5 % α-(1,6)-glucosidic linkages [2]. Despite this chemical simplicity, the α-glucan chains are arranged radially in a supramolecular assembly forming water insoluble granules varying in size (1–150 µm), morphology, crystal-type packing, and crystallinity (15–45 %) [3, 4]. This organization renders starch, especially its crystalline regions, into a relatively challenging substrate for complete enzymatic deconstruction. The extent of starch resistance to enzymatic hydrolysis is correlated to botanical origin and processing, both factors having influence on crystal packing, crystallinity, and supramolecular structure of the starch granule.

The high utilization potential of starch as a renewable biological resource and an industrial feedstock motivates efforts to improve starch hydrolysis yields, particularly from more resistant starch types and for shortening process times. Gain in yields of hydrolysis would have a considerable impact on efficiency and cost of industrial starch processing as well as the reduction of environmental impact of this process [5].

The classical paradigm of polysaccharide degradation by hydrolytic enzymes, which has been valid for decades, was recently re-visited by the discovery of the copper-dependent lytic polysaccharide monooxygenases (LPMOs) [6, 7]. LPMOs catalyze the oxidative cleavage of glycosidic linkages of polysaccharides in the presence of molecular oxygen and an external electron donor, by hydroxylation of either the C1 or C4 carbon of the glycosidic bond [8, 9]. These enzymes play an instrumental role in degradation of recalcitrant crystalline polysaccharides such as cellulose [7, 10] and chitin [6], which renders LPMOs into pivotal tools in industrial biomass conversion [11]. Recently it has been established that some LPMOs oxidize non-crystalline hemicelluloses and soluble cello-oligosaccharides [1215]. Moreover, the recent discovery of starch-active LPMOs [1618], assigned into auxiliary activity family 13 (AA13) in the CAZy database [19], suggest that LPMOs play a role in starch degradation together with amylolytic hydrolases.

Filamentous fungi secrete impressive amounts of hydrolytic enzymes targeting polysaccharides. In addition to glycoside hydrolyses (GHs), several fungi also produce a multitude of LPMOs [20, 21]. In contrast to lignocellulose matrices, little is known regarding the type and composition of enzyme cocktails deployed by fungi for the degradation of starches differing in botanical origin and properties. Notably, the involvement of oxidative enzymes including LPMOs in starch degradation has not been demonstrated in vivo. Here, we provide new insight into the excellent starch-degrading capabilities of the well-studied saprophytic ascomycete Aspergillus nidulans that is taxonomically related to well-established industrial cell factory species such as Aspergillus niger and Aspergillus oryzae [22]. By integrating secretomics and enzyme activity assays, we analyzed temporal changes of the enzymes secreted by A. nidulans to sustain growth on three different starches in the course of 5 days. The data demonstrate differences in growth and secretomes on the selected starch substrates. A common feature of growth on starch was that two AA13 LPMOs including the modular starch-specific enzyme joint to a starch-binding domain of family 20 (CBM20) were among the most abundant CAZymes together with a variety of LPMOs and other oxidative enzymes. This finding suggests that oxidative cleavage of α-glucosidic bonds plays a significant role in starch breakdown. Altogether, the novel insight into enzymatic activities secreted by A. nidulans and related fungi for efficient starch breakdown is relevant for design of enzyme mixtures with enhanced bioconversion efficiencies of starches especially those resistant to hydrolytic degradation.

Results

Starch substrates and fungal growth

To assess the ability of A. nidulans to sense differences in the origin and structure of the starch substrates and to fine-tune the composition of secreted enzymes accordingly, this fungus was grown on wheat, high-amylose (HA) maize, and pea starches, and the resulting secretomes were analyzed. Aspergillus nidulans grew efficiently on wheat and HA maize starch and no intact starch granules were distinguished from the fungal biomass after 5 days, suggesting extensive degradation of both starches. By contrast, growth was poor on pea starch, leaving significant amounts of intact starch granules at culture harvest, which clearly demonstrates important differences due to the botanical origin and properties of the starch on enzymatic deconstruction and growth.

Enzymatic analysis of amylolytic activities

The α-amylase and α-glucosidase activities were measured in the filtered culture supernatants. The average activities of the biological triplicates in different starch media at days 1–5 are shown in Fig. 1. Enzymatic activities were growth-substrate dependent and the highest α-amylase and α-glucosidase activities were measured in the wheat and maize starch culture supernatants, respectively. The α-amylase activity increased to a maximum in 3–4 days and decreased thereafter, with activity maximum (0.21 U/ml) after 4 days in wheat starch (Fig. 1a). By contrast, the α-amylase activity in the pea starch culture supernatants was barely detectable, consistent with the poor growth on this substrate.

Fig. 1
figure 1

Activity of amylolytic hydrolases. Measurement of secreted α-amylase (a) and α-glucosidase (b) activities from Aspergillus nidulans grown on wheat (black), high-amylose maize (gray) or pea starch (white) for 5 days (see “Methods” section). The data are presented as means ± standard deviations of three biological replicates

The α-glucosidase activity increased over time in all samples and showed the highest activity (53 U/ml) in HA maize starch medium on day 5 (Fig. 1b). The glucosidase activity was roughly similar between days 1 and 3 in the pea starch culture, whereas a one-fold increase was observed at days 4–5. By contrast, the α-glucosidase activities in the HA maize and wheat cultures increased steadily until day 5 (Fig. 1b).

Survey of secreted A. nidulans proteins

Filtered supernatants from A. nidulans cultures grown on wheat, HA maize, and pea starch media were analyzed using liquid chromatography combined with tandem mass spectrometry (LC-MS/MS). The analysis of the data set (Additional file 1: Table S1) revealed dynamic secreted protein profiles over the course of 5 days. The theoretical complete proteome of A. nidulans contains 10,556 sequences of which 9.7 % are predicted to be secreted using a combination of three different algorithms. Of the 937 identified proteins in this study, 33 % were predicted to be secreted, which approximately represents 30 % of the theoretical secretome. The identified secreted proteins on days 3, 4, and 5 were assigned to different functional categories including proteases and various carbohydrate-active proteins and clustered both according to abundance and trend related to increase/decrease over time (Additional file 2: Figure S1, Additional file 3: Figure S2, respectively). The number of secreted proteins detected in each culture supernatant varied between 174 (pea starch, day 5) and 221 (HA maize starch, day 1) and generally the number of identified proteins decreased at day 5, as compared to days 3 and 4, but less so in pea (≈4.4 %) followed by wheat (≈7.5 %) and maize (≈9 %) starches (Fig. 2). Approximately 20 % of the secreted proteins were assigned as uncharacterized (lacking characterized homologues). For the remaining secretome, CAZymes (carbohydrate-active enzymes and proteins assigned into the CAZy database, http://www.cazy.org) represented the largest category (Fig. 3) amounting to 44 % of the secretome (Fig. 2).

Fig. 2
figure 2

Overview of the secreted proteins in Aspergillus nidulans starch cultures. The figure depicts the distribution of the secreted proteins from culture supernatants of A. nidulans grown on wheat, high-amylose maize or pea starch on culture days 3–5. Identified carbohydrate-active enzymes (CAZymes) are assigned in the following categories: auxiliary activity (AA), carbohydrate esterase (CE), glycoside hydrolase (GH), polysaccharide lyase (PL). Secreted proteins lacking functionally characterized homologs were assigned as uncharacterized and proteases are shown to highlight their abundance. The remaining proteins belong to a variety of functional categories and are represented as “other” for clarity

Fig. 3
figure 3

Heat map comparison of expression patterns of CAZymes detected after 3–5 days growth of Aspergillus nidulans on minimal media supplemented with: high-amylose (HA) maize, wheat, or pea starch. The colors in the heat map indicate the label-free quantification (LFQ) intensity reported by MaxQuant ranging from 2 × 106 (light green) to 4 × 1011 (light red). Missing values were imputed from a normal distribution located at the quantification limit. Glycoside hydrolases (GHs) are in red, auxiliary activities (AAs) in blue, polysaccharide lyases (PLs) in light brown, carbohydrate esterases (CEs) in yellow, and other proteins in gray. Secretion prediction is a combination of SignalP, Phobius, and WolfPSort where at least two algorithms had to agree

Approximately 30 % of the 478 annotated CAZymes were identified in the secretomes despite growth on purified starches and not crude plant material. The data revealed that most of the identified CAZymes are GHs or AAs. Generally, the total number of identified CAZymes was lower in the pea starch cultures (Fig. 2). Aspergillus nidulans secreted multiple enzymes (4–10 enzymes) from certain CAZy families, e.g., GH3, GH16 and GH43, and CE10 as well as AA3, AA7, and AA9, with the three latter families harboring oxidative activities (Additional file 4: Table S2). Some ancillary carbohydrate module (CBM) families were identified including four CBM20s which mediate starch-binding in addition to cellulose-specific (CBM1 and CBM6), different β-glucans (CBM6, CBM24, and CBM43), chitin (CBM18), and arabinofuranosyl (CBM42)-binding modules.

Cell wall-degrading activities from GH5, GH6, GH7, GH53, and GH62 were present only in wheat and HA maize cultures. Beyond this, only a few families were not identified in all samples, including GH35 which was only observed in HA maize and pea starches, and GH24 and GH95 only identified in HA maize and pea starch-induced samples, respectively (Additional file 4: Table S2). Notably, a GH13 α-amylase (AmyF, Q5B7U2), a GH31 α-glucosidase (AgdB, G5EB11), and an AA13 LPMO (AnLPMO13B, Q5B027) were the most abundant CAZymes in all samples (Fig. 3). The analysis revealed the secretion of six different putative LPMOs belonging to families AA9 and AA13 in all starch samples, and one AA11, which was only found in the pea starch culture on day 3 (Table 1; Additional file 4: Table S2).

Table 1 Predicted secreted lytic polysaccharide monooxygenases encoded in the Aspergillus nidulans genome

Starch-degrading enzymes

α-Amylases of GH13, glucoamylases of GH15, and α-glucosidases of GH31 are candidates for starch degradation. Indeed, seven of the twelve predicted extracellular enzymes of these families were identified in the secretome (Table 2; Fig. 4). Three putative GH13 α-amylases (G5EB45, G5EAT0, and Q5B7U2), and three GH31 α-glucosidases (G5EB03, G5EB11, and Q5BET9) were present in all samples on all days. A GH15 α-glucoamylase (C8VLL3) was identified in all wheat and HA maize starch medium and only on day 4 and 5 in pea. A second glucoamylase was also identified (GlaB, Q5AWC8) (Table 2) and included based on a previous study showing its purification from A. oryzae cultures [23], despite its lack of a canonical secretion signal. Moreover, the recently functionally characterized starch-degrading AA13 LPMO [17] was identified in addition to a second AA13 enzyme (Table 1; Fig. 4). Interestingly, the abundance of the GH13 α-amylase is having a CBM20 (G5EAT0) decreases over time (Fig. 4; Additional file 3: Figure S2). Otherwise, all amylolytic proteins were highly abundant at all time points or showed an increase over time (Fig. 4; Additional file 3: Figure S2).

Fig. 4
figure 4

Heat map identical to Fig. 3, but filtered to contain amylolytic hydrolases and AA enzymes detected after 3–5 days growth of Aspergillus nidulans on minimal media supplemented with: high-amylose (HA) maize, wheat, or pea starch. Highlighted in green are AA13s, in purple—AA9s and AA11s, in blue—AA3s and AA7s, in red—hydrolytic enzymes associated with starch degradation, and in gray—other proteins

Table 2 Predicted secreted enzymes assigned in CAZy families implicated in starch degradation by the Aspergillus nidulans genome

Redox-active enzymes

The A. nidulans genome encodes 13 LPMOs: 9 AA9s, 2 AA11s, and 2 AA13s (Table 1). Two AA9 are assigned in clade LPMO1, one in LPMO2, five in clade LPMO3, and one in an unclassified cluster, as phylogenetically defined by Vu et al. [8]. None of the five clade 3 LPMOs were identified, whereas the remaining four enzymes were observed, albeit at varying abundance (Table 1; Fig. 4). The two AA13 LPMOs were among the most abundant proteins in the culture supernatants, AnLPMO13B (Q5B027) seemingly at higher abundance than AnLPMO13A (Q5B1W7). The N-terminal histidine residue of both these AA13 enzymes seemed to be methylated in virtually all samples (only exception being the day 2,3 replicates from the pea culture) based on the identification of the “H(Me)GYLTVPASR” peptide (best MASCOT score = 78) that is identical in both enzymes, while the non-methylated counterpart could not be detected. Notably, the abundance of both proteins increased slightly over time (Fig. 4; Additional file 3: Figure S2).

In addition to LPMOs, the secretomes also contain a variety of other redox-active enzymes, including a catalase, laccase, and thioredoxin reductase as well as members of AA3 and AA7 (Fig. 4; Additional file 3: Figure S2). The six (out of 18) putative AA3 enzymes secreted by A. nidulans belong to sub group 2 that accommodates aryl-alcohol and glucose oxidases as well as dehydrogenases [24]. Of these, two (Q5AV48 and Q5AUN2) were highly abundant, whereas the initially low abundance C8VDT4 shows a substantial increase over time (Fig. 4; Additional file 3: Figure S2). The remaining AA3s are only present at low abundance. A. nidulans secretomes show the presence of 10 AA7 putative glyco-oligosaccharide oxidases (GOOs), three of which (Q5AY23, Q5B7X9, and C8VCU1) were highly abundant in most samples, the two former and the latter showing a decreasing and increasing trend over time, respectively (Fig. 4; Additional file 3: Figure S2).

The secretomes contain several non-CAZyme redox-active enzymes. As a matter of fact, the most abundant protein in all collected supernatants (except day 5 maize and day 3 pea) is Catalase B (CatB, P78619), which also shows a slight decrease in abundance over time in wheat and maize, and an opposite trend in pea (Fig. 4; Additional file 5: Tables S3–S5). Thioredoxin reductase is abundant in most samples (among top 20 proteins in wheat and pea starch samples, Additional file 5: Tables S3 and S5) and shows a decreasing trend in the wheat and maize samples like CatB. A laccase (Q5BEX9) is also present in the secretome at low abundance, and more so over time.

Plant cell wall-degrading non-oxidative enzymes

Ten putative β-glucoside and cellulose-degrading enzymes were identified. A GH1 (Q5AR97) and five putative GH3 β-glucosidases (Q5B5S8, Q5B6C6, Q5AWD4, Q5B7X0, and Q5B9F2) were identified in all samples. Two of these enzymes (Q5AR97, Q5B7X0) appeared only after 4 days in the pea starch. Notably, cellobiohydrolases (GH6) and endoglucanases (GH5_7) active on cellulose were not identified in the pea starch culture supernatant. Two GH5 (Q5AUG2, Q5BDU5) were only present in the maize, whereas the GH6 (C8VSG6) and GH7 (Q5B7R2) were also identified in wheat starch cultures. Twelve putative hemicellulose-degrading enzymes were identified including a GH10 (Q00177) and a GH11 (P55332) β-1,4-endoxylanases and four putative GH43 β-xylosidases/α-l-arabinofuranosidases (C8VCT5, C8VKG9, Q5AUM3, and Q5AV99). α-Arabinofuranosidases of GH54 (O74288) and GH62 (Q5AUX2) were also identified, the latter only in wheat and HA maize starch media together with a GH53 β-1,4-endogalactanase (Q5B153). Moreover, a putative GH93 exo-arabinanase (Q5BBM0), a GH27 (Q5AX28) and a GH36 (Q5AU92), both the latter being α-1,6-galactosidases were identified. Putative pectin-degrading enzymes including a GH2 β-glucuronidase (Q5BAN5), a GH105 unsaturated rhamnogalacturonan hydrolase (Q5AQP7), and two putative GH28 exo-polygalacturonases (Q5ASG9 and Q873X6) were detected, the latter only in HA maize starch cultures. Finally, a GH35 putative β-galactosidase (Q5BFC4) was detected in HA maize and pea starch cultures.

Fungal cell wall active enzymes

Fifteen secreted enzymes putatively involved in fungal cell wall modification and degradation were identified including six GH16 glycanases (C8VUN8, Q5AY11, Q5AYL0, Q5BAP5, Q5BGT5, and Q5B4L5), a GH17 β-1,3-endoglucanase (Q5AUT0), two putative GH55 β-1,3-exoglucanase (Q5B3Q5 and C8VQV2), three GH72 β-1,3-glucanosyltransferases (C8VSK8, Q5AVM3, and Q5AW19), a GH81 endo-β-1,3-endoglucanase (C8VT57), and a putative GH20 N-acetylglucosaminidase (G5EB27). The putative GH71 α-1,3-glucanase (C8VKQ6) was detected only on days 4 and 5 in all cultures. In addition, three α-mannan-degrading enzymes, one GH47 (Q5BF93), and two GH92 (C8VSU5 and Q5BAV5) were identified.

Purification of starch-binding enzymes from wheat starch cultures

β-Cyclodextrin affinity chromatography was used to capture the fraction of A. nidulans secreted enzymes that possess affinity for starch. The elution fraction from this purification was analyzed by SDS-PAGE (Fig. 5) and dominant protein bands were identified using MALDI-TOF mass spectrometry (Table 3). The common feature to all proteins identified on the gel is the presence of CBM20 [25]. The three most intense bands contain the GH15 α-glucoamylase (GlaB, Q5AWC8), the GH13_1 α-amylase (G5EAT0), and the biochemically characterized starch-active LPMO AnAA13A (Q5B1W7, Table 1), that appeared as the most abundant protein in this starch-associated secretome fraction. The CBM20 of the α-amylase was also identified in two lower molecular mass bands on the gel, which is likely due to proteolytic cleavage (Fig. 5).

Fig. 5
figure 5

Purification of starch-binding enzymes from the secretome. β-cyclodextrin-sepharose affinity chromatography purification of the Aspergillus nidulans wheat starch culture supernatant harvested after 5 days (lane II). Three dominant bands were identified using mass spectrometry (MALDI-TOF MS, see “Methods” section). The bands were identified as glucoamylase (1), α-amylase (2 and 2*-CBM20 domain), and AnLPMO13A (3). Molecular mass marker in kDa (lane I)

Table 3 Identification of β-cyclodextrin-sepharose affinity chromatography purified secreted proteins in Aspergillus nidulans culture grown on wheat starch for 5 days

Discussion

The abundance of starch from terrestrial biomass is only surpassed by cellulose. An excess of two billion tons of starch are annually produced from cereals and coarse grains and another 700 million tons from roots and tubers (The food and agriculture organization of the United Nations [FAO]; http://www.fao.org/). The bulk of extracted starch from these crops is used for non-food industrial applications including first generation biofuel that still contributes more than two thirds of the total biofuel production [26]. Technical and economic issues continue to hamper a larger transition to second generation lignocellulose-based feedstocks [27], which together with the large investments made in first generation production facilities is likely to maintain the demand for starch as a feedstock in the time to come. Therefore, improving the efficiency of starch depolymerization in industrial applications has significant economic and environmental potential.

The synergistic action of LPMOs with hydrolases [10, 28] in the breakdown of recalcitrant polysaccharides has spurred interest in the exploitation of these enzymes in biofuel production [11]. Several proteomics studies revealed that LPMOs are highly co-secreted with hydrolases during fungal growth on plant-derived biomass or purified polysaccharides [20, 29]. The oxidative activity of starch-specific LPMOs from Neurospora crassa [16] and A. nidulans [17] toward resistant retrograded starch was recently demonstrated and the presence of the A. nidulans LPMO elicited a 100-fold increase in the release of maltose by a commercial β-amylase on the same substrate [17]. This in vitro activity suggests that starch-active AA13 LPMOs may play an important role in starch depolymerization by distinct fungi. Currently, there are no data on the involvement and biological significance of LPMOs in starch degradation. In this study we investigate the protein inventory secreted by A. nidulans during growth on three different starches differing in botanical origin and properties and show that starch-specific LPMOs are highly abundant in all starch secretomes together with LPMOs of different specificities and other oxidative enzymes.

A. nidulans secretomes are fine-tuned to the different starch matrices

Fungal growth was far better on cereal starches from wheat and maize, while the pea starch was a poor growth substrate. The overall secretion profiles differed between the three substrates, but were more similar between the botanically closer cereal starches as compared to pea (Fig. 3; Additional file 2: Figure S1). This was also supported by enzyme assays that revealed much higher α-glucosidase and α-amylase activities in the cereal starch cultures and a substantial increase of α-amylase activity at day 3, which was not observed in the pea starch culture (Fig. 1). This is also consistent with large differences in the 20 most abundant proteins (Additional file 5: Tables S3–S5). The diversity of identified CAZymes was the highest in HA maize starch and the lowest in pea starch secretomes that contained a larger relative proportion of unknown enzymes (Fig. 2). These data highlight the fine-tuning of the secretomes to each starch substrate. The three used starches are different in their fine structures and composition. The HA maize starch has the highest amylose content with approximately 70 % (w/w), followed by pea starch (>45 % w/w) and wheat (25 % w/w). HA maize starch granules are the smallest (5–25 µm), followed by pea (5–45 µm) and wheat (20–35 µm) [3, 4]. Given the chemical simplicity of starch, it is unlikely that the secretomes merely reflect the differences in starch structure per se, but also the presence of other non-starch components, including polysaccharides, lipids, and proteins, which may impede the action of amylolytic enzymes [30]. The starch granules and protein storage vacuoles in monocot cereals like wheat and maize reside in the endosperm, while in pea the cotyledon tissue serves the same storage function [1, 31]. Pea starch isolation has been reported to be more difficult due to the presence of fine fiber and insoluble flocculent proteins, significant amounts of which remain associated with the starch [32]. Indeed, numerous proteases are identified in the secretomes and more so in the pea than the cereal starch cultures (Fig. 2; Additional file 2: Figure S1). Pea starch is generally less well studied with respect to microbial degradation and the reasons for the poorer growth on this substrate are unclear. Starch from HA maize has been shown to be the least efficiently degraded by an α-amylase from Aspergillus fumigatus, followed by starch from pea [33]. The excellent growth of A. nidulans on the more resistant HA maize starch supports the possibility that non-starch components in the pea starch matrix rather than the recalcitrance of the starch is responsible for the poor growth. Indeed, legume cotyledons including pea are known to harbor a range of protease and hydrolase inhibitors [34], which may represent more challenging growth conditions and delay the deconstruction of the starch.

The major components in the cereal endosperm cell walls are arabinoxylans and mixed linkage glucans, in addition to small amounts of cellulose, xyloglucan, and pectins [35]. By contrast, xyloglucan is the main structural polysaccharide in cotyledon cell walls [36]. This is consistent with the presence of cellulases of GH5, GH6, and GH7, in addition to GH62 arabinoxylan arabinofuranosidase exclusively in the cereal starch cultures (Fig. 3; Additional file 4: Table S2). Xylose has been suggested to be the inducer of several polysaccharide-degrading enzymes in aspergilli [37]. Possibly the co-secretion of a range of hemicellulose-degrading enzymes including a GH10 xylanase with amylolytic enzymes (Fig. 3; Additional file 3: Figure S2) is suggestive of the presence of xylose derived either from arabinoxylan in cereals or xyloglucan remaining in the legume pea starch preparations. This co-secretion may contribute to the degradation of cell wall components and render starch more accessible to the action of amylolytic enzymes. An “anti-caging” effect of xylanases and other cell wall hydrolases has been observed in vitro starch digestion studies [38, 39]. Also the secretion of five LPMOs of AA9, which feature in the breakdown of cellulose and hemicellulose (including an enzyme with a putative cellulose binding module, Table 1) supports the presence of cell wall components in the starch matrices [7, 13, 14]. These findings motivate further studies on the effect of low concentrations of key hemicellulose polysaccharide-degrading enzymes on the yields of starch degradation.

The hydrolytic starch-degrading machinery

The genomic and enzymatic capabilities of A. nidulans with regard to complex polysaccharide degradation [19, 40, 41] manifest the saprophytic lifestyle of this fungus. The starch degradation model in fungi has been based on the biochemically well-described amylolytic hydrolases from aspergilli [33, 42]. The present study is the first to report high resolution secretome data in response to growth on different starches, providing insight on the impact of starch botanical origin on growth and protein secretion profiles.

The data revealed the secretion of a core of starch hydrolytic battery of A. nidulans comprising three α-amylases, three α-glucosidases, and two glucoamylases (Table 2). The absence of AmyD and AmyE of GH13_1 from the secretome (Table 2) is in excellent agreement with the fungal cell wall remodeling role of these enzymes [43]. An intriguing observation is the decrease in abundance of the GH13-CBM20 (AmyB, G5EAT0) over time (Fig. 4; Additional file 3: Figure S2). The capture of the CBM20 proteolytic degradation product of this enzyme from the starch-binding fraction of the secretome (Fig. 5) proves proteolytic degradation has occurred, which is also in agreement with the decrease of α-amylase activity on day 5 in cereal starches, despite the slight increase in abundance of the other α-amylase (AmyF) not possessing a CBM20. By contrast, a pattern of increased abundance is observed for the CBM20-containing LPMO (AnLPMO13A, Q5B1W7) and GH15 glucoamylases (Fig. 4; Additional file 3: Figure S2). This expression pattern suggests that the starch-binding α-amylase AmyB (G5EAT0) plays an instrumental role in the initial attack on granular starch. It is tempting to speculate that the endo-acting activity of AmyB creates sites for attack by both glucoamylases and α-glucosidases that release glucose from the non-reducing ends of starch and solubilized fragments thereof. The high abundance and relative increase in AA13 LPMOs with time, observed here for the first time, indicate that these activities are of great importance in starch deconstruction. The possible increase of resistant starch structures that escape hydrolytic degradations during growth may justify maintained or increased levels of these enzymes that have been shown to boost the degradation of resistant starch, toward the later stages of growth, which will be further discussed below.

Abundant deployment of LPMOs in starch degradation by A. nidulans

Several studies investigating the secretomes of plant biomass-degrading microorganisms have reported the abundant presence of LPMOs, indicating their importance in biomass conversion [20, 29]. However, the significance of LPMOs in starch utilization in vivo has not been addressed. Two studies recently reported oxidative cleavage of starch by fungal AA13 LPMOs [16, 17]. Indeed, LPMOs from both family AA9 and AA13 were identified in the starch secretomes, the latter family being the most abundant. As a matter of fact, one of the AA13 LPMOs, AnLPMO13B (Q5B027) is among the top four most abundant proteins at all time points in the wheat and maize starch cultures (Fig. 4; Additional file 5: Tables S3, S4), indicating a prominent role in starch degradation. The other AA13 LPMO (AnLPMO13A), which has a CBM20 appended, is also abundant in the secretomes of cereal starches (always among the top 10 %). The abundance of AnLPMO13A may be underestimated as compared to AnLPMO13B, due to the preferential CBM20-mediated binding of the former onto granular starch (concluded form β-cyclodextrin capture, Fig. 5). AnLPMO13A has already been demonstrated to cleave amylose chains and to boost hydrolysis of retrograded starch [17], whereas no activity data are available for AnLPMO13B. However, the two AA13 LPMOs are highly similar (catalytic modules share 81 % sequence identity), indicating that they also share substrate specificity. A notable difference is that AnLPMO13B possesses a 75 amino acid C-terminal extension of unknown function as compared to the long linker and the C-terminal CBM20 of AnLPMO13A. The absence of AnLPMO13B in secretome fraction captured by β-cyclodextrin affinity suggests that the affinity of this enzyme to insoluble starch is much lower than AnLPMO13A that possesses a CBM20. These differences indicate that the AA13 LPMOs may target different parts of the starch structure, yielding complementary activities. This has been observed for the conserved cellulose-targeting LPMO pair encoded by Streptomyces coelicolor, which contributed synergistically when combined in cellulose degradation [44]. Interestingly no activity of the AA13 catalytic module lacking the CBM20 was detected [17]. This highlights the role of the CBM20 in targeting the enzyme to the insoluble starch structures.

Although less abundant than the family AA13 LPMOs, four AA9 LPMOs and one AA11 LPMOs were identified in the cultures (Table 1; Additional file 4: Table S2). So far, AA9 LPMOs have only been shown to depolymerize cellulose and hemicellulose substrates. Thus, it is reasonable to speculate that these enzymes are secreted by the fungus to aid gaining access to the starch granules that in vivo are shielded by a protective lignocellulose layer, e.g., grain bran or to disrupt cell walls in the starchy endosperm tissue. Interestingly, the five clade 3 AA9 LPMOs of the A. nidulans genome were not represented in the secretomes (Table 1). Clade 3 LPMOs are proposed to have mixed activity, yielding both C1- and C4-oxidized products, whereas clade 1 members oxidize the C1 carbon and clade 2 the C4 carbon of the substrate. The significance of this expression pattern is not known, but may be related to the type of substrate the fungus is exposed to.

One of the two AA11 LPMOs is observed in the secretomes, albeit at low abundance. The enzyme shows only moderate similarity to the chitin-degrading AoAA11 (51 %) [45], precluding reliable functional assignment. Assuming chitin specificity (like AoAA11), it is possible that the function is related to fungal cell wall remodeling.

The role of non-LPMO redox-active enzymes in starch degradation

It is well established that the sub-family 1 members of the glucose, methanol, choline (GMC) oxidoreductases (AA3s), i.e., the cellobiose dehydrogenases (CDHs), can provide electrons to LPMOs [28, 46]. Maltodextrins are known as poor substrates for CDHs [47], thus it may be that other redox-active enzymes play this role during starch degradation. Indeed, multiple enzymes from the AA3 family and the AA7 family (glyco-oligosaccharide oxidases with one enzyme shown to be active on maltose other malto-oligosaccharides) are identified in the culture supernatant (Fig. 4).

Both aryl-alcohol oxidases and glucose dehydrogenases, both assigned into AA3_2, have been, very recently, shown to be able to transfer electrons to LPMOs [24, 48]. Two (C8VDT4 and Q5B9S6) of the six identified putative AA3_2 enzymes in the present study are not currently assigned in CAZy, whereas the remaining four sequences populate clusters lacking any biochemically characterized enzymes in a recent phylogenetic analysis of AA3_2 [49], which precludes a reliable assignment of their activity. Notably, the two identified AA3_2 putative enzymes Q5AUN2 and Q5B624 populate a single cluster adjacent to characterized glucose dehydrogenases in the same phylogenetic analysis [49]. The identified AA7 enzymes (and possibly some of the AA3_2) are predicted oxidases (not dehydrogenases as the CDHs), indicating that one of their main products is H2O2. A possible role for enzymes generating H2O2 is to provide lignin-degrading enzymes like laccases and peroxidases with hydrogen peroxide, required for lignin depolymerization. On the other hand, no or little lignin is present in the substrates provided in the present experiments, possibly indicating a different role of these enzymes. Biochemical characterization of the identified AA7 and AA3 enzymes is warranted to verify their roles.

The H2O2 generated from the activity of AA3 and AA7 oxidases on products derived from depolymerization of starch (or other components in the substrate matrices) provides a rationale for the dominance of catalase B (CatB) in most secretomes. Catalase catalyzes the conversion of H2O2 to O2 and H2O, thus protecting the fungus from H2O2 toxic effects (i.e., hydroxyl radicals emerging from Fenton chemistry). Indeed, the expression of CatB has previously been demonstrated to be induced by H2O2 or H2O2-generating conditions [50], and the enzyme activity was suggested to protect the fungus from the toxic side-products encountered in aerobic growth. It should be noted that LPMOs, which represent some of the most abundant proteins in the cultures, are capable of substantial H2O2 production [51]. A second highly abundant protein observed in the fungal cultures is thioredoxin reductase (Fig. 4; Additional file 5: Tables S3–S5), a protein also related to detoxification of reactive oxygen species (ROS) [52]. Thus, it seems that the fungus is actively protecting itself from ROS during aerobic degradation of biomass. These findings could inspire future research on design of enzyme cocktails for biomass depolymerization, i.e., protecting enzymes from toxic byproducts formed by oxidative enzymes. As a matter of fact, one study has already shown the advantageous use of catalase in enzyme cocktails containing high levels of LPMOs [53].

AmyR regulation of amylolytic hydrolases and LPMOs

The AmyR transcription factor, that is activated by binding isomaltose, has been shown to regulate the expression of the α-glucosidases agdA, agdB, and agdE agdF, the α-amylases amyA, amyB, and amyF, and the glucoamylase glaB genes in A. nidulans [54]. All these enzymes were indeed identified in the present study (Table 2). Analysis of the A. nidulans genome revealed the presence of the AmyR consensus sequence 5′-CGGN8CGG-3′ at around −300 basepairs (−318 to −305) in the promotor region of the AnLPMO13A and around −600 basepairs (−612 to −609) for AnLPMO13B. This is the first time a common regulatory link has been identified between GHs and LPMOs. Interestingly, hierarchical clustering analysis revealed that both AA13 enzymes have very similar secretion patterns as distinct amylolytic enzymes under the control of AmyR (Fig. 4). These first data on the temporal distribution of starch-degrading enzymes suggest a sophisticated regulatory mechanism, whereby the fungus deploys specific enzymes at different stages of the starch growth. As mentioned above, it is tempting to speculate that the increased secretion of starch-specific AA13s correlates with the later stages of the culture, where the more resistant structures in the starch substrate accumulate. Further studies are needed to unveil the precise preference of the two types of AA13 enzymes observed in this study.

Conclusions

This is the first study that uses high resolution secretomics analysis to probe the temporal changes of protein profiles of the model saprophytic ascomycete A. nidulans growing on different starchy substrates. The data, for the first time, unveil a conspicuous abundance of AA13 LPMOs in the secretomes, which is suggestive of an instrumental role of these enzymes in starch degradation. Another novel finding is the identification of binding sites of the transcriptional regulator AmyR, which established a co-regulatory link between GHs and LPMOs featuring in starch degradation. Beyond a core amylolytic machinery, the secretomes were clearly correlated to the starch type used for growth and numerous CAZymes, proteases, and oxidoreductases were secreted. A possible rationale for this is targeting non-starch component in the substrate matrix. The abundance and co-secretion patterns of LPMOs, AA3, AA7, and other oxidative enzymes is particularly intriguing and merits further work to understand the role of these enzymes. This novel insight promotes our understanding of the degradation of different starches, including those of more resistant nature, and inspires formulation of better commercial enzyme cocktails for more efficient exploration of this important biomass resource.

Methods

Carbohydrate substrates and assay chemicals

The wheat starch was from Sigma-Aldrich (S5127, unmodified), pea and high-amylose (HA) maize starches were from KMC (Brande, Denmark). The starches for growth experiments were not autoclaved to minimize changes of the native starch granule structure, but washed twice in 70 % ethanol and water and subsequently pelleted by centrifugation (14,000×g, 5 min) before resuspension in autoclaved growth media or as substrates for enzymatic assays. p-Nitrophenyl-α-d-glucopyranoside (PNPG) and p-nitrophenol (PNP) were from Sigma-Aldrich, and insoluble Blue Starch was a custom preparation from Pharmacia (Uppsala, Sweden).

Fungal strain and culture conditions

A. nidulans strain FGSC A4, (FGSC, Kansas City, MO), was pre-grown on malt extract agar plates containing 1 % wheat, HA maize or pea starches for 5 days at 30 °C until new mycelia were formed. Mycelial plugs were used to inoculate 1 L of minimal medium containing 1 % (w/v) carbon source with a start pH of 6.5 in 3 L baffled shake flasks. The minimal medium contained, per liter, 6 g NaNO3, 0.52 g KCl, 0.52 g MgSO4·7H2O, 1.52 g KH2PO4, and 2 ml of Hutner’s trace elements. Hutner’s trace elements contained, per liter, 2.2 g ZnSO4·7H2O, 1.1 g H3BO3, 0.5 g MnCl2·4H2O, 0.5 g FeSO4·7H2O, 0.16 g CoCl2·6H2O, 0.16 g CuSO4·5H2O, 0.11 g (NH4)Mo7O24·4H2O, and 5 g EDTA. The cultures were incubated on a rotary shaker (150 rpm) at 30 °C for 5 days. All experiments were performed in biological triplicates. The medium was supplemented with antibiotics (100 µg/ml ampicillin and 34 µg/ml chloramphenicol) to inhibit growth of possible bacterial contamination from the starch substrates.

α-amylase and α-glucosidase activities

Both enzymatic activities were measured and averaged from the biological triplicates.

α-amylase activity of the fungal cultures was assayed toward insoluble Blue Starch (iBS). The activity was measured using iBS (6.25 mg/ml) suspended in 10 mM MES, pH 6.5. The reaction mixture (900 µl) was incubated for 15 min at 37 °C after addition of the culture filtrate (100 µl). The reaction was stopped by addition of 0.5 M NaOH (200 µl). After centrifugation (4000×g, 3 min) 200 µl supernatant was transferred to a 96 well microtiter plate and A 620 values were used to determine enzyme activity. One activity unit (U) was defined as the amount of enzyme that leads to an increase in A 620 of 1 absorbance unit in the reaction mixture under these experimental conditions.

The α-glucosidase activity was determined using p-nitrophenyl α-d-glucopyranoside (PNPG) as the substrate. The activity in culture filtrates was assayed toward 2 mM PNPG in 10 mM MES, 0.005 % Triton X-100, pH 6.5 and 10 μl culture supernatant in 50 μl reactions, at 37 °C for 60 min. The reaction was stopped by adding 1 M Na2CO3 (200 μl). The amount of released p-nitrophenol (PNP) was measured spectrophotometrically at A 410 using PNP (0–2 mM) as a standard. One unit of activity (U) was defined as the amount of enzyme that released 1 μmol/min of PNP at the given conditions.

Preparation of secretome samples and mass spectrometry analysis

Culture filtrates containing the fungal secretomes were separated from mycelia and residual insoluble starch by centrifugation (15,000×g, 15 min, 4 °C) and filtration (0.45 µm hydrophilic PVDF membranes; Millipore) and stored at −20 °C until further use. Samples (10 ml) from three biological replicate cultures were collected from day 3 to 5 under sterile conditions. Protein concentrations of the culture filtrates were determined using a Bradford assay (Bio-Rad Laboratories, Hercules, CA) with a BSA standard according to the manufacturer’s instructions.

Proteins of culture filtrates (2 ml) were precipitated by direct addition of 500 µl cold 50 % trichloroacetic acid (TCA), incubated at 4 °C for 8 h, centrifuged (15,000×g, 15 min, 4 °C) to pellet precipitated protein and washed with 300 µl ice-cold 10 mM HCl in 90 % acetone. After centrifugation (15,000×g, 15 min, 4 °C), pellets were dried, redissolved in 100 µl 20 mM Tris–HCl at pH 8, reduced with 10 mM (final concentration) dithiothreitol (DTT) and incubated at room temperature for 30 min. Alkylation was performed by adding 15 mM (final concentration) iodoacetamide (IAA) and incubating at room temperature for 30 min in the dark. The proteins were digested with 20 µl 12.5 ng/µl sequencing-grade modified trypsin (Promega, Madison, WI) and incubated at 37 °C for 16 h. Trypsination was stopped with 0.5 % (final concentration) trifluoroacetic acid (TFA). Samples were dried in SpeedVac and purified with ZipTip C18 pipette tips (Merck Millipore, Cork, Ireland) according to manufacturer’s instructions. Purified samples were dried and dissolved in 2 % acetonitrile (ACN) and 0.1 % TFA mix.

The peptides were analyzed in two technical replicates using a nanoHPLC-MS/MS system consisting of a Dionex Ultimate 3000 RSLCnano (Thermo Scientific, Bremen, Germany) connected to a Q-Exactive hybrid quadrupole-orbitrap mass spectrometer (Thermo Scientific) equipped with a nano-electrospray ion source. Samples were loaded onto a trap column (Acclaim PepMap100, C18, 5 µm, 100 Å, 300 µm i.d. ×5 mm, Thermo Scientific) and backflushed onto a 50 cm analytical column (Acclaim PepMap RSLC C18, 2 µm, 100 Å, 75 µm i.d., Thermo Scientific). At the start, the columns were in 96 % solution A (0.1 % (v/v) formic acid), 4 % solution B (80 % (v/v) ACN, 0.1 % (v/v) formic acid). Peptides were eluted using a 90 min gradient from 4 to 13 % (v/v) solution B in 2 min, 13–45 % (v/v) in 70 min and finally to 55 % B in 5 min before the wash phase at 90 % B, all at a flow rate of 300 nl/min. To isolate and fragment the ten most intense peptide precursor ions at any given time throughout the chromatographic elution, the Q-Exactive mass spectrometer was operated in data-dependent mode to switch automatically between orbitrap-MS and higher-energy collisional dissociation (HCD) orbitrap-MS/MS acquisition. The selected precursor ions were then excluded for repeated fragmentation for 20 s. The resolution was set to R = 70,000 and R = 35,000 for MS and MS/MS, respectively. For optimal acquisition of MS/MS spectra, automatic gain control (AGC) target values were set to 1,000,000 charges and a maximum injection time of 128 ms.

MS raw files were imported into MaxQuant [55, 56] version 1.4.1.2 and proteins were identified and quantified using the MaxLFQ algorithm [57]. The samples were searched against a database containing the total theoretical proteome of A. nidulans downloaded from UniProt (10,557 sequences) [58] supplemented with common contaminants such as keratins, trypsin, and BSA. In addition, reversed sequences of all protein entries were concatenated to the database to estimate the false discovery rate (FDR). Protein N-terminal acetylation, oxidation of methionine, conversion of glutamine to pyroglutamic acid, and deamidation of asparagine and glutamine were used as variable modifications, while carbamidomethylation of cysteine residues was used as a fixed modification. Trypsin was used as digestion enzyme and two missed cleavages were allowed. The ‘match between runs’ feature of MaxQuant was enabled with default parameters, to transfer identifications between samples based on accurate mass and retention time [57]. This was done to increase the number of identified peptides and set so that transfer of peptides was only allowed between samples from the same substrate. All identifications were filtered to achieve a protein FDR of 1 % and two ratio counts were required for a valid protein quantification. To evaluate the active site N-terminal histidine residue in both AA13 enzymes methylated, we performed a separate search using Mascot [59] including a variable modification of histidine-methylation as well as semitryptic cleavage pattern to obtain N-terminal matches after signal peptide cleavage.

Post-processing was done using Perseus version 1.5.0.31. Proteins categorized as only identified by site and matches to reversed sequences or contaminants were removed. Furthermore, the proteins were filtered so that a valid quantification existed for at least two of the three replicates on at least one substrate. Intensities were log-transformed and missing values imputed based on a downshifted normal distribution. Hierarchical clustering and heat map generation were done with Euclidean distance measure and average linkage. For visualization of trending proteins, ANOVA (permutation-based FDR, p < 0.10) was used to filter out proteins that showed no significant change over time. Protein quantitative values were z-scored within substrate and imputed values removed prior to clustering and heat map generation. Secretion prediction was a combination of SignalP [60], Phobius [61] and WoLF PSort [62] where at least two prediction algorithms had to agree. Secreted CAZymes were annotated using dbCAN (http://csbl.bmb.uga.edu/dbCAN/) [63].

β-Cyclodextrin affinity chromatography purification of secreted proteome and identification of most abundant proteins

The supernatant of the culture grown on wheat starch for 5 days was supplemented with (NH4)2SO4 to 0.5 M final concentration and agitated until the salt was completely dissolved. Thereafter, the sample was filtrated (0.45 µm membrane filters; Merck, Millipore) and applied to an XK 16/20 column packed with a 20 ml β-cyclodextrin (β-CD) sepharose affinity resin and pre-equilibrated with 4 column volumes of 10 mM sodium acetate buffer, pH 5.5, with 500 mM NaCl at 1 ml/min [64]. After sample loading at 0.75 ml/min, the column was washed with eight column volumes of the above buffer at 1 ml/min. Bound proteins were eluted with four column volumes of 20 mM sodium acetate buffer, 7 mM β-CD, pH 5.5. The purification was carried out using an ÄKTA Avant chromatograph interfaced by UNICORN 5.0 control software (GE Healthcare).

The eluted sample was analyzed by SDS-PAGE using Novex NuPAGE® 4–12 % Bis-Tris Gels (Invitrogen, Carlsbad, CA) in the XCell SureLock Mini-Cell system (invitrogen), according to manufacturer’s instructions and visualized by staining with InstantBlue solution (Expedeon, Cambridgeshire, UK). Spots were excised from the visually most prominent protein bands from the stained gel, washed in 300 µl 40 % ethanol at 50 °C for 15 min and in 100 µl 100 % ACN at room temperature for 10 min. The proteins in the gel were reduced with 50 µl 10 mM DTT in 100 mM NH4HCO3 at 56 °C for 45 min. Cysteine alkylation was performed by adding 100 µl 55 mM IAA in 100 mM NH4HCO3 and incubating at room temperature for 30 min in the dark. The gels were washed with 100 µl 50 % ACN, then 100 µl 100 % ACN and dried. The proteins in the gel were digested with 10 µl 12.5 ng/µl sequencing-grade modified trypsin (Promega) dissolved in 10 mM NH4HCO3 and incubated at 4 °C for 45 min, followed by addition of 10 µl 10 mM NH4HCO3 and incubation at 37 °C for 16 h. Samples (1 or 2 µl) were spotted directly onto a MTP AnchorChip target plate (Bruker Daltonics, Bremen, Germany), allowed to dry, and overlaid with 1 µl 0.5 µg/µl α-cyano-4-hydroxycinnamic acid (CHCA) matrix in 90 % ACN, 0.1 % TFA. The MS analyses were performed using an Ultraflex II MALDI-TOF/TOF MS (Bruker Daltonics). The obtained mass spectra were processed with FlexAnalysis and BioTools software both provided by the instrument manufacturer. Combination of MS and MS/MS data was used as input for databases searching for the spectra from MALDI-TOF/TOF using an in-house-licensed Mascot search engine (Matrix Science, London, UK). Proteins were identified using NCBInr database. The following parameters were set for searching: allowed global modification, carbamidomethyl cysteine; variable modification, oxidation of methionine; missed cleavages, 1; peptide tolerance, 80 ppm; MS/MS tolerance, ± 0.5 Da. The protein identification was considered valid if it matched more than two peptides and the significance threshold for protein identifications was p < 0.05.

Identification of the consensus binding for the AmyR transcription factor

The AmyR consensus binding site sequence 5′-CGGN8CGG-3′ was identified by searching 1000 bp upstream of the start codon of the A. nidulans genes encoding AnLPMO13A and AnLPMO13B. The positions of the 5′ and 3′ ends of the AmyR binding site with putative translational start site as +1 were calculated.

Abbreviations

AA:

auxiliary activity enzyme

CAZyme:

carbohydrate-active enzyme

CBM:

carbohydrate-binding module

CDH:

cellobiose dehydrogenase

β-CD:

β-cyclodextrin

FGSC:

fungal genetic stock center

GH:

glycoside hydrolase

HA:

high amylose

LC–MS/MS:

liquid chromatography-tandem mass spectroscopy

LPMO:

lytic polysaccharide monooxygenase

PL:

pectin lyase

References

  1. Zeeman SC, Kossmann J, Smith AM. Starch: its metabolism, evolution, and biotechnological modification in plants. Annu Rev Plant Biol. 2010;61:209–34.

    Article  CAS  Google Scholar 

  2. Pérez S, Bertoft E. The molecular structures of starch components and their contribution to the architecture of starch granules: a comprehensive review. Starch/Stӓrke. 2010;62(8):389–420.

    Article  Google Scholar 

  3. Buléon A, Colonna P, Planchot V, Ball S. Starch granules: structure and biosynthesis. Int J Biol Macromol. 1998;23(2):85–112.

    Article  Google Scholar 

  4. Tester RF, Karkalas J, Qi X. Starch - composition, fine structure and architecture. J Cereal Sci. 2004;39(2):151–65.

    Article  CAS  Google Scholar 

  5. Sun H, Zhao P, Ge X, Xia Y, Hao Z, Liu J, Peng M. Recent advances in microbial raw starch degrading enzymes. Appl Biochem Biotechnol. 2010;160(4):988–1003.

    Article  CAS  Google Scholar 

  6. Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sørlie M, Eijsink VGH. An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science. 2010;330(6001):219–22.

    Article  CAS  Google Scholar 

  7. Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JCN, Johansen KS, Krogh KBRM, Jørgensen CI, Tovborg M, Anthonsen A, et al. Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci. 2011;108(37):15079–84.

    Article  CAS  Google Scholar 

  8. Vu VV, Beeson WT, Phillips CM, Cate JHD, Marletta MA. Determinants of regioselective hydroxylation in the fungal polysaccharide monooxygenases. J Am Chem Soc. 2014;136(2):562–5.

    Article  CAS  Google Scholar 

  9. Hemsworth GR, Davies GJ, Walton PH. Recent insights into copper-containing lytic polysaccharide mono-oxygenases. Curr Opin Struct Biol. 2013;23(5):660–8.

    Article  CAS  Google Scholar 

  10. Forsberg Z, Vaaje-Kolstad G, Westereng B, Bunæs AC, Stenstrøm Y, Mackenzie A, Sørlie M, Horn SJ, Eijsink VGH. Cleavage of cellulose by a CBM33 protein. Protein Sci. 2011;20(9):1479–83.

    Article  CAS  Google Scholar 

  11. Johansen KS. Discovery and industrial applications of lytic polysaccharide mono-oxygenases. Biochem Soc Trans. 2016;44(1):143–9.

    Article  CAS  Google Scholar 

  12. Isaksen T, Westereng B, Aachmann FL, Agger JW, Kracher D, Kittl R, Ludwig R, Haltrich D, Eijsink VGH, Horn SJ. A C4-oxidizing lytic polysaccharide monooxygenase cleaving both cellulose and cello-oligosaccharides. J Biol Chem. 2014;289(5):2632–42.

    Article  CAS  Google Scholar 

  13. Agger JW, Isaksen T, Varnai A, Vidal-Melgosa S, Willats WGT, Ludwig R, Horn SJ, Eijsink VGH, Westereng B. Discovery of LPMO activity on hemicelluloses shows the importance of oxidative processes in plant cell wall degradation. Proc Natl Acad Sci. 2014;111(17):6287–92.

    Article  CAS  Google Scholar 

  14. Frommhagen M, Sforza S, Westphal AH, Visser J, Hinz SWA, Koetsier MJ, van Berkel WJH, Gruppen H, Kabel MA. Discovery of the combined oxidative cleavage of plant xylan and cellulose by a new fungal polysaccharide monooxygenase. Biotechnol Biofuels. 2015;8(1):101.

    Article  Google Scholar 

  15. Bennati-Granier C, Garajova S, Champion C, Grisel S, Haon M, Zhou S, Fanuel M, Ropartz D, Rogniaux H, Gimbert I, Record E, Berrin JG. Substrate specificity and regioselectivity of fungal AA9 lytic polysaccharide monooxygenases secreted by Podospora anserina. Biotechnol Biofuels. 2015;8(1):90.

    Article  Google Scholar 

  16. Vu VV, Beeson WT, Span EA, Farquhar ER, Marletta MA. A family of starch-active polysaccharide monooxygenases. Proc Natl Acad Sci. 2014;111(38):13822–7.

    Article  CAS  Google Scholar 

  17. Lo Leggio L, Simmons TJ, Poulsen JCN, Frandsen KEH, Hemsworth GR, Stringer MA, von Freiesleben P, Tovborg M, Johansen KS, De Maria L, et al. Structure and boosting activity of a starch-degrading lytic polysaccharide monooxygenase. Nat Commun. 2015;6:5961.

    Article  CAS  Google Scholar 

  18. Vu VV, Marletta MA. Starch-degrading polysaccharide monooxygenases. Cell Mol Life Sci. 2016;73(14):2809–19.

    Article  CAS  Google Scholar 

  19. Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol Biofuels. 2013;6(1):41.

    Article  CAS  Google Scholar 

  20. Navarro D, Rosso MN, Haon M, Olivé C, Bonnin E, Lesage-Meessen L, Chevret D, Coutinho PM, Henrissat B, Berrin JG. Fast solubilization of recalcitrant cellulosic biomass by the basidiomycete fungus Laetisaria arvalisinvolves successive secretion of oxidative and hydrolytic enzymes. Biotechnol Biofuels. 2014;7(1):143.

    Article  Google Scholar 

  21. Saykhedkar S, Ray A, Ayoubi-Canaan P, Hartson SD, Prade R, Mort AJ. A time course analysis of the extracellular proteome of Aspergillus nidulans growing on sorghum stover. Biotechnol Biofuels. 2012;5(1):52.

    Article  CAS  Google Scholar 

  22. Wortman JR, Gilsenan JM, Joardar V, Deegan J, Clutterbuck J, Andersen MR, Archer D, Bencina M, Braus G, Coutinho P, et al. The 2008 update of the Aspergillus nidulans genome annotation: a community effort. Fungal Genet Biol. 2009;46(1):S2–13.

    Article  CAS  Google Scholar 

  23. Hata Y, Ishida H, Kojima Y, Ichikawa E, Kawato A, Suginami K, IMayasu S. Comparison of two glucoamylases produced by Aspergillus oryzae in solid-state culture (koji) and in submerged culture. J Ferment Bioeng. 1997;84(6):532–7.

    Article  CAS  Google Scholar 

  24. Garajova S, Mathieu Y, Beccia MR, Bennati-Granier C, Biaso F, Fanuel M, Ropartz D, Guigliarelli B, Record E, Rogniaux H, Henrissat B, Berrin JG. Single-domain flavoenzymes trigger lytic polysaccharide monooxygenases for oxidative degradation of cellulose. Sci Rep. 2016;6:28276.

    Article  CAS  Google Scholar 

  25. Christiansen C. Abou Hachem M, Janecek S, Viksø-Nielsen A, Blennow A, Svensson B. The carbohydrate-binding module family 20—diversity, structure, and function. FEBS J. 2009;276(18):5006–29.

    Article  CAS  Google Scholar 

  26. Ho DP, Ngo HH, Guo W. A mini review on renewable sources for biofuel. Bioresour Technol. 2014;169:742–9.

    Article  CAS  Google Scholar 

  27. Yue D, You F, Snyder SW. Biomass-to-bioenergy and biofuel supply chain optimization: overview, key issues and challenges. Comput Chem Eng. 2014;66:36–56.

    Article  CAS  Google Scholar 

  28. Langston JA, Shaghasi T, Abbate E, Xu F, Vlasenko E, Sweeney MD. Oxidoreductive cellulose depolymerization by the enzymes cellobiose dehydrogenase and glycoside hydrolase 61. Appl Environ Microbiol. 2011;77(19):7007–15.

    Article  CAS  Google Scholar 

  29. Hori C, Ishida T, Igarashi K, Samejima M, Suzuki H, Master E, Ferreira P, Ruiz-Duenas FJ, Held B, Canessa P, et al. Analysis of the Phlebiopsis gigantea genome, transcriptome and secretome provides insight into its pioneer colonization strategies of wood. PLoS Genet. 2014;10(12):e1004759.

    Article  Google Scholar 

  30. Ring SG, Gee JM, Whittam M, Orford P, Johnson IT. Resistant starch: its chemical form in foodstuffs and effect on digestibility in vitro. Food Chem. 1988;28(2):97–109.

    Article  CAS  Google Scholar 

  31. Smith AM. The biosynthesis of starch granules. Biomacromolecules. 2001;2(2):335–41.

    Article  CAS  Google Scholar 

  32. Ratnayake WS, Hoover R, Warkentin T. Pea starch: composition, structure and properties—a review. Starch/Stӓrke. 2002;54(6):217–34.

    Article  CAS  Google Scholar 

  33. Planchot V, Colonna P, Gallant DJJ, Bouchet B. Extensive degradation of native starch granules by α-amylase from Aspergillus fumigatus. J Cereal Sci. 1995;21(2):163–71.

    Article  CAS  Google Scholar 

  34. Salgado P, Freire JPB, Ferreira RB, Teixeira A, Bento O, Abreu MC, Toullec R, Lalles JP. Immunodetection of legume proteins resistant to small intestinal digestion in weaned piglets. J Sci Food Agric. 2003;83(15):1571–80.

    Article  CAS  Google Scholar 

  35. Palmer R, Cornuault V, Marcus SE, Knox JP, Shewry PR, Tosi P. Comparative in situ analyses of cell wall matrix polysaccharide dynamics in developing rice and wheat grain. Planta. 2014;241(3):669–85.

    Article  Google Scholar 

  36. Knudsen KEB. Fiber and nonstarch polysaccharide content and variation in common crops used in broiler diets. Poult Sci. 2014;93(9):2380–93.

    Article  CAS  Google Scholar 

  37. Andersen MR, Vongsangnak W, Panagiotou G, Salazar MP, Lehmann L, Nielsen J. A trispecies Aspergillus microarray: comparative transcriptomics of three Aspergillus species. Proc Natl Acad Sci. 2008;105(11):4387–92.

    Article  CAS  Google Scholar 

  38. Le DM, Fojan P, Azem E, Pettersson D, Pedersen NR. Visualization of the anticaging effect of Ronozyme WX xylanase on wheat substrates. Cereal Chem. 2013;90(5):439–44.

    Article  Google Scholar 

  39. Choct M. Feed non-starch polysaccharides: chemical structures and nutritional significance. Feed Milling Int. 1997;191:13–26.

    Google Scholar 

  40. Coutinho PM, Andersen MR, Kolenova K, VanKuyk PA, Benoit I, Gruben BS, Trejo-Aguilar B, Visser H, Van Solingen P, Pakula T. Post-genomic insights into the plant polysaccharide degradation potential of Aspergillus nidulans and comparison to Aspergillus niger and Aspergillus oryzae. Fungal Genet Biol. 2009;46(1):S161–9.

    Article  CAS  Google Scholar 

  41. Bauer S, Vasu P, Persson S, Mort AJ, Somerville CR. Development and application of a suite of polysaccharide-degrading enzymes for analyzing plant cell walls. Proc Natl Acad Sci. 2006;103(30):11417–22.

    Article  CAS  Google Scholar 

  42. Matsubara T, Ammar YB, Anindyawati T, Yamamoto S, Ito K, Iizuka M, Minamiura N. Degradation of raw starch granules by α-amylase purified from culture of Aspergillus awamori KT-11. J Biochem Mol Biol. 2004;37(4):422–8.

    Article  CAS  Google Scholar 

  43. van der Kaaij RM, Janecek S, van der Maarel MJEC, Dijkhuizen L. Phylogenetic and biochemical characterization of a novel cluster of intracellular fungal α-amylase enzymes. Microbiology. 2007;153(12):4003–15.

    Article  Google Scholar 

  44. Forsberg Z, Mackenzie AK, Sorlie M, Rohr ÅK, Helland R, Arvai AS, Vaaje-Kolstad G, Eijsink VGH. Structural and functional characterization of a conserved pair of bacterial cellulose-oxidizing lytic polysaccharide monooxygenases. Proc Natl Acad Sci. 2014;111(23):8446–51.

    Article  CAS  Google Scholar 

  45. Hemsworth GR, Henrissat B, Davies GJ, Walton PH. Discovery and characterization of a new family of lytic polysaccharide monooxygenases. Nat Chem Biol. 2014;10(2):122–6.

    Article  CAS  Google Scholar 

  46. Tan T, Kracher D, Gandini R, Sygmund C, Kittl R, Haltrich D, Hällberg BM, Ludwig R, Divne C. Structural basis for cellobiose dehydrogenase action during oxidative cellulose degradation. Nat Commun. 2015;6:7542.

    Article  Google Scholar 

  47. Henriksson G, Johansson G, Pettersson G. A critical review of cellobiose dehydrogenases. J Biotechnol. 2000;78(2):93–113.

    Article  CAS  Google Scholar 

  48. Kracher D, Scheiblbrandner S, Felice AKG, Breslmayr E, Preims M, Haltrich D, Eijsink VGH, Ludwig R. Extracellular electron transfer systems fuel cellulose oxidative degradation. Science. 2016;352(6289):1098–101.

    Article  CAS  Google Scholar 

  49. Couturier M, Mathieu Y, Li A, Navarro D, Drula E, Haon M, Grisel S, Ludwig R, Berrin JG. Characterization of a new aryl-alcohol oxidase secreted by the phytopathogenic fungus Ustilago maydis. Appl Microbiol Biotechnol. 2016;100(2):697–706.

    Article  CAS  Google Scholar 

  50. Calera JA, Sánchez-Weatherby J, López-Medrano R, Leal F. Distinctive properties of the catalase B of Aspergillus nidulans. FEBS Lett. 2000;475(2):117–20.

    Article  CAS  Google Scholar 

  51. Kittl R, Kracher D, Burgstaller D, Haltrich D, Ludwig R. Production of four Neurospora crassa lytic polysaccharide monooxygenases in Pichia pastoris monitored by a fluorimetric assay. Biotechnol Biofuels. 2012;5(1):79.

    Article  CAS  Google Scholar 

  52. Thon M, Al-Abdallah Q, Hortschansky P, Brakhage AA. The Thioredoxin system of the filamentous fungus Aspergillus nidulans: impact on development and oxidative stress response. J Biol Chem. 2007;282(37):27259–69.

    Article  Google Scholar 

  53. Scott BR, Huang HZ, Frickman J, Halvorsen R, Johansen KS. Catalase improves saccharification of lignocellulose by reducing lytic polysaccharide monooxygenase-associated enzyme inactivation. Biotechnol Lett. 2016;38:425–34.

    Article  CAS  Google Scholar 

  54. Nakamura T, Maeda Y, Tanoue N, Makita T, Kato M, Kobayashi T. Expression profile of amylolytic genes in Aspergillus nidulans. Biosci Biotechnol Biochem. 2006;70(10):2363–70.

    Article  CAS  Google Scholar 

  55. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26(12):1367–72.

    Article  CAS  Google Scholar 

  56. Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M. Andromeda: a peptide search engine integrated into the MaxQuant environment. J Proteome Res. 2011;10(4):1794–805.

    Article  CAS  Google Scholar 

  57. Cox J, Hein MY, Luber CA, Paron I. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol Cell Proteomics. 2014;13(9):2513–26.

    Article  CAS  Google Scholar 

  58. Consortium TU. UniProt: a hub for protein information. Nucleic Acids Res. 2014;43(D):204–12.

    Google Scholar 

  59. Perkins DN, Pappin DJC, Creasy DM, Cottrell JS. Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. 1999;20(18):3551–67.

    Article  CAS  Google Scholar 

  60. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.

    Article  CAS  Google Scholar 

  61. Käll L, Krogh A, Sonnhammer ELL. Advantages of combined transmembrane topology and signal peptide prediction-the Phobius web server. Nucleic Acids Res. 2007;35(W):429–32.

    Article  Google Scholar 

  62. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35(W):585–7.

    Article  Google Scholar 

  63. Yin Y, Mao X, Yang J, Chen X, Mao F, Xu Y. DbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40(W):445–51.

    Article  Google Scholar 

  64. Vester-Christensen MB. Abou Hachem M, Naested H, Svensson B. Secretory expression of functional barley limit dextrinase by Pichia pastoris using high cell-density fermentation. Protein Expr Purif. 2010;69(1):112–9.

    Article  CAS  Google Scholar 

Download references

Authors’ contributions

LN, GVK, and MAH designed all the experiments. BS helped in the design of the starch growth experiments. LN performed all the experiments, except secretome analysis. MØA performed secretome analysis. LN and MØA analyzed secretomics data. LN, MAH, and GVK wrote the first manuscript draft. BS helped in revising the manuscript. All authors participated in finalizing the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of supporting data

All supporting data related to this work are included as additional supporting files.

Consent for publication

All the authors consent to publish this research work.

Funding

This work was supported by a Grant (NNF12OC0000769) to MAH from the Novo Nordisk foundation in the area “Biotechnology-based Synthesis and Production”. GVK was funded by Grant 214138 provided by the Norwegian Research Council.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Gustav Vaaje-Kolstad or Maher Abou Hachem.

Additional information

Laura Nekiunaite and Magnus Ø. Arntzen contributed equally to this work

Additional files

13068_2016_604_MOESM1_ESM.xlsx

Additional file 1: Table S1. Proteins identified in the secretome of Aspergillus nidulans. The table shows the MaxQuant output of identified proteins and the quantitative LFQ intensities associated with all replicates and conditions. Prediction of secretion was done as a combination of SignalP, Phobius and WoLF PSORT. CAZy annotation was done with dbCAN.

13068_2016_604_MOESM2_ESM.pdf

Additional file 2: Figure S1. Heat map comparison of expression patterns of 312 secreted proteins detected after 3–5 days growth of A. nidulans on minimal media supplemented with: high-amylose (HA) maize, wheat, or pea starch. The colors in the heat map indicate the label-free quantification (LFQ) intensity reported by MaxQuant ranging from 2 × 106 (light green) to 4 × 1011 (light red). Missing values were imputed from a normal distribution located at the quantification limit. Highlighted in green are AA13s, in purple—AA9s and AA11s, in blue—AA3s and AA7s, in red—hydrolytic enzymes associated with starch degradation, in gray—proteases, and in brown—cell wall-degrading enzymes. Secretion prediction is a combination of SignalP, Phobius, and WolfPSort where at least two algorithms had to agree.

13068_2016_604_MOESM3_ESM.pdf

Additional file 3: Figure S2. A trending heat map of proteins predicted to be secreted during growth of Aspergillus nidulans on high-amylose (HA) maize, pea or wheat starch. Two main clusters are indicated in the graph, proteins with decreasing levels over time (I) and proteins with increasing levels over time (II). The intensities are z-score normalized per row, separately for each substrate, to emphasize the trend of protein amount during growth independent of total level and substrate-level differences. Proteins showing no change over time (ANOVA, permutation-based FDR >0.1) were removed from the plot. Missing values (protein level less than detection limit of mass spectrometer) are colored white. Highlighted in green are AA13s, in purple—AA9s and AA11s, in blue—AA3s and AA7s, in red—hydrolytic enzymes associated with starch degradation, in gray—proteases, and in brown—cell wall-degrading enzymes.

13068_2016_604_MOESM4_ESM.docx

Additional file 4: Table S2. Comparison of detected CAZymes in the secretome of Aspergillus nidulans during growth on wheat, high-amylose maize and pea starch at day 3, 4 and 5.

13068_2016_604_MOESM5_ESM.docx

Additional file 5: Tables S3–S5. Table S3. Top 20 detected proteins in the secretome of Aspergillus nidulans during growth on wheat starch at day 3, 4 and 5. Table S4. Top 20 detected proteins in the secretome of Aspergillus nidulans during growth on high-amylose maize starch at day 3, 4 and 5. Table S5. Top 20 detected proteins in the secretome of Aspergillus nidulans during growth on pea starch at day 3, 4 and 5.

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nekiunaite, L., Arntzen, M.Ø., Svensson, B. et al. Lytic polysaccharide monooxygenases and other oxidative enzymes are abundantly secreted by Aspergillus nidulans grown on different starches. Biotechnol Biofuels 9, 187 (2016). https://doi.org/10.1186/s13068-016-0604-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-016-0604-0

Keywords