Comparative analysis of differential gene expression reveals novel insights into the heteroblastic foliage functional traits of Pinus massoniana seedlings.

Pinus massoniana needles, rich in medicinal polysaccharides and flavonoids, undergo heteroblastic foliage, transitioning from primary needles (PN) to secondary needles (SN) during growth, resulting in altered functional traits. Despite its significance, the molecular regulatory mechanisms governing these traits remain unclear. This study employs Iso-Seq and RNA-Seq analyses to explore differentially expressed genes (DEGs) associated with functional traits throughout the main growth season of heteroblastic foliage. Co-expression network analysis identified 34 hub genes and 17 key transcription factors (TFs) influencing light-harvesting antenna, photosystem I and II, crucial in photosynthesis regulation. Additionally, 14 genes involved in polysaccharide metabolism pathways, synthesizing sucrose, glucose, UDP sugars, and xylan, along with four genes in flavonoid biosynthesis pathways, regulating p-coumaroyl-CoA, quercetin, galangin, and myricetin production, exhibited differential expression between PN and SN. Further analysis unveils a highly interconnected network among these genes, forming a pivotal cascade of TFs and DEGs. Therefore, heteroblastic changes significantly impact needle functional traits, potentially affecting the pharmacological properties of PN and SN. Thus, these genomic insights into understanding the molecular-level differences of heteroblastic foliage, thereby establishing a foundation for advancements in the pharmaceutical industry related to needle-derived products.