Secondary cell wall (SCW) thickening is a major cellular developmental stage determine wood structure and properties. Although the molecular regulation of cell wall deposition during tracheary element differentiation has been well established in primary growth systems, less is known about the gene regulatory processes involved in the multi-layered SCW thickening of mature trees. Using third-generation (long-read single-molecule real-time, SMRT) and second-generation (short-read sequencing by synthesis, SBS) sequencing methods, we established a Pinus bungeana transcriptome resource with comprehensive functional and structural annotation for the first time. Using these approaches, we generated high spatial resolution datasets for the vascular cambium, xylem expansion regions, early SCW thickening, late SCW thickening, and mature xylem tissues of 71-year-old Pinus bungeana trees. A total of 79,390 non-redundant transcripts, 31,808 long non-coding RNAs, and 5,147 transcription factors were annotated and quantified in different xylem tissues at all growth and differentiation stages. Furthermore, using this high spatial resolution dataset, we established a comprehensive transcriptomic profile and found that members of the NAC, WRKY, SUS, CESA, and LAC gene families are major players in early SCWs formation in tracheids, whereas members of the MYB and LBD transcription factor families are highly expressed during late SCWs thickening. Our results provide new molecular insights into the regulation of multi-layered SCW thickening in conifers. The high spatial resolution datasets provided can serve as important gene resources for improving softwoods.

Ray parenchyma cells are involved in the initiation of heartwood formation. The position within a ray influences the timing of ray parenchyma cell differentiation and function; however, there is little information concerning the positional influence on the cellular changes of ray parenchyma cells from sapwood and heartwood. In this study, radial variations in morphology, size, and ultrastructure of ray parenchyma cells were studied by combined transmission electron microscopy and optical microscopy. Results showed that cellular traits of ray parenchyma cells in Populus tomentosa were all affected by both radial position in the secondary xylem and position within a ray. Specifically, radial variations in cellular traits were more evident in isolation cells, which were not adjacent to vessel elements. Both cell length and cell width/length ratio of isolation cells were bigger than contact cells, which contacted adjacent vessel elements via pits. Moreover, the secondary wall thickening and lignification of contact cells developed in the current-year xylem, much earlier than isolation cells. Secondary walls in contact cells were in a polylamellate structure with a protective layer on the inner side. No alteration in the ultrastructure of contact cells occurred in the sapwood-heartwood transition zone, except that most contact cells died. By contrast, in the transition zone, isolation cells still lived. A thin secondary wall began to deposit on the thick primary wall of isolation cells, with two isotropic layers on the inner side of the primary wall and secondary wall respectively being characteristic. Meanwhile, starch grains in isolation cells were depleted, and dark polyphenolic droplets lost their spherical shape and flowed together. Furthermore, the intercellular spaces of isolation cells became densified in the transition zone. Overall, cellular changes suggested that the positional information of ray parenchyma cells appeared to be an important factor in the transformation from sapwood to heartwood. Unlike contact cells, isolation cells were more elongated, specialized in radial transport, had a delayed formation of secondary walls, and were involved in the synthesis of heartwood substances. Our result promotes the elucidation of the involvement of xylem rays in heartwood formation.