Ancestral sequence reconstruction as a tool to study the evolution of wood decaying fungi
Iván Ayuso-Fernández1*† 
Gonzalo Molpeceres2† 
Susana Camarero2 
Francisco Javier Ruiz-Dueñas2 
Angel T. Martínez2*- 1Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences (NMBU), Ås, Norway
- 2Centro de Investigaciones Biológicas “Margarita Salas” (CIB), CSIC, Madrid, Spain
The study of evolution is limited by the techniques available to do so. Aside from the use of the fossil record, molecular phylogenetics can provide a detailed characterization of evolutionary histories using genes, genomes and proteins. However, these tools provide scarce biochemical information of the organisms and systems of interest and are therefore very limited when they come to explain protein evolution. In the past decade, this limitation has been overcome by the development of ancestral sequence reconstruction (ASR) methods. ASR allows the subsequent resurrection in the laboratory of inferred proteins from now extinct organisms, becoming an outstanding tool to study enzyme evolution. Here we review the recent advances in ASR methods and their application to study fungal evolution, with special focus on wood-decay fungi as essential organisms in the global carbon cycling.
1 Introduction
1.1 Molecular paleogenetics concept
Ancestral protein resurrection opens fascinating ways to test evolutionary hypotheses otherwise impossible to be addressed. The concept of molecular resurrection stems from the seminal work by Linus Pauling and Emile Zuckerkandl, where they coined the idea of molecular paleogenetics – or how to mathematically infer the ancestral sequences of existing proteins and genes (Pauling et al., 1963). In their work, they proposed a method based on comparing existing sequences to calculate the probability of amino acids in ancestral nodes from protein phylogenies (Figure 1). In this way, the most probable sequence at any ancestral node of a phylogeny could be obtained, and even resurrected in the laboratory to evaluate its biochemical properties. Unfortunately, sequence information in their time was scarce, and it was not until the 90s, with the accumulation of gene (and genome) sequencing information in databases, that the first examples of ancestral proteins resurrected in the laboratory appeared (Malcolm et al., 1990; Stackhouse et al., 1990). From then on, the ideas of Pauling and Zuckerkandl were brought to life, and the field of paleogenetics started a beautiful and fascinating journey. There has been a continuous development and refinement of algorithms (and their implementation as software), together with the identification of the main caveats of ASR and ancestral protein characterization efforts in many protein families. The theoretical and experimental aspects of ASR continue to develop in parallel, and the detailed history behind ASR can be found in dedicated reviews and books (Chang and Donoghue, 2000; Liberles, 2007; Merkl and Sterner, 2016).
Today the examples of resurrected proteins cover dozens of protein families and range all fields of life sciences (Gumulya and Gillam, 2016). Some studies aim to describe the origin of life, pointing to hot ancestral environments as the broth where the first organisms appeared (Gaucher et al., 2003). Physiology and behavior of extinct species can also benefit greatly from ASR studies, as shown with the reconstruction of dinosaur vision (Chang et al., 2002) and the origins of alcohol metabolism in hominids (Carrigan et al., 2015). Analyzing the evolution of enzyme families is appealing, and outstanding works have shown how specific point mutations can restrict substrate specificity (Voordeckers et al., 2012), how plants diversified the production of secondary metabolites (Huang et al., 2012) or even how cancer-related kinases acquired resistance to drugs (Wilson et al., 2015).
But ASR is not only used in evolutionary biochemistry, it also has clear applications in industry and medicine. In this sense, outstanding examples show how ASR has been used to optimize vaccines using computational approaches, such as those targeting HIV by calculation of centralized genes that minimize the genetic distance to circulating strains (Nickle et al., 2003; Arenas and Posada, 2010b). ASR has also been proposed as a central tool for protein engineers (Spence et al., 2021), and could be used together with other tools (e.g. molecular directed evolution) to optimize industrial processes such as those related to the enzymatic degradation and valorization of biomass in the biorefinery context (Alcalde, 2015). Of course, the above is a small selection of a continuously growing number of works where ASR methods are used to study protein evolution or to take advantage of ancestral gene properties.
In this review we will focus on fungi, which stand out as one of the most diverse kingdoms of life. Among them, we will analyze how ASR has been useful to study the wood-decay fungal enzymes as essential components in the global carbon cycle. Wood-decay fungi secrete a highly diverse enzyme consortium to attack the main components of the plant cell wall, being essentially lignin, cellulose and hemicelluloses. Analyzing how they evolved and how their plant cell wall degrading enzymes (PCWDE) acquired the properties they have today is important from a fundamental and an applied point of view.
1.2 ASR methodology
Before diving into fungal evolution and how ASR can help its study, we will briefly describe the methodology to reconstruct ancestral sequences. Since the ultimate goal of ASR is getting proteins resurrected in the laboratory to evaluate an evolutionary hypothesis or characterize interesting properties, an ancestral resurrection experiment is only as good as the inferred sequences. Therefore, special care must be taken to obtain ancestral sequences, and a good understanding of the methods and the sources of uncertainty and ambiguity will improve the overall quality of the reconstruction. Figure 1 shows the simplified pipeline required for enzyme resurrection. Typically, a number of protein sequences up to a few hundred are aligned, and a phylogeny is built. The phylogeny, the aligned extant sequences, and an evolutionary model describing the probabilities of change from one amino acid to another are then used to infer the ancestral sequences. Once these sequences are manually curated, the genes will be synthesized and the proteins resurrected in the laboratory, measuring the properties of interest to assess their modification through evolution.
Figure 1 Simplified pipeline to resurrect ancestral proteins, including computational and experimental work. The sequences of the family of interest are aligned using accurate algorithms. The phylogeny is then obtained, and together with an evolutionary model and the multiple sequence alignment (MSA), the ancestral sequences at each node can be obtained. Once the ancestral sequence is selected, it can be resurrected (produced) in the laboratory to obtain the experimental data of interest (e.g. UV-vis or CD spectral properties, enzyme activity or thermal stability).
The inference of ancestral sequences can be done with three different approaches: maximum parsimony (MP), maximum likelihood (ML) or Bayesian (often called hierarchical Bayesian, HB) methods, the two latter being the so-called probabilistic methods (Selberg et al., 2021). All three require a multiple sequence alignment (MSA) of the proteins of interest and a phylogenetic tree explaining relationships between the orthologues selected, and ML and HB methods allow the use of different models of evolution (often called substitution models) (Arenas, 2015). Taking this into account, when doing ASR, it is crucial to obtain accurate MSAs and phylogenetic trees independent of the method used to infer the ancestral sequences.
Obtaining an accurate MSA of the proteins of interest is far from trivial. The high number of algorithms available make this a difficult choice. It is known that the alignment method has an impact in the ancestral reconstruction and there is potential to introduce biases in the inferred sequences (Vialle et al., 2018), with some methods [MAFFT (Katoh et al., 2002), PRANK (Löytynoja, 2014)] performing better than others with simulated data sets. Trying to systematically analyze how important this choice is, Aadland and Kolaczkowski (2020) tested the accuracy of some of the most popular tools for sequence alignment. Their study shows that to improve accuracy, the best choice is to integrate the results of different algorithms combining the information from different sequence alignments, diminishing this way an important source for errors in the ancestral reconstruction.
The selection of protein substitution models (models of evolution) is a common practice in phylogenetic reconstruction and ASR, often done by obtaining the best-fitting model scored using likelihood-based methods (Darriba et al., 2011). However, there is debate in the field since it has been said that model selection has little importance in the final reconstructions (Abadi et al., 2019; Tao et al., 2020). To evaluate the impact of model selection in ASR, computational studies using real and simulated data have recently shown that the best-fitting model yields the most accurate sequences, and if the best-fitting model cannot be applied, the most similar models are preferred aiming to obtain the best reconstructions possible (Del Amparo and Arenas, 2022).
The phylogeny for ASR is mainly obtained with ML or Bayesian methods. The true phylogeny of the family of interest is rarely known, but some studies have shown that its uncertainty has little impact in the reconstruction (Hanson-Smith et al., 2010). In this sense, ML methods for ASR ignore the uncertainty in the phylogeny, while Bayesian methods incorporate it by sampling the distribution of ancestral states. However, using simulated data, Hanson-Smith et al. (2010) showed that Bayesian approaches to integrate over different topologies of the phylogenetic tree do not improve the accuracy of the reconstruction over ML methods. This proves that ASR is robust to uncertainty since the conditions that generate such uncertainty are also making the ancestral state identical between different trees.
Once an accurate MSA and a phylogeny are obtained, ASR can be performed with three different approaches as described above. MP methods infer the ancestral sequence that explains the minimum number of changes leading to extant proteins. Since they were the first methods employed to obtain ancestral states (Fitch, 1971), they are the less sophisticated and barely used today. ML and HB methods have been continuously improved since the early applications of ASR, and a detailed description of them can be found in other reviews (Merkl and Sterner, 2016; Selberg et al., 2021). Among the advantages of Bayesian methods, it has been argued that they integrate over the uncertainty of the reconstruction and despite this having little impact compared to ML methods, it has also been shown that ML can bias the reconstruction to overrepresentation of common amino acids at each site (Williams et al., 2006; Hobbs et al., 2012). This can lead to less accurate but more stable and active enzymes, reason why ML methods can be preferred if ASR is used with protein engineering purposes (Spence et al., 2021).
The accessibility of ASR methods has increased in the last years with the implementation of algorithms as software and the implementation of automated pipelines on online servers, as summarized in Table 1. In fact, many methods used to perform phylogenetic analyses allow ancestral sequence reconstruction. Among these methods, PAML (latest version 4.7) stands out as one of the most used within ASR (Yang, 2007). The implementation of ASR methods on servers can be automated even to include as input only the unaligned sequences of interest, but the user is encouraged to revise reliability of the results and accuracy of the reconstructed sequences.
Table 1 Relevant methods to perform ASR.
Most of the methods rely on empirical amino acid substitution models and assume that an entire sequence evolves at the same rate. ProtASR2 is an exception, as it allows inclusion of structural information to perform ASR (Arenas et al., 2017; Arenas and Bastolla, 2020). This method is based on substitution models that consider stability of the native, unfolded and misfolded states of a protein to avoid high and low hydrophobicity predictions, reconstructing proteins less biased towards higher stabilities (such as the ML methods based on empirical substitution models) and closer to the folding stability of simulated proteins (Arenas et al., 2015). Another important caveat when performing ASR is that the methods employed assume one unique phylogeny and that recombination events did not happen through evolution, despite recombination being a common and widespread genetic event. Pointing this out and addressing the issue, Arenas and Posada (2010a) analyzed the effect of recombination events in ASR studies. They showed that independently of the ASR method used, not considering recombination can bias ASR performed with nucleotides, codons and proteins, even at low simulated recombination events.
Finally, there is a handful of databases useful to collect the sequence information needed for ASR. Public databases constantly updated include those from NCBI (https://www.ncbi.nlm.nih.gov/) and UniProt (https://www.uniprot.org/), and in the case of biomass degradation two more specialized are JGI-Mycocosm (https://mycocosm.jgi.doe.gov/mycocosm/home) and CAZy (http://www.cazy.org/). There is also one recent database called Revenant, with hand curated information of resurrected proteins, that is linked to many other databases (https://revenant.bioinformatica.org/).
2 Fungal evolution
Given the high diversity of fungi, their evolution is sometimes difficult to study. There is a strong consensus that fungi (and nucleariids) are the sister group of holozoans within the clade Ophistokonta (Baldauf, 2008), and their divergence is usually estimated in the Mesoproterozoic/early Neoproterozoic era (~1000 mya) (Heckman et al., 2001; Parfrey et al., 2011). The last common ancestor of fungi is considered to have been non-filamentous and aquatic, with flagellated spores (James et al., 2006). However, some studies point to an earlier and different origin of fungi that would imply a revision of this phylogeny (Butterfield, 2005; Bengtson et al., 2017). The work by Bengtson and coworkers is especially controversial given the fact that they found fungus-like filamentous fossils originated in the early Paleoproterozoic era (~2400 mya), an origin considerably older (and different) than that previously thought. Additionally, it has been demonstrated that complex multicellularity in fungi appeared in different taxa as a convergent evolutionary adaptation, so the roots of these organisms have to be further investigated in detail (Nagy et al., 2018).
The origin and further evolution of Dikarya (including Ascomycota and Basidiomycota) is less troublesome (Figure 2). The time-calibration of different phylogenetic trees using fungal fossils establishes the origin of Dikarya in the Neoproterozoic era (~750 mya), and the divergence of Ascomycota and Basidiomycota in the Cambric period (~500 mya) (Parfrey et al., 2011; Floudas et al., 2012). The evolution of Basidiomycota shows a rapid and continuous diversification of the different groups that form the phylum (Zhao et al., 2017). Within them, the appearance of the class Agaricomycetes, where wood-rotting fungi are included, is estimated to occur at the end of the Carboniferous period (~300 mya) (Floudas et al., 2012; Zhao et al., 2017) with global geological consequences. The work of Floudas and coworkers postulates that the appearance of wood-rotting fungi, together with their production of the first lignin degrading enzymes, contributed to the end of biomass accumulation in form of coal at the end of the Carboniferous period. However, geoclimatic factors would have also significantly contributed to coal formation under ever-wet tropical conditions, and its decline could also be related to climatic shifts toward drier conditions (Hibbett et al., 2016; Nelsen et al., 2016).
Figure 2 Timeline of the main diversification events in the evolution of wood-decay fungi in the subkingdom Dikaria, phyla Ascomycota and Basidiomycota, subphyla Saccharomycotina, Agaricomycotyna and Pezizomycotina, classes Dacrymycetes, Agaricomycetes and Tremellomycetes and orders Agaricales and Polyporales; Camb., Cambric period; Ordov., Ordovician; Sil., Sillurian; Devon., Devonian; Carbon., Carboniferous; Cretac., Cretaceous; Cen., Cenozoic era; N, Neogene period. Agaric, Agaricales; Polyp, Polyporales. The diversification times were taken from Floudas et al. (2012).
Either way, after Agaricomycetes appeared there was a huge diversification of these fungi, generating the great number of species existing today (~21000). During this speciation, there was a complex evolution in terms of gene duplication/loss events, especially concerning PCWDE genes. In this sense, white-rot fungi (the most efficient lignin degraders) evolved with a higher ratio of gene duplication related to lignin-degrading enzymes, while brown-rot (with a preferential degradation of cellulose and limited degradation of lignin) suffered a clear loss of these genes (Ruiz-Dueñas et al., 2013; Nagy et al., 2017). Although the work of Floudas and coworkers (2012) suggested the common ancestor of Agaricomycetes as a white-rot fungus, a recent work involving the genome sequences of undersampled species shows that white-rot fungi evolved later in the Agaricomycetes (Nagy et al., 2016). Thus, the ancestor of both white- and brown-rot fungi could have been a poor wood degrader ("soft-rot") fungus, with an earlier diversification of wood carbohydrate active enzymes compared to ligninolytic peroxidases (although sampling more fungal genomes is necessary to clarify this point) (Nagy et al., 2016). It is important to note that the classic dichotomy white/brown-rot is useful to compare fungi, but it might be not representative of all degradation strategies. Aside from typical soft-rot decay, there are poor wood-degrading fungi displaying a white-rot like phenotype without ligninolytic genes, which could be a transition between the two classic phenotypes (Riley et al., 2014; Schilling et al., 2020).
2.1 Evolution of agaricomycetes and their PCWDE
Given the importance of lignin removal in the biorefinery context (Ragauskas et al., 2014), many efforts in the past decade were dedicated to sequencing Agaricomycetes genomes and to analyzing their lignocellulolytic enzyme machinery, together with the evaluation of the phylogenetic relationships of species and enzyme families involved in plant biomass degradation.
Especially relevant with respect to the evolution of wood-rotting fungi was the work of Floudas and colleagues (2012), who performed the first comprehensive study of Agaricomycetes genomes by analyzing the enzyme content through evolution, and the special role of ligninolytic peroxidases that will be discussed in detail below. Similar studies focused on Polyporales highlighted the differences between brown-rot and white-rot fungi in evolution, suggesting that igninolytic and related genes were lost in the brown-rot lineage (Ruiz-Dueñas et al., 2013; Ferreira et al., 2015). While the loss of ligninolytic heme peroxidases and the reduction of laccases seem clear in the brown-rot lineage, the enzymes involved in H2O2 production such as glucose–methanol–choline (GMC) oxidoreductases and copper-radical oxidases are widely distributed in the two lineages.
Recent works focused on the order Agaricales demonstrated that the expansion of ligninolytic genes is important not only for wood white-rot fungi, but also for fungi with other lignocellulose-decaying lifestyles (Ruiz-Dueñas et al., 2021). Moreover, one of the peroxidases of Agrocybe pediades, an example of grass litter fungus, displays lignin-degrading capabilities matching those of the white-rot ligninolytic enzymes (Sánchez-Ruiz et al., 2021), and one of the laccases of the same fungus, secreted during first days of solid-state fermentation of wheat straw, shows similar kinetics with lignin-derived phenols as those shown by laccases from white-rot Polyporales (Aza et al., 2021). Given that efficient lignin-decay capabilities appear in different nutritional modes in fungi with a varied PCWDE portfolio, the intriguing evolutionary relationships between distant taxonomical orders are still subject of intense investigation.
2.2 PCWDE families involved in lignocellulose degradation
Understanding the evolution of wood-rotting fungi requires the investigation of the enzyme families involved in lignocellulose degradation. Despite ligninolytic peroxidases having been extensively studied in the past decades due to their role in lignin degradation, wood-decay fungi secrete a consortium of different types of PCWDE. A rough classification can be made according to the main plant-cell wall component they act on. Here we will briefly describe the main families involved in biomass decay based on their role in lignin or cellulose/hemicellulose degradation (Figure 3), and we will outline the evolution of PCWDE in terms of ASR and enzyme resurrection works published to date. Given ASR is a relatively new technique, there is a limited number of publications so far, and in some cases, we extend the scope to related enzyme families from bacteria to show how ASR could help to better understand the evolution PCWDE. The works analyzed here, together with the methods for ASR employed by the respective authors, are summarized in Table 2.
Figure 3 Classic model of enzymatic plant cell wall degradation. Hydrolysis of cellulose involves cellobiohydrolases (CBH) that hydrolyze the cellulose chains from reducing and non-reducing ends, endoglucanases (EG) that cut randomly inside the chains, and β-glucosidases (βGL) that convert the cellobiose units produced by CBHs to glucose. The oxidation of the crystalline cellulose surface by LPMOs results in an internal chain cleavage. The wide set of enzymes involved in the degradation of hemicellulose are represented by a single designation (HC ENZs). Laccases (LAC) directly oxidize the phenolic lignin moiety, or the nonphenolic lignin via redox mediators. Within peroxidases, lignin peroxidases (LiP) are able to oxidize directly the major nonphenolic lignin moiety, manganese peroxidases (MnP) oxidize Mn2+ to Mn3+ (the diffusible Mn3+ chelates oxidize the phenolic moiety of lignin) and versatile peroxidases (VP) catalyze the reactions described for both LiPs and MnPs. Brown-rot fungi also degrade cellulose and hemicellulose via in vivo Fenton chemistry. The glucose-methanol-choline oxidase/dehydrogenase superfamily (GMC) are exemplified by the aryl-alcohol oxidases (AAO) that generate hydrogen peroxide in the redox cycling of aromatic metabolites. Cellulose chains are depicted as linked olive-green hexagons, hemicelluloses in blue and lignin in dark green. Released lignin units and aromatic metabolites are represented as green hexagons.
Table 2 Summary of PCWDE works using ASR, and the methods used by the authors.
2.2.1 Enzymes involved in lignin degradation
Lignin, the most abundant aromatic polymer in nature, is formed by phenyl-propanoid units derived from the oxidative radical coupling of cinnamyl alcohols catalyzed by peroxidases and laccases in the plant cell wall (Vanholme et al., 2010). It gives rigidity to plants while protecting cellulose against pathogen attacks. During plant evolution, its functions have been also linked to UV protection and water-proofing vessels (Weng and Chapple, 2010). The main fungal enzyme families related to lignin depolymerization are heme peroxidases and laccases, the former including class-II ligninolytic peroxidases (PODs) and to some extent the so-called dye-decolorizing peroxidases (DyPs). The GMC superfamily of oxidoreductases should be also considered here given their complementary role, even if they do not modify lignin directly.
Ligninolytic peroxidases are class-II heme peroxidases from the peroxidase-catalase superfamily (Zámocký et al., 2015), classified as AA2 in CAZy database (Lombard et al., 2014) and also known as PODs. They use H2O2 as the electron acceptor to start the catalytic cycle, and they can be classified into three families based on their catalytic properties: i) lignin peroxidases (LiP), enzymes able to oxidize directly the major nonphenolic moiety of lignin (Mester et al., 2001; Sáez-Jiménez et al., 2016); ii) short and long manganese peroxidases (MnP), performing the Mn2+ to Mn3+ oxidation, allowing the formation of diffusible Mn3+ chelates that oxidize the minor phenolic lignin moiety (Fernández-Fueyo et al., 2014); and iii) versatile peroxidases (VP), combining the catalytic properties of LiP and MnP (Ruiz-Dueñas et al., 2009). Fungal ligninolytic peroxidases have a crucial role in lignin degradation and carbon recycling in nature.
The first ancestral character reconstruction studies on fungal PCWDE were done with these enzymes, showing that ancestral MnPs appeared in the Carboniferous period (around 400 mya) providing ancient fungi with new enzymatic tools to degrade recalcitrant polymers, and therefore contributing to the end of biomass accumulation in the form of coal (Floudas et al., 2012). According to the reconstruction of only specific amino acid positions involved in their catalysis, MnPs would have incorporated a solvent exposed catalytic tryptophan midway in the evolution of wood-rotting fungi, generating VPs. Later, ancestral VPs would lose the Mn-binding site, leading to the most efficient LiPs observed today (Floudas et al., 2012; Ruiz-Dueñas et al., 2013). This evolutionary hypothesis involving consecutive changes in the oxidation sites of ligninolytic peroxidases was experimentally proven by Ayuso-Fernández et al. (2017) performing ASR (and ancestor resurrection) in Polyporales, where most wood white-rot fungi are included. The authors sampled 10 Polyporales genomes obtaining a total of 113 curated ligninolytic peroxidase protein sequences to find key enzymes in the evolutionary trajectory towards the most efficient LiPs existing today. Their structural models showed the changes in oxidation sites proposed by Floudas et al. (2012), but ASR allowed the resurrection and characterization of the ancestral peroxidases in the laboratory to confirm the putative activities inferred in silico. The activity assays performed by the authors showed that the common ancestor of Polyporales peroxidases was a MnP, able to oxidize Mn2+ and low redox potential phenols, implying that ancestral Polyporales were in fact able to oxidize lignin, mainly using diffusible Mn3+-chelates. Ancestral MnPs diversified later into the types of ligninolytic peroxidases observed today, and the lineage leading to LiPs included an exposed catalytic tryptophan, generating ancestral VPs. These ancestral versatile enzymes would later lose the Mn-binding site, originating ancestral LiPs and the LiPs existing today, with only the surface exposed tryptophan as oxidation site. The authors showed that there was a progressive increase in the efficiency of Mn2+ oxidation until the incorporation of the surface tryptophan in VPs, indicating an initial preference for the cation diffusion strategy to degrade lignin. However, when the tryptophan appeared in evolution, direct oxidation of the major non-phenolic lignin moiety was possible, and that activity was improved (selected) in the evolution towards the most efficient LiPs, concomitant with the loss of their capability to oxidize Mn2+. This type of evolution incorporating a solvent exposed catalytic tryptophan happened several times and is a convergent trait in the evolution of white-rot fungi, as indicated by the convergence to a catalytic tryptophan independently in two different peroxidase lineages (LiP and VP) in Polyporales (Ayuso-Fernández et al., 2018). Direct oxidation of the major non-phenolic lignin was improved in both lineages, leading to extant LiPs or VPs, but through different changes in evolution and with different catalytic properties. Also, when the catalytic tryptophan appeared, enzyme stability in acidic pH, which is relevant to lignin oxidation in nature, was improved in both lineages. Together with the progressive changes in oxidation sites related to oxidative capabilities, the heme site was also reshaped through evolution. Using protein NMR and redox-potential measurement (Figure 4, right), the authors showed that the redox potential of ligninolytic peroxidases increased with time and was correlated with a subtle rearrangement in the coordination between the so-called proximal histidine and the iron of the heme cofactor (Ayuso-Fernández et al., 2019a).
Figure 4 Main events in the evolution of ligninolytic peroxidases (right), in parallel with lignin evolution in gymnosperms and angiosperms (left). Initially, Mn2+ to Mn3+ oxidation was the main fungal strategy to oxidize lignin from gymnosperms, with only G units. This oxidation would be performed in the three acidic residues forming the Mn-binding site in ancestral MnPs (AncMnP). Once the solvent-exposed tryptophan appeared in the first ancestral VPs (AncVP), the direct oxidation of the major non-phenolic moiety of lignin would be possible and improved through the evolution of peroxidases up to today’s LiPs. Concomitant with the change in oxidation sites, the redox potential and the stability towards acidic pHs increased. In parallel, lignin became more complex and less phenolic with the appearance of angiosperms, including S lignin in its structure, and the above changes in ligninolytic enzymes were a fungal response to this increased complexity. The time-calibration was taken from (Ayuso-Fernández et al., 2019b) and (Morris et al., 2018).
The evolution of Polyporales peroxidases has also been related to plant evolution using lignins from different natural origins (Ayuso-Fernández et al., 2019b). Taking advantage of stopped-flow spectrophotometry and 2D-NMR, it was shown that ancestral ligninolytic peroxidases oxidized lignin from Gymnosperms better than lignin from Angiosperms while, after the incorporation of the solvent-exposed catalytic tryptophan, the preference changed to lignin from Angiosperms. Angiosperm lignin, whose appearance and diversification is more recent than Gymnosperm lignin (Morris et al., 2018), is more complex and less phenolic because it incorporates syringil (S) units, in addition to the guaiacyl (G) units synthesized by Gymnosperms. A time-calibrated peroxidase phylogeny indicates that the solvent-exposed tryptophan, which allows oxidization of the major non-phenolic lignin moiety directly, appeared near the time of Angiosperms diversification (Figure 4, left), pointing to co-evolution between a better defense against some pathogens, and other physiological/anatomical advantages in plants, and a better attack by more efficient lignin degradation in fungi (Morris et al., 2018; Ayuso-Fernández et al., 2019b).
More recently, peroxidase evolution was addressed in a wider array of ligninolytic fungi. 336 class-II peroxidases were identified in 42 genomes of Agaricomycetes, including not only Polyporales but also species of Agaricales, Russulales, Boletales and Amylocorticiales (Ruiz-Dueñas et al., 2021). ASR conducted with these enzymes corroborated the reconstruction studies described above for the Polyporales enzymes. In addition, two independent evolutionary pathways were identified leading to the appearance of the surface catalytic tryptophan in VPs and LiPs from Agaricales, with different transition enzymes compared to those reconstructed in Polyporales. Moreover, additional evolutionary pathways explaining the emergence of different novel MnP subfamilies in Agaricales and Russulales were identified. Although none of these novel MnPs or their reconstructed ancestors have yet been characterized, it has been suggested that their evolution in Agaricales most likely responds to an evolutionary adaptation to the different ligninocellulosic substrates (wood, decayed wood, grass litter or forest litter) on which these fungi grow in nature.
Laccases (EC 1.10.3.2, p-diphenol:dioxygen oxidoreductases, AA1_1 in CAZy), are multicopper oxidases (MCOs) that require O2 as electron acceptor and can directly oxidize a wide arrange of phenolic substrates, as well as aryl amines, N-heterocycles and benzenethiols (Xu, 1996; Baldrian, 2006). During their catalytic cycle, O2 is reduced to H2O as byproduct of the reaction. Given their low catalytic requirements, laccases have a great biotechnological potential. In fact, fungal laccases are the oxidoreductases with the highest number of applications reported to date in a range of industrial sectors, organic chemistry, pulp and paper, food and textile industries or detoxification of pollutants (Rodríguez Couto and Toca Herrera, 2006).
In fungi, MCOs play diverse physiological roles in morphogenesis, spore resistance and pigment formation, stress defense, fungal pathogen/host plant interaction, humus turnover and, as described previously, lignin degradation (Giardina et al., 2010; Janusz et al., 2020). In contrast to ligninolytic peroxidases, that can directly oxidize non-phenolic lignin, laccases have a relatively low redox potential that in principle restricts their oxidation action to the minor phenolic lignin moiety (around 20% of lignin polymer). However, this limitation is overcome in the presence of low molecular-weight compounds such as lignin-derived phenols that act as redox mediators of laccases (Camarero et al., 2005). Once oxidized by laccase, the mediator radical can oxidize recalcitrant substrates such as the nonphenolic subunits of lignin (Cañas and Camarero, 2010; Hilgers et al., 2018). Moreover, in the presence of lignin phenols laccase may contribute to fungal lipid peroxidation in vivo, thus expanding its role in the biodegradation of lignin and other recalcitrant aromatic compounds (Srebotnik and Boisson, 2005; Camarero et al., 2008) acting synergistically with MnPs.
Like ligninolytic peroxidases, laccases are found in high-copy numbers in white-rot fungi. An analysis of laccases sensu stricto in Polyporales suggested that there was only a single gene in the common ancestor of these fungi that appeared in the end of the Jurassic (Savinova et al., 2019). The authors hypothesized that the expansion of laccase genes happened in the second half of the Cretaceous era, relating it to the rise of a lignin more resistant to degradation with the expansion of Angiosperms, and in turn a need to increase the tools that fungi had to degrade plant biomass. By contrast, the results of a recent work covering diverse species of different orders from Agaricomycetes including Polyporales (Ruiz-Dueñas et al., 2021) showed that their ancient common ancestor (also dated in the late Jurassic) already possessed several sensu-stricto laccase genes from which the current laccases diversified.
The combination of ASR and directed evolution methods were shown to be an effective approach for lacasse engineering. The available evolutionary information helps to push the boundaries of enzyme design and promote the development of biocatalysts more suitable for industrial applications. In this way, ASR generates new starting points for enzyme design that are broadly different from the protein sequences of the extant enzymes. The resurrected enzymes can be excellent starting points for directed evolution to rescue promiscuous activities lost during natural evolution that can be used today to confer new enzyme functionalities for biotechnological purposes (Alcalde, 2015). A case in point was carried out in an ASR study using a final set of 87 laccase sequences from Basidiomycota fungi where three ancestral nodes were reconstructed (although only two of them could be biochemically characterized) (Gomez-Fernandez et al., 2020). As a proof of concept, the authors performed a preliminary directed evolution experiment using LacAnc100 ancestor, their most recent reconstructed sequence, towards the oxidation of β-diketones. The final laccase variant produced oxidation rates for 1,3-cyclopentanedione, the model compound used in the high-throughput screening assay, that were 160% higher than those of parental LacAnc100. In addition, the authors analyzed the S224G substitution in the OB-1 laccase obtained by directed evolution of PM1 basidiomycete laccase, which notably improved the activity detected for OB-1 compared to the native enzyme (Maté et al., 2010). This mutation was also an ancestral mutation shared by all resurrected ancestors. Therefore, the knowledge obtained by ASR can also be implemented in the engineering of contemporary enzymes by the incorporation of ancestral mutations.
Dye-decolorizing peroxidases are heme-containing peroxidases present in fungi, bacteria and archaea. They do not belong to the peroxidase-catalase superfamily, having a clear and distinct phylogenetic origin within the chlorite-dismutase superfamily (Zámocký et al., 2015). Fungal DyPs share a structural fold, and since their heme site and catalytic cycle are similar to the above ligninolytic peroxidases, a convergent evolution has been proposed (Linde et al., 2015). Their role in nature remains essentially unknown, and as peroxidases they are overall inefficient, indicating that this might not be their physiological role (Singh and Eltis, 2015), which correlates with their little impact in Agaricales lifestyle evolution (Ruiz-Dueñas et al., 2021). However, related to lignin oxidation, it seems that some fungal DyPs can modify lignin (Linde et al., 2021), but there is an ongoing debate whether other types of DyPs are able to truly oxidize this recalcitrant polymer given their poor performance or the amounts of enzyme needed to observe any effect (Min et al., 2015; Rai et al., 2021). Traditionally, DyPs have been classified in classes A, B, C and D according to primary structural homologies, but a more recent classification considering tertiary structure was proposed, resulting in the alternative designations P, I and V. The details of DyP classification are beyond the scope of this review, but for clarity, the majority of fungal DyPs were classified in family D belonging to the recent type V.
There is only one work using ASR to study fungal DyPs (Zitare et al., 2021), where the authors obtained an ancestral D-type enzyme with the main purpose of establishing a good system to characterize DyPs better. They used 641 sequences from Basidiomycota fungi and reconstructed two ancestral nodes from their phylogeny. Unfortunately, only one of those ancestors (together with an alternative variant to cover ambiguity in the reconstruction) was expressed as a soluble protein. The full characterization of their ancestral DyP showed that its capability to oxidize typical DyP substrates is diminished. Intriguingly, the ancestral DyP is able to oxidize Mn2+, despite apparently having no Mn2+-binding site, which could be linked to an ancestral role in lignin degradation. However, since time calibration was not performed, the correlation of this activity with the ancestral MnPs discussed above is difficult. Interestingly, the authors used their ancestral setup to perform structure-function mutagenesis studies focusing on conserved amino acids in the distal heme site. Their results showed the first direct evidence of the role of two conserved residues in the heme site, related to H2O2 reduction during the catalytic cycle of DyPs. This work shows how ancestral enzymes are useful not only to evaluate evolutionary hypotheses, but also to characterize unknown structural features of enzyme families.
GMC oxidoreductases, forming the glucose-methanol-choline oxidase/dehydrogenase superfamily, assist in the degradation of lignin (and crystalline carbohydrates) by generating the H2O2 used by peroxidases (and lytic polysaccharide monooxygenases, LPMOs), and also by reducing oxidized lignin products to avoid repolymerization (Marzullo et al., 1995). It is a varied group including aryl-alcohol oxidase (AAO), alcohol (methanol) oxidase, cellobiose dehydrogenase, glucose oxidase, glucose dehydrogenase, pyranose oxidase and pyranose dehydrogenase, and all together form the AA3 family in CAZy (Sützl et al., 2019). This large and diverse family of enzymes share a common structural fold that possesses a flavin binding motif for the adenine dinucleotide (FAD) cofactor. There are no ASR works related to fungal GMC oxidoreductases, but ancestors from mammalian flavin-containing monooxygenases (FMOs), proteins that also hold a FAD-binding motif, were resurrected by ASR (Nicoll et al., 2020). These FMO ancestors possessed a well-conserved FAD binding domain and were active, pointing out that ancestral proteins containing this co-factor can be obtained and making future reconstruction of ancestral GMC oxidoreductases plausible.
2.2.2 Enzymes involved in cellulose and hemicellulose degradation
Cellulose is the main component of the plant cell wall and the most abundant polymer on Earth, being an essential renewable source in the biorefinery context (Payne et al., 2015). Despite its simple chemical composition (β-1,4 linked anhydroglucose units), the degradation of cellulose in nature requires the concerted action of multiple enzymes, including both hydrolases and oxidoreductases. The three classical types of glycoside hydrolases (GH) (outlined in detail below) act synergistically on amorphous regions of the cellulose fibers (Figure 3). However, the action of lytic polysaccharide monooxygenases (LPMOs) is required to break cellulose crystallinity, boosting with it the action of GHs. In addition, brown-rot fungi also degrade cellulose and hemicellulose via in vivo Fenton chemistry (Figure 3).
Endoglucanases (EG, endo-1,4-β-D-glucanases, EC 3.2.1.4) randomly cleave β-1,4 bonds in internal amorphous regions of cellulose, generating new non-reducing ends. Their fold can be varied, but all of them display a large cleft with the catalytic amino acids in order to accommodate cellulose fibrils (Davies and Henrissat, 1995). In CAZy, they are classified in structural families GH5, GH6, GH7, GH9, GH12, GH44, GH45, and GH74 (Couturier et al., 2016).
Cellobiohydrolases (CBH, cellulose 1,4-β-cellobiosidases, EC 3.2.1.91) release cellobiose (the glucose disaccharide) from cellulose fragments released by EGs. They can act on reducing or non-reducing ends (Couturier et al., 2016), and are processive enzymes meaning that they can slide through the cellulose fiber to continue cleavage. They are included in families GH6 and GH7 in CAZy.
β-glucosidases (BGL, EC 3.2.1.21) cleave the cellobiose dimer into glucose monomers. They are not processive enzymes, and with the cleavage of cellobiose they cause product inhibition on CBHs, a bottleneck in cellulose degradation (Sørensen et al., 2013). In CAZy they are included in families GH1 to GH3.
There are no examples of fungal ancestral GHs, but a reconstruction of ancestral EGs was carried out using 32 sequences of the EG Cel5A family from bacteria (Barruetabeña et al., 2019). Since the carbohydrate-binding modules (CBM) of the sequences are very heterogeneous and poorly conserved, only the catalytic region was considered for the ASR. The results indicate that the LFCA_EG ancestor, belonging to the oldest node of Firmicutes phyla, is a promising biocatalyst for biotechnological applications since it displays high activity over a broad range of temperature and pH values, and it acts on different substrates and displays high heterologous expression yields. Its thermostability and good activity at high temperatures are important features for commercial applications in the pulp and paper industry or the conversion of cellulosic biomass into fermentable sugars for biofuels production (Bayer et al., 2004; Patel et al., 2019). Interestingly, LFCA_EG also shows processive endoglucanase and exoglucanase activity, although no sequences coding hydrolases with this activity were included in the phylogeny. Unlike the fungal cellulases, the bacterial cellulases can be involved in assembling multi-enzyme complexes, called cellulosomes, together with other enzymes related to the degradation of lignocellulose (Bayer et al., 2004). In this sense, LFCA_EG enhances its activity when it is assembled into a minicellulosome.
As an example of direct biotechnological application, the LFCA_EG ancestor was tested in the transformation of cellulosic resources into nanocellulose (nano-sized form of cellulose), which is a high value material suitable for high-performance applications such as tissue engineering. Surprisingly, this enzyme can single-handedly produce pure nanocellulose particles bellow 500 nm with elevated stability, crystallinity, and controlled aspect ratio. This was reported to be the first attempt to generate nanocellulose by a single enzyme in an efficient manner (Alonso-Lerma et al., 2020).
Lytic polysaccharide monooxygenases (LPMOs) are mono-copper enzymes with a characteristic flat surface capable of breaking the crystallinity of cellulose (and chitin) by hydroxylation of the 1 or 4 position of the β-1,4 glycosidic bond. Their discovery is recent (Vaaje-Kolstad et al., 2010), but since LPMO discovery many types with different properties and roles in nature have been described, forming families AA9-AA11 and AA13-AA17 in CAZy. Fungal LPMOs belong classically to family AA9, but some AA13 and AA16 fungal LPMOs have been described. The high diversity of AA9s (even within the same species) has been related to their distinct action on cellulose and complex and varied hemicelluloses (Monika et al., 2022). Together with the recent discovery showing their use of H2O2 as the true co-substrate in nature (Bissaro et al., 2017), and their strong connections in co-expression matrix analyses (Arntzen et al., 2020), LPMOs became a central piece in enzymatic biomass degradation. A first study on LPMO ASR has been recently reported, showing how the reshape of their surfaces could be important for substrate specificity (Ayuso-Fernández et al., 2022).
Hemicelluloses degradation. Hemicelluloses are the third major component of plant cell walls, after cellulose and lignin. They possess a complex chemical structure, being a mixture of branched polysaccharides with three types of backbones and many types of substitutions (Saha, 2003; van den Brink and de Vries, 2011). Degradation of hemicelluloses requires a myriad of enzymes, both for the cleavage of the backbones and the release of the substitutions in those backbones. This multienzymatic and complex degradation is out of the scope of this work, but has been reviewed in detail elsewhere (van den Brink and de Vries, 2011; Li et al., 2022).
Carbohydrate-binding modules. All the aforementioned enzymes acting on cellulose or hemicellulose frequently contain an additional domain or domains that guide the binding to the polymer, known as carbohydrate-binding modules (CBM) (Boraston et al., 2004). CBMs are usually small in size, flexible and stable, and they are widely distributed in nature (Liu et al., 2022). There is a massive number of CBMs annotated in CAZy, currently classified in 91 families. The evolution of these noncatalytic domains is complex, with many such domains appearing across different enzyme families, which could indicate domain transferences between genes and species.
2.3 Lipases
Lipolytic enzymes (EC 3.1.1) are ubiquitously produced in nature and they encompass a diverse group of hydrolases catalyzing the cleavage and formation of ester bonds (Ali et al., 2012). Their classification is complex, and several systems have been used over time (Sarmah et al., 2018). Although they do not directly modify cellulose/hemicellulose or lignin, their main role in biomass degradation has been related to the synergistic effect of feruloyl esterases with cellulases, xylanases and pectinases to degrade complex carbohydrates (Wong, 2006; Arntzen et al., 2020). Fungal lipases also interact with cuticle, the cutin polyesters and waxes barrier that covers plants aerial surfaces (Arya et al., 2021).
Roughly, lipases catalyze the hydrolysis or synthesis of a broad range of different carboxylic esters, showing high specificity towards glyceridic substrates. Sterol esterases (EC 3.1.1.13) show the same capability, but acting on sterol esters as their natural substrates (Hasan et al., 2006; Vaquero et al., 2016). However, in fungi it has been reported that a group of enzymes included in the Candida rugosa (recently designated Diutina rugosa) like lipase family (abH03.01, homologous family in the Lipase Engineering Database) show a broad substrate specificity combining both lipases and sterol esterases properties acting on acylglycerols and sterol esters. Due to the wide specificity of members of this group, it has been proposed to reclassify them as “Versatile Lipases” (Barriuso et al., 2016). These enzymes have been related to the initial degradation of epicuticular waxes and cuticle (Heredia, 2003). Degradation of these external layers makes the plant polysaccharides more accessible to the different CAZymes.
Structurally, all enzymes mentioned above possess a substrate-binding pocket consisting of a large hydrophobic cavity covered by a mobile amphipathic α-helix (named lid or flap). The lid stays closed in an aqueous solution under physiological conditions. However, when the enzyme is in the presence of hydrophobic substrates, the lid rearranges its position (creating an open gate) and the catalytic site becomes accessible (Rodríguez-Salarichs et al., 2021).
To study the evolution of the versatile lipases in fungi, 127 sequences from Agaricales were selected for ASR (Barriuso and Martínez, 2017). Unfortunately, in this study the predicted ancestral proteins were not expressed, so experimental characterization was not available, and conclusions had to be drawn from in silico predictions related to amino acid sequences and 3D homology models. The inferred sequences corresponding to ancestral nodes selected as key intermediates in the evolution showed that both the number of amino acids and hydrophobicity increase throughout the versatile lipases’ lineage in the lid region. In fact, homology models from intermediate nodes showed that the lid of the oldest ancestors held a loop instead of the α-helix present in the lid of extant lipases. Thus, the gradual extension of this region might have led to the formation of the α-helix. It is important to note that the lid is responsible for substrate recognition, so changes in the hydrophobicity and length of this region may be involved in substrate specificity (Mancheño et al., 2003). Similarly, the size of the putative substrate-binding pocket as well as the length of the intramolecular tunnel also changed through evolution (Barriuso and Martínez, 2017). The substrate-binding site is located in an internal and hydrophobic tunnel, whose shape and amino acid composition could be related to the different substrate specificities and catalytic properties of these enzymes (Mancheño et al., 2003). A hypothesis based on similar enzymes is that a longer intramolecular tunnel might enable the communication of the substrate-binding pocket with the outside of the protein, creating a putative novel exit of the reaction products.
Although not directly related to the scope of this review, we switch kingdoms to focus on lipolytic enzymes from bacteria to discuss a couple examples of ancestral lipases characterized experimentally. Bacterial lipases are divided into an increasing number of families with the identification and characterization of novel enzymes. In particular, family I (true lipases) is the largest family and is divided into several subfamilies (Arpigny and Jaeger, 1999; Jaeger and Eggert, 2002; Kovacic et al., 2018). An ASR study was carried out using 83 sequences of bacterial lipases, and the oldest node was selected for resurrection (Rozi et al., 2022). Family 1.3 was selected for ASR due to their singularity in structure, amino acid sequences and secretion mechanism. The last common ancestor of family 1.3, which the authors named “LUCA”, showed remarkable properties related to its substrate profiling with a wider promiscuity for medium and long chain substrates (C10-C16 substrates) compared to present-day 1.3 lipases. It was hypothesized that evolution could have directed modern 1.3 lipases to prefer medium chain substrates. On the other hand, LUCA displayed an optimum temperature at 70°C, which is one of the highest optimum temperatures observed for other modern 1.3 lipases. The LUCA ancestral enzyme possessed other biochemical characteristics such as thermostability, operability at high pH values, and tolerance towards different organic solvents, making it a potentially valuable biocatalyst for industry. Again, this result suggests that the optimum temperature of modern 1.3 lipases decreased through evolution to adapt to environmental conditions. An extreme example of this evolutionary adaptation would be the cold-adapted lipase AMS8 isolated from an Antarctic Pseudomonas that reached the optimum activity at 30°C (Ganasen et al., 2016).
3 Applicability of ASR in protein engineering
Biomass pretreatment and enzymatic hydrolysis are the most expensive steps in biomass upcycling in biorefineries (Payne et al., 2015). For that reason, significant efforts have been made to not only discover new and more efficient enzymes, but also to improve the existing enzymes through protein engineering (Percival Zhang et al., 2006; Alcalde, 2015). ASR becomes then an excellent tool to explore an increased sequence space with non-existing proteins, mining for enzymes with interesting catalytic properties.
Aside from its use in developing evolutionary hypotheses, ASR potential in protein engineering is well documented (Risso et al., 2018; Spence et al., 2021). One of the typical features usually observed in ancestral enzymes is their higher thermostability, with improved melting temperatures of up to +30°C compared with extant enzymes (Trudeau et al., 2016). Commonly, this feature is noticed in proteins predicted for up to a billion years ago (Risso et al., 2014), but it is not perfectly clear if this is an artifact of reconstructions based in ML methods or a true trait of ancestral proteins (Wheeler et al., 2016). A good example of the ASR methods employed to obtain stable variants is the production of a hyper-thermostable L-amino acid oxidase used to perform deracemization to D-amino acids. With a Tm > 95°C and long-term thermal stability, this enzyme is a perfect candidate for industrial application (Ishida et al., 2021).
ASR is used also to study structure/function relationships. Sampling intermediates in evolution can give essential clues to understand how active sites are shaped, which is crucial to engineer new activities or properties in enzymes (Voordeckers et al., 2012; Ayuso-Fernández et al., 2018; Spence et al., 2021). Also, ASR can aid in better understanding of enzyme function to guide protein engineering efforts. In this case ASR methods have been used to show unprecedent allosteric mechanisms (Gamiz-Arco et al., 2021), and to study the evolution of protein-protein interactions (Jemth et al., 2018). Ancestral enzymes are useful also to expand the reactions observed in nature with artificial enzymes. Catalytic promiscuity through evolution is an ongoing debate (Copley, 2015; Glasner et al., 2020), but if the family of interest displays an array of activities, ASR can be used to obtain versatile enzymes or even catalysts with unique reactivities (Huang et al., 2012; Voordeckers et al., 2012; Babkova et al., 2017).
4 Conclusions
Since Pauling and Zuckerkandl conceived the idea of paleogenetics, the development of ASR methods has led to a continuously increasing number of ancestral proteins resurrected in the laboratory, both for protein engineering purposes and to explain evolutionary histories. This is the case also for the limited number of ASR works studying wood-rotting fungi. But this must be seen as an opportunity: ASR has only scratched the surface of possibilities in these fascinating organisms. Most of PCWDE families are still not studied with ASR, and in those already studied, there is still room to resurrect more ancestors or to analyze different subfamilies and study them in the laboratory. Here we have reviewed the current knowledge of ancestral enzymes of wood-rotting fungi, together with a summary of the methods available to do ASR. We foresee a growing number of works in the upcoming years consolidating ASR as an important tool to study fungal evolution and to improve the applicability of fungal enzymes.
Author contributions
IA-F and GM reviewed the literature and wrote the initial draft. FR-D, SC, and ATM performed the critical revision of the article. All the authors participated in the final editing and approved the submitted version.
Funding
This research was funded by the Genobioref (BIO2017-86559-R) and Lig2Plast (PID2021-126384OB-I00) projects of the Spanish Ministry of Science and Innovation (co-financed by FEDER funds); by the WoodZymes (H2020-BBI-JTI-2017-792070) and CuBE – ERC Synergy Project (H2020-ERC-2019-SyG-856446) EU projects, and by the Consejo Superior de Investigaciones Científicas project PIE-202120E019 and grant 2021AEP106, SusPlast platform, and program for the Spanish Recovery, Transformation and Resilience Plan funded by the Recovery and Resilience Facility of the European Union, established by the Regulation (EU) 2020/2094. G.M. acknowledges The Tatiana Pérez de Guzmán el Bueno Foundation for his predoctoral Environment grant. Finally, we acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
Acknowledgments
We thank the editor and the referees for their time and effort spent to improve our review.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Aadland K., Kolaczkowski B. (2020). Alignment-integrated reconstruction of ancestral sequences improves accuracy. Genome Biol. Evol. 12 (9), 1549–1565. doi: 10.1093/gbe/evaa164
Abadi S., Azouri D., Pupko T., Mayrose I. (2019). Model selection may not be a mandatory step for phylogeny reconstruction. Nat. Commun. 10 (1), 934. doi: 10.1038/s41467-019-08822-w
Alcalde M. (2015). Engineering the ligninolytic enzyme consortium. Trends Biotechnol. 33 (3), 155–162. doi: 10.1016/j.tibtech.2014.12.007
Ali Y., Verger R., Abousalham A. (2012). “Lipases or esterases: Does it really matter?,” in Toward a new bio-Physico-Chemical classification BT - lipases and phospholipases: Methods and protocols, vol. pp . Ed. Sandoval G. (Totowa, USA:Humana Press), 31–51). doi: 10.1007/978-1-61779-600-5_2
Alonso-Lerma B., Barandiaran L., Ugarte L., Larraza I., Reifs A., Olmos-Juste R., et al. (2020). High performance crystalline nanocellulose using an ancestral endoglucanase. Commun. Materials. 1 (1), 57. doi: 10.1038/s43246-020-00055-5
Arenas M. (2015). Trends in substitution models of molecular evolution. Frontiers in genetics. 6, 319. doi: 10.3389/fgene.2015.00319
Arenas M., Bastolla U. (2020). ProtASR2: Ancestral reconstruction of protein sequences accounting for folding stability. Methods Ecol. Evol. 11 (2), 248–257. doi: 10.1111/2041-210X.13341
Arenas M., Posada D. (2010a). The effect of recombination on the reconstruction of ancestral sequences. Genetics 184 (4), 1133–1139. doi: 10.1534/genetics.109.113423
Arenas M., Posada D. (2010b). Computational design of centralized HIV-1 genes. Curr. HIV Res. 8 (8), 613–621. doi: 10.2174/157016210794088263
Arenas M., Sánchez-Cobos A., Bastolla U. (2015). Maximum-likelihood phylogenetic inference with selection on protein folding stability. Mol. Biol. Evol. 32 (8), 2195–2207. doi: 10.1093/molbev/msv085
Arenas M., Weber C. C., Liberles D. A., Bastolla U. (2017). ProtASR: An evolutionary framework for ancestral protein reconstruction with selection on folding stability. Systematic. Biol. 66 (6), 1054–1064. doi: 10.1093/sysbio/syw121
Arntzen M.Ø., Bengtsson O., Várnai A., Delogu F., Mathiesen G., Eijsink V. G. H. (2020). Quantitative comparison of the biomass-degrading enzyme repertoires of five filamentous fungi. Sci. Rep. 10 (1), 20267. doi: 10.1038/s41598-020-75217-z
Arpigny J. L., Jaeger K. E. (1999). Bacterial lipolytic enzymes: Classification and properties. Biochem. Journal. 343. Pt. 1(Pt. 1), 177–183. doi: 10.1042/bj3430177
Arya G. C., Sarkar S., Manasherova E., Aharoni A., Cohen H. (2021). The plant cuticle: An ancient guardian barrier set against long-standing rivals. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.663165
Ashkenazy H., Penn O., Doron-Faigenboim A., Cohen O., Cannarozzi G., Zomer O., et al. (2012). FastML: A web server for probabilistic reconstruction of ancestral sequences. Nucleic Acids Res. 40, W580–W584. doi: 10.1093/nar/gks498
Ayuso-Fernández I., De Lacey A. L., Cañada F. J., Ruiz-Dueñas F. J., Martínez A. T. (2019a). Increase of redox potential during the evolution of enzymes degrading recalcitrant lignin. Chem. – A. Eur. J. 25 (11), 2708–2712. doi: 10.1002/chem.201805679
Ayuso-Fernández I., Golten O., Hall K., Sørlie M., Kjendseth Røhr Å., Eijsink V. G. H. (2022). “Ancestral sequence reconstruction of lytic polysaccharide monooxygenases,” in Università di Siena. Ed. Oxizymes (Siena, Italy:Università di Siena), 73–74.
Ayuso-Fernández I., Martínez A. T., Ruiz-Dueñas F. J. (2017). Experimental recreation of the evolution of lignin-degrading enzymes from the Jurassic to date. Biotechnol. Biofuels 10 (1), 67. doi: 10.1186/s13068-017-0744-x
Ayuso-Fernández I., Rencoret J., Gutiérrez A., Ruiz-Dueñas F. J., Martínez A. T. (2019b). Peroxidase evolution in white-rot fungi follows wood lignin evolution in plants. Proc. Natl. Acad. Sci. 116 (36), 17900 LP–17905. doi: 10.1073/pnas.1905040116
Ayuso-Fernández I., Ruiz-Dueñas F. J., Martínez A. T. (2018). Evolutionary convergence in lignin-degrading enzymes. Proc. Natl. Acad. Sciences. 115(25). 6428. LP. –, 6433. doi: 10.1073/pnas.1802555115
Aza P., Molpeceres G., Ruiz-Dueñas F. J., Camarero S. (2021). “Heterologous expression, engineering and characterization of a novel laccase of agrocybe pediades with promising properties as biocatalyst,” in Journal of fungi (Basel, Switzerland:MDPI) (5). doi: 10.3390/jof7050359
Babkova P., Sebestova E., Brezovsky J., Chaloupkova R., Damborsky J. (2017). Ancestral haloalkane dehalogenases show robustness and unique substrate specificity. Chembiochem : A. Eur. J. Chem. Biol. 18 (14), 1448–1456. doi: 10.1002/cbic.201700197
Baldauf S. L. (2008). An overview of the phytogeny and diversity of eukaryotes. In. overview. phytogeny. Diversity eukaryotes. (Vol. 46. Issue. 3. pp, 263–273). doi: 10.3724/SP.J.1002.2008.08060
Baldrian P. (2006). Fungal laccases – occurrence and properties. FEMS Microbiol. Rev. 30 (2), 215–242. doi: 10.1111/j.1574-4976.2005.00010.x
Barriuso J., Martínez M. J. (2017). Evolutionary history of versatile-lipases from agaricales through reconstruction of ancestral structures. BMC Genomics 18 (1), 12. doi: 10.1186/s12864-016-3419-2
Barriuso J., Vaquero M. E., Prieto A., Martínez M. J. (2016). Structural traits and catalytic versatility of the lipases from the Candida rugosa-like family: A review. Biotechnol. Adv. 34 (5), 874–885. doi: 10.1016/j.biotechadv.2016.05.004
Barruetabeña N., Alonso-Lerma B., Galera-Prat A., Joudeh N., Barandiaran L., Aldazabal L., et al. (2019). Resurrection of efficient precambrian endoglucanases for lignocellulosic biomass hydrolysis. Commun. Chem. 2 (1), 76. doi: 10.1038/s42004-019-0176-6
Bayer E. A., Belaich J.-P., Shoham Y., Lamed R. (2004). The cellulosomes: multienzyme machines for degradation of plant cell wall polysaccharides. Annu. Rev. Microbiol. 58, 521–554. doi: 10.1146/annurev.micro.57.030502.091022
Bengtson S., Rasmussen B., Ivarsson M., Muhling J., Broman C., Marone F., et al. (2017). Fungus-like mycelial fossils in 2.4-billion-year-old vesicular basalt. Nat. Ecol. Evol. 1 (6), 141. doi: 10.1038/s41559-017-0141
Bissaro B., Røhr Å.K., Müller G., Chylenski P., Skaugen M., Forsberg Z., et al. (2017). Oxidative cleavage of polysaccharides by monocopper enzymes depends on H2O2. Nat. Chem. Biol. 13 (10), 1123–1128. doi: 10.1038/nchembio.2470
Boraston A. B., Bolam D. N., Gilbert H. J., Davies G. J. (2004). Carbohydrate-binding modules: Fine-tuning polysaccharide recognition. Biochem. Journal. 382(Pt. 3), 769–781. doi: 10.1042/BJ20040892
Bouckaert R., Vaughan T. G., Barido-Sottani J., Duchêne S., Fourment M., Gavryushkina A., et al. (2019). BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis. PloS Comput. Biol. 15, 1–28. doi: 10.1371/journal.pcbi.1006650
Butterfield N. J. (2005). Probable proterozoic fungi. Paleobiology 31 (1), 165–182. doi: 10.1666/0094-8373(2005)031<0165:PPF>2.0.CO;2
Cai W., Pei J., Grishin N. V. (2004). Reconstruction of ancestral protein sequences and its applications. BMC Evolutionary. Biol. 4, 33. doi: 10.1186/1471-2148-4-33
Camarero S., CaÑas A. I., Nousiainen P., Record E., Lomascolo A., MartÍnez M. J., et al. (2008). p-hydroxycinnamic acids as natural mediators for laccase oxidation of recalcitrant compounds. Environ. Sci. Technol. 42 (17), 6703–6709. doi: 10.1021/es8008979
Camarero S., Ibarra D., Martínez M. J., Martínez A. T. (2005). Lignin-derived compounds as efficient laccase mediators for decolorization of different types of recalcitrant dyes. Appl. Environ. Microbiol. 71 (4), 1775–1784. doi: 10.1128/AEM.71.4.1775-1784.2005
Cañas A. I., Camarero S. (2010). Laccases and their natural mediators: Biotechnological tools for sustainable eco-friendly processes. Biotechnol. Adv. 28 (6), 694–705. doi: 10.1016/j.biotechadv.2010.05.002
Carrigan M. A., Uryasev O., Frye C. B., Eckman B. L., Myers C. R., Hurley T. D., et al. (2015). Hominids adapted to metabolize ethanol long before human-directed fermentation. Proc. Natl. Acad. Sci. 112 (2), 458 LP–463. doi: 10.1073/pnas.1404167111
Chang B. S., Donoghue M. J. (2000). Recreating ancestral proteins. Trends Ecol. Evol. 15 (3), 109–114. doi: 10.1016/s0169-5347(99)01778-4
Chang B. S. W., Jönsson K., Kazmi M. A., Donoghue M. J., Sakmar T. P. (2002). Recreating a functional ancestral archosaur visual pigment. Mol. Biol. Evol. 19 (9), 1483–1489. doi: 10.1093/oxfordjournals.molbev.a004211
Copley S. D. (2015). An evolutionary biochemist’s perspective on promiscuity. Trends Biochem. Sci. 40 (2), 72–78. doi: 10.1016/j.tibs.2014.12.004
Couturier M., Bennati-Granier C., Urio M. B., Ramos L. P., Berrin J.-G. (2016). Fungal Enzymatic Degradation of Cellulose BT - Green Fuels Technology. Biofuels, pp 133–146). doi: 10.1007/978-3-319-30205-8_6
Darriba D., Taboada G. L., Doallo R., Posada D. (2011). ProtTest 3: Fast selection of best-fit models of protein evolution. Bioinf. (Oxford. England). 27 (8), 1164–1165. doi: 10.1093/bioinformatics/btr088
Davies G., Henrissat B. (1995). Structures and mechanisms of glycosyl hydrolases. Structure. (London. England.: 1993). 3 (9), 853–859. doi: 10.1016/S0969-2126(01)00220-9
Del Amparo R., Arenas M. (2022). Consequences of substitution model selection on protein ancestral sequence reconstruction. Mol. Biol. Evol. 39 (7), msac144. doi: 10.1093/molbev/msac144
Fernández-Fueyo E., Acebes S., Ruiz-Dueñas F. J., Martínez M. J., Romero A., Medrano F. J., et al. (2014). Structural implications of the c-terminal tail in the catalytic and stability properties of manganese peroxidases from ligninolytic fungi. Acta Crystallographica. Section. D Biol. Crystallography. 70(Pt. 12), 3253–3265. doi: 10.1107/S1399004714022755
Ferreira P., Carro J., Serrano A., Martínez A. T. (2015). A survey of genes encoding H2O2-producing GMC oxidoreductases in 10 polyporales genomes. Mycologia 107 (6), 1105–1119. doi: 10.3852/15-027
Fitch W. M. (1971). Toward defining the course of evolution: Minimum change for a specific tree topology. Systematic. Biol. 20 (4), 406–416. doi: 10.1093/sysbio/20.4.406
Floudas D., Binder M., Riley R., Barry K., Blanchette R. A., Henrissat B., et al. (2012). The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. Sci. (New. York. N.Y.) 336, 1715–19. doi: 10.1126/science.1221748
Gamiz-Arco G., Gutierrez-Rus L. I., Risso V. A., Ibarra-Molero B., Hoshino Y., Petrović D., et al. (2021). Heme-binding enables allosteric modulation in an ancient TIM-barrel glycosidase. Nat. Commun. 12 (1), 380. doi: 10.1038/s41467-020-20630-1
Ganasen M., Yaacob N., Rahman R. N. Z. R. A., Leow A. T. C., Basri M., Salleh A. B., et al. (2016). Cold-adapted organic solvent tolerant alkalophilic family I.3 lipase from an Antarctic. Pseudomonas. Int. J. Biol. Macromol. (92), 1266–1276. doi: 10.1016/j.ijbiomac.2016.06.095
Gaucher E. A., Thomson J. M., Burgan M. F., Benner S. A. (2003). Inferring the palaeoenvironment of ancient bacteria on the basis of resurrected proteins. Nature 425 (6955), 285–288. doi: 10.1038/nature01977
Giardina P., Faraco V., Pezzella C., Piscitelli A., Vanhulle S., Sannia G. (2010). Laccases: A never-ending story. Cell. Mol. Life Sci. 67 (3), 369–385. doi: 10.1007/s00018-009-0169-1
Glasner M. E., Truong D. P., Morse B. C. (2020). How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation. FEBS J. 287 (7), 1323–1342. doi: 10.1111/febs.15185
Gomez-Fernandez B. J., Risso V. A., Rueda A., Sanchez-Ruiz J. M., Alcalde M. (2020). Ancestral resurrection and directed evolution of fungal mesozoic laccases. Appl. Environ. Microbiol. 86 (14), 1–15. doi: 10.1128/AEM.00778-20
Gumulya Y., Gillam E. M. J. (2016). Exploring the past and the future of protein evolution with ancestral sequence reconstruction: the ‘retro’ approach to protein engineering. Biochem. J. 474 (1), 1–19. doi: 10.1042/BCJ20160507
Hanson-Smith V., Kolaczkowski B., Thornton J. W. (2010). Robustness of ancestral sequence reconstruction to phylogenetic uncertainty. Mol. Biol. Evol. 27 (9), 1988–1999. doi: 10.1093/molbev/msq081
Hasan F., Shah A. A., Hameed A. (2006). Industrial applications of microbial lipases. Enzyme Microbial. Technol. 39 (2), 235–251. doi: 10.1016/j.enzmictec.2005.10.016
Heckman D. S., Geiser D. M., Eidell B. R., Stauffer R. L., Kardos N. L., Hedges S. B. (2001). Molecular evidence for the early colonization of land by fungi and plants. Sci. (New. York. N.Y.). 293 (5532), 1129–1133. doi: 10.1126/science.1061457
Heredia A. (2003). Biophysical and biochemical characteristics of cutin, a plant barrier biopolymer. Biochim. Biophys. Acta (BBA). - Gen. Subj. 1620 (1), 1–7. doi: 10.1016/S0304-4165(02)00510-X
Hibbett D., Blanchette R., Kenrick P., Mills B. (2016). Climate, decay, and the death of the coal forests. Curr. Biol 26, R536–R567. doi: 10.1016/j.cub.2016.01.014
Hilgers R., Vincken J.-P., Gruppen H., Kabel M. A. (2018). Laccase/mediator systems: Their reactivity toward phenolic lignin structures. ACS Sustain. Chem. Eng. 6 (2), 2037–2046. doi: 10.1021/acssuschemeng.7b03451
Hobbs J. K., Shepherd C., Saul D. J., Demetras N. J., Haaning S., Monk C. R., et al. (2012). On the origin and evolution of thermophily: Reconstruction of functional precambrian enzymes from ancestors of Bacillus. Mol. Biol. Evol. 29 (2), 825–835. doi: 10.1093/molbev/msr253
Huang R., Hippauf F., Rohrbeck D., Haustein M., Wenke K., Feike J., et al. (2012). Enzyme functional evolution through improved catalysis of ancestrally nonpreferred substrates. Proc. Natl. Acad. Sci. 109 (8), 2966 LP–2971. doi: 10.1073/pnas.1019605109
Ishida C., Miyata R., Hasebe F., Miyata A., Kumazawa S., Ito S., et al. (2021). Reconstruction of hyper-thermostable ancestral l-amino acid oxidase to perform deracemization to d-amino acids. ChemCatChem 13 (24), 5228–5235. doi: 10.1002/cctc.202101296
Jaeger K.-E., Eggert T. (2002). Lipases for biotechnology. Curr. Opin. Biotechnol. 13 (4), 390–397. doi: 10.1016/s0958-1669(02)00341-5
James T. Y., Kauff F., Schoch C. L., Matheny P. B., Hofstetter V., Cox C. J., et al. (2006). Reconstructing the early evolution of fungi using a six-gene phylogeny. Nature 443 (7113), 818–822. doi: 10.1038/nature05110
Janusz G., Pawlik A., Świderska-Burek U., Polak J., Sulej J., Jarosz-Wilkołazka A., et al. (2020). Laccase properties, physiological functions, and evolution. Int. J. Mol. Sci. 21 (3), 966. doi: 10.3390/ijms21030966
Jemth P., Karlsson E., Vögeli B., Guzovsky B., Andersson E., Hultqvist G., et al. (2018). Structure and dynamics conspire in the evolution of affinity between intrinsically disordered proteins. Sci. Adv 4, 1–13. doi: 10.1126/sciadv.aau4130
Katoh K., Misawa K., Kuma K., Miyata T. (2002). MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30 (14), 3059–3066. doi: 10.1093/nar/gkf436
Kovacic F., Babic N., Krauss U., Jaeger K.-E. (2018). Classification of lipolytic enzymes from bacteria BT - aerobic utilization of hydrocarbons, oils and lipids. Ed. Rojo F. (Cham, Switzerland: Springer International Publishing), 1–35. doi: 10.1007/978-3-319-39782-5_39-1
Liberles D. A. (2007). Ancestral sequence reconstruction (Oxford, England: Oxford University Press on Demand).
Li X., Dilokpimol A., Kabel M. A., de Vries R. P. (2022). Fungal xylanolytic enzymes: Diversity and applications. Bioresour. Technol. 344, 126290. doi: 10.1016/j.biortech.2021.126290
Linde D., Ayuso-Fernández I., Laloux M., Aguiar-Cervera J. E., de Lacey A. L., Ruiz-Dueñas F. J., et al. (2021). Comparing ligninolytic capabilities of bacterial and fungal dye-decolorizing peroxidases and class-II peroxidase-catalases. Int. J. Mol. Sci. 22 (5), 1–23. doi: 10.3390/ijms22052629
Linde D., Ruiz-Dueñas F. J., Fernández-Fueyo E., Guallar V., Hammel K. E., Pogni R., et al. (2015). Basidiomycete DyPs: Genomic diversity, structural–functional aspects, reaction mechanism and environmental significance. Arch. Biochem. Biophysics. 574, 66–74. doi: 10.1016/j.abb.2015.01.018
Liu Y., Wang P., Tian J., Seidi F., Guo J., Zhu W., et al. (2022). Carbohydrate-binding modules of potential resources: Occurrence in nature, function, and application in fiber recognition and treatment. Polymers, 14(9). doi: 10.3390/polym14091806
Lombard V., Golaconda Ramulu H., Drula E., Coutinho P. M., Henrissat B. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 42 (Database issue), D490–D495. doi: 10.1093/nar/gkt1178
Löytynoja A. (2014). Phylogeny-aware alignment with PRANK. Methods Mol. Biol. (Clifton. N.J.). (1079), 155–170. doi: 10.1007/978-1-62703-646-7_10
Malcolm B. A., Wilson K. P., Matthews B. W., Kirsch J. F., Wilson A. C. (1990). Ancestral lysozymes reconstructed, neutrality tested, and thermostability linked to hydrocarbon packing. Nature 345 (6270), 86–89. doi: 10.1038/345086a0
Mancheño J. M., Pernas M. A., Martínez M. J., Ochoa B., Rúa M. L., Hermoso J. A. (2003). Structural insights into the lipase/esterase behavior in the Candida rugosa lipases family: Crystal structure of the lipase 2 isoenzyme at 1.97Å resolution. J. Mol. Biol. 332 (5), 1059–1069. doi: 10.1016/j.jmb.2003.08.005
Marzullo L., Cannio R., Giardina P., Santini M. T., Sannia G. (1995). Veratryl alcohol oxidase from Pleurotus ostreatus participates in lignin biodegradation and prevents polymerization of laccase-oxidized substrates. J. Biol. Chem. 270 (8), 3823–3827. doi: 10.1074/jbc.270.8.3823
Maté D., García-Burgos C., García-Ruiz E., Ballesteros A. O., Camarero S., Alcalde M. (2010). Laboratory evolution of high-redox potential laccases. Chem. Biol. 17 (9), 1030–1041. doi: 10.1016/j.chembiol.2010.07.010
Merkl R., Sterner R. (2016). Ancestral protein reconstruction: Techniques and applications. Biol. Chem. 397 (1), 1–21. doi: 10.1515/hsz-2015-0158
Mester T., Ambert-Balay K., Ciofi-Baffoni S., Banci L., Jones A. D., Tien M. (2001). Oxidation of a tetrameric nonphenolic lignin model compound by lignin peroxidase. J. Biol. Chem. 276 (25), 22985–22990. doi: 10.1074/jbc.M010739200
Min K., Gong G., Woo H. M., Kim Y., Um Y. (2015). A dye-decolorizing peroxidase from Bacillus subtilis exhibiting substrate-dependent optimum temperature for dyes and β-ether lignin dimer. Sci. Rep. 5 (1), 8245. doi: 10.1038/srep08245
Monika T., A. H. O., Heidi Ø., Anikó V., Francisco V., H. E. V. G., et al. (2022). Comparison of six lytic polysaccharide monooxygenases from Thermothielavioides terrestris shows that functional variation underlies the multiplicity of LPMO genes in filamentous fungi. Appl. Environ. Microbiol. 88 (6), e00096–e00022. doi: 10.1128/aem.00096-22
Morris J. L., Puttick M. N., Clark J. W., Edwards D., Kenrick P., Pressel S., et al. (2018). The timescale of early land plant evolution. Proc. Natl. Acad. Sci. United. States America 115 (10), E2274–E2283. doi: 10.1073/pnas.1719588115
Musil M., Khan R. T., Beier A., Stourac J., Konegger H., Damborsky J., et al. (2021). FireProtASR: A web server for fully automated ancestral sequence reconstruction. Briefings Bioinf. 22 (4), bbaa337. doi: 10.1093/bib/bbaa337
Nagy L. G., Kovács G. M., Krizsán K. (2018). Complex multicellularity in fungi: evolutionary convergence, single origin, or both? Biol. Rev. 93 (4), 1778–1794. doi: 10.1111/brv.12418
Nagy L. G., Riley R., Bergmann P. J., Krizsán K., Martin F. M., Grigoriev I. V., et al. (2017). Genetic bases of fungal white-rot wood decay predicted by phylogenomic analysis of correlated gene-phenotype evolution. Mol. Biol. Evol. 34 (1), 35–44. doi: 10.1093/molbev/msw238
Nagy L. G., Riley R., Tritt A., Adam C., Daum C., Floudas D., et al. (2016). Comparative genomics of early-diverging mushroom-forming fungi provides insights into the origins of lignocellulose decay capabilities. Mol. Biol. Evol. 33 (4), 959–970. doi: 10.1093/molbev/msv337
Nelsen M. P., DiMichele W. A., Peters S. E., Boyce C. K. (2016). Delayed fungal evolution did not cause the Paleozoic peak in coal production. Proc. Natl. Acad. Sci. 113 (9), 2442 LP–2447. doi: 10.1073/pnas.1517943113
Nickle D. C., Jensen M. A., Gottlieb G. S., Shriner D., Learn G. H., Rodrigo A. G., et al. (2003). “Consensus and ancestral state HIV vaccines,” in Science Vol. 299, Issue 5612. (New York, N.Y: Science), 1515–1518). doi: 10.1126/science.299.5612.1515c
Nicoll C. R., Bailleul G., Fiorentini F., Mascotti M. L., Fraaije M. W., Mattevi A. (2020). Ancestral-sequence reconstruction unveils the structural basis of function in mammalian FMOs. Nat. Struct. Mol. Biol. 27 (1), 14–24. doi: 10.1038/s41594-019-0347-2
Parfrey L. W., Lahr D. J. G., Knoll A. H., Katz L. A. (2011). Estimating the timing of early eukaryotic diversification with multigene molecular clocks. Proc. Natl. Acad. Sciences. 108(33). 13624. LP. –, 13629. doi: 10.1073/pnas.1110633108
Patel A. K., Singhania R. R., Sim S. J., Pandey A. (2019). Thermostable cellulases: Current status and perspectives. Bioresour. Technol. 279, 385–392. doi: 10.1016/j.biortech.2019.01.049
Pauling L., Zuckerkandl E., Henriksen T., Lövstad R. (1963). Chemical paleogenetics. molecular restoration studies of extinct forms of life. Acta chemica scandinavica. 17 supl., 9–16. doi: 10.3891/acta.chem.scand.17s-0009
Payne C. M., Knott B. C., Mayes H. B., Hansson H., Himmel M. E., Sandgren M., et al. (2015). Fungal cellulases. Chem. Rev. 115 (3), 1308–1448. doi: 10.1021/cr500351c
Percival Zhang Y.-H., Himmel M. E., Mielenz J. R. (2006). Outlook for cellulase improvement: screening and selection strategies. Biotechnol. Adv. 24 (5), 452–481. doi: 10.1016/j.biotechadv.2006.03.003
Ragauskas A. J., Beckham G. T., Biddy M. J., Chandra R., Chen F., Davis M. F., et al. (2014). Lignin valorization: Improving lignin processing in the biorefinery. Sci. (New. York. N.Y.). 344 (6185), 1246843. doi: 10.1126/science.1246843
Rai A., Klare J. P., Reinke P. Y. A., Englmaier F., Fohrer J., Fedorov R., et al. (2021). “Structural and biochemical characterization of a dye-decolorizing peroxidase from dictyostelium discoideum,” in International journal of molecular sciences Vol. 22 Issue 12. doi: 10.3390/ijms22126265
Riley R., Salamov A., Brown D. W., Nagy L. G., Dimitrios F., Held B. W., et al. (2014). Extensive sampling of basidiomycete genomes demonstrates inadequacy of the white-rot/brown-rot paradigm for wood decay fungi. Proc. Natl. Acad. Sci. 111 (27), 9923–9928. doi: 10.1073/pnas.1400592111
Risso V. A., Gavira J. A., Sanchez-Ruiz J. M. (2014). Thermostable and promiscuous precambrian proteins. Environ. Microbiol. 16 (6), 1485–1489. doi: 10.1111/1462-2920.12319
Risso V. A., Sanchez-Ruiz J. M., Ozkan S. B. (2018). Biotechnological and protein-engineering implications of ancestral protein resurrection. Curr. Opin. Struct. Biol. 51, 106–115. doi: 10.1016/j.sbi.2018.02.007
Rodríguez Couto S., Toca Herrera J. L. (2006). Industrial and biotechnological applications of laccases: A review. Biotechnol. Adv. 24 (5), 500–513. doi: 10.1016/j.biotechadv.2006.04.003
Rodríguez-Salarichs J., García de Lacoba M., Prieto A., Martínez M. J., Barriuso J. (2021). Versatile lipases from the Candida rugosa-like family: A mechanistic insight using computational approaches. J. Chem. Inf. Modeling. 61 (2), 913–920. doi: 10.1021/acs.jcim.0c01151
Ronquist F., Teslenko M., van der Mark P., Ayres D. L., Darling A., Höhna S., et al. (2012). MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Systematic. Biol. 61 (3), 539–542. doi: 10.1093/sysbio/sys029
Ross C. M., Foley G., Boden M., Gillam E. M. J. (2022). Using the evolutionary history of proteins to engineer insertion-deletion mutants from robust, ancestral templates using graphical representation of ancestral sequence predictions (GRASP). Methods Mol. Biol. (Clifton. N.J.). 2397, 85–110. doi: 10.1007/978-1-0716-1826-4_6
Rozi M. F. A. M., Rahman R. N. Z. R. A., Leow A. T. C., Ali M. S. M. (2022). Ancestral sequence reconstruction of ancient lipase from family I.3 bacterial lipolytic enzymes. Mol. Phylogenet. Evol. 168, 107381. doi: 10.1016/j.ympev.2021.107381
Ruiz-Dueñas F. J., Barrasa J. M., Sánchez-García M., Camarero S., Miyauchi S., Serrano A., et al. (2021). Genomic analysis enlightens agaricales lifestyle evolution and increasing peroxidase diversity. Mol. Biol. Evol. 38 (4), 1428–1446. doi: 10.1093/molbev/msaa301
Ruiz-Dueñas F. J., Lundell T., Floudas D., Nagy L. G., Barrasa J. M., Hibbett D. S., et al. (2013). Lignin-degrading peroxidases in polyporales: An evolutionary survey based on 10 sequenced genomes. Mycologia 105 (6), 1428–1444. doi: 10.3852/13-059
Ruiz-Dueñas F. J., Morales M., García E., Miki Y., Martínez M. J., Martínez A. T. (2009). Substrate oxidation sites in versatile peroxidase and other basidiomycete peroxidases. J. Exp. Bot. 60 (2), 441–452. doi: 10.1093/jxb/ern261
Sørensen A., Lübeck M., Lübeck P. S., Ahring B. K. (2013). Fungal β-glucosidases: A bottleneck in industrial use of lignocellulosic materials. Biomolecules 3 (3), 612–631. doi: 10.3390/biom3030612
Sáez-Jiménez V., Rencoret J., Rodríguez-Carvajal M. A., Gutiérrez A., Ruiz-Dueñas F. J., Martínez A. T. (2016). Role of surface tryptophan for peroxidase oxidation of nonphenolic lignin. Biotechnol. Biofuels 9 (1), 198. doi: 10.1186/s13068-016-0615-x
Saha B. C. (2003). Hemicellulose bioconversion. J. Ind. Microbiol. Biotechnol. 30 (5), 279–291. doi: 10.1007/s10295-003-0049-x
Sánchez-Ruiz M. I., Ayuso-Fernández I., Rencoret J., González-Ramírez A. M., Linde D., Davó-Siguero I., et al. (2021). “Agaricales mushroom lignin peroxidase: From structure–function to degradative capabilities,” in Antioxidants (Basel, Switzerland:MDPI) Vol. 10, Issue 9). doi: 10.3390/antiox10091446
Sarmah N., Revathi D., Sheelu G., Yamuna Rani K., Sridhar S., Mehtab V., et al. (2018). Recent advances on sources and industrial applications of lipases. Biotechnol. Prog. 34 (1), 5–28. doi: 10.1002/btpr.2581
Savinova O. S., Moiseenko K. V., Vavilova E. A., Chulkin A. M., Fedorova T. V., Tyazhelova T. V., et al. (2019). “Evolutionary relationships between the laccase genes of polyporales: Orthology-based classification of laccase isozymes and functional insight from trametes hirsuta,” in Frontiers in microbiology (Lausanne, Switzerland:Frontiers Media SA) Vol. 10. doi: 10.3389/fmicb.2019.00152
Schilling J. S., Kaffenberger J. T., Held B. W., Ortiz R., Blanchette R. A. (2020). “Using wood rot phenotypes to illuminate the “gray” among decomposer fungi,” in Frontiers in microbiology (Lausanne, Switzerland: Frontiers Media SA) Vol. 11. doi: 10.3389/fmicb.2020.01288
Selberg A. G. A., Gaucher E. A., Liberles D. A. (2021). Ancestral sequence reconstruction: From chemical paleogenetics to maximum likelihood algorithms and beyond. J. Mol. Evol. 89 (3), 157–164. doi: 10.1007/s00239-021-09993-1
Singh R., Eltis L. D. (2015). The multihued palette of dye-decolorizing peroxidases. Arch. Biochem. Biophysics. 574, 56–65. doi: 10.1016/j.abb.2015.01.014
Spence M. A., Kaczmarski J. A., Saunders J. W., Jackson C. J. (2021). Ancestral sequence reconstruction for protein engineers. Current opinion in structural biology. 69, 131–141. doi: 10.1016/j.sbi.2021.04.001
Srebotnik E., Boisson J.-N. (2005). Peroxidation of linoleic acid during the oxidation of phenols by fungal laccase. Enzyme Microbial. Technol. 36 (5), 785–789. doi: 10.1016/j.enzmictec.2005.01.004
Stackhouse J., Presnell S. R., McGeehan G. M., Nambiar K. P., Benner S. A. (1990). The ribonuclease from an extinct bovid ruminant. FEBS Lett. 262 (1), 104–106. doi: 10.1016/0014-5793(90)80164-E
Sützl L., Foley G., Gillam E. M. J., Bodén M., Haltrich D. (2019). The GMC superfamily of oxidoreductases revisited: Analysis and evolution of fungal GMC oxidoreductases. Biotechnol. Biofuels 12 (1), 118. doi: 10.1186/s13068-019-1457-0
Tamura K., Stecher G., Kumar S. (2021). MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 38 (7), 3022–3027. doi: 10.1093/molbev/msab120
Tao Q., Barba-Montoya J., Huuki L. A., Durnan M. K., Kumar S. (2020). Relative efficiencies of simple and complex substitution models in estimating divergence times in phylogenomics. Mol. Biol. Evol. 37 (6), 1819–1831. doi: 10.1093/molbev/msaa049
Trudeau D. L., Kaltenbach M., Tawfik D. S. (2016). On the potential origins of the high stability of reconstructed ancestral proteins. Mol. Biol. Evol. 33 (10), 2633–2641. doi: 10.1093/molbev/msw138
Vaaje-Kolstad G., Westereng B., Horn S. J., Liu Z., Zhai H., Sørlie M., et al. (2010). An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Sci. (New. York. N.Y.). 330 (6001), 219–222. doi: 10.1126/science.1192231
van den Brink J., de Vries R. P. (2011). Fungal enzyme sets for plant polysaccharide degradation. Appl. Microbiol. Biotechnol. 91 (6), 1477–1492. doi: 10.1007/s00253-011-3473-2
Vanholme R., Demedts B., Morreel K., Ralph J., Boerjan W. (2010). Lignin biosynthesis and structure. Plant Physiol. 153 (3), 895–905. doi: 10.1104/pp.110.155119
Vaquero M. E., Barriuso J., Martínez M. J., Prieto A. (2016). Properties, structure, and applications of microbial sterol esterases. Appl. Microbiol. Biotechnol. 100 (5), 2047–2061. doi: 10.1007/s00253-015-7258-x
Vialle R. A., Tamuri A. U., Goldman N. (2018). Alignment modulates ancestral sequence reconstruction accuracy. Mol. Biol. Evol. 35 (7), 1783–1797. doi: 10.1093/molbev/msy055
Voordeckers K., Brown C. A., Vanneste K., van der Zande E., Voet A., Maere S., et al. (2012). Reconstruction of ancestral metabolic enzymes reveals molecular mechanisms underlying evolutionary innovation through gene duplication. PloS Biol. 10 (12), e1001446. doi: 10.1371/journal.pbio.1001446
Weng J.-K., Chapple C. (2010). The origin and evolution of lignin biosynthesis. New Phytol. 187 (2), 273–285. doi: 10.1111/j.1469-8137.2010.03327.x
Wheeler L. C., Lim S. A., Marqusee S., Harms M. J. (2016). The thermostability and specificity of ancient proteins. Curr. Opin. Struct. Biol. 38, 37–43. doi: 10.1016/j.sbi.2016.05.015
Williams P. D., Pollock D. D., Blackburne B. P., Goldstein R. A. (2006). Assessing the accuracy of ancestral protein reconstruction methods. PloS Comput. Biol. 2 (6), e69. doi: 10.1371/journal.pcbi.0020069
Wilson C., Agafonov R. V., Hoemberger M., Kutter S., Zorba A., Halpin J., et al. (2015). Kinase dynamics. using ancient protein kinases to unravel a modern cancer drug’s mechanism. Science 347 (6224), 882–886. doi: 10.1126/science.aaa1823
Wong D. W. S. (2006). Feruloyl esterase. Appl. Biochem. Biotechnol. 133 (2), 87–112. doi: 10.1385/ABAB:133:2:87
Xu F. (1996). Oxidation of phenols, anilines, and benzenethiols by fungal laccases: correlation between activity and redox potentials as well as halide inhibition. Biochemistry 35 (23), 7608–7614. doi: 10.1021/bi952971a
Yang Z. (2007). PAML 4: Phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24 (8), 1586–1591. doi: 10.1093/molbev/msm088
Zámocký M., Hofbauer S., Schaffner I., Gasselhuber B., Nicolussi A., Soudi M., et al. (2015). Independent evolution of four heme peroxidase superfamilies. Arch. Biochem. Biophys. 574, 108–119. doi: 10.1016/j.abb.2014.12.025
Zhao R.-L., Li G.-J., Sánchez-Ramírez S., Stata M., Yang Z.-L., Wu G., et al. (2017). A six-gene phylogenetic overview of basidiomycota and allied phyla with estimated divergence times of higher taxa and a phyloproteomics perspective. Fungal Diversity 84 (1), 43–74. doi: 10.1007/s13225-017-0381-5
Keywords: ancestral sequence reconstruction, wood decay fungi, lignocellulosic biomass, plant cell-wall degrading enzymes (PCWDE), evolution
Citation: Ayuso-Fernández I, Molpeceres G, Camarero S, Ruiz-Dueñas FJ and Martínez AT (2022) Ancestral sequence reconstruction as a tool to study the evolution of wood decaying fungi. Front. Fungal Biol. 3:1003489. doi: 10.3389/ffunb.2022.1003489
Received: 26 July 2022; Accepted: 22 September 2022;
Published: 14 October 2022.
Edited by:
Donald O. Natvig, University of New Mexico, United StatesReviewed by:
Carolina Elena Girometta, University of Pavia, ItalyMiguel Arenas, University of Vigo, Spain
Copyright © 2022 Ayuso-Fernández, Molpeceres, Camarero, Ruiz-Dueñas and Martínez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Iván Ayuso-Fernández, ivan.ayuso-fernandez@nmbu.no; Angel T. Martínez, atmartinez@cib.csic.es
†These authors share first authorship