Abstract
Celiac disease remains underdiagnosed despite an established diagnostic pathway, reflecting phenotypic heterogeneity, delayed recognition, and dependence on resource-intensive testing. This narrative review synthesizes evidence on artificial intelligence and phenomics in celiac disease (CD) across four operational layers: electronic health record phenotyping, Human Phenotype Ontology-based semantic encoding, machine-learning pre-screening from routine clinical data, and deep-learning-assisted histopathology. A targeted literature search of PubMed/MEDLINE, Embase, and Google Scholar covered publications from January 2010 through December 2025, using combinations of CD/coeliac disease with artificial intelligence, machine learning, deep learning, phenomics, Human Phenotype Ontology, computable phenotype, electronic health records, and natural language processing, supplemented by targeted searches and citation chaining. We present a curated minimum viable CD phenome and discuss clinical actionability, age-specific considerations, and current pediatric validation gaps. Selected models have demonstrated promising CD detection or pre-screening performance in specific datasets, but prospective, multicenter evidence with external validation remains insufficient to establish clinical utility. Artificial intelligence should augment expert clinical care rather than replace it.
Keywords
Celiac disease,
Artificial intelligence,
Phenomics,
Human Phenotype Ontology,
Deep learning,
Computable phenotype,
Electronic health records,
Precision medicine.
Introduction
Celiac disease (CD) is a chronic immune-mediated enteropathy triggered by dietary gluten in genetically susceptible individuals, affecting an estimated 1.4% of the global population by serology and 0.7% by biopsy confirmation.1 Despite its prevalence, more than 80% of cases remain undiagnosed.2,3 Diagnostic delays have been reported across healthcare settings and may be prolonged in both high- and low-resource environments.4-8
Diagnostic delay has been associated with clinical consequences, including osteoporosis, iron-deficiency anemia, and neurological complications.5,6 Persistent underdiagnosis reflects the phenotypic heterogeneity of CD, spanning classical malabsorption to silent or extraintestinal presentations, and the resource-intensive nature of diagnostic evaluation.9,10 The challenge is particularly relevant in pediatric patients, who may present with short stature (HP:0004322), failure to thrive (HP:0001508), and dental enamel defects, while adults may present predominantly with extraintestinal or minimally symptomatic disease.11 Potential CD is best described as persistently positive CD-specific serology in the presence of normal or non-atrophic small-intestinal mucosa; precise diagnostic and management criteria depend on the applicable guideline and clinical context.11
Artificial intelligence (AI) and phenomics may help address these bottlenecks. Machine learning (ML) may identify pre-diagnostic patterns in routine clinical data,12 while deep learning (DL) has shown promising performance for automated CD detection from digitized biopsy images.13 Standardized ontologies such as the Human Phenotype Ontology (HPO) enable computable, interoperable phenotype representations that can bridge clinical records, research databases, and AI pipelines.14,15 This review synthesizes current evidence across these domains and proposes a four-layer implementation framework for AI-assisted CD diagnostics, with explicit attention to clinical actionability, age-specific phenotypic considerations, and equitable applicability.
A targeted narrative literature search was conducted across PubMed/MEDLINE, Embase, and Google Scholar for publications from January 2010 through December 2025 addressing AI, ML, DL, phenomics, or computable phenotyping in CD; earlier publications were considered without date restriction for foundational HPO and EHR-phenotyping methodology. Search terms included combinations of “CD” or “coeliac disease” with “artificial intelligence,” “machine learning,” “deep learning,” “phenomics,” “Human Phenotype Ontology,” “computable phenotype,” “electronic health records,” and “natural language processing.” Additional targeted searches addressed extreme gradient boosting (XGBoost), convolutional neural networks, vision transformers, histopathology, potential CD, and pre-diagnostic screening. Citation chaining from high-yield publications supplemented the database searches. Studies were considered when they reported original AI-model performance data relevant to CD or methodological approaches directly applicable to CD phenotyping. Conference abstracts without full-text data, case reports, and editorials without original quantitative findings were excluded. The literature was reviewed by the authors for relevance and applicability; no formal risk-of-bias or methodological quality-assessment tool was applied.
Computable phenotyping of CD
Limitations of ICD-based case ascertainment
The International Classification of Diseases (ICD) has historically been used to identify CD in administrative and research databases. ICD-10 code K90.0 has shown high specificity in some settings but may have limited sensitivity when diagnoses are recorded under related codes or when extraintestinal manifestations predominate.16 Evidence on coding performance is heterogeneous, and a specific 15–30% misclassification rate for CD cohorts is not sufficiently supported by the cited literature; accordingly, we do not assign a universal misclassification rate.
Phecodes, which aggregate related ICD codes into clinically meaningful phenotype categories, offer a partial solution. The phecode for CD (phecode 557.1) combines K90.0 with related codes for malabsorption and gluten sensitivity.17,18 However, phecodes remain dependent on the completeness and accuracy of underlying ICD coding and do not capture laboratory, pathological, or narrative data.
Multimodal electronic health record phenotyping
More sophisticated computable phenotypes for CD integrate multiple electronic health record (EHR) data streams: structured diagnosis codes, laboratory values (tissue transglutaminase immunoglobulin A (tTG-IgA), hemoglobin, ferritin, albumin), medication records, pathology reports, and unstructured clinical notes processed by natural language processing (NLP).19-21 Algorithmic phenotyping frameworks such as PheNorm and PheVis have been applied to autoimmune and other conditions, 22,23 with approaches that may be transferable to CD phenotyping. The CALIBER platform in the United Kingdom has demonstrated high-fidelity EHR phenotyping across numerous conditions using linked primary care, hospital, and registry data.24 Ludvigsson et al.25 demonstrated that a computerized algorithm applied to electronic medical record data could identify individuals requiring CD testing with high specificity. The integration of pathology-report NLP with structured EHR data may capture histological and longitudinal dimensions that are not available from ICD codes alone.21,26
The Human Phenotype Ontology as a semantic framework for CD
HPO architecture and the 2024 update
The Human Phenotype Ontology (HPO) is a structured, hierarchical vocabulary of human phenotypic abnormalities, originally described by Robinson et al.27 Successive updates have expanded its coverage substantially; the 2024 release comprises more than 18,000 terms and over 170,000 disease–phenotype annotations across thousands of diseases.14 HPO terms are organized as a directed acyclic graph with major subhierarchies covering abnormalities of the digestive, hematological, metabolic, and nervous systems.28,29 The 2024 update also expanded multilingual and global resources, supporting broader use in diverse populations in which CD epidemiology and human leukocyte antigen (HLA) distributions may differ.14,30,31
HPO-based phenotyping pipelines and age-specific considerations
Several computational pipelines exploit HPO for clinical decision support and differential diagnosis. Encoding clinical data with HPO enables similarity-based matching between patient phenotype profiles and disease-specific HPO annotations, a strategy established in rare-disease diagnostics and potentially applicable to CD.32,33 For CD, HPO can represent a broad phenotypic spectrum, including villous atrophy (HP:0011473), malabsorption (HP:0002024), and extraintestinal manifestations. The proposed CD phenome should therefore be regarded as a curated conceptual feature set rather than a validated disease-specific frequency model.
An important consideration for HPO-based AI models is the distinction between pediatric and adult phenotypic presentations. In children, CD may manifest through growth failure, failure to thrive, and dental enamel defects. In adults, extraintestinal or minimally symptomatic presentations are common, including iron-deficiency anemia, elevated transaminases, and peripheral neuropathy. This phenotypic divergence suggests that HPO feature weighting, training-cohort composition, and risk thresholds may require age-stratified calibration. The cited ML pre-screening study by Dreyfuss et al.12 included adolescents and adults, and its performance has not been independently validated in an exclusively pediatric population.
The PhenoBrain system illustrates the potential of HPO-based AI in rare-disease differential diagnosis.34 Evaluated across 2,271 cases covering 431 rare diseases, the system reported top-3 recall of 0.613 and top-10 recall of 0.813 in the reported physician-comparison setting.34 These findings were obtained in a non-celiac rare-disease cohort and therefore provide indirect methodological support only; they do not establish diagnostic performance in CD and should not be interpreted as evidence that HPO-based AI outperforms specialists in CD.
Curated HPO terms for the CD phenome
Table 1 presents a curated set of HPO terms intended as a minimum viable CD phenome across major clinical domains. The terms were selected conceptually to represent gastrointestinal, hematological, metabolic/growth, neurological, dermatological, and reproductive features relevant to CD. Because the cited literature does not provide a validated CD-specific frequency estimate or machine-learning assessment of feature utility for this exact 18-term set, the table deliberately omits a frequency column and avoids claims that any term is obligatory, universal, or independently validated as a predictive feature. The set should therefore be considered a hypothesis-generating semantic layer requiring formal validation before clinical deployment.
| HPO term | HPO ID | Clinical domain | Clinical relevance |
|---|
| Villous atrophy | HP:0011473 | Gastrointestinal | Important histological feature of established CD; not universal because diagnosis may be established without biopsy in selected pediatric pathways |
| Malabsorption | HP:0002024 | Gastrointestinal | Represents impaired intestinal nutrient absorption and may accompany symptomatic disease |
| Chronic diarrhea | HP:0002028 | Gastrointestinal | Common classical presentation but not required and may be absent in atypical disease |
| Abdominal pain | HP:0002027 | Gastrointestinal | Common but nonspecific gastrointestinal symptom |
| Abdominal distention | HP:0003270 | Gastrointestinal | May accompany gastrointestinal symptoms and malabsorption |
| Iron deficiency anemia | HP:0001891 | Hematological | Important extraintestinal manifestation and a clinically relevant trigger for testing |
| Decreased circulating folate concentration | HP:0100507 | Hematological | May occur with nutritional deficiency or impaired intestinal absorption |
| Thrombocytosis | HP:0001894 | Hematological | May occur as a reactive hematological abnormality |
| Elevated circulating hepatic transaminase concentration | HP:0002910 | Metabolic/hepatic | Can accompany CD and may prompt investigation when otherwise unexplained |
| Osteoporosis | HP:0000939 | Metabolic/skeletal | May occur in association with chronic malabsorption and altered bone health |
| Short stature | HP:0004322 | Metabolic/growth | Relevant pediatric growth phenotype that may prompt CD evaluation |
| Failure to thrive | HP:0001508 | Metabolic/growth | Relevant pediatric phenotype, particularly in younger children |
| Peripheral neuropathy | HP:0009830 | Neurological | Recognized extraintestinal neurological manifestation |
| Ataxia | HP:0001251 | Neurological | Uncommon neurological presentation associated with gluten-related disease |
| Cognitive impairment | HP:0100543 | Neurological | Nonspecific neurological feature reported in association with CD |
| Abnormal blistering of the skin | HP:0008066 | Dermatological | Characteristic blistering and vesiculobullous cutaneous manifestation of gluten-sensitive enteropathy; include when clinically documented |
| Infertility | HP:0000789 | Reproductive | Reported reproductive manifestation, particularly relevant in adult presentations |
| Recurrent spontaneous abortion | HP:0200067 | Reproductive | Reported reproductive manifestation associated with untreated CD |
Artificial intelligence in CD histopathology
The Marsh–Oberhuber system and its limitations
The Marsh–Oberhuber classification remains widely used for describing histological changes associated with CD, grading mucosal abnormalities from increased intraepithelial lymphocytes through villous atrophy.11 However, interobserver variability can affect intermediate-grade assessment. Quantitative digital pathology offers an approach beyond subjective grading. Gruver et al.35 developed pathologist-trained machine-learning classifiers that quantify villous and crypt measurements and intraepithelial lymphocyte features, generating quantitative scores that correlate with modified Marsh assessment and dietary intervention response. Griffin et al.36 reported the feasibility of automated cell-type and tissue-region classification across multisite whole-slide images (WSIs).
Deep learning models for CD histopathology and potential CD
Jaeckle et al.13 evaluated a machine-learning model for binary detection of CD from hematoxylin and eosin-stained duodenal biopsy WSIs. The model was trained using 3,383 WSIs from four hospitals and evaluated on 644 previously unseen scans from a geographically separate regional National Health Service (NHS) Trust, constituting an independent test dataset. The model achieved accuracy, sensitivity, and specificity exceeding 95%, with an area under the receiver operating characteristic curve (AUC) exceeding 0.99. Model–pathologist agreement was comparable with interpathologist agreement.13 Importantly, this study addressed binary CD detection rather than automated Marsh grading.
A particularly challenging diagnostic scenario is potential CD, characterized by persistently positive CD-specific serology in the presence of normal or non-atrophic small-intestinal mucosa; precise definition and management depend on guideline and clinical context.37 Piccialli et al.37 applied ML to predict biopsy-positive conversion in this population, identifying clinical and laboratory features associated with progression. Such work suggests a possible future role for AI-assisted risk stratification, but it does not establish a validated deep-learning method for continuous mucosal risk assessment.
Explainability remains important for clinical adoption. Attention-map visualization and class activation mapping may help highlight image regions contributing to a model’s prediction, potentially supporting review by a pathologist. In a clinical workflow, such visualizations should be treated as aids to interpretation rather than proof of biological causality.
Carreras reported convolutional neural network (CNN)-based CD classification with accuracy of 99.7%, precision of 99.6%, recall of 99.3%, and an F1 score of 99.5% on duodenal biopsy images, although single-center validation limits generalizability.38 Koh et al.39 and Stoleru et al.40 reported the feasibility of automated biopsy and capsule-endoscopy interpretation, respectively. Soylu and Bozkir applied EfficientNet to IgA-class endomysial antibody immunofluorescence pattern recognition.41 A systematic review by Hartmann Tolić et al.42 summarized promising performance across AI-assisted CD image-analysis tasks, while emphasizing heterogeneity and the need for broader validation. Key studies across these applications are summarized in Table 2.12,13,34-41
| Study | Algorithm | Dataset/task | Key result and interpretation |
|---|
| Jaeckle et al. (2025)13 | ML ensemble | 3,383 training/cross-validation WSIs; 644 WSI independent test set from geographically separate NHS Trust. Binary CD vs non-CD detection | Accuracy, sensitivity and specificity > 95%; AUC > 0.99. Promising independent-test performance; not automated Marsh grading |
| Dreyfuss et al. (2024)12 | XGBoost and comparators | Development: 677 highly seropositive cases/176,293 controls; test: 153 cases/41,087 controls; Maccabi, Israel | XGBoost AUC 0.86. Adolescents/adults; seropositivity/autoimmunity outcome, not biopsy-confirmed CD |
| Gruver et al. (2023)35 | Pathologist-trained ML | Multisite biopsy data/WSIs; quantitative feature assessment | Scores correlated with modified Marsh assessment and dietary response |
| Griffin et al. (2024)36 | Automated DL | Multisite WSIs; cell-type and tissue-region classification | Feasibility demonstrated; broader external validation remains needed |
| Piccialli et al. (2021)37 | ML | Potential CD cohort; biopsy-positive conversion prediction | Potential risk-stratification role; not a validated DL mucosal-risk model |
| Carreras (2024)38 | CNN | Single-center duodenal biopsy images; CD vs non-CD | Accuracy 99.7%; broader external validation needed |
| Koh et al. (2021)39 | CNN/ML | Duodenal biopsy images; automated interpretation | Supports feasibility of automated biopsy-image interpretation |
| Stoleru et al. (2022)40 | ML | Capsule endoscopy images; CD detection | Accuracy 94.1%; F1 score 94% |
| Soylu and Bozkir (2025)41 | EfficientNet | EMA immunofluorescence images | Emerging AI application to serological image interpretation |
Mao et al./ PhenoBrain (2025)34 | HPO-based AI | 2,271 cases, 431 rare diseases | Non-celiac rare-disease cohort; indirect methodological support only |
Machine learning for pre-diagnostic screening
Pre-diagnostic signatures in routine laboratory data
A pre-diagnostic phase may precede CD recognition, during which laboratory abnormalities such as iron-deficiency anemia and elevated transaminases can occur.9,10 Dreyfuss et al.12 developed five ML models using deidentified electronic medical record data. The model-development dataset included 677 highly seropositive cases and 176,293 controls; a distinct test dataset included 153 highly seropositive cases and 41,087 controls. The study population comprised adolescents and adults, and the outcome was identification of previously unrecognised CD autoimmunity or seropositivity rather than biopsy-confirmed CD. XGBoost achieved an AUC of 0.86, and performance remained above 0.80 at a four-year pre-index interval in the reported analyses. Iron-deficiency anemia, transaminitis, and low high-density lipoprotein levels were among the leading features identified by model interpretation.12
Regarding clinical actionability, an ML model could conceptually be used to prioritize patients for established CD serological testing. Any model-generated threshold would require local calibration and should not be treated as a validated universal cut-off. Once CD-specific serology is obtained, subsequent diagnostic evaluation should follow current age-appropriate clinical guidelines rather than a fixed threshold embedded in this framework. HLA-DQ2/DQ8 testing should be considered an adjunctive test in selected situations of diagnostic uncertainty, not a routine confirmatory step. Any proposed monitoring interval should likewise be regarded as conceptual and unvalidated.
B-cell receptor repertoire and multiomics approaches
Beyond conventional laboratory data, emerging ML approaches exploit deeper biological signals. Shemesh et al.43 applied ML to naïve B-cell receptor repertoire sequencing data, showing that B-cell receptor features distinguished patients with CD from controls in the study dataset. Abraham et al.44 reported genomic prediction of CD risk using statistical learning, and Piroozkhah et al.45 reviewed emerging multiomics approaches that combine genomics, proteomics, metabolomics, microbiome data, and AI. These approaches remain primarily investigational and require validation before they can be incorporated into routine screening.
A four-layer integration framework
Framework architecture and clinical feasibility
The convergence of EHR phenotyping, HPO-based semantic encoding, ML pre-screening, and DL-assisted histopathology can be conceptualized as a four-layer implementation framework for AI-assisted CD diagnostics (Fig. 1).13 The layers are conceptualized as sequential, with outputs potentially informing the next layer. The framework is not proposed as a replacement for current practice but as a modular architecture that could be implemented incrementally according to institutional infrastructure and validation status.
Layers 1 and 2 may be the most readily deployable components in settings with structured EHR systems. Layer 3 could support targeted serological testing where appropriate data are available, whereas Layer 4 requires digitized WSI infrastructure. The principal implementation challenges include data interoperability, representative training and validation datasets, clinical workflow integration, clinician oversight, and appropriate governance. Claims of clinical readiness should remain proportionate to the available evidence.
Layer 1 - EHR phenotyping: Structured and unstructured EHR data can be processed using NLP and rule-based algorithms to extract computable phenotypes. ICD codes can be augmented with laboratory values, medication records, and pathology reports to construct multimodal patient representations.19-21,26
Layer 2 - HPO semantic encoding: Patient phenotype profiles can be mapped to HPO terms using automated NLP pipelines or structured data extraction.32,33,46 The resulting HPO-encoded profiles can support similarity-based matching, interoperability, and research data integration. This layer draws on the curated minimum viable CD phenome in Table 1, which is a conceptual feature set requiring formal validation.
Layer 3 - ML pre-screening: HPO-encoded phenotype profiles, augmented with longitudinal laboratory data, can be processed by ML models to generate a CD risk score.12,25,37 In a validated implementation, model outputs could help prioritize patients for serological testing. Clinician review of model-derived feature contributions may provide additional transparency.
Layer 4 - DL histopathology: For patients undergoing endoscopic biopsy, digitized WSIs can be processed by DL models for automated CD detection or quantitative histological analysis.13,35,38-40 The Jaeckle et al.13 model performs binary CD detection, and its AUC > 0.99 was obtained on an independent test set of 644 WSIs from a geographically separate NHS Trust. It should not be characterized as an automated Marsh-grading system. Other studies have explored quantitative Marsh-related features, but these approaches should be distinguished from binary CD detection. Discordant or low-confidence cases should remain under expert pathologist review.
HPO-to-phecode-to-ICD translation
A potential enabling component of the framework is cross-mapping among HPO terms, phecodes, and ICD codes. McArthur et al.32 developed a systematic mapping between HPO terms and phecodes, enabling phenotype data generated in research contexts to be related to clinical coding systems. The Open Annotations for Rare Diseases resource provides real-world-data-derived HPO and rare-disease phenotype annotations extracted from EHRs using ontology mapping and NLP.47 This rare-disease resource is not a CD-specific validation resource, but it illustrates how real-world EHR data can complement expert-curated phenotype annotations.
Limitations
Data quality and standardization
The performance of AI models for CD diagnosis is constrained by the quality, representativeness, and standardization of training data. Histopathological datasets vary in staining protocols, scanner hardware, sampling, and annotation practices, limiting cross-site generalizability.42 Federated learning and harmonized data standards may facilitate broader validation. The AUC > 0.99 reported by Jaeckle et al.13 was obtained on an independent test set from a geographically separate NHS Trust, which is an important strength; nevertheless, prospective, multicenter studies assessing clinical utility and validation in broader populations remain necessary. Binary CD detection should not be conflated with automated Marsh grading.
Algorithmic bias and health equity
AI models trained predominantly on data from high-income populations may perform differently in underrepresented groups. CD has distinct epidemiological patterns and HLA distributions across populations,30,31 and AI tools should therefore be evaluated in diverse cohorts before deployment. The current evidence base remains concentrated in selected academic and healthcare systems, underscoring the need for inclusive dataset curation and prospective bias assessment.
Explainability and clinical trust
The adoption of AI in clinical practice requires interpretable and clinically reviewable outputs. Shapley-value-based explanations and attention-map visualizations may assist clinicians in understanding model behavior and identifying potential errors.12,13 Regulatory requirements vary by jurisdiction and intended use; therefore, developers should address applicable regulatory, validation, documentation, monitoring, and post-market requirements rather than assuming a single universal pathway.
Age-specific validation gaps
As noted in the discussion of HPO-based phenotyping and age-specific considerations, the principal ML pre-screening evidence currently includes adolescents and adults rather than exclusively pediatric cohorts.12 The lack of prospective validation in dedicated pediatric populations is a critical limitation. Pediatric CD registries capable of providing adequately powered, age-stratified datasets are therefore needed.
Future directions
Emerging technologies
Foundation models pretrained on large pathology image corpora may improve generalization for histopathological tasks when CD-specific training data are limited, but their performance in CD remains to be established prospectively. Large language models integrated with HPO-based pipelines may further automate phenotype extraction from clinical notes.33,46 Multiomics integration, combining genomics, proteomics, microbiome data, and EHR phenotypes, represents a potential future layer for precision medicine.45
Research priorities
Future research should prioritize: (1) construction of large, diverse, multicenter CD histopathology and EHR datasets that include underrepresented populations; (2) prospective real-world validation of ML pre-screening algorithms in primary care and pediatric-specific cohorts; (3) international harmonisation of CD-relevant HPO annotation sets with explicit age-stratified feature weighting; (4) development and external validation of DL models for subtle mucosal changes; (5) prospective bias auditing and inclusive AI development; and (6) implementation and regulatory science studies to support safe translation of AI-based CD tools into clinical workflows.
Conclusions
Celiac disease presents a compelling setting for the application of AI and phenomics to a common, underdiagnosed, and phenotypically heterogeneous condition. The convergence of computable EHR phenotyping, HPO-based semantic encoding, ML pre-screening, and DL-assisted histopathology provides a coherent four-layer framework for future precision diagnostics. Selected models have demonstrated promising CD detection performance, including performance comparable with that of pathologists in specific datasets, while HPO-based AI has shown promise in non-celiac rare-disease differential diagnosis. These findings support further investigation but do not establish clinical superiority or routine clinical utility in CD.
Realizing the potential of this framework requires attention to dataset diversity and standardization, algorithmic bias, age-specific validation, clinical integration, and governance. The four-layer architecture is modular and may be implemented incrementally as individual components achieve adequate validation. Prospective, multicenter studies with external validation are needed before routine clinical deployment can be recommended.
Declarations
Acknowledgments
The authors thank the open-access scientific community for making the literature and ontological resources cited in this review freely available. In accordance with the journal’s policy on AI tool disclosure, the authors declare the use of BIOMNI as an AI-assisted tool to support literature synthesis and language editing during manuscript preparation. Figure 1 was designed by the authors specifically for this manuscript; its visual elements, including icons and schematic components, were created by the authors rather than copied or adapted from third-party copyrighted sources. The figure was subsequently refined with the assistance of artificial intelligence during manuscript preparation. The authors remain responsible for verifying the scientific accuracy and licensing status of all visual components.
Funding
This work received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. No grant numbers apply, as no funding was received in support of this research.
Conflict of interest
The authors declare no conflict of interest related to this publication.
Author contributions
Conceptualization (HR, NB, IB); literature review (HR, NB); writing—original draft (HR); writing—review and editing (HR, NB, IB); informatics and technical review (IB); figure and table design (HR). All authors contributed to interpretation, approved the final manuscript, and agree to be accountable for all aspects of the work.