Molecular Evolution and Diversity of Lignin Degrading Heme Peroxidases in the Agaricomycetes

The plant and microbial peroxidase superfamily encompasses three classes of related protein families. Class I includes intracellular peroxidases of prokaryotic origin, class II includes secretory fungal peroxidases, including the lignin degrading enzymes manganese peroxidase (MnP), lignin peroxidase (LiP), and versatile peroxidase (VP), and class III includes the secretory plant peroxidases. Here, we present phylogenetic analyses using maximum parsimony and Bayesian methods that address the origin and diversification of class II peroxidases. Higher-level analyses used published full-length sequences from all members of the plant and microbial peroxidase superfamily, while lower-level analyses used class II sequences only, including 43 new sequences generated from Agaricomycetes (mushroom-forming fungi and relatives). The distribution of confirmed and proposed catalytic sites for manganese and aromatic compounds in class II peroxidases, including residues supposedly involved in three different long range electron transfer pathways, was interpreted in the context of phylogenies from the lower-level analyses. The higher-level analyses suggest that class II sequences constitute a monophyletic gene family within the plant and microbial peroxidase superfamily, and that they have diversified extensively in the basidiomycetes. Peroxidases of unknown function from the ascomycete Magnaporthe grisea were found to be the closest relatives of class II sequences and were selected to root class II sequences in the lower-level analyses. LiPs evidently arose only once in the Polyporales, which harbors many white-rot taxa, whereas MnPs and VPs are more widespread and may have multiple origins. Our study includes the first reports of partial sequences for MnPs in the Hymenochaetales and Corticiales.     

Fungal biodegradation and enzymatic modification of lignin

Lignin, the most abundant aromatic biopolymer on Earth, is extremely recalcitrant to degradation. By linking to both hemicellulose and cellulose, it creates a barrier to any solutions or enzymes and prevents the penetration of lignocellulolytic enzymes into the interior lignocellulosic structure. Some basidiomycetes white-rot fungi are able to degrade lignin efficiently using a combination of extracellular ligninolytic enzymes, organic acids, mediators and accessory enzymes. This review describes ligninolytic enzyme families produced by these fungi that are involved in wood decay processes, their molecular structures, biochemical properties and the mechanisms of action which render them attractive candidates in biotechnological applications. These enzymes include phenol oxidase (laccase) and heme peroxidases [lignin peroxidase (LiP), manganese peroxidase (MnP) and versatile peroxidase (VP)]. Accessory enzymes such as H₂O₂-generating oxidases and degradation mechanisms of plant cell-wall components in a non-enzymatic manner by production of free hydroxyl radicals (·OH) are also discussed.     

Biodegradation of lignocellulosics: microbial, chemical, and enzymatic aspects of the fungal attack of lignin.

Wood is the main renewable material on Earth and is largely used as building material and in paper-pulp manufacturing. This review describes the composition of lignocellulosic materials, the different processes by which fungi are able to alter wood, including decay patterns caused by white, brown, and soft-rot fungi, and fungal staining of wood. The chemical, enzymatic, and molecular aspects of the fungal attack of lignin, which represents the key step in wood decay, are also discussed. Modern analytical techniques to investigate fungal degradation and modification of the lignin polymer are reviewed, as are the different oxidative enzymes (oxidoreductases) involved in lignin degradation. These include laccases, high redox potential ligninolytic peroxidases (lignin peroxidase, manganese peroxidase, and versatile peroxidase), and oxidases. Special emphasis is given to the reactions catalyzed, their synergistic action on lignin, and the structural bases for their unique catalytic properties. Broadening our knowledge of lignocellulose biodegradation processes should contribute to better control of wood-decaying fungi, as well as to the development of new biocatalysts of industrial interest based on these organisms and their enzymes.     

Continuous production of ligin peroxidase by phanerochaete chrysosporium immobilized on a sintered glass carrier

Phanerochaete chrysosporium cells were immobilized on the sintered porous glass support. Such a biocatalizer was used as a bed of the enzymatic reactor system for the continuous production of lignin peroxidase. From the after culture fluid the lignin peroxidase enzymatic activity was recovered and purified applying anion exchangers. Additionally, some physico-chemical properties of lignin peroxidases were determined.     

Role of laccases and peroxidases in saprotrophic activities in the lichen Usnea undulata

