Estimating Residence Times of Lymphocytes in Ovine Lymph Nodes

The ability of lymphocytes to recirculate between blood and secondary lymphoid tissues such as lymph nodes (LNs) and spleen is well established. Sheep have been used as an experimental system to study lymphocyte recirculation for decades and multiple studies document accumulation and loss of intravenously (i.v.) transferred lymphocytes in efferent lymph of various ovine LNs. Yet, surprisingly little work has been done to accurately quantify the dynamics of lymphocyte exit from the LNs and to estimate the average residence times of lymphocytes in ovine LNs. In this work we developed a series of mathematical models based on fundamental principles of lymphocyte recirculation in the body under non-inflammatory (resting) conditions. Our analysis suggested that in sheep, recirculating lymphocytes spend on average 3 h in the spleen and 20 h in skin or gut-draining LNs with a distribution of residence times in LNs following a skewed gamma (lognormal-like) distribution. Our mathematical models also suggested an explanation for a puzzling observation of the long-term persistence of i.v. transferred lymphocytes in the efferent lymph of the prescapular LN (pLN); the model predicted that this is a natural consequence of long-term persistence of the transferred lymphocytes in circulation. We also found that lymphocytes isolated from the skin-draining pLN have a 2-fold increased entry rate into the pLN as opposed to the mesenteric (gut-draining) LN (mLN). Likewise, lymphocytes from mLN had a 3-fold increased entry rate into the mLN as opposed to entry rate into pLN. In contrast, these cannulation data could not be explained by preferential retention of cells in LNs of their origin. Taken together, our work illustrates the power of mathematical modeling in describing the kinetics of lymphocyte migration in sheep and provides quantitative estimates of lymphocyte residence times in ovine LNs.


