Back to Authors

Gubatan, John

https://orcid.org/0000-0001-6037-2883

Mayo Clinic

Mayo Clinic

Citations 6

Single-cell immune repertoire atlas maps coordinated circulating adaptive immune states in inflammatory bowel disease
Gubatan, John , Ye, Jiayu, Lund‐Andersen, Christin , Canas, Jorge, Zhou, Yuxi , Boye, Theresa, Hoang, Jacqueline, Sojwal, Raoul, Fardeen, Touran , Tran, Tracy, Sokzini, Peacha, Koh, Ally H, Gibson, Marielle, Huang, Yuting , Peterson, Kathryn , Sinha, Sidhartha, Rogalla, S, Nielsen, Ole Haagen , Rosen, Michael J, Zhou, Yuxi , Canas, Jorge , Huang, Yuting , Ye, Jiayu, Boye, Theresa Louise , Hoang, Jacqueline, Sojwal, Raoul S. , Tracy, Tran, Rosen, Michael J., Sandve, G K , Gubatan, John Mark

No container title available | September 23, 2026 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.64898/2026.09.20.753043 Get It

Single-cell studies have defined immune states in inflammatory bowel disease (IBD), but how adaptive receptor histories organize circulating immunity remains unclear. We generated a single-cell transcriptomic atlas of peripheral blood from 249 participants with Crohn's disease, ulcerative colitis, or non-IBD control status, including 182 with productive TCR and BCR recovery. Expanded TCR clonotypes marked inflammatory-memory and cytotoxic states, while distinct but similar paired TCRs shared inflammatory programs across participants. BCR lineage maturation linked IgA-associated mucosal and plasma B cell programs to somatic mutation and class switching, distinguishing maturation-associated biology from clonal expansion. Helper, regulatory, and cytotoxic T-cell programs covaried with B-cell states, and inferred interactions nominated reciprocal antigen-presentation and helper pathways. Repertoire-based machine learning distinguished diagnosis, inflammation, and contemporaneous six-month treatment-response status. Together, this atlas connects receptor architecture to coordinated systemic immune remodeling, establishes a foundation for repertoire-informed patient stratification, and prioritizes candidate mechanisms of IBD pathogenesis.
Single-cell profiling reveals stromal niche and α5β1 regulation of neutrophil phenotypes in ulcerative colitis
Zhou, Yuxi , Canas, Jorge , Ye, Jiayu , Hoang, Jacqueline, Ferkel-Soltau, Sonia, Holman, Derek, Sojwal, Raoul, Fardeen, Touran , Vazques, Gustavo, Huang, Yuting , Peterson, Kathryn , Nolan, Garry, Sinha, Sidhartha, Rogalla, S, Gubatan, John , Huang, Yuting , Hoang, Jacqueline, Sojwal, Raoul S. , Nolan, Garry P. , Gubatan, John Mark

No container title available | September 23, 2026 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.64898/2026.09.21.753318 Get It