Lichens produce various oxidoreductases including heme-containing peroxidases and the copper-containing phenol oxidases tyrosinase and laccase. Our earlier findings suggested that significant oxidoreductase activity occurs mainly in lichens from the order Peltigerales. Here we show that the non-Peltigeralean lichen Usnea can display significant activities of peroxidases and laccases. Strong evidence for the involvement of peroxidases and laccases in saprotrophic activities comes from the observation that their activities are induced by “starvation” due to prolonged dark storage, and also by treatment with soluble cellulose and lignin breakdown products. We also show that, given a quinone and chelated Fe, Usnea can produce hydroxyl radicals; these radicals contribute to the break down of carbohydrates or lignin. However, hydroxyl radical production is independent of laccase and peroxidase activity. Laccases and peroxidases are involved in other aspects of lichen biology; here we show that peroxidases, but not laccases, can break down lichen substances. Reduction in the amounts of lichen substances will reduce photoprotection, which will increase the photosynthetic capacity of thalli during winter when light intensities are low.     

Computational analysis of the Phanerochaete chrysosporium v2.0 genome database and mass spectrometry identification of peptides in ligninolytic cultures reveal complex mixtures of secreted proteins

The white-rot basidiomycete Phanerochaete chrysosporium employs extracellular enzymes to completely degrade the major polymers of wood: cellulose, hemicellulose, and lignin. Analysis of a total of 10,048 v2.1 gene models predicts 769 secreted proteins, a substantial increase over the 268 models identified in the earlier database (v1.0). Within the v2.1 'computational secretome,' 43% showed no significant similarity to known proteins, but were structurally related to other hypothetical protein sequences. In contrast, 53% showed significant similarity to known protein sequences including 87 models assigned to 33 glycoside hydrolase families and 52 sequences distributed among 13 peptidase families. When grown under standard ligninolytic conditions, peptides corresponding to 11 peptidase genes were identified in culture filtrates by mass spectrometry (LS-MS/MS). Five peptidases were members of a large family of aspartyl proteases, many of which were localized to gene clusters. Consistent with a role in dephosphorylation of lignin peroxidase, a mannose-6-phosphatase (M6Pase) was also identified in carbon-starved cultures. Beyond proteases and M6Pase, 28 specific gene products were identified including several representatives of gene families. These included 4 lignin peroxidases, 3 lipases, 2 carboxylesterases, and 8 glycosyl hydrolases. The results underscore the rich genetic diversity and complexity of P. chrysosporium's extracellular enzyme systems.     

Asynchronous origins of ectomycorrhizal clades of Agaricales

The ectomycorrhizal (ECM) symbiosis is the most widespread biotrophic nutritional mode in mushroom-forming fungi. ECM fungi include, though are not limited to, about 5000 described species of Agaricales from numerous, independently evolved lineages. Two central hypotheses suggest different explanations for the origin of ECM fungal diversity: (i) dual origins, initially with the Pinaceae in the Jurassic and later with angiosperms during the Late Cretaceous, and (ii) a simultaneous and convergent radiation of ECM lineages in response to cooling climate during the Palaeogene and advancing temperate ECM plant communities. Neither of these hypotheses is supported here. While we demonstrate support for asynchronous origins of ECM Agaricales, the timing of such events appears to have occurred more recently than suggested by the first hypothesis, first during the Cretaceous and later during the Palaeogene. We are also unable to reject models of rate constancy, which suggests that the diversity of ECM Agaricales is not a consequence of convergent rapid radiations following evolutionary transitions from saprotrophic to ECM habits. ECM lineages of Agaricales differ not only in age, but also in rates of diversification and rate of substitution at nuclear ribosomal RNA loci. These results question the biological uniformity of the ECM guild.     

Phylogenetic relationships in class I of the superfamily of bacterial, fungal, and plant peroxidases.

Molecular phylogeny among catalase-peroxidases, cytochrome c peroxidases, and ascorbate peroxidases was analysed. Sixty representative sequences covering all known subgroups of class I of the superfamily of bacterial, fungal, and plant heme peroxidases were selected. Each sequence analysed contained the typical peroxidase motifs evolved to bind effectively the prosthetic heme group, enabling peroxidatic activity. The N-terminal and C-terminal domains of catalase-peroxidases matching the ancestral tandem gene duplication event were treated separately in the phylogenetic analysis to reveal their specific evolutionary history. The inferred unrooted phylogenetic tree obtained by three different methods revealed the existence of four clearly separated clades (C-terminal and N-terminal domains of catalase-peroxidases, ascorbate peroxidases, and cytochrome c peroxidases) which were segregated early in the evolution of this superfamily. From the results, it is obvious that the duplication event in the gene for catalase-peroxidase occurred in the later phase of evolution, in which the individual specificities of the peroxidase families distinguished were already formed. Evidence is presented that class I of the heme peroxidase superfamily is spread among prokaryotes and eukaryotes, obeying the birth-and-death process of multigene family evolution.     

