# Fungal dye-decolorizing peroxidase diversity: roles in either intra- or extracellular processes

Martino Adamo, Enrico Büttner, Harald Kellner, Patricia Luis, Sophie Comtet-Marre, Laurent Vallon, Rocio Prego, Martin Hofrichter, Mariangela Girlanda, Pierre Peyret, et al.

► To cite this version:

Martino Adamo, Enrico Büttner, Harald Kellner, Patricia Luis, Sophie Comtet-Marre, et al. Fungal dye-decolorizing peroxidase diversity: roles in either intra- or extracellular processes. Applied Microbiology and Biotechnology, 2022, 106 (8), pp.2993-3007. 10.1007/s00253-022-11923-0. hal-03650390

HAL Id: hal-03650390

https://hal.inrae.fr/hal-03650390

Submitted on 25 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Fungal dye-decolorizing peroxidase diversity: roles in either intra- or extracellular processes

Martino Adamo<fnref n="1" />,<fnref n="2" /> ID · Sophie Comtet-Marre<fnref n="3" /> · Enrico Büttner<fnref n="4" /> · Harald Kellner<fnref n="4" /> · Patricia Luis<fnref n="2" /> · Laurent Vallon<fnref n="2" /> · Rocio Prego<fnref n="3" /> · Martin Hofrichter<fnref n="4" /> · Mariangela Girlanda<fnref n="1" /> · Pierre Peyret<fnref n="3" /> · Roland Marmeisse<fnref n="1" />,<fnref n="2" />,<fnref n="5" />,<fnref n="6" />

Received: 28 February 2022 / Revised: 6 April 2022 / Accepted: 8 April 2022
© The Author(s) 2022

## Abstract

Fungal dye-decolorizing peroxidases (DyPs) have found applications in the treatment of dye-contaminated industrial wastes or to improve biomass digestibility. Their roles in fungal biology are uncertain, although it has been repeatedly suggested that they could participate in lignin degradation and/or modification. Using a comprehensive set of 162 fully sequenced fungal species, we defined seven distinct fungal DyP clades on basis of a sequence similarity network. Sequences from one of these clades clearly diverged from all others, having on average the lower isoelectric points and hydropathy indices, the highest number of N-glycosylation sites, and N-terminal sequence peptides for secretion. Putative proteins from this clade are absent from brown-rot and ectomycorrhizal species that have lost the capability of degrading lignin enzymatically. They are almost exclusively present in white-rot and other saprotrophic *Basidiomycota* that digest lignin enzymatically, thus lending support for a specific role of DyPs from this clade in biochemical lignin modification. Additional nearly full-length fungal DyP genes were isolated from the environment by sequence capture by hybridization; they all belonged to the clade of the presumably secreted DyPs and to another related clade. We suggest focusing our attention on the presumably intracellular DyPs from the other clades, which have not been characterized thus far and could represent enzyme proteins with novel catalytic properties.

## Key points

* A fungal DyP phylogeny delineates seven main sequence clades.

* Putative extracellular DyPs form a single clade of *Basidiomycota* sequences.

* Extracellular DyPs are associated to white-rot fungi.

**Keywords** DyP-type peroxidase · Lignin degradation · Fungi · Sequence capture

Martino Adamo and Sophie Comtet-Marre are Equal contributors

<figure><img src="image_2.png" /></figure>

<footnote>Martino Adamo
martino.adamo@unito.it</footnote>

<footnote n="1">¹ Department of Life Sciences and Systems Biology (DBIOS), Università degli Studi Di Torino, 25 Viale P.A. Mattioli, 10125 Torino, Italy</footnote>

<footnote n="2">² Univ Lyon, Université Claude Bernard Lyon 1, CNRS, INRAE, UMR Ecologie Microbienne, VetAgro Sup43 Boulevard du 11 Novembre 1918, 69622 Villeurbanne Cedex, France</footnote>

<footnote n="3">³ Université Clermont Auvergne, INRAE, MEDISt, 63000 Clermont-Ferrand, France</footnote>

<footnote n="4">⁴ Department of Bio- and Environmental Sciences, International Institute Zittau, Technische Universität Dresden, Zittau, Germany</footnote>

<footnote n="5">⁵ Institut de Systématique, Evolution, Biodiversité (ISYEB), Muséum National d'Histoire Naturelle, CNRS, Sorbonne Université, EPHE, Université des Antilles, CP39, 57 rue Cuvier, 75005 Paris, France</footnote>

<footnote n="6">⁶ Institute for Sustainable Plant Protection (IPSP), National Research Council (CNR), 25 Viale P.A. Mattioli, 10125 Torino, Italy</footnote>
**Introduction**

Lignin is a complex plant cell-wall polymer made of differ-
ent cross-linked phenylpropanoid monomers (Schmidt 2006)
that represents about 30% of the plant-derived organic carbon
available in the biosphere (Boerjan et al. 2003). Besides pro-
viding structural rigidity in living plants, another important
function is to protect cell wall polysaccharides from degra-
dation by hydrolytic enzymes secreted by microorganisms
(lignin barrier). Control of lignin degradation represents an
important field of research either for the use of this polymer
as a source of aromatic molecules or simply for its removal
from plant biomass to access cell-wall polysaccharides for
biorefinery and biofuel production (Duval and Lawoko 2014;
Chio et al. 2019).

White-rot (WR) fungi, a specific functional group of
saprotrophic *Basidiomycota*, are considered the main group
of microorganisms able to perform enzymatic ligninolysis.
Phylogenomic studies support the assumption that the WR
phenotype is an ancestral character that arose once early in
the evolution of the *Basidiomycota* (Floudas et al. 2012).
Subsequently, the WR phenotype was lost several times
independently during the course of evolution to give rise
to different lineages of either brown-rot (BR) saprotrophic
or mutualistic ectomycorrhizal (EM) *Basidiomycota* that
associate to plant roots (Kohler et al. 2015; Martin et al.
2016; Nagy et al. 2017; Miyauchi et al. 2020). Both BR
and EM fungi have lost the capacity of performing enzy-
matic ligninolysis, although being capable of modifying the
lignin structure to some extent using oxidizing mechanisms
(Arantes et al. 2012; Rineau et al. 2012).

From a mechanistic point of view, biochemical and genetic
studies have demonstrated that several secreted enzymes
belonging to the so-called heme-containing, fungal-specific
class II peroxidases (E.C. 1.11.1.-) participate directly or
indirectly in lignin decomposition (Hammel and Cullen 2008;
Salame et al. 2014; Ayuso-Fernández et al. 2018). Among
these class II peroxidases, lignin (LiP, E.C. 1.11.1.14) and
versatile (VP, E.C. 1.11.1.16) peroxidases interact directly
with lignin to cleave non-phenolic ether bonds that link phe-
nylpropanoid monomers. A further class II enzyme, manga-
nese peroxidase (MnP, E.C. 1.11.1.13), but also VPs, oxidize
Mn²⁺ to Mn³⁺, which is then chelated with organic acids
and cleaves as low-molecular mass and diffusible redox-
mediator phenolic structures of lignin (Ruiz-Dueñas et al.
2021). Besides biochemical studies, the essential role of class
II peroxidases in lignin decomposition was also supported by
phylogenomic studies. Diversification of this protein family
correlates with diversification of WR fungi and the several
transitions from a WR to BR or EM trophic modes almost
always went along with a loss or a sharp contraction of class
II peroxidase gene family in the fungal genomes (Kohler

et al. 2015; Martin et al. 2016; Nagy et al. 2017; Miyauchi
et al. 2020; Ruiz-Dueñas et al. 2021).

Despite the undoubted contribution of class II peroxi-
dases to lignin decomposition, their exclusive role in this
process is debated. Indeed, substantial lignin degradation
seemingly occurs in the absence of organisms producing
these enzymes, as in the case of the termite gut (Li et al.
2017). On the other hand, other microbial enzymes have
been shown to be able to cleave phenolic and/or non-
phenolic bonds of artificial lignin substrates (e.g. β-O-4-
dimers = adlerol) and/or oxidize Mn²⁺ to Mn³⁺ and there-
fore behave similarly to either class II LiP/VPs or MnPs
(Fernández-Fueyo et al. 2015). One such enzyme family
is that of “dye-decolorizing peroxidase” (DyP) (pfam no.
PF04261, E.C. 1.11.1.19) (Celis and Dubois 2015). These
heme-containing peroxidases are unrelated to class II peroxi-
dases and were identified for their unique ability to oxidize
high-redox potential anthraquinone dyes (Kim et al. 1995),
a property which makes them good candidates for the reme-
diation of dye-polluted waste-water (Scheibner et al. 2008).
DyPs present a ferredoxin-like fold (Singh and Eltis 2015),
which is unique among peroxidases, and are present in all
three domains of life (Bacteria, Archaea, and Eukarya) as
opposed to class II peroxidases that are restricted to the
Eumycota (true fungi) (Zámocký et al. 2012; Celis and
Dubois 2015). Global phylogenies of the DyP (sub)family
show that, with very few exceptions, fungal sequences are
monophyletic and form a well-supported clade that does not
include sequences from other taxonomic groups (Sugano
2009; Zámocký et al. 2015).

Different levels of biochemical information exist for at least
13 fungal DyPs (12 from *Basidiomycota*, one from *Ascomy-
cota*), including 3D structures (based on protein crystals) for
two of them (Yoshida et al. 2011; Strittmatter et al. 2013),
identification of amino acids participating to catalysis through
site-directed mutagenesis and biophysical measurements
(Linde et al. 2015b) as well as reconstruction and characteri-
zation of ancestral fungal sequences (Zitare et al. 2021).

