Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes’ intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

24Mushroom-forming fungi in the order Agaricales represent an independent origin of 25 bioluminescence in the tree of life, yet the diversity, evolutionary history, and timing of 26 the origin of fungal luciferases remain elusive. We sequenced the genomes and 27transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage 28comprising the majority of bioluminescent fungi. We show that bioluminescence 29 evolved in the common ancestor of Mycena spp. and the marasmioid clade of 30Agaricales and was maintained through at least 160 million years of evolution. We 31revealed Mycena exhibit two-speed genomes and resolve how the luciferase cluster was 32 derived by duplication and translocation, frequently rearranged and lost in most Mycena 33 species, but conserved in the Armillaria lineage. Luciferase cluster members were co-34 expressed across developmental stages, with highest expression in fruiting body caps 35and stipes, suggesting fruiting-related adaptive functions. Our results contribute to 36 understanding a de novo origin of bioluminescence and the corresponding gene cluster 37 in a diverse group of enigmatic fungal species. 38