Comparative Genomics and Transcriptomics Analyses Reveal Divergent Plant Biomass-Degrading Strategies in Fungi

Plant biomass is one of the most abundant renewable carbon sources, which holds great potential for replacing current fossil-based production of fuels and chemicals. In nature, fungi can efficiently degrade plant polysaccharides by secreting a broad range of carbohydrate-active enzymes (CAZymes), such as cellulases, hemicellulases, and pectinases. Due to the crucial role of plant biomass-degrading (PBD) CAZymes in fungal growth and related biotechnology applications, investigation of their genomic diversity and transcriptional dynamics has attracted increasing attention. In this project, we systematically compared the genome content of PBD CAZymes in six taxonomically distant species, Aspergillus niger, Aspergillus nidulans, Penicillium subrubescens, Trichoderma reesei, Phanerochaete chrysosporium, and Dichomitus squalens, as well as their transcriptome profiles during growth on nine monosaccharides. Considerable genomic variation and remarkable transcriptomic diversity of CAZymes were identified, implying the preferred carbon source of these fungi and their different methods of transcription regulation. In addition, the specific carbon utilization ability inferred from genomics and transcriptomics was compared with fungal growth profiles on corresponding sugars, to improve our understanding of the conversion process. This study enhances our understanding of genomic and transcriptomic diversity of fungal plant polysaccharide-degrading enzymes and provides new insights into designing enzyme mixtures and metabolic engineering of fungi for related industrial applications.