Richness,species composition and functional groups in Agaricomycetes communities along a vegetation and elevational gradient in the Andean Yungas of Argentina

The Neotropics are among the least explored regions from a mycological perspective. A few recent molecular studies in South America have shown high fungal diversity as well as numerous groups of mostly undescribed taxa. Through soil metabarcoding analysis we compared richness and species composition among macrofungal communities, belonging to Agaricales, Russulales, Boletales and Phallomycetidae groups, in three elevational forests types in the subtropical Yungas of Northwestern Argentina (Piedmont forest; Montane forest, Montane cloud forest). The aims of this study were to assess richness of taxonomic and functional groups along the elevation gradient and to assess the relationships between environmental variables and species composition in the studied fungal communities. The results have shown rich Agaricomycetes communities, diversely structured among forests habitats. The elevation gradient differentially affected the richness and distribution of Agaricales, Russulales, Boletales and Phallomycetidae. Based on fungal trophic modes and guilds, the gradient also affected the ectomycorrhizal taxa distribution. When considering the basidiomata growth forms (agaricoid, boletoid, gasteroid, etc.), only the secotioid type showed significant elevational differences. Additional analyses indicated that saprotrophic nutritional mode was dominant along the entire gradient, being partially replaced by biotrophic modes at higher elevations. Fungal communities in the Montane cloud forests are most dissimilar when compared with communities at the Piedmont forest and Montane forest, which is consistent with the different biogeographic origins of these forests. DNA metabarcoding sequence analysis provided detailed information on the diversity and taxonomic and functional composition of macrofungal communities.     

Delignification of Wood Chips and Pulps by Using Natural and Synthetic Porphyrins: Models of Fungal Decay

Kraft pulps, prepared from softwoods, and small chips of birch wood were treated with heme and tert-butyl hydroperoxide in aqueous solutions at reflux temperature. Analyses of treated pulps showed decreases in kappa number (a measure of lignin content) from about 36 to less than 2, with concomitant increases in brightness (80% increase in the better samples). Analyses of treated wood chips revealed selective delignification and removal of hemicelluloses. After 48 h of treatment, lignin losses from the wood chips approached 40%, and xylose/mannose (hemicellulose) losses approached 70%, while glucose (cellulose) losses were less than 10%. Examination of delignified chips by transmission electron microscopy showed that the removal of lignin occurred in a manner virtually indistinguishable from that seen after decay by white rot fungi. Various metalloporphyrins, which act as biomimetic catalysts, were compared to horseradish peroxidase and fungal manganese peroxidase in their abilities to oxidize syringaldazine in an organic solvent, dioxane. The metalloporphyrins and peroxidases behaved similarly, and it appeared that the activities of the peroxidases resulted from the extraction of heme into the organic phase, rather than from the activities of the enzymes themselves. We concluded that heme-tert-butyl hydroperoxide systems in the absence of a protein carrier mimic the decay of lignified tissues by white rot fungi.     

Evolution of structure and function of Class I peroxidases

The phylogenetics of Class I of the heme peroxidase-catalase superfamily currently representing over 940 known sequences in all available genomes of prokaryotes and eukaryotes has been analysed. The robust reconstructed tree for 193 Class I peroxidases with 6 selected Class II representatives reveals all main trends of molecular evolution. It suggests how the ancestral peroxidase gene might have been transferred from prokaryotic into eukaryotic genomes. Besides well known families of catalase-peroxidases, cytochrome c peroxidases and ascorbate peroxidases, the phylogenetic analysis shows for the first time the presence of two new well separated clades of hybrid-type peroxidases that might represent evolutionary bridges between catalase-peroxidases and cytochrome c peroxidases (type A) as well as between ascorbate peroxidases and Class II peroxidases (type B). Established structure-function relationships are summarized. Presented data give useful hints on the origin and evolution of catalytic promiscuity and specificity and will be a valuable basis for future functional analysis of Class I enzymes as well as for de novo design.     

Patterns of lignin degradation and oxidative enzyme secretion by different wood- and litter-colonizing basidiomycetes and ascomycetes grown on beech-wood