INTRODUCTION
One of the peculiar properties of the mammalian adaptive immune system is the ability of its lymphocytes to recirculate between multiple tissues in the body; that is lymphocytes in the blood are able to enter the tissues and after some residence times in the tissues, they return to circulation (1). The pattern of lymphocyte recirculation in general depends on the lymphocyte type (e.g., B or T cell), status of the lymphocyte (resting vs. activated), and perhaps tissues via which lymphocytes are migrating. Naive, antigen-unexperienced lymphocytes, primarily recirculate between secondary lymphoid tissues such as lymph nodes, spleen, and Peyer's patches (2)(3)(4). Following activation after exposure to an antigen, naive lymphocytes become activated and differentiate into effector lymphocytes which are able to access non-lymphoid tissues such as the skin and gut epithelium (2,3,(5)(6)(7).
Why lymphocytes recirculate is not entirely clear (8,9). Because the frequency of lymphocytes specific to any given antigen is in general low and the place of entry of any pathogen is unknown by the naive host, recirculation of lymphocytes may increase the chance to that pathogen-specific cell will encounter its antigen (10,11). Experimental evidence that impairing lymphocyte recirculation influences the ability of the host to respond to infections is very limited. For example, the use of the drug FTY720 (fingolimod) in humans has been associated with a higher incidence of severe infections (12,13). FTY720 prevents lymphocyte exit from lymph nodes, thus, reducing their ability to recirculate (14)(15)(16). However, whether the side-effects of FTY720 is exclusively due to its impact on lymphocyte recirculation is unknown.
The ability of lymphocytes to recirculate between blood and lymph have been nicely demonstrated in now classical experiments by Gowans in rats and later by Hall and Morris in sheep (17,18). Interestingly, subpopulations of lymphocytes can migrate preferentially to different regions of the body, based on their origin as well as their type (19)(20)(21)(22)(23)(24). Molecular interactions between receptors and associated ligands corresponding to the selective entry of lymphocytes to both lymphoid and nonlymphoid tissue have been relatively well-characterized (25)(26)(27)(28). However, the actual kinetics of lymphocyte recirculation have been characterized mostly qualitatively, and we still do not fully understand how long lymphocytes reside in the spleen, LNs, and Peyer's patches, and how such residence times depend on lymphocyte type and animal species.
Understanding lymphocyte migration via lymph nodes (LNs) may be of particular importance for larger animals, including humans, where LNs constitute the majority of the secondary lymphoid tissues (29,30). Lymphocytes may enter the LNs via two routes: from the blood via high endothelial venules (HEVs) or from the afferent lymph draining interstitial fluids from surrounding tissues (1). Experimental measurements suggest that under resting, noninflammatory conditions most lymphocytes (about 80-90%) enter the LNs via HEVs (31). Lymphocytes in LNs exit the nodes with efferent lymphatics which either passes to the next LN in the chain of LNs, or via right or left lymphatic ducts return to the circulation (1,9,32). In contrast, cells can only enter the spleen from the blood and cells exiting the spleen directly return to circulation (33,34).
In the past, to study lymphocyte recirculation via individual LNs sheep or cattle have been used (18,(35)(36)(37)(38)(39)(40). In such experiments, lymphocytes are collected from specific tissue of the animal, e.g., blood, a removed LN, or efferent lymph of a LN, labeled with a radioactive or fluorescent label, and then reinfused back into the same animal (e.g., intravenously, i.v.). The dynamics of the labeled cells is then monitored in the blood, or more commonly, in the efferent lymph of different LNs over time (37, e.g., see Figure 1A). The dynamics of the labeled cells in the efferent lymph of various LNs follow a nearly universal pattern-the number of labeled cells increases initially, reaches a peak and then slowly declines over time (20,31,42,43, and see Figures 2C,D). Given that in many such experiments, the peak of labeled cells in the efferent lymph is reached in 24 h and declines slowly (e.g., Figure 2C), it can be interpreted that the average residence time of lymphocytes in the ovine lymph nodes is about 48 h. To our knowledge, the actual residence time in ovine LNs has not been regularly reported. Interestingly, with the use of mathematical modeling an accurate quantification of how long lymphocytes spend in LNs in sheep has been recently performed (44).
In their pioneering study Thomas et al. (44) analyzed data on migration of lymphocytes via individual ovine LNs. To quantify these dynamics, the authors developed a mathematical model which considers cell migration via a LN as a random walk between multiple sub-compartments in the LN. In the model cells entering the LN start in the first sub-compartment, progress via the series of subcompartments by a random walk and eventually exit the node by leaving the last subcompartment (44). As their model allows for the possibility for cells to spend variable lengths of time in the LN, it can naturally explain the long duration of labeled lymphocytes exiting the LN. By fitting the model to several different sets of data, the authors concluded that the average residence time of lymphocytes in ovine LNs is about 31 h (44).
In addition to the average residence time the distribution of residence times may be important. In particular, if lymphocytes that just entered the LN have the same chance of exiting it as lymphocytes that already spent some time in it (i.e., distribution of residence times is exponential), this could suggest that exit of lymphocytes from LN is a simple Poisson-like stochastic process. Indeed, one recent study suggested that residency time of naive CD4 and CD8 T cells in LNs of mice is exponentially distributed (45). In contrast, if lymphocytes require some time to be spent in the LN, for example, to acquire the ability to exit the node, then the distribution of residence time cannot be exponential. Recent work from our group suggests that residence time of thoracic duct lymphocytes in LNs of rats is not exponential and is best described by a gamma distribution with the shape parameter k = 2 or 3 (41). Our analysis of data on kinetics of lymphocyte exit from inguinal LNs of photoconvertable Kaede mice did not allow to firmly establish the shape of the residence time distribution (9). Whether the distribution of residence time of lymphocytes in ovine LNs is exponential or more complex is unknown.
In this paper, we formulated a series of mathematical models aimed at describing the kinetics of recirculation of lymphocytes in sheep. All models take into account basic physiological constrains on lymphocyte migration, for example, that lymphocytes enter the LNs continuously from the blood and lymphocytes that exit LNs return back to circulation. The models were fitted to a series of experimental data from previously published studies on lymphocyte migration in sheep. Our results suggest that the distribution of residence times of lymphocytes in ovine LNs is best described by a nonexponential distribution with estimated average residence times being 12-22 h depending on the type of lymphocytes used in experiments. The long-term presence of labeled lymphocytes in efferent lymph of cannulated LNs was explained by a continuous entry of new cells from the blood to the LN, thus, simplifying a previous modeling result (44). Overall, our analysis provides a quantitative framework to estimate kinetics of lymphocyte recirculation using measurements of lymphocyte numbers in the blood and efferent lymph of ovine LNs. Such a framework may be useful to understand the efficacy and potential limitations of immunotherapies involving adoptive transfer of T cells in humans (46-48).

Experimental Data
For our analyses we digitized the data from three publications (20,42,49). We describe in short design of these experiments and how the data have been collected. For more detail the reader is referred to the original publications.

Lymphocyte Dynamics in Blood and Efferent Lymph (Dataset #1)
Experiments of Frost et al. (42) have been performed with 4-24 month old Alpenschalf or Schwarzkopf sheep ( Figure 1A). Lymphocytes were collected from the efferent lymph of different LNs (e.g., prescapular) then labeled with 51 Cr and resuspended in buffer. A surgical procedure to install an indwelling cannula into the jugular vein was performed to allow for i.v. infusions and collection of blood samples. The efferent duct of specific lymph nodes was also cannulated to allow for collection of lymph that passed through the node. The number of labeled lymphocytes in efferent lymph was determined by washing cells in a buffer and counting with a gamma scintillation counter. Labeled lymphocytes in venous or peripheral blood were counted similarly by a gamma scintillation counter then compared to the activity of plasma from the same volume of blood.
Labeled lymphocytes were injected i.v. and the efferent lymph of prescapular LN (pLN) was collected. Lymph was collected as described above at 20-min intervals for 3 h or daily for about 2 weeks. The number of labeled lymphocytes in the lymph node was expressed as cpm per 10 7 cells collected. Peripheral blood samples were only collected for certain experiments and they were taken at 10 min, 30 min, 1 h, and 3-h intervals thereafter. The number of labeled lymphocytes in peripheral blood was expressed as cpm/ml of whole blood minus the activity in the plasma. Cells exiting the lymph node were measured for 120 h, while cells in the blood were measured for 90 h, then once more at the 120th h. For this reason, we discuss lymphocyte migration kinetics in terms of short-and long-term migration experiments. The short-term experiments include the dynamics of the labeled lymphocytes in both blood and efferent lymph for the first 90 h, and the long-term dataset considers all available data (measurement up to 120 h).
According to "Blood Volume of Farm Animals, " Hampshire sheep less than a year to 3 years old have on average 6.3-5.8 ml blood per 100 g weight and an average weight of 92-156 lbs (50). This results in about 2.63 L of blood for Hampshire sheep less than a year old and about 4.1 L for Hampshire sheep 2-3 years old. The breed of sheep used in this study weight from 60 kg (Alpenschaf) to 100 kg (Schwarzkopf) at maturity, so it is reasonable to assume that these total blood estimates are less than those calculated for Hampshire sheep. We assume the volume of blood in the sheep is an average of these given volumes, specifically V = 3.4 L. This estimate was used to convert the estimate of the number of labeled lymphocytes per ml of blood to the total number of lymphocytes in the whole blood.
In one set of experiments [ Figure 1 in (42)] 2.5 × 10 9 labeled cells (representing 12.6 × 10 6 cpm) were injected i.v. into sheep corresponding on average RL o = 2.5 × 10 9 cells/(12.6 × 10 6 cpm) = 200 cells/cpm. Because in the original data, the cells in the blood at time t (RL t ) were measured in cpm/ml of blood, the total number of cells in the blood B at time t was calculated using the formula where RL o = 200 cells/cpm, V = 3.4 L and RL t was digitized from Figure 1 of Frost et al. (42). In the same experiments labeled cells exiting the pLN were measured as cpm/10 7 cell with each data point summing 3 h of cell collected every 20 min. Therefore, the amount of labeled cells exiting the pLN per hour (C t ) is one third of this number. The total output of the pLN is given in Figure 3 in Frost et al. (42) and was estimated to be f = 10.85 × 10 7 cell/hr. We calculate total cells exiting the pLN per hour at time t as where RL o = 200 cells/cpm. The final data (dataset1.csv, given as Supplement to the paper) includes changes in the total number of labeled lymphocytes in the peripheral blood and the number of labeled lymphocytes exiting pLN per hour over time.

Migration of Lymphocytes From Afferent to Efferent Lymph (Dataset #2)
To further investigate how lymphocytes migrate via ovine LNs we digitized data shown in Figure 4 of Young et al. (49).
In these experiments different subsets of lymphocytes (CD4 T cells, CD8 T cells, γ δ T cells, and B cells) were collected from efferent prescapular lymph and labeled with PKH-26 or CFSE to distinguish between different cell subsets. A maximum of 2 × 10 6 cells of each type were infused into two popliteal afferent lymphatics over 1 h. Cells were collected as they exited the cannulated efferent lymph of the popliteal LN (poLN), and phenotyped. The final data (dataset2.csv, given as Supplement to the paper) includes the percent of labeled cells found in the efferent lymph of the poLN over time.

Migration of T Lymphocytes via Skin-Draining and Gut-Draining Lymph Nodes (Dataset #3)
The final set of experimental data we used come from recirculation experiments of Reynolds et al. (20) with young sheep (Figure 6). Lymphocytes were isolated by cannulating the skin-draining pLN or the ileal end of the gut-draining mesenteric lymph node chain (mLN). Collected lymphocytes were enriched for T cells and labeled with FITC or RITC based on origin. Labeled cells were re-infused into the same animal i.v. and cell frequency (given in labeled cells per 10 4 cells) was reported for 240 h [ Figure 1 in (20)]. Because the authors did not report the overall cell output in efferent lymph of pLN and mLN, we used the provided numbers of the frequency of labeled cells in the lymph for fitting the models. Existing data suggest similar output of lymphocytes from pLN (1 − 5 × 10 8 cells/h) or mLN [1 − 10 × 10 8 cells/h, (37)] which in part justifies our approach. The final data (dataset3.csv, given as Supplement to the paper) includes the number of labeled cells from skin or intestinal lymph per 10 4 cells found in efferent lymph of pLN or mLN over time during cannulation. The raw data from the cited papers were extracted using digitizing software Engauge Digitizer (digitizer.sourceforge.net).

Basic Assumptions
In our models we assume that blood is the main supplier of lymphocytes to other tissues and that when exiting these tissues, lymphocytes return to the blood (41,51). We also assume that lymphocyte have the same exit kinetics for all LNs (except in the analysis of the dataset #3), and that cell infusion does not impact lymphocyte migration via individual LNs. There is a disagreement whether residency time of lymphocytes depends on the LN type which may be due to differences in animal species or lymphocyte subsets used, or in experimental techniques (41,45).

Recirculation model
To predict the dynamics of i.v. transferred labeled lymphocytes in the blood and efferent lymph of pLN we extended a previously proposed compartmental model describing recirculation kinetics of lymphocytes in the whole body (41). The model ( Figure 1B) predicts the number of i.v. transferred lymphocytes in the blood (B), spleen (S), lymph nodes (L), and other non-lymphoid tissues (T). In the model we ignored migration of lymphocytes via vasculature of the lung and liver since previous work suggested that resting lymphocytes pass via these tissues, at least in rats, within a minute (41). This is different from migration of activated lymphocytes via these tissues which could take hours (9). Yet, it should be emphasized that because in experiments that we have analyzed migration of lymphocytes via lung and liver vasculature has not been measured, it is possible that lung/liver may represent "spleen" or other "non-lymphoid tissues" in our model. We explore the impact of including lung/liver in our recirculation model on estimates of lymphocyte residency time in the LNs in section 4.
In the model, cells in the blood can migrate to the spleen, lymph nodes, or other tissues at rates m Bi and cells can return to circulation from these tissues at rates m iB where i = S, L, T. When exiting spleen or non-lymphoid tissues lymphocyte follow the first order kinetics, so the decline of cells in the tissues in the absence of any input is given by an exponential function. This in part is based on our previous work suggesting of migration of thoracic duct lymphocytes via spleen can be described as first order kinetics (41). In contrast, migration of lymphocytes via LNs may not follow the first order kinetics [e.g., (41)] and thus was modelled by assuming k sub-compartments in the nodes with equal transit rates m LB . Such sub-compartments may represent different areas in the LNs, for example, paracortex and medulla. The "recirculation" model was adopted from our previous study (41) to describe the migration kinetics of labeled lymphocytes in the animal. Lymphocytes in the blood B may migrate to three (n = 3) major tissue compartments: spleen S (at a rate m BS ), lymph nodes L (at a rate m BL ), or other peripheral tissues T (at a rate m BT ). A fraction λ of cells migrating to lymph nodes migrate to the cannulated LN (Lo), thus can be measured, while remaining cells (1 − λ) migrate to other lymph nodes (L). Exit of lymphocytes from the spleen or peripheral tissues follows first order kinetics at rates m SB and m TB , respectively. In contrast, residency times in LNs are gamma distributed and are modeled by assuming k sub-compartments (in the figure k = 3) where exit from each sub-compartment is given by the rate m LB . (C) The alternative "blood-LN dynamics" model allows for cells in the blood (B) to enter the cannulated LN at a rate λm BL and exit the LN by passing via k sub-compartments at a rate m LB . Cells in the blood also leave the blood at a rate α (in case of a single exponential decline).
Mathematically, we used the sub-compartments to model nonexponential residency time of lymphocytes in the LNs.
To describe accumulation and loss of labeled lymphocytes in the cannulated lymph nodes we assume that a fraction of lymphocytes λ migrating to lymph nodes migrate to the cannulated node Lo 1 (e.g., pLN, Figure 1B) while 1 − λ cells migrate to other LNs (L 1 ). We assume that cells do not die but the process of migration to non-lymphoid tissues with no return back to circulation is equivalent to cell death. We did consider several alternative models in which death rate was added to the model (see section 4). Taken together, with these assumptions the mathematical model for the kinetics of lymphocyte recirculation in sheep is given by equations: where cells exiting the cannulated lymph node, m LB Lo k do not return to the blood because they are sampled, and m LB Lo k is the rate of labeled lymphocyte exit from the sampled lymph node (which is compared to experimental data, e.g., column 2 in the dataset #1, and see Figures 2C,D).
Our experimental data were restricted only to labeled lymphocytes found in the blood and efferent lymph (or only efferent lymph). It may be argued that a minimal recirculation model to describe such data should only involve blood and LNs (i.e., n = 1). Also, the number of sub-compartments k in LNs is also unknown. Therefore, in our analyses we fitted a series of models assuming different values for k and n to the data from (42) (dataset #1) and compared the quality of the model fit to data using AIC (see section 3).
The average residence time of lymphocytes in the spleen or non-lymphoid tissues is 1/m SB and 1/m TB , respectively. The residence time of lymphocytes in the LNs is given as RT = k/m LB . The initial number of labeled cells in the blood varied by experiment and is indicated in individual graphs. In some fits parameter λ could not be identified from the data and thus was fixed (indicated by absent predicted confidence intervals). The rest of the parameters were fit.

Blood-LN dynamics model
Many studies of lymphocyte recirculation that reported kinetics of accumulation and loss of labeled lymphocytes in efferent  Figure 1C) with j = 2 and k = 3 resulted in the best fit. The estimated average RT of lymphocytes in the pLN is shown. Fits of models that assume different numbers of tissue compartments, different number of sub-compartments in LNs, or different numbers of exponential functions are shown in Table 1 and Table S1. Parameters for the best fits of these models are given in Table 2. lymph of ovine LNs did not report dynamics of these cells in the blood which is the major limitation of such studies. Therefore, to gain insight into whether LN cannulation data alone can be used to infer lymphocyte residence times in the LNs we propose an alternative model which only considers lymphocyte dynamics in the blood, one cannulated LN, and the efferent lymph of the cannulated LN. In this "blood-LN dynamics" model ( Figure 1C) the dynamics of lymphocytes in the blood is described by a phenomenological function B given as a sum of j declining exponentials. The rationale to use such a model stems from the kinetics of labeled lymphocyte dynamics in the blood observed in (42) and other studies [Figures 2A,B and (44)] even though actual dynamics in other experimental systems may not follow the same pattern. The dynamics of labeled cells in the sampled lymph node is thus driven by the continuous entry of labeled cells from the blood into the LN. The model is then given by following equations X i and α i are the initial values and the rate of decline in the i th exponential function, λm BL is the overall rate at which lymphocytes from the blood enter the LN, and m LB is the rate at which lymphocytes move between k sub-compartments in the LN and exit the LN. The rate at which lymphocytes exit the LN and thus are sampled in the LN efferent lymph is m LB Lo k . It should be noted that parameters λ, m BL and X i in many cases are not identified from the data on lymphocyte dynamics in the efferent lymph, thus, the main parameter that we are interested in is the average residence time of lymphocytes in the LNs given by RT = k/m LB . Because the dynamics of cells in the blood is generally unknown when fitting the model predictions to data, we varied the number of exponential functions j = 1, 2, 3 and compared the quality of fits of different models using AIC.

Migration When Cells Are Injected Into Afferent Lymph
In one set of experiments migration of labeled lymphocytes via LNs was measured by directly injecting lymphocytes into the afferent lymph of a LN and observing accumulation and loss of these cells in the efferent lymph of the LN. To use these data to estimate the lymphocyte residence time in the LN we assume that cells injected into afferent lymph A migrate into the lymph node at a rate m A , and then the cells migrate via each of k subcompartments in the LN at a rate m LB . With these assumptions the dynamics of cells in the afferent lymph and the LN are given by equations: where initially all labeled cells were in the afferent lymph. As previously stated, the rate of lymphocyte exit from the LN via efferent lymph is given by m LB Lo k . The average residence time of lymphocytes in the LN is then RT = k/m LB .

Homing to Different Lymph Nodes
In the final set of experiments Reynolds et al. (20) collected lymphocytes from efferent lymph of pLN or mLN, labeled and then re-infused the collected cells i.v. into the same animal. The labeled cells were then collected in the efferent lymph of the pLN and mLN. Because the authors did not report the dynamics of labeled cells in the blood, we extended the "blood-LN dynamics" model (see  to describe cell migration from the blood to the efferent lymph of two LNs. The number of labeled lymphocytes found in the i th sub-compartment of the pLN and mLN are given by Lo 1,i and Lo 2,i , respectively: where m L1B and m L2B are the rate of lymphocyte exit from the pLN and mLN, respectively, m BL1 and m BL2 are the rates of lymphocyte entry from the blood to pLN and mLN, respectively, and j = 1, 2 in fitting models to data. Because the data clearly showed the difference in accumulation of lymphocytes in different LNs, we considered two alternative explanation for this difference. In one model we assume that the difference in kinetics is due to differences in the rate of lymphocyte entry into specific LNs (m BL1 = m BL2 ) while residence times are identical in the two LNs (m L1B = m L2B ). In the alternative model, the rate of entry into the LNs are the same but residence times may differ (m BL1 = m BL2 and m L1B = m L2B ).

Statistics
The models were fitted to data in R (version 3.1.0) using modFit routine in FME package (version 1.3.5) by log-transforming the data (single or two different measurements) and model predictions and by minimizing the sum of squared residuals. For example, when fitting the recirculation model to the data on lymphocyte numbers in the blood and efferent lymph from Frost et al. (42) the SSR was calculated in the following way: Frontiers in Immunology | www.frontiersin.org where experimental measurements are given in Equations (1) and (2) and model predictions are from Equations (3) and (8), respectively, for l measurements at times t i , i = 1 . . . l. Numerical solutions of the system of equations were obtained using ODE solver lsoda (from the deSolve package) with default absolute and relative error tolerance. Different algorithms such as BFGS, L-BFGS-B, or Marquart in the modFit routine were used to find parameter estimates. Discrimination between alternative models was done using corrected Akaike Information Criterion, AIC (52) where SSR is the sum of squared residuals, N is the number of data points, and p is the number of model parameters fitted to data. The model with the minimal AIC score among all tested models was viewed as the best fit model, but a difference of AIC score of 1-3 between best fit and second best fit models was generally viewed as not significant (52). Predicted 95% confidence intervals for estimated parameters were calculated as ±2σ with standard deviation σ provided for each parameter by the modFit routine.

Estimating Lymphocyte Residency Time in the LNs Using Lymphocyte Dynamics in Efferent Lymph
To gain quantitative insights into the kinetics lymphocyte migration via sheep lymph nodes we first fitted our recirculation model (given by Equations 3-9) to the "short-term migration" data on lymphocyte dynamics in the blood and efferent lymph of the prescapular LN from Frost et al. (42) (see section 2 for more detail on the data). While the overall structure of the recirculation model was defined by the number n of different compartments through which lymphocytes could recirculate (Figure 1B), we first investigated how many such compartments are in fact necessary to describe the experimental data by varying n between 1 and 4. Additionally, we tested how many sub-compartments k in the LNs are needed for best description of the data (see . The analysis revealed that n = 3 tissue compartments and k = 3 sub-compartments in the LNs are needed to adequately describe the dynamics of labeled cells in the blood and efferent lymph ( Table 1). Such a model could accurately describe simultaneously the loss of labeled cells in the blood and accumulation and loss of labeled cells in the efferent lymph (Figures 2A,C). The model predicted the existence of two recirculation compartments with average residence times of 2.4 and 19.5 h, with the latter compartment corresponding to LNs in the sheep. The nature of the first compartment is unclear but given the estimated residence time it is likely that it represents the spleen [e.g., see (41)]. The final third compartment was needed to explain the long-term loss of labeled cells from the blood at a rate of about d = 0.02/h. The predicted rate of lymphocyte migration to tissues (m BS + m BL + m BT ≈ 1.4 − 1.6/h ≫ d, see Table 2) is higher than the observed rate d because of the return of lymphocytes that had migrated to lymphoid tissues back to circulation. The model also naturally explains the long-term decline in the number of labeled lymphocytes found in the efferent lymph of the cannulated pLN which is simply driven by the decline of labeled cells in the blood. The analysis also suggests that none of the obvious characteristics of the distribution of the lymphocyte exit rate from the LN such as the time of the peak or the average of the overall distribution (e.g., see Figure 2C) accurately represent the average residence time. This result strongly suggests that to accurately estimate lymphocyte residence times from LN cannulation experiments it is critical to use appropriate mathematical models.
Assuming a smaller (k = 1) or a larger (k = 4) number of sub-compartments in the LNs resulted in poorer fits of the data ( Table 1). The intuitive reason of why the model in which lymphocyte residence times in LNs are exponentially distributed (k = 1) does not fit the data well follows from the rapid loss of labeled lymphocytes in the blood within the first hours after lymphocyte transfer (Figure 2A). Rapid decline in the number of labeled lymphocytes in the blood reduces the rate at which new labeled cells enter the pLN which would have resulted in a relatively rapid exit of cells from the pLN for exponentially distributed residence time. Similarly, the model in which there are too many sub-compartments would force the distribution of cells in the efferent lymph to be even broader, thus, also resulting in poorer fit. Thus, this analysis suggests that migration of lymphocytes via LNs is not described by a simple exponential function and there is a requirement for lymphocytes to spend some minimal time in LNs before exiting into circulation.
It is interesting to note how the dynamics of labeled lymphocytes in the blood may be used to infer recirculation kinetics of cells. Indeed, the initial rapid decline of the number of labeled lymphocytes is explained in the model by migration to secondary lymphoid tissues and change in the decline rate at 2-3 h after lymphocyte transfer is naturally explained by the exit of initially migrated cells from one of the compartments (most likely spleen) back to the blood. Thus, lymphocyte kinetics in the blood suggests residence time in first compartment of about 2-3 h ( Table 2).
The recirculation model makes a strong assumption that the dynamics of labeled lymphocytes in the blood and efferent lymph are due to migration of lymphocytes into and out of different tissues (Figure 1B). When the experimental data is provided for the dynamics of labeled cells in the blood as in dataset #1, we are able to use Equations (3)-(9) to accurately describe the data and estimate lymphocyte residence times in various tissues. However, when experimental data do not contain measurements of the dynamics of labeled cells in the blood, predictions of the recirculation model remain speculative. Therefore, to estimate residence times of lymphocytes in LNs in the absence of such data, we developed an alternative mathematical model. This model involves a smaller number of assumptions, the major of which is the physiological constraint that most lymphocytes enter lymph nodes from blood (2). Intuitively, however, this model allows us to estimate the time taken by cells to migrate from the blood to the efferent lymph of a LN irrespective of the specific route of this migration.
In the"blood-LN dynamics" model, the dynamics of labeled lymphocytes in the blood is described phenomenologically as a sum of several exponential functions, and by fitting a series of such models we found that the dynamics of labeled cells in the first 90 h after cell transfer is best described by a sum of two exponentials (Table S1 and Figure 2B). The model predicted a rapid initial loss of lymphocytes in the blood at a rate of 2/h (halflife time of about 21 min) and a slower loss rate of 0.02/h after the first 4 h (half-life time of 35 h, Figure 2B and Table 2).
By fitting a series of mathematical models in which the number of sub-compartments in the pLN was varied, to the data on dynamics of labeled lymphocytes in the efferent lymph we found that k = 3 sub-compartments provided fits of the best quality (Table S1 and Figure 2D). Importantly, the model predicted the average residence time of lymphocytes in the pLN of 19.6 h which is nearly identical to the value found by fitting recirculation model to the same data ( Table 2).
Both recirculation and blood-LN dynamics models were then fitted to the long-term migration data in which the dynamics of labeled cells in efferent lymph was measured continuously for 120 h while in the blood there was an extra measurement at 120 h (Figures 3A-D). The recirculation model with n = 3 tissue compartments and k = 3 sub-compartments was able to accurately describe the data (Figures 3A,C) although the model underestimated the number of labeled lymphocytes in the blood at 120 h post-transfer ( Figure 3A and Table S2). Interestingly, the model required a positive rate of lymphocyte return from the "third" tissue back to circulation ( Figure 3A and Table 2). In the absence of lymphocyte return from the tissue compartment In the recirculation model the three tissue compartments are suggested to be spleen, LNs, and other non-lymphoid tissues and the migration rates from the blood to these compartments are denoted as m ij with i, j = B, S, L, T. In the blood-LN dynamics model it was not possible to estimate accurately the initial number of labeled lymphocytes in the blood (X 1 ), so that parameter was fixed to X 1 = 2.5×10 9

cells. Residence times in LNs were calculated as RT = k/m LB and as 1 m iB for other compartments (i = S, T).
the model poorly matches the number of labeled cells in the blood (results not shown). Importantly, the recirculation model predicted similar average residence times of lymphocytes in the first two compartments (representing spleen and LNs) to that of the model fitted to short-term dataset ( Table 2).
Perhaps unsurprisingly, to describe the dynamics of labeled lymphocytes in the blood over 120 h the sum of three different exponential functions was required (Table S1). Furthermore, the model with k = 3 sub-compartments in the pLN was able to describe the dynamics of labeled lymphocytes in efferent  (42). In these experiments, there were additional measurements for labeled lymphocytes in the blood and efferent lymph (compare to Figure 2A). The best fit recirculation model assumed n = 3 and k = 3. The average residency times RT i is given for each compartment. The best fit blood-LN dynamics model had j = 3 exponential decay functions and k = 3 sub-compartments for the LNs. The best fit models were determined via a series of trials (see Tables S1, S2). Parameter estimates and 95% CIs are shown in Table 2.
lymph with best quality (Table S1 and Figures 3B,D). The model predicted the average residence time of lymphocytes in LNs to be 19.6 h which is consistent with results from the recirculation model fitted to the same data or the models fitted to short-term migration data.
Taken together, analysis of data from

Migration of Lymphocytes From Afferent to Efferent Lymph Suggests Non-exponentially Distributed Residence Time in LN
Analysis on the dynamics of labeled lymphocytes transferred i.v. into sheep suggested that migration of lymphocytes via LN follows a multi-stage process which can be described as cell migration via identical sub-compartments (Figures 1B,C). Since the time it takes for lymphocytes to cross the endothelial barrier and enter LNs is very short [few minutes, (53)], the finding that distribution of lymphocyte residence times are not exponential could still be due to some unknown processes. Therefore, to further investigate the issue of the distribution of residence times of lymphocytes in ovine LNs we analyzed experimental dataset #2 (49). In these experiments, Young et al. (49) isolated lymphocytes from the efferent lymph of the pLN, labeled and injected the cells into the afferent lymph of the popliteal LN (poLN), and then measured exit of the labeled cells the efferent lymph of the poLN (Figure 4 and see section 2 for experimental data detail). Cells, injected into the afferent lymph, cannot move to any other tissue but the draining LN, and thus, such data allow to directly evaluate kinetics of lymphocyte migration via individual LN.
To describe these data, we adapted the blood-LN dynamics model to include migration of labeled lymphocytes from the afferent lymph to the LN and then to the efferent lymph (Equations 13-15). The model has 3 unknown parameters that must be estimated from the data (A(0), m A , m LB ). Unfortunately, the original data for cell dynamics for individual animals were not available, and the digitized data only included 3 time points which does not allow accurate estimation of all model parameters (results not shown). Therefore, to investigate the dynamics of labeled cells in the efferent lymph we fitted a series of mathematical models with a varying number of subcompartments k in the LN and average residence times RT = k/m LB fixed to several different values to the experimental data (Table S3). Analysis revealed that several sub-compartments are needed for accurate description of the data and the actual number of sub-compartments varied for different cell subtypes, but was never less than k = 3 ( Table S3). The expected residence times also varied with the cell type but overall were within 18-20 h range which is consistent with the previous analysis of Frost et al. (42) data (Figure 4).
By fitting the data with the model in which the number of sub-compartments k was varied we found that the estimated residence time of lymphocytes in the poLN was dependent on the assumed k (Table S4). This is consistent with our recent result on estimating residence time of T and B lymphocytes in LNs of mice using the data from photoconvertable Kaede mice (9). Interestingly, the model fit predicted a relatively slow movement of lymphocytes from the afferent lymph to the LN, which is determined by the parameter m A (1/m A ≈ 5 h) and was dependent on the number of sub-compartments. These results also support our conclusion that residence times of lymphocytes in ovine poLN are not exponentially distributed and the average residence time for different lymphocyte subsets is around 20 h.

Impact of Lymphocyte Kinetics in the Blood on Estimates of the Lymphocyte Residency Time in the LNs
Our analysis of the Frost et al. (42) data demonstrated the usefulness of having measurements of the dynamics of labeled lymphocytes both in the blood and efferent lymph of a specific cannulated LN. Unfortunately, many published studies that we reviewed lacked measurements of lymphocyte counts in the blood and only recorded cell numbers in the efferent lymph, often as the percent of labeled cells in the overall population. An important question is whether the estimates of the residency time of lymphocytes in LNs, found by fitting mathematical models to the data on lymphocyte counts in efferent lymph only, depend on the assumed lymphocyte dynamics in the blood.
In several different experiments the decline of i.v. injected labeled lymphocytes in the blood is bi-exponential with the rapid decline in cell numbers within a few hours and slower decline in the next days (44, see Figure 2B). By fitting the blood-LN dynamics model (with k = 3) to the data on the dynamics of labeled lymphocytes in efferent lymph (shown in Figure 2D) we found that fits of similar quality could be obtained independently whether the dynamics of labeled cells in the blood follow either a single or bi-exponential decline (results not shown). However, the estimates of the average residence time of lymphocytes in LNs was dependent on the assumed model of lymphocyte dynamics in the blood; namely, assuming a bi-exponential decline resulted in longer average residence times (results not shown).
To investigate this issue further we fitted the blood-LN dynamics model (with k = 3) assuming that the number of labeled cells in the blood follows an exponential decline, to the data on labeled cell dynamics in efferent lymph. In this analysis we either fitted the rate of cell decline in the blood (α 1 ) or fixed it to different values (Figure 5). We found that the decline rate of labeled cells in the blood has a dramatic impact on the quality of the model fit of the data as well as on estimates of the average residence times (Figure 5B). In particular, assuming that the number of labeled cells remains constant in the blood (α 1 = 0) predicts a constant output of labeled cells in efferent lymph, and as the result, failed to accurately describe the data (Figure 5B). Similarly, assuming that the loss of labeled cells occurs relatively rapidly (α 1 = 0.05/h) also results in poor fits of the data and longer average residence time of lymphocytes in the LN (Figure 5B). However, allowing the rate of lymphocyte loss α 1 to be fitted resulted in good fits of the data further suggesting that the long-term dynamics of labeled cells in the efferent lymph is the consequence of cell dynamics in the blood.

Lymphocytes Migrate More Rapidly to LNs Which the Cells Recently Exited
All data analyzed so far have been for lymphocytes isolated from pLNs which migrate back to pLN or poLNs. An important question is whether the average residence time of lymphocytes varies across types of LN. Indeed, previous analysis of migration of naive T cells in mice suggest that lymphocytes spend less time in gut-draining mesenteric LNs (mLNs) than in skin-draining pLNs (45). In contrast, another study suggested similar residency times of thoracic duct lymphocytes in pLN and mLN of rats (41). To address this issue we analyzed experimental dataset #3 (20).
In their studies, Reynolds et al. (20) isolated T lymphocytes from efferent lymph of pLN or mLN, labeled them with different fluorescent dyes, re-injected the cells, and measured their exit in the efferent lymph of pLN and mLN (Figure 6 and see section 2 for experimental data detail). The data showed that T cells isolated from pLN accumulate to higher numbers in the efferent lymph of pLN as compared to cells from mLN and vice versa (Figure 7). There could be at least two alternative explanations for such differential accumulation of cells in LNs of their origin: preferential migration or preferential retention. According to the preferential migration hypothesis, cells from pLN have a higher rate of entry into pLN than the rate at which cells from mLN enter the pLN (and vice versa). In contrast, in the preferential retention hypothesis, cells from pLN have a longer residence time in pLN as compared to cells from mLN (and vice versa).
To discriminate between these alternative hypotheses, we fitted the blood-LN dynamics model to these data. The blood-LN dynamics model was chosen because of its relative simplicity and because Reynolds et al. (20) did not report dynamics of transferred lymphocytes in the blood, which was needed for accurate estimation of parameters of the recirculation model. Specifically, we assumed that the T cell kinetics follow biexponential decline (j = 2 in Equation 10) and that lymphocytes must traverse via k = 3 sub-compartments in the LNs. In the "preferential migration" model we fixed the average residence times of lymphocytes in LNs for cells from pLN and mLN (determined by the parameter m LB ) but allowed different rates of entry into the LN from the blood (determined by the parameter m BL , see . This model can accurately describe experimental data (Figures 7A,B). Interestingly, the model predicted 2 fold higher entry rate into pLN by cells of pLN origin as compared to cells of mLN origin, and 3 fold higher entry rate into mLN by cells of mLN origin, as compared to  Table S4.   (Table S5). Importantly, assuming identical average residence time of T cells from pLN in skin-draining or gut-draining LNs resulted in fits of excellent quality suggesting the average residence time of T cells from pLN does not depend on the LN type. However, T cells from mLN migrated via LNs nearly 2 fold faster than T cells from pLN suggesting that the average residence time does depend on the origin of T cells.
In the alternative "preferential retention" model we fixed the rate of lymphocyte entry from the blood to the LNs and allowed the residence times (or more precisely, the rate of exit of T cells from the LNs) to vary depending on the LN type. This model failed to accurately describe the data (Figures 7C,D) suggesting that the data cannot be explained solely by increased retention of cells in the LN of their origin. Importantly, allowing both entry

DISCUSSION
It is well understood that some lymphocytes are able to recirculate between blood and secondary lymphoid tissues such as lymph nodes. In part, this understanding came from multiple experiments on lymphocyte migration from the blood to efferent lymph of various LNs in sheep. Yet, while the data on the kinetics of lymphocyte migration via individual LNs have been published, quantitative interpretation of these data has been lacking until recently. In particular, the average residence time of lymphocytes in the ovine LNs remained largely unknown and there was incomplete understanding of why labeled lymphocytes persisted in the efferent lymph of cannulated LNs.
The first attempt known to us to explain lymphocyte dynamics in efferent lymph during cannulation experiments in sheep was by Thomas et al. (44) who modeled lymphocyte migration via the LN as a random walk. The model suggested that longterm detection of labeled lymphocytes in efferent lymph was due to inability of some lymphocytes to exit the LN. Here we formulated several alternative mathematical models, based on the basic understanding of lymphocyte recirculation in mammals which accurately explain the cannulation data, and proposed an alternative explanation for long term detection of labeled cells in efferent lymph. Namely, because labeled cells persist in the blood, continuous reentry of such cells into the LN can naturally explain long-term persistence of labeled cells in the efferent lymph.
Our mathematical modeling approach allowed to provide novel estimates of lymphocyte residence time in ovine lymph nodes which vary between 12 and 20 h depending on lymphocyte type and being approximately independent of the type of LN (e.g., skin-or gut-draining, Figures 2, 7). Furthermore, the combination of data and mathematical model predicted an existence of a compartment with a shorter residence time, about 2-3 h, which we propose is likely to be the spleen. Indeed, we recently published a similar estimate of lymphocyte residence time in the spleen of rats (41). However, this prediction remains to be tested experimentally, and it remains possible that the unobserved compartment with the residency time of 2-3 h represents vasculature of the lung or liver. We consider this alternative to be unlikely given that we previously estimated that resting lymphocytes spend a rather short time in the lung/liver vasculature of rats [less than 1 min, (41)].
Parameter estimates also suggest a relatively short residence time of lymphocytes in the blood (e.g., for long-term migration data 1/(m BS + m BL + m BT ) ≈ 1.6/h or RT B = 28 min). This is relatively similar to previous estimates (9,41).
While we did not specifically model lymphocyte migration via lung or liver vasculature, it is well understood that lymphocytes in the blood do pass via these tissues (41,54,55). To determine whether inclusion of the lung/liver vasculature impacts our estimates of the lymphocyte residence time in LNs, we extended the recirculation model (given in Equations 3-9) by adding an equation for lymphocytes in the lung/liver vasculature V (which would be identical to Equation 4). The migration of lymphocytes from the blood to the vasculature was given by rate m BV and exit of lymphocytes from the vasculature into the circulation was given by rate m VB . Then we varied the rates m BV = m VB from 0/h to 30/h and fitted other parameters of the model to the short-term recirculation data. Analysis showed that the estimate of the lymphocyte residence time in the LNs varied relatively little with changes in lymphocyte migration rate via the lung/liver vasculature (from 19.5 to 17.6 h) suggesting that the estimate of the lymphocyte residency time in LNs is robust to exclusion of the lung/liver vasculature from the model. Interestingly, however, increasing the rate of lymphocyte migration via the lung/liver vasculature reduced the quality of the model fit to data (measured as SSR or AIC) which is in line with our result that the models with more than n = 3 tissue compartments describe the cannulation data with poorer quality.
Another important conclusion from our analyses is that the data on lymphocyte dynamics in efferent lymph is not described well by a model in which residence times of lymphocytes in LNs are exponentially distributed ( Table 1). In part, this is because of the wide distribution in the exit rates of labeled lymphocytes in efferent lymph over time. However, describing cell migration via LNs as a simple one directional process and ignoring the ability of lymphocytes to remain in the LN for longer [e.g., by including a "backflow" in cell movement as was done by Thomas et al. (44)] may be an over-simplification. Yet, because the model in which lymphocyte residence times are gamma distributed describes the experimental data with acceptable quality (e.g., see Figure 2), introducing additional details/parameters contradicts the fundamental "Occam's razor" principle.
Our analysis suggests difficulty with interpreting data from ovine LN cannulation experiments in which the dynamics of transferred lymphocytes is not tracked in the blood. In particular, we found that estimates of lymphocyte residence times in LNs do depend on the assumed model for lymphocyte dynamics  [16][17][18][19][20] were fit to dataset #3 from Reynolds et al. (20) assuming that difference in lymphocyte dynamics is due to differential entry of lymphocytes into pLN and mLN. (C,D) Models fit assuming that there is differential residence time of lymphocytes in the LNs. The number of cells exiting pLN (boxes) or mLN (triangles) was recorded for lymphocytes collected originally from pLN (A,C) or from mLN (B,D). Model parameters are shown in Table S5. The average residence times of lymphocytes in pLN (RT p ) and mLN (RT m ) predicted by the best fit model are shown.
in the blood (e.g., single vs. double exponential decline and slow or rapid decline). Therefore, future studies on lymphocyte recirculation kinetics in sheep should always attempt to measure and report concentration and total numbers of transferred cells in the blood.
One of the fundamental questions of lymphocyte recirculation is whether lymphocytes in the blood have some "memory" of the specific LN they recently came from, and if such memory exists, whether it comes from preferential entry into a specific LN or from preferential retention in the LN. Several experimental studies have addressed the question qualitatively. For example, activated lymphocytes, or lymphoblasts, collected from the intestinal lymph of sheep were shown to accumulate preferentially in tissues associated with the gut (56), and a similar finding was reported for lymphoblasts isolated from intestinal lymph of rats (57,58). In contrast, lymphoblasts isolated from peripheral lymph preferentially accumulated in peripheral lymph nodes (57). There has also been a distinction in migratory preference based on cellular subset as it has been observed that small lymphocytes accumulate in mucosal sites such as Peyer's patches (59,60).
We used mathematical modeling to investigate whether preferential accumulation of lymphocytes in the LN of their origin is due to preferential entry or preferential retention for one specific dataset (20). Our analysis showed that a model with preferential retention was not able to accurately describe the experimental data, while the model in which cells could preferentially enter a LN was able to describe the data well (Figure 7). Intuitively, this may be because the earliest increase in the number of cells found in the efferent lymph seems to be driven by rate of cell entry into the node and the data clearly indicate difference in cell accumulation in the efferent lymph depending on the cell's origin. While we have not addressed it formally, it is possible that the overall distribution of residence times (determined by the parameter k) may be different in LNs of different types.
There are a number of limitations with experimental data and our modeling analyses that need to be highlighted. In particular, in all of our experimental data, the dynamics of labeled lymphocytes in the efferent lymph was reported as a frequency of total cells, which required recalculation to determine the total number of cells exiting a specific LN per unit of time [e.g., (42)]. Similarly, calculation of the total number of lymphocytes in the blood requires the knowledge of the total blood volume of animals which was not reported. The required recalculations may introduce errors (e.g., due to incorrectly assumed blood volume in the animals) and thus may influence the values for some estimated model parameters. For example, a smaller assumed blood volume in animals would naturally lead to a lower number of transferred lymphocytes in Frost et al. (42) experiments detected in the blood which should directly impact the estimate of the rate of lymphocyte migration from the blood to the LN. Thus, the absolute values of estimated rates at which lymphocytes are predicted to migrate to LNs from the blood should be treated as approximate. However, estimates of lymphocyte residency time in LNs should be robust to changes in the scaling of lymphocyte numbers in the lymph.
One of the major assumptions we made in the models was that all labeled cells have identical migratory characteristics, e.g., all cells are capable of entering and exiting the LNs and do so at the same rates. In many previous studies the types of lymphocytes used in recirculation experiments (e.g., naive or memory lymphocytes, B or T cells) were not specified, and it is very well possible that migratory properties vary by cell type [e.g., see (9)]. It is clear however that including multiple cell subpopulations will increase complexity of the models, making them unidentifiable from the data we have available. Also, the models based on kinetically homogeneous cell populations could describe the data reasonably well, which suggests that there is no need to introduce a more complex model for such data. Yet, comparison of predictions found by the best fit models for the data did indicate some discrepancy, for example, the best fit model was not fully capable of capturing the peak of the exit rate of labeled lymphocytes in efferent lymph (e.g., Figure 5B). It is possible that including slow and fast recirculating cell sub-populations may be able to fully capture the peak in labeled cells even though we were not able to improve fit of these data by extending the model to two sub-populations with different migration kinetics (results not shown). Additional data that includes variability in lymphocyte dynamics in efferent lymph between different animals may be useful to further show the need for more complex models. Another major assumption of our modeling approach is that lymphocytes in circulation enter LNs via HEVs and not via afferent lymph of the tissues. While there is some experimental support for this assumption for lymphocytes migrating in noninflammatory conditions (2), it is clear that during inflammation in the skin, many cells may enter the skin-draining LNs via afferent lymph (19).
At its core, the combination of experimental data and our models allowed us to estimate the time it takes for lymphocytes to migrate from the blood to the efferent lymph of specific LNsand our results suggest that this time is gamma distributed. For lymphocytes migrating to LNs via HEVs, this distribution is likely to be related to lymphocyte residence time in LNs as lymphocytes pass via HEV rather quickly (53). However, if lymphocytes migrate from the blood to efferent lymph by first entering non-lymphoid tissues (e.g., skin), then exiting the tissue into afferent lymphatics, and then passing via the LN-then our estimates of the average residence time of lymphocytes in LNs are upper bound values. It is also possible that the rate at which lymphocytes leave the final, k th sub-compartment in the LN may be different from that for other sub-compartments. However, there was no need to increase model complexity as the model with a constant rate m LB was sufficient to accurately describe the data.
In most of our models we ignored the possibility of cell death. When describing labeled lymphocyte dynamics during short-term (< 90 h) migration experiments with the recirculation model (Figure 2) we found the need to have a tissue compartment which acts as a sink and thus may represent a death process ( Table 2). However, there appears to be an equilibrium reached by recirculating lymphocytes in the blood by 120 h of the experiment ( Figure 3A) suggesting a limited role of death process in determining overall dynamics of labeled lymphocytes. Still, we performed some additional analyses by adding death rate to all tissue compartments and found that the best fit is found when such death rate is small or non-existent (Table S6). It is possible, however, that the early loss of lymphocytes in the blood as was observed in Frost et al. (42) data (Figure 2A) may be due cell death. It is important to highlight that high rates of cell death may influence interpretation of the data and estimates of the model parameters so future studies should attempt to quantify total cell numbers in as many tissues as possible.
Even with all limitations in the data and assumptions of the models, we provided a quantitative framework to analyze data from LN cannulation experiments in sheep. The models, developed in the paper, may need to be tailored to explain kinetics of lymphocyte recirculation in specific experiments. As illustrated in this work, greater insights into mechanisms regulating lymphocyte migration in large animals such as sheep and humans may thus be obtained by combining the use of quantitative experiments and mathematical modeling.

DATA AVAILABILITY
All datasets analyzed for this study are included in the Supplementary Files.

AUTHOR CONTRIBUTIONS
The study was originally designed by VG. Data were digitized from cited publications by MM. All major analyses were performed by MM. The paper was written by MM and VG.