Abstract Background and Aims Fibroblast–neutrophil interactions are implicated in ulcerative colitis (UC) treatment resistance, but whether fibroblasts modify neutrophil inflammatory activity remains unclear. We investigated responses to fibroblast exposure, stromal depletion, and α5β1-directed treatment. Methods We integrated human and mouse colonic single-cell RNA sequencing (scRNA-seq), spatial protein profiling (CODEX), fibroblast–neutrophil coculture, and dextran sulfate sodium (DSS) mouse colitis with fibroblast activation protein (FAP) ablation or α5β1-directed treatment. Readouts included cellular representation and protein fluorescence via flow cytometry, neutrophil extracellular trap (NET)-associated elastase activity, RNA programs and epithelial responses using scRNA-seq. Results Inflamed UC tissue contained higher fractions of inflammatory fibroblasts and oncostatin M (OSM)- and CXCR4-associated neutrophils. Inflammatory fibroblasts had higher FAP–α5β1 module scores, and proximity to FAP-high fibroblasts was associated with higher neutrophil OSM and CXCR4 fluorescence. UC fibroblasts increased neutrophil OSM, CXCR4 and myeloperoxidase (MPO) fluorescence and NET-associated elastase activity relative to control fibroblasts, while α5β1-directed treatment attenuated these responses. Both mouse FAP ablation and α5β1 blockade in DSS colitis reduced histologic inflammation and overall neutrophil frequency. Ablation broadly reduced recovered neutrophil representation, whereas blockade increased OSM-positive and PADI4-positive percentages within neutrophils despite lower MPO fluorescence within these subsets. Both interventions were associated with lower epithelial chemokine scores but distinct absorptive and mucus/secretory responses. Conclusions Fibroblast exposure modifies neutrophil inflammatory phenotype and effector-associated activity, identifying candidate therapeutic pathways in UC. Evaluating these responses alongside neutrophil representation and tissue inflammation could inform pharmacodynamic assessment of stromal-directed therapies. Graphical Abstract What You Need to Know Background and Context Fibroblast–neutrophil interactions are associated with treatment resistance in ulcerative colitis. Whether fibroblast exposure modifies neutrophil inflammatory phenotype and effector activity, beyond its established recruitment role, remains unclear. New Findings UC fibroblasts increased neutrophil inflammatory proteins and NET-associated elastase activity; α5β1-directed treatment attenuated these responses. FAP ablation and α5β1-directed treatment reduced colitis and overall neutrophil frequency. α5β1-directed treatment lowered MPO within OSM-positive and PADI4-positive neutrophils despite higher marker-positive percentages. Epithelial single-cell RNA sequencing revealed reduced chemokine programs with both interventions but distinct barrier-associated and regenerative expression patterns. Limitations Cellular drug targets remain unresolved. Spatial associations and elastase activity do not establish direct contact, NET formation, or clinical efficacy. Clinical Research Relevance Fibroblast-associated pathways offer candidate therapeutic opportunities to modify neutrophil inflammatory phenotype and effector activity in UC. Coculture responses were pharmacologically modifiable, and stromal interventions reduced experimental colitis. Measuring neutrophil frequency, subset proteins, and effector activity alongside tissue inflammation could guide pharmacodynamic assessment. Clinical efficacy and treatment-selection utility require validation. Basic Research Relevance Human coculture provides experimental evidence that fibroblast exposure modifies neutrophil inflammatory proteins and NET-associated elastase activity. FAP+ stromal ablation and α5β1-directed treatment produced distinct patterns of neutrophil representation and subset proteins in colitis. Epithelial single-cell RNA sequencing revealed reduced chemokine programs with both interventions but differing barrier-associated and regenerative expression patterns. Together, these findings motivate mechanistic studies of how stromal recruitment cues, soluble signals, and matrix-associated interactions regulate neutrophil responses, and whether these changes contribute to reduced epithelial inflammation and tissue repair. Lay Summary Fibroblasts, support cells in the intestine, increased inflammatory activity in immune cells called neutrophils. Blocking a pathway involved in cell attachment reduced these responses in cultures. Targeting fibroblast-associated pathways also reduced experimental colitis, with changes in intestinal lining gene activity related to inflammation, barrier function, and repair. These findings suggest potential approaches to treating ulcerative colitis.
Unsupervised Machine Learning for Adaptive Immune Receptors with immuneML
Pavlović, Milena , Würtzen, Charlotte , Kanduri, Chakravarthi , Mamica, Maria , Scheffer, Lonneke , Lund‐Andersen, Christin , Gubatan, John Mark , Ullmann, Theresa , Greiff, Victor , Sandve, G K , Gubatan, John

No container title available | April 18, 2026 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.64898/2026.04.15.718648 Get It

Abstract Machine learning (ML) enables adaptive immune receptor repertoires (AIRRs) analyses for biomarker identification and therapeutic development. With the majority of AIRR data partially or imperfectly labeled, unsupervised ML is essential for motif discovery, biologically meaningful clustering, and generation of novel receptor sequences. However, no unified framework for unsupervised ML exists in the AIRR field, hindering the assessment of model robustness and generalizability. Here, we present an immuneML release advancing unsupervised ML in the AIRR field through unified clustering workflows, interpretable generative modeling, integration with protein language model embeddings, dimensionality reduction, and visualization. We demonstrate immuneML’s utility in three use cases: (i) benchmarking generative models for epitope-specific sequence generation, assessing specificity and novelty, (ii) systematic evaluation of clustering approaches on experimental receptor sequences against biological properties, such as epitope specificity and MHC, and (iii) unsupervised analysis of an experimental AIRR dataset to examine potential confounding, a practice widespread in related fields but unexplored in AIRR analyses.
P0009 Machine Learning Enables Deep Immunophenotyping of Inflammatory Bowel Disease Using Single-Cell Transcriptomic and Immune Repertoire Signatures in Peripheral Blood Mononuclear Cells
Gubatan, J M, Jiayu, Y, Christin, L A, Hoang, Jacqueline, Sojwal, Raoul S. , Fardeen, Touran , Tran, T, Sokzini, P, Koh, Ai Leen , Gibson, M, Sinha, S, Rogalla, S, Rosen, M, Nielsen, O H, Nielsen, Ole Haagen , Gubatan, John , Ye, Jiayu, Lund‐Andersen, Christin , Boye, Theresa Louise , Sokzini, Peacha, Koh, A, Gibson, Michelle , Sinha, Sidhartha R. , Sandve, G K , Gubatan, John Mark

