Abstract
Pleurotus tuoliensis (Pt) and P. eryngii var. eryngii (Pe) are important edible mushrooms. The epigenetic and gene expression signatures characterizing major developmental transitions in these two mushrooms remain largely unknown. Here, we report global analyses of DNA methylation and gene expression in both mushrooms across three major developmental transitions, from mycelium to primordium and to fruit body, by whole-genome bisulfite sequencing (WGBS) and RNA-seq-based transcriptome profiling. Our results revealed that in both Pt and Pe the landscapes of methylome are largely stable irrespective of genomic features, e.g., in both protein-coding genes and transposable elements (TEs), across the developmental transitions. The repressive impact of DNA methylation on expression of a small subset of genes is likely due to TE-associated effects rather than their own developmental dynamics. Global expression of gene orthologs was also broadly conserved between Pt and Pe, but discernible interspecific differences exist especially at the fruit body formation stage, and which are primarily due to differences in trans-acting factors. The methylome and transcriptome repertories we established for the two mushroom species may facilitate further studies of the epigenetic and transcriptional regulatory mechanisms underpinning gene during development in Pleurotus and related genera.