Transcriptome Profiling Reveals Candidate Genes Related to Stipe Gradient Elongation of Flammulina filiformis

Stipe gradient elongation is an important and remarkable feature in the development of most mushroom fruiting bodies. However, its molecular mechanism has rarely been described. Here, the decreasing trend of stipe elongation and increasing trend of cell length in a gradient from the top to the base of the stipe were determined in a model basidiomycete mushroom: Flammulina filiformis. According to RNA-seq results, 1409 differentially expressed genes (DEGs) were identified among elongation region (ER), transition region (TR), and stable region (SR) samples, including 26 transcription factors (TFs). Based on Short Time-series Expression Miner (STEM) clustering of DEGs, clusters 1 and 3, with obvious expression trends that were consistent with or in contrast to the elongation rate, were screened. The cluster 1 DEGs were mainly involved in the GO cellular component category and KEGG genetic information processing class; however, the cluster 3 DEGs were mainly involved in metabolic processes. Furthermore, qRT-PCR confirmed that key genes of the long-chain fatty acid synthesis pathway were involved in stipe gradient elongation and regulated by NADPH oxidase-derived ROS signaling molecules. These findings provide an essential basis for understanding the molecular mechanism of stipe gradient elongation.