Ectomycorrhiza (ECM) is a symbiotic relation between plant and fungi that is essential for nutrient uptake of many stand forming trees. There are two conflicting views about the evolution of ECM in fungi suggesting (1) relatively few transitions to ECM followed by reversals to non-ECM, or (2) many independent origins of ECM and no reversals. In this study, we compare these, and other, hypotheses and test the impact of different models on inference. We assembled a dataset of five marker gene sequences (nuc58, nucLSU, nucSSU, rpb1, and rpb2) and 2,174 fungal taxa covering the three subphyla: Agaricomycotina, Mucoromycotina and Pezizomycotina. The fit of different models, including models with variable rates in clades or through time, to the pattern of ECM fungal taxa was tested in a Bayesian framework, and using AIC and simulations. We find that models implementing variable rates are a better fit than models without rate shift, and that the conclusion about the relative rate between ECM and non-ECM depend largely on whether rate shifts are allowed or not. We conclude that standard constant-rate ancestral state reconstruction models are not adequate for the analysis of the evolution of ECM fungi, and may give contradictory results to more extensive analyses.

The genomes of ectomycorrhizal (ECM) fungi have a reduced number of genes encoding Carbohydrate-Active EnZymes (CAZymes), expansions in transposable elements (TEs) and small secreted proteins (SSPs) compared with saprotrophs. Fewer genes for specific peptidases and lipases in ECM fungi are also reported. It is unclear whether these changes occur at the shift to the ECM habit or are more gradual throughout the evolution of ECM lineages.

We generated a genomic dataset of 20 species in the ECM lineage Inocybaceae and compared them with six saprotrophic species.

Inocybaceae genomes have fewer CAZymes, peptidases, lipases, secondary metabolite clusters and SSPs and higher TE content than their saprotrophic relatives. There was an increase in the rate of gene family evolution along the branch with the transition to the ECM lifestyle. This branch had very high rate of evolution in CAZymes and had the largest number of contractions. Other significant changes along this branch included expansions in transporters, transposons-related genes and communication genes such as fungal kinases.

There is a high concentration of changes in proximity to the transition to the ECM lifestyle, which correspond to the identified key changes for the gain of this lifestyle.