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.

In recent years, eukaryotic long non-coding RNAs (lncRNAs) have been identified as important factors involved in a wide variety of biological processes, including histone modification, alternative splicing and transcription enhancement. The expression of lncRNAs is highly tissue-specific and is regulated by environmental stresses. Recently, a large number of plant lncRNAs have been identified, but very few of them have been studied in detail. Furthermore, the mechanism of lncRNA expression regulation remains largely unknown. Arabidopsis HISTONE DEACETYLASE 6 (HDA6) and LSD1-LIKE 1/2 (LDL1/2) can repress gene expression synergistically by regulating H3Ac/H3K4me. In this research, we performed RNA-seq and ChIP-seq analyses to further clarify the function of HDA6-LDL1/2. Our results indicated that the global expression of lncRNAs is increased in hda6/ldl1/2 and that this increased lncRNA expression is particularly associated with H3Ac/H3K4me2 changes. In addition, we found that HDA6-LDL1/2 is important for repressing lncRNAs that are non-expressed or show low-expression, which may be strongly associated with plant development. GO-enrichment analysis also revealed that the neighboring genes of the lncRNAs that are upregulated in hda6/ldl1/2 are associated with various developmental processes. Collectively, our results revealed that the expression of lncRNAs is associated with H3Ac/H3K4me2 changes regulated by the HDA6-LDL1/2 histone modification complex.

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