Recent advances in imaging technologies now allow for real-time tracking of fast-moving immune cells as they search for targets such as pathogens and tumor cells through complex three-dimensional tissues. Cytotoxic T cells are specialized immune cells that continually scan tissues for such targets to engage and kill, and have emerged as the principle mediators of breakthrough immunotherapies against cancers. Modeling the way these T cells move is of great value in furthering our understanding of their collective search efficiency. T-cell motility is characterized by heterogeneity at two levels: (a) Individual cells display different distributions of translational speeds and turning angles, and (b) each cell can during a given track, its motility, switch between local search and directional motion. Despite a likely considerable influence on a motile population’s search performance, statistical models that accurately capture both such heterogeneities in a distinguishing manner are lacking. Here, we model three-dimensional T-cell trajectories through a spherical representation of their incremental steps and compare model outputs to real-world motility data from primary T cells navigating physiological environments. T cells in a population are clustered based on their directional persistence and characteristic “step lengths” therein capturing between-cell heterogeneity. The motility dynamics of cells within each cluster are individually modeled through hidden Markov model to capture within-cell transitions between local and more extensive search patterns. We explore the importance of explicitly capturing altered motility patterns when cells lie in close proximity to one another, through a non-homogenous hidden Markov model.
MillerMJWeiSHParkerI, et al. Two-photon imaging of lymphocyte motility and antigen response in intact lymph node. Science2002; 296: 1869–1873.
2.
MasopustDSchenkelJM. The integration of t cell migration, differentiation and function. Nat Rev Immunol2013; 13: 309.
3.
KrummelMFBartumeusFGérardA. T cell migration, search strategies and mechanisms. Nat Rev Immunol2016; 16: 193.
4.
KienleKLämmermannT. Neutrophil swarming: an essential process of the neutrophil tissue response. Immunol Rev2016; 273: 76–93.
5.
WeningerWBiroMJainR. Leukocyte migration in the interstitial space of non-lymphoid organs. Nat Rev Immunol2014; 14: 232.
6.
GaleanoNLKwanRYQWeningerW, et al. Antigen-specific T cells fully conserve antitumour function following cryopreservation. Immunol Cell Biology2016; 94: 411–418.
7.
TextorJPeixotoAHenricksonSE, et al. Defining the quantitative limits of intravital two-photon lymphocyte tracking. Proc Natl Acad Sci USA2011; 108: 12401–12406.
8.
SchienbeinMGrulerH. Langevin equation, Fokker-Planck equation and cell migration. Bull Math Biol1993; 55: 585–608.
9.
NiculescuITextorJDe BoerRJ. Crawling and gliding: a computational model for shape-driven cell migration. PLoS Comput Biol2015; 11: e1004280.
10.
ManolopoulouIMatheuMPCahalanMD, et al. Bayesian spatio-dynamic modeling in cell motility studies: Learning nonlinear taxic fields guiding the immune response. J Am Stat Assoc2012; 107: 855–865.
11.
TextorJSinnMDe BoerRJ. Analytical results on the beauchemin model of lymphocyte migration. BMC Bioinform2013; 14(6): S10.
12.
BeaucheminCDixitNMPerelsonAS. Characterizing T cell movement within lymph nodes in the absence of antigen. J Immunol2007; 178: 5505–5512.
13.
BaniganEJHarrisTHChristianDA, et al. Heterogeneous cd8+ T cell migration in the lymph node in the absence of inflammation revealed by quantitative migration analysis. PLoS Comput Biology2015; 11: e1004058.
14.
MokhtariZMechFZitzmannC, et al. Automated characterization and parameter-free classification of cell tracks based on local migration behavior. PLoS One2013; 8: e80808.
15.
EsterMKriegelHPSanderJ, et al. A density-based algorithm for discovering clusters in large spatial databases with noise. In: Kdd. vol. 96; 1996. p. 226–231.
16.
KalnisPMamoulisNBakirasS. On discovering moving clusters in spatio-temporal data. In: International Symposium on Spatial and Temporal Databases. Springer; 2005. pp. 364–381.
17.
LeeJGHanJWhangKY. Trajectory clustering: a partition-and-group framework. In: Proceedings of the 2007 ACM SIGMOD international conference on Management of data. ACM; 2007. pp. 593–604.
18.
JebaraTSongYThadaniK. Spectral clustering and embedding with hidden markov models. In: European Conference on Machine Learning. Springer; 2007. pp. 164–175.
19.
ZhangPDengMVan de WegheN. Clustering spatio-temporal trajectories based on kernel density estimation. In: International Conference on Computational Science and Its Applications. Springer; 2014. pp. 298–311.
20.
EsterMKriegelHPSanderJ, et al. Incremental clustering for mining in a data warehousing environment. In: VLDB. vol. 98. Citeseer; 1998. p. 323–333. Sydney.
21.
WangMWangALiA. Mining spatial-temporal clusters from geo-databases. In International Conference on Advanced Data Mining and Applications 2006. pp. 263–270. Xi'an, China: Springer.
22.
WuHRYehMYChenMS. Profiling moving objects by dividing and clustering trajectories spatiotemporally. IEEE Trans Knowl Data Eng2013; 25: 2615–2628.
23.
PanXHeYWangH, et al. Mining regular behaviors based on multidimensional trajectories. Expert Syst Appl2016; 66: 106–113.
24.
PlankMJAuger-MéthéMCodlingEA. Lévy or not? analysing positional data from animal movement paths. In: Dispersal, individual movement and spatial ecology. Springer; 2013. p. 33–52.
25.
DustinML. Stop and go traffic to tune T cell responses. Immunity2004; 21: 305–314.
26.
MéhesEVicsekT. Collective motion of cells: from experiments to models. Integr Biol (Camb)2014; 6: 831–854.
27.
MetznerCMarkCSteinwachsJ, et al. Superstatistical analysis and modelling of heterogeneous random walks. Nat Commun2015; 6: 7516.
28.
ReadMNBaileyJTimmisJ, et al. Leukocyte motility models assessed through simulation and multi-objective optimization-based model selection. PLoS Comput Biology2016; 12: e1005082.
29.
BartolucciFFarcomeniA. A discrete time event-history approach to informative drop-out in mixed latent markov models with covariates. Biometrics2015; 71: 80–89.
30.
RaffaJDDubinJA. Multivariate longitudinal data analysis with mixed effects hidden markov models. Biometrics2015; 71: 821–831.
31.
PattersonTABassonMBravingtonMV, et al. Classifying movement behaviour in relation to environmental conditions using hidden markov models. J Anim Ecol2009; 78: 1113–1123.
32.
LangrockRKingRMatthiopoulosJ, et al. Flexible and practical modeling of animal telemetry data: hidden markov models and extensions. Ecology2012; 93: 2336–2342.
33.
PattersonTAPartonALangrockR, et al. Statistical modelling of individual animal movement: an overview of key methods and a discussion of practical challenges. ArXiv e-prints 2016;1603.07511.
34.
PedersenMWRightonDThygesenU, et al. Geolocation of north sea cod (gadus morhua) using hidden markov models and behavioural switching. Can J Fish Aquat Sci2008; 65: 2367–2377.
35.
ThygesenUHPedersenMWMadsenH. Geolocating fish using hidden markov models and data storage tags. In: Tagging and tracking of marine animals with electronic devices. Springer; 2009. p. 277–293.
36.
VermardYRivotEMahévasS, et al. Identifying fishing trip behaviour and estimating fishing effort from vms data using bayesian hidden markov models. Ecol Modell2010; 221: 1757–1769.
37.
MitaniYSatoKItoS, et al. A method for reconstructing three-dimensional dive profiles of marine mammals using geomagnetic intensity data: results from two lactating weddell seals. Polar Biol2003; 26: 311–317.
38.
LiuYZidekJVTritesAW, et al. Bayesian data fusion approaches to predicting spatial tracks: Application to marine mammals. Ann Appl Stat2016; 10: 1517–1546.
39.
LaplancheCMarquesTAThomasL. Tracking marine mammals in 3d using electronic tag data. Methods in Ecology and Evolution2015; 6: 987–996.
40.
OrellanaMVMatraiPALeckC, et al. Marine microgels as a source of cloud condensation nuclei in the high arctic. Proc Natl Acad Sci USA2011; 108: 13612–13617.
41.
SimsDWSouthallEJHumphriesNE, et al. Scaling laws of marine predator search behaviour. Nature2008; 451: 1098.
42.
HumphriesNEQueirozNDyerJRM, et al. Environmental context explains lévy and brownian movement patterns of marine predators. Nature2010; 465: 1066.
43.
SpeziaL. Bayesian variable selection in non-homogeneous hidden markov models through an evolutionary monte carlo method. Comput Stat Data Anal2020; 143: 106840.
44.
KokiCMeligkotsidouLVrontosI. Forecasting under model uncertainty: Non-homogeneous hidden markov models with pòlya-gamma data augmentation. J Forecast2020; 39: 580–598.
45.
MeligkotsidouLDellaportasP. Forecasting with non-homogeneous hidden markov models. Stat Comput2011; 21: 439–449.
46.
HolsclawTGreeneAMRobertsonAW, et al. Bayesian nonhomogeneous markov models via pólya-gamma data augmentation with applications to rainfall modeling. Ann Appl Stat2017; 11: 393–426.
47.
MaruottiARocciR. A mixed non-homogeneous hidden markov model for categorical data, with application to alcohol consumption. Stat Med2012; 31: 871–886.
48.
BeltmanJBMaréeAFMLynchJN, et al. Lymph node topology dictates T cell migration behavior. J Exp Med2007; 204: 771–780.
49.
HuismanRKoedijkKGKoolCJM, et al. Tail-index estimates in small samples. J Bus Econ Stat2001; 19: 208–216.
50.
EfronBTibshiraniRJ. An introduction to the bootstrap. CRC press, 1993.
51.
HastieTTibshiraniRFriedmanJH, et al. The elements of statistical learning: data mining, inference, and prediction. vol. 2. Springer, 2009: 485–579.
52.
TorkashvandE. Supplement to “modeling three dimensional T cell motility using clustering and hidden markov models”. Stat Methods Med Res2023.
53.
NealRM, et al. Mcmc using hamiltonian dynamics. Handbook of Markov Chain Monte Carlo2011; 2.
54.
BergerJOMorenoEPericchiLR, et al. An overview of robust Bayesian analysis. Test1994; 3: 5–124.
55.
TorkashvandEJozaniMJTorabiM. Constrained Bayes estimation in small area models with functional measurement error. Test2016; 25: 710–730.
56.
SchadDJBetancourtMVasishthS. Toward a principled bayesian workflow in cognitive science. arXiv preprint arXiv:190412765 2019.
57.
ZucchiniWMacDonaldILLangrockR. Hidden Markov models for time series: an introduction using R. 150. CRC press, 2016.
58.
Stan Development Team. RStan: the R interface to Stan 2018;R package version 2.17.3.
59.
KimJOgdenRTKimH. A method to identify differential expression profiles of time-course gene data with fourier transformation. BMC bioinformatics2013; 14: 310.
60.
TsengGCWongWH. Tight clustering: a resampling-based approach for identifying stable and tight patterns in data. Biometrics2005; 61: 10–16.
61.
GelmanACarlinJBSternHS, et al. Bayesian data analysis. 2. FL: CRC press Boca Raton, 2014.
62.
GelmanAHwangJVehtariA. Understanding predictive information criteria for Bayesian models. Stat Comput2014; 24: 997–1016.
63.
TextorJDannenbergKBerryJ, et al. Motilitylab: Quantitative analysis of motion 2016;R package version 0.2-5.
64.
ChoHKBowmanKPNorthGR. A comparison of gamma and lognormal distributions for characterizing satellite rain rates from the tropical rainfall measuring mission. J Appl Meteorol2004; 43: 1586–1597.
65.
JonesCARCibulaMFengJ, et al. Micromechanics of cellularized biopolymer networks. Proc Natl Acad Sci USA2015; 112: E5117–E5122.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