Journal of Crohn’s and Colitis | January 21, 2026 | Volume: 20 | Issue: Supplement_1 | Pages: No page available

10.1093/ecco-jcc/jjaf231.190 Get It

Abstract Background Inflammatory bowel disease (IBD), including ulcerative colitis (UC) and Crohn’s disease (CD), is characterized by immune dysregulation1,2. Blood-based biomarkers in diagnosing patients with IBD are lacking. We tested the hypothesis that circulating immune cell types and immune repertoire could distinguish patients with IBD from controls. Methods We performed deep immunophenotyping with scRNA-seq with T-cell receptor (TCR) and B-cell receptor (BCR) sequencing on peripheral blood mononuclear cells (PBMCs) from patients with IBD and healthy controls (87 UC, 127 CD, 35 controls). We applied machine learning to quantify cell abundance (miloR3), construct diagnostic classifiers (immuneML4), and infer cell-cell ligand-receptor interactions (MultiNicheNet5) between patients with IBD and controls. Results A scRNA-seq atlas of over 1.1 million PBMCs consisting of 53 distinct cell types was generated (Fig1A). In classifying UC vs controls, support vector machine model had a balanced accuracy of 0.90, consisting of increased IgA and IgG plasma B cells, switched memory B cells, TRegs, and decreased CD16 monocytes, cDC1, non-switched B cells, and MAIT (Fig.1B). In classifying CD vs controls, K-nearest neighbors model had a balanced accuracy of 0.87, consisting of increased IgA plasma B cells, non-switched B cells, CD4 Th22, TReg cytotoxic, CD4 Th1/Th17, and pDC and decreased CD16 monocytes, gdT, CD4 Th17, CD4 Th2, and MAIT cells.Fig.1C). Inflammation was associated with increased IgG plasma B in patients with both UC and CD. Using 4-mer amino acid motif of the BCR complementary-determining region 3 (CDR3) logistic regression model, immuneML had a balanced accuracy of 0.82 (0.75-0.92) in classifying CD vs controls (Fig.1D with models showing top 25 BCR predictive features). Using 4-mer amino acid motif of the BCR CDR3 logistic regression model, immuneML had a balanced accuracy of 0.90 (0.70-1.00) in classifying UC vs controls (Fig.1E with models showing top 25 BCR predictive features). Using 3-mer amino acid motif of the TCR CDR3 logistic regression model, immuneML had a balanced accuracy of 0.76 (0.72-081) in classifying CD vs UC. (Fig.1F with models showing top 25 TCR predictive features). In UC vs controls, MultiNicheNet demonstrated increased signaling in ADAM10-NOTCH2 between TRegs and naïve B cells and IL15-IL2RA between CD14 monocytes and CD4 Th17 cells. In CD vs controls, there was increased signaling in IL15-IL2RA between CD14 monocytes and CD4 Th17 cells and ADAM10-TSPAN5 between CD14 monocytes and TRegs (Fig.1G). Conclusion Applying machine learning models to scRNA-seq and immune repertoire datasets enable blood-based deep immunophenotyping with diagnostic potential and reveals novel insights into cell signaling in IBD. References: 1.Kong L, Pokatayev V, Lefkovith A, Carter GT, Creasey EA, Krishna C, Subramanian S, Kochar B, Ashenberg O, Lau H, Ananthakrishnan AN, Graham DB, Deguine J, Xavier RJ. The landscape of immune dysregulation in Crohn’s disease revealed through single-cell transcriptomic profiling in the ileum and colon. Immunity. 2023 Feb 14;56(2):444-458.e5. doi: 10.1016/j.immuni.2023.01.002. Epub 2023 Jan 30. Erratum in: Immunity. 2023 Dec 12;56(12):2855. doi: 10.1016/j.immuni.2023.10.017. PMID: 36720220; PMCID: PMC9957882. 2.Smillie CS, Biton M, Ordovas-Montanes J, Sullivan KM, Burgin G, Graham DB, Herbst RH, Rogel N, Slyper M, Waldman J, Sud M, Andrews E, Velonias G, Haber AL, Jagadeesh K, Vickovic S, Yao J, Stevens C, Dionne D, Nguyen LT, Villani AC, Hofree M, Creasey EA, Huang H, Rozenblatt-Rosen O, Garber JJ, Khalili H, Desch AN, Daly MJ, Ananthakrishnan AN, Shalek AK, Xavier RJ, Regev A. Intra- and Inter-cellular Rewiring of the Human Colon during Ulcerative Colitis. Cell. 2019 Jul 25;178(3):714-730.e22. doi: 10.1016/j.cell.2019.06.029. PMID: 31348891; PMCID: PMC6662628. 3.Dann E, Henderson NC, Teichmann SA, Morgan MD, Marioni JC. Differential abundance testing on single-cell data using k-nearest neighbor graphs. Nat Biotechnol. 2022 Feb;40(2):245-253. doi: 10.1038/s41587-021-01033-z. Epub 2021 Sep 30. PMID: 34594043; PMCID: PMC7617075. 4.Pavlović M, Scheffer L, Motwani K, Kanduri C, Kompova R, Vazov N, Waagan K, Bernal FLM, Costa AA, Corrie B, Akbar R, Al Hajj GS, Balaban G, Brusko TM, Chernigovskaya M, Christley S, Cowell LG, Frank R, Grytten I, Gundersen S, Haff IH, Hovig E, Hsieh PH, Klambauer G, Kuijjer ML, Lund-Andersen C, Martini A, Minotto T, Pensar J, Rand K, Riccardi E, Robert PA, Rocha A, Slabodkin A, Snapkov I, Sollid LM, Titov D, Weber CR, Widrich M, Yaari G, Greiff V, Sandve GK. The immuneML ecosystem for machine learning analysis of adaptive immune receptor repertoires. Nat Mach Intell. 2021 Nov;3(11):936-944. doi: 10.1038/s42256-021-00413-z. Epub 2021 Nov 16. PMID: 37396030; PMCID: PMC10312379. 5.Browaeys, R., Gilis, J., Sang-Aram, C., De Bleser, P., Hoste, L., Tavernier, S., Lambrechts, D., Seurinck, R. and Saeys, Y., 2023. MultiNicheNet: a flexible framework for differential cell-cell communication analysis from multi-sample multi-condition single-cell transcriptomics data. BioRxiv, pp.2023-06. Conflict of interest: Dr. Gubatan, John Mark: No conflict of interest Jiayu, Ye: No conflict of interest Christin, Lund-Andersen: No conflict of interest Boye, Theresa: No conflict of interest Hoang, Jacqueline: No conflict of interest Sojwal, Raoul: No conflict of interest Fardeen, Touran: No conflict of interest Tran, Tracy: No conflict of interest Sokzini, Peacha: No conflict of interest Koh, Ally: No conflict of interest Gibson, Marielle: No conflict of interest Sinha, Sidhartha: No conflict of interest Rogalla, Stephan: No conflict of interest Rosen, Michael: No conflict of interest Nielsen, Ole Haagen: No conflict of interest Sandve, Geir Kjetil: No conflict of interest
Long-term follow-up of gastrointestinal CAR T-cell lymphoma: homing, clonal expansion, and response to cyclosporine
Hosoya, Hitomi , Torres, Armando N. Bastidas , Fernandez‐Pol, Sebastian , Gubatan, John , Durán, George E. , Lohman, Claire , Wright, Charles V., Sahaf, Bita , Mackall, Crystal L. , Sidana, Surbhi , Khodadoust, Michael S. , Mikkilineni, Lekha, Duran, George E , Lohman, Claire , Mackall, Crystal , Miklos, David B , Sidana, Surbhi , Kurtz, David M , Khodadoust, Michael S , Hosoya, Hitomi , Najidh, Safa , Dong, Fei , Ehlinger, Zachary J. , Wright, Charles, Miklos, David B. , Kurtz, David M. , Mikkilineni, Lekha, Gubatan, John Mark

