Abstract
Oral mucosa exhibits region-specific keratinization, essential for periodontal health, yet the spatial and molecular mechanisms driving these differences remain poorly understood. This study aimed to generate a high-resolution spatial transcriptomic atlas of the human oral mucosa around the mucogingival junction, to reveal stromal-epithelial interactions, that distinguish keratinized from non-keratinized programs. Formalin-fixed paraffin-embedded specimens from the mucogingival junction area of two healthy donors were analyzed with the 10 × Genomics Visium HD platform, yielding two keratinized and two non-keratinized regions. Spatial clustering, pseudotime trajectory inference, cell-type integration with a single-cell reference, and ligand-receptor network analysis were applied to delineate epithelial and stromal compartments. Sixteen reproducible clusters, recapitulating tissue architecture, were identified and revealed distinct transcriptional signatures, distinguishing gingiva from lining mucosa. Pseudotime analysis revealed bifurcating epithelial lineages, originating from a shared basal progenitor layer into keratinized and non-keratinized programs. Gingival keratinization was driven by stromal collagen ligands (COL1A1, COL1A2, COL6A1, COL6A2) engaging epithelial receptors (CD44, SDC1), further reinforced within the epithelium by desmosomal adhesion via DSG1-DSC2/3. Gingival keratinization emerges from integrated stromal collagen signaling and epithelial adhesion. This spatially resolved framework advances understanding of oral mucosal specialization and provides a foundation for biologically guided regenerative therapies.