Skip to main content
Download PDF

Boosting LPMO-driven lignocellulose degradation by polyphenol oxidase-activated lignin building blocks

Biotechnology for Biofuels volume 10, Article number: 121 (2017) Cite this article

Abstract

Background

Many fungi boost the deconstruction of lignocellulosic plant biomass via oxidation using lytic polysaccharide monooxygenases (LPMOs). The application of LPMOs is expected to contribute to ecologically friendly conversion of biomass into fuels and chemicals. Moreover, applications of LPMO-modified cellulose-based products may be envisaged within the food or material industry.

Results

Here, we show an up to 75-fold improvement in LPMO-driven cellulose degradation using polyphenol oxidase-activated lignin building blocks. This concerted enzymatic process involves the initial conversion of monophenols into diphenols by the polyphenol oxidase MtPPO7 from Myceliophthora thermophila C1 and the subsequent oxidation of cellulose by MtLPMO9B. Interestingly, MtPPO7 shows preference towards lignin-derived methoxylated monophenols. Sequence analysis of genomes of 336 Ascomycota and 208 Basidiomycota reveals a high correlation between MtPPO7 and AA9 LPMO genes.

Conclusions

The activity towards methoxylated phenolic compounds distinguishes MtPPO7 from well-known PPOs, such as tyrosinases, and ensures that MtPPO7 is an excellent redox partner of LPMOs. The correlation between MtPPO7 and AA9 LPMO genes is indicative for the importance of the coupled action of different monooxygenases in the concerted degradation of lignocellulosic biomass. These results will contribute to a better understanding in both lignin deconstruction and enzymatic lignocellulose oxidation and potentially improve the exploration of eco-friendly routes for biomass utilization in a circular economy.

Background

Fungal carbohydrate converting enzymes are considered important for eco-friendly application in plant biomass degradation. The resulting carbohydrates are sources for the production of biochemicals or biofuels and new enzymatically modified cellulose-based products are envisaged for the future.

Next to carbohydrates, phenolic compounds are also major components of lignocellulosic plant biomass. Phenolic compounds are present either as free molecules or in conjugated form as part of lignin or bound to carbohydrates. Lignin is one of the major constituents of (secondary) plant cell walls, together with the polysaccharides cellulose and hemicellulose. Lignin is composed of the three aromatic monolignol units: coniferyl, sinapyl, and para-coumaryl alcohol. Cellulose consists of β-(1→4)-linked glucosyl chains that interact with each other via hydrogen bonds and van der Waals forces, which results in the formation of crystalline cellulose fibrils. Unlike cellulose, hemicellulose is a heteropolymer varying in its monosaccharide composition and linkages between the monosaccharides. Examples are xylan, mannan, or β-(1→3, 1→4)-linked β-glucan. Hemicellulose interacts with lignin through ester and ether linkages, thereby forming a network that embeds the cellulose microfibrils [1, 2].

Recent studies focused on the function of lytic polysaccharide monooxygenases (LPMOs) have confirmed that these enzymes drive the oxidative degradation of cellulose, and they are considered important for the enzymatic degradation of plant biomass [3]. LPMOs are classified based on their sequence in the Carbohydrate Active enzyme (CAZy; [4]) database as auxiliary activity (AA) families AA9, AA10, AA11, and AA13. In brief, LPMOs have been reported to oxidize β-(1→4)-linked glucan chains at either the C1- or C4-carbon position or both, resulting in the cleavage of glucan chains [3, 5, 6]. Other LPMOs of these AA families have been described to oxidize the (1→4)-linkage of chitin, xylan, hemicellulose, such as xyloglucan and glucomannan, soluble cellodextrins, and starch [7,8,9,10,11].

In order to oxidize polysaccharides, LPMOs demand electrons to activate molecular oxygen in their copper-containing active site [3, 5, 6]. The electrons can be donated by reducing agents, like low molecular weight compounds (ascorbic acid, gallic acid) or the macromolecule lignin [12,13,14]. Other ways of providing LPMOs with electrons have been reported, such as flavocytochrome-dependent cellobiose dehydrogenases (CDHs), light-induced pigments, light-driven chemical oxidation of water, or diphenol-regenerating GMC (glucose-methanol-choline-oxidase/dehydrogenase)-oxidoreductases [5, 12, 15, 16]. The mechanistic understanding of electron donation systems is highly relevant in order to enable optimization of LPMO activity and, thereby, plant biomass degradation.

Phenolic compounds, including small molecular weight compounds that serve as lignin building blocks and lignin, are intrinsically present in plant biomass and are natural electron donors for LPMO activity. However, monophenols are not optimal electron donors for LPMO activity, because of their relatively high redox potential [12]. Compounds with a 1,2-benzenediol or a 1,2,3-benzenetriol moiety have, compared to monophenols, a lower redox potential [13]. Their low redox potential enables them to reduce the copper ion in the active site of LPMOs and enhance the LPMO activity [12].

The enzymatic oxidation of phenolic compounds is a well-known reaction in many natural environments. For example, it causes browning of food products, contributes to taste in tea fermentation, and plays a role in plant biomass decomposition [17,18,19]. These oxidative systems involve the activity of laccases (EC.1.10.3.2), peroxidases (EC.1.11.1), tyrosinases (EC.1.14.18.1), and catechol oxidases (EC 1.10.3.1) [18,19,20]. In particular tyrosinases, also often referred to as polyphenol oxidases, are of interest due to their ability to hydroxylate phenolic compounds [21]. This so-called monophenolase activity typically involves the ortho-hydroxylation of monophenols into ortho-diphenols, compounds that comprise a 1,2-benzenediol moiety. Tyrosinases also exhibit diphenolase activity, which is characterized by the oxidation of these ortho-diphenols into ortho-quinones [22]. With respect to biomass degradation, the diphenolase activity of polyphenol oxidases is not conductive to polysaccharides oxidation since ortho-quinones cannot be utilized by LPMOs [12]. Based on their monophenolase activity, tyrosinases have been shown to use non-methoxylated monophenols as substrates, rather than using methoxylated monophenols, that are the predominant structural units of lignin [23].

In this study, we investigated if polyphenol oxidases, in particular tyrosinases, can enhance LPMO activity. Therefore, the tyrosinase MtPPO7 from Myceliophthora thermophila C1 was used, which showed activity towards phenolic compounds and was obtained from DuPont Industrial Biosciences. In addition, we used the commercially available tyrosinase AbPPO from the white button mushroom Agaricus bisporus. Both tyrosinases were each incubated with the previously described MtLPMO9B in the presence of various plant phenolic compounds and cellulose [13, 24]. We found that MtPPO7 is highly active towards methoxylated monophenols, including monomeric lignin building blocks. This activity can strongly boost the release of oxidized gluco-oligosaccharides as catalyzed by MtLPMO9B and thereby the degradation of cellulose. In addition, we found a strong correlation between genes encoding MtPPO like proteins and AA9 LPMOs in fungal genomes of 336 Ascomycota and 208 Basidiomycota.

Results

MtPPO7 addition improves cellulose oxidation by MtLPMO9B

Our previous results have shown that reducing agents with a 1,2-benzenediol or 1,2,3-benzenetriol moiety gave the highest release of non-oxidized and C1-oxidized gluco-oligosaccharides from regenerated amorphous cellulose (RAC) incubated with three MtLPMOs compared to the incubation with compounds comprising only a single hydroxyl group [13]. Hence, in the current research we hypothesized that MtLPMO activity can benefit from enzymes which have hydroxylating capacity, such as polyphenol oxidases (Fig. 1) [21]. We choose MtPPO7, which originates like the MtLPMO9B employed here, from the thermophilic filamentous fungus Myceliophthora thermophila C1. As a reference, the well-studied tyrosinase AbPPO from the edible mushroom Agaricus bisporus was used [19, 25,26,27,28].

Fig. 1
figure 1

Schematic presentation of the concerted action of polyphenol oxidases and LPMOs. Monophenols with a 1-hydroxy, 2-methoxy moiety are hydroxylated by MtPPO7 from M. thermophila C1. The resulting methoxylated catechols are excellent electron donors for MtLPMO9B activity. In comparison to MtPPO7, the mushroom tyrosinase AbPPO is able to convert non-methoxylated monophenols into compounds comprising a 1,2-benzenediol moiety. The released compounds comprising a 1,2-benzenediol moiety are able to donate electrons for MtLPMO9B activity. However, AbPPO exhibits a strong diphenolase activity which reduces the available amount of these compounds comprising a 1,2-benzenediol moiety for MtLPMO9B due to further oxidation of these compounds into ortho-quinones (indicated by dashed arrow)

Full size image

The activity of MtLPMO9B towards RAC, with and without addition of MtPPO7, was determined using 21 phenolic compounds as potential electron donors. Many of these phenolic compounds are methoxylated as these are the predominant structural units of lignin. Results are presented in Table 1 and Fig. 2. The 21 different phenolic compounds used were classified into three groups and further divided into subgroups, as previously described [13]. In short, group I represents compounds comprising one hydroxyl group, group II compounds comprising a 1,2-benzenediol moiety, and group III compounds comprising a 1,2,3-benzenetriol moiety (Fig. 2; Table 1).