Blood | November 25, 2025 | Volume: 147 | Issue: 11 | Pages: 1191-1198

10.1182/blood.2025031423 Get It

Abstract Chimeric antigen receptor (CAR) T-cell therapy has emerged as a transformative treatment for hematological malignancies, yet its potential to drive lymphomagenesis poses significant clinical concerns. In this study, we investigated the mechanisms underlying CAR T-cell–associated lymphomagenesis in the gastrointestinal (GI) tract on a single case, focusing specifically on the role of integrin α4β7 expression and a predisposing somatic SH2B3 mutation. We observed oligoclonal CAR T cells homing to, and clonally expanding in, the GI tract, with the dominant expanded clone harboring both a pathogenic SH2B3 mutation and a CAR transgene integration within a TFCP2 locus. The clonal CAR T cells subsequently transitioned beyond the GI tract into the peripheral blood, suggesting a potential pathway for systemic dissemination. We found clinical, histological, and molecular evidence demonstrating the efficacy of cyclosporine in reducing the expanded malignant clone and achieving durable clinical remission for more than a year. Our findings highlight the complex interplay between CAR T-cell therapy, preexisting genetic vulnerabilities, and the GI microenvironment, emphasizing the need for vigilant monitoring and tailored therapeutic strategies to address the risks associated with CAR T-cell lymphomagenesis.
Immune Effector Cell-Associated Enterocolitis following Chimeric Antigen Receptor T-cell Therapy in Multiple Myeloma
Song, Jinming , Lekakis, Lazaros J. , Samraj, Annie N. , Naresh, Kikkeri N. , Vazquez‐Martinez, Mariola , Mikkilineni, Lekha , Corraes, Andre De Menezes Silva , Kalariya, Nilesh, Patel, Krina K. , Ferreri, Christopher J. , Voorhees, Peter M. , Valdes, Cesar Rodriguez , Cowan, Andrew J., Sborov, Douglas W. , Lin, Yi, Naresh, Kikkeri N , Fortuna, Gliceida Galarza , Banerjee, Rahul , Fernandez‐Pol, Sebastian , Baz, Rachid , Spiegel, Jay Y. , Gubatan, John , Sidana, Surbhi , Shim, Kevin G. , Fonseca, Rafaël , Richard, Shambavi , Hansen, Doris K. , Wang, Yinghong

