Cumulative disease-course phenotypes in Crohn's disease: a retrospective single center FAMD clustering study

Background and aimsCrohns disease (CD) is a highly heterogeneous disease. Patients with similar characteristics at baseline may end up on very different therapeutic and surgical paths. We aimed to identify disease-course trajectory phenotypes from a mix of clinical, treatment, surgical and comorbidity variables, and evaluate whether information already available at diagnosis could anticipate which trajectory a patient would follow. MethodsWe analysed a retrospective cohort of 98 unique CD patients enrolled at the IRCCS Ca Granda Policlinico Hospital, Milan. Factor analysis of mixed data (FAMD) was applied to 30 disease-course variables spanning demographic and Montreal descriptors, extra-intestinal manifestation burden, individual advanced-therapy drug exposures, surgical burden and comorbidity/symptom burden; clusters were then ranked by treatment burden using a line-weighted therapy score built separately from drug line positions and therapeutic-class weights. Clustering itself was carried out in the retained FAMD coordinate space, and cluster robustness was checked with silhouette metrics and bootstrap Jaccard stability. Baseline-associated variables were assessed separately, while advanced-therapy episodes were reconstructed from the therapy timeline. ResultsThree trajectory phenotypes emerged: C1 (n=56), C2 (n=31) and C3 (n=11). C3 was characterised by younger age at diagnosis, a heavier extra-intestinal manifestation burden, more intensive and prolonged therapy, greater biologic exposure and a higher cumulative number of surgeries, with comorbidity burden rising across the gradient. Montreal behaviour and comorbidity burden emerged as the strongest baseline correlate of trajectory membership, ahead of age at diagnosis and extra-intestinal manifestation burden; penetrating (B3) behaviour, notably, was actually more common in C1 than in C3. A model built on the diagnosis-available variables independently linked to trajectory membership (age at diagnosis, disease location and Montreal behaviour) showed moderate discrimination for the C3 trajectory (area under the curve, AUC=0.80, 95% CI 0.67-0.93). A sensitivity analysis that also added extra-intestinal manifestation count and comorbidity burden, whose timing relative to diagnosis we could not confirm, pushed discrimination to AUC=0.90. ConclusionsFAMD-space clustering identified stable, clinically interpretable CD-course trajectory phenotypes, including a small, high-burden C3 subgroup. This C3 trajectory, small, early-onset and treatment-refractory, shares clinical features that overlap with the difficult-to-treat CD phenotype, though C3 was derived by unsupervised clustering rather than validated against that frameworks own criteria; variables available at or near diagnosis showed exploratory discrimination for the high-burden C3 trajectory, though external validation and careful timing verification are needed before any clinical use.