A single-nucleus and spatial transcriptomic atlas of poplar leaves reveals the regulation of leaf polarity and cuticle deposition.
Leaf adaxial-abaxial polarity is fundamental for plant morphogenesis and environmental adaptation through asymmetric cell differentiation. Emerging evidence reveals dorsoventral metabolic gradients act downstream of transcriptional networks to fine-tune cellular specialization. While conserved transcription factors (e.g., HD-ZIP III and KANADI) establish initial polarity, the molecular networks driving position-specific cellular differentiation and their integration with metabolic adaptation remain unclear. Leveraging single-nucleus and spatial transcriptomics, we resolve major cell classes (mesophyll, epidermal, and vascular-associated) and their adaxial-abaxial subtypes, revealing dorsoventral polarity in transcriptional profiles and metabolic pathways. Adaxial cells are enriched in phenylpropanoid/flavonoid biosynthesis, while abaxial cells show preferential activation of stress and hormone signaling. Notably, we identify MYC2 as a key regulator of adaxial cuticle biosynthesis, binding to promoters of lipid biosynthetic and transport genes (e.g., CER10 and LTPG1) and promoting cuticle thickening. Our study uncovers how positional identity shapes transcriptional and metabolic polarity in leaves, with MYC2 emerging as a central regulator coordinating organ-specific adaptations. These findings provide insights into the spatial regulation of plant development and stress resilience, offering potential strategies for engineering stress-tolerant woody crops.