# Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species

Yi-Fei Sun, Annie Lebreton, Jia-Hui Xing, Yu-Xuan Fang, Jing Si, Emmanuelle Morin, Shingo Miyauchi, Elodie Drula, Steven Ahrendt, Kelly Cobaugh, et al.

► To cite this version:

Yi-Fei Sun, Annie Lebreton, Jia-Hui Xing, Yu-Xuan Fang, Jing Si, et al.. Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species. Journal of Fungi, 2022, 8 (3), pp.311. 10.3390/jof8030311 . hal-03627194

HAL Id: hal-03627194

https://amu.hal.science/hal-03627194

Submitted on 28 Jun 2023

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in *Ganoderma* Species

Yi-Fei Sun <fnref n="1" />,<fnref n="2" />, Annie Lebreton <fnref n="2" />,<fnref n="3" /> ID, Jia-Hui Xing <fnref n="1" />, Yu-Xuan Fang <fnref n="1" />, Jing Si <fnref n="1" />, Emmanuelle Morin <fnref n="2" /> ID,
Shingo Miyauchi <fnref n="2" />,<fnref n="4" />, Elodie Drula <fnref n="5" /> ID, Steven Ahrendt <fnref n="6" />, Kelly Cobaugh <fnref n="6" />, Anna Lipzen <fnref n="6" />, Maxim Koriabine <fnref n="6" />,
Robert Riley <fnref n="6" />, Annegret Kohler <fnref n="2" />, Kerrie Barry <fnref n="6" /> ID, Bernard Henrissat <fnref n="7" />,<fnref n="8" /> ID, Igor V. Grigoriev <fnref n="6" />,<fnref n="9" /> ID,
Francis M. Martin <fnref n="2" />,<fnref n="3" />,* and Bao-Kai Cui <fnref n="1" />,<fnref n="3" />,*

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

check for
updates

**Citation:** Sun, Y.-F.; Lebreton, A.; Xing, J.-H.; Fang, Y.-X.; Si, J.; Morin, E.; Miyauchi, S.; Drula, E.; Ahrendt, S.; Cobaugh, K.; et al. Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in *Ganoderma* Species. *J. Fungi* **2022**, *8*, 311. https://doi.org/10.3390/jof8030311

Academic Editors: Yang Xiao and
László Nagy

Received: 8 February 2022

Accepted: 16 March 2022

Published: 18 March 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The *Ganoderma* species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven *Ganoderma* genomes, including three new genomes, *G. australe*, *G. leucocontextum*, and *G. lingzhi*. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that *Ganoderma* species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the *Ganoderma* species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, *G. australe* has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. *G. australe* also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of *G. australe* secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the *Ganoderma* species.

**Keywords:** *Ganoderma*; genomics; secondary metabolism; secretome; terpenes; wood decay

<footnote n="1">¹ Institute of Microbiology, School of Ecology and Nature Conservation, Beijing Forestry University, Beijing 100083, China; yifeisun2016@163.com (Y.-F.S.); jhxing321@sina.com (J.-H.X.); yuxuanf2000@163.com (Y.-X.F.); jingsi1788@126.com (J.S.)</footnote>

<footnote n="2">² Université de Lorraine, INRAE, UMR Interactions Arbres/Microorganismes (IAM), Centre INRAE Grand Est-Nancy, 54280 Champenoux, France; annie.lebreton@inrae.fr (A.L.); emmanuelle.morin@inrae.fr (E.M.); smiyauchi@mpipz.mpg.de (S.M.); annegret.kohler@inrae.fr (A.K.)</footnote>

<footnote n="3">³ Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, Beijing Forestry University, Beijing 100083, China</footnote>

<footnote n="4">⁴ Max Planck Institute for Plant Breeding Research, Department of Plant Microbe Interactions, 50829 Cologne, Germany</footnote>

<footnote n="5">5 INRAE, Aix Marseille University, UMR1163 Biodiversité et Biotechnologie Fongiques, 13009 Marseille, France; elodie.drula@inrae.fr</footnote>

<footnote n="6">6 U.S. Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA; sahrendt0@lbl.gov (S.A.); klcobaugh@lbl.gov (K.C.); alipzen@lbl.gov (A.L.); mkoriabine@lbl.gov (M.K.); rwriley@lbl.gov (R.R.); kwbarry@lbl.gov (K.B.); ivgrigoriev@lbl.gov (I.V.G.)</footnote>

<footnote n="7">7 DTU Bioengineering, Technical University of Denmark, 2800 Kongens Lyngby, Denmark; bernard.henrissat@gmail.com</footnote>

<footnote n="8">8 Department of Biological Sciences, King Abdulaziz University, Jeddah 999088, Saudi Arabia</footnote>

<footnote n="9">⁹ Department of Microbial and Plant Biology, University of California Berkeley, Berkeley, CA 94720, USA</footnote>

<footnote>* Correspondence: francis.martin@inrae.fr (F.M.M.); cuibaokai@bjfu.edu.cn (B.-K.C.); Tel.: +33-383394080 (F.M.M.); +86-1062336309 (B.-K.C.)</footnote>
# 1. Introduction

*Ganoderma* species (Ganodermataceae, Polyporales, Basidiomycota) are both ecologically and economically relevant fungi found in forest ecosystems. As wood-decay fungi, *Ganoderma* species efficiently decompose the components of plant cell walls, including lignin, cellulose and hemicellulose [1]. *G. boninense* Pat. causes a basal stem rot (BSR) on oil palm trees [2,3], while *G. lucidum* and *G. australe* have been used for degradation of environmental pollutants [4,5]. As medicines, *Ganoderma* fruiting bodies have been used over the last 2000 years to produce drugs used for improving immunity, and in anti-aging and anti-cancer treatments in humans [6–8]. For example, spore powder or basidiocarp slices of *G. lingzhi* Sheng H. Wu, Y. Cao, and Y.C. Dai and *G. tsugae* Murrill are used in Asian traditional medicine to improve health. On the other hand, polysaccharides and glycans extracted from *G. sinense* J.D. Zhao, L.W. Hsu and X.Q. Zhang have been used as clinical drugs [9].