DyPs can form surface protein radicals (surface-exposed
tryptophan radicals) via a long-range electron transfer
(LRET) and are hereby able to oxidize phenolic compounds
and synthetic dyes. In *Pleurotus ostreatus* PosDyP4, a
Mn²⁺-oxidation activity located on a different surface site,
surrounded by acidic amino acids, was identified (Fernán-
dez-Fueyo et al. 2015). This activity is analogous to that
accomplished by MnPs and VPs (Fernández-Fueyo et al.
2018), and Krahe et al. (2020) were able to use a closely
related DyP of *Pleurotus sapidus* with Mn²⁺-oxidation activ-
ity to cleave alkenes (aryl alkenes (E)-methyl isoeugenol,
α-methylstyrene, and *trans*-anethole) in the presence of
Mn²⁺, which could be relevant in flavor production. A fur-
ther Mn²⁺-oxidizing DyP activity has been recently reported
for the first characterized ascomycetous DyP from *Xylaria grammica* (Kimani et al. 2021).

The roles of DyPs in fungal biology have never been addressed by forward genetic approaches, and they have been proposed to contribute to the detoxification of naturally occurring soil-borne or plant metabolites (Valette et al. 2017), including natural anthraquinone derivatives (Sugawara et al. 2019). They have also repeatedly been suspected to participate in lignin decomposition and in the oxidation of lignin-related compounds (Liers et al. 2010; Celis and Dubois 2015; Cragg et al. 2015; de Gonzalo et al. 2016; Catucci et al. 2020; de Eugenio et al. 2021), despite the fact that fungal class II peroxidases perform far better on "lignin substrates" than bacterial or fungal DyPs (Linde et al. 2021). In support of their putative roles in lignin oxidation, fungal DyPs have been also often reported as secreted enzymes, i.e. extracellularly acting biocatalysts. Several of them were purified from fungal culture filtrates, and they have been identified as well in the secretomes of fungi growing on lignin-rich substrates (Salvachúa et al. 2016). From an evolutionary perspective, as for class II peroxidases and other gene families, the fungal DyP (sub)family diversification parallels the evolution of WR fungi and regresses in BR and EM lineages (Nagy et al. 2017).

In the absence of forward genetic data and in the presence of uncertain biochemical activities, the objective of the present study has been to further evaluate the potential biological functions of fungal DyPs, including their putative contribution to lignin decomposition using phylogenetic and correlative approaches. Two strategies were followed. One was to establish a comprehensive view on the evolution of the DyP (sub)family in the fungal kingdom, using available genomic data, and to confront these data to our current knowledge on DyPs' catalytic mechanisms and evolution of trophic modes in fungi. A second strategy was to isolate new fungal DyP genes expressed in different environments (soils and decomposing wood) in an attempt to characterize potentially divergent DyP sequences distantly related to sequences present in public databases. The latter approach could potentially identify new DyPs with potentially useful catalytic properties for novel applications in biotechnology.

# Materials and methods

## Species and sequence datasets

