Refractory coeliac disease (RCD) is a severe complication of coeliac disease despite strict gluten avoidance. RCD type 2 (RCD2) is characterized by aberrant intraepithelial lymphocytes (IELs) [1][2] and carries a high lymphoma risk, while RCD type 1 (RCD1) remains a diagnosis of exclusion because of the lack of aberrant IELs despite having villous atrophy. We asked whether RCD1 represents a distinct biological state and whether persistent inflammation is accompanied by coordinated immune and epithelial remodelling.
Duodenal biopsies from 24 participants including 6 healthy controls, 7 active coeliac disease (ACD), 9 RCD1 and 2 RCD2 patients were profiled using single-cell RNA sequencing, CITE-seq and immune-receptor sequencing of CD45+ cells, with single-cell RNA sequencing of EPCAM+ epithelial cells. Spatial transcriptomics was performed in 8 biopsies.
Aberrant IELs that carry lymphoma-driver mutations occupied multiple transcriptionally distinct states. We discovered large molecular heterogeneity in IEL populations, with expansion of both αβ and γδ clones increasing with disease severity. IEL populations carried a shared remodelling programme in RCD1, including GZMK-expressing memory-like and cytotoxic states. Within epithelial cells, we found highly heterogeneous enterocyte, and transit-amplifying (TA) clusters. Disease-associated epithelial cells displayed progressive regenerative and immune-recognition programmes, including REG3A and HLA-E upregulation. Enteroendocrine (EEC) population increased with disease severity, with NEUROD1 activating regulon activity and upregulation of nociceptive channel in RCD1. Pseudotime analysis identified distinct epithelial differentiation trajectory in RCD1. Spatial analysis further confirmed the presence of REG3A+ and HLA-E+ enterocytes and NK-like lymphocyte clones.
Single-cell multi-omics analysis revealed a distinct molecular signature of RCD1 compared to ACD characterized by coordinated immune, epithelial reprogramming and clonal expansion. These findings provide a biological framework for distinguishing RCD1 from ACD. Analysis of larger, longitudinal RCD1 and RCD2 cohorts will be required to determine whether these conditions represent distinct disease states or stages along a shared pathogenic continuum.