*Ganoderma* is the largest genus in Ganodermataceae including 461 taxa recorded in Index Fungorum (http://www.indexfungorum.org/, accessed on 16 January 2022) as of 15 September 2021. However, only four complete genomes of *Ganoderma* species have been reported at this date. The genome of the well-known medicinal polypore, *G. lucidum*, was published by Chen et al. [10]. This study focused on the analyses of genes encoding cytochrome P450s (CYPs), transporters and regulatory proteins which are related to secondary metabolism (SM) and wood degradation. Then, Binder et al. [11] used the genome of an unknown North American *Ganoderma* species to unravel phylogenetic relationships within the Polyporales. The sequencing and analyses of *G. sinense*, another well-known medicinal fungus, provided new highlights on the genome methylation patterns, small RNA transcriptome, SM and defense processes [12]. Utomo et al. [13] sequenced and analyzed a pathogenic strain of *G. boninense* isolated from an oil palm tree with severe symptoms of BSR disease. Additional genomes are required to explore the full diversity of *Ganoderma* gene sets involved secondary metabolism and wood white-rot decay and identify the genetic traits specific to the different species. Therefore, we sequenced and analyzed three additional genomes of *Ganoderma* species and compared them to other previously sequenced taxa, i.e., *G. australe*, *G. leucocontextum* and *G. lingzhi*.

*Ganoderma australe* is widely distributed on deciduous trees in central, eastern and southern China, and *G. lingzhi* is widely cultivated thanks to its high medicinal value, whereas *G. leucocontextum* was recently described from southwestern China and is now cultivated in Yunnan and Tibet for medicinal production. In addition to these three new genomes, we included in our analysis seven publicly available genomes of *Ganoderma* species. Finally, we compared the *Ganoderma* gene sets to those obtained on 51 Basidiomycota species, including white-rot and brown-rot wood decayers, and biotrophic fungi. Here, we provided new insights on the evolutionary relationships between *Ganoderma* species and other Polyporales, and we identified idiosyncrasies in the gene repertoire involved in SM and secreted enzymes.

# 2. Materials and Methods

## 2.1. Strain and Culture Conditions

Fruiting bodies of *Ganoderma leucocontextum* Dai 12418 was collected on *Quercus*, and the strain was cultured on Potato Dextrose Agar (Potato Extract 200 g, Agar 20 g, Dextrose 20 g, water 1 L); *G. australe* Cui 17254 and *G. lingzhi* Cui 9166 were collected on other angiosperm logs, and the strains were cultured on Malt Extract Agar (Malt Extract 20 g, Agar 18 g, $KH_2PO_4$ 3 g, Glucose 10 g, water 1 L). The strains were cultured for 7–14 days in the dark at 25 °C. The taxonomic affiliation of the three strains was confirmed by phylogenetic analyses based on ITS sequences. Isolates were deposited in the culture collection of the Institute of Microbiology, Beijing Forestry University and are available upon request.
## 2.2. DNA Extractions

For genomic DNA extraction, fresh 7-day-old vegetative agar mycelial cultures grown on cellophane sheets were harvested, snap frozen in liquid nitrogen and stored at -80 °C. High molecular weight genomic DNA of *Ganoderma australe* was extracted by using the sodium dodecylsulfate (SDS) method [14], whereas a modified cetyltrimethyl ammonium bromide (CTAB) method [15] was used for genomic DNA extraction of *G. leucocontextum* and *G. lingzhi*.

## 2.3. Genome Sequencing and Assembly

The genome of the strain Dai 12418 from *Ganoderma leucocontextum* was sequenced using the Pacific Biosciences platform (PacBio, Menlo Park, CA, USA); 1.5 µg of genomic DNA was sheared to 10 kb using Covaris g-Tube or Diagenode megaruptor tube. The sheared DNA was treated with DNA prep to remove single-stranded ends and DNA damage repair mix followed by end repair/A Tail and ligation of barcoded overhang adapters using SMRTbell Express Template Prep 2.0 Kit (PacBio, Menlo Park, CA, USA). The library was purified with AMPure PB beads and libraries with different barcodes were pooled at equimolar for up to four maximum fungal genomes with a total sum of genome sizes of 200 Mb. A BluePippin size selection (Sage Science, Beverly, MA, USA) was then performed on the pool to remove the shorter peaks physically. PacBio Sequencing primer was then annealed to the SMRTbell template library and sequencing polymerase was bound to them using Sequel II Binding kit 2.0. The prepared SMRTbell template libraries were then sequenced on a Pacific Biosystems' Sequel II sequencer using 8 M v1 SMRT cells and Version 2.0 sequencing chemistry with 1 × 1800 sequencing movie run times. Filtered subread data was processed to remove artifacts. Mitochondria was assembled separately with the CCS reads using an in-house tool (assemblemito.py), used to filter the CCS reads, and polished with gcpp—algorithm arrow version SMRTLink v8.0.0.80529 (https://www.pacb.com/support/software-downloads, accessed on 16 January 2022). The mitochondria-filtered CCS reads were then assembled with Flye version 2.7.1-b1590 (https://github.com/fenderglass/Flye, accessed on 16 January 2022) (—g 40 M—asm-coverage 50—pacbio-corr) and polished with gcpp—algorithm arrow version SMRTLINK v8.0.0.80529. Contigs less than 1000 bp were excluded.

The *Ganoderma australe* Cui 17254 and *G. lingzhi* Cui 9166 genomes were sequenced using PacBio and Illumina platform. Sequencing libraries for Illumina were generated using NEBNext® Ultra™ DNA Library Prep Kit (NEB, Ipswich, MA, USA). One µg of genomic DNA was fragmented by sonication to 350 bp (400 bp for *G. lingzhi* DNA) using Covaris g-Tube. The DNA fragments were end-polished, A-tailed, and ligated with the full-length adaptor for Illumina sequencing with further PCR amplification. The libraries were purified by AMPure XP system (Beckman Coulter, Brea, IN, USA) and were analyzed for size distribution by Agilent2100 Bioanalyzer and quantified using real-time PCR. The libraries for single-molecule real-time (SMRT) on PacBio platform was constructed with an insert size of 20 kb using the SMRT bell TM Template kit (version 1.0, Pacific Biosciences, Menlo Park, CA, USA). The DNA fragments were repaired DNA damage and ends and prepared blunt ligation reaction. The library was purified with 0.45X AMPure PB beads and size-selection using the BluePippin System (Sage Science, Beverly, MA, USA). The libraries were analyzed for size distribution by Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). The whole genome of *G. australe* was sequenced using PacBio Sequel and Illumina NovaSeq PE150 at Beijing Novogene Bioinformatics Technology Co., Ltd. (Beijing, China). The whole genome of *G. lingzhi* was sequenced using PacBio RS II platform and Illumina MiSeq platform at Shanghai Personalbio Technology Co., Ltd. (Shanghai, China). The low-quality reads were filtered (less than 500 bp) to obtain clean data. Preliminary assembly was conducted with SMRTLink v5.0.1 (https://www.pacb.com/support/software-downloads, accessed on 16 January 2022), and long reads (more than 6000 bp) were selected. By the variant Caller module of the SMRT Link software, the arrow algorithm was used to correct and count the variant sites in the preliminary assembly results.
## 2.4. Genome Annotation and Quality Check

The genome of *Ganoderma leucocontextum* was annotated using the JGI Annotation Pipeline with the support of their corresponding Trinity transcriptomes. Both assembly and annotations are available from JGI genome portal in MycoCosm (https://mycocosm.jgi.doe.gov/, accessed on 16 January 2022) [16,17]. The de novo annotation of *G. australe* and *G. lingzhi* were conducted by Augustus v3.03 [18], Genewise v2.4.1 [19] and Evidence-Modeler [20].

The quality of the predicted proteomes was evaluated by using the tool Benchmarking Universal Single-Copy Orthologs (BUSCO v.4.1.3) [21] with the Basidiomycota gene set downloaded from https://busco-data.ezlab.org/v4/data/lineages/basidiomycota_odb10.2020-09-10.tar.gz, accessed on 16 January 2022.

## 2.5. Synteny Analyses

The synteny analyses were conducted on the 10 largest scaffolds of each of the selected *Ganoderma* genomes. Pair-wise comparisons and identification of syntenic blocks were performed by using the R package DECIPHER [22] with default parameters. The synteny blocks between every two species were visualized with the R package Circlize [23]. Data management, integration, and visualization were as described in Hage et al. [24].

## 2.6. Phylogenetic Analyses

The phylogenetic analyses of 12 *Ganoderma* species were conducted by Maximum Likelihood (ML) with 23 ITS sequences. *Sanguinoderma* sp. Cui 17238 was used as the outgroup. The ML analyses were performed in RAxML-HPC v. 8.2.3 [25] involving 1000 ML searches under the GTRGAMMA model and 1000 rapid bootstrap replicates with the GTRCAT model to obtain the best tree and ML bootstrap. All trees were viewed in FigTree 1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/, accessed on 16 January 2022). The ML bootstrap ≥ 50% were presented on topologies from ML analyses.

## 2.7. Phylogenomic Analyses

A total of 58 Basidiomycota species was used for the phylogenomic analyses, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, a single endophyte and two pathogens. *Melampsora larici-populina* Kleb. and *Ustilago maydis* (DC.) Corda were used as outgroup species (Table 1).

<table><caption>**Table 1.** List of the 58 genomes used for phylogenomic and comparative analyses. The new genomes are shown in bold.</caption><thead><tr><th>Species Name</th><th>Strain</th><th>Reference</th><th>Lifestyle</th><th>Order</th><th>Clade</th></tr></thead><tbody><tr><td>Abortiporus biennis</td><td>CIRM-BRFM1778</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>residual polyporoid clade</td></tr><tr><td>Bjerkandera adusta</td><td>HHB-12826-SP</td><td>[11]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Ceriporiopsis subvermispora</td><td>B</td><td>[26]</td><td>White-rot</td><td>Polyporales</td><td>gelatoporia clade</td></tr><tr><td>Dichomitus squalens</td><td>LYAD-421 SS1</td><td>[27]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Earliella scabrosa</td><td>CIRM-BRFM 1817</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Epithele typhae</td><td>CBS 203.58</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Fomes fomentarius</td><td>CIRM-BRFM 1821</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Fomitiporia mediterranea</td><td>MF3/22 #7</td><td>[27]</td><td>White-rot</td><td>Hymenochaetales</td><td>-</td></tr><tr><td><strong>Ganoderma australe</strong></td><td>Cui 17254</td><td><strong>this study</strong></td><td><strong>White-rot</strong></td><td><strong>Polyporales</strong></td><td><strong>core polyporoid clade</strong></td></tr><tr><td>Ganoderma boninense</td><td>G3</td><td>[13]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td><strong>Ganoderma leucocontextum</strong></td><td>Dai 12418</td><td><strong>this study</strong></td><td><strong>White-rot</strong></td><td><strong>Polyporales</strong></td><td><strong>core polyporoid clade</strong></td></tr><tr><td><strong>Ganoderma lingzhi</strong></td><td>Cui 9166</td><td><strong>this study</strong></td><td><strong>White-rot</strong></td><td><strong>Polyporales</strong></td><td><strong>core polyporoid clade</strong></td></tr><tr><td>Ganoderma lucidum</td><td>G.260125-1</td><td>[10]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Ganoderma sinense</td><td>ZZ0214-1</td><td>[12]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Ganoderma sp.</td><td>10597 SS1</td><td>[11]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr></tbody></table>
<table><caption>**Table 1. Cont.**</caption><thead><tr><th>Species Name</th><th>Strain</th><th>Reference</th><th>Lifestyle</th><th>Order</th><th>Clade</th></tr></thead><tbody><tr><td>Heterobasidion annosum</td><td>TC32-1</td><td>[28]</td><td>White-rot</td><td>Russulales</td><td>-</td></tr><tr><td>Hexagonia nitida</td><td>CIRM-BRFM 1802</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Hydnopolyporus fimbriatus</td><td>CBS384.51</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Irpe x lacteus</td><td>CCBAS Fr. 238 617/93</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Leiotrametes lactinea</td><td>CIRM-BRFM 1664</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Leiotrametes menziesii</td><td>CIRM-BRFM 1781</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Lentinus tigrinus</td><td>ALCF2SS1-7</td><td>[29]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Obba rivulosa</td><td>3A-2</td><td>[30]</td><td>White-rot</td><td>Polyporales</td><td>gelatoporia clade</td></tr><tr><td>Panus rudis</td><td>PR-1116 ss-1</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>residual polyporoid clade</td></tr><tr><td>Phanerochaete carnosa</td><td>HHB-10118-Sp</td><td>[31]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Phanerochaete chrysosporium</td><td>RP-78</td><td>[32]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Phlebia brevispora</td><td>HHB-7030 SS6</td><td>[11]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Phlebia centrifuga</td><td>FBCC195</td><td>[33]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Phlebia radiata</td><td>FBCC0043-79</td><td>[34]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Phlebiopsis gigantea</td><td>5-6</td><td>[35]</td><td>White-rot</td><td>Polyporales</td><td>phlebioid clade</td></tr><tr><td>Polyporus arcuarius</td><td>HBB13444</td><td>[36]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Polyporus brumalis</td><td>BRFM 1820</td><td>[37]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Polyporus squamosus</td><td>CCBS 676</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Pycnoporus cinnabarinus</td><td>BRFM 137</td><td>[38]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Pycnoporus sanguineus</td><td>BRFM 1264</td><td>[39]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Rigidoporus microporus</td><td>ED310</td><td>[40]</td><td>White-rot</td><td>Polyporales</td><td>-</td></tr><tr><td>Stereum hirsutum</td><td>FP-91666 SS1</td><td>[27]</td><td>White-rot</td><td>Russulales</td><td>-</td></tr><tr><td>Trametes pubescens</td><td>FBCC735</td><td>[41]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Trametes versicolor</td><td>FP-101664 SS1</td><td>[27]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Trametopsis cervina</td><td>CIRM-BRFM 1824</td><td>[24]</td><td>White-rot</td><td>Polyporales</td><td>core polyporoid clade</td></tr><tr><td>Amylocystis lapponica</td><td>SKaAmylap13</td><td>[24]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Antrodia serialis</td><td>Sig1Antser10</td><td>[24]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Antrodia sinuosa</td><td>LB1</td><td>PRJNA196036</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Daedalea quercina</td><td>L15889 ss-12</td><td>[42]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Fibroporia radiculosa</td><td>TFFH 294</td><td>[43]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Fomitopsis betulina</td><td>CIRM-BRFM 1772</td><td>[24]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Fomitopsis pinicola</td><td>FP-58527 SS1</td><td>[27]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Laetiporus sulphureus</td><td>93-53</td><td>[42]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Postia placenta</td><td>MAD-698-R-SB12</td><td>[44]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Rhodofomes roseus</td><td>CIRM-BRFM 1785</td><td>[24]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Wolfiporia cocos</td><td>MD-104 SS10</td><td>[27]</td><td>Brown-rot</td><td>Polyporales</td><td>antrodia clade</td></tr><tr><td>Amanita muscaria</td><td>Koide</td><td>[45]</td><td>Ectomycorrhizal</td><td>Agaricales</td><td>-</td></tr><tr><td>Laccaria bicolor</td><td>S238N-H82</td><td>[46]</td><td>Ectomycorrhizal</td><td>Agaricales</td><td>-</td></tr><tr><td>Pisolithus microcarpus</td><td>441</td><td>[45]</td><td>Ectomycorrhizal</td><td>Boletales</td><td>-</td></tr><tr><td>Rhizopogon vinicolor</td><td>AM-OR11-026</td><td>[47]</td><td>Ectomycorrhizal</td><td>Boletales</td><td>-</td></tr><tr><td>Piriformospora indica</td><td>DSM 11827</td><td>[48]</td><td>Endophyte</td><td>Sebacinales</td><td>-</td></tr><tr><td>Melampsora larici-populina</td><td>98AG31</td><td>[49]</td><td>Pathogen</td><td>Pucciniales</td><td>-</td></tr><tr><td>Ustilago maydis</td><td>521</td><td>[50]</td><td>Pathogen</td><td>Ustilaginales</td><td>-</td></tr></tbody></table>

The orthologous protein clusters of the 58 proteomes were identified with OrthoFinder v2.4.0 [51]. The sequences of each of 143 single copy, conserved orthologs were aligned using MAFFT v7.471 [52]. Poorly aligned regions were removed with Trimal v1.4.1 [53]. Based on these alignments, a maximum likelihood phylogenomic tree was constructed by RAxML-NG v0.9.0 [54] using partitions corresponding to an orthologous group, and their associated best-fit model for each partition of the concatenate protein alignments were estimated by Modeltest-NG v0.1.6 (Table S1) [55].

## 2.8. Annotation of Transposable Elements

Transposable elements (TEs) were identified as described in Payen et al. [56]. Briefly, de novo repeat sequences were predicted in unmasked genome assemblies of 58 genomes,
using RepeatScout 1.0.6 [57]. Sequences ≥ 100 bp and ≥10 occurrences were filtered out. The selected sequences were annotated by searching homologous sequences against the fungal references in REPBASE v.22.08 (http://www.girinst.org/server/RepBase/index.php, accessed on 16 January 2022) using tBLASTx [58]. The coverage of TEs in the genomes, including unknown TEs, was estimated by REPEATMASKER open 4.1.1 (http://www.repeatmasker.org, accessed on 16 January 2022).

## 2.9. Protein Functional Annotation

The gene clusters related to SM biosynthetic pathways were predicted by antiSMASH 4.2.0 [59] and visualized along with the species phylogenetic relationship on iTOLv5 (https://itol.embl.de/, accessed on 16 January 2022) [60]. PFAM domain searches were performed with HMMER [61]. The most abundant PFAM protein domains (abundance > 100) were visualized along with the species phylogenetic relationship on iTOLv5. Secreted proteins were predicted as described in Pellegrin et al. [62], and proteins with a size < 300 amino acids were identified as small secreted proteins (SSP). Carbohydrate-active enzymes (CAZymes) were manually curated by the CAZy team (http://www.cazy.org, accessed on 16 January 2022). The annotation of secreted proteases and lipases was performed by BLASTP search (E-value = 10⁻⁵) against MEROPS (http://merops.sanger.ac.uk/, accessed on 16 January 2022) and Lipase Engineering Database (http://www.led.uni-stuttgart.de/, accessed on 16 January 2022).

# 3. Results

## 3.1. Main Genome Features

Three *Ganoderma* species were newly sequenced in this study: *Ganoderma australe* strain Cui 17254, *G. leucocontextum* strain Dai 12418 and *G. lingzhi* strain Cui 9166. The size of the *G. australe* assembly was 84.27 Mb and 20,460 protein-coding genes were predicted (Table 2). This is the largest *Ganoderma* genome sequenced to date. The size of the genome assemblies for *G. leucocontextum* and *G. lingzhi* was lower at 60.34 and 60.56 Mb, and 15,007 and 16,592 protein-coding genes were predicted on the assemblies, respectively. Between 73.1% to 99.8% of a benchmark set of conserved fungal genes (BUSCO) were found in genome assemblies, indicating that assembled genomes captured most of the coding gene space, although the gene annotation for *G. lingzhi* appeared to be more fragmented (Table 2).

<table><caption>**Table 2.** Genome features of the three newly sequenced *Ganoderma* genomes.</caption><thead><tr><th>Genome Feature</th><th>G. australe</th><th>G. leucocontextum</th><th>G. lingzhi</th></tr></thead><tbody><tr><td>Genome size (Mb)</td><td>84.27</td><td>60.34</td><td>60.56</td></tr><tr><td>Number of scaffolds</td><td>93</td><td>843</td><td>342</td></tr><tr><td>Scaffold L50 (bp)</td><td>1,745,385</td><td>205,166</td><td>402,014</td></tr><tr><td>Scaffold N50</td><td>17</td><td>66</td><td>37</td></tr><tr><td>Longest scaffold (bp)</td><td>4,455,856</td><td>1,715,371</td><td>2,154,085</td></tr><tr><td>Shortest scaffold (bp)</td><td>33,016</td><td>1008</td><td>980</td></tr><tr><td>GC content (%)</td><td>55.48</td><td>55.95</td><td>55.88</td></tr><tr><td>Protein-coding genes</td><td>20,460</td><td>16,952</td><td>15,007</td></tr><tr><td>Average gene length (bp)</td><td>1582</td><td>1846</td><td>1605</td></tr><tr><td>Complete BUSCOs (%)</td><td>84.6</td><td>99.8</td><td>73.1</td></tr></tbody></table>

We compared the three newly sequenced *Ganoderma* genomes to four published *Ganoderma* genomes and 51 other Basidiomycota genomes. The genomes size of the 58 investigated species ranged from 19.66 to 109.88 Mb (Figure 1), with 6785 to 26,226 predicted genes (Table S2). Except for *Melampsora larici-populina*, *Ganoderma* species displayed the largest genomes, i.e., *G. australe* and *G. boninense* were two-to three-fold larger than other wood decayers. No significant differences (*p* > 0.05) in the average genome size were found between white-rot fungi, brown-rot fungi and ectomycorrhizal fungi. The number of
duplicated BUSCO genes was higher (23.4%) in *G. australe*. The synteny analyses showed no evidence for whole-genome or segmental duplications (Figure S1), it may indicate a polymorphic dikaryon (Table S2).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Comparison of the genome size among the 58 selected species. Bars represent size of 58 genomes and boxplots show proportion of genome size in different lifestyles. The new genomes are shown in red.</figcaption></figure>

## 3.2. Macrosynteny between *Ganoderma* Genomes

The top 10 scaffolds of each *Ganoderma* genomes, covering 13% to 75% of the whole assemblies, were selected to perform a macrosynteny analysis (Figure S2, Table S3A). We observed the highest percentage of syntenic segments (72%) between *G. lingzhi* and *G. lucidum*, while *G. boninense* and *G. leucocontextum* showed a lower rate when compared to other species (Table S3B), reflecting a higher genome divergence.

## 3.3. Phylogenetic Analyses of *Ganoderma* Strains

The phylogenetic relationship of 12 *Ganoderma* species was conducted based on 23 ITS sequences and one *Sanguinoderma* sp. as outgroup (Figure 2). Two *G. lucidum* strains were not clustered together indicating that the identification of *G. lucidum* G.260125-1 should be considered as *G. lingzhi* actually. *Ganoderma* sp. 10597 SS1 clustered with *G. sessile* suggesting that this strain may pertain to this sessile species in this phylogeny.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** ML analyses of 12 *Ganoderma* strains based on ITS sequences. ML bootstrap values higher than 50% are shown. Stains with sequenced genomes are shown in bold.</figcaption></figure>

3.4. *Phylogenomic Analyses of the Ganoderma Species and Related Polyporales*

Our phylogenetic analyses, based on 143 single-copy conserved protein sequences, was in agreement with previous phylogenetic analyses using single or multi-locus approaches [24,63]. We also identified five major clades in Polyporales: the core polyporoid clade, the antrodia clade, the gelatoporia clade, the phlebioid clade, and the residual polyporoid clade. *Ganoderma* species clustered in the core polyporoid clade with a high bootstrap value (Figure 3, 100% ML bootstrap). The phylogenomic status of *Ganoderma* species was consistent with the macro-synteny conservation results.

The taxonomic status of *Rigidoporus microporus* has been changed to Hymenochaetales [64], and here, it was confirmed again by phylogenomic analyses. Among the sampled species, the brown-rot fungi and the ectomycorrhizal fungi formed monophyletic groups. The present analyses also confirmed that the white-rot lifestyle is evolutionary polyphyletic.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Phylogenetic relationship within the *Ganoderma* clade and 58 selected Basidiomycota species. The 'best tree' resulting from maximum likelihood (ML) analyses based on 143 single-copy genes and ML bootstrap values are shown. Colors are coded for the five lifestyles.</figcaption></figure>

### 3.5. Transposable Element Identification

The dominant TEs in *Ganoderma* species belonged to the *Gypsy* and *Copia* families of long terminal repeats (LTR) retrotransposons. Besides, the proportion of unknown TEs in *Ganoderma* species was large, especially in *G. australe* (10.5% of the total assembly) and *G. lingzhi* (11.1% of the total assembly), indicating they likely played a key role in genome rearrangements. *G. leucocontextum* displayed the largest TE coverage (18.32% of the total assembly), and the most diverse TE distribution, including simple repeats (*Repetitive*), *IS3EU* and *Dada* of DNA transposon, *RTEX* and *L1* of non-LTR retrotransposon which are unique repeat elements in *Ganoderma* species (Figure 4A). The number of *IS3EU* and *Dada* sequences in *G. leucocontextum* was significantly larger comparing with other species (Figure 4B). Although occurring at moderate copy numbers, the *G. boninense* genome contained unique TE families, such as the *DIRS* LTR retrotransposons and several DNA repeated elements. *G. australe* also contained more copy of *Tad1* non-LTR retrotransposon and *Helitron* DNA transposons by comparing to other *Ganoderma* species.

TE coverage in the 58 analyzed genomes ranged from 0.37% (*Phlebiopsis gigantea*) to 41.69% (*Fomitiporia mediterranea*) (Table S4). Comparing to other lifestyle fungi, white-rot
fungi have low TE coverage except *F. mediterranea* (*p* < 0.01), among which *Ganoderma* species have relatively higher TE coverage (*p* < 0.05). Ectomycorrhizal fungi showed the larger repeat element coverage, in which unknown TE elements were the most abundant (>12%).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** (A) Distribution and coverage of transposable elements (TEs) identified in the 58 selected genomes. The bubble size is proportional to the coverage of each of transposable elements (shown inside the bubbles). The bars on top show the total coverage per genome. (B) The copy number of transposable elements (TEs) identified in the 58 genomes. The bubble size is proportional to TE copy number (shown inside the bubbles). The bars on top show the total copy number per genome. Color codes for the five fungal lifestyles are shown at the bottom of the figure.</figcaption></figure>
### 3.6. Biosynthetic Gene Clusters

We identified a total of 16 types of SM biosynthetic clusters (Figure 5, Table S5). Except for the *Rhizopogon vinicolor*, the pathogens, ectomycorrhizal fungi and endophytes generally had a lower number of genes involved in SM than saprotrophic species. This difference was mainly driven by a lower content in t1pks, terpenes related genes and miscellaneous genes tagged as "putative" and "other" by the antiSMASH software. The antrodia clade composed of brown-rot species was segregated from other Polyporales clades by an enrichment in t1pks, indole, and fatty acid associated genes, and a depletion in NRPS. Among white-rot fungi in Polyporales, the phlebioid clade displayed a higher content in t3pks and a lower content in terpenes associated genes. The core polyporoid clade contained the species with diverse sets of biosynthetic genes clusters. Gene clusters involved in terpene synthesis and unknown metabolites (i.e., putative SM gene clusters) were enriched in *Ganoderma* species, especially in *G. australe* which contained 31 terpene associated gene clusters and 80 putative biosynthetic gene clusters. Noteworthy, *G. lingzhi* and *G. lucidum* as the main medically relevant species showed the lowest number of terpenes associated gene clusters, whereas *G. australe* and the pathogenic *G. boninense* had the highest content. *G. lucidum*, *G. leucocontextum*, *Trametes versicolor*, and *T. pubescens* genomes encoded a cluster related to lantipeptide production, and the first three were identified as associated with a t1pks.

### 3.7. Pfam Protein Domains Found in the Genomes

More than 5000 Pfam protein domains were identified in the 58 selected genomes. A total of 32 Pfam categories were sorted according to their gene copy number (>100) (Figure 6, Table S6). *Ganoderma australe* showed the higher number of Pfam protein domains, i.e., five-fold larger than *Ustilago maydis* (with 864 protein domains). In *G. australe*, 13 protein domains were prominent, including cytochrome P450s (PF00067: p450) involved in SM, heterokaryon incompatibility protein (PF06985: HET) often related to vegetative incompatibility (VI) and membrane transporters of the major facilitator superfamily (PF07690: MFS_1). Besides, other protein domains putatively playing a role in epigenetic regulation (PF00078: RVT_1, PF00385: Chromo, PF00665: rve) and protein-protein interactions (PF12937: F-box-like) were also enriched in *G. australe*. Comparing to other *Ganoderma* species, *G. boninense* and *G. leucocontextum* encoded additional protein domains (e.g., PF17667: Pkinase_fungal, PF18758: KDZ, PF18759: Plavaka, PF18803: CxC2, PF20149: DUF6532, PF20151: DUF6533, and PF20152: DUF6534), which were only found in parasitic and symbiotic fungi. Protein kinases domains (PF00069: Pkinase, PF07714: Pkinase_Tyr) and protein-protein interactions domains (PF00400: WD40, PF12894: A0PC4_WD40) were especially enriched in ectomycorrhizal fungi.
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Presence and abundance of the gene clusters involved in secondary metabolite biosynthesis along with the species phylogenetic relationship between the 58 selected fungal species. The heatmap depicts the number of the gene clusters according to the color scale from blue to red. Color codes for the five fungal lifestyles are shown at the top left of the figure.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Distribution and abundance of the top 100 Pfam protein domains in the selected 58 species. The heatmap depicts the protein domain copy number according to the color scale from purple to green. Color codes for the five fungal lifestyles are shown at the top left of the figure.</figcaption></figure>

### 3.8. *The Predicted Secretome*

*Ganoderma* species have the largest secretome among the 58 selected fungal species (Figure 7, Table S7). Small secreted proteins (SSPs) represented 41% to 81% of the total secreted proteins and 4% to 28% SSPs were annotated as CAZymes, lipases, or proteases. Most known SSPs were annotated as CAZymes, especially enriched in *G. australe* and *G. lingzhi* (Figure 8, Table S8).

As secreted CAZymes plays a key role in wood degradation, the evolution of their content through Polyporales was investigated. The different clades in Polyporales have a distinctive CAZyme content and the *Ganoderma* species also contained a distinctive CAZyme repertoire within the core polyporoid clade. The antrodia clade composed of brown-rot fungi have a two-fold lower content of secreted CAZymes compared to the white-rot species in Polyporales (141 ± 24 vs. 237 ± 71 respectively, Table S9). Among the 148 CAZymes sub-families (grouped in 101 families) with secreted genes representatives, 48 (grouped in 38 families) exhibited a lower number of genes in brown-rots compared to white-rots in Polyporales (BM test, FDR padj < 0.01). It included cellulolytic enzymes (GH6, GH7), LPMOs (AA9), ligninolytic PODs (AA1_1), heme-associated PODs (AA2), the carbohydrate-binding module CBM1, CAZymes involved in bacterial cell wall degradation (GH25, GH79), fungal cell wall degradation (GH20, GH76, GH92, GH128, GH135, GH152),
and other CAZymes involved in plant cell wall degradation. CAZymes was globally enriched in the core polyporoid clade compared to the phlebioid clade. This enrichment was associated to the expansion of AA1, AA14, CE1, EXPN, GH16, GH17, GH18, GH25, and GH30 families. Of note, the AA7 family was absent from the phlebioid clade and AA12 was absent from the core polyporoid clade. Among the core polyporoid clade, the *Ganoderma* species have enriched CAZymes, and 11 families were expanded, including GH18, GH16, AA1, GH43, CE16, GH3, GH128, GH47, GH115, GH25, and GH1 families. Whereas only PL4 was enriched (and only found) in the other species of the core polyporoid clade. Compared to other species, *G. australe* had the largest set of secreted CAZymes, in which GH18 and EXPN families were highly enriched in this species. It also contained the largest repertoire of PCWDEs, MCWDEs, and enzymes acting on pectin, peptidoglycans, and chitin (Figure 9, Table S10).

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Distribution of predicted secreted proteins (i.e., secretome) in *Ganoderma* species and other selected fungi. Bars represent the gene copy number for CAZymes, lipases, proteases, and SSPs. Color codes for the five fungal lifestyles are shown at the top left of the figure.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Distribution of small secreted proteins (SSPs) in *Ganoderma* species and other selected fungi. Bars represent the number of annotated (CAZymes, lipases, proteases) and other SSPs. All species are annotated by five lifestyles.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure 9.** Distribution of secreted CAZymes involved in plant and microbial cell wall degradation in *Ganoderma* species and other selected fungal genomes. The bubble size is proportional to the number of secreted CAZymes grouped for 11 categories. Colors are coded by five lifestyles. The bar plots show the count of genes involved in PCWDE and MCWDE (**left**), and the ratio of PCWDE to MCWDE (**right**).</figcaption></figure>
A total of 3880 genes encoding for secreted proteases (in 73 MEROPS families) were identified in 58 fungal genomes (Table S11). Eight protease families (A01A, G01, M28E, M36, S08A, S10, S28, and S53) have more than 100 members, and the largest represented protease subfamily was A01A with 1332 proteins. No significant differences were found between saprotrophic and symbiotic species (*p* > 0.05), but five protease families (C40, M24X, M57, M77, S08B) were only found in white-rot fungi, although with low gene copy numbers. Among these families, the M57 protease family including four proteins, was found in *Ganoderma australe*, *G. leucocontextum*, *G. lingzhi*, and *G. lucidum*. In addition, *Ganoderma* species have more G01, M28E, and M35 proteases.

Comparing to other lifestyles, white-rot fungi have the higher set of secreted lipases (*p* < 0.001) and among the 58 analyzed fungal secretomes, *Ganoderma australe* contained the highest number of secreted lipases, the most abundant being GGGX lipases (Table S12). Most genomes displayed less than seven GX lipases, except for *G. boninense* and *G. leucocontextum* with ten and nine GX lipases respectively. They are prominent enzymes catalyzing a wide range of reactions on various cellular substrates [65,66].

**4. Discussion**

In this study, we provided three newly sequenced genomes of ecologically and economically relevant *Ganoderma* species. A phylogeny, based on 23 rDNA ITS sequences from 12 *Ganoderma* species, allowed us to determine the phylogenetic status of the newly sequenced species in *Ganoderma*. This phylogeny concurred with previous studies which were carried out by using a smaller set of *Ganoderma* species [11,63,67], except for *G. lucidum*. Indeed, this species displayed inconsistency between the two available strains, suggesting a misidentification of *G. lucidum* G.260125-1 strain which was purchased from a company. According to our study, this strain may pertain to the *G. lingzhi* which is widely cultivated for its medicinal usage in China [7]. This phylogeny also indicated that *Ganoderma* sp. 10597-SS1 probably pertained to *G. sessile*. The phylogenomic analyses based on 58 genomes, including seven *Ganoderma* species, supported the ITS phylogenetic result and confirmed that *Ganoderma* belong to the core polyporoid clade of Polyporales as defined by Justo et al. [63].

*Ganoderma* species have been used for centuries in traditional medicine thanks to their well-known arsenal of antimicrobial, anti-aging, antioxidant, anti-inflammatory, and immunomodulating compounds (e.g., polysaccharides, triterpenoids, and peptides) [68,69]. The drastic reduction of sequencing cost in the last decades allowed genome-wide mining of medicinal compound. As a result, the number of new genomes rapidly increased over the year, such as genomes of *G. tsugae* CCMJ4178 [67] and *G. leucocontextum* DH-8 [70]. The strain DH-8 of *G. leucocontextum* was newly sequenced and its genome was 50.05 Mb with 58 scaffolds [70], the genome of strain Dai 12418 of *G. leucocontextum* sequenced in our study was 60.34 Mb with 843 scaffolds. In our study, 21 terpene gene clusters were predicted in strain Dai 12418 against the 10 terpene gene clusters predicted in strain DH-8. The two-fold enrichment in strain Dai 12418 compared to strain DH-8 was unexpected, despite taking into differences in annotation tools and the potential bias of assembly quality. It may suggest a substantial intraspecific genome polymorphism. Direct comparisons of other functional gene categories were hampered by major differences in annotation methodologies. For example, Liu et al. [70] predicted 614 CAZymes genes using HMM, while we only found 291 CAZymes genes, using the in-house pipeline from the CAZymes database followed by expert manual curation.

One of the main targets for medically-relevant products are the secondary metabolites associated genes, especially the terpenes in *Ganoderma* species. Compared to other species, the terpene genes of *Ganoderma* species were indeed expanded, however, substantial differences were observed among this genus. *G. lingzhi* and *G. lucidum* as the main medically species have the lowest number of terpenes associated clusters whereas the widespread *G. australe* and the pathogenic *G. boninense* have the highest content. Terpene related gene expansion in *G. australe* and *G. boninense* might explain their ecological ability to develop
on a broader set of substrates. On the other hand, a shift of the terpenome composition in *G. lingzhi*, *G. tsugae*, and *G. leucocontextum* might have led to the production of terpenes with beneficial properties for human health. This study also allowed us to identify a polyketide synthase complex involved in the synthesis of the antibacterial lantipeptide [71] in *G. lucidum*, *G. leucocontextum*, *Trametes versicolor*, and *T. pubescens*, likely explaining the use of formulations based on these fungi as antibiotics [10,72–74]. Besides, the large content of gene clusters encoding unknown biochemical function(s) were found in Polyporales species, especially *G. australe*, suggested an outstanding ability to synthesize a large set of secondary metabolites of yet unknown function. This confirmed that the economically and medically relevant secondary metabolites of these fungi represent an untapped resource.

As Polyporales is an important group of wood decayers, we investigated the decay ability of these species by analyzing their secretome. The current evolution was posited that early Agaricomycetes were saprotrophic and different lifestyles were derived from it [75,76]. It also supported that the white-rot lifestyle is linked to the ability to degrade lignin. The acquisition of lignin degradation in the Agaricomycetes was estimated in the late Carboniferous and further evolved with the evolution of lignin complexity in Polyporales and Agaricales taxa [77]. The monophyletic antrodia clade, only composed of brow-rot species, was derived from white-rot lineages in current phylogenomic analyses. According to Baldrian and Valášková [75], brown-rot fungi, despite having several independent origins, were associated with the loss of ligninolytic PODs, heme dye-decolorizing PODs, heme-thiolate POD/peroxygenases (HTPs), cellulolytic enzymes (GH6, GH7, LPMO), and the carbohydrate-binding module CBM1 reduced their ability to degrade lignin into its partial modification by releasing Fenton-generated hydroxyl radicals in the colonized material. This statement was confirmed in the antrodia clade, in which those losses were also accompanied by losses in other PCW-degrading enzymes, and FCW- and BCW-degrading enzymes. The enhanced ability of white-rot fungi to decompose substrates was also marked by a secreted lipase enrichment. White-rot specific proteases (C40, M24X, M57, M77, S08B) were identified, and they were related to genetic regulation, bacterial cell-wall modification, nutrient transformation, and other cellular physiological functions.

The various white-rot Polyporales clades also displayed divergent wood degrading abilities. The potential ligninolytic ability and more generally PCW degradation ability of the core polyporoid clade was the highest among Polyporales due to the expansion of ligninolytic PODs (AA1) and hemicellulose degrading enzyme (AA14, CE1, GH30). The difference among white-rot clades also showed in the enzyme arsenal to compete with other microbes for substrates and copy with biotic threats. The specific expansion in BCW-degrading enzyme (GH25) and FCW-degrading enzyme (GH16, GH17, GH18) was also observed in the core polyporoid clade. *Ganoderma* species deepened their lignocellulose degradation ability with further expansion of the ligninolytic PODs (AA1), hemicellulose (GH43, CE16, GH115), BCW- (GH25), and FCW- (GH16, GH18, GH128, GH47) degrading enzymes. Aminopeptidase Ap1 (M28E), deuterolysin (M35) and scytalidoglutamic peptidase (G01) were also enriched in *Ganoderma* species, but further functional analyses are needed to clarify their role(s). Noteworthy, *G. australe* showed an even wider degrading gene repertoire including higher PCW-degrading CAZyme, lipase and protease content. Similarly, numerous cytochrome P450 monooxygenases (PF00067: p450) were found in *G. australe*. These enzymes are known for their role in lignin and xenobiotic degradation [78]. This pattern supported the known ability of *G. australe* to more efficiently decompose lignocellulose than other *Ganoderma* species. An increased CAZyme repertoire in BCW- and FCW-degrading enzymes was also found in *G. australe*, which could perform the higher competitiveness than other microorganism for the substrate during colonization on the hosts. Associated with its enriched number of biosynthetic genes and clusters, such as lytic transglycolase (PF03330) and cerato-platanins (PF07249), this repertoire could be used to combat with miscellaneous biotic threats.

The abundance of HET genes related to heterokaryon incompatibility, MFS transporters, and epigenetic regulation in *Ganoderma australe* indicated that these mechanisms
play a key role in the stable reproduction and evolution of this wood decayer [79–81]. With 168 protein domains detected among seven *Ganoderma* species (PF00249, PF08914, PF11831, PF12776, PF13837, PF13921, PF15963), MYB transcription factors were also abundant in *Ganoderma* species. These transcription factors are the largest transcription factor families in eukaryotic organisms and play key role in variable development and physiological activities [82,83]. Wang et al. identified 75 MYB transcription factors in five *Ganoderma* species after manually curation, and of the gene copy number found in each species was lower than the detected results in our study. This is likely due to the differences between strains and technical methods used here. Further exploration of MYB genes can help to clarify its potential function during the growth and development of *Ganoderma* fungi.

The present comparative analysis of the publicly available *Ganoderma* genomes re-
vealed a series of genetic features specific to this lineage of wood decayers. This study
provided foundational information to characterize further ecological traits of this important
group of decomposers. In addition, this information will be used to characterize the regula-
tion of genes involved in SM biosynthesis pathways at the transcriptomic level, including
antimicrobial compounds and medicinal drugs.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof8030311/s1, Figure S1: Synteny analyses of *Ganoderma australe* itself. The scatter plot shows correlation between all contigs, Figure S2: Macrosynteny between *Ganoderma* species. The percentage of syntenic hits occurring in the same order on the compared block of sequences between two genomes are shown inside each Circos plot, Table S1: Model list for each orthologous group of concatenate protein alignments, Table S2: Main features of the 58 genomes used in this study, Table S3: (A) Size and genome coverage of the 10 largest scaffolds in *Ganoderma* genomes, (B) synteny summary statistic of the 10 largest scaffolds in *Ganoderma* species, Table S4: Summary statistics of components and coverage for transposable elements (TEs) in the 58 genomes, Table S5: Number of the gene clusters involved in secondary metabolites biosynthesis identified in the 58 genomes, Table S6: Number of the filtered protein domains (abundance > 100) in Pfam database among the 58 genomes, Table S7: Count of each category (CAZymes, lipases, SSP, proteases) in 58 fungal secretomes, Table S8: Number of small secreted proteins in 58 secretomes, Table S9: Number of secreted CAZymes in 58 secretomes, Table S10: Composition and ratios of secreted CAZymes for plant and microbial cell wall degradation in 58 genomes, Table S11: Number of secreted proteases in 58 secretomes, Table S12: Number of secreted lipases in 58 secretomes.

**Author Contributions:** Supervision: B.-K.C. and F.M.M.; data curation: Y.-F.S., B.-K.C. and F.M.M.; formal analysis: Y.-F.S. and A.L. (Annie Lebreton); methodology: F.M.M., A.L. (Annie Lebreton), S.M., E.M. and E.D.; resources: Y.-F.S., J.-H.X., Y.-X.F., J.S., S.A., K.C., A.L. (Anna Lipzen), M.K., R.R., A.K., K.B., B.H. and I.V.G.; software: A.L. (Annie Lebreton) and E.M.; visualization: Y.-F.S. and A.L. (Annie Lebreton); writing—original draft: Y.-F.S.; writing—review and editing: A.L. (Annie Lebreton), B.-K.C. and F.M.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the National Natural Science Foundation of China (Nos. U2003211, 31870008), the Scientific and Technological Tackling Plan for the Key Fields of Xinjiang Production and Construction Corps (No. 2021AB004), Beijing Forestry University Outstanding Young Talent Cultivation Project (No. 2019JQ03016), the Laboratory of Excellence Advanced Research on the Biology of Tree and Forest Ecosystems (ARBRE; ANR-11-LABX 0002 01, project POLOMICS) and the China Scholarship Council. The work conducted by the U.S. Department of Energy Joint Genome Institute, a DOE Office of Science User Facility, was supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data generated or analyzed during this study are included in this published article and supplementary materials. Genomic data used in this study (Table 1) are available from NCBI (https://www.ncbi.nlm.nih.gov/genome, accessed on 16 January 2022) and the JGI Genome Portal (http://genome.jgi.doe.gov, accessed on 16 January 2022). The raw sequencing data of *Ganoderma australe*, *G. lingzhi* and *G. leucocontextum* is deposited on NCBI linked to BioProject:
PRJNA775667 and PRJNA791677, Biosample: SAMN22612695, SAMN22612696 and SAMN24337934, and the accession numbers are SRR16605592, SRR16605591 and JAKETP000000000.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Rodríguez-Couto, S. Industrial an environmental application of white-rot fungi. *Mycosphere* **2017**, *8*, 456-466. [CrossRef]

2. Pilotti, C.A. Stem rots of oil palm caused by *Ganoderma boninense*: Pathogen biology and epidemiology. *Mycopathologia* **2005**, *159*, 129-137. [CrossRef] [PubMed]

3. Seman, I.B. R&D on Biology, Detection and Management of *Ganoderma* Disease in Oil Palm. Ph.D. Thesis, Universiti Putra Malaysia, Selangor, Malaysia, 2018.

4. Si, J.; Wu, Y.; Ma, H.F.; Cao, Y.J.; Sun, Y.F.; Cui, B.K. Selection of a pH- and temperature-stable laccase from *Ganoderma australe* and its application for bioremediation of textile dyes. *J. Environ. Manag.* **2021**, *299*, 113619. [CrossRef] [PubMed]

5. Wang, H.; Deng, W.; Shen, M.; Yan, G.; Zhao, W.; Yang, Y. A laccase Gl-LAC-4 purified from white-rot fungus *Ganoderma lucidum* had a strong ability to degrade and detoxify the alkylphenol pollutants 4-n-octylphenol and 2-phenylphenol. *J. Hazard. Mater.* **2021**, *408*, 124775. [CrossRef] [PubMed]

6. Boh, B.; Berovic, M.; Zhang, J.S.; Lin, Z.B. *Ganoderma lucidum* and its pharmaceutically active compounds. *Biotechnol. Annu. Rev.* **2007**, *13*, 265-301. [CrossRef] [PubMed]

7. Cao, Y.; Wu, S.H.; Dai, Y.C. Species clarification of the prize medicinal *Ganoderma* mushroom 'Lingzhi'. *Fungal Divers.* **2012**, *56*, 49-62. [CrossRef]

8. Zhang, Y.R.; Jiang, Y.F.; Zhang, M.; Zhang, L.J. *Ganoderma sinense* polysaccharide: An adjunctive drug used for cancer treatment. *Prog. Mol. Biol. Transl. Sci.* **2019**, *163*, 165-177. [CrossRef] [PubMed]

9. Chan, K.M.; Yue, G.G.L.; Li, P.; Wong, E.C.W.; Lee, J.K.M.; Kennelly, E.J.; Lau, C.B.S. Screening and analysis of potential anti-tumor components from the stipe of *Ganoderma sinense* using high-performance liquid chromatography/time-of-flight mass spectrometry with multivariate statistical tool. *J. Chromatogr. A* **2017**, *1487*, 162-167. [CrossRef]

10. Chen, S.L.; Xu, J.; Liu, C.; Zhu, Y.J.; Nelson, D.R.; Zhou, S.G.; Li, C.F.; Wang, L.Z.; Guo, X.; Sun, Y.Z.; et al. Genome sequence of the model medicinal mushroom *Ganoderma lucidum*. *Nat. Commun.* **2012**, *3*, 913. [CrossRef] [PubMed]

11. Binder, M.; Justo, A.; Riley, R.; Salamov, A.; Lopez-Giraldez, F.; Sjökvist, E.; Copeland, A.; Foster, B.; Sun, H.; Larsson, E.; et al. Phylogenetic and phylogenomic overview of the Polyporales. *Mycologia* **2013**, *105*, 1350-1373. [CrossRef] [PubMed]

12. Zhu, Y.J.; Xu, J.; Sun, C.; Zhou, S.G.; Xu, H.B.; Nelson, D.R.; Qian, J.; Song, J.Y.; Luo, H.M.; Xiang, L.; et al. Chromosome-level genome map provides insights into diverse defense mechanisms in the medicinal fungus *Ganoderma sinense*. *Sci. Rep.* **2015**, *5*, 11087. [CrossRef] [PubMed]

13. Utomo, C.; Tanjung, Z.A.; Aditama, R.; Buana, R.F.N.; Pratomo, A.D.M.; Tryono, R.; Liwang, T. Draft genome sequence of the phytopathogenic fungus *Ganoderma boninense*, the causal agent of basal stem rot disease on oil palm. *Genome Announc.* **2018**, *6*, e00122-18. [CrossRef] [PubMed]

14. Lim, H.J.; Lee, E.H.; Yoon, Y.; Chua, B.; Son, A. Portable lysis apparatus for rapid single-step DNA extraction of *Bacillus subtilis*. *J. Appl. Microbiol.* **2016**, *120*, 379-387. [CrossRef] [PubMed]

15. Fulton, T.M.; Chunwongse, J.; Tanksley, S.D. Microprep protocol for extraction of DNA from tomato and other herbaceous plants. *Plant Mol. Biol. Report.* **1995**, *13*, 207-209. [CrossRef]

16. Grigoriev, I.V.; Nikitin, R.; Haridas, S.; Kuo, A.; Ohm, R.; Otillar, R.; Riley, R.; Salamov, A.; Zhao, X.L.; Korzeniewski, F.; et al. MycoCosm portal: Gearing up for 1000 fungal genomes. *Nucleic Acids Res.* **2014**, *42*, 699-704. [CrossRef] [PubMed]

17. Kuo, A.; Bushnell, B.; Grigoriev, I.V. Fungal genomics: Sequencing and annotation. In *Fungi*; Martin, F., Ed.; Elsevier Academic Press: Cambridge, UK, 2014; pp. 1-52.

18. Stanke, M.; Diekhans, M.; Baertsch, R.; Haussler, D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. *Bioinformatics* **2008**, *24*, 637-644. [CrossRef] [PubMed]

19. Birney, E.; Clamp, M.; Durbin, R. GeneWise and Genomewise. *Genome Res.* **2004**, *14*, 988-995. [CrossRef]

20. Haas, B.J.; Salzberg, S.L.; Zhu, W.; Pertea, M.; Allen, J.E.; Orvis, J.; White, O.; Buell, C.R.; Wortman, J.R. Automated eukaryotic gene structure annotation using eEvidenceModeler and the program to assemble spliced alignments. *Genome Biol.* **2008**, *9*, R7. [CrossRef] [PubMed]

21. Simão, F.A.; Waterhouse, R.M.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. *Bioinformatics* **2015**, *31*, 3210-3212. [CrossRef] [PubMed]

22. Wright, E.S. DECIPHER: Harnessing local sequence context to improve protein multiple sequence alignment. *BMC Bioinform.* **2015**, *16*, 322. [CrossRef] [PubMed]

23. Gu, Z.G.; Gu, L.; Eils, R.; Schlesner, M.; Brors, B. Circlize Implements and enhances circular visualization in R. *Bioinformatics* **2014**, *30*, 2811-2812. [CrossRef]

24. Hage, H.; Miyauchi, S.; Virágh, M.; Drula, E.; Min, B.; Chaduli, D.; Navarro, D.; Favel, A.; Norest, M.; Lesage-Meessen, L.; et al. Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay. *Environ. Microbiol.* **2021**, *23*, 5716-5732. [CrossRef] [PubMed]
25. Stamatakis, A. RAxML Version 8: A tool for phylogenetic analyses and post analysis of large phylogenies. *Bioinformatics* **2014**, *30*, 1312–1313. [CrossRef] [PubMed]

26. Fernandez-Fueyoa, E.; Ruiz-Dueñasa, F.J.; Ferreirab, P.; Floudas, D.; Hibbett, D.S.; Canessa, P.; Larrondo, L.F.; James, T.Y.; Seelenfreund, D.; Lobos, S.; et al. Comparative genomics of *Ceriporiopsis subvermispora* and *Phanerochaete chrysosporium* provide insight into selective ligninolysis. *Proc. Natl. Acad. Sci. USA* **2011**, *109*, 5458–5463. [CrossRef] [PubMed]

27. Floudas, D.; Binder, M.; Riley, R.; Barry, K.; Blanchette, R.A.; Henrissat, B.; Martínez, A.T.; Otillar, R.; Spatafora, J.W.; Yadav, J.S.; et al. The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes. *Science* **2012**, *336*, 1715–1719. [CrossRef] [PubMed]

28. Olson, Å.; Aerts, A.; Asiegbu, F.; Belbahri, L.; Bouzid, O.; Broberg, A.; Canbäck, B.; Coutinho, P.M.; Cullen, D.; Dalman, K.; et al. Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen. *New Phytol.* **2012**, *194*, 1001–1013. [CrossRef] [PubMed]

29. Wu, B.J.; Xu, Z.Y.; Knudson, A.; Carlson, A.; Chen, N.Y.; Kovaka, S.; LaButti, K.; Lipzen, A.; Pennachio, C.; Riley, R.; et al. Genomics and development of *Lentinus tigrinus*: A white-rot wood-decaying mushroom with dimorphic fruiting bodies. *Genome Biol. Evol.* **2018**, *10*, 3250–3261. [CrossRef] [PubMed]

30. Miettinen, O.; Riley, R.; Barry, K.; Cullen, D.; de Vries, R.P.; Hainaut, M.; Hatakka, A.; Henrissat, B.; Hildén, K.; Kuo, R.; et al. Draft genome sequence of the white-rot fungus *Obba rivulosa* 3A-2. *Genome Announc.* **2016**, *4*, e00976-16. [CrossRef]

31. Suzuki, H.; MacDonald, J.; Syed, K.; Salamov, A.; Hori, C.; Aerts, A.; Henrissat, B.; Wiebenga, A.D.; VanKuyk, P.A.; Barry, K.; et al. Comparative genomics of the white-rot fungi, *Phanerochaete carnosa* and *P. chrysosporium*, to elucidate the genetic basis of the distinct wood types they colonize. *BMC Genom.* **2012**, *13*, 444. [CrossRef] [PubMed]

32. Ohm, R.A.; Riley, R.; Salamov, A.; Min, B.; Choi, I.; Grigoriev, I.V. Genomics of wood-degrading fungi. *Fungal Genet. Biol.* **2014**, *72*, 82–90. [CrossRef] [PubMed]

33. Mäkelä, M.R.; Peng, M.; Granchi, Z.; Chin-A-Woeng, T.; Hegi, R.; Pelt, S.I.V.; Ahrendt, S.; Riley, R.; Hainaut, M.; Henrissat, B.; et al. Draft genome sequence of the basidiomycete white-rot fungus *Phlebia centrifuga*. *Genome Announc.* **2018**, *6*, e01414-17. [CrossRef]

34. Kuuskeri, J.; Häkkinen, M.; Laine, P.; Smolander, O.; Tamene, F.; Miettinen, S.; Nousiainen, P.; Kemell, M.; Auvinen, P.; Lundell, T. Time-scale dynamics of proteome and transcriptome of the white-rot fungus *Phlebia radiata*: Growth on spruce wood and decay effect on lignocellulose. *Biotechnol. Biofuels* **2016**, *9*, 192. [CrossRef]

35. Hori, C.; Ishida, T.; Igarashi, K.; Samejima, M.; Suzuki, H.; Master, E.; Ferreira, P.; Ruiz-Dueñas, F.J.; Held, B.; Canessa, P.; et al. Analysis of the *Phlebiopsis gigantea* genome, transcriptome and secretome provides insight into its pioneer colonization strategies of wood. *PLoS Genet.* **2014**, *10*, e1004759. [CrossRef]

36. Varga, T.; Krizsán, K.; Földi, C.; Dima, B.; Sánchez-García, M.; Sánchez-Ramírez, S.; Szöllősi, G.J.; Szarkándi, J.G.; Papp, V.; Albert, L.; et al. Megaphylogeny resolves global patterns of mushroom evolution. *Nat. Ecol. Evol.* **2019**, *3*, 668–678. [CrossRef]

37. Miyauchi, S.; Rancon, A.; Drula, E.; Hage, H.; Chaduli, D.; Favel, A.; Grisel, S.; Henrissat, B.; Herpoël-Gimbert, I.; Ruiz-Dueñas, F.J.; et al. Integrative visual omics of the white-rot fungus *Polyporus brumalis* exposes the biotechnological potential of its oxidative enzymes for delignifying raw plant biomass. *Biotechnol. Biofuels* **2018**, *11*, 201. [CrossRef]

38. Levasseur, A.; Lomascolo, A.; Chabrol, O.; Ruiz-Dueñas, F.J.; Boukhris-Uzan, E.; Piumi, F.; Kües, U.; Ram, A.F.J.; Murat, C.; Haon, M.; et al. The genome of the white-rot fungus *Pycnoporus cinnabarinus*: A basidiomycete model with a versatile arsenal for lignocellulosic biomass breakdown. *BMC Genom.* **2014**, *15*, 486. [CrossRef] [PubMed]

39. Miyauchi, S.; Hage, H.; Drula, E.; Lesage-Meessen, L.; Berrin, J.G.; Navarro, D.; Favel, A.; Chaduli, D.; Grisel, S.; Haon, M.; et al. Conserved white-rot enzymatic mechanism for wood decay in the Basidiomycota genus *Pycnoporus*. *DNA Res.* **2020**, *27*, dsaa011. [CrossRef]

40. Oghenekaro, A.O.; Kovalchuk, A.; Raffaello, T.; Camarero, S.; Gressler, M.; Henrissat, B.; Lee, J.; Liu, M.; Martínez, A.T.; Miettinen, O.; et al. Genome sequencing of *Rigidoporus microporus* provides insights on genes important for wood decay, latex tolerance and interspecific fungal interactions. *Sci. Rep.* **2020**, *10*, 5250. [CrossRef] [PubMed]

41. Granchi, Z.; Peng, M.; Chi-A-Woeng, T.; de Vries, R.P.; Hildén, K.; Mäkelä, M.R. Genome sequence of the basidiomycete white-rot fungus *Trametes pubescens* FBCC735. *Genome Announc.* **2017**, *5*, e01643-16. [CrossRef]

42. Nagy, L.G.; Riley, R.; Tritt, A.; Adam, C.; Daum, C.; Floudas, D.; Sun, H.; Yadav, J.S.; Pangilinan, J.; Larsson, K.H.; et al. Comparative genomics of early-diverging mushroom-forming fungi provides insights into the origins of lignocellulose decay capabilities. *Mol. Biol. Evol.* **2016**, *33*, 959–970. [CrossRef] [PubMed]

43. Tang, J.D.; Perkins, A.D.; Sonstegard, T.S.; Schroeder, S.G.; Burgess, S.C.; Diehla, S.V. Short-read sequencing for genomic analysis of the brown rot fungus *Fibroporia radiculosa*. *Appl. Environ. Microbiol.* **2012**, *78*, 2272–2281. [CrossRef] [PubMed]

44. Gaskell, J.; Kersten, P.; Larrondo, L.F.; Canessa, P.; Martinez, D.; Hibbett, D.; Schmoll, M.; Kubicek, C.P.; Martinez, A.T.; Yadav, I.; et al. Draft genome sequence of a monokaryotic model brown-rot fungus *Postia (Rhodonia) placenta* SB12. *Genom. Data* **2017**, *14*, 21–23. [CrossRef] [PubMed]

45. Kohler, A.; Kuo, A.; Nagy, L.; Morin, E.; Barry, K.; Buscot, F.; Canbäck, B.; Choi, C.; Cichocki, N.; Clum, A.; et al. Convergent losses of decay mechanisms and rapid turnover of symbiosis genes in mycorrhizal mutualists. *Nat. Genet.* **2015**, *47*, 410–415. [CrossRef] [PubMed]

46. Martin, F.; Aerts, A.; Ahrén, D.; Brun, A.; Danchin, E.G.; Duchaussoy, F.; Gibon, J.; Kohler, A.; Lindquist, E.; Pereda, V.; et al. The genome of *Laccaria bicolor* provides insights into mycorrhizal symbiosis. *Nature* **2008**, *452*, 88–92. [CrossRef] [PubMed]
47. Mujic, A.B.; Kuo, A.; Tritt, A.; Lipzen, A.; Chen, C.; Johnson, J.; Sharma, A.; Barry, K.; Grigoriev, I.V.; Spatafora, J.W. Comparative genomics of the ectomycorrhizal sister species *Rhizopogon vinicolor* and *Rhizopogon vesiculosus* (Basidiomycota: Boletales) reveals a divergence of the mating type B Locus. *Genes Genomes Genet.* **2017**, *7*, 1775–1789. [CrossRef] [PubMed]

48. Zuccaro, A.; Lahrmann, U.; Güldener, U.; Langen, G.; Pfiffi, S.; Biedenkopf, D.; Wong, P.; Samans, B.; Grimm, C.; Basiewicz, M.; et al. Endophytic life strategies decoded by genome and transcriptome analyses of the mutualistic root symbiont *Piriformospora indica*. *PLoS Pathog.* **2011**, *7*, e1002290. [CrossRef]

49. Duplessisa, S.; Cuomob, C.A.; Lin, Y.C.; Aerts, A.; Tisserant, E.; Veneault-Fourrey, C.; Joly, D.L.; Hacquard, S.; Amselem, J.; Cantarel, B.L.; et al. Obligate biotrophy features unraveled by the genomic analysis of rust fungi. *Proc. Natl. Acad. Sci. USA* **2011**, *108*, 9166–9171. [CrossRef]

50. Kämper, J.; Kahmann, R.; Bölker, M.; Ma, L.J.; Brefort, T.; Saville, B.J.; Banuett, F.; Kronstad, J.W.; Gold, S.E.; Müller, O.; et al. Insights from the genome of the biotrophic fungal plant pathogen *Ustilago maydis*. *Nature* **2006**, *444*, 97–101. [CrossRef] [PubMed]

51. Emmas, D.; Kelly, S. OrthoFinder: Phylogenetic orthology inference for comparative genomics. *Genome Biol.* **2019**, *20*, 238. [CrossRef] [PubMed]

52. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. *Mol. Biol. Evol.* **2013**, *30*, 772–780. [CrossRef]

53. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. *Bioinformatics* **2009**, *25*, 1972–1973. [CrossRef]

54. Kozlov, A.M.; Darriba, D.; Flouri, T.; Morel, B.; Stamatakis, A. RAxML-NG: A fast, scalable, and user-friendly tool for maximum likelihood phylogenetic inference. *Bioinformatics* **2019**, *35*, 4453–4455. [CrossRef] [PubMed]

55. Darriba, D.; Posada, D.; Kozlov, A.M.; Stamatakis, A.; Morel, B.; Flouri, T. ModelTest-NG: A new and scalable tool for the selection of DNA and protein evolutionary models. *Mol. Biol. Evol.* **2019**, *37*, 291–294. [CrossRef] [PubMed]

56. Payen, T.; Murat, C.; Martin, F. Reconstructing the evolutionary history of gypsy retrotransposons in the Périgord black truffle (Vittad.). *Mycorrhiza* **2016**, *26*, 553–563. [CrossRef] [PubMed]

57. Price, A.L.; Jones, N.C.; Pevzner, P.A. De novo identification of repeat families in large genomes. *Bioinformatics* **2005**, *21*, 351–358. [CrossRef]

58. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

59. Blin, K.; Wolf, T.; Chevrette, M.G.; Lu, X.W.; Schwalen, C.J.; Kautsar, S.A.; Duran, H.G.S.; de Los Santos, E.L.C.; Kim, H.U.K.; Nave, M.; et al. antiSMASH 4.0—Improvements in chemistry prediction and gene cluster boundary identification. *Nucleic Acids Res.* **2017**, *45*, 36–41. [CrossRef] [PubMed]

60. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. *Nucleic Acids Res.* **2021**, *49*, 293–296. [CrossRef] [PubMed]

61. Eddy, S.R. Accelerated profile HMM searches. *PLoS Comput. Biol.* **2011**, *7*, e1002195. [CrossRef] [PubMed]

62. Pellegrin, C.; Morin, E.; Martin, F.M.; Veneault-Fourrey, C. Comparative analysis of secretomes from ectomycorrhizal fungi with an emphasis on small-secreted proteins. *Front. Microbiol.* **2015**, *6*, 1278. [CrossRef]

63. Justo, A.; Miettinen, O.; Floudas, D.; Ortiz-Santana, B.; Sjökvist, E.; Lindner, D.; Nakasone, K.; Niemelä, T.; Larsson, K.H.; Ryvarden, L.; et al. A revised family-level classification of the Polyporales (Basidiomycota). *Fungal Biol.* **2017**, *121*, 798–824. [CrossRef] [PubMed]

64. Wu, F.; Chen, J.J.; Ji, X.H.; Vlasák, J.; Dai, Y.C. Phylogeny and diversity of the morphologically similar polypore genera *Rigidoporus*, *Physisporinus*, *Oxyporus*, and *Leucophellinus*. *Mycologia* **2017**, *109*, 749–765. [CrossRef] [PubMed]

65. Mendes, D.B.; Da Sliva, F.F.; Guarda, P.M.; Almeida, A.F.; de Oliveira, D.P.; Morais, P.B.; Guarda, E.A. Lipolytic enzymes with hydrolytic and esterification activities produced by filamentous fungi isolated from decomposition leaves in an aquatic environment. *Enzym. Res.* **2019**, *8182425*. [CrossRef] [PubMed]

66. Swarbrick, C.M.D.; Nanson, J.D.; Patterson, E.I.; Forwood, J.K. Structure, function, and regulation of thioesterases. *Prog. Lipid Res.* **2020**, *79*, 101036. [CrossRef] [PubMed]

67. Jiang, N.; Hu, S.; Peng, B.; Li, Z.H.; Yuan, X.H.; Xiao, S.J.; Fu, Y.P. Genome of *Ganoderma* species provides insights into the evolution, conifers substrate utilization, and terpene synthesis for *Ganoderma tsugae*. *Front. Microbiol.* **2021**, *12*, 724451. [CrossRef] [PubMed]

68. Baby, S.; Johnson, A.J.; Govindan, B. Secondary metabolites from *Ganoderma*. *Phytochemistry* **2015**, *114*, 66–101. [CrossRef]

69. Pattanayak, S.; Das, S.; Biswal, G. *Ganoderma*: The wild mushroom with wonderful health benefits. *J. Pharmacogn. Phytochem.* **2020**, *9*, 313–316. [CrossRef]

70. Liu, Y.C.; Huang, L.H.; Hu, H.P.; Cai, M.J.; Liang, X.W.; Li, X.M.; Zhang, Z.; Xie, Y.Z.; Xiao, C.; Chen, S.D.; et al. Whole-genome assembly of *Ganoderma leucocontextum* (Ganodermataceae, Fungi) discovered from the Tibetan Plateau of China. *G3 Genes Genomes Genet.* **2021**, *11*, jkab337. [CrossRef] [PubMed]

71. Schnell, N.; Entian, K.D.; Schneider, U.; Götz, F.; Zähner, H.; Kellner, R.; Jung, G. Prepeptide sequence of epidermin, a ribosomally synthesized antibiotic with four sulphide-rings. *Nature* **1988**, *333*, 276–278. [CrossRef] [PubMed]

72. Benson, K.F.; Stamets, P.; Davis, R.; Nally, R.; Taylor, A.; Slater, S.; Jensen, G.S. The mycelium of the *Trametes versicolor* (Turkey tail) mushroom and its fermented substrate each show potent and complementary immune activating properties in vitro. *BMC Complement. Altern. Med.* **2019**, *19*, 342. [CrossRef]
73. Gao, X.; Qi, J.Y.; Ho, C.T.; Li, B.; Mu, J.J.; Zhang, Y.T.; Hu, H.P.; Mo, W.P.; Chen, Z.Z.; Xie, Y.Z. Structural characterization and immunomodulatory activity of a water-soluble polysaccharide from *Ganoderma leucocontextum* fruiting bodies. *Carbohydr. Polym.* **2020**, *249*, 116874. [CrossRef] [PubMed]

74. Im, K.H.; Nguyen, T.K.; Choi, J.; Lee, T.S. In vitro antioxidant, anti-diabetes, anti-dementia, and inflammation inhibitory effect of *Trametes pubescens* fruiting body extracts. *Molecules* **2016**, *21*, 639. [CrossRef] [PubMed]

75. Baldrian, P.; Valášková, V. Degradation of cellulose by basidiomycetous fungi. *FEMS Microbiol. Rev.* **2008**, *32*, 501–521. [CrossRef] [PubMed]

76. Lebreton, A.; Bonnardel, F.; Dai, Y.C.; Imberty, A.; Martin, F.M.; Lisacek, F. A comprehensive phylogenetic and bioinformatics survey of lectins in the fungal kingdom. *J. Fungi* **2021**, *7*, 453. [CrossRef] [PubMed]

77. Hage, H.; Rosso, M.N. Evolution of fungal carbohydrate-active enzyme portfolios and adaptation to plant cell-wall polymers. *J. Fungi* **2021**, *7*, 185. [CrossRef] [PubMed]

78. Syed, K.; Nelson, D.R.; Riley, R.; Yadav, J.S. Genomewide annotation and comparative genomics of cytochrome P450 monooxygenases (P450s) in the polypore species *Bjerkandera adusta*, *Ganoderma* sp. and *Phlebia brevispora*. *Mycologia* **2013**, *105*, 1445–1455. [CrossRef] [PubMed]

79. Debets, F.; Yang, X.; Griffiths, A.J. Vegetative incompatibility in *Neurospora*: Its effect on horizontal transfer of mitochondrial plasmids and senescence in natural populations. *Curr. Genet.* **1994**, *26*, 113–119. [CrossRef]

80. Nuss, D. Hypovirulence: Mycoviruses at the fungal–plant interface. *Nat. Rev. Microbiol.* **2005**, *3*, 632–642. [CrossRef]

81. Madej, M.G. Comparative sequence–function analysis of the major facilitator superfamily: The “Mix-and-Match”. *Methods Enzymol.* **2015**, *557*, 521–549. [CrossRef]

82. Riechmann, J.L.; Heard, J.; Martin, G.; Reuber, L.; Jiang, C.; Keddie, J.; Adam, L.; Pineda, O.; Ratcliffe, O.J.; Samaha, R.R.; et al. *Arabidopsis* transcription factors: Genome-wide comparative analysis among eukaryotes. *Science* **2000**, *290*, 2105–2110. [CrossRef] [PubMed]

83. Wang, L.; Huang, Q.H.; Zhang, L.L.; Wang, Q.F.; Liang, L.; Liao, B.S. Genome-wide characterization and comparative analysis of MYB transcription factors in *Ganoderma* species. *Genes Genomes Genet.* **2020**, *10*, 2653–2660. [CrossRef] [PubMed]