Publicly available DyP sequences were extracted from GenBank and the JGI Mycocosm portal (https://genome.jgi.doe.gov/programs/fungi/index.jsf; Grigoriev et al. 2014) databases as well as from our own unpublished whole-genome sequences. DyP protein sequences were identified in public databases using BlastP searches (Altschul et al. 1990) using as queries the protein sequences of two biochemically

characterized DyPs from *Auricularia auricula-judae* (accession No. JQ650250) and *Bjerkandera adusta* (CDN40127.1). A threshold E-value of $10^{-5}$ was retained for the inclusion of the sequences in the dataset. DyP sequences extraction from the unpublished genomes was performed using HMMER (http://hmmer.org/). A HMMER profile was first created starting from a Clustal Omega alignment (Sievers and Higgins 2014) of 150 DyP sequences from the public databases. Phylogenetically-related bacterial DyP sequences, all belonging to the DyP clade C (Zámocký et al. 2015), were identified following blastp searches using as queries different sequences representative of the phylogenetic diversity of the fungal DyPs.

## Phylogenetic and network analyses

Only putative full-length protein sequences, starting with methionine and ending with a stop codon, corresponding to CDS longer than 1000 bp and without gaps in the regions where amino acids essential for catalysis have been identified (as described in Linde et al. 2015a), were used for a phylogenetic reconstruction of the fungal DyP (sub)family. In the case of the *Aspergillus* and *Penicillium*, genera for which the genomes of numerous species are available, only sequences from *Aspergillus oryzae* and *Penicillium italicum* were included in the phylogenetic analysis. Sequences were aligned using MUSCLE (Edgar 2004). ModelFinder (Kalyaanamoorthy et al. 2017), as implemented on the IQ-Tree web Server (http://iqtree.cibiv.univie.ac.at/; (Trifinopoulos et al. 2016)), was used to find the best substitution model (LG+R). IQ-TREE Web Server was also used to perform a maximum likelihood phylogenetic analysis (Nguyen et al. 2015). One thousand bootstrap replicates were performed to assess tree topology.

A fungal species phylogeny that encompassed all examined fully-sequenced basidiomycetous species was obtained by modifying a publicly available phylogenetic tree (#Tr106951, Zhao et al. 2017), built using an alignment of six genes, downloaded from TreeBase (www.treebase.org/). Species whose genomes have not been sequenced were removed from the original six-genes alignment while the corresponding six genes from newly sequenced species were added to it to construct a new species tree following the approach presented in Zhao et al. (2017). Species taxonomy was assigned using the MycoBank database (www.mycobank.org). Both gene and species trees were edited using the online software iTOL (https://itol.embl.de; (Letunic and Bork 2016)). We inferred gene loss and duplication along the species phylogenetic tree using a comparative genomics method with NOTUNG 2.9 as described in the original article of Chen et al. (2000).

To visualize similarities between fungal DyP protein sequences, a sequence similarity network (SSN) was
computed using the EFI-Enzyme Similarity Tool (https://efi.igb.illinois.edu/efi-est/; (Gerlt et al. 2015)) starting from the complete set of DyP protein sequences. The lowest pairwise alignment score limit for the output file was set to 66. The full-network output file was visualized with Cytoscape 3.5.1 (Shannon et al. 2003). The network was represented with an “organic” layout, and only edges with a pairwise alignment score above 110 were visualized. To further evaluate the potential contribution of agaricomycetous DyP genes to organic matter (e.g. wood) and more specifically, to lignin degradation, we compared the prevalence of the different DyP clades in polyphyletic species groups of *Agaricomycota*, which differed with respect to their trophic modes (“trophic guilds”). Using the FUNguild (Nguyen et al. 2016) and experts’ opinions, we defined four trophic modes. “White-rot” (WR), sensu stricto encompasses saprotrophic or facultative pathogenic, lignicolous species (living on wood and coarse woody debris), which degrade lignin enzymatically to produce a “white-rot” phenotype in compact wood (where finally whitish cellulose fibers prevail). All of these species possess mostly numerous class-II peroxidase genes (Mn, versatile, and/or lignin peroxidases) in their genomes (Floudas et al. 2012). “Brown-rot” (BR) species, although also lignicolous, produce a brown-rot phenotype that results from a selective degradation of cell-wall polysaccharides (cellulose and xylan/hemicelluloses) along with a nonenzymatic modification of the lignin polymer. “Ectomyccorrhizal” (EM) species are mutualistic symbionts of mainly trees (larger woody plants) that have (almost) lost the ability to degrade plant cell walls (Shah et al. 2016). Finally, a fourth group named “other saprotrophs” (OS) encompasses all other species that are either litter-decomposing saprotrophs (e.g. *Agaricus bisporus*), coprophilic species (e.g. *Coprinopsis cinerea*), or endophytic ones (e.g. *Sebacina vermifera*). Since these OS species do not grow on compact wood, their modes of degrading the different plant polymers do not exactly fit into the white-rot or brown-rot definitions, although several authors classify several of these species as primarily white-rotters (Eastwood et al. 2011; Nagy et al. 2016).

## Prediction of sequence features

Prediction of N-terminal signal peptides for eukaryotic secreted proteins was performed with SignalP (https://services.healthtech.dtu.dk/service.php?SignalP-5.0; (Nielsen, 2017)) with default settings except for the D-cutoff value, which was set at 0.34 for both the noTM (transmembrane domain) and the TM networks. SignalP raw results were reported on the protein phylogenetic tree. The number of putative N-glycosylation sites was predicted using NetNGlyc (https://services.healthtech.dtu.dk/service.php?NetNGlyc-1.0) (Gupta et al. 2002). Other protein features, such

as protein length, absolute numbers of specific amino acids, predicted protein isoelectric point (pI), and global protein hydropathy using the “grand average of hydropathy” (GRAVY), were predicted using tools of the sequence manipulation suite (https://www.bioinformatics.org/sms2/; (Stothard 2000)).

Principal component analysis (PCA) was performed to cluster DyP protein sequences according to their protein features in the R environment using the R 3.6.3 (R Core Team 2019). Principal components were extracted with the prcomp function of the R base package and visualized with the ggbiplot package v0.55. Kmeans analysis was performed with the kmeans function of the R base package. Variable’s (i.e. protein features) contribution to principal components (PC) was determined using the fviz_contrib function from factoextra v1.0.6 R package. Statistical differences were tested using the nonparametric Kruskal-Wallis test and Dunn’s post hoc tests with dunn.test package v1.3.4 (Dinno 2015). DyP abundance comparisons among different fungal guilds were performed with the Yates-corrected Chi-square test.

## Probes design for gene capture by hybridization

From 1267 publicly available fungal DyP DNA-coding sequences, we designed 69 (70 bp long) degenerated probes (Supplemental Table S1) with the KASpOD software (https://g2im.u-clermont1.fr/kaspod/index.php, (Parisot et al. 2012)). In silico probe coverage evaluation, on the 1267 fungal DyP sequences detects 94.2% of the ascomycetous sequences and 81.5% of the basidiomycetous ones (with four allowed mismatches). Oligonucleotides corresponding to probe sequences were synthesized with two flanking adaptor sequences for their PCR amplification and conversion to biotinylated RNA probes using the T7 RNA polymerase, as described in Bragalini et al. (2014).

## Environmental samples, RNA extraction, and cDNA synthesis

Grassland and forest soil, as well as deadwood samples from six different contrasted geographic sites in Italy and France, were included in this study (Supplemental Table S2). Four of them had already been described in Adamo et al. (2020) and one in Bragalini et al. (2014).

RNA from soil samples was extracted from 2 g of material using the RNA Power Soil extraction kit from MOBIO Laboratories (Carlsbad, CA, USA) according to the manufacturer’s instructions. Deadwood RNA was extracted from 100 mg of decaying wood following the protocol described in Adamo et al. (2017). Eukaryotic cDNAs were
synthesized from 2 µg of total environmental RNA using the “Mint-2 cDNA synthesis and amplification kit” according to the manufacturer’s instructions (Evrogen, Moskow, Russian Federation). First-strand synthesis was initiated at the mRNA poly-A end using a modified poly-dT primer (CDS-4 M). All resulting cDNAs were bordered at their 5’ and 3’ extremities by the same M1 sequence (5’-AAGCAG TGGTATCAACGCAGAGT-3’) that can be used for priming cDNA PCR amplification.

## Gene capture by hybridization

DyP cDNA sequence capture was performed as already described in Denonfoux et al. (2013) and Bragalini et al. (2014). Briefly, 2 µg of heat-denatured cDNA were hybridized to 500 ng of a mix of biotinylated RNA probes for 24 h at 65 °C in microcentrifuge tubes. The probes/DyP cDNA hybrids were then captured on streptavidin-coated paramagnetic beads (Dynabeads M-280 Streptavidin, Invitrogen, Waltham, MA USA). After several washing steps to eliminate unbound cDNAs, the captured cDNAs were detached from the beads using 0.1 M NaOH and purified using the “Qiaquick PCR purification kit” (Qiagen, Hilden, Germany). Captured cDNAs were further amplified using “Encyclo Plus PCR kit” and the M1 primer and were purified using the “Qiaquick PCR purification kit.” Captured cDNA was then subjected to the second round of capture identical to the first one. The enrichment in DyP sequences along the capture protocol was evaluated via a semi-quantitative PCR using DyP-specific primers (Bragalini et al. 2014; Kellner et al. 2014).

## Sequencing of captured cDNAs and bioinformatics analysis of the sequencing data

Since known fungal DyP CDSs range in size from 1500 to 2000 bp, captured cDNA was the first size fractionated to sequence only putatively full-length cDNAs and to eliminate truncated and contaminating sequences. Fractionation of 1 µg amplified cDNA (between 1400 and 2000 bp) was performed using a “BluePippin” instrument (Sage Science, Beverly, MA, USA). Fractionated cDNA was quantified by fluorimetry using the “Qubit dsDNA HS Assay Kit” and “Qubit Fluorimeter 2.0” (Thermo Fisher Scientific, Waltham, MA USA) and re-amplified by PCR using the M1 primer and “KAPA HiFi Taq polymerase” (KAPA Biosystems, Wilmington, MA USA). The 50 µl PCR mix contained 10 ng of cDNA, 10 µl of 5X KAPA HiFi Fidelity Buffer, 1.5 µl of 10 mM dNTPs, 10 µl of the 10 µM concentrated M1 primer, and 1 U of KAPA HiFi Taq polymerase. After an initial denaturation at 94 °C for 3 min, cDNA fragments were amplified for 25 cycles

comprising 20 s at 98 °C, 15 s at 66 °C, and 1 min at 72 °C. Amplification was terminated by a final elongation at 72 °C for 3 min. The “Agencourt® AMPure® XP kit” (Beckman Coulter, Brea, CA USA) was used to perform a final purification and partial size fractionation of the captured cDNA. A 0.5X concentration of beads was used in order to remove small unspecific fragments from the eluted material. Purified cDNA quantity was assessed using the “Qubit dsDNA HS Assay Kit” and “Qubit Fluorimeter 2.0”, while DNA purity was assessed by measuring the OD 260/OD 280 and OD 260/OD 230 ratios using a NanoDrop TM 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA USA).

All captured cDNA samples were sequenced using the Illumina HiSeq 2000 2×250 bp technology (I.G.A. Technologies, Udine, Italy) using the Illumina MiSeq 2×250 bp standard protocol for samples preparation.

Adapter sequences were eliminated using Cutadapt (Martin 2011). Sequence quality was evaluated with Trimmomatic (Bolger et al. 2014); bases with quality lower than 20 were eliminated, and only sequences longer than 60 bases were kept. Trimmed sequences were assembled using IDBA-UD (Peng et al. 2012) and the default parameter. To obtain longer contigs, the resulting contigs were then assembled using CAP3 and default parameters (Huang 1999). For the identification of DyP sequences, a similarity search between known fungal DyP sequences and contigs was performed using DIAMOND v0.9.22.123 (Buchfink et al. 2015), with the BLASTx command in “sensitive” mode (i.e. a maximal e-value of $1.\times10^{-5}$ and a minimal identity of 50%). Matching cDNA sequences were translated using the Expasy (Gasteiger 2003) translate tool and further analyzed for the presence of the DyP peroxidase domain using ScanProSite and the Prosite database (de Castro et al. 2006; Sigrist et al. 2012). Identified DyP proteins were clustered using a 90% identity threshold using CD-HIT v4.7 (Fu et al. 2012). Sequences are available on NCBI with the accession numbers OM674899-OM674912.

## Results

### Global distribution of DyP genes in the fungal kingdom (Eumycota)

Among the 227 fungal genomes inspected, 454 DyP-coding genes were identified in 162 species belonging exclusively to the *Dikarya* (*Ascomycota* and *Basidiomycota*; Supplemental Table S3) and not in any other fungal taxonomic groups. In the *Ascomycota*, 134 DyP genes were identified in the genomes of 84 species, all from the *Pezizomycotina* and none from the *Saccharomycotina* and *Taphrinomycotina*. For *Basidiomycota*, DyP genes were absent from the
early-branching clades of *Ustilagomycotina* (subphylum of smut fungi), and present in the *Pucciniomycotina* (rust fungi) and *Agaricomycotina*. In these latter groups, 306 sequences were identified in the genomes of 78 species (34% of inspected species). Predicted numbers of putative DyP-encoding genes per genome differed greatly, ranging from 0 to 1 up to 64 genes in the case of the basidiomycetous species *Sphaerobolus stellatus* (which is remarkable since this fungus contains also 78 genes of unspecific peroxygenases/UPOs, EC 1.11.2.1; Hofrichter et al. 2020). In *Ascomycota*, the number of predicted DyP genes never exceeded six per genome (e.g. *Paraconiothyrium sporulosum*). A detailed list of DyP genes identified in the inspected species is given in Supplemental Table S4.

## Evolution of fungal DyP proteins

DyP amino acid sequences were used to compute a sequence similarity network, which individualized seven well-separated clusters of sequences (Fig. 1) using a pairwise alignment score of 110. Following a phylogenetic analysis, each cluster defined one single monophyletic clade, with the exception of cluster VI-1, from which emerged clade VI-2 (Fig. 2). From the taxonomic point of view, while clade VI-1 encompassed sequences from both *Basidiomycota* and *Ascomycota*, other clades had a narrower distribution pattern as illustrated for clade II present in the plant pathogens belonging to the class *Pucciniomycotina* (with the exception of one sequence originating from the *Agaricales Termitomyces* sp.).

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Sequence-similarity clustering of DyP protein sequences retrieved from fully sequenced fungal genomes and of biochemically-characterized fungal DyPs. Clustering using the EFI-EST server identified seven separate clusters (I–VI-1 and VI-2). Symbol shape indicates the presence (diamond) or absence (circle) of a predicted signal peptide. Symbol color indicates the trophic mode of the corresponding fungal species; red, white rot; green, brown rot or ectomycorrhizal; grey, other. Symbols with a yellow contour indicate biochemically-characterized enzymes</figcaption></figure>

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Phylogenetic analysis of the fungal DyP peroxidase family based on the alignment of 301 protein sequences retrieved from fully-sequenced genomes, corresponding to biochemically characterized enzymes, or encoded by environmental cDNAs isolated by gene capture by hybridization. Branch colors correspond to each of the sequence clusters illustrated in Fig. 1. The clade with black branches corresponds to bacterial sequences with the highest similarity values to fungal protein sequences. Grey dots indicate the bootstrap support (1000 replicates) of the corresponding branch; smallest dots, 75%; largest ones, 100%. The inner circle, the taxonomic origin of the sequences; blue, *Basidiomycota*; red, *Ascomycota*; grey, captured environmental sequences. Red stars, the presence of a putative signal peptide; red triangles, biochemically-characterized enzymes. Outer circle (orange bars), theoretical hydropathy (GRAVY) indices of the polypeptides (mostly negatives, a few positive values in cluster VI-2). Inner circle (light blue bars), pI values (expressed as pI-7)</figcaption></figure>

When representative sequences of each of the seven fungal DyP clades were individually used in blast searches to identify their closest relatives among bacterial DyP sequences, the same set of bacterial sequences was identified. These bacterial sequences all belonged to group C of DyP sequences as defined by (Zámocký et al. 2015). A selection of these bacterial sequences was included in the fungal DyP phylogeny. They all clustered in a single basal clade (Fig. 2).

For each fungal protein sequence, we predicted a number of features that either reflected the global protein structure (GRAVY hydropathy index, isoelectric point (pI), sequence length, number of putative N-glycosylation sites, number of acidic amino acid residues at the putative Mn²⁺ oxidation site), which could be related to catalytic activities as in the case of the number of aromatic residues (tyrosine, tryptophan, and phenylalanine), which can participate

in long-range electron transfers (LRETs) (Supplemental Table S4). When incorporated in a PCA analysis, these different variables separated several of the proteins according to the clades they belonged to (Fig. 3). The first two axes of the PCA altogether explained 52% of the variation in the dataset. Clade VI-2 was particularly well separated from all other clades with significantly lower average pI values and significantly higher GRAVY indices (Fig. 2), as well as the number of putative N-glycosylation sites per polypeptide (Dunn's post hoc test, *p* < 0.05) (Supplemental Fig. S1). Furthermore, 69% of the proteins affiliated with clade VI-2 also possessed a predicted N-terminal signal peptide for secretion (Fig. 2). Outside of clade VI-2, only four proteins from *Puccinia* spp., belonging to clade II, also processed a putative signal peptide. Interestingly, 13 of the 16 purified and characterized fungal DyPs identified in the literature belong to clade VI-2; the three others were affiliated with clade
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Principal component analysis (PCA) separates fungal DyP protein sequences (dots) belonging to different sequence clusters (same color code as in Fig. 2) according to their structural characteristics (vectors: hydropathy index (GRAVY), theoretical isoelectric point (pI), number of predicted N-glycosylation sites (N-glyc), polypeptide length (aa), number of specific amino acid residues (acidic, F, W, Y))</figcaption></figure>

VI-1. Most of these characterized proteins indeed appeared to be extracellular as they were purified from fungal culture filtrates presenting a dye-decolorizing activity (toward Reactive Blue5).

Inspection of amino acid residues present at positions essential for catalysis in the AauDyPI enzyme from *A. auricula-judae* (Strittmatter et al. 2013; Linde et al. 2015b) highlighted differences between the different clades (Supplemental Fig. S2). Positions in the distal side of the heme pocket are well conserved across all the clades, except in the case of Clade III, where the aspartate residue is replaced by a glycine one (position 168 in AauDyPI), and phenylalanine frequently replaced by valine (position 359). In the case of the proximal side of the heme pocket, while the histidine and aspartate residues, implicated in Fe³⁺ ligation as well as in peroxide-binding and cleavage, are well conserved across the clades, the central valine (position 253) residue present in the *AauDyPI* protein is frequently replaced by isoleucine, methionine, or phenylalanine residues, the latter two being predominant in clades I and IV, and III, respectively (Supplemental Fig. S2). Finally, the radical-forming residues tyrosine and tryptophan, probably involved in LRETs, present at positions 337 and 377 of the AauDyPI protein are conserved in all clades, except in clades VI-1 and VI-2, where they are frequently replaced by phenylalanine, another aromatic amino acid that can participate in LRETs (Acebes et al. 2017).

## DyPs and the evolution of trophic modes in the Agaricomycotina

We constructed a species phylogenetic tree that included 128 *Agaricomycotina* species and two species belonging to *Ustilaginomycotina* (as an outgroup) with fully sequenced genomes (Fig. 4), based on the tree published by Zhao et al. (2017). According to Zhao et al. (2017), the topology of this tree reflected the time-dependent divergence between the main agaricomycetous orders. Mapping the occurrence of DyP genes and clusters on this consensus species tree showed that most major agaricomycetous DyP clades (II, III, V, and VI, see Fig. 2) diverged early during the evolution of this taxon, being already present in all the most basal taxa of the *Agaricomycotina*, such as *Cantharellales*, *Sebacinales*, *Auriculariales*, and *Tremellales* (Hibbett et al. 2007). The tree also illustrates the recurrent loss, but also duplication, of several of the DyP clades with a differential distribution among the main agaricomycetous orders. For example, current genomic data suggest the loss of DyP genes belonging to clades VI-1 and VI-2 in the ancestor of the *Boletales*, which seemingly only possess members of clade V.

DyP phylogeny did not match fungal species phylogeny. One single species could host genes belonging to up to three clades (as in the case of the agaric species *Gymnopus luxurians* or the polypore *B. adusta* that both possess DyPs from clades VI-1, V, and III), and the presence
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Phylogenetic tree of 128 fully-sequenced *Agaricomycotina* species and two *Ustilaginomicotina* (used as outgroup) computed based on a six genes sequence alignment according to Zhao et al. (2017). Species names are colored according to the order they belong to. Colored dots associated with species names indicate their trophic mode, as reported in the legend. The number of predicted</figcaption></figure>

DyPs in each sequenced genome and their distribution in the different sequence clusters (Fig. 1 and 2) are presented in the bar chart. Arrows pointing downwards or upwards indicate significant expansion or reduction in the number of DyP genes affiliated to clusters VI-1 and VI-2 at the leaves of the tree (according to NOTUNG).

of genes encoding DyPs representative of one clade never excluded the presence of a gene of another clade. We predicted the number of gene loss or gene gain events, using the NOTUNG software, separately for the most

sequence-rich clades VI-1 and VI-2. Diversification at the level of the tree leaves (species-level) of clades VI-1 and VI-2 was not systematically coordinated between DyP clades (Fig. 4).
To further evaluate the potential contribution of agarico-
mycetous DyP genes to organic matter conversion, and more
specifically, to lignin degradation, we compared the prevalence
of the different DyP clades in *Agaricomycotina* with
respect to their trophic modes. The first observation from the
cluster analysis is that clade VI-2 (putatively secreted extracellular DyPs) encompassed sequences exclusively from WR
and OS species and none from BR or EM ones (Fig. 1).
This distribution is specific to clade VI-2. Indeed, all other clades with agaricomycetous sequences (clades II, III, V, and VI-1) contained sequences from WR, OS, BR, and EM species (the latter being absent from clade III). Although clade VI-2 sequences are absent from nearly two-thirds of the species classified as WR, and more WR species have been sequenced compared to BR and EM ones (Fig. 4), Chi² tests strongly support (p < 0.01) the higher prevalence of clade VI-2 sequences in WR (or WR and OS) species than in BR and EM ones. On the contrary, we did not find evidence for a higher prevalence of other clades in groups of species presenting a specific trophic mode (WR, OS, BR, or EM; Chi² test, p < 0.5).

Diversity of fungal DyPs expressed
in the environment

A set of 69 biotinylated degenerated RNA probes, whose
sequences spanned the known phylogenetic diversity of
fungal DyP sequences, was designed and used to capture
homologous cDNAs synthesized from mRNA directly
extracted from different environmental matrices (decaying
wood, forest and grassland leaf-litter, and soils). DNA recov-
ered after two successive rounds of capture (enrichment)
was randomly sequenced, and the reads were successfully
assembled to recover 14 full-length or nearly full-length DyP
cDNAs. All of these environmental sequences encode dif-
ferent DyP not present in the database. They were clearly
affiliated with either clade VI-2 (5 sequences) or clade VI-1
(9 sequences) (Fig. 2). In clade VI-1, three environmental
DyP sequences were related to reference sequences from
*Ascomycota* and six to sequences from *Basidiomycota*. We
affiliated “cluster 4” to *Basidiomycota*, although also close
to *Ascomycota* sequences in the phylogeny. In clade VI-2,
three of the five protein sequences had a predicted N-termi-
nal signal peptide.

Discussion

We provide an exhaustive and comprehensive view of DyP-
encoding genes distribution during the evolution of the
*Dikarya* fungi, using all fungal genomes publicly available at
the start of the study. Fungal (*Eumycota*) DyP coding genes
of the clade D, as defined by Zamocky et al. (2015), are

present in some, but not all, *Dikarya* species. Search for the
most similar bacterial sequences and their placement in the
fungal DyP phylogenetic tree suggest that fungal sequences
may have originated from a single horizontal transfer of a
bacterial class C gene to an ancestor of the *Dikarya* more
than 400 My ago (Taylor et al. 2004), followed by multiple
independent events of gene loss, duplication and diversifica-
tion during the course of evolution of this taxonomic group.

The access to a large (more than 450) number of fungal DyP sequences allowed the identification of distinct clades, defined on the basis of sequence similarities that diverged from each other at different periods during the evolution of the *Dikarya*. In the *Agaricomycotina*, and to a far lesser extent in the *Pezzizomycotina*, DyP diversification did not parallel taxa diversification since DyP genes belonging to different (up to three) clades can co-exist in the genome of a single extant fungal species. The four main DyP clades (III, V, VI-1, and VI-2) present in extant *Agaricomycotina* have thus followed different and apparently independent evolutionary trajectories, being independently eliminated, recruited, or amplified in specific lineages. It seems, for example, that clade VI-2 experienced more frequent episodes of both gene diversification and loss than clade VI-1. Consequently, clade VI-1 is more prevalent in extant agaricomycetous species compared to clade VI-2, but with an average lower number of copies per genome (1.5 versus 3.9). These independent evolutionary trajectories suggest that enzymes from each of the clades may participate in different functions and are under different selective pressures.

DyP clades, which form distinct clusters in network anal-
yses, can also be, to some extent, distinguished from each
other based on “global” protein characteristics (hydropathy
indices, isoelectric points, absolute numbers of aromatic
residues, and abundance of *N*-glycosylation sites or the pres-
ence/absence of an *N*-terminal signal peptide for secretion).
As observed for other enzyme families, including enzymes
implicated in biomass degradation (Aspeborg et al. 2012;
Viborg et al. 2019), these observations suggest that the
divergence between DyPs could be associated with signifi-
cant changes in catalytic properties, substrate range, and cel-
lular functions (Zallot et al. 2021). It is, however, difficult to
speculate on the nature of these changes since current struc-
tural (e.g. 3D structure) and functional information on fungal
DyPs are fragmentary and with regard to related enzymes
belonging exclusively to the phylogenetically-related clades
VI-1 and VI-2. Enzymes belonging to this latter clade share
a number of very distinctive structural features that clearly
distinguish them from those of other clades. The most strik-
ing one is the presence of a sequence peptide that suggests
secretion through the conventional endoplasmic reticulum
pathway. Although SignalP may sometimes fail to predict
signal peptides, we are confident that in our case, this should
be a marginal problem as the presence of a signal peptide
was predicted for almost all of the sequences of clade VI-2. This hypothesis is further supported by the presence in clade VI-2 proteins of a higher number of putative N-glycosylation sites, a common feature of endoplasmic reticulum-secreted proteins (Medus et al. 2017). Although it is tempting to hypothesize that the emergence of clade VI-2 corresponded to a transition from an intracellular to an extracellular status, this hypothesis would need further experimental validation. On the one hand, a majority of the biochemically characterized fungal DyPs have been purified from culture filtrates that precisely belong to clade VI-2 (see Fig. 2). On the other hand, novel DyPs also purified from culture filtrates belong to clade VI-1, from which emerged clade VI-2. Clade VI-1 is the only clade that encompasses both ascomycetous and basidiomycetous sequences, and one of the characterized enzymes was from the wood-rot *Ascomycota X. grammica*. It could thus be hypothesized that clade VI-1 DyPs are either intracellular or extracellular, secreted through non-conventional secretion systems as proposed for a rhamnosi-dase of *Xylaria polymorpha* and other fungal proteins (Nghi et al. 2012; Prudovsky et al. 2003). Acquisition or loss of structural elements favoring secretion of DyPs is also docu-mented in the bacteria, where type A DyPs possess a Tat secretion signal not found in the other bacterial clades such as clade C from which fungal DyPs may derive (de Gonzalo et al. 2016). Secreted type A bacterial DyPs are suspected to be involved in chemical lignin modification (Sugano and Yoshida 2021).

Besides clade VI (VI-1 and VI-2) DyPs, no functional information exists for DyPs belonging to other clades, including their putative cellular localization, although several, but not all, DyPs of rust fungi (*Puccinomycotina*) possess signal peptides. We suggest that future functional studies should focus on these other clades, not only to precise their cellular localization but also on their catalytic activities and substrate range to understand their cellular functions and evaluate their potential use as catalysts in biotechnology.

Regarding the most debated role of fungal DyPs, i.e. their contribution to lignin degradation, our correlative analyses suggest a facultative contribution of a subset of these enzymes in this process, as suggested in Nagy et al. (2017), who considered the fungal DyP family as a whole without splitting it in distinct clades with different evolutionary histories. We show a strong association between clade VI-2 and species of *Agaricales* classified as white-rotters or litter-decomposers (soil saprotrophs) that, for many of them, can degrade lignin and humic substances (Steffen et al. 2000, 2002). Intuitively, enzymes active on bulky extracellular polymers need to be secreted as we hypothesized that they are in the case of clade VI-2 DyPs. This strong association contrasts with the complete absence of this clade from the genomes of either brown rot or ectomycorrhizal species sequenced thus far. In the

case of the *Boletales*, a taxonomic group composed of BR and ECM species, the loss of clade VI-2 DyPs could even represent an ancestral characteristic. BR and ECM species can, however, possess DyPs belonging to other clades. If clade VI-2 DyPs participate in lignin degradation/modification, this role may be rather accessory as (i) many WR species lack these enzymes and (ii) species that possess clade VI-2 DyPs also possess class II peroxidases, including MnPs, VPs, and/or LiPs. Besides *Agaricales*, clade VI-2 DyPs are also present in saprotrophic or plant symbionts/endophytes of the orders *Cantharellales/Sebacinales* that are not considered WR species (Nagy et al. 2016) but possess a rich set of genes encoding plant cell-wall degrading enzymes in their genomes (Kohler et al. 2015). Among these species is *Sistotrema* sp. (*Cantharellales*), sometimes considered a white-rot fungus although not possessing class-II peroxidases (Riley et al. 2014).

Contribution of clade VI-1 enzymes to lignin degradation by the wood-rot ascomycete *X. grammica* has also been suggested as one corresponding enzyme purified from culture filtrates of this species possesses a Mn²⁺-oxidizing activity (Kimani et al. 2021). Incidentally, we observe that several other ascomycetous species possessing clade VI-1 DyPs are associated with decaying wood (e.g. *Ascofurine* sp., *Diaporthe ampelina*, *Eutypa lata*, *Kretzschmaria deusta*, *Rosellinia necatrix*, and *Xylaria* spp.). These genera and species could be targeted to study the transcription profile of their DyPs when grown on lignocellulosic substrates, their (extra)cellular localization, and substrate specificity. Other ascomycetous species possessing clade VI-2 DyPs were however not associated with rotting wood, this is also the case for all ascomycetous species processing clade I or IV DyPs, the two other clades present in this phylum.

Regarding DyPs expressed in the environment, we successfully assembled 14 complete or nearly complete DyP sequences following sequence capture of cDNAs retro-transcribed from soil and decaying wood RNAs. All DyP sequences were affiliated to either clades VI-1 or VI-2 despite the fact that capture probes also targeted sequences from the other fungal DyP clades. Clades VI-1 and VI-2 sequences represent altogether about two-thirds of all DyP present in databases. The exclusive presence of clades VI-1 or VI-2 DyPs among environmental sequences suggests that fungal species harboring DyPs belonging to these clades are either overrepresented and/or that the corresponding genes are overexpressed in the studied environments. These observations are comforted by the data of Kellner et al. (2014), who used PCR to obtain partial DyP sequences from different forest soils. Affiliation of these 61 partial sequences to the different clades showed that 81% originated from clades VI-1 (42%) and VI-2 (24%), while a minority of the sequences came from clades I (4 sequences), III (6), and V(1). The absence of environmental sequences that define
novel DyP clades may signify that fungal genome sequenc-
ing may have already captured the entire phylogenetic diver-
sity of the DyP D subfamily or those additional clades are
too divergent to be captured by the designed probes or too
rare to be easily identified using this approach. The occur-
rence of expressed clades VI-1 and VI-2 DyPs, from both
basidiomycetous and ascomycetous taxa in organic matter
(OM)-rich environments (leaf-litter, uppermost soil) rein-
forces our hypothesis for their involvement in OM degrada-
tion and at the same time should incite us to characterize the
catalytic properties of the more enigmatic and somewhat
structurally divergent members of clades I to V.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00253-022-11923-0.

**Acknowledgements** We would like to acknowledge Dr. Francesco
Venice (University of Turin) for the suggestions in bioinformatics.
The study was performed in the framework of the PeroxiDiv project
co-funded by the ANR (Agence Nationale de la Recherche, France;
ANR-RF-2015-01) and the DFG (Deutsche Forschungsgemeinschaft,
Germany). Martino Adamo was supported by a grant from the Ital-
ian Ministry of Education, University and Research (MIUR) and the
CMIRA (Coopération et mobilité internationales Rhône-Alpes) pro-
gram of region Région Auvergne-Rhône-Alpes. We are grateful to
the Mésocentre Clermont Auvergne University and AuBi platform for
providing computing resources.

**Author contribution** RM, MA, HK, MG, PL, and MH conceived and designed the research. MA, EB, and LV conducted experiments. MA, PP, SCM, and RP analyzed data. MA, MG, and RM wrote the manuscript. All authors read and approved the manuscript.

**Data availability** Sequences generated during the current study are available in the NCBI-nr repository, with accession numbers OM674899-OM674912. All accession numbers of the sequences analyzed during this study are included in this published article (Supplementary Table S4).

This article is licensed under a Creative Commons Attribution 4.0 Inter-
national License, which permits use, sharing, adaptation, distribution
and reproduction in any medium or format, as long as you give appro-
priate credit to the original author(s) and the source, provide a link to
the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in
the article's Creative Commons licence, unless indicated otherwise in
a credit line to the material. If material is not included in the article's
Creative Commons licence and your intended use is not permitted by
statutory regulation or exceeds the permitted use, you will need to
obtain permission directly from the copyright holder. To view a copy of
this licence, visit http://creativecommons.org/licenses/by/4.0/.

References

Acebes S, Ruiz-Dueñas FJ, Toubes M, Sáez-Jiménez V, Pérez-Boada
M, Lucas MF, Martínez AT, Guallar V (2017) Mapping the long-
range electron transfer route in ligninolytic peroxidases. J Phys
Chem B 121:3946–3954. https://doi.org/10.1021/acs.jpcb.7b008
35

Adamo M, Voyron S, Girlanda M, Marmeisse R (2017) RNA extrac-
tion from decaying wood for (Meta)transcriptomic analyses. Can
J Microbiol 63:841–850. https://doi.org/10.1139/cjm-2017-0230

Adamo M, Voyron S, Chialva M, Marmeisse R, Girlanda M (2020) Metabarcoding on both environmental DNA and RNA highlights differences between fungal communities sampled in different habitats. PLoS ONE 15:e0244682. https://doi.org/10.1371/journal.pone.0244682

Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215:403–410. https://doi.org/10.1016/S0022-2836(05)80360-2

Arantes V, Jellison J, Goodell B (2012) Peculiarities of brown-rot fungi and biochemical Fenton reaction with regard to their potential as a model for bioprocessing biomass. Appl Microbiol Biotechnol 94:323–338. https://doi.org/10.1007/s00253-012-3954-y

Aspeborg H, Coutinho PM, Wang Y, Brumer H, Henrissat B (2012) Evolution, substrate specificity and subfamily classification of glycoside hydrolase family 5 (GH5). BMC Evol Biol 12:186. https://doi.org/10.1186/1471-2148-12-186

Ayuso-Fernández I, Ruiz-Dueñas FJ, Martínez AT (2018) Evolutionary convergence in lignin-degrading enzymes. Proc Natl Acad Sci 115:6428–6433. https://doi.org/10.1073/pnas.1802555115

Boerjan W, Ralph J, Baucher M (2003) Lignin biosynthesis. Annu Rev Plant Biol 54:519–546. https://doi.org/10.1146/annurev.arplant.54.031902.134938

Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. https://doi.org/10.1093/bioinformatics/btu170

Bragalini C, Ribière C, Parisot N, Vallon L, Prudent E, Peyretaillade E, Girlanda M, Peyret P, Marmeisse R, Luis P (2014) Solution hybrid selection capture for the recovery of functional full-length eukaryotic cDNAs from complex environmental samples. DNA Res 21:685–694. https://doi.org/10.1093/dnares/dsu030

Buchfink B, Xie C, Huson DH (2015) Fast and sensitive protein alignment using DIAMOND. Nat Methods 12:59–60. https://doi.org/10.1038/nmeth.3176

Catucci G, Valetti F, Sadeghi SJ, Gilardi G (2020) Biochemical features of dye-decolorizing peroxidases: current impact on lignin degradation. Biotechnol Appl Biochem 67:751–759. https://doi.org/10.1002/bab.2015

Celis AI, Dubois JL (2015) Substrate, product, and cofactor: the extraordinarily flexible relationship between the CDE superfamily and heme. Arch Biochem Biophys 574:3–17. https://doi.org/10.1016/j.abb.2015.03.004

Chen K, Durand D, Farach-Colton M (2000) NOTUNG: a program for dating gene duplications and optimizing gene family trees. J Comput Biol 7:429–447. https://doi.org/10.1089/106652700750050871

Chio C, Sain M, Qin W (2019) Lignin utilization: a review of lignin depolymerization from various aspects. Renew Sustain Energy Rev 107:232–249. https://doi.org/10.1016/j.rser.2019.03.008

Cragg SM, Beckham GT, Bruce NC, Bugg TDH, Distel DL, Dupree P, Etxabe AG, Goodell BS, Jellison J, McGeehan JE, McQueen-Mason SJ, Schnorr K, Walton PH, Watts JEM, Zimmer M (2015) Lignocellulose degradation mechanisms across the tree of life. Curr Opin Chem Biol 29:108–119. https://doi.org/10.1016/j.cbpa.2015.10.018

de Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N (2006) ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res 34:W362–W365. https://doi.org/10.1093/nar/gkl124

de Eugenio LI, Peces-Pérez R, Linde D, Prieto A, Barriuso J, Ruiz-Dueñas FJ, Martínez MJ (2021) Characterization of a dye-decolorizing peroxidase from *Irpex lacteus* expressed in *Escherichia coli*: an enzyme with wide substrate specificity able to transform

<figure><img src="image_8.png" /></figure>
lignosulfonates. *J Fungi* 7:325. https://doi.org/10.3390/jof7050325

de Gonzalo G, Colpa DI, Habib MHM, Fraaije MW (2016) Bacterial enzymes involved in lignin degradation. *J Biotechnol* 236:110–119

Denonfoux J, Parisot N, Dugat-Bony E, Biderre-Petit C, Boucher D, Morgavi DP, Le PD, Peyretaillade E, Peyret P (2013) Gene capture coupled to high-throughput sequencing as a strategy for targeted metagenome exploration. *DNA Res* 20:185–196. https://doi.org/10.1093/dnares/dst001

Dinno A (2015) Nonparametric pairwise multiple comparisons in independent groups using Dunn's test. *Stata J* 15:292–300

Duval A, Lawoko M (2014) A review on lignin-based polymeric, micro- and nano-structured materials. *React Funct Polym* 85:78–96. https://doi.org/10.1016/j.reactfunctpolym.2014.09.017

Eastwood DC, Floudas D, Binder M, Majcherczyk A, Schneider P, Aerts A, Asiegbu FO, Baker SE, Barry K, Bendiksby M, Blumen-tritt M, Coutinho PM, Cullen D, de Vries RP, Gathman A, Good-ell B, Henrissat B, Ihrmark K, Kauserud H, Kohler A, LaButti K, Lapidus A, Lavin JL, Lee Y-H, Lindquist E, Lilly W, Lucas S, Morin E, Murat C, Oguiza JA, Park J, Pisabarro AG, Riley R, Rosling A, Salamov A, Schmidt O, Schmutz J, Skrede I, Stenlid J, Wiebenga A, Xie X, Kües U, Hibbett DS, Hoffmeister D, Högberg N, Martin F, Grigoriev IV, Watkinson SC (2011) The plant cell wall-decomposing machinery underlies the functional diversity of forest fungi. *Science* 333:762–765. https://doi.org/10.1126/science.1205411

Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res* 32:1792–1797. https://doi.org/10.1093/nar/gkh340

Fernández-Fueyo E, Linde D, Almendral D, López-Lucendo MF, Ruiz-Duenas FJ, Martínez AT (2015) Description of the first fungal dye-decolorizing peroxidase oxidizing manganese(II). *Appl Microbiol Biotechnol* 99:8927–8942. https://doi.org/10.1007/s00253-015-6665-3

Fernández-Fueyo E, Davó-Siguero I, Almendral D, Linde D, Baratto MC, Pogni R, Romero A, Guallar V, Martínez AT (2018) Description of a non-canonical Mn(II)-oxidation site in peroxidases. *ACS Catal* 8:8386–8395. https://doi.org/10.1021/acscatal.8b02306

Floudas D, Binder M, Riley R, Barry K, Blanchette RA, Henrissat B, Martinez AT, Otillar R, Spatafora JW, Yadav JS, Aerts A, Benoit I, Boyd A, Carlson A, Copeland A, Coutinho PM, de Vries RP, Ferreira P, Findley K, Foster B, Gaskell J, Glotzer D, Gorecki P, Heitman J, Hesse C, Hori C, Igarashi K, Jürgens JA, Kallen N, Kersten P, Kohler A, Kües U, Kumar TKA, Kuo A, LaButti K, Larrondo LF, Lindquist E, Ling A, Lombard V, Lucas S, Lundell T, Martin R, McLaughlin DJ, Morgenstern I, Morin E, Murat C, Nagy LG, Nolan M, Ohm RA, Patyshakuliyeva A, Rokas A, Ruiz-Duenas FJ, Sabat G, Salamov A, Samejima M, Schmutz J, Slot JC, St. John F, Stenlid J, Sun H, Sun S, Syed K, Tsang A, Wiebenga A, Young D, Pisabarro A, Eastwood DC, Martin F, Cullen D, Grigoriev IV, Hibbett DS (2012) The paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. *Science* 336:1715–1719. https://doi.org/10.1126/science.1221748

Fu L, Niu B, Zhu Z, Wu S, Li W (2012) CD-HIT: accelerated for clustering the next-generation sequencing data. *Bioinformatics* 28:3150–3152. https://doi.org/10.1093/bioinformatics/bts565

Gasteiger E (2003) ExPASy: the proteomics server for in-depth protein knowledge and analysis. *Nucleic Acids Res* 31:3784–3788. https://doi.org/10.1093/nar/gkg563

Gerlt JA, Bouvier JT, Davidson DB, Imker HJ, Sadkhin B, Slater DR, Whalen KL (2015) Enzyme function initiative-enzyme similarity tool (EFI-EST): a web tool for generating protein sequence similarity networks. *Biochim Biophys Acta - Proteins Proteomics* 1854:1019–1037. https://doi.org/10.1016/j.bkapap.2015.04.015

Grigoriev IV, Nikitin R, Haridas S, Kuo A, Ohm R, Otillar R, Riley R, Salamov A, Zhao X, Korzeniewski F, Smirnova T, Nordberg H, Dubchak I, Shabalov I (2014) MycoCosm portal: gearing up for 1000 fungal genomes. *Nucleic Acids Res* 42:D699–D704. https://doi.org/10.1093/nar/gkt1183

Gupta R, Jensen LJ, Brunak S (2002) Orphan Protein Function and Its Relation to Glycosylation. In: Bioinformatics and Genome Analysis. Springer Berlin Heidelberg, Berlin, Heidelberg, pp 275–294

Hammel KE, Cullen D (2008) Role of fungal peroxidases in biological ligninolysis. *Curr Opin Plant Biol* 11:349–355. https://doi.org/10.1016/j.pbi.2008.02.003

Hibbett DS, Binder M, Bischoff JF, Blackwell M, Cannon PF, Eriksson OE, Huhndorf S, James T, Kirk PM, Lücking R, Thorsten Lumbsch H, Lutzoni F, Matheny PB, McLaughlin DJ, Powell MJ, Redhead S, Schoch CL, Spatafora JW, Stalpers JA, Vilgalys R, Aime MC, Aptroot A, Bauer R, Begerow D, Benny GL, Castlebury LA, Crous PW, Dai YC, Gams W, Geiser DM, Griffith GW, Gueidan C, Hawksworth DL, Hestmark G, Hosaka K, Humber RA, Hyde KD, Ironside JE, Köljalg U, Kurtzman CP, Larsson KH, Lichtwardt R, Longcore J, Miadlikowska J, Miller A, Moncalvo JM, Mozley-Standridge S, Oberwinkler F, Parmasto E, Reeb V, Rogers JD, Roux C, Ryvarden L, Sampaio JP, Schüßler A, Sugiyama J, Thorn RG, Tibell L, Untereiner WA, Walker C, Wang Z, Weir A, Weiss M, White MM, Winka K, Yao YJ, Zhang N (2007) A higher-level phylogenetic classification of the fungi. *Mycol Res* 111:509–547. https://doi.org/10.1016/j.mycres.2007.03.004

Hofrichter M, Kellner H, Herzog R, Karich A, Liers C, Scheibner K, Kimani VW, Ullrich R (2020) Fungal peroxygenases: a phylogenetically old superfamily of heme enzymes with promiscuity for oxygen transfer reactions. In: Nevalainen H. (eds) Grand Challenges in Fungal Biotechnology. Grand Challenges in Biology and Biotechnology. Springer, Cham, 369–403. Springer, Cham

Huang X (1999) CAP3: a DNA Sequence assembly program. *Genome Res* 9:868–877. https://doi.org/10.1101/gr.9.9.868

Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS (2017) ModelFinder: fast model selection for accurate phylogenetic estimates. *Nat Methods* 14:587–589. https://doi.org/10.1038/nmeth.4285

Kellner H, Luis P, Pecyna MJ, Barbi F, Kapturska D, Krußer D, Zak DR, Marmeisse R, Vandenbol M, Hofrichter M (2014) Widespread occurrence of expressed fungal secretory peroxidases in forest soils. *PLoS One* 9:e95557. https://doi.org/10.1371/journal.pone.0095557

Kim SJ, Ishikawa K, Hirai M, Shoda M (1995) Characteristics of a newly isolated fungus, *Geotrichum candidum* Dec 1, which decolorizes various dyes. *J Ferment Bioeng* 79:601–607. https://doi.org/10.1016/0922-338X(95)94755-6

Kimani V, Ullrich R, Büttner E, Herzog R, Kellner H, Jehmlich N, Hofrichter M, Liers C (2021) First dye-decolorizing peroxidase from an ascomycetous fungus secreted by *Xylaria grammica*. *Bio-molecules* 11:1–17. https://doi.org/10.3390/biom11091391

Kohler A, Kuo A, Nagy LG, Morin E, Barry KW, Buscot F, Canbäck B, Choi C, Cichocki N, Clum A, Colpaert J, Copeland A, Costa MD, Doré J, Floudas D, Gay G, Girlanda M, Henrissat B, Herrmann S, Hess J, Högberg N, Johansson T, Khouja H-R, LaButti K, Lahrmann U, Levasseur A, Lindquist EA, Lipzen A, Marmeisse R, Martino E, Murat C, Ngan CY, Nehls U, Plett JM, Pringle A, Ohm RA, Perotto S, Peter M, Riley R, Rineau F, Ruytinx J, Salamov A, Shah F, Sun H, Tarkka M, Tritt A, Veneault-Fourrey C, Zuccaro A, Tunlid A, Grigoriev IV, Hibbett DS, Martin F (2015) Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. *Nat Genet* 47:410–415. https://doi.org/10.1038/ng.3223

Krahe N-K, Berger RG, Ersoy F (2020) A dyP-type peroxidase of *Pleurotus sapidus* with alkene cleaving activity. *Molecules* 25:1536. https://doi.org/10.3390/molecules25071536
Letunic I, Bork P (2016) Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. *Nucleic Acids Res* 44:W242–W245. https://doi.org/10.1093/nar/gkw290

Li H, Yelle DJ, Li C, Yang M, Ke J, Zhang R, Liu Y, Zhu N, Liang S, Mo X, Ralph J, Currie CR, Mo J (2017) Lignocellulose pretreatment in a fungus-cultivating termite. Proc Natl Acad Sci U S A 114:4709–4714. https://doi.org/10.1073/pnas.1618360114

Liers C, Bobeth C, Pecyna M, Ullrich R, Hofrichter M (2010) DyP-like peroxidases of the jelly fungus *Auricularia auricula-judae* oxidize nonphenolic lignin model compounds and high-redox potential dyes. Appl Microbiol Biotechnol 85:1869–1879. https://doi.org/10.1007/s00253-009-2173-7

Linde D, Pogni R, Cañellas M, Lucas F, Guallar V, Baratto MC, Sinicropi A, Sáez-Jiménez V, Coscolín C, Romero A, Medrano FJ, Ruiz-Dueñas FJ, Martínez AT (2015a) Catalytic surface radical in dye-decolorizing peroxidase: a computational, spectroscopic and site-directed mutagenesis study. Biochem J 466:253–262. https://doi.org/10.1042/BJ20141211

Linde D, Ruiz-Dueñas FJ, Fernández-Fueyo E, Guallar V, Hammel KE, Pogni R, Martínez AT (2015) Basidiomycete DyPs: genomic diversity, structural-functional aspects, reaction mechanism and environmental significance. Arch Biochem Biophys 574:66–74

Linde D, Ayuso-Fernández I, Laloux M, Aguiar-Cervera JE, de Lacey AL, Ruiz-Dueñas FJ, Martínez AT (2021) Comparing ligninolytic capabilities of bacterial and fungal dye-decolorizing peroxidases and class-II peroxidase-catalases. Int J Mol Sci 22:2629. https://doi.org/10.3390/ijms22052629

Martin M (2011) Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 17:10. https://doi.org/10.14806/ej.17.1.200

Martin F, Kohler A, Murat C, Veneault-Fourrey C, Hibbett DS (2016) Unearthing the roots of ectomycorrhizal symbioses. Nat Rev Microbiol 14:760–773

Medus ML, Gomez GE, Zacchi LF, Couto PM, Labriola CA, Labanda MS, Bielsa RC, Clérico EM, Schulz BL, Caramelo JJ (2017) N-glycosylation triggers a dual selection pressure in eukaryotic secretory proteins. Sci Rep 7:8788. https://doi.org/10.1038/s41598-017-09173-6

Miyauchi S, Kiss E, Kuo A, Drula E, Kohler A, Sánchez-García M, Morin E, Andreopoulos B, Barry KW, Bonito G, Buée M, Carver A, Chen C, Cichocki N, Clum A, Culley D, Crous PW, Fauchery L, Girlanda M, Hayes RD, Kéri Z, LaButti K, Lipzen A, Lombard V, Magnuson J, Maillard F, Murat C, Nolan M, Ohm RA, Pangilinan J, Pereira MF, Perotto S, Peter M, Pfister S, Riley R, Sitrit Y, Stielow JB, Szöllősi G, Žifčáková L, Štursová M, Spatafora JW, Tedersoo L, Vaario L-M, Yamada A, Yan M, Wang P, Xu J, Bruns T, Baldrian P, Vilgalys R, Dunand C, Henrissat B, Grigoriev IV, Hibbett D, Nagy LG, Martin FM (2020) Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits. Nat Commun 11:5125. https://doi.org/10.1038/s41467-020-18795-w

Nagy LG, Riley R, Tritt A, Adam C, Daum C, Floudas D, Sun H, Yadav JS, Pangilinan J, Larsson KH, Matsuura K, Barry K, Labutti K, Kuo R, Ohm RA, Bhattacharya SS, Shirouzu T, Yoshinaga Y, Martin FM, Grigoriev IV, Hibbett DS (2016) Comparative genomics of early-diverging mushroom-forming fungi provides insights into the origins of lignocellulose decay capabilities. Mol Biol Evol 33:959–970. https://doi.org/10.1093/molbev/msv337

Nagy LG, Riley R, Bergmann PJ, Krizsan K, Martin FM, Grigoriev IV, Cullen D, Hibbett DS (2017) Genetic bases of fungal white rot wood decay predicted by phylogenomic analysis of correlated gene-phenotype evolution. Mol Biol Evol 34:35–44. https://doi.org/10.1093/molbev/msw238

Nguyen LT, Schmidt HA, Von Haeseler A, Minh BQ (2015) IQ-TREE: a fast and effective stochastic algorithm for estimating

maximum-likelihood phylogenies. Mol Biol Evol 32:268–274. https://doi.org/10.1093/molbev/msu300

Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, Schilling JS, Kennedy PG (2016) FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol 20:241–248. https://doi.org/10.1016/J.FUNECO.2015.06.006

Nielsen H (2017) Predicting secretory proteins with SignalP. Methods Mol Biol 1611:59–73. https://doi.org/10.1007/978-1-4939-7015-5_6

Nghi DH, Bittner B, Kellner H, Jehmlich N, Ullrich R, Pecyna MJ, Nousiainen P, Sipilä J, Huong LM, Hofrichter M, Liers C (2012) Appl Environ Microbiol 78(14):4893–4901. https://doi.org/10.1128/AEM.07588-11

Parisot N, Denonfoux J, Dugat-Bony E, Peyret P, Peyretaillade E (2012) KASpOD-a web service for highly specific and explorative oligonucleotide design. Bioinformatics 28:3161–3162. https://doi.org/10.1093/bioinformatics/bts597

Peng Y, Leung HCM, Yiu SM, Chin FYL (2012) IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 28:1420–1428. https://doi.org/10.1093/bioinformatics/bts174

Prudovsky I, Mandinova A, Soldi R, Bagala C, Graziani I, Landriscina M, Tarantini F, Duarte M, Bellum S, Doherty H, Maciag T (2003) The non-classical export routes: FGF1 and IL-1α point the way. J Cell Sci 116(24):4871–4881. https://doi.org/10.1242/jcs.00872

Qin X, Su X, Luo H, Ma R, Yao B, Ma F (2018) Deciphering lignocellulose deconstruction by the white rot fungus *Irpex lacteus* based on genomic and transcriptomic analyses. Biotechnol Biofuels 11:1–14. https://doi.org/10.1186/s13068-018-1060-9

R Core Team (2019) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria

Riley R, Salamov AA, Brown DW, Nagy LG, Floudas D, Held BW, Levasseur A, Lombard V, Morin E, Otillar R, Lindquist EA, Sun H, LaButti KM, Schmutz J, Jabbour D, Luo H, Baker SE, Pisabarro AG, Walton JD, Blanchette RA, Henrissat B, Martin F, Cullen D, Hibbett DS, Grigoriev IV (2014) Extensive sampling of basidiomycete genomes demonstrates inadequacy of the white-rot/brown-rot paradigm for wood decay fungi. Proc Natl Acad Sci 111:9923–9928. https://doi.org/10.1073/pnas.1400592111

Rineau F, Roth D, Shah F, Smits M, Johansson T, Canbäck B, Olsen PB, Persson P, Grell MN, Lindquist E, Grigoriev IV, Lange L, Tunlid A (2012) The ectomycorrhizal fungus *Paxillus involutus* converts organic matter in plant litter using a trimmed brown-rot mechanism involving Fenton chemistry. Environ Microbiol 14:1477–1487. https://doi.org/10.1111/j.1462-2920.2012.02736.x

Ruiz-Dueñas FJ, Barrasa JM, Sánchez-García M, Camarero S, Miyauchi S, Serrano A, Linde D, Babiker R, Drula E, Ayuso-Fernández I, Pacheco R, Padilla G, Ferreira P, Barriuso J, Kellner H, Castanera R, Alfaro M, Ramírez L, Pisabarro AG, Riley R, Kuo A, Andreopoulos W, LaButti K, Pangilinan J, Tritt A, Lipzen A, He G, Yan M, Ng V, Grigoriev IV, Cullen D, Martin F, Rosso M-N, Henrissat B, Hibbett D, Martínez AT (2021) Genomic analysis enlightens agaricales lifestyle evolution and increasing peroxidase diversity. Mol Biol Evol 38:1428–1446. https://doi.org/10.1093/molbev/msaa301

Salame TM, Knop D, Levinson D, Mabjeesh SJ, Yarden O, Hadar Y (2014) Inactivation of a *Pleurotus ostreatus* versatile peroxidase-encoding gene (*mnp2*) results in reduced lignin degradation. Environ Microbiol 16:265–277. https://doi.org/10.1111/1462-2920.12279

Salvachúa D, Katahira R, Cleveland NS, Khanna P, Resch MG, Black BA, Purvine SO, Zink EM, Prieto A, Martínez MJ, Martínez AT, Simmons BA, Gladden JM, Beckham GT (2016) Lignin
depolymerization by fungal secretomes and a microbial sink. Green Chem 18:6046–6062. https://doi.org/10.1039/C6GC01531J

Savelli B, Li Q, Webber M, Jemmat AM, Robitaille A, Zamocky M, Mathé C, Dunand C (2019) RedoxiBase: a database for ROS homeostasis regulated proteins. Redox Biol 26:101247. https://doi.org/10.1016/j.redox.2019.101247

Scheibner M, Hülsdau B, Zelena K, Nimtz M, De Boer L, Berger RG, Zorn H (2008) Novel peroxidases of *Marasmius scorodonius* degrade β-carotene. Appl Microbiol Biotechnol 77:1241–1250. https://doi.org/10.1007/s00253-007-1261-9

Schmidt O (2006) Wood and tree fungi. Springer-Verlag, Berlin Heidelberg, p 2006

Shah F, Nicolás C, Bentzer J, Ellström M, Smits M, Rineau F, Canbäck B, Floudas D, Carleer R, Lackner G, Braesel J, Hoffmeister D, Henrissat B, Ahrén D, Johansson T, Hibbett DS, Martin F, Persson P, Tunlid A (2016) Ectomycorrhizal fungi decompose soil organic matter using oxidative mechanisms adapted from saprotrophic ancestors. New Phytol 209:1705–1719. https://doi.org/10.1111/nph.13722

Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T (2003) Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13:2498–2504. https://doi.org/10.1101/gr.1239303

Sievers F, Higgins DG (2014) Clustal Omega, accurate alignment of very large numbers of sequences. In: Russell D. (eds) Multiple Sequence Alignment Methods. Methods in Molecular Biology (Methods and Protocols), vol 1079, Humana Press, Totowa, NJ, 105–116

Sigrist CJA, de Castro E, Cerutti L, Cuche BA, Hulo N, Bridge A, Bougueleret L, Xenarios I (2012) New and continuing developments at PROSITE. Nucleic Acids Res 41:D344–D347. https://doi.org/10.1093/nar/gks1067

Singh R, Eltis LD (2015) The multihued palette of dye-decolorizing peroxidases. Arch Biochem Biophys 574:56–65

Steffen KT, Hofrichter M, Hatakka A (2000) Mineralisation of 14 C-labelled synthetic lignin and ligninolytic enzyme activities of litter-decomposing basidiomycetous fungi. Appl Microbiol Biotechnol 54:819–825. https://doi.org/10.1007/s002530000473

Steffen KT, Hatakka A, Hofrichter M (2002) Degradation of humic acids by the litter-decomposing basidiomycete *Collybia dryophila*. Appl Environ Microbiol 68:3442–3448. https://doi.org/10.1128/AEM.68.7.3442-3448.2002

Stothard P (2000) The sequence manipulation suite: JavaScript programs for analyzing and formatting protein and DNA sequences. Biotechniques 28:1102. https://doi.org/10.2144/00286ir01

Strittmatter E, Liers C, Ullrich R, Wachter S, Hofrichter M, Plattner DA, Piontek K (2013) First crystal structure of a fungal high-redox potential dye-decolorizing peroxidase substrate interaction sites and long-range electron transfer. J Biol Chem 288:4095–4102. https://doi.org/10.1074/jbc.M112.400176

Sugano Y (2009) DyP-type peroxidases comprise a novel heme peroxidase family. Cell Mol Life Sci 66:1387–1403

Sugano Y, Yoshida T (2021) DyP-Type Peroxidases: Recent Advances and Perspectives. Int J Mol Sci 22(11):5556. https://doi.org/10.3390/ijms22115556

Sugawara K, Igeta E, Amano Y, Hyuga M, Sugano Y (2019) Degradation of antifungal anthraquinone compounds is a probable physiological role of DyP secreted by *Bjerkandera adusta*. AMB Express 9:56. https://doi.org/10.1186/s13568-019-0779-4

Taylor TN, Hass H, Kerp H, Krings M, Hanlin RT (2004) Perithecial ascomycetes from the 400 million year old rhynie chert: an example of ancestral polymorphism. Mycologia 96:1403. https://doi.org/10.2307/3762157

Trifinopoulos J, Nguyen LT, von Haeseler A, Minh BQ (2016) W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res 44:W232–W235. https://doi.org/10.1093/nar/gkw256

Valette N, Perrot T, Sormani R, Gelhaye E, Morel-Rouhier M (2017) Antifungal activities of wood extractives. Fungal Biol Rev 31(3):113–123. https://doi.org/10.1016/j.fbr.2017.01.002

Viborg AH, Terrapon N, Lombard V, Michel G, Czjzek M, Henrissat B, Brumer H (2019) A subfamily roadmap of the evolutionarily diverse glycoside hydrolase family 16 (GH16). J Biol Chem 294:15973–15986. https://doi.org/10.1074/jbc.RA119.010619

Yoshida T, Tsuge H, Konno H, Hisabori T, Sugano Y (2011) The catalytic mechanism of dye-decolorizing peroxidase DyP may require the swinging movement of an aspartic acid residue. FEBS J 278:2387–2394. https://doi.org/10.1111/j.1742-4658.2011.08161.x

Zallot R, Oberg N, Gerlt JA (2021) Discovery of new enzymatic functions and metabolic pathways using genomic enzymology web tools. Curr Opin Biotechnol 69:77–90. https://doi.org/10.1016/j.copbio.2020.12.004

Zámocký M, Gasselhuber B, Furtmüller PG, Obinger C (2012) Molecular evolution of hydrogen peroxide degrading enzymes. Arch Biochem Biophys 525:131–144

Zámocký M, Hofbauer S, Schaffner I, Gasselhuber B, Nicolussi A, Soudi M, Pirker KF, Furtmüller PG, Obinger C (2015) Independent evolution of four heme peroxidase superfamilies. Arch Biochem Biophys 574:108–119. https://doi.org/10.1016/j.abb.2014.12.025

Zhao RL, Li GJ, Sánchez-Ramírez S, Stata M, Yang ZL, Wu G, Dai YC, He SH, Cui BK, Zhou JL, Wu F, He MQ, Moncalvo JM, Hyde KD (2017) A six-gene phylogenetic overview of basidiomycota and allied phyla with estimated divergence times of higher taxa and a phyloproteomics perspective. Fungal Divers 84:43–74. https://doi.org/10.1007/s13225-017-0381-5

Zitare UA, Habib MH, Rozeboom H, Mascotti ML, Todorovic S, Fraaije MW (2021) Mutational and structural analysis of an ancestral fungal dye-decolorizing peroxidase. FEBS J 288:3602–3618. https://doi.org/10.1111/febs.15687

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.