No container title available | August 10, 2024 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.21203/rs.3.rs-4727434/v1 Get It

Abstract Background: Several immune effector cell (IEC) therapies have been approved for multiple myeloma (MM), including chimeric antigen receptor T-cell (CAR-T) therapies such as idecabtagene vicleucel (ide-cel) and ciltacabtagene autoleucel (cilta-cel). We report 14 cases of IEC-associated enterocolitis, in five cases fatal, following CAR-T therapy in MM. Methods: We performed a retrospective analysis of cases of diarrhea, enterocolitis, or colitis following commercial ide-cel or cilta-cel (including out-of-specification products) at 10 centers within the United States Myeloma Immunotherapy Consortium. Patients were treated according to the local standard of care at each participating center. Results: We identified 14 cases of IEC-associated enterocolitis out of 1159 CAR-T infusions, corresponding to a 1.2% incidence (0.2% for ide-cel and 2.3% for cilta-cel). This proportion excludes 3 cases of potential or confirmed T-cell lymphoproliferative disorders of the bowel. Patients developed acute-onset symptoms (typically non-bloody Grade 3+ diarrhea) with negative infectious workup beginning a median of 79 days (range: 22-210 days) after CAR-T therapy and a median of 74 days after cytokine release syndrome resolution (range: 2-202 days). Biopsies uniformly demonstrated inflammation, with specific findings including intra-epithelial lymphocytosis, crypt dropout on colon biopsies, and villous blunting on duodenal biopsies. Systemic corticosteroids were initiated in 10 patients (71%) a median of 26 days (range: 1-80 days) following symptom onset. Only 40% of steroid-treated patients had symptomatic improvement; subsequent infliximab or vedolizumab led to improvement in 50% and 33% of corticosteroid-refractory patients, respectively. With a median follow-up of 224 days, n = 4 have had diarrhea resolution, n = 5 have had improved or stable symptoms, and n = 5 (36%) have died from IEC-associated enterocolitis or treatment-emergent infections. Conclusions: IEC-associated enterocolitis is a distinct but rare complication of CAR-T therapy characterized by non-bloody diarrhea typically beginning 1-3 months after infusion. Thorough diagnostic workup (including evaluation for potential T-cell malignancies) and appropriate management with an expert gastroenterologist are essential. Further research is critical to understand the pathogenesis of this novel toxicity. In the interim, the early use of infliximab or vedolizumab may potentially hasten symptom resolution and lower reliance on high-dose corticosteroids during the post-CAR-T period.