Abstract
<jats:p>Background and aims: Perianal fistulizing Crohn's disease (pCD) affects 20% of patients with Crohn's disease (CD) and severely impacts quality of life. Current therapies fail to provide sustained relief, and the molecular underpinnings of pCD remain poorly understood. This study aimed to elucidate the molecular landscape of perianal fistulas through multiomic profiling. Methods: Paired fistula-tract and adjacent rectal mucosal biopsies were collected from 63 patients (48 with pCD, 15 with cryptoglandular fistula [CPTGL]). Longitudinal sampling generated 101 unique molecular profiles, comprising of both RNA sequencing (RNA-seq) and 16S rRNA sequencing from both tissue sites, followed by integrative multiomics and gene network analyses. Results: Unsupervised clustering revealed three patient clusters primarily defined by host gene expression, with minimal contribution from microbial profiles. Fistulae in clusters 1 and 2 showed strong immune activation and epithelial-mesenchymal transition (EMT). In contrast, cluster 3 fistulae displayed keratinization and metabolic reprogramming resembling psoriatic skin, together with reduced JAK-STAT signalling. Cluster 3 was enriched for patients classified as TOpClass:2a, who are more suitable for surgical repair (p = 0.02). Conversely, cluster 2, characterized by rectal keratinization and EMT in both fistula and rectum, showed the highest MRI inflammatory-mass score (p = 0.02) and a greater risk of subsequent ileostomy (Kaplan-Meier; p = 0.008). An independent RNAseq dataset validated the keratinization signature in a subset of pCD fistulae. Conclusion: Integrated multiomic analysis identified distinct molecular subtypes of pCD with surgical and therapeutic relevance. These findings refine the molecular understanding of pCD and support a precision medicine approach.</jats:p>