- Open Access
Characteristic pancreatic and splenic immune cell infiltration patterns in mouse acute pancreatitis
Cell & Bioscience volume 11, Article number: 28 (2021)
A systemic evaluation of immune cell infiltration patterns in experimental acute pancreatitis (AP) is lacking. Using multi-dimensional flow cytometry, this study profiled infiltrating immune cell types in multiple AP mouse models.
Three AP models were generated in C57BL/6 mice via cerulein (CAE) injection, alcohol and palmitoleic acid (EtOH + POA) injection, and alcohol diet feeding and cerulein (EtOH + CAE) injection. Primary pancreatic cells and splenocytes were prepared, and multi-dimensional flow cytometry was performed and analyzed by manual gating and computerized PhenoGraph, followed by visualization with t-distributed stochastic neighbor embedding (t-SNE).
CAE treatment induced a time-dependent increase of major innate immune cells and a decrease of follicular B cells, and TCD4+ cells and the subtypes in the pancreas, whereas elicited a reversed pattern in the spleen. EtOH + POA treatment resulted in weaker effects than CAE treatment. EtOH feeding enhanced CAE-induced amylase secretion, but unexpectedly attenuated CAE-induced immune cell regulation. In comparison with manual gating analysis, computerized analysis demonstrated a remarkable time efficiency and reproducibility on the innate immune cells and B cells.
The reverse pattern of increased innate and decreased adaptive immune cells was consistent in the pancreas in CAE and EtOH + POA treatments. Alcohol feeding opposed the CAE effect on immune cell regulation. Together, the immune profiling approach utilized in this study provides a better understanding of overall immune responses in AP, which may facilitate the identification of intervention windows and new therapeutic strategies. Computerized analysis is superior to manual gating by dramatically reducing analysis time.
Acute pancreatitis (AP) is an inflammatory disorder associated with substantial morbidity and mortality . Alcohol consumption serves as a major etiologic factor in AP, with alcohol-induced or alcoholic pancreatitis accounting for up to 50% of AP cases [2,3,4]. While the majority of AP cases resolve without significant sequelae, approximately 16 % progress to chronic pancreatitis (CP), primarily in patients with strong risk factors such as alcohol abuse [5,6,7]. CP is the number one cause of type 3c diabetes  and a major risk factor for pancreatic cancer, one of the most lethal cancer types . Unfortunately, the standard care for pancreatitis lacks specific pharmacological therapies but remains primarily supportive.
The pathogenesis of AP is thought to originate from sterile autodigestion , with subsequent inflammatory responses playing a critical role in the disease [11, 12]. Previous studies have shown that infiltrating immune cells, particularly innate immune cells such as neutrophils, macrophages/monocytes, and dendritic cells, are elevated during AP development [13,14,15,16,17], while adaptive or acquired immune cells, including T and B cells, have been reported to be differentially regulated [18,19,20]. However, previously published works on this topic in experimental AP mouse models used pooled samples from multiple mouse pancreata in order to obtain sufficient cell counts for the measurement of specific immune cells [13, 21]. For these reasons, a comprehensive profiling of the infiltrating immune cells during AP is currently lacking and therefore desired. The outcomes of such a study will lead to a better understanding of the overall immune responses in AP, as well as identify key infiltrating immune cells for the intervention windows and the development of innovative therapeutic strategies, ultimately improving outcomes of AP patients and reducing the risk of developing CP, diabetes, and pancreatic cancer.
For over 50 years, flow cytometry has been used to identify populations of cells including immune cells . Although long-established, this traditional technique has several limitations. For instance, only a few fluorescent-labeled antibodies can be used for detection and the required manual gating analysis is time-consuming with unavoidable biases among individual investigators . To circumvent these limitations, multi-dimensional flow cytometry, with 8–12 fluorescent-labeled antibodies, was developed . Later on, the expansion of computerized analysis led to reduced workload and bias that inherently encompassed the prior manual gating technique, as well as better visualization. For instance, t-stochastic neighbor embedding (t-SNE) emerged, providing a better way to visualize multi-dimensional cytometry data through nonlinear dimensionality reduction according to the similarities of cells, creating an easily understandable 2D or 3D scatter plot [25, 26]. PhenoGraph was then developed with the ability to cluster cells from high dimensional data by converting single-cell data into a graph representing the connections and relationships between the cells, displaying the data on a t-SNE map . Thus, the application of PhenoGraph and t-SNE in multi-dimensional flow cytometry provide powerful tools for a deeper understanding of cell populations and the potential to replace manual gating .
In this study, we utilized multi-dimensional flow cytometry to profile infiltrating immune cell types in the pancreas and spleen from three AP mouse models that recapitulate a range of AP pathogenesis. We then analyzed and compared the results by manual gating and computerized analysis. Our study identified characteristic regulation of immune cells in the pancreas and spleen in several experiment AP models, and revealed remarkable time-efficiency of the computerized analysis.
Characteristic increase of innate immune cell types but decrease of adaptive immune cell types in the pancreas in cerulein (CAE)-induced acute pancreatitis (AP)
To assess the regulation of infiltrating immune cells during AP development, we utilized the widely used and well-defined CAE mouse model in C57BL/6 mice. AP was confirmed by serum amylase level (Fig. 1a) and histopathological scores (Fig. 1b), both increased at 16 h and declined to the control levels at 30 h onward to day 7.
To profile the infiltrating immune cell types in the pancreas during AP development, we isolated pancreatic immune cells from the mice and performed multi-dimensional flow cytometry using two mouse immune phenotype panels (Table 1). By traditional manual gating, we found that the infiltrating immune cells increased during AP development were primarily innate immune cells (Fig. 1c). Compared to the baseline controls (CON), neutrophils (CD11b + Ly6G+) increased by 14.54-fold at 16 h and 4.60-fold at 30 h (p < 0.05), monocytes (or macrophages, CD11b + Ly6C+) increased by 18.36-fold at 16 h (p < 0.05), and eosinophils (CD11b + Ly6G-Ly6C-SSChigh) increased by 13.08-fold at 16 h and 4.16-fold at 30 h (p < 0.05). Conventional dendritic cells (cDCs, IA/IE + CD19-CD11c+) and red pulp macrophages (RP macs, CD11b+/-IA/IE + CD11c-) increased at 16 h (2.45- and 1.92-fold, respectively, p < 0.05) and at 30 h (2.27- and 2.19-fold, respectively, p < 0.05). All cell types returned to levels similar to those of CON by day 7, except for RP macs, which were still elevated at day 7.
Regarding the changes in adaptive immune cells, compared to CON, total B cells (IA/IE + CD19+) did not change significantly, but follicular B cells (IA/IE + CD19 + CD23 + CD21/35high) decreased 0.65-fold at 16 h (Fig. 1d, p < 0.05). Total T cells (CD5+) did not change significantly. TCD4+, regulatory TCD4+CD25+, naïve TCD4+CD25−, effector memory TCD4+CD62L−CD44+, and naïve TCD8+CD62L+CD44− cells decreased at 16 h (0.69-, 0.70-, 0.69-, 0.62-, and 0.70-fold, respectively, p < 0.05), and returned to the levels of CON at 30 h (except for naïve TCD8+CD62L+CD44− cells) and day 7 (except for effector memory TCD4+CD62L−CD44+ cells). Double negative TCD4−CD8− cells increased at 30 h and persisted to day 7 (1.22- and 1.18-fold, p < 0.05, Fig. 1e).
To better visualize the multi-dimensional cytometry results, the data from each mouse for each time point were organized by t-SNE, which forms a 2D map showing clusters of the specific cell populations, for the innate immune cells and the adaptive immune cells, B cells, from Panel 1 (Fig. 1f), and for the adaptive immune cells, T cells, from Panel 2 (Fig. 1g).
We further conducted computerized analysis directly on the CD45 + cells from the multi-dimensional flow data. PhenoGraph was applied and a heatmap was generated based on the surface marker expression. Cell populations were designated according to the combination of the markers. The regulation pattern of the innate immune cells and B cells generated by PhenoGraph was consistent with that by manual gating (Fig. 2a, b). Further, the regulation pattern of T subtypes generated by PhenoGraph was consistent with that by manual gating for TCD4+, TCD8+ and, TCD4−CD8− cells, but different for total TCD5+ cells. The regulation of TCD4+CD25+, TCD4+CD25−, TCD4+CD62L−CD44+, and TCD8+CD62L+CD44− cells was not identified by PhenoGraph analysis. Additional regulated T cell subtypes, TCD4+CD62L−CD44−, TCD4+CD8+, TCD4+CD8+CD62L−CD44−, and TCD4−CD8−CD62L−CD44− cells were identified by PhenoGraph analysis (Fig. 2c). The data were visualized by t-SNE maps from Panel 1 (Fig. 2d) and Panel 2 (Fig. 2e). Overall, the computerized approach took significantly less time and effort that was spent on the manual gating method.
Differentially regulated innate and adaptive immune cells in the spleen, compared to the immune alterations in the pancreas in CAE-induced AP
The inflammatory profiles of the spleen were also investigated. Compared to the controls, out of innate immune cells, only monocytes decreased (0.60-fold, p < 0.05), while natural killer T (NKT) cells (CD161 + CD5+) increased at 16 h (1.31-fold, p < 0.05) and returned to the control levels at 30 h and day 7 (Fig. 3a). Regarding the adaptive immune cells, B cells did not change (Fig. 3b). Regulatory TCD4+CD25+, TCD8+, and naïve TCD8+CD62L+CD44− cells increased, while TCD4+CD62L−CD44−, effector memory TCD8+CD62L−CD44+, TCD8+CD62L−CD44−, and double positive TCD4+CD8+ cells decreased (p < 0.05, Fig. 3c). t-SNE was conducted to generate a 2D map for better visualization from Panel 1 (Fig. 3d) and Panel 2 (Fig. 3e).
Computerized analysis was once again conducted on the CD45 + cells from the multi-dimensional flow data of the spleen, and revealing similar observations for monocytes (Fig. 4a) and B cells (Fig. 4b), but different for regulations on NKT cells (Fig. 4a) and T cell subtypes (Fig. 4c). t-SNE was conducted from Panel 1 (Fig. 4d) and Panel 2 (Fig. 4e).
Contrasting infiltrating immune cell types in the pancreas and the spleen of alcohol-induced AP models
Alcohol consumption is a major etiologic factor of AP. However, a critical lack of access to human pancreatic tissue has made it difficult to obtain a full understanding of the pathogenic immune mechanisms involved . Therefore, several alcohol animal models have been established, recapitulating human alcoholic pancreatitis for the study of the disease pathogenesis [29,30,31]. To investigate the regulation of infiltrating inflammatory cells in alcohol-induced AP, we took a similar approach as described in the above-mentioned CAE-induced AP model, and profiled the infiltrating inflammatory cells at 16 h in two alcohol-induced AP models, EtOH + POA and EtOH + CAE.
In the EtOH + POA model (Fig. 5), EtOH and POA injections induced higher levels of amylase secretion (Fig. 5a), but similar histopathological scores (Fig. 5b), compared with the controls. Manual gating of the multi-dimensional flow cytometry revealed a similar pattern of innate immune cells and T cells in the pancreas at 16 h, although weaker than that observed in the CAE model (Fig. 5c); neutrophils increased and double positive TCD4+CD8+ cells decreased in the spleen (Fig. 5d).
In the EtOH + CAE model (Fig. 6), mice fed with EtOH liquid diet (dEtOH) for 2 weeks gained less bodyweight than the mice pair-fed with control diet (dCON, Fig. 6a). CAE injections (3 hourly) combined with EtOH feeding induced higher amylase secretion compared to the baseline controls (CON) and dCON feeding groups (Fig. 6b, *p < 0.05), suggesting that EtOH feeding may synergize the CAE effect on amylase secretion. However, the increase of histopathological scores observed for the combined CAE injection and EtOH feeding group did not reach significant levels (Fig. 6c). In the pancreas, manual gating of the multi-dimensional flow cytometry demonstrated increased monocytes, decreased B cells, and decreased T cells and their respective subtypes in the dCON + CAE group, the effect of which was attenuated in the dEtOH + CAE group (Fig. 6d). In the spleen, increased eosinophils were observed in the dEtOH + CAE group. Increased monocytes and increased regulatory TCD4+CD25+, effector memory TCD4+CD62L−CD44+, effector memory TCD8+CD62L−CD44+, TCD8+CD62L−CD44−, and TCD4+CD8+CD62L−CD44+ cells were observed in the dCON + CAE group, the effect of which was attenuated in the dEtOH + CAE group, except for TCD8+CD62L−CD44− cells (Fig. 6e).
In this study, we demonstrated a complex and variable regulation of innate and adaptive immune cells in the pancreas and spleen during AP progression. Specifically, we found that in the pancreas, several innate immune cells increased during the acute phase (16 h) and then declined to baseline levels at recovery phases (30 h–day 7), while adaptive immune cells decreased during the acute phase and recovered to baseline levels subsequently. In the spleen, different patterns were observed compared to those in the pancreas, in that most innate immune cells populations did not change in response to AP induction, particularly T cell subtypes were differentially regulated.
Previous studies investigating the immune cell regulation in AP report individual or a limited number of immune cell types independently in human pancreas studies [14, 19], animal CAE model [13, 16, 32] and other animal models [16,17,18, 20, 33]. Regarding innate immune cell types in literature, neutrophils have been reported to increase in the early phase of AP , while monocytes/macrophages increase in the early and late phases of AP [13, 14, 17, 33]. Increased eosinophils have been detected in human eosinophilic pancreatitis [34, 35], the CAE-induced CP model , and the partial duct ligation mouse AP model . DCs have been shown to increase after AP induction [16, 17]. Our data, consistently and systematically, demonstrated the increase of neutrophils, monocytes/macrophages, eosinophils, and cDCs at the early phase of AP and the recovery at day 7 in the CAE-AP model (Fig. 1a). NK cells and NKT cells were reported to be increased in the pancreas in adenovirus-mediated mouse AP , however we did not observe an increase of NK cells and NKT cells in the pancreas. This variance could be related to pancreatitis induction methods. RP macs are named due to being principally localized in the red pulp of the spleen. They are necessary for maintaining blood homeostasis by phagocytosis of injured and senescent erythrocytes and blood-borne particulates . However, the corresponding population of RP macs in the pancreas is unclear. Here we report an increase of RP macs in the pancreas during AP. Therefore, further study on this special population of cells in the pancreas is warranted. Regarding adaptive immune cells, we observed decreased B cell subtypes in AP, consistent with literature reported in human pancreas study . Variable regulation patterns of T cells and the subtypes have been documented in human pancreas study [14, 19] and animal models [16, 33]. We observed a trend of decreased T cells and the subtypes in AP, and also identified the regulation of several additional T cell subtypes. Overall, our data revealed that the immune cell regulation persisted after serum amylase levels and histopathological changes of the pancreas recovered to baseline levels.
In the spleen, we observed decreased monocytes and increased NKT cells at 16 h, consistent with literature reports [20, 33]. Furthermore, increased TCD4+ cells were reported in the blood, liver, and spleen of sodium-taurocholate-induced necrotizing acute pancreatitis in mice . We did not observe increased TCD4+ cells, but demonstrated differentially regulated T cell subtypes in the spleen of the CAE-induced AP model.
Multi-dimensional flow cytometry serves as an effective systems biology tool to adequately evaluate complex cellular systems and clarify cellular identity and function . Although supplying a more robust analysis with increased fluorescent-labeled antibody capabilities, the manual gating analysis can be quite taxing on the individual researcher and is inherently prone to bias. Newly developed computerized analysis techniques, such as PhenoGraph and t-SNE, dramatically reduce the time to analyze data while providing parallel, if not greater, quality. Our study demonstrated that computerized analysis techniques took a fraction of the time for analysis and achieved comparable results for the identification of innate immune cells and B cells as that by manual gating. However, further optimization on analyzing T cell subtypes is warranted.
Comparing computerized analysis with manual gating in this study, our study supports the use of computerized analysis techniques directly without manual gating in further studies for innate immune cells and B cells. As high-dimensional flow cytometry (18 to 30 parameters) and mass cytometry (> 40 parameters) are the future direction in the field of immunology research [39, 40], for instance, combining the two mouse immune panels used in the current study would further improve efficiency.
Since there are fewer studies on the infiltrating immune cells in alcohol-induced AP animal models, we performed immune profiling in two alcohol-induced AP models to explore this aspect. We observed a similar immune profile pattern although with a weaker effect, in the pancreas of the EtOH + POA model compared to the CAE model. In the EtOH + CAE model, we only observed increased monocytes out of the tested innate immune cells and decreased B and T cell subtypes in the pancreas in the dCON + CAE group. Interestingly, most of the effects were attenuated in the dEtOH + CAE group (Fig. 6). The opposing effects of EtOH feeding on CAE-induced immune responses appear to be contradicting to its enhancing effect on CAE-induced amylase secretion. Further investigation is warranted to uncover the underlying mechanism.
There are several aspects of novelty in this study. First, we utilized three AP models to recapitulate a range of AP pathogenesis in order to demonstrate a “big picture” of immune cell infiltration and regulation patterns in the pancreas and spleen. Second, we demonstrated a systemic analysis of the immune cell regulation in AP. Third, we discovered an unexpected opposing effect of EtOH on CAE-induced immune cell regulation. In addition, we validated the time efficiency and data reproducibility using computerized analysis by directly comparing with traditional manual gating analysis. There are several limitations in this study. For feasibility reasons, only pancreas and spleen were analyzed in current studies. Peripheral blood and other co-morbid organs should be included, as they have been independently studied in clinical research and animal studies. The alcohol models utilized in this study only induced amylase secretion, with no significant pancreatic histopathological damages, resulting in weaker immune responses than that observed in the CAE model. Therefore, alcohol-induced animal AP models should be optimized to achieve stronger AP induction. Although certain inconsistent patterns of the immune cell regulation were derived from three different AP animal models, upregulation of innate immune cells and downregulation of adaptive immune cells in the pancreas are consistent across these different models, and, importantly, validate the findings derived from clinical studies. Although reproducibility of each experiment is fair, it is still necessary to include baseline controls for each repeat experiment in order to normalize and combine the results within the same treatment. Use of untreated mice as baseline controls is an acceptable option in the pancreas research field, particularly in comparison of multiple time points post-injury . Taking a similar approach, in this study, untreated mice in individual experiments were set as baseline controls in an effort to compare the time-dependent responses post-AP induction, also, to keep consistency to compare magnitude of the responses among different AP models.
In conclusion, by using multi-dimension flow cytometry, we profiled infiltrating immune cell types in the pancreas and spleen of CAE- and alcohol-induced AP mouse models. Specifically, we found that the majority of infiltrating immune cells in the AP pancreas were innate immune cells, accompanied with downregulated adaptive immune cells. We observed differential regulation patterns in the spleen, where monocyte/macrophages decreased and NKT cells increased, associated with differentially regulated T cell subtypes. The EtOH models demonstrated similar but weaker trends on infiltrating immune cell types when compared to the CAE model. The results in this study further attested that computerized analysis techniques are superior to manual gating in ease and time efficiency, and data reproducibility particularly for innate immune cells and B cells. Overall, the results from this comprehensive profiling demonstrate the time-dependent immune responses in AP, which may lead to identification of cellular targets and time windows for further studies and therapeutic advancements.
Cerulein (CAE), the decapeptide analog of the potent pancreatic secretagogue cholecystokinin, was purchased from Bachem Americas, Inc. (Torrance, CA); Phadebas Amylase Test from Magle Life Sciences (Cambridge, MA); alcohol (EtOH, 200 proof, 100%) from Fisher Scientific (Hampton, NH); palmitoleic acid (POA, a non-oxidative ethanol metabolite) from Sigma-Aldrich (St. Louis, MO); and rodent liquid diets from Bio-Serv (Flemington, NJ, Catalog # F1258SP for EtOH feeding, Catalog # F1259SP for pair-fed control).
Animals and AP induction
All animal experiments were performed according to the protocols approved by the Animal Welfare Committee of the University of Texas Health Science Center at Houston. C57BL/6 mice (6–8 weeks old, male) were purchased from Jackson Laboratories, Inc. (Bar Harbor, ME). The animals were housed in a climate-controlled room with an ambient temperature of 23 °C and a 12:12-h light-dark cycle. Animals were fed standard laboratory chow (LabDiet 5053, LabDiet, St. Louis, MO), given water ad libitum, and randomly assigned to control or experimental groups.
Three AP mouse models were generated: (1) CAE model. AP was induced by CAE (50 µg/kg, 9 injections at hourly intervals, intraperitoneal (ip)) as described previously . The mice were euthanized at 16 h, 30 h and day 7 after the last CAE injection (n = 6 mice/group). (2) EtOH-POA model. The mice were injected (ip) with EtOH (1.32 g/kg) followed by POA (4 mg/kg) 1 h apart [29, 30]. The mice were euthanized 16 h after the POA injection (n = 3 mice/group). (3) EtOH + CAE model. The mice were fed with Liebere-DeCarli EtOH liquid diet (5 % v/v EtOH) or pair-fed with control diet for 2 weeks . All mice received CAE (50 µg/kg, 3 injections at hourly intervals, ip) on the last day. The mice were euthanized 16 h after the last CAE injection (n = 6 mice/group). The mouse blood, pancreas, and spleen were harvested for analysis. For each set of experiment, untreated mice were included as baseline controls (n = 3 mice) to keep consistency for comparison between different sets of experiments.
Serum amylase levels were determined using the Phadebeas Amylase test tablets as instructed by the manufacturer and as previously described .
Pancreatic morphological evaluation
Paraffin-embedded pancreas samples were sectioned (5 µm), stained with hematoxylin & eosin (H&E), and examined by an experienced pathologist blinded to the sample identity. Histopathological changes of AP were evaluated according to the criteria from previously published works [42, 43], and scored from 0 to 4 as absence to the most severe injury on the extent of edema, inflammation, and acinar necrosis.
Single cell suspension from the pancreas was prepared as previously described [13, 21]. Briefly, pancreas tissue was minced and digested in digestion buffer (Hank’s balanced salt solution (HBSS) without Mg2+, without Ca2+, containing 8 % fetal bovine serum (FBS), 2 mg/ml collagenase IV, 20 µg/ml DNase I) at 37 °C and shaking for 15 min. The digested tissue was filtered through a 100 µm cell strainer and washed with fluorescence-activated cell sorter (FACS) buffer (HBSS containing 2 % FBS and 4 % bovine serum albumin). The suspended cells were centrifuged at 400g for 5 min. The cell pellet was treated with red blood cell (RBC) lysis buffer and EasyStep dead cell removal kit (StemCell Technologies, Vancouver, Canada), centrifuged at 400g for 5 min, and suspended for flow cytometry.
Single cell suspension from the spleen was prepared by the cytometry and cell sorting core (CCSC) of Baylor College of Medicine (Houston, TX). Briefly, spleen tissue was dissociated and digested in digestion buffer (RPMI 1640 containing 2 mg/ml DNase I, 10 mg/ml collagenase II) in the tubes, placed on gentleMACS Dissociator (Miltenyi Biotec, Bergisch Gladbach, Germany), and incubated at room temperature for 10 min, twice. The stopping buffer (1x PBS, 0.1 M Ethylenediaminetetraacetic acid (EDTA)) was added to the digestion mix, which was then filtered through a mesh filter cap and centrifuged at 400g for 5 min. The cell pellet was suspended in FACs buffer (1x phosphate-buffered saline (PBS) without Mg2+, without Ca2+, 2 % FBS, 25 mM HEPES, 2 mM EDTA) and incubated with RBC Lysis buffer, washed again and suspended in FACS buffer for flow cytometry.
Flow cytometry assays
All flow cytometry assays were performed by CCSC of Baylor College of Medicine. Briefly, the isolated cells were incubated with the block solution (FACS buffer with 1:125 dilution of CD16/CD32), and then stained with the mouse immune phenotype Panel 1 or 2 antibody cocktail (Table 1) at 4 °C for 20 min in dark. The cells isolated from the spleen were also used as the Fluorescence Minus One (FMO) control. The cells were centrifuged at 400g for 5 min, and the cell pellet was resuspended in running buffer (FACS Buffer with one drop/ml Nuc Blue Fixed DAPI Stain. Life Technologies, Carlsbad, CA). The flow cytometry reaction was run on LSRII Cell Analyzers (BD Biosciences, San Jose, CA); 500,000 live singlet events for Panel 1 and 350,000 live singlet events for Panel 2 were recorded.
Cytometry data were analyzed with Flowjo 10 (Treestar Inc, Ashland, OR). All FCS files were cleaned by Flowjo plugin FlowAI . Conventional manual gating was first conducted for identification of cell populations and statistical analysis. Each cell type/population with mean event number more than 100 in any group were adopted to have coefficient of variation ≥ 10 %  and analyzed as % of total leukocytes (CD45+). The results of manual gating were visualized in 2D t-SNE maps. Briefly, the CD45 + cells of all samples were downsampled to 3000 (or all for events less than 3000) and concatenated. t-SNE was run on the global concatenated data with iteration of 1000, perplexity of 30, learning rate of 200 and theta of 0.5 to obtain a 2D map by flowjo plugin t-SNE . Group gating by sample ID was conducted on the global concatenated data to obtain group concatenated data as well as t-SNE maps for each group. Cell populations of each group were superposed on t-SNE maps respectively. For computerized analysis, PhenoGraph was run on the concatenated data with the parameter k = 30 to define subpopulations by Flowjo plugin PhenoGraph , and the clusters generated by PhenoGraph were then applied to each group and superposed onto t-SNE maps. A heatmap of surface marker expression was generated by Flowjo plugin ClusterExplorer for interpretation of the cell populations clustered by PhenoGraph.
Statistics were conducted with GraphPad Prism 8 (GraphPad Software. La Jolla, CA) and SPSS 24.0 (IBM Corp., New York, NY). Data are presented as mean ± standard error (SE). One-way ANOVA followed by LSD post hoc analysis was used to determine the differences among multiple groups. Student’s t test was used to compare the difference between two groups. p value less than 0.05 was considered significant.
Availability of supporting data
All data generated or analyzed during this study are included in this article.
t-distributed stochastic neighbor embedding
- NK cells:
Natural killer cells
Conventional dendritic cells
- RP macs:
Red pulp macrophages
Frossard JL, Steer ML, Pastor CM. Acute pancreatitis. Lancet. 2008;371(9607):143–52.
Go V, Everhart J, editors. Pancreatitis. Washington: US Government Printing Office; 1994.
Herreros-Villanueva M, Hijona E, Banales JM, Cosme A, Bujanda L. Alcohol consumption on pancreatic diseases. World J Gastroenterol. 2013;19(5):638–47.
Murtaugh LC, Keefe MD. Regeneration and repair of the exocrine pancreas. Annu Rev Physiol. 2015;77:229–49.
Yadav D, Whitcomb DC. The role of alcohol and smoking in pancreatitis. Nat Rev Gastroenterol Hepatol. 2010;7(3):131–45.
Lankisch PG, Breuer N, Bruns A, Weber-Dany B, Lowenfels AB, Maisonneuve P. Natural history of acute pancreatitis: a long-term population-based study. Am J Gastroenterol. 2009;104(11):2797–805.
Kleeff J, Whitcomb DC, Shimosegawa T, Esposito I, Lerch MM, Gress T, et al. Chronic pancreatitis. Nat Rev Dis Primers. 2017;3:17060.
Hart PA, Bellin MD, Andersen DK, Bradley D, Cruz-Monserrate Z, Forsmark CE, et al. Type 3c (pancreatogenic) diabetes mellitus secondary to chronic pancreatitis and pancreatic cancer. Lancet Gastroenterol Hepatol. 2016;1(3):226–37.
Ryan DP, Hong TS, Bardeesy N. Pancreatic adenocarcinoma. N Engl J Med. 2014;371(11):1039–49.
Saluja AK, Lerch MM, Phillips PA, Dudeja V. Why does pancreatic overstimulation cause pancreatitis? Annu Rev Physiol. 2007;69:249–69.
Al Mofleh IA. Severe acute pancreatitis: pathogenetic aspects and prognostic factors. World J Gastroenterol. 2008;14(5):675–84.
Jensen JN, Cameron E, Garay MV, Starkey TW, Gianani R, Jensen J. Recapitulation of elements of embryonic development in adult mouse pancreatic regeneration. Gastroenterology. 2005;128(3):728–41.
Zhao Q, Wei Y, Pandol SJ, Li L, Habtezion A. STING signaling promotes inflammation in experimental acute pancreatitis. Gastroenterology. 2018;154(6):1822–35.
Bhatnagar A, Wig JD, Majumdar S. Expression of activation, adhesion molecules and intracellular cytokines in acute pancreatitis. Immunol Lett. 2001;77(3):133–41.
Hatano N, Sugiyama M, Watanabe T, Atomi Y. Opsonin receptor expression on peritoneal exudative and circulatory neutrophils in murine acute pancreatitis. Pancreas. 2001;23(1):55–61.
Bedrosian AS, Nguyen AH, Hackman M, Connolly MK, Malhotra A, Ibrahim J, et al. Dendritic cells promote pancreatic viability in mice with acute pancreatitis. Gastroenterology. 2011;141(5):1915-26 e1–14.
Van Gassen N, Van Overmeire E, Leuckx G, Heremans Y, De Groef S, Cai Y, et al. Macrophage dynamics are regulated by local macrophage proliferation and monocyte recruitment in injured pancreas. Eur J Immunol. 2015;45(5):1482–93.
Schmidt AI, Kuhlbrey C, Lauch R, Wolff-Vorbeck G, Chikhladze S, Hopt UT, et al. The predominance of a naive T helper cell subset in the immune response of experimental acute pancreatitis. Pancreatology. 2017;17(2):209–18.
Caluianu EI, Alexandru DO, Tartea EA, Georgescu M, Mercu TD, Trasca ET, et al. Assessment of T- and B-lymphocytes and VEGF-A in acute pancreatitis. Roman J Morphol Embryol. 2017;58(2):481–6.
Schmidt AI, Seifert GJ, Lauch R, Wolff-Vorbeck G, Chikhladze S, Hopt UT, et al. Organ-specific monocyte activation in necrotizing pancreatitis in mice. J Surg Res. 2015;197(2):374–81.
Xue J, Sharma V, Hsieh MH, Chawla A, Murali R, Pandol SJ, et al. Alternatively activated macrophages promote pancreatic fibrosis in chronic pancreatitis. Nat Commun. 2015;6:7158.
Fulwyler MJ. Electronic separation of biological cells by volume. Science. 1965;150(3698):910–1.
Mair F, Hartmann FJ, Mrdjen D, Tosevski V, Krieg C, Becher B. The end of gating? An introduction to automated analysis of high dimensional cytometry data. Eur J Immunol. 2016;46(1):34–43.
Roederer M, De Rosa S, Gerstein R, Anderson M, Bigos M, Stovel R, et al. 8 color, 10-parameter flow cytometry to elucidate complex leukocyte heterogeneity. Cytometry. 1997;29(4):328–39.
Maaten. Lvd HG. Visualizing Data using t-SNE. J Mach Learn Res. 2008;9:2579–605.
Amir el AD, Davis KL, Tadmor MD, Simonds EF, Levine JH, Bendall SC, et al. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat Biotechnol. 2013;31(6):545–52.
Levine JH, Simonds EF, Bendall SC, Davis KL, Amir el AD, Tadmor MD, et al. Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis. Cell. 2015;162(1):184–97.
Habtezion A. Inflammation in acute and chronic pancreatitis. Curr Opin Gastroenterol. 2015;31(5):395–9.
Vigna SR, Shahid RA, Liddle RA. Ethanol contributes to neurogenic pancreatitis by activation of TRPV1. 2014;28(2):891–6.
Barlass U, Dutta R, Cheema H, George J, Sareen A, Dixit A, et al. Morphine worsens the severity and prevents pancreatic regeneration in mouse models of acute pancreatitis. Gut. 2018;67(4):600–2.
Ren Z, Yang F, Wang X, Wang Y, Xu M, Frank JA, et al. Chronic plus binge ethanol exposure causes more severe pancreatic injury and inflammation. Toxicol Appl Pharmcol. 2016;308:11–9.
Gukovskaya AS, Vaquero E, Zaninovic V, Gorelick FS, Lusis AJ, Brennan ML, et al. Neutrophils and NADPH oxidase mediate intrapancreatic trypsin activation in murine experimental acute pancreatitis. Gastroenterology. 2002;122(4):974–84.
Shifrin AL, Chirmule N, Gao GP, Wilson JM, Raper SE. Innate immune responses to adenoviral vector-mediated acute pancreatitis. Pancreas. 2005;30(2):122–9.
Abraham SC, Leach S, Yeo CJ, Cameron JL, Murakata LA, Boitnott JK, et al. Eosinophilic pancreatitis and increased eosinophils in the pancreas. Am J Surg Pathol. 2003;27(3):334–42.
Tian L, Fu P, Dong X, Qi J, Zhu H. Eosinophilic pancreatitis: three case reports and literature review. Mol Clin Oncol. 2016;4(4):559–62.
Manohar M, Verma AK, Venkateshaiah SU, Mishra A. Role of eosinophils in the initiation and progression of pancreatitis pathogenesis. Am J Physiol Gastrointest Liver Physiol. 2018;314(2):G211–22.
Kurotaki D, Uede T, Tamura T. Functions and development of red pulp macrophages. Microbiol Immunol. 2015;59(2):55–62.
O’Donnell EA, Ernst DN, Hingorani R. Multiparameter flow cytometry: advances in high resolution analysis. Immune Netw. 2013;13(2):43–54.
Spitzer MH, Nolan GP. Mass cytometry: single cells, many features. Cell. 2016;165(4):780–91.
Weber LM, Robinson MD. Comparison of clustering methods for high-dimensional single-cell flow and mass cytometry data. Cytom Part A. 2016;89(12):1084–96.
Boggs K, Wang T, Orabi AI, Mukherjee A, Eisses JF, Sun T, et al. Pancreatic gene expression during recovery after pancreatitis reveals unique transcriptome profiles. Sci Rep. 2018;8(1):1406.
Cao Y, Yang W, Tyler MA, Gao X, Duan C, Kim SO, et al. Noggin attenuates cerulein-induced acute pancreatitis and impaired autophagy. Pancreas. 2013;42(2):301–7.
Schmidt J, Lewandrowski K, Fernandez-del Castillo C, Mandavilli U, Compton CC, Warshaw AL, et al. Histopathologic correlates of serum amylase activity in acute experimental pancreatitis. Dig Dis Sci. 1992;37(9):1426–33.
Monaco G, Chen H, Poidinger M, Chen J, de Magalhaes JP, Larbi A. flowAI: automatic and interactive anomaly discerning tools for flow cytometry data. Bioinformatics. 2016;32(16):2473–80.
Allan AL, Keeney M. Circulating tumor cell analysis: technical and statistical considerations for application to the clinic. J Oncol. 2010;2010:426218.
Jimenez-Carretero D, Ligos JM, Martínez-López M, Sancho D, Montoya MC. Flow cytometry data preparation guidelines for improved automated phenotypic analysis. J Immunol. 2018;200(10):3319–31.
The authors thank Jiajing Li at Department of Surgery, and Histology Laboratory at Department of Pathology and Laboratory Medicine, UTHSC-Houston for technical support. This project was supported by the Cytometry and Cell Sorting Core at Baylor College of Medicine with funding from the CPRIT Core Facility Support Award (CPRIT-RP180672), the NIH (CA125123 and RR024574), and the assistance of Dr. Paul Porter and Mr. Joel M. Sederstrom.
This study was supported by the National Institute of Health (NIH) grant 1 R21 AA027014-01A1 (T.C.K), Jack H Mayfield M.D. Distinguished Professorship in Surgery (T.C.K), Dean’s fund for Summer Research Program (J.M.D).
Ethical approval and consent to participate
All animal experiments were performed according to the protocols approved by the Animal Welfare Committee of the University of Texas Health Science Center at Houston.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yang, B., Davis, J.M., Gomez, T.H. et al. Characteristic pancreatic and splenic immune cell infiltration patterns in mouse acute pancreatitis. Cell Biosci 11, 28 (2021). https://doi.org/10.1186/s13578-021-00544-1
- Mouse acute pancreatitis
- Immune profiling
- t-distributed stochastic neighbor embedding (t-SNE)