Transcriptome and miRNAome analysis reveals components regulating tissue differentiation of bamboo shoots

Primary thickening determines bamboo yield and wood property. However, little is known about the regulatory networks involved in this process. This study identified a total of 58,652 genes and 150 miRNAs via transcriptome and small RNA sequencing using the underground thickening shoot samples of wild-type (WT) Moso bamboo (Phyllostachys edulis) and a thick wall (TW) variant (P. edulis “Pachyloen”) at five developmental stages (WTS1/TWS1-WTS5/TWS5). A total of 14,029 (65.17%) differentially expressed genes and 68 (45.33%) differentially expressed miRNAs were identified from the WT, TW, and WTTW groups. The first two groups were composed of four pairwise combinations, each between two successive stages (WTS2/TWS2_versus_WTS1/TWS1, WTS3/TWS3_versus_WTS2/TWS2, WTS4/TWS4_versus_WTS3/TWS3, and WTS5/TWS5_versus_WTS4/TWS4), and the WTTW group was composed of five combinations, each between two relative stages (TWS1–5_versus_WTS1–5). Additionally, among the phytohormones, zeatin showed more remarkable changes in concentrations than indole-3-acetic acid, gibberellic acid, and abscisic acid throughout the five stages in the WT and the TW groups. Moreover, 125 cleavage sites were identified for 387 miRNA–mRNA pairs via degradome sequencing (P