Table 1 RAC incubated with MtLPMO9B only and with AbPPO or MtPPO7 addition in the presence of 21 reducing agents
Full size table
Fig. 2
figure 2

RAC incubated with MtLPMO9B with or without MtPPO7 or AbPPO addition in the presence of 21 reducing agents. The numbers are total sums of the integrated peak areas of released non-oxidized (shaded red, blue, and yellow) and C1-oxidized (red, blue, and yellow) gluco-oligosaccharides after incubation of regenerated amorphous cellulose (RAC; 1.5 mg mL−1) with MtLPMO9B only (red bars) (5.0 μg mL−1), MtLPMO9B together with AbPPO (blue bars, 2.5 µL mL−1), or MtLPMO9B together with MtPPO7 (yellow bars) (5.0 μg mL−1) based on HPAEC. The reducing agents (2 mM) are numbered (X-axis) and specified in Table 1. Vertical dotted lines separate the reducing agents into groups (I, II, II) and subgroups (a, b, c). Group I, monophenols (Ia), compounds with a 1-hydroxy,2-methoxy moiety (Ib) or a 1-hydroxy-2,6-dimethoxy moiety (Ic). Group II, compounds with 1,2-benzenediols moiety (IIa) and compounds with a 1,2-dihydroxy-3-methoxy moiety (IIb). Group III, reducing agents with a 1,2,3-benzenetriol moiety (III). All incubations were performed in duplicate and the standard deviations are represented by error bars. See “Methods” for details

Full size image

As expected, from our previous research [13], monophenols (group Ia) are less efficient electron donors for MtLPMO9B, with the exception of 4-hydroxybenzoic acid (no. 1). The low MtLPMO9B activity towards RAC resulted in a very small amount of released non-oxidized and C1-oxidized gluco-oligosaccharides (Fig. 2; Table 1; Additional file 1: Figure S1). The addition of MtPPO7 to RAC incubated with MtLPMO9B did not increase the electron donor capacities of this group Ia compounds, except for para-coumaric acid (no. 3). In contrast, the addition of MtPPO7 to RAC incubated with MtLPMO9B in the presence of methoxylated monophenols (group Ib, compounds with a 1-hydroxy, 2-methoxy moiety) led in all cases to a high increase (up to ±75 times, no. 5) in the release of non-oxidized and C1-oxidized gluco-oligosaccharides (Fig. 2; Additional file 1: Figure S1). Similarly, the oxidative cleavage of RAC incubated with MtLPMO9B increased, when MtPPO7 was added to group Ic compounds, which comprise a 1-hydroxy-2,6-dimethoxy moiety. For group IIa, IIb, and III phenolic compounds, the addition of MtPPO7 to incubations of RAC with MtLPMO9B led to a moderate increase of non-oxidized and C1-oxidized gluco-oligosaccharides (4–39%), with one exception (no. 19) (Fig. 2; Table 1).

Incubations performed with AbPPO showed deviating results from corresponding experiments with MtPPO7. First, unlike MtPPO7, addition of AbPPO to non-methoxylated monophenols (group Ia) resulted in increased levels of non-oxidized and C1-oxidized oligosaccharides during the incubation of RAC with MtLPMO9B. This different activity was expected on the basis of the known hydroxylating capacity of AbPPO towards non-methoxylated monophenolic compounds. Especially, the addition of AbPPO in the presence of para-coumaric acid (no. 3) increased MtLPMO9B-catalyzed degradation of RAC compared to the incubation without AbPPO (Fig. 2; Table 1). In contrast to MtPPO7, AbPPO addition to group Ib, Ic, II, and III compounds improved the MtLPMO9B-catalyzed RAC degradation relatively moderate (5–361%, Table 1) or even decreased the RAC degradation by 1–44% (Fig. 2; Table 1).

MtLPMO9B-mediated cellulose oxidation in the presence of MtPPO7 using time-dependent incubations

The ability of MtPPO7 to convert methoxylated phenolic compounds and to increase the oxidative degradation of RAC by MtLPMO9B was further investigated using time-dependent incubations (24 h). These incubations were performed in the presence of two methoxylated monophenols (guaiacol and ferulic acid, group Ib), as well as in the presence of a non-methoxylated ortho-diphenol (3-methylcatechol, group IIa) (Fig. 3; Additional file 2: Figure S2). After 4 h of incubation, a steady MtPPO7-induced oxidation of guaiacol (no. 9) was observed (Additional file 3: Figure S3). At the same time, the amounts of released non-oxidized and C1-oxidized gluco-oligosaccharides by MtLPMO9B from RAC also increased steadily (Fig. 3c). A similar trend for the release of non-oxidized and C1-oxidized gluco-oligosaccharides was shown for the incubation with ferulic acid (no. 8, Additional file 2: Figure S2). MtPPO7 showed a relatively low efficiency towards 3-methylcatechol (no. 17) compared to guaiacol (Table 1; Additional file 3: Figure S3). As reported above, no significant increase of MtLPMO9B-released non-oxidized and C1-oxidized gluco-oligosaccharides from RAC was observed when comparing incubations with and without MtPPO7 (Fig. 3d). In presence of MtPPO7 though, the initial rate (0–6 h) of released non-oxidized and C1-oxidized gluco-oligosaccharides was lower compared to the incubation without MtPPO7 (Fig. 3d).

Fig. 3
figure 3

