{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,25]],"date-time":"2026-09-25T17:59:17Z","timestamp":1790359157309,"version":"4.1.0"},"reference-count":34,"publisher":"Oxford University Press (OUP)","issue":"9","license":[{"start":{"date-parts":[[2026,9,2]],"date-time":"2026-09-02T00:00:00Z","timestamp":1788307200000},"content-version":"vor","delay-in-days":2,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000054","name":"National Cancer Institute","doi-asserted-by":"publisher","award":["P30 CA016672"],"award-info":[{"award-number":["P30 CA016672"]}],"id":[{"id":"10.13039\/100000054","id-type":"DOI","asserted-by":"publisher"}]},{"name":"Cancer Prevention and Research Institute of Texas","award":["Rp180140"],"award-info":[{"award-number":["Rp180140"]}],"id":[{"id":"https:\/\/ror.org\/0003xa228","id-type":"ROR","asserted-by":"crossref"}]},{"name":"Andrew Sabin Family Fellowship"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2026,8,31]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>Examination of T cell receptor (TCR) clonality has become a way of understanding immunologic response to cancer and its interventions in recent years. An aspect of these analyses is determining which receptors expand or contract statistically significantly as a function of an exogenous perturbation such as therapeutic intervention.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>We characterize the commonly used Fisher\u2019s exact test approach for such analyses and propose an alternative formulation that does not necessitate pairwise, within-patient comparisons. We develop this flexible Bayesian longitudinal mixture model that accommodates variable length patient followup and handles missingness where present, not omitting data in estimation because of structural practicalities. Once clones are partitioned by the model into dynamic (expanding or contracting) and static categories, one can associate their counts or other characteristics with disease state, interventions, baseline biomarkers, and patient prognosis. We apply these developments to a cohort of prostate cancer patients who underwent randomized metastasis-directed therapy or not. Our analyses reveal a significant increase in clonal expansions among metastasis-directed therapy (MDT) patients and their association with later progressions both independent and within strata of MDT. Analysis of receptor motifs and VJ gene enrichment combinations using a high-dimensional penalized log-linear model we develop also suggests distinct biological characteristics of expanding clones, with and without inducement by MDT.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>An example model implementation in R\/STAN is available at doi.org\/10.5281\/zenodo.21209546<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btag659","type":"journal-article","created":{"date-parts":[[2026,9,2]],"date-time":"2026-09-02T12:47:25Z","timestamp":1788353245000},"source":"Crossref","is-referenced-by-count":0,"title":["Identifying expanding TCR clonotypes with a longitudinal Bayesian mixture model and their associations with cancer patient prognosis, metastasis-directed therapy, and VJ gene enrichment"],"prefix":"10.1093","volume":"42","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-3174-1656","authenticated-orcid":false,"given":"David","family":"Swanson","sequence":"first","affiliation":[{"name":"Department of Biostatistics, University of Texas MD Anderson Cancer Center, Houston, TX, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Alexander","family":"Sherry","sequence":"additional","affiliation":[{"name":"Department of Radiation Oncology, Mayo Clinic , Rochester, MN,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Cara","family":"Haymaker","sequence":"additional","affiliation":[{"name":"Department of Translational Molecular Pathology, University of Texas MD Anderson Cancer Center, Houston, TX, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-4510-0382","authenticated-orcid":false,"given":"Alexandre","family":"Reuben","sequence":"additional","affiliation":[{"name":"Departments of Thoracic, Head and Neck Medical Oncology, University of Texas MD Anderson Cancer Center, Houston, TX, United States"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Chad","family":"Tang","sequence":"additional","affiliation":[{"name":"Department of Radiation Oncology, University of Texas MD Anderson Cancer Center , Houston, TX,","place":["United States"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"286","published-online":{"date-parts":[[2026,9,2]]},"reference":[{"key":"2026092513142921300_btag659-B1","volume-title":"Categorical Data Analysis","author":"Agresti","year":"2013"},{"key":"2026092513142921300_btag659-B2","doi-asserted-by":"publisher","first-page":"e2207516120","DOI":"10.1073\/pnas.2207516120","article-title":"Inferring the T cell repertoire dynamics of healthy individuals","volume":"120","author":"Bensouda Koraichi","year":"2023","journal-title":"Proc Natl Acad Sci U S A"},{"key":"2026092513142921300_btag659-B3","first-page":"2","article-title":"Hamiltonian Monte Carlo for hierarchical models","volume":"79","author":"Betancourt","year":"2015","journal-title":"Curr Trends Bayesian Methodol Appl"},{"key":"2026092513142921300_btag659-B4","doi-asserted-by":"crossref","DOI":"10.1201\/9781003456285","volume-title":"Statistical Inference","author":"Casella","year":"2024"},{"key":"2026092513142921300_btag659-B5","doi-asserted-by":"crossref","first-page":"1816068","DOI":"10.3389\/fimmu.2026.1816068","article-title":"HLA alleles shape distinct biases in the usage preferences of tcr v \u03b2 segments","volume":"17","author":"Castorina","year":"2026","journal-title":"Front Immunol"},{"key":"2026092513142921300_btag659-B6","doi-asserted-by":"publisher","first-page":"e70164","DOI":"10.1002\/cam4.70164","article-title":"Human TCR repertoire in cancer","volume":"13","author":"Chen","year":"2024","journal-title":"Cancer Med"},{"key":"2026092513142921300_btag659-B7","doi-asserted-by":"publisher","first-page":"107023","DOI":"10.1016\/j.imlet.2025.107023","article-title":"Characteristics and specific differential gene analysis of the tcr immune repertoire in secondary adult HLH lymphocytes","volume":"275","author":"Chen","year":"2025","journal-title":"Immunol Lett"},{"key":"2026092513142921300_btag659-B8","doi-asserted-by":"publisher","first-page":"743448","DOI":"10.3389\/fimmu.2021.743448","article-title":"Longitudinal high-throughput sequencing of the T-cell receptor repertoire reveals dynamic change and prognostic significance of peripheral blood TCR diversity in metastatic colorectal cancer during chemotherapy","volume":"12","author":"Chen","year":"2021","journal-title":"Front Immunol"},{"key":"2026092513142921300_btag659-B9","doi-asserted-by":"publisher","first-page":"187","DOI":"10.1111\/j.2517-6161.1972.tb00899.x","article-title":"Regression models and life-tables","volume":"34","author":"Cox","year":"1972","journal-title":"J R Stat Soc Ser B (Methodological)"},{"key":"2026092513142921300_btag659-B10","doi-asserted-by":"publisher","first-page":"4517","DOI":"10.1128\/JVI.03474-14","article-title":"Dynamics of the cytotoxic T cell response to a model of acute viral infection","volume":"89","author":"DeWitt","year":"2015","journal-title":"J Virol"},{"key":"2026092513142921300_btag659-B11","doi-asserted-by":"publisher","first-page":"87","DOI":"10.2307\/2340521","article-title":"On the interpretation of x2 from contingency tables, and the calculation of P","volume":"85","author":"Fisher","year":"1922","journal-title":"J R Stat Soc"},{"key":"2026092513142921300_btag659-B12","volume-title":"Statistical Methods for Research Workers. Number 5","author":"Fisher","year":"1928"},{"key":"2026092513142921300_btag659-B13","doi-asserted-by":"publisher","first-page":"123","DOI":"10.1111\/j.1467-9868.2010.00765.x","article-title":"Riemann manifold Langevin and Hamiltonian Monte Carlo methods","volume":"73","author":"Girolami","year":"2011","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026092513142921300_btag659-B14","doi-asserted-by":"publisher","first-page":"26","DOI":"10.1002\/path.5396","article-title":"Characteristics, dynamic changes, and prognostic significance of TCR repertoire profiling in patients with renal cell carcinoma","volume":"251","author":"Guo","year":"2020","journal-title":"J Pathol"},{"key":"2026092513142921300_btag659-B15","first-page":"1593","article-title":"The No-U-Turn sampler: adaptively setting path lengths in Hamiltonian Monte Carlo","volume":"15","author":"Hoffman","year":"2014","journal-title":"J Mach Learn Res"},{"key":"2026092513142921300_btag659-B16","doi-asserted-by":"publisher","first-page":"e1008814","DOI":"10.1371\/journal.pcbi.1008814","article-title":"Predicting recognition between T cell receptors and epitopes with TCRGP","volume":"17","author":"Jokinen","year":"2021","journal-title":"PLoS Comput Biol"},{"key":"2026092513142921300_btag659-B17","doi-asserted-by":"publisher","first-page":"587014","DOI":"10.3389\/fimmu.2020.587014","article-title":"Characteristics of TCR repertoire associated with successful immune checkpoint therapy responses","volume":"11","author":"Kidman","year":"2020","journal-title":"Front Immunol"},{"key":"2026092513142921300_btag659-B18","doi-asserted-by":"publisher","first-page":"963","DOI":"10.2307\/2529876","article-title":"Random-effects models for longitudinal data","volume":"38","author":"Laird","year":"1982","journal-title":"Biometrics"},{"key":"2026092513142921300_btag659-B19","doi-asserted-by":"publisher","first-page":"627","DOI":"10.1080\/10618600.2014.938812","article-title":"Learning interactions via hierarchical group-lasso regularization","volume":"24","author":"Lim","year":"2015","journal-title":"J Comput Graph Stat"},{"key":"2026092513142921300_btag659-B20","doi-asserted-by":"publisher","first-page":"149","DOI":"10.1146\/annurev-clinpsy-080822-051727","article-title":"Missing data analysis","volume":"20","author":"Little","year":"2024","journal-title":"Annu Rev Clin Psychol"},{"key":"2026092513142921300_btag659-B21","volume-title":"Statistical Analysis with Missing Data","author":"Little","year":"2019"},{"key":"2026092513142921300_btag659-B22","doi-asserted-by":"publisher","first-page":"603","DOI":"10.1200\/JCO.2024.42.3_suppl.603","article-title":"Addition of metastasis-directed therapy to systemic therapy for oligometastatic pancreatic ductal adenocarcinoma (EXTEND): a multicenter, randomized phase II trial","volume":"42","author":"Ludmir","year":"2024","journal-title":"JCO"},{"key":"2026092513142921300_btag659-B23","doi-asserted-by":"publisher","first-page":"29544","DOI":"10.1038\/srep29544","article-title":"Analyzing the CDR3 repertoire with respect to TCR\u2014beta chain V-D-J and V-J rearrangements in peripheral T cells using HTS","volume":"6","author":"Ma","year":"2016","journal-title":"Sci Rep"},{"key":"2026092513142921300_btag659-B24","doi-asserted-by":"publisher","first-page":"417","DOI":"10.1111\/j.1467-9868.2010.00740.x","article-title":"Stability selection","volume":"72","author":"Meinshausen","year":"2010","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026092513142921300_btag659-B25","doi-asserted-by":"publisher","first-page":"eadx1751","DOI":"10.1126\/sciadv.adx1751","article-title":"Quantitative and large-scale investigation of human TCR-HLA cross-reactivity","volume":"11","author":"Pan","year":"2025","journal-title":"Sci Adv"},{"key":"2026092513142921300_btag659-B26","volume-title":"R: A Language and Environment for Statistical Computing","author":"R Core Team","year":"2022"},{"key":"2026092513142921300_btag659-B27","volume-title":"Fundamentals of Biostatistics","author":"Rosner","year":"2006"},{"key":"2026092513142921300_btag659-B28","doi-asserted-by":"publisher","first-page":"e0213684","DOI":"10.1371\/journal.pone.0213684","article-title":"Model to improve specificity for identification of clinically-relevant expanded T cells in peripheral blood","volume":"14","author":"Rytlewski","year":"2019","journal-title":"PLoS One"},{"key":"2026092513142921300_btag659-B29","doi-asserted-by":"publisher","first-page":"910","DOI":"10.1016\/j.ijrobp.2022.05.023","article-title":"Definitive local consolidative therapy for oligometastatic solid tumors: results from the lead-in phase of the randomized basket trial EXTEND","volume":"114","author":"Sherry","year":"2022","journal-title":"Int J Radiat Oncol Biol Phys"},{"key":"2026092513142921300_btag659-B30","doi-asserted-by":"publisher","first-page":"D419","DOI":"10.1093\/nar\/gkx760","article-title":"VDJdb: a curated database of T-cell receptor sequences with known antigen specificity","volume":"46","author":"Shugay","year":"2018","journal-title":"Nucleic Acids Res"},{"key":"2026092513142921300_btag659-B31","author":"Stan Development Team","year":"2024"},{"key":"2026092513142921300_btag659-B32","doi-asserted-by":"publisher","first-page":"267","DOI":"10.1111\/j.2517-6161.1996.tb02080.x","article-title":"Regression shrinkage and selection via the lasso","volume":"58","author":"Tibshirani","year":"1996","journal-title":"J R Stat Soc Ser B Stat Methodol"},{"key":"2026092513142921300_btag659-B33","doi-asserted-by":"publisher","first-page":"210","DOI":"10.1038\/s43018-019-0022-x","article-title":"Immune awakening revealed by peripheral T cell dynamics after one cycle of immunotherapy","volume":"1","author":"Valpione","year":"2020","journal-title":"Nat Cancer"},{"key":"2026092513142921300_btag659-B34","doi-asserted-by":"publisher","first-page":"1461","DOI":"10.1080\/08820139.2025.2550374","article-title":"Biased usage of V\/D\/J genes and clonal diversity in IgG repertoires correlates with disease activity and clinical features in systemic autoimmune diseases","volume":"54","author":"Wang","year":"2025","journal-title":"Immunol Invest"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btag659\/70922587\/btag659.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/9\/btag659\/70922587\/btag659.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/42\/9\/btag659\/70922587\/btag659.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,9,25]],"date-time":"2026-09-25T17:14:40Z","timestamp":1790356480000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/doi\/10.1093\/bioinformatics\/btag659\/8779979"}},"subtitle":[],"editor":[{"given":"Ercument","family":"Cicek","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"editor"}]}],"short-title":[],"issued":{"date-parts":[[2026,8,31]]},"references-count":34,"journal-issue":{"issue":"9","published-print":{"date-parts":[[2026,8,31]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btag659","relation":{},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,9]]},"published":{"date-parts":[[2026,8,31]]},"article-number":"btag659"}}