Introduction
Plant biomass is a magnificent renewable carbon source on Earth, which is of major importance for ecology and is an attractive substrate for the production of biofuels and biochemicals. Plant cell walls mainly consist of polysaccharides (e.g., cellulose, hemicellulose, and pectin), lignin, and small amounts of protein. In addition, starch and inulin are commonly found as storage polysaccharides.
In nature, fungi play a central role in the degradation of plant biomass due to their unique ability to secrete a broad range of extracellular enzymes for decomposition of recalcitrant plant polysaccharides into mono-or oligosaccharides that can be metabolized for their growth and reproduction. Based on sequence and structural similarity, carbohydrate-active enzymes are catalogued into five major families in the Carbohydrate-Active enZymes (CAZymes) database (http://www.cazy.org/ (accessed on 24 February 2023)) [1,2], including glycoside hydrolases (GHs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), glycosyltransferases (GTs), and auxiliary activities (AAs). Fungal plant biomass

Phylogenetic Analyses
Protein sequences of each selected CAZy family (PBD-related) were aligned using MAFFT v7.310 [29,30]. Unusually long and incomplete sequences were corrected manually based on NCBI BlastX searching results, whereas duplicate and ambiguous sequences were discarded. Phylogenetic trees generated using FastTree [23] were used to distinguish true orthologs across six fungi. The resulting trees were visualized using iTOL [31].
The culture conditions were maintained as similarly as possible to allow a good comparison of the species. In detail, A. niger NRRL3, A. nidulans FGSC A4, and P. subrubescens FBCC1632/CBS132785 were pre-cultured in complete medium [32] with 2% D-fructose, and mycelial aliquots were then transferred to minimal medium [32] with 25 mM D-glucose, D-fructose, D-galactose, D-mannose, L-rhamnose, D-xylose, L-arabinose, D-galacturonic acid, and D-glucuronic acid, respectively, and cultivated for 2 h. The same cultivation approach was used for T. reesei QM6a and P. chrysosporium PR-78, but with media optimized for these species [33,34] and cultivated for 4 h. The cultivation approach for D. squalens can be found in our previous study [18]. So, only the basic salt medium is specific for the species, but the carbon source is the same in all species. Mycelial samples were harvested and immediately frozen in liquid nitrogen.
Total RNA was extracted from ground mycelial samples using TRIzol reagent (Invitrogen). Purification of mRNA, synthesis of cDNA library, and sequencing were conducted at DOE Joint Genome Institute (JGI). The details of transcriptome data generation and related downstream analysis were described previously [24].
In addition, the CAZy genes with strong sugar specificity were selected using a sugar specificity index (SSI), which was calculated by following the equation that has been widely used in studies of gene expression specificity [35,36]: For the above equation, x i is the expression of the gene in each specific monosaccharide condition i, and n is the number of total tested monosaccharide conditions. We define the CAZyme genes with SSI > 0.7 as sugar-specific genes (SSGs). In addition, we filtered out the extremely low-expressed genes whose maximum expression levels (FPKM) were less than 10 in all tested conditions. For each SSG, we defined the corresponding sugar growth condition in which the gene showed highest expression level among all tested sugars as its inducing sugar.

Growth Profiling on 18 Plant Biomass-Related Substrates
For growth phenotype analyses, all strains were grown on minimal medium (MM) [32] with monosaccharides, disaccharides at 25 mM, and polysaccharides at 1% final concentration. Growth was performed at 30 • C for two Aspergillus spp. and 25 • C for the other species. Media with no carbon source was used as a control. If growth on a specific carbon source is the same as on no carbon source, it is considered as no growth. Growth on D-glucose was used as an internal reference for growth. Substrates were obtained from Sigma-Aldrich (

Genomic Potential of PBD CAZymes Shows Diversity across Six Fungi
In total, the genome repertoire of 65 CAZy (sub-)families and feruloyl esterases (FAEs) involved in degradation of plant polysaccharides was comparatively analyzed in six fungi (Figures 1 and S1). Specifically, the selected PBD-related enzymes studied included four auxiliary activities (AAs) families (AA3, AA9, AA13, and AA16), six carbohydrate esterase families (CE1, CE5, CE8, CE12, CE15, and CE16), five polysaccharide lyase families (PL1, PL3, PL4, PL9, and PL11), 39 glycoside hydrolase (GH) families, and FAEs ( Figure 1). Notably, ten enzyme families (GH3, GH12, CE16, GH2, GH35, GH31, GH43, GH51, GH54, and FAE) are multifunctional and involved in the degradation of more than one plant polysaccharide ( Figure S1). The potential for PBD as inferred from the genomic repertoire of PBD CAZymes shows clear diversity across the studied fungi ( Figure 1 and Table S2). Overall, the CA-Zyme repertoire differences correlate well with the taxonomic distance of the studied The potential for PBD as inferred from the genomic repertoire of PBD CAZymes shows clear diversity across the studied fungi ( Figure 1 and Table S2). Overall, the CAZyme repertoire differences correlate well with the taxonomic distance of the studied fungi. The three Eurotiomycetes (A. niger, A. nidulans, and P. subrubescens) possess a relatively high number of PBD-related CAZymes compared to T. reesei and the Basidiomycetes (P. chrysosporium and D. squalens), especially with respect to enzymes from GH and PL families ( Figure S2 and Table S3). T. reesei contains fewer CAZymes involved in each plant polysaccharide than the other three Ascomycetes. The two Basidiomycetes have fewer hemicellulolytic, pectinolytic, and starch-degrading enzymes, while the total numbers of cellulolytic enzymes are comparable to Eurotiomycetes. Notably, a higher number of AA9 genes were identified in the Basidiomycetes. In the following sections, we discuss in more detail the diversity of genomic potential with respect to degradation of each polysaccharide.
EGLs act mostly on amorphous regions of cellulose, releasing reducing and nonreducing chain ends [38,39]. They are mainly classified in GH12, GH131, and GH45 and three subfamilies of GH5 (GH5_4, GH5_5, and GH5_22) ( Table S3). Many of these EGL families showed noticeable differences between both evolutionarily close and distant species, especially the GH5 (sub-)family. GH5_4 has a single copy in A. nidulans and P. subrubescens but is missing in all other species. GH5_22 was not identified in the two Aspergillus species and T. reesei but is present as two copies in P. subrubescens and the two Basidiomycetes. P. subrubescens has more copies of GH5_5 and GH12 than the other species. In contrast to most EGLs with comparable or more copies in Ascomycetes than Basidiomycetes, GH45 and GH131 have multiple copies in both Basidiomycetes, but are present as a single copy in three Ascomycetes (except for the absence of GH131 in T. reesei).
CBHs release cellobiose from the reducing (CBH I) or non-reducing (CBH II) ends of the cellulose. Fungal CBH I and CBH II genes are classified in GH7 and GH6, respectively. Gene copies of GH7 are notably larger (eight) in P. chrysosporium and vary among the other species (two to four). GH6 was present as a single copy in T. reesei and the Basidiomycetes and with two genes in the three Eurotiomycetes.
BGLs are involved in the cleavage of cellobiose into glucose and are mainly grouped into the GH1 and GH3 families. GH3 BGLs are greatly expanded in the three Eurotiomycetes, while GH1 is only significantly expanded in P. subrubescens.
LPMOs and cellobiose dehydrogenases (CDHs) have been suggested to be involved in the oxidative cleavage of cellulose and boosted cellulose degradation [26,40,41]. AA9 is a widely distributed LPMO in fungi, which was present as multiple copies in all studied fungi and was significantly expanded in the Basidiomycetes compared to the Ascomycetes. In contrast, AA16 LPMO was only present as a single copy in the Eurotiomycetes and was missing in the other species, while AA3_1 CDH showed less variation among the studied fungi.

Hemicellulose Degradation
In contrast to cellulose, hemicellulose is a heterogeneous polysaccharide, which is composed of different residues that form diverse backbones and branches. Degradation of plant hemicelluloses requires a variety of hemicellulolytic enzymes, which are divided over at least twenty GH families, four CE families, and FAE (Figures 1 and S1 and Table S1).
Overall, the genomes of A. nidulans and P. subrubescens contain more copies of mannanases than the other fungi ( Figure 1). The differences are that A. nidulans has more copies of MAN (GH26, GH134, and GH5_7), whereas P. subrubescens harbors large numbers of AGL (GH27 and GH36)-and MND (GH2)-encoding genes. In general, the other four species have fewer mannanases, except that T. reesei and A. niger contain considerable numbers of AGLs (especially GH27) and MNDs, and the Basidiomycetes have a higher number of CE16 copies. Notably, GH36, GH26, and GH134 were absent in both Basidiomycetes (Figure 1), and GH26 and GH134 were absent in the Basidiomycetes and T. reesei.

Pectin Degradation
Pectin is a very complex polymer, composed of a family of galacturonic acid-rich polysaccharides, including: homogalacturonan (HG), xylogalacturonan (XGA), and two types of rhamnogalacturonan (RG-I and RG-II) [50]. Its degradation requires a variety of enzymes (Table S1), which are classified into more than 20 CAZy families, mainly including pectin hydrolases (GH2, GH5_16, GH28, GH35, GH43, GH51, GH54, GH78, GH88, GH93, and GH105), pectin lyases (PL1, PL3, PL4, PL9, and PL11) [51], pectin methylesterases (CE8), and pectin/rhamnogalacturonan acetyl esterase (CE12) [52]. The gene content of pectinolytic genes showed huge differences among the six studied species. In general, the Eurotiomycete fungi have a higher number and more diverse set of pectinolytic genes than the other three studied species (Figure 1). Several enzyme families involved in pectin degradation are reduced or absent in T. reesei, including the absence of four GH families (GH51, GH53, GH88, and GH93) and many fewer genes encoding for GH28 and GH78 compared to the other three Ascomycetes. PLs are exclusively present and largely varied in three Eurotiomycota species, of which A. nidulans has more types of PL families (e.g., PL1, PL3, PL4, PL9, and PL11). The other studied species lost most of the PLs, except D. squalens, which has one gene of PL4_3. Regarding pectin methylesterases, CE8 was found in multiple copies in most species except in T. reesei, while CE12 was present in the three Eurotiomycota species and D. squalens and was absent in T. reesei and P. chrysosporium.
Inulin consists of a branched β-2,1-linked chain of D-fructose with a terminal Dglucose residue [56,57]. GH32 is the only CAZy family involved in inulin degradation, and it mainly acts as endo-inulinase, exo-inulinase, and invertase. GH32 enzymes were absent in T. reesei and P. chrysosporium but were expanded in A. niger and P. subrubescens with five and nine copies, respectively (Figure 1).

Expression Profile of CAZy Genes during Fungi Grown on Different Monosaccharides
PBD CAZymes mediate the release of monomeric/oligo-sugars from plant polysaccharides. In turn, the presence of these small sugars can induce or repress the expression of CAZyme-encoding genes. Here, we investigated the transcriptome diversity of CAZymes in response to different plant-derived monosaccharides of studied fungi.
Overall, the expression profiles of CAZy genes varied across the tested species and enzyme families. The inducing pattern of pectinolytic genes showed considerable consistency across the studied species (Figure 2). For most of the studied Ascomycetes, relative high abundance of pectinolytic genes were observed on L-arabinose, D-galacturonic acid, D-glucuronic acid, or L-rhamnose, which are the most abundant sugars present in the polymers of pectin. A. nidulans is an exception, in which pectinolytic genes mainly had higher expression on L-rhamnose and were not clearly induced by D-galacturonic acid ( Figure 2 and Table S4). Additionally, the genes involved in the degradation of xylan and xyloglucan showed relatively high expression on L-arabinose in the Ascomycetes, particularly P. subrubescens, which had more highly expressed genes on L-arabinose, while they were also induced by D-xylose in A. niger and induced by L-rhamnose in A. nidulans. Notably, no clear sugar-inducing pattern was observed for the studied Basidiomycetes, except that most PBD-related CAZymes showed relatively high expression levels on D-galacturonic acid compared to the other tested sugars for P. chrysosporium.
In addition, we also analyzed the sugar specificity index of each CAZy gene. Sugar specificity is often described based on gene expression levels, and we define genes with a high sugar specificity index as a group of genes whose function and expression are specific for one certain sugar condition. This index can be used to evaluate if the specific CAZy genes involved in certain polysaccharide degradation showed corresponding inducing patterns on the corresponding monosaccharides ( Figure S3 and Table S5). As expected, in the three Eurotiomycetes, the CAZy genes involved in pectin degradation showed clear D-galacturonic acid-and L-rhamnose-inducing specificity (except mainly L-rhamnoseinducing patterns in A. nidulans), and the CAZy genes involved in xylan or xyloglucan degradation revealed L-arabinose-and D-xylose-inducing patterns. Notably, parts of CAZy genes involved in the degradation of the above polysaccharides also showed high expression levels on D-galactose in A. niger and A. nidulans, but not in P. subrubescens, which indicates metabolic cross talk between sugar-sensing pathways and/or overlap in regulatory networks. In addition, ten and eighteen CAZymes involved in pectin, xylan, or mannan degradation showed sugar-specific induction on D-glucuronic acid in A. nidulans and P. subrubescens, respectively. In contrast, T. reesei only had a small number of CAZy genes with a high sugar specificity index, such as six CAZy genes involved in xylan and xyloglucan that showed L-arabinose induction and three pectinolytic genes with Dgalacturonic acid induction specificity (Table S5). We detected more than twenty CAZymes highly expressed on D-galacturonic acid in P. chrysosporium, but only nine genes belong to pectinolytic genes. In D. squalens, due to the limited transcriptome data, we only observed a very small number of genes showing sugar specificity ( Figure S3 and Table S5).
PBD CAZymes mediate the release of monomeric/oligo-sugars from plant polysaccharides. In turn, the presence of these small sugars can induce or repress the expression of CAZyme-encoding genes. Here, we investigated the transcriptome diversity of CA-Zymes in response to different plant-derived monosaccharides of studied fungi.
Overall, the expression profiles of CAZy genes varied across the tested species and enzyme families. The inducing pattern of pectinolytic genes showed considerable consistency across the studied species (Figure 2). For most of the studied Ascomycetes, relative high abundance of pectinolytic genes were observed on L-arabinose, D-galacturonic acid, D-glucuronic acid, or L-rhamnose, which are the most abundant sugars present in the polymers of pectin. A. nidulans is an exception, in which pectinolytic genes mainly had higher expression on L-rhamnose and were not clearly induced by D-galacturonic acid ( Figure 2 and Table S4). Additionally, the genes involved in the degradation of xylan and xyloglucan showed relatively high expression on L-arabinose in the Ascomycetes, particularly P. subrubescens, which had more highly expressed genes on L-arabinose, while they were also induced by D-xylose in A. niger and induced by L-rhamnose in A. nidulans. Notably, no clear sugar-inducing pattern was observed for the studied Basidiomycetes, except that most PBD-related CAZymes showed relatively high expression levels on Dgalacturonic acid compared to the other tested sugars for P. chrysosporium.

Growth Profiles Do Not Always Correlate with the (Post-)Genomics Profiles of CAZy Genes
The noticeable diversity of genomics and expressional profiles of CAZy genes identified above in the tested fungi strongly suggests that they employ different PBD strategies. To evaluate the correlation between (post-)genomic profiles and actual carbon utilization ability, we analyzed the growth profile of these fungi on nine monosaccharides, two disaccharides, and seven polysaccharides ( Figure 3A). Growth on D-glucose and no carbon source were used as an internal reference for growth. The relative difference between growth on different carbon sources was examined between the species. In general, the Eurotiomycetes showed better growth on most of the tested polysaccharides than the other species, which is consistent with a higher gene content of the corresponding degrading genes ( Figure 3B). In contrast, possibly due to the absence and reduction of many CAZy genes in the genome, T. reesei grew relatively poorly on most of the tested polysaccharides (except starch). Different species showed huge variation of growth on inulin. A. niger and P. subrubescens grew well on inulin, whereas the other four species grew very poorly. This correlates well with the diversity of inulin-degrading genes in different species. A. niger and P. subrubescens have many copies of inulinolytic genes, in contrast to the other species that are impaired in their inulin-degrading machinery. Compared to the other two analyzed Eurotiomycetes, A. nidulans had fewer copies of GH32, and it was particularly missing copies of endo-inulinase (Table S3), which reduces its ability to efficiently utilize this substrate as a carbon source. In T. reesei and P. chrysosporium, GH32 was fully absent. In D. squalens, two copies of GH32 were identified, but their expression was extremely low (Table S4).
Besides these consistent observations for genomic repertoire, expression, and growth, inconsistencies were also observed. Good growth on citrus pectin was observed for T. reesei and P. chrysosporium, while their genomes encode fewer pectinolytic genes In general, the Eurotiomycetes showed better growth on most of the tested polysaccharides than the other species, which is consistent with a higher gene content of the corresponding degrading genes ( Figure 3B). In contrast, possibly due to the absence and reduction of many CAZy genes in the genome, T. reesei grew relatively poorly on most of the tested polysaccharides (except starch). Different species showed huge variation of growth on inulin. A. niger and P. subrubescens grew well on inulin, whereas the other four species grew very poorly. This correlates well with the diversity of inulin-degrading genes in different species. A. niger and P. subrubescens have many copies of inulinolytic genes, in contrast to the other species that are impaired in their inulin-degrading machinery. Compared to the other two analyzed Eurotiomycetes, A. nidulans had fewer copies of GH32, and it was particularly missing copies of endo-inulinase (Table S3), which reduces its ability to efficiently utilize this substrate as a carbon source. In T. reesei and P. chrysosporium, GH32 was fully absent. In D. squalens, two copies of GH32 were identified, but their expression was extremely low (Table S4).
Besides these consistent observations for genomic repertoire, expression, and growth, inconsistencies were also observed. Good growth on citrus pectin was observed for T. reesei and P. chrysosporium, while their genomes encode fewer pectinolytic genes than the Eurotiomycetes. Additionally, we observed active gene expression of pectinolytic genes for these two species on D-galacturonic acid, the major sugar component of pectin (Figures 2 and 3C). In contrast, D. squalens grew poorly on pectin. In line with the crucial role of pectin lyases (PLs) on degrading apple pectin [58], the absence of PLs on the genome of T. reesei and P. chrysosporium could explain their relatively poor growth on apple pectin. All tested fungi could grow relatively well on starch, while none of the studied species grew well on cellulose as the sole carbon source despite a significant number of cellulases encoded in their genomes. The Eurotiomycetes grew well on xylan, while the other three tested species showed only minimal growth on xylan, which is consistent with the fact that fewer xylanase genes were identified in their genomes and actively expressed during growth of these fungi on xylan-related monosaccharides (Figure 3). In contrast, we observed poor correlation between (galacto-)mannan (guar gum) and mannan-degrading CAZyme content in the fungal genomes. The Eurotiomycetes grew well on guar gum, even though A. niger has a very small number of mannan-degrading genes. Compared to A. niger, T. reesei and the Basidiomycetes have a comparable number of mannan-degrading genes, but they all showed impaired growth. Unexpectedly, none of the fungi showed expression specificity on mannose ( Figure 3C and Table S5), but some showed high sugar specificity on other monosaccharides, like A. niger and P. subrubescens on L-arabinose. For A. niger, it has been shown that mannobiose, rather than mannose, is the inducer of the mannanolytic enzyme system [59], which could explain why there is no mannose-specific induction observed.

Discussion
The genome content of CAZymes showed considerable differences among the tested species, especially among the taxonomically distant species. The CAZymes' repertoire and their expression under different carbon sources of each specific fungus imply the possible different PBD approaches. In general, the studied Eurotiomycetes possess more diverse sets of PBD genes that enable them to grow well on most tested polysaccharides, while T. reesei and the two Basidiomycetes grew poorly on most polysaccharides due to a relatively small repository or less actively expressed PBD genes. Particularly, T. reesei is the only tested Ascomycete that showed impaired growth on maltose and starch, which correlates with fewer starch-degrading enzymes identified in the genome of this species.
Among the three tested Eurotiomycota, only A. niger and P. subrubescens can grow well on both sucrose and inulin, while A. nidulans showed impaired growth on inulin. The divergence could be explained by the lower number of copies of GH32, and the low expression level of intracellular invertase SucB (Aspnid1|2301) [60].
The detailed structures and compositions of polysaccharides significantly affect their degradation by specific fungi, especially for species with less variety of enzyme families encoded in their genomes, which was demonstrated by the growth profile of T. reesei and P. chrysosporium on pectin extracted from two different plant sources.
Although crucial PBD CAZy genes were present in their genomes, the two tested wood-degrading Basidiomycetes showed relatively poor growth on most tested polysaccharides and showed clear different gene expression profiles of PBD genes than the tested Ascomycetes. The possible reason could be that these fungi mainly grow on wood in nature and have adapted their induction system to grow on complex polysaccharides as carbon sources instead of being optimized for individual polysaccharides as we tested here. In addition, other observations have suggested that inducers of CAZy genes in Basidiomycetes may be disaccharides, rather than the monosaccharides tested here (M. Mäkelä & R.P. de Vries, unpublished results).
Poor growth on cellulose was observed for all tested species, although cellulolytic genes were identified in their genomes. One possible explanation is that the necessary enzymes were not fully activated during fungal growth on cellulose. For instance, many cellulolytic genes in A. niger are controlled by the xylanolytic transcriptional activator XlnR, which was mainly activated in the presence of xylose [61,62]. In line with this, visible growths on the cellulose-derived oligosaccharide, cellobiose, were observed for all six species (Figure 3A), which require many fewer degrading enzymes to fully decompose than intact cellulose. In addition, commercial cellulose is not structurally identical to natural cellulose, due to the method of extraction.
We identified several CAZy families that are uniquely present in certain species or phyla, which could indicate a special physiological role and potential industrial applications. For example, the AA16 family only exists in the three Eurotiomycetes, while the AA9 family is broadly present in all six species. A recent study suggested these two fungal LPMO families are counterparts and target the same substrates [26]. Contradicting this, a recent study found that AA16 can only boost the cellulase activity of AA9 but cannot break down cellulose independently [63]. Clearly, more detailed studies in this area are warranted. In addition, GH44 enzymes were only identified in D. squalens [47], which indicates the scarcity and unique function of this enzyme. Additionally, we found CE15 only in the two Basidiomycetes and T. reesei, which was proposed to cleave ester linkages between lignin and glucuronoxylan [64]. GH141, whose active site pocket of the hydrolase is likely to be significantly more open than other fucosidases (e.g., GH95 and GH29), was only present in P. subrubescens, revealing that it has a more variable PBD mode [65]. Only A. nidulans contains genes encoding PL3 and PL9 pectate lyases, as well as PL11 rhamnogalacturonan lyases, which indicates that A. nidulans may have great potential for commercial applications in pectin degradation [66].
In addition, some genes encoding a specific family present phyla specificity. GH54 is only present in the four Ascomycetes and has been reported to remove decorations including L-arabinose and D-galactose in the presence of metal cofactors, which may suggest that it could be introduced in enzymatic cocktails to assist in the hydrolysis of xylan and xyloglucan. Similarly, GH62 was only found in the four Ascomycetes, and it is the sole GH family that contains only arabinofuranosidase (ABF) specificities [67]. However, high variability was also observed in the number of genes within the same phyla. For example, multiple gene copies of GH51-, GH54-, and GH62-encoding ABFs were identified in P. subrubescens, which implies a potential enhancement in the ability of P. subrubescens to release L-arabinose from plant biomass in a more specific manner [68]. Moreover, this is in line with our finding that CAZy genes involved in xylan, pectin, and xyloglucan in P. subrubescens show evident L-arabinose specificity ( Figure 3C).
Additionally, the expression profiles of PBD genes in six tested species showed dramatic differences during fungal growth on nine monosaccharides (Figures 2 and 3C). For the closely related Eurotiomycetes, although they have comparable content of pectinolytic genes, the varying expression patterns of these genes suggest diversity in the underlying complex regulatory networks. In A. niger, a comparable set of pectinolytic genes were induced on L-arabinose, L-rhamnose, and D-galacturonic acid, while the same set of genes in A. nidulans were largely induced by L-rhamnose alone, and in P. subrubescens they were mainly induced by L-arabinose and D-galacturonic acid (Table S5). Similarly, sugar-specific CAZymes involved in xylan or xyloglucan degradation in different fungi also displayed different inducing patterns on L-arabinose and D-xylose.

Conclusions
In this study, we revealed a high diversity of fungal plant polysaccharide degrading enzymes at both the genomic and transcriptomic levels. Considering the importance of the ecological role and increasing industrial potential of fungal plant biomass degradation, the findings revealed in this study will improve our understanding of fungal diversity and provide new insights into designing enzyme mixtures and metabolic engineering of fungi for related industrial applications.
Supplementary Materials: The following supporting information can be downloaded at: https://ww w.mdpi.com/article/10.3390/jof9080860/s1. Figure S1: CAZy families related to plant biomass degradation. Families in bold are involved in the degradation of multiple polysaccharides; Figure S2: Pie chart presenting the proportion of CAZy genes involved in plant biomass degradation in the genomes of the six fungi. GH: glycoside hydrolase, AA: auxiliary activities, CE: carbohydrate esterases, PL: polysaccharide lyases, FAE: feruloyl esterases; Figure S3: Heatmap visualization of expression profiles of CAZy genes with sugar-specific induction in the six fungi during their growth on nine monosaccharides; Table S1: CAZymes involved in plant polysaccharide degradation, and the number of putative CAZy-encoding genes in the six fungi related to the degradation of specific polysaccharides; Table S2: Orthologs of CAZymes detected among six studied species; Table S3: Summary of CAZy genes involved in plant polysaccharide degradation across the studied six fungi (including annotations from JGI); Table S4: Expression profiling of CAZyme-encoding genes involved in plant polysaccharide degradation across the studied six fungi; Table S5: Summary of CAZy genes with high sugar specificity in six fungi.