Activity of MtLPMO9B towards amorphous cellulose in the presence and absence of MtPPO7. HPAEC elution pattern of regenerated amorphous cellulose (RAC; 1.5 mg mL−1) incubated with MtLPMO9B (5.0 μg mL−1) with (yellow) and without (red) MtPPO7 (5.0 μg mL−1) addition in the presence of a guaiacol (no. 9 specified in Table 1, 2 mM) and b 3-methylcatechol (no. 17 specified in Table 1, 2 mM) after 24 h. The incubation of RAC with MtLPMO9B results in the formation of non-oxidized gluco-oligosaccharides (GlcOSn) and C1-oxidized gluco-oligosaccharides (GlcOS #n ). Based on HPAEC, integrated peak areas are shown as the total sum of released non-oxidized (shaded red and yellow) and C1-oxidized (red and yellow) gluco-oligosaccharides after incubation of RAC (1.5 mg mL−1) with MtLPMO9B only (red bars, 5.0 μg mL−1) and MtLPMO9B together with MtPPO7 (yellow bars, 5.0 μg mL−1) in the presence of c guaiacol and d 3-methylcatechol. All incubations were performed in duplicate, and the standard deviations are presented as error bars. See “Methods” for data analysis and details

Full size image

Activity of MtPPO7 towards methoxylated phenolic compounds

The most striking observation from the experiments described above is the conversion of methoxylated phenolic compounds by MtPPO7 into products, which enhance the oxidative degradation of cellulose by MtLPMO9B. The fate of these phenolic compounds upon MtPPO7 incubation was further studied in the absence of MtLPMO9B using UHPLC-UV-MSn.

Based on UV measurements, non-methoxylated phenolic compounds of group Ia, group IIa, and group III remained constant in concentration during MtPPO7 incubation, showing that MtPPO7 has a low activity towards compounds comprising a 1,2-benzenediol moiety (Fig. 4; Additional file 4: Table S1). Different from MtPPO7, in previous research AbPPO has been reported to be highly active towards these types compounds [19, 24, 25].

Fig. 4
figure 4

UV response areas of 21 reducing agents incubated in the presence and absence of MtPPO7. The total sum is shown of integrated peak areas of 21 reducing agents (2 mM) with (yellow bar) and without (red bar) addition of MtPPO7 (5.0 μg mL−1). The reducing agents are numbered (bold on the left) and specified in Table 1. Samples were incubated in a 50 mM potassium phosphate buffer (pH = 6.0) containing 2.5 µM copper(II)-chloride for 24 h at 50 °C and measured by UHPLC-UV (280 nm). Bold numbers on the right reducing agents conversion (%) by MtPPO7, which is based on the difference of integrated peak areas (UV 280 nm) of the reducing agents incubated with MtPPO7 compared to the incubation of reducing agents only. All incubations were performed in duplicate. See “Methods” for details

Full size image

MtPPO7 was active towards all eight methoxylated monophenols (group Ib) tested. Based on UV, the concentrations of group Ib compounds decreased during MtPPO7 incubation between 28 and 67% compared to the same compounds incubated without MtPPO7 (Fig. 4). Three types of MtPPO7 reactions were observed based on mass differences between substrate and products formed (Additional file 4: Table S1): (A) hydroxylation (m/z +16; e.g., no. 5 and 7), (B) decarboxylation (m/z −44, e.g., no. 8), and (C) demethylation (m/z −14, e.g., no. 11). Especially the hydroxylation of methoxylated phenolic compounds is a key reaction, since the products formed comprise a second hydroxyl group. The formed ortho-diphenols are known to be efficient electron donors for LPMOs. Decarboxylation and demethylation occurred either in the presence of MtPPO7 or as a result of polymerization reactions in both the presence and absence of MtPPO7. Several group Ib compounds spontaneously formed dimers in the absence of MtPPO7 (no. 8 or no. 11). These dimers were almost absent when group Ib compounds were incubated with MtPPO7 (Additional file 4: Table S1). Similar to group Ib compounds, the concentrations of compounds comprising a 1-hydroxyl-2,6-dimethoxy moiety (group Ic) incubated with MtPPO7 decreased, based on UV, between 6% and 64% compared to the incubation without MtPPO7 (Fig. 4; Additional file 4: Table S1). Again, masses indicating decarboxylation (−44 Da) and demethylation (−14 Da) reactions were formed. MtPPO7 was also able to convert group IIb compounds (up to 87% substrate conversion, no. 20) that comprise a 1,2-dihydroxy-3-methoxy moiety (Fig. 4). The reactions observed were decarboxylation (m/z −44, e.g., no. 19) and dimerization (e.g., no. 20) of group IIb compounds based on masses formed (Additional file 4: Table S1).

In general, most phenolic compounds (such as no. 17 and no. 21) that were incubated with MtPPO7 formed insoluble complexes, which likely result from polymerization reactions of ortho-quinones formed. These insoluble complexes were not determined by UHPLC-UV-MSn. The complexation of ortho-quinones resulted, obviously, in a decrease of the UV signal of the substrates (Fig. 4; Additional file 5: Figure S4).

Conversion of guaiacol and 3-methylcatechol by MtPPO7

For a better discrimination between monophenolase and diphenolase activity, MtPPO7 conversion of guaiacol and 3-methylcatechol was monitored over a period of 24 h by UHPLC-UV-MSn. The conversion of guaiacol by MtPPO7 resulted in the initial formation of 3-methoxycatechol (Additional file 4: Table S1; Additional file 5: Figure S4, Additional file 6: Figure S5). Further reactions resulted in the formation of brown pigments indicating that 3-methoxycatechol was oxidized into ortho-quinones, which are likely to polymerize and form insoluble complexes. Other masses determined by UHPLC-UV-MSn indicated the polymerization of guaiacol and compounds originating from the oxidation of guaiacol by MtPPO7 into trimers (m/z 399, 401, 415). Based on the masses detected (Additional file 4: Table S1; Additional file 5: Figure S4), a scheme is presented of possible reaction pathways of guaiacol occurring during MtPPO7 incubation (Additional file 6: Figure S5).

The oxidation of guaiacol by MtPPO7 mainly occurred between 4 and 8 h of incubation, whereas 3-methylcatechol showed to be oxidized by MtPPO7 within the first 2 h. During the incubation of 3-methylcatechol with MtPPO7 pink pigments are formed, which precipitate after sample centrifugation. Based on MSn, masses of 3-methylcatechol (m/z 123) and masses indicating the dimerization of 3-methylcatechol (m/z 245) were detected upon incubation of both 3-methylcatechol with MtPPO7 and 3-methylcatechol alone (Additional file 4: Table S1). No detectable amounts of new products were formed during the incubation of 3-methylcatechol with MtPPO7 compared to the incubation of 3-methylcatechol only (Additional file 4: Table S1).

Structural model of MtPPO7

A structural model of MtPPO7 was generated based on the available structure of a catechol oxidase from Aspergillus oryzae (AoCO4, Protein Data Bank entry: 4j3p) (Fig. 5). MtPPO7 and AoCO4 share 38% amino acid sequence identity. The structural model of MtPPO7 shows a four-helix bundle fold with the presence of six conserved histidines coordinating the two copper ions in the active site, which is typical for PPO-like tyrosinases and catechol oxidases. Six of the seven cysteines are involved in conserved disulfide bridges (Cys47–Cys393, Cys75–Cys134, and Cys196–Cys234) and expected to be relevant for the thermo-tolerance of MtPPO7 (Fig. 5). Based on the model, the large distance (10.5 Å) between the sulfur atom of Cys302 and Cɛ-atom of His110 prevents formation of a thioether bond, which is present in other PPOs such as AbPPO3 and AbPPO4 [25]. Characteristics of the amino acid sequence and structural model of MtPPO7 are further described in the “Discussion” section.

Fig. 5
figure 5

Structural model of MtPPO7. The structural model of MtPPO7 was generated based on the crystal structure of a catechol oxidase from Aspergillus oryzae (AoCO4, Protein Data Bank entry: 4j3p) [42]. MtPPO7 and AoCO4 share an amino acid sequence identity of 38%. The MtPPO7 model is predominantly α-helical with the catalytic copper site situated in the four-helix bundle. The coordination of the two copper atoms (blue) by six histidine residues (orange) is strictly conserved. The three disulfide bridges Cys47-Cys398, Cys75-Cys134, and Cys196-Cys234 (yellow) are also conserved

Full size image

Genome-wide analysis of AA9 LPMO, AbPPO, and MtPPO7

Sequence annotations of 336 Ascomycota and 208 Basidiomycota genomes [29] were used in this analysis. Some of the numbers of genes obtained encoding AA9 LPMOs were verified by comparison with published data. For example, we have identified 22 and 18 genes encoding AA9 LPMOs in M. thermophila and T. terrestris, respectively, which is similar to previously reported results [30]. In total, 277 Ascomycota genomes and 178 Basidiomycota genomes contained AA9 LPMOs encoding genes (Fig. 6). The two PPOs used in this work are regarded as two different proteins due to their low sequence identity of 12%. The two PPO classes used for the genome analysis were not further divided into short and long tyrosinases [31]. Over 90% of the Ascomycota that comprise genes encoding AA9 LPMOs showed also the presence of genes encoding either AbPPO-like, MtPPO7-like proteins or both (Fig. 6a). In contrast, only around 40% of the AA9 LPMOs encoding Basidiomycota contained genes that encode AbPPO-like, MtPPO7-like proteins or both. The percentage of Ascomycota and Basidiomycota genomes studied that contained neither genes encoding AA9 LPMOs, AbPPOs, nor MtPPO7s were 15 and 8%, respectively (Fig. 6a). The average number of genes found per genome encoding AA9 LPMOs was higher in Ascomycota (13.0) than in Basidiomycota (10.6) (Fig. 6b, c). Both the total and the average number of genes encoding AbPPO-like and MtPPO7-like proteins were also higher in Ascomycota compared to Basidiomycota (Fig. 6b, c).

Fig. 6
figure 6

Genome analysis of Ascomycota and Basidiomycota. The protein sequence annotations of 336 Ascomycota and 208 Basidiomycota genomes were used for genome analysis. In total, a 277 Ascomycota and b 178 Basidiomycota comprise AA9 LPMOs encoding genes (blue column). Next to this occurrence of AA9 LPMOs encoding genes in Ascomycota and Basidiomycota, the co-occurrence of both AbPPO-like and MtPPO7-like, or either AbPPO-like or MtPPO7-like enzymes is indicated on the right of each column. The number of Ascomycota and Basidiomycota without co-occurrence of AA9 LPMOs encoding genes was 59 and 30, respectively (red column). In total, 50 Ascomycota and 17 Basidiomycota species were determined that do not carry any genes encoding AA9 LPMOs, AbPPOs, and MtPPO7s. c Total number of AA9 LPMOs, AbPPOs, and MtPPO7s encoding genes annotated in 336 Ascomycota and 208 Basidiomycota genomes. d Distribution of annotated AA9 LPMOs, AbPPOs, and MtPPO7s encoding genes within 336 Ascomycota and 208 Basidiomycota genomes. See “Methods” for details

Full size image

Principal component analysis (PCA) was performed on all 336 Ascomycota and 208 Basidiomycota. In addition, we used the numbers which describe the presence of AA9 LPMO, MtPPO7s, and AbPPOs encoding genes in each fungal species.

The first two components of the PCA explained 69.8 and 19.7% variation in the data, respectively. Correlations of 0.75 and 0.5 were observed between the presence of genes encoding AA9 LPMOs and MtPPO7s in Ascomycota and Basidiomycota, respectively (Fig. 7b). Lower correlations of 0.59 and 0.34 were observed in the genes encoding AA9 LPMOs and AbPPOs. No correlation (correlation ≤0.25) was observed between AbPPOs and MtPPO7, for both Ascomycota and Basidiomycota. Based on the presence of at least 10 annotated genes that encode cellulose-degrading enzymes per fungus, 27 Ascomycota and 23 Basidiomycota genomes were selected (Additional file 7: Table S2, Additional file 8: Table S3) [32]. The selected 27 Ascomycota genomes showed a higher correlation (0.60) between genes encoding AA9 LPMOs and MtPPO7s compared to the 23 selected Basidiomycota (0.53) (Fig. 7c). Also, no correlation was found between genes encoding MtPPO7- and AbPPO-like genes in the selected Basidiomycota species (Fig. 7).

Fig. 7
figure 7

Statistical genome analysis of Ascomycota and Basidiomycota. a Principal component analysis (PCA) of genes encoding AA9 LPMOs, AbPPOs, and MtPPO7s of 336 Ascomycota and 208 Basidiomycota genomes. The numbers which describe the presence of genes encoding AA9 LPMOs, MtPPO7s, and AbPPOs in each fungal species were used as variables and are presented as principal components of Ascomycota (red dots) and Basidiomycota (blue dots). The first two components of the PCA explained 69.8 and 19.6% variation in the data, respectively. Vectors (red arrows) with a similar orientation illustrate a high correlation between the gene families (b). The axis legends indicate the overall contribution of the three gene families encoding AA9 LPMOs, AbPPOs, and MtPPO7s. b Correlation between the three gene families encoding AA9 LPMOs, AbPPOs, and MtPPO7s of 336 Ascomycota (upper triangle) and 208 Basidiomycota (lower triangle). c Correlation between genes encoding AA9 LPMOs, AbPPOs, and MtPPO7s of 27 Ascomycota (upper triangle) and 23 Basidiomycota (lower triangle), which have at least ten annotated genes encoding cellulose-degrading enzymes per fungus. Species of both fungal classes are listed in Additional file 7: Table S2 and Additional file 8: Table S3. Numerical values of the correlations of b and c are presented in Additional file 9: Figure S6. See “Methods” for details

Full size image

Discussion

MtPPO7 enhances cellulose oxidation by MtLPMO9B

Many PPOs, tyrosinases in particular, have been described to oxidize monophenols and various phenolic compounds comprising a 1,2-benzenediol or 1,2,3-benzenetriol moiety. However, none of these tyrosinases show a high activity towards methoxylated compounds (group Ib) such as ferulic acid (no. 8) [33,34,35,36]. Importantly, MtPPO7 hydroxylates methoxylated phenolic compounds (group Ib) and thereby improves the activity of MtLPMO9B up to 75 times (no. 5, Table 1). On the other hand, AbPPO addition to RAC incubated with MtLPMO9B in the presence non-methoxylated monophenols also resulted in a significant increase of the MtLPMO9B activity (up to 99 times). Especially the activity of MtPPO7 towards methoxylated phenolic compounds is of high relevance, as these compounds are abundant as structural lignin units and, therefore, intrinsically present during plant biomass degradation. Here, we only determined the effect of polyphenol oxidases on the LPMO-mediated cellulose degradation in the presence of phenolic compounds that are known as lignin building blocks but we did not, for example, investigate the impact of these enzymes in the presence of washed or unwashed pretreated biomass. In comparison to the substrate used in this thesis, we expect a lower effect of MtPPO7-like enzymes on the LPMO-mediated substrate oxidation using pretreated biomass, since this biomass contains a mixture of multiple potential electron-donating compounds, which were described in an earlier study [13].

A recent study described that H2O2, in addition to or instead of O2 (Fig. 1), acts as a co-substrate for LPMOs, while reducing agents are still needed to activate the active site copper [37]. Once copper is activated, LPMOs oxidize substrates under the use of H2O2. The latter can be formed by the reduction of O2 through reducing agents [37]. Considering these recent findings, it can only be hypothesized how MtPPO7 increased the MtLPMO9B-mediated cellulose oxidation, which was investigated in our research. It is possible that the oxidation of the phenolic compounds led to the formation of compounds that have either an enhanced reducing efficiency on MtLPMO9B or led to an enhanced H2O2 generation which increased the MtLPMO9B activity.

Activity of MtPPO7 towards phenolic compounds

Based on our results, we conclude that MtPPO7 improves the activity of MtLPMO9B by two main reactions. First, MtPPO7 hydroxylates methoxylated monophenols at the ortho-position (monophenolase activity) and forms compounds comprising a 1,2-benzenediol moiety. These compounds have, compared to methoxylated monophenols, a lower redox potential and are known to be good electron donors for LPMOs [12,13,14]. Secondly, MtPPO7 inhibits the formation of non-enzymatic coupling reactions, which occurred during the incubation of several cinnamic acid derivatives in the absence of MtPPO7, but were not formed when MtPPO7 was present during the incubations (Additional file 4: Table S1). Non-enzymatic coupling reactions of cinnamic acid derivatives of group I, II, or III result in the formation of bulky phenolic polymers. The aliphatic acrylic acid (prop-2-enoic acid) group present in cinnamic acid derivatives constitutes an elongation of the conjugated aromatic ring. This group takes part in polymerization reactions, which can be caused by radical formation due to the presence of copper (2.5 µM) [18, 38]. However, it remains unclear, how MtPPO7 prevents (reduces) the polymerization of cinnamic acid derivatives.

Importantly, MtPPO7 showed only a low efficiency towards non-methoxylated compounds comprising a 1,2-benzenediol or 1,2,3-benzenetriol moiety. In contrast to AbPPO, MtPPO7 had no inhibitory effect on the MtLPMO9B activity towards RAC in the presence of these phenolic compounds (Table 1). Although we did not determine the formation of ortho-quinones, it is possible that MtPPO7 further oxidizes methoxylated compounds comprising a 1,2-benzenediol moiety into ortho-quinones due to the observed pigment formation (Additional file 5: Figure S4, Additional file 6: Figure S5). However, this pigment formation could also partly result from non-enzymatic polymerization reactions or formation of metal-catechol complexes, as indicated in Fig. 8 [39]. Still, based on the monophenolase activity we would consider MtPPO7 to be a tyrosinase. The low efficiency of MtPPO7 to convert compounds comprising a 1,2-benzenediol moiety is expected to be an advantage. MtLPMO9B can utilize these compounds as electron donors and improve the oxidative activity towards cellulose. In the case of a strong diphenolase activity, however, the available amount of these compounds comprising a 1,2-benzenediol moiety would rapidly decrease due to further oxidation of these compounds into ortho-quinones. Hence, we expect for some other PPOs, such as AbPPO, a negative effect on the LPMO-mediated cellulose oxidation. This effect results from the strong activity of these PPOs towards compounds comprising a 1,2-benzenediol or 1,2,3-benzenetriol moiety, which are potential reducing agents for LPMOs [33, 40, 41]. The ortho-quinones formed (no. 17 and 18) (Fig. 2; Table 1) have been shown to be less efficient electron donors for LPMOs than their ortho-diphenol precursors [12].

Fig. 8
figure 8

Proposed reaction pathway of MtPPO7- and AbPPO-mediated oxidation of phenolic compounds. a Phenolic compounds comprising a 2-methoxy moiety are hydroxylated by MtPPO7 from M. thermophila C1 at the 6-position. Different from MtPPO7, AbPPO from A. bisporus shows a low efficiency (dotted arrows) towards these methoxylated monophenols. The formed compounds comprising a 1,2-dihydroxy-3-methoxy moiety have, compared to the methoxylated monophenols, a lower electron potential and an increased electron-donating capacity for LPMOs. b Non-methoxylated monophenols are hydroxylated at the ortho-position by AbPPO into compounds comprising a 1,2-benzenediol moiety, whereas MtPPO7 shows a low efficiency (dotted arrows) towards non-methoxylated monophenols. The formed compounds comprising a a 1,2-dihydroxy-3-methoxy and b 1,2-benzenediol moiety are expected to be further oxidized into ortho-quinones by either the MtPPO7- and AbPPO-mediated or catalyst-mediated oxidation, such as copper which was present during the incubation. As indicated in Additional file 6: Figure S5, these ortho-quinones, including intermediate products formed by the oxidation of phenolic compounds, are expected to polymerize and form insoluble complexes

Full size image

Structural characteristics of MtPPO7

The structural model of MtPPO7 shows a four-helix bundle architecture with the presence of six histidine residues coordinating the two copper ions in the active site, which is typical for PPO-like tyrosinases and catechol oxidases (Fig. 5). We clearly showed that the preference of MtPPO7 towards methoxylated phenolic compounds and the low activity towards diphenols is different from reported fungal or plant catechol oxidases and tyrosinases (Fig. 4). Here, we further discuss the structure–function relationship of these PPOs in comparison with the generated model of MtPPO7. However, it should be noted that the structural model of MtPPO7 was generated based on a catechol oxidase from Aspergillus oryzae (AoCO4) and both proteins share an overall amino acid sequence identity of only 38% [42]. It has already been shown that marginal structural differences determine whether PPOs exhibit a stronger monophenolase or diphenolase activity [43]. Based on our model, it is possible that structural features discussed below deviate compared to the crystal structure which is, however, not available yet.

Firstly, the presence of a thioether bridge between a cysteine and histidine is assumed to be important for the monophenolase and diphenolase activity. This bond has been described to be present in many fungal PPOs, such as AbPPO3/4 from A. bisporus, NcPPO from Neurospora crassa, and also in plant catechol oxidases, such as IbCO from Ipomoea batatas [26, 44, 45]. However, similar to AoCO4 from A. oryzae [42], MtPPO7 does not form this thioether bond (Fig. 5). Secondly, MtPPO7 shows the presence of the PPO typical ‘gate residue’ (Phe300 in MtPPO7, Phe261 in IbCO, and Val299 in AoCO4). It has been hypothesized that a bulky ‘gate residue,’ such as a phenylalanine, limits the accessibility of monophenols to the copper A ion (CuA) and, therefore, might be a limiting factor for the monophenolase activity [44, 46]. However, in AoCO4 the bulky phenylalanine is replaced by Val299 and still AoCO4 shows a rather low activity towards monophenols [31, 42]. Our data suggest that Phe300 in MtPPO7 does not limit the accessibility for (methoxylated) monophenols. Thirdly, it also has recently been proposed that the presence of certain amino acid residues next to the CuB-coordinating histidines indicates whether a tyrosinase or catechol oxidase exhibits a stronger monophenolase and a weak diphenolase activity or a weak monophenolase and a strong diphenolase activity [43]. However, none of the described amino acid residues at that specific position (Pro, Asn, Glu, and Gln) is present in MtPPO7, which has a tyrosine (Tyr268) at the corresponding position. In MtPPO7 a tyrosine (Tyr268) occupies. Whether Tyr268 is responsible for the observed specificity towards methoxylated monophenols can only be hypothesized at this moment. Fourthly, MtPPO7 shares with tyrosinases and catechol oxidases the characteristic ‘tyrosine motif’ (Tyr-X-Tyr/Phe or Tyr/Phe-X-Tyr; in MtPPO7 residues 399-401). This motif is highly conserved among plants and fungi [31, 42, 47]. Fifthly, MtPPO7 does not contain the ‘YG motif’ (Gly-Tyr motif), which is a typical feature of fungal tyrosinases and catechol oxidases [47]. In conclusion, based on substrate specificity (Fig. 1) and modeled structure (Fig. 5), we do not classify MtPPO7 as a typical tyrosinase or catechol oxidase.

Genome analysis

Based on genome analysis, we found a positive correlation between genes encoding AA9 LPMOs and MtPPO7-like enzymes in Ascomycota and Basidiomycota, which has not been shown previously. Interestingly, this correlation strengthens the evidence that fungi benefit from the concerted activity of AA9 LPMOs and MtPPO7-like enzymes in nature, which is in agreement with major findings obtained from the experimental data of this work. The proportion of annotated AA9 LPMOs within Ascomycota and Basidiomycota (82 and 86%, respectively) is about 5–10% lower with the findings of a recent study (92% in Ascomycota and Basidiomycota), which may result from the different species chosen (Fig. 6) [12]. The two PPOs (MtPPO7 and AbPPO) used in this work share only a low sequence identity and, therefore, are considered to belong to different subgroups of polyphenol oxidases. It is important to realize that active tyrosinases, such as the mushroom tyrosinase AbPPO3, have been described to consist of large (H) and small (L) subunits [45, 48]. We focused for the genome analysis on genes encoding the protein of the catalytic H subunit, which contains the highly conserved binuclear copper site. The function of the L subunit is so far not known [45]. All purified enzymes used in this work are secreted. However, we did not differentiate between secreted and non-secreted fungal AA9 LPMOs and PPOs in the genome analysis. It may be that most of the enzymes detected here take part in the external fungal metabolism, as can be deducted from the high abundance of the annotated MtPPO7s and AbPPOs throughout Ascomycota and Basidiomycota. We, therefore, specifically selected genomes of plant biomass degrading fungi that contain at least ten genes encoding known cellulose-degrading enzymes per fungus and found again a similar correlation between genes encoding AA9 LPMOs, AbPPOs, and MtPPO7s (Fig. 7) [32, 49]. In contrast to MtPPO7s, genes encoding AbPPOs showed a weak correlation with AA9 LPMOs in Ascomycota and Basidiomycota (Fig. 7), especially for the selected cellulose-degrading Basidiomycota (Fig. 7c). One possible explanation is that a cellulose-rich environment is associated with the abundance of lignin as seen in soft- and hardwoods. Still, there may be considerable variation in the substrate specificity and specific activity of PPOs within the MtPPO7- or AbPPO-like gene families. The latter seems to be consistent with the low sequence identity among PPOs and, in addition, by the diverse substrate specificities of PPOs that have only marginal differences in their amino acid sequence [31, 43, 50].

Conclusions

For the first time, we demonstrate the importance of the coupled action of different monooxygenases in the concerted degradation of plant biomass. We demonstrated that MtPPO7 is particularly active towards methoxylated phenolic compounds that are the predominant structural units of lignin. This feature distinguishes MtPPO7 from the well-known mushroom tyrosinase AbPPO, which stimulates the LPMO activity via its ability to hydroxylate non-methoxylated monophenols. However, its strong diphenolase activity limits the applicability of this tyrosinase for producing electron-donating capacity for LPMOs. In addition, we established that genes encoding MtPPO7-like enzymes and AA9 LPMOs are highly correlated throughout Ascomycota and Basidiomycota, suggesting that AA9 LPMOs benefit from the activity of MtPPO7-like enzymes in nature. Further understanding in both lignin deconstruction and enzymatic lignocellulose oxidation will lead to more eco-friendly routes for biomass utilization, which is seen as a prerequisite for a circular economy.

Methods

Enzyme expression, production, and purification

MtPPO7 (UniProt: KX772412) was over-expressed in the homologous Myceliophthora thermophila C1 strain. A low protease/low (hemi-)cellulase-producing M. thermophila C1 strain was used to produce MtPPO7 [51, 52]. The MtPPO7-containing culture broth was fractionated to obtain a pure MtPPO7 preparation. This preparation was provided by DuPont Industrial Biosciences. MtLPMO9B was expressed and purified as described in [13].

Protein identification

The pure MtPPO7 fraction was analyzed by LC–mass spectrometry confirming the presence of the MtPPO7 by ‘The Scripps Research Institute’ (San Diego, CA, USA).

Purification and identification of mushroom tyrosinase

Tyrosinase from the edible button mushroom Agaricus bisporus was purified from a commercial enzyme preparation (Sigma-Aldrich, Steinheim, Germany) as described previously [25]. The purified enzyme preparation (referred to as AbPPO) was shown to contain the isoforms PPO3 and PPO4 [25].

Cellulose substrate and reducing agents

Regenerated amorphous cellulose (RAC) was prepared from Avicel PH-101 as described previously [7, 53]. All reducing agents used throughout this study (Table 1) were purchased from Sigma-Aldrich (Steinheim, Germany).

Incubation conditions for MtLPMO9B, MtPPO7, and AbPPO

Regenerated amorphous cellulose (1–2 mg mL−1, see Figure captions) was dissolved in 50 mM ammonium acetate buffer (pH = 5.0), with or without addition of reducing agents (final concentration of 2 mM). The standard enzyme concentrations of MtLPMO9B, MtPPO7, and AbPPO used in this work were 5.0, 5.0, and 0.7 µg protein mg−1 substrate, respectively. All samples were incubated for 20 h at 50 °C in a head-over-tail rotator in portions of 1 mL total volume (Stuart rotator, Bibby Scientific, Stone, UK) at 20 rpm. Supernatants of all incubations were analyzed by HPAEC and RP-UHPLC-UV-ESI–MSn.

Oligosaccharide analysis

Oligosaccharides were analyzed by high-performance anion-exchange chromatography (HPAEC) with pulsed amperometric detection (PAD) using a HPAEC system (ICS-5000, Dionex, Sunnyvale, CA, USA) as described previously [7].

RP-UHPLC-UV-ESI-MSn analysis

Supernatants of all incubations were subjected to an Accela reversed-phase high-performance liquid chromatography (RP-UHPLC) system coupled to electron spray ionization mass spectrometry (Thermo Scientific, San Jose, CA, USA). Injected samples (5 μL) were separated using an Acquity C18 column (2.1 × 150 mm, 1.7 μm particle size) with an Acquity UHPLC Shield RP18 Vanguard guard column (2.1 × 5 mm, 1.7 μm particle size). Both columns were purchased from Waters (Milford, MA, USA). Gradient elution with eluent A (H2O + 1% (v/v) acetonitrile + 0.1% (v/v) HOAc) and eluent B (acetonitrile + 0.1% (v/v) HOAc) was performed according to the following steps: From 0 to 17.7 min a linear gradient from 5 to 60% B; from 17.7 to 21.7 min, isocratic 100% B, and from 21.7 to 26 min, isocratic 5% B. The flow rate and the injection volume were 0.4 mL min−1 and 5 μL, respectively. The column temperature was set to 30 °C and the photodiode array detector was operated in the range of 200–400 nm.

Samples were further analyzed using an LTQ-Velos mass spectrometer (Thermo Scientific) equipped with a ESI–MS. Data were collected over a m/z range of 90–1500 in both negative (NI) and positive (PI) modes. The collision energy was set to 35%.

Structural modeling

An alignment was made of the amino acid sequence of MtPPO7 and the amino acid sequence of catechol oxidase from Aspergillus oryzae (AoCO4), which scored highest in a Blast search using the MtPPO7 sequence against the Protein Data Bank (38% amino acid identity). Using this alignment and the available structure of AoCO4 (PDB id: 4J3P) as template, structural models were obtained for MtPPO7 using the Modeller program version 9.16. One hundred comparative models were generated, after which the model with lowest corresponding DOPE score was selected for inspection and image generation using Pymol (Pymol, The PyMOL Molecular Graphics System, version 1.5.0.4 Schrödinger, LLC, New York, NY, USA).

Genome-wide analysis

Fungal genomes were obtained from the JGI MycoCosm portal [29]. In total, protein sequence annotations of 336 Ascomycota and 208 Basidiomycota genomes were used. BLAST databases for those protein sequences were created using ‘makeblastdb’ program in BLAST + v2.2.30 [54]. The protein BLAST was performed separately using ten MtPPO7 sequences [PDB id: 4j3p and closely related MtPPO7s: Q2UNF9, Aspergillus oryzae (strain ATCC 42149); G2QC95, Myceliophthora thermophila (ATCC 42464); Q2H7I7, Chaetomium globosum (ATCC 6205); G0SFX8, Chaetomium thermophilum (DSM 1495); L7IAQ4, Magnaporthe oryzae (strain Y34); L7JMT9, Magnaporthe oryzae (strain P131); G4N2I5, Magnaporthe oryzae (strain 70-15); A0A084GCK1, Scedosporium apiospermum; A0A0C4DYF2, Magnaporthiopsis poae (ATCC 64411); J3P591, Gaeumannomyces graminis var. tritici (strain R3-111a-1)], ten AA9 LPMO sequences (PDB id: 4d7u and closely related LPMOs: Q7SHI8, Neurospora crassa (strain ATCC 24698); G2QCJ3, Myceliophthora thermophila (strain ATCC 42464); F7W1P4, Sordaria macrospora (strain ATCC MYA-333); G2RB73, Thielavia terrestris (strain ATCC 38088); Q2H8N9, Chaetomium globosum (strain ATCC 6205); G0S408, Chaetomium thermophilum (strain DSM 1495); F8MLY8, Neurospora tetrasperma (strain FGSC 2508); T0L448, Colletotrichum gloeosporioides (strain Cg-14); A0A0H4K9X4 and A0A1C9CXI0, Myceliophthora thermophila C1); and four AbPPOs (AbPPO1, AbPPO2, AbPPO3, and AbPPO4 [55, 56]) as query sequences. Resulting sequences below E-value cut-off of 0.001 with query coverage above 60% for AA9 LPMOs, 65% for MtPPO7s, and 85% for AbPPOs were considered for further analysis. Selection of cellulase-rich Ascomycota and Basidiomycota was based on the presence of at least 10 genes encoding cellulose-degrading enzymes, which are classified in the CAZy database as glycosyl hydrolase families GH1, GH3, GH5, GH6, GH7, GH12, GH45. The GH gene families were selected based on Kubicek et al. [57]. Previous data [32] were used to determine the number of annotated genes encoding cellulose-degrading enzymes. Based on this selection, 27 Ascomycota and 23 Basidiomycota species were selected for the Pearson correlation analysis (Fig. 7; Additional file 7: Table S2, Additional file 8: Table S3). All the statistical analyses were performed in R [58].

References

  1. Lam TBT, Kadoya K, Iiyama K. Bonding of hydroxycinnamic acids to lignin: ferulic and p-coumaric acids are predominantly linked at the benzyl position of lignin, not the β-position, in grass cell walls. Phytochemistry. 2001;57(6):987–92. doi:10.1016/S0031-9422(01)00052-8.

    Article  CAS  Google Scholar 

  2. Takahashi N, Koshijima T. Ester linkages between lignin and glucuronoxylan in a lignin-carbohydrate complex from beech (Fagus crenata) wood. Wood Sci Technol. 1988;22(231–241):231–41. doi:10.1007/BF00386018.

    Article  CAS  Google Scholar 

  3. Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sorlie M, et al. An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science. 2010;330(6001):219–22. doi:10.1126/science.1192231.

    Article  CAS  Google Scholar 

  4. Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014:42;D490–5. doi:10.1093/nar/gkt1178.

    Article  CAS  Google Scholar 

  5. Phillips CM, Beeson WT, Cate JH, Marletta MA. Cellobiose dehydrogenase and a copper-dependent polysaccharide monooxygenase potentiate cellulose degradation by Neurospora crassa. ACS Chem Biol. 2011;6(12):1399–406. doi:10.1021/cb200351y.

    Article  CAS  Google Scholar 

  6. Quinlan RJ, Sweeney MD, Lo Leggio L, Otten H, Poulsen JC, Johansen KS, et al. Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci USA. 2011;108(37):15079–84. doi:10.1073/pnas.1105776108.

    Article  CAS  Google Scholar 

  7. Frommhagen M, Sforza S, Westphal AH, Visser J, Hinz SW, Koetsier MJ, et al. Discovery of the combined oxidative cleavage of plant xylan and cellulose by a new fungal polysaccharide monooxygenase. Biotechnol Biofuels. 2015;8:101. doi:10.1186/s13068-015-0284-1.

    Article  Google Scholar 

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

    Article  CAS  Google Scholar 

  9. Isaksen T, Westereng B, Aachmann FL, Agger JW, Kracher D, Kittl R, et al. A C4-oxidizing lytic polysaccharide monooxygenase cleaving both cellulose and cello-oligosaccharides. J Biol Chem. 2014;289(5):2632–42. doi:10.1074/jbc.M113.530196.

    Article  CAS  Google Scholar 

  10. Forsberg Z, Rohr AK, Mekasha S, Andersson KK, Eijsink VG, Vaaje-Kolstad G, et al. Comparative study of two chitin-active and two cellulose-active AA10-type lytic polysaccharide monooxygenases. Biochemistry. 2014;53(10):1647–56. doi:10.1021/bi5000433.

    Article  CAS  Google Scholar 

  11. Agger JW, Isaksen T, Varnai A, Vidal-Melgosa S, Willats WG, Ludwig R, et al. Discovery of LPMO activity on hemicelluloses shows the importance of oxidative processes in plant cell wall degradation. Proc Natl Acad Sci USA. 2014;111(17):6287–92. doi:10.1073/pnas.1323629111.

    Article  CAS  Google Scholar 

  12. Kracher D, Scheiblbrandner S, Felice AKG, Breslmayr E, Preims M, Ludwicka K, et al. Extracellular electron transfer systems fuel cellulose oxidative degradation. Science. 2016. doi:10.1126/science.aaf3165.

    Google Scholar 

  13. Frommhagen M, Koetsier MJ, Westphal AH, Visser J, Hinz SWA, Vincken J-P, et al. Lytic polysaccharide monooxygenases from Myceliophthora thermophila C1 differ in substrate preference and reducing agent specificity. Biotechnol Biofuels. 2016;9(1):1–17. doi:10.1186/s13068-016-0594-y.

    Article  Google Scholar 

  14. Westereng B, Cannella D, Wittrup Agger J, Jørgensen H, Larsen Andersen M, Eijsink VGH, et al. Enzymatic cellulose oxidation is linked to lignin by long-range electron transfer. Scientific Rep. 2015;5:18561. doi:10.1038/srep18561.

    Article  CAS  Google Scholar 

  15. Bissaro B, Forsberg Z, Ni Y, Hollmann F, Vaaje-Kolstad G, Eijsink VGH. Fueling biomass-degrading oxidative enzymes by light-driven water oxidation. Green Chem. 2016;18:5357–66. doi:10.1039/c6gc01666a.

    Article  CAS  Google Scholar 

  16. Cannella D, Mollers KB, Frigaard NU, Jensen PE, Bjerrum MJ, Johansen KS, et al. Light-driven oxidation of polysaccharides by photosynthetic pigments and a metalloenzyme. Nat Commun. 2016;7:11134. doi:10.1038/ncomms11134.

    Article  CAS  Google Scholar 

  17. Gao Z-J, Liu J-B, Xiao X-G. Purification and characterisation of polyphenol oxidase from leaves of Cleome gynandra L. Food Chem. 2011;129(3):1012–8. doi:10.1016/j.foodchem.2011.05.062.

    Article  CAS  Google Scholar 

  18. Dec J, Haider K, Bollag J-M. Decarboxylation and demethylation of naturally occurring phenols during coupling reactions and polymerization. Soil Sci. 2001;166(10):660–71.

    Article  CAS  Google Scholar 

  19. Kuijpers TFM, van Herk T, Vincken J-P, Janssen RH, Narh DL, van Berkel WJH, et al. Potato and mushroom polyphenol oxidase activities are differently modulated by natural plant extracts. J Agric Food Chem. 2014;62(1):214–21. doi:10.1021/jf4043375.

    Article  CAS  Google Scholar 

  20. Beloqui A, Pita M, Polaina J, Martínez-Arias A, Golyshina OV, Zumárraga M, et al. Novel polyphenol oxidase mined from a metagenome expression library of bovine rumen: biochemical properties, structural analysis, and phylogenetic relationships. J Biol Chem. 2006;281(32):22933–42. doi:10.1074/jbc.M600577200.

    Article  CAS  Google Scholar 

  21. Lerch K. Neurospora tyrosinase: structural, spectroscopic and catalytic properties. Mol Cell Biochem. 1983;52(2):125–38. doi:10.1007/bf00224921.

    Article  CAS  Google Scholar 

  22. Konrad L. Tyrosinase: molecular and active-site structure. Enzymatic browning and its prevention. ACS Symposium Series: American Chemical Society; 1995. p. 64–80.

  23. Vanholme R, Demedts B, Morreel K, Ralph J, Boerjan W. Lignin biosynthesis and structure. Plant Physiol. 2010;153(3):895–905. doi:10.1104/pp.110.155119.

    Article  CAS  Google Scholar 

  24. Kuijpers TFM, Narváez-Cuenca C-E, Vincken J-P, Verloop AJW, van Berkel WJH, Gruppen H. Inhibition of enzymatic browning of chlorogenic acid by sulfur-containing compounds. J Agric Food Chem. 2012;60(13):3507–14. doi:10.1021/jf205290w.

    Article  CAS  Google Scholar 

  25. Kuijpers TFM, Gruppen H, Sforza S, van Berkel WJH, Vincken J-P. The antibrowning agent sulfite inactivates Agaricus bisporus tyrosinase through covalent modification of the copper-B site. FEBS J. 2013;280(23):6184–95. doi:10.1111/febs.12539.

    Article  CAS  Google Scholar 

  26. Mauracher SG, Molitor C, Al-Oweini R, Kortz U, Rompel A. Latent and active abPPO4 mushroom tyrosinase cocrystallized with hexatungstotellurate(VI) in a single crystal. Acta Crystallographica Sect D. 2014;70(9):2301–15. doi:10.1107/S1399004714013777.

    Article  CAS  Google Scholar 

  27. Mauracher SG, Molitor C, Al-Oweini R, Kortz U, Rompel A. Crystallization and preliminary X-ray crystallographic analysis of latent isoform PPO4 mushroom (Agaricus bisporus) tyrosinase. Acta Crystallographica Secti F Struct Biol Commun. 2014;70(Pt 2):263–6. doi:10.1107/s2053230x14000582.

    Article  CAS  Google Scholar 

  28. Wu J, Chen H, Gao J, Liu X, Cheng W, Ma X. Cloning, characterization and expression of two new polyphenol oxidase cDNAs from Agaricus bisporus. Biotechnol Lett. 2010;32(10):1439–47. doi:10.1007/s10529-010-0329-2.

    Article  CAS  Google Scholar 

  29. Grigoriev IV, Nikitin R, Haridas S, Kuo A, Ohm R, Otillar R, et al. MycoCosm portal: gearing up for 1000 fungal genomes. Nucleic Acids Res. 2014;42(Database issue):D699–704. doi:10.1093/nar/gkt1183.

    Article  CAS  Google Scholar 

  30. Berka RM, Grigoriev IV, Otillar R, Salamov A, Grimwood J, Reid I, et al. Comparative genomic analysis of the thermophilic biomass-degrading fungi Myceliophthora thermophila and Thielavia terrestris. Nat Biotechnol. 2011;29(10):922–7. doi:10.1038/nbt.1976.

    Article  CAS  Google Scholar 

  31. Gasparetti C, Faccio G, Arvas M, Buchert J, Saloheimo M, Kruus K. Discovery of a new tyrosinase-like enzyme family lacking a C-terminally processed domain: production and characterization of an Aspergillus oryzae catechol oxidase. Appl Microbiol Biotechnol. 2010;86(1):213–26. doi:10.1007/s00253-009-2258-3.

    Article  CAS  Google Scholar 

  32. Zhao Z, Liu H, Wang C, Xu J-R. Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genom. 2013;14(1):1–15. doi:10.1186/1471-2164-14-274.

    Article  Google Scholar 

  33. Selinheimo E, NiEidhin D, Steffensen C, Nielsen J, Lomascolo A, Halaouli S, et al. Comparison of the characteristics of fungal and plant tyrosinases. J Biotechnol. 2007;130(4):471–80. doi:10.1016/j.jbiotec.2007.05.018.

    Article  CAS  Google Scholar 

  34. Zhou P, Smith NL, Lee CY. Potential purification and some properties of Monroe apple peel polyphenol oxidase. J Agric Food Chem. 1993;41(4):532–6. doi:10.1021/jf00028a004.

    Article  CAS  Google Scholar 

  35. Paul B, Gowda LR. Purification and characterization of a polyphenol oxidase from the seeds of field bean (Dolichos lablab). J Agric Food Chem. 2000;48(9):3839–46. doi:10.1021/jf000296s.

    Article  CAS  Google Scholar 

  36. Zou Y, Hu W, Jiang A, Ma K. Partial purification and characterization of a novel extracellular tyrosinase from Auricularia auricula. Appl Biochem Biotechnol. 2013;172(3):1460–9. doi:10.1007/s12010-013-0638-8.

    Article  Google Scholar 

  37. Bissaro B, Rohr AK, Skaugen M, Forsberg Z, Horn SJ, Vaaje-Kolstad G, et al. Fenton-type chemistry by a copper enzyme: molecular mechanism of polysaccharide oxidative cleavage. bioRxiv. 2016. doi:10.1101/097022.

    Google Scholar 

  38. Rosazza JPN, Huang Z, Dostal L, Volm T, Rousseau B. Review: biocatalytic transformations of ferulic acid: an abundant aromatic natural product. J Ind Microbiol. 1995;15(6):457–71. doi:10.1007/bf01570016.

    Article  CAS  Google Scholar 

  39. Schweigert N, Zehnder AJB, Eggen RIL. Chemical properties of catechols and their molecular modes of toxic action in cells, from microorganisms to mammals. Environ Microbiol. 2001;3(2):81–91. doi:10.1046/j.1462-2920.2001.00176.x.

    Article  CAS  Google Scholar 

  40. Raymond J, Rakariyatham N, Azanza JL. Purification and some properties of polyphenoloxidase from sunflower seeds. Phytochemistry. 1993;34(4):927–31. doi:10.1016/S0031-9422(00)90689-7.

    Article  CAS  Google Scholar 

  41. Lee CY, Smith NL, Pennesi AP. Polyphenoloxidase from DeChaunac grapes. J Sci Food Agric. 1983;34(9):987–91. doi:10.1002/jsfa.2740340915.

    Article  CAS  Google Scholar 

  42. Hakulinen N, Gasparetti C, Kaljunen H, Kruus K, Rouvinen J. The crystal structure of an extracellular catechol oxidase from the ascomycete fungus Aspergillus oryzae. JBIC, J Biol Inorg Chem. 2013;18(8):917–29. doi:10.1007/s00775-013-1038-9.

    Article  CAS  Google Scholar 

  43. Solem E, Tuczek F, Decker H. Tyrosinase versus catechol oxidase: one asparagine makes the difference. Angewandte Chemie Int Edit. 2016;55(8):2884–8. doi:10.1002/anie.201508534.

    Article  CAS  Google Scholar 

  44. Klabunde T, Eicken C, Sacchettini JC, Krebs B. Crystal structure of a plant catechol oxidase containing a dicopper center. Nat Struct Biol. 1998;5(12):1084–90. doi:10.1038/4193.

    Article  CAS  Google Scholar 

  45. Ismaya WT, Rozeboom HJ, Weijn A, Mes JJ, Fusetti F, Wichers HJ, et al. Crystal structure of Agaricus bisporus mushroom tyrosinase: identity of the tetramer subunits and interaction with tropolone. Biochemistry. 2011;50(24):5477–86. doi:10.1021/bi200395t.

    Article  CAS  Google Scholar 

  46. Decker H, Schweikardt T, Tuczek F. The first crystal structure of tyrosinase: all questions answered? Angewandte Chemie Int Edit. 2006;45(28):4546–50. doi:10.1002/anie.200601255.

    Article  CAS  Google Scholar 

  47. Marusek CM, Trobaugh NM, Flurkey WH, Inlow JK. Comparative analysis of polyphenol oxidase from plant and fungal species. J Inorg Biochem. 2006;100(1):108–23. doi:10.1016/j.jinorgbio.2005.10.008.

    Article  CAS  Google Scholar 

  48. Strothkamp KG, Jolley RL, Mason HS. Quaternary structure of mushroom tyrosinase. Biochem Biophys Res Commun. 1976;70(2):519–24. doi:10.1016/0006-291X(76)91077-9.

    Article  CAS  Google Scholar 

  49. Kubicek CP, Kubicek EM. Enzymatic deconstruction of plant biomass by fungal enzymes. Curr Opin Chem Biol. 2016;35:51–7. doi:10.1016/j.cbpa.2016.08.028.

    Article  CAS  Google Scholar 

  50. Goldfeder M, Kanteev M, Adir N, Fishman A. Influencing the monophenolase/diphenolase activity ratio in tyrosinase. Biochimica et Biophysica Acta Proteins Proteom. 2013;1834(3):629–33. doi:10.1016/j.bbapap.2012.12.021.

    Article  CAS  Google Scholar 

  51. Visser H, Joosten V, Punt PJ, Gusakov AV, Olson PT, Joosten R, et al. Development of a mature fungal technology and production platform for industrial enzymes based on a Myceliophthora thermophila isolate, previously known as Chrysosporium lucknowense C1. Ind Biotechnol. 2011;7:214–23. doi:10.1089/ind.2011.0003.

    Article  CAS  Google Scholar 

  52. Punt PJ, Burlingame RP, Pynnonen CM, Olson PT, Wery J, Visser J, Heinrich, et al. Chrysosporium lucknowense protein production system. Patent WO/2010/107303; 2010.

  53. Zhang YHP, Cui J, Lynd LR, Kuang LR. A transition from cellulose swelling to cellulose dissolution by o-phosphoric acid: evidence from enzymatic hydrolysis and supramolecular structure. Biomacromol. 2006;7(2):644–8. doi:10.1021/bm050799c.

    Article  CAS  Google Scholar 

  54. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST + : architecture and applications. BMC Bioinform. 2009;10(1):1–9. doi:10.1186/1471-2105-10-421.

    Article  Google Scholar 

  55. Wichers HJ, Recourt K, Hendriks M, Ebbelaar CEM, Biancone G, Hoeberichts FA, et al. Cloning, expression and characterisation of two tyrosinase cDNAs from Agaricus bisporus. Appl Microbiol Biotechnol. 2003;61(4):336–41. doi:10.1007/s00253-002-1194-2.

    Article  CAS  Google Scholar 

  56. N-y Li, W-m Cai, Q-l Jin, Q-p Qin, F-l Ran. Molecular cloning and expression of polyphenoloxidase genes from the mushroom, Agaricus bisporus. Agric Sci China. 2011;10(2):185–94. doi:10.1016/S1671-2927(09)60305-9.

    Article  Google Scholar 

  57. Kubicek CP, Starr TL, Glass NL. Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Ann Rev Phytopathol. 2014;52(1):427–51. doi:10.1146/annurev-phyto-102313-045831.

    Article  Google Scholar 

  58. Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016. https://www.R-project.org.

Download references

Authors’ contributions

MF, MJK, JV, SWAH, DW, J-PV, and MAK designed research; MF performed research and analyzed data; SKM and MF performed genome analysis; AHW and MF performed structural modeling; and MF, SKM, AHW, WJHB, HG, and MAK wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

This research was funded by The Graduate School VLAG (Advanced Studies in Food Technology, Agrobiotechnology, Nutrition and Health Sciences) of Wageningen University & Research (Wageningen, The Netherlands) and by a VICI Grant (865.14.001) of the Netherlands Organization for Scientific Research (NWO) to D.W.

Competing interests

The authors declare that they have no competing interests.

Availability of supporting data

Supporting data will be provided as Additional files 1, 2, 3, 4, 5, 6, 7, 8, 9.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

  1. Laboratory of Food Chemistry, Wageningen University & Research, Bornse Weilanden 9, 6708 WG, Wageningen, The Netherlands

    Matthias Frommhagen, Jean-Paul Vincken, Harry Gruppen & Mirjam A. Kabel

  2. Laboratory of Biochemistry, Wageningen University & Research, Stippeneng 4, 6708 WE, Wageningen, The Netherlands

    Sumanth Kumar Mutte, Adrie H. Westphal, Dolf Weijers & Willem J. H. van Berkel

  3. DuPont Industrial Biosciences, Nieuwe Kanaal 7-S, 6709 PA, Wageningen, The Netherlands

    Martijn J. Koetsier & Sandra W. A. Hinz

  4. Fungal Genetics & Technology Consultancy, P.O. Box 39b, 6700 AJ, Wageningen, The Netherlands

    Jaap Visser

Authors
  1. Matthias Frommhagen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Sumanth Kumar Mutte
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Adrie H. Westphal
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Martijn J. Koetsier
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Sandra W. A. Hinz
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jaap Visser
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jean-Paul Vincken
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Dolf Weijers
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Willem J. H. van Berkel
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Harry Gruppen
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Mirjam A. Kabel
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Mirjam A. Kabel.

Additional files

13068_2017_810_MOESM1_ESM.png

Additional file 1: Figure S1. Activity of MtLPMO9B towards amorphous cellulose in the presence and absence of MtPPO7 or AbPPO. HPAEC elution pattern of regenerated amorphous cellulose (RAC; 1.5 mg mL−1) incubated with MtLPMO9B (red, 5.0 μg mL−1) only, or with either AbPPO (blue, 2.5 µL mL−1) or MtPPO7 (yellow, 5.0 μg mL−1) in the presence of (a) para-coumaric acid (no. 3 specified in Table 1, 2 mM) and (b) 3-hydroxy-4-methoxycinnamic acid (no. 5 specified in Table 1, 2 mM). The incubation of RAC with MtLPMO9B results in the formation of non-oxidized gluco-oligosaccharides (GlcOSn) and C1-oxidized gluco-oligosaccharides (GlcOS #n ). See “Methods” for details.

13068_2017_810_MOESM2_ESM.png

Additional file 2: Figure S2. Release of oligosaccharides from RAC incubated with MtLPMO9B in the presence and absence of MtPPO7 throughout 24 h. Samples were incubated in the presence of ferulic acid (no. 8 specified in Table 1). The total sum is shown as integrated peak areas of released non-oxidized (shaded red and shaded yellow) and C1-oxidized (red and yellow) gluco-oligosaccharides after incubation of regenerated amorphous cellulose (RAC; 1.5 mg mL−1) with MtLPMO9B only (red bars, 5 mg mL−1) and MtLPMO9B together with MtPPO7 (yellow bars, 5 mg mL−1) based on HPAEC. All incubations were performed in duplicate, and the standard deviations are presented as error bars. See “Methods” for details.

13068_2017_810_MOESM3_ESM.png

Additional file 3: Figure S3. Concentration of phenolic compounds incubated in the presence and absence of MtPPO7. (a) guaiacol (no. 9 specified in Table 1, 2 mM) and (c) 3-methylcatechol (no. 17 specified in Table 1, 2 mM) were incubated with MtPPO7 (yellow bar, 5 μg mL−1) or without (red bar). Samples were incubated in a 50 mM potassium phosphate buffer (pH = 6.0) containing 2.5 µM copper(II)-chloride for 24 h at 50 °C. The conversion of guaiacol and 3-methylcatechol by MtPPO7 was calculated by subtracting the determined concentration of the incubation of guaiacol or 3-methylcatechol in the presence of MtPPO7 from the concentration that was determined by the incubation of guaiacol and 3-methylcatechol alone. This conversion was expressed as the relative decrease of the guaiacol and 3-methylcatechol concentration and is shown in (b) and (d), respectively. See “Methods” for details.

13068_2017_810_MOESM4_ESM.xlsx

Additional file 4: Table S1. Compounds detected after incubation of 21 reducing agents without or with MtPPO7 by UHPLC/UV-MS1.

13068_2017_810_MOESM5_ESM.png

Additional file 5: Figure S4. UHPLC-UV-MSn elution profile of guaiacol incubated (a) in the presence and (b) absence of MtPPO7. Guaiacol (no. 9 specified in Table 1, 2 mM) was incubated with (5 μg mL−1) or without MtPPO7. Samples were incubated in a 50 mM potassium phosphate buffer (pH = 6.0) containing 2.5 µM copper(II)-chloride for 24 h at 50 °C. Annotation of the peaks based on UV was done by using mass spectrometry (Additional file 6: Figure S5). See “Methods” for details.

13068_2017_810_MOESM6_ESM.png

Additional file 6: Figure S5. Schematic presentation of possible reaction pathways of guaiacol incubated in the presence of MtPPO7. In short, MtPPO7 hydroxylates guaiacol (no. 9 specified in Table 1) into 3-methoxycatechol (monophenolase activity). Although not determined, it is likely that 3-methoxycatechol is further oxidized by MtPPO7 into the corresponding ortho-quinone (diphenolase activity, dashed arrow). These ortho-quinones are expected to polymerize and form insoluble complexes. Guaiacol itself forms insoluble complexes via auto-oxidation, which results from the presence of copper during the incubation for 24 h at 50 °C. The decrease in guaiacol concentration during the incubation of guaiacol without MtPPO7 is also shown in Additional file 3: Figure S3a. The determined masses indicate the presence of multiple trimers (399.06, 401.06 and 415.04) consisting of polymerized 3-methylcatechol and guaiacol (Additional file 5: Figure S4). As described above, the polymerization reactions are expected to be catalyzed by copper during the incubation conditions applied. All masses were determined by UHPLC-UV-MSn after incubation of guaiacol (2 mM) with MtPPO7 (5.0 μg mL−1). Samples were incubated for 24 h at 50 °C in 50 mM potassium phosphate, pH 6.0, containing 2.5 µM copper(II)-chloride.

Additional file 7: Table S2. Selected cellulase-rich Ascomycota from the JGI database1.

Additional file 8: Table S3. Selected cellulase-rich Basidiomycota from the JGI database1.

13068_2017_810_MOESM9_ESM.png

Additional file 9: Figure S6. Correlation of AA9 LPMOs, AbPPOs and MtPPO7s encoding genes in Ascomycota and Basidiomycota. Correlation between the three gene families encoding AA9 LPMOs, AbPPOs and MtPPO7s of (a) 336 Ascomycota and (b) 208 Basidiomycota. Correlation between genes encoding AA9 LPMOs, AbPPOs and MtPPO7s of (c) 27 Ascomycota and (d) 23 Basidiomycota, which have at least ten annotated genes encoding cellulose-degrading enzymes. Species of selected fungal classes are listed in Additional file 7: Table S2 and Additional file 8: Table S3. Graphical presentations of the correlations of (a) together with (b) and (c) and (d) are shown in Fig. 7. See “Methods” for details.

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

Frommhagen, M., Mutte, S.K., Westphal, A.H. et al. Boosting LPMO-driven lignocellulose degradation by polyphenol oxidase-activated lignin building blocks. Biotechnol Biofuels 10, 121 (2017). https://doi.org/10.1186/s13068-017-0810-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-017-0810-4

Keywords

Biotechnology for Biofuels and Bioproducts

ISSN: 2731-3654

Contact us