Secondary xylem produced by stem secondary growth is the main source of tree biomass and possesses great economic and ecological value in papermaking, construction, biofuels, and the global carbon cycle. The secondary xylem formation is a complex developmental process, and the underlying regulatory networks and potential mechanisms are still under exploration. In this study, using hybrid poplar (Populus alba x Populus glandulosa clone 84K) as a model system, we first ascertained three representative stages of stem secondary growth and then investigated the regulatory network of secondary xylem formation by joint analysis of transcriptome and miRNAs. Notably, 7507 differentially expressed genes (DEGs) and 55 differentially expressed miRNAs (DEMs) were identified from stage 1 without initiating secondary growth to stage 2 with just initiating secondary growth, which was much more than those identified from stage 2 to stage 3 with obvious secondary growth. DEGs encoding transcription factors and lignin biosynthetic enzymes and those associated with plant hormones were found to participate in the secondary xylem formation. MiRNA-target analysis revealed that a total of 85 DEMs were predicted to have 2948 putative targets. Among them, PagmiR396d-PagGRFs, PagmiR395c-PagGA2ox1/PagLHW/PagSULTR2/PagPolyubiquitin 1, PagmiR482d-PagLAC4, PagmiR167e-PagbHLH62, and PagmiR167f/g/h-PagbHLH110 modules were involved in the regulating cambial activity and its differentiation into secondary xylem, cell expansion, secondary cell wall deposition, and programmed cell death. Our results give new insights into the regulatory network and mechanism of secondary xylem formation.
