Comparative transcriptome and WGCNA reveal key genes involved in lignocellulose degradation in Sarcomyxa edulis

The developmental transcriptomes of Sarcomyxa edulis were assessed to explore the molecular mechanisms underlying lignocellulose degradation. Six stages were analyzed, spanning the entire developmental process: growth of mycelium until occupying half the bag (B1), mycelium under low-temperature stimulation after occupying the entire bag (B2), appearance of mycelium in primordia (B3), primordia (B4), mycelium at the harvest stage (B5), and mature fruiting body (B6). Samples from all six developmental stages were used for transcriptome sequencing, with three biological replicates for all experiments. A co-expression network of weighted genes associated with extracellular enzyme physiological traits was constructed using weighted gene co-expression network analysis (WGCNA). We obtained 19 gene co-expression modules significantly associated with lignocellulose degradation. In addition, 12 key genes and 8 kinds of TF families involved in lignocellulose degradation pathways were discovered from the four modules that exhibited the highest correlation with the target traits. These results provide new insights that advance our understanding of the molecular genetic mechanisms of lignocellulose degradation in S. edulis to facilitate its utilization by the edible mushroom industry.