The degradation of lignocellulose and the secretion of extracellular oxidoreductases were investigated in beech-wood (Fagus sylvatica) microcosms using 11 representative fungi of four different ecophysiological and taxonomic groups causing: (1) classic white rot of wood (e.g. Phlebia radiata), (2) 'nonspecific' wood rot (e.g. Agrocybe aegerita), (3) white rot of leaf litter (Stropharia rugosoannulata) or (4) soft rot of wood (e.g. Xylaria polymorpha). All strong white rotters produced manganese-oxidizing peroxidases as the key enzymes of ligninolysis (75-2200 mU g(-1)), whereas lignin peroxidase activity was not detectable in the wood extracts. Interestingly, activities of two recently discovered peroxidases - aromatic peroxygenase and a manganese-independent peroxidase of the DyP-type - were detected in the culture extracts of A. aegerita (up to 125 mU g(-1)) and Auricularia auricula-judae (up to 400 mU g(-1)), respectively. The activity of classic peroxidases correlated to some extent with the removal of wood components (e.g. Klason lignin) and the release of small water-soluble fragments (0.5-1.0 kDa) characterized by aromatic constituents. In contrast, laccase activity correlated with the formation of high-molecular mass fragments (30-200 kDa). The differences observed in the degradation patterns allow to distinguish the rot types caused by basidiomycetes and ascomycetes and may be suitable for following the effects of oxidative key enzymes (ligninolytic peroxidases vs. laccases, role of novel peroxidases) during wood decay.     

Heterogeneity and regulation of manganese peroxidases from Phanerochaete chrysosporium.

Lignin and Mn peroxidases are two families of isozymes produced by the lignin-degrading fungus Phanerochaete chrysosporium under nutrient nitrogen or carbon limitation. We purified to homogeneity the three major Mn peroxidase isozymes, H3 (pI = 4.9), H4 (pI = 4.5), and H5 (pI = 4.2). Amino-terminal sequencing of these isozymes demonstrates that they are encoded by different genes. We also analyzed the regulation of these isozymes in carbon- and nitrogen-limited cultures and found not only that the lignin and Mn peroxidases are differentially regulated but also that differential regulation occurs within the Mn peroxidase isozyme family. The isozyme profile and the time at which each isozyme appears in secondary metabolism differ in both nitrogen- and carbon-limited cultures. Each isozyme also responded differently to the addition of a putative inducer, divalent Mn. The stability of the Mn peroxidases in carbon- and nitrogen-limited cultures was also characterized after cycloheximide addition. The Mn peroxidases are more stable in carbon-limited cultures than in nitrogen-limited cultures. They are also more stable than the lignin peroxidases. These data collectively suggest that the Mn peroxidase isozymes serve different functions in lignin biodegradation.     

Anthracenediethanol inhibits lignin degradation by Phanerochaete chrysosporium by competing for oxidation by lignin peroxidase,and not by trapping singlet oxygen

The biodegradation of anthracene-9, 10-diethanol by the ligninolytic fungus Phanerochaete chrysosporium, previously though to involve singlet oxygen, is shown to be catalyzed by lignin peroxidases. Veratryl alcohol stimulated the enzymatic degradation of anthracenediethanol, and anthracenediethanol inhibited enzymatic oxidation of veratryl alcohol. Competition for oxidation by lignin peroxidase is suggested as the mechanism of the inhibition of lignin biodegradation by anthracenediethanol and related anthracene derivatives.Abbreviations ADE anthracene-9,10-diethanol - AES anthracene-9,10-bisethanesulfonic acid - DHP dehydrogenative polymerizate - DMF N,N-dimethylformamide - EPX 9,10-endoperoxide of ADE - PMR proton magnetic resonance     

A first-principles model of early evolution: emergence of gene families, species, and preferred protein folds

In this work we develop a microscopic physical model of early evolution where phenotype—organism life expectancy—is directly related to genotype—the stability of its proteins in their native conformations—which can be determined exactly in the model. Simulating the model on a computer, we consistently observe the “Big Bang” scenario whereby exponential population growth ensues as soon as favorable sequence–structure combinations (precursors of stable proteins) are discovered. Upon that, random diversity of the structural space abruptly collapses into a small set of preferred proteins. We observe that protein folds remain stable and abundant in the population at timescales much greater than mutation or organism lifetime, and the distribution of the lifetimes of dominant folds in a population approximately follows a power law. The separation of evolutionary timescales between discovery of new folds and generation of new sequences gives rise to emergence of protein families and superfamilies whose sizes are power-law distributed, closely matching the same distributions for real proteins. On the population level we observe emergence of species—subpopulations that carry similar genomes. Further, we present a simple theory that relates stability of evolving proteins to the sizes of emerging genomes. Together, these results provide a microscopic first-principles picture of how first-gene families developed in the course of early evolution.     

Analysis and expression of the class III peroxidase large gene family in Arabidopsis thaliana

Higher plants possess a large set of the classical guaiacol peroxidases (class III peroxidases, E.C. 1.11.1.7). These enzymes have been implicated in a wide array of physiological processes such as H(2)O(2) detoxification, auxin catabolism and lignin biosynthesis and stress response (wounding, pathogen attack, etc.). During the last 10 years, molecular cloning has allowed the isolation and characterization of several genes encoding peroxidases in plants. The achievement of the large scale Arabidopsis genome sequencing, combined with the DNA complementary to RNA (cDNA) expressed sequence tags projects, provided the opportunity to draw up the first comprehensive list of peroxidases in a plant. By screening the available databases, we have identified 73 peroxidase genes throughout the Arabidopsis genome. The evolution of the peroxidase multigene family has been investigated by analyzing the gene structure (intron/exon) in correlation with the phylogenetic relationships between the isoperoxidases. An evolutionary pattern of extensive gene duplications can be inferred and is discussed. Using a cDNA array procedure, the expression pattern of 23 peroxidases was established in the different organs of the plant. All the tested peroxidases were expressed at various levels in roots, while several were also detected in stems, leaves and flowers. The specific functions of these genes remain to be determined.     

DyP-like peroxidases of the jelly fungus Auricularia auricula-judae oxidize nonphenolic lignin model compounds and high-redox potential dyes

The jelly fungus Auricularia auricula-judae produced an enzyme with manganese-independent peroxidase activity during growth on beech wood (∼300 U l⁻¹). The same enzymatic activity was detected and produced at larger scale in agitated cultures comprising of liquid, plant-based media (e.g. tomato juice suspensions) at levels up to 8,000 U l⁻¹. Two pure peroxidase forms (A. auricula-judae peroxidase (AjP I and AjP II) could be obtained from respective culture liquids by three chromatographic steps. Spectroscopic and electrophoretic analyses of the purified proteins revealed their heme and peroxidase nature. The N-terminal amino acid sequence of AjP matched well with sequences of fungal enzymes known as “dye-decolorizing peroxidases”. Homology was found to the N-termini of peroxidases from Marasmius scorodonius (up to 86%), Thanatephorus cucumeris (60%), and Termitomyces albuminosus (60%). Both enzyme forms catalyzed not only the conversion of typical peroxidase substrates such as 2,6-dimethoxyphenol and 2,2′-azino-bis(3-ethylthiazoline-6-sulfonate) but also the decolorization of the high-redox potential dyes Reactive Blue 5 and Reactive Black 5, whereas manganese(II) ions (Mn²⁺) were not oxidized. Most remarkable, however, is the finding that both AjPs oxidized nonphenolic lignin model compounds (veratryl alcohol; adlerol, a nonphenolic β-O-4 lignin model dimer) at low pH (maximum activity at pH 1.4), which indicates a certain ligninolytic activity of dye-decolorizing peroxidases.     

Prokaryotic origins of the non-animal peroxidase superfamily and organelle-mediated transmission to eukaryotes

Members of the superfamily of plant, fungal, and bacterial peroxidases are known to be present in a wide variety of living organisms. Extensive searching within sequencing projects identified organisms containing sequences of this superfamily. Class I peroxidases, cytochrome c peroxidase (CcP), ascorbate peroxidase (APx), and catalase peroxidase (CP), are known to be present in bacteria, fungi, and plants, but have now been found in various protists. CcP sequences were detected in most mitochondria-possessing organisms except for green plants, which possess only ascorbate peroxidases. APx sequences had previously been observed only in green plants but were also found in chloroplastic protists, which acquired chloroplasts by secondary endosymbiosis. CP sequences that are known to be present in prokaryotes and in Ascomycetes were also detected in some Basidiomycetes and occasionally in some protists. Class II peroxidases are involved in lignin biodegradation and are found only in the Homobasidiomycetes. In fact class II peroxidases were identified in only three orders, although degenerate forms were found in different Pezizomycota orders. Class III peroxidases are specific for higher plants, and their evolution is thought to be related to the emergence of the land plants. We have found, however, that class III peroxidases are present in some green algae, which predate land colonization. The presence of peroxidases in all major phyla (except vertebrates) makes them powerful marker genes for understanding the early evolutionary events that led to the appearance of the ancestors of each eukaryotic group.