Abstract
Transcriptional and proteomic profiling of individual cells have revolutionized interpretation of biological phenomena by providing cellular landscapes of healthy and diseased tissues1,2. These approaches, however, do not describe dynamic scenarios in which cells continuously change their biochemical properties and downstream ‘behavioural’ outputs3,4,5. Here we used 4D live imaging to record tens to hundreds of morpho-kinetic parameters describing the dynamics of individual leukocytes at sites of active inflammation. By analysing more than 100,000 reconstructions of cell shapes and tracks over time, we obtained behavioural descriptors of individual cells and used these high-dimensional datasets to build behavioural landscapes. These landscapes recognized leukocyte identities in the inflamed skin and trachea, and uncovered a continuum of neutrophil states inside blood vessels, including a large, sessile state that was embraced by the underlying endothelium and associated with pathogenic inflammation. Behavioural screening in 24 mouse mutants identified the kinase Fgr as a driver of this pathogenic state, and interference with Fgr protected mice from inflammatory injury. Thus, behavioural landscapes report distinct properties of dynamic environments at high cellular resolution.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
All data and materials used in the study are available to any researcher for purposes of reproducing or extending these analyses. Source data are provided with this paper.
Code availability
All code used are available to any researcher for purposes of reproducing or extending these analyses. The code for ACME is available at https://doi.org/10.5281/zenodo.5638537.
References
Kwok, I. et al. Combinatorial single-cell analyses of granulocyte–monocyte progenitor heterogeneity reveals an early uni-potent neutrophil progenitor. Immunity 53, 303–318.e5 (2020).
Amit, I., Regev, A. & Hacohen, N. Strategies to discover regulatory circuits of the mammalian immune system. Nat. Rev. Immunol. 11, 873–880 (2011).
Lämmermann, T. et al. Neutrophil swarms require LTB4 and integrins at sites of cell death in vivo. Nature 498, 371–375 (2013).
Sreeramkumar, V. et al. Neutrophils scan for activated platelets to initiate inflammation. Science 346, 1234–1238 (2014).
Woodfin, A. et al. The junctional adhesion molecule JAM-C regulates polarized transendothelial migration of neutrophils in vivo. Nat. Immunol. 12, 761–769 (2011).
Papalexi, E. & Satija, R. Single-cell RNA sequencing to explore immune cell heterogeneity. Nat. Rev. Immunol. 18, 35–45 (2018).
Sumen, C., Mempel, T. R., Bazo, I. B. & von Andrian, U. H. Intravital microscopy: visualizing immunity in context. Immunity 21, 315–329 (2004).
Lelkes, E., Headley, M. B., Thornton, E. E., Looney, M. R. & Krummel, M. F. The spatiotemporal cellular dynamics of lung immunity. Trends Immunol. 35, 379–386 (2014).
Yamada, K. M. & Sixt, M. Mechanisms of 3D cell migration. Nat. Rev. Mol. Cell Biol. 20, 738–752 (2019).
Goh, C. C. et al. Real-time imaging of dendritic cell responses to sterile tissue injury. J. Invest. Dermatol. 135, 1181–1184 (2015).
Goh, C. C. et al. The impact of ischemia–reperfusion injuries on skin resident murine dendritic cells. Eur. J. Immunol. 48, 1014–1019 (2018).
Adrover, J. M. et al. Programmed ‘disarming’ of the neutrophil proteome reduces the magnitude of inflammation. Nat. Immunol. 21, 135–144 (2020).
Adrover, J. M. et al. A neutrophil timer coordinates immune defense and vascular protection. Immunity 50, 390–402.e10 (2019).
Hasenberg, A. et al. Catchup: a mouse model for imaging-based tracking and modulation of neutrophil granulocytes. Nat. Methods 12, 445–452 (2015).
Hidalgo, A. et al. Heterotypic interactions enabled by polarized neutrophil microdomains mediate thromboinflammatory injury. Nat. Med. 15, 384–391 (2009).
Mócsai, A., Walzog, B. & Lowell, C. A. Intracellular signalling during neutrophil recruitment. Cardiovasc. Res. 107, 373–385 (2015).
Weir, M. C. et al. Selective inhibition of the myeloid Src-family kinase Fgr potently suppresses AML cell growth in vitro and in vivo. ACS Chem. Biol. 13, 1551–1559 (2018).
García-Prieto, J. et al. Neutrophil stunning by metoprolol reduces infarct size. Nat. Commun. 8, 14780 (2017).
Hirahashi, J. et al. Mac-1 (CD11b/CD18) links inflammation and thrombosis after glomerular injury. Circulation 120, 1255–1265 (2009).
McArdle, S. et al. Migratory and dancing macrophage subsets in atherosclerotic lesions. Circ. Res. 125, 1038–1051 (2019).
Mempel, T. R., Henrickson, S. E. & von Andrian, U. H. T-cell priming by dendritic cells in lymph nodes occurs in three distinct phases. Nature 427, 154–159 (2004).
Boulch, M. et al. A cross-talk between CAR T cell subsets and the tumor microenvironment is essential for sustained cytotoxic activity. Sci. Immunol. 6, eabd4344 (2021).
Ivanovitch, K., Temiño, S. & Torres, M. Live imaging of heart tube development in mouse reveals alternating phases of cardiac differentiation and morphogenesis. eLife 6, e30668 (2017).
Lindquist, R. L. et al. Visualizing dendritic cell networks in vivo. Nat. Immunol. 5, 1243–1250 (2004).
Faust, N., Varas, F., Kelly, L. M., Heck, S. & Graf, T. Insertion of enhanced green fluorescent protein into the lysozyme gene creates mice with green fluorescent granulocytes and macrophages. Blood 96, 719–726 (2000).
Palomino-Segura, M., Virgilio, T., Morone, D., Pizzagalli, D. U. & Gonzalez, S. F. Imaging cell interaction in tracheal mucosa during influenza virus infection using two-photon intravital microscopy. J. Vis. Exp. 58355 (2018).
Bauer, C. A. et al. Dynamic Treg interactions with intratumoral APCs promote local CTL dysfunction. J. Clin. Invest. 124, 2425–2440 (2014).
Tran Cao, H. S. et al. Development of the transgenic cyan fluorescent protein (CFP)-expressing nude mouse for ‘technicolor’ cancer imaging. J. Cell. Biochem. 107, 328–334 (2009).
Devi, S. et al. Neutrophil mobilization via plerixafor mediated CXCR4 inhibition arises from lung demargination and blockade of neutrophil homing to the bone marrow. J. Exp. Med. 210, 2321–2336 (2013).
Çiçek, Ö., Abdulkadir, A., Lienkamp, S. S., Brox, T. & Ronneberger, O. In MICCAI 2016. Lecture Notes in Computer Science Vol. 9901 (eds Ourselin, S. et al.) (Springer, 2016).
Welch, G. & Bishop, G. An introduction to the Kalman Filter. In Pract. (2006).
Hao, Y. et al. Integrated analysis of multimodal single-cell data. Cell 184, 3573–3587.e29 (2021).
Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S 4th edn (Springer, 2002).
Bates, D., Mächler, M., Bolker, B. M. & Walker, S. C. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Lloyd, S. P. Least squares quantization in PCM. IEEE Trans. Inf. Theory 28, 129–137 (1982).
Henique, C. et al. Genetic and pharmacological inhibition of microRNA-92a maintains podocyte cell cycle quiescence and limits crescentic glomerulonephritis. Nat. Commun. 8, 1829 (2017).
Acknowledgements
We thank all members of the Hidalgo laboratory and M. Desco for discussion; P. Frenette for inspiring this study; C. C. Goh and E. Y. Kim for seeding imaging experiments; the electron microscopy unit from the faculty of Medicine of Universidad Autonoma de Madrid for help with experiments; E. Marín, L. Cabezuela, E. Santos, R. Mota and the animal facility at CNIC for animal husbandry, animal procedures and histology; J. Rossaint, M. Gunzer, J.A. Enriquez, A. Mocsai, R.W. Hendricks, G. Sabio, M. Sperandio, E. Hirsch and B. Walzog for the generous gift of mutant mice; and C. Torroja, D. Jiménez and M. Desco for technical advice. This study was supported by RTI2018-095497-B-I00 from Ministerio de Ciencia e Innovación (MCIN), HR17_00527 from Fundación La Caixa, Transatlantic Network of Excellence (TNE-18CVD04) from the Leducq Foundation, and FET-OPEN (no. 861878) from the European Commission to A.H. M.P-S. is supported by a Federation of European Biochemical Societies and the EMBO ALTF (no. 1142-2020) long-term fellowships. J.S. is supported by a fellowship (PRE2019-089130) from MICINN and A.A.-C. is supported by fellowship CF/BQ/DR19/11740022 from La Caixa Foundation. J.L.Y.L. was supported by A*STAR and a Juan de la Cierva JCI-2017-33136 Fellowship from MICINN. S.D.C. is a recipient of a Marie Sklodowska-Curie fellowship (749731). M.G. is supported by SAF2017-89116R-P from MCIN and HR18_00120 from la Fundación La Caixa. T.R.M. is supported by grant NIH AI163223, L.G.N. is supported by SIgN core funding from A*STAR, and G.F.C. is supported by MCIN/AEI/10.13039/501100011033 (grant PID2019-110895RB-I00) and by Junta de Comunidades de Castilla-La Mancha (SBPLY/19/180501/000211). F.S.-C. is supported by MCIN (grant RTI2018-102084-B-I00), O.S. is supported by the Leducq Foundation (TNE-18CVD04), F.D.-d.-M. is supported by MCIN (TEC2017-84395-P), and T.E.S. is supported by the National Cancer Institute, NIH grant CA233576. The CNIC is supported by the MCIN and the Pro-CNIC Foundation.
Author information
Authors and Affiliations
Contributions
Conceptualization: A.H. Methodology: M.P.-S., F.S.-C., A.Z., T.E.S., P.-L.T., H.C.E.W., O.S., M.M.-M., G.F.C., I.G.-D., F.D.-d.-M. and A.H. Investigation: G.C., M.P.-S., M.M.-M., J.S., D.G.A., G.F.C., J.L.Y.L., R.M., J.M.A., A.A.-C., S.M.-S., A.S.d.V., S.D.C. and M.D.P. Visualization and intravital microscopy: G.C., M.P.-S., M.M.-M., J.S. and D.G.A. Image data analysis: M.P.-S., I.G.-D., M.M.-M., J.L.Y.L. and M.D.P. Myocardial infarction and glomerulonephritis experiments: G.C. and S.M.-S. Funding acquisition: A.H., F.S.-C., P.-L.T., L.G.N. and G.F.C. Supervision: P.-L.T., S.F.G., T.R.M., A.A.-S., L.G.N., G.F.C., I.G.-D., F.D.-d.-M. and A.H. Writing, original draft: M.P.-S. and A.H. Writing, review and editing: All authors. Contribution note: M.M.-M. and J.S. contributed equally as secondary authors in this paper.
Corresponding author
Ethics declarations
Competing interests
A.H. is a consultant for Flagship Pioneering. The other authors declare no competing interests.
Peer review information
Nature thanks Klaus Ley and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Selection of parameters for behavioral analyses.
(A) Application of size filters. Left, representative images showing use of >40 voxel filter to eliminate subcellular objects in the trachea experiments. Comparison of the raw vs. surface reconstructed objects (see merged) eliminates fragment-like objects, as shown in the insets. Right, tSNE representation of the trachea dataset in which the filter for cell size was set to 0 (no filter) or 160 voxels, showing objects with sizes below 40 voxels (threshold used in our experiments). The number of objects for each representation are shown in brackets. Note that segregation of neutrophils and DCs into different visual clusters is compromised in the absence of filter. (B) Workflow for parameter selection. 4D images were analyzed to extract morphometric and kinetic parameters (118 in our experiments using Imaris software). We performed supervised selection of the best characteristics following criteria of redundancy, biological features of interest in the specific biological setup, or removal of non-biological parameters such as arbitrary position. In parallel we generated correlation networks for all parameters and each experiment (118 in the experiments reported here), and we visualized the distribution of the selected parameters in the correlation networks (see below). Finally, we reduced dimensionality using the selected parameters to identify behavioral clusters for further validation. For the “training” experiments shown here, where cell identities were known, we determined LRI/ARI to complement our correlation networks with the power of each parameter to classify cells correctly. The workflow is fully adaptable to other image analysis tools, as well as algorithms to establish correlation between parameters and for dimensional reduction, including elastic net regression methods. (C) Correlation networks used for the experiments shown in Fig. 1 with Imaris image analysis. Networks on the left column highlight the specific 31 parameters selected, which are identified by the number code shown in Supplementary Table 2. Correlation networks on the right column correspond to the three datasets shown in Fig. 1 (influenza infection in the trachea, ischemia-reperfusion in skin and laser injury in skin), showing parameters as circles whose diameters are proportional to their LRI, as well as the positive (red) and negative (blue) correlations between each pair of parameters. The thickness of the links is proportional to the absolute value of the Pearson correlation coefficient for each pair, and the distance between parameters reflect the similarity of the Pearson coefficients with the rest of parameters. (D) Violin plots showing LRI/ARI values for all 118 vs. the selected 31 parameters. Lines represent medians. Data compared by Mann-Whitney non-parametric test. (E) Heatmap showing the LRI/ARI values for each of the 31 selected parameters for each experiment, as well as the geometric mean for the three experiments combined, reflecting the average power of each parameter in our experiments. (F) Quality of parameter selection. Top, number of parameters selected, and predictive power (LRI or ARI in red) of the parameters selected by Lasso regression compared with our list of 31 selected parameters. Bottom, comparison of the distribution of the 25 parameters selected by Lasso regression and our 31 selected parameters across the correlation network for the trachea experiment. (G) tSNE plots generated by considering only the morphometric, or only the kinetic parameters, or both combined. Donut plots show the distribution of the analyzed cell types (macrophages, DCs and neutrophils) in each cluster. Note that the accuracy in identifying specific cell types for each cluster is always highest when both classes of parameters are combined. (H) tSNE plots generated by considering all 118 parameters or only the selected 31 across all three experimental setups (influenza infection in the trachea, ischemia-reperfusion in skin and laser injury in skin). Cell classification per cluster was better for the selected 31 parameters. (I) tSNE plots showing the classification of cells into clusters by using all 118 parameters and a standard single cell analytical pipeline (Seurat_v4). Donut plots indicate the distribution of the analyzed cell types (macrophages, DCs and neutrophils) in each cluster.
Extended Data Fig. 2 Behavioral landscape of the infected trachea.
(A) Heatmap of the 31 behavioral parameters used for the trachea infection analysis. For a full list of all extracted parameters please refer to Supplementary Table 2. (B) Pearson correlation matrix for all 118 parameters extracted from the trachea imaging experiment, with the selected parameters marked in red font. (C) Segmentation of cells from the trachea by four different combinations of morpho-kinetic parameters. We randomly chose 1, 6, 15, 25 and 31 parameters (list of parameter code numbers shown at right) and used them to represent the separation of neutrophils and DCs using tSNE. The original set of parameters used in Fig. 1c is at the top. Parameters are ordered from higher to lower LRI (left to right) to better visualize the classification value of each parameter used in the plot analyses. (D) LRI (score of cell identities) are proportional to the number of parameters extracted from the imaging experiments and combined to infer identities. Violin plots show the distribution of LRIs when using 1-5 parameters to classify cells in the trachea experiment, assuming that only sets of 5, 25, 31, 50 or 118 parameters are available for analysis. Note that the LRIs shown here are for the full 118 parameter set, and are not comparable with the 31 subset of selected parameters, which feature higher LRI values, as shown in the violin plots on the far right. (E) Individual analyses of the behavior of DCs and neutrophils from the original dataset, shown as tSNE plots for each population. Each behavioral parameter can be visualized and compared across cell subsets and parameters to infer positive or negative correlations, as shown for Distance to DC which negatively correlates with cell speeds in the Pearson correlation matrix of the 31 parameters used in the final analysis (F).
Extended Data Fig. 3 Behavioral landscape of the skin under ischemia-reperfusion.
(A) Representative image of I/R injury of the skin (original image on top; reconstruction of volumes and tracks at bottom). Below, examples of a typical GFPlo macrophage, a GFPhi neutrophil and an YFP+ DC used to classify the cells post-analysis. (B) Heatmap of all parameters and classified by cluster (0, 1 y 2) from the plot in Fig. 1h, and further divided into subclusters shown in (C). Below, expression plots of selected parameters. (C) tSNE plot showing all subclusters identified in the heatmap in (B). Donut plots indicate the fraction of neutrophils, DCs and macrophages in each cluster. Bottom panels show the behavioral maps generated by back-gating each cluster into the original position for each cell so that maps show the position of cells with the same behavioral profile.
Extended Data Fig. 4 Behavioral landscape of laser injury in the skin.
(A) Representative image of laser burn injury (original image left; reconstruction of volumes and tracks at right), (B) Heatmap of the all scored parameters, showing DCs and neutrophils. Expression tSNE plots of selected parameters are shown at bottom. (C) Individual analyses of the behavior of DCs and neutrophils from the original dataset, shown as tSNE plots for each population. Each behavioral parameter can be visualized and compared across cell subsets and parameters to infer random or gradient distribution for each population. For example, the location of the laser injury can be extracted as a parameter (left, yellow arrowhead) that shows graded behaviors of neutrophils relative to their distance to the wound, but not for DCs. (D) Examples of behavioral maps generated by projecting the intensity of specific parameters onto the XY location of individual cells at all time points. Actual image, plot-map by cell type and behavioral maps are shown. (E) Sub-clustering identifies two behavioral clusters of neutrophils and one for DCs (top), which were projected back onto their corresponding xyz position thus giving a profile of the distribution of behavioral clusters in the skin anatomy (middle). The neutrophil clusters feature differences in various parameters, as shown in the expression plots (arrowheads in the bottom tSNE plots). (F) Representative image of regulatory T cells (Treg) and cytotoxic T cells (CTL) in a CT26 carcinoma (red outline) in the skin, and tSNE plots of the cells classified by behavioral phenotype and by cell type. Donut plots show the match between both classifications. (G) Heatmap of the differentially scored parameters discriminating CTLs and Tregs. (H–J) Behavioral landscapes and maps of CTLs in carcinoma-bearing mice (H), neutrophils inside or outside inflamed vessels (I), and bone marrow neutrophils before and after administration of the mobilizing chemokine CXCL1 (J). Donut plots and expression plots illustrate the correlation between behavioral patterns or parameters and their localization in tissues. Dashed lines in the behavioral maps in (H–I) delineate tumor-stroma or vessel-parenchyma borders, respectively. Data are from one experiment per condition to visualize the distribution of cells in a single anatomical area.
Extended Data Fig. 5 Neutrophil states inside inflamed venules.
(A) Analysis of the cremaster dataset using Imaris software and UMAP representation show less defined behavioral clusters than using ACME (compare with Fig. 2d–f). Donut plots show the distribution of clusters in control and platelet-depleted mice. (B) Anomalous morphometric reconstructions of fast rolling cells, shown in top and side 3D views of cells moving at different speeds. Firmly adherent B1 neutrophils are shown for reference. (C) Rapid changes in morphology for neutrophils in the B2 group, following inchworm-type crawling during a 90 s recording; scale bar, 10 μm. (D) Membrane extensions (yellow arrowheads) forming around large oblate neutrophils in the B3 group, but not from B1 or B2; scale bar, 5 μm. (E) Representative micrograph of an inflamed vessel from Ly6GCre; Rosa26tdTom mouse with several neutrophils exhibiting B2 and B3 behavioral profiles (arrowheads), and “footprints” beneath B3 cells; scale bar, 10 μm. The presence of the footprints for each behavior is quantified in (F), where n is number of neutrophils analyzed. (G) Micrographs and quantification of CD11b expression measured by in vivo imaging across the different behavioral groups, with rolling neutrophils included as reference cells; data is from the indicated number of cells (in brackets), from 6 mice per group. (H) Micrographs and quantification of the number of beads phagocytosed by neutrophils from each behavioral group, including rolling cells; n is the number of cells (in brackets) from 6 mice analyzed per group. Scale bar, 5 μm. (I) Representative 3D image of an inflamed cremaster vessel showing examples of B2 and B3 neutrophils (left image), which were examined for extravasation across the endothelial wall over time (arrowheads in insets, right). (J) Percent of B3-type neutrophils that localize in junctional vs. non-junctional areas, and (K) the frequency of transendothelial migration (TEM) for each behavioral group of neutrophils; n is 5 mice per group, with the indicated number of analyzed cells (brackets). All bar graphs show mean ± SEM and data were analyzed by one-way ANOVA with Tukey’s multigroup comparison test (H, K) or unpaired two-tailed t-test (J). Number of analyzed cells per group from 3-5 mice each are indicated in brackets.
Extended Data Fig. 6 Transitional states of neutrophils in vessels.
(A) Heatmap of all parameters across all behaviors, including the three sub-groups in B2. (B) UMAP based on hierarchical clustering to identify two additional behavioral clusters within B2. (C) Distribution of cells in each sub-cluster B2.1, B2.2 and B2.3 for the indicated parameters, showing for example that cells B2.3 feature sizes and distances to the vessel wall similar to those of B3. Data analyzed by one-way ANOVA. (D) Transitions between behavioral clusters shown graphically in the UMAP (left) and quantified at right. (E) Scheme illustrating the most common transitions typically involving passage through B2, suggesting that this is an obligate transitional stage for neutrophils in inflamed vessels. Drawings in each group represent the silhouettes of representative cells at different times as in Fig. 2h.
Extended Data Fig. 7 Track parameters in the behavioral screening.
(A) Heatmap of the differentially scored behaviors among the three main behavioral groups (B1, B2 and B3). Outlined in red are the specific behaviors chosen for our screening in Fig. 3. Note that “tortuosity” is an inverse measure of “directionality”, which was used in our screening. (B) Speed and directionality obtained by epifluorescence (2D) analysis of cremasteric venules in mice with mixed chimeric bone marrow of wild-typeDsRed and non-fluorescent mutant donors, which provided internal controls for each group. Thick lines show means; The number of analyzed cells per group is shown in brackets as (control, mutant), and were obtained from at least 3 mice per group. Data analyzed by unpaired two-tailed t-test.
Extended Data Fig. 8 Morphometric parameters in the behavioral screening.
Ellipticity prolate and H/L ratios measured for individual cells in static 3D reconstructions from 24 mutant and 3 control groups as summarized in Fig. 3. Values are from cremasteric venules in mice with mixed chimeric bone marrow of wild-typeDsRed and non-fluorescent mutant donors, which provided internal controls for each group. Thick lines show means; The number of analyzed cells per group is shown in brackets as (control, mutant), and were obtained from at least 3 mice per group. Data analyzed by unpaired two-tailed t-test analysis.
Extended Data Fig. 9 Protection from myocardial injury by targeting Fgr.
(A) Micrographs of NETs (positive for citH3 and MPO; red) and vessels (blue) in cremasteric venules of wild-type subjected or not to I/R, and Fgr−/− mice subjected to I/R. Right, quantification of NETs per tissue volume; Data shown as mean ± SEM and n are number of mice per group. (B) Competitive recruitment of wild-type and Fgr−/− neutrophils to the peritoneal cavity after zymosan injection, or to the bronchoalveolar space of lungs after LPS instillation in mixed chimeric mice; n are numbers of mice analyzed. Selplg−/− neutrophils are shown for comparison of impaired migration. Values are normalized to reference wild-type competitors across the different groups and given as migration efficiencies. Data shown as mean ± SEM and n are number of mice analyzed. (C) Micrographs of Weibel-Palade bodies (WPB) and vacuoles in myocardial vessels after sham or ischemic challenge, which are quantified in (D). These measures of vascular damage are dependent on neutrophils, as shown after experimental depletion with 1A8 antibody (E). Data shown as mean ± SEM, and n is the number of micrographs analyzed, from 2 mice. (F) Effect of the Fgr antagonist TL02-59 in myocardial death upon ischemia-reperfusion, when given after ischemia at the time of reperfusion. Micrographs of heart sections at left illustrate the protective effect on myocardial death (outlined whitish regions). Data are normalized to the area at risk (AAR) and shown as mean from 4 mice per group. (G) Combination of neutrophil depletion with 1A8 antibody, and Fgr deficiency in transplanted mice. The infarcted areas are normalized with the areas at risk; n is number of mice per group. (H) Combination of neutrophil depletion with the Fgr agonist TL02-59. Data shown as mean ± SEM; n are mice per group. (I) Combination of the Fgr inhibitor in hematopoietic Fgr−/− mice, with no effect in further protecting from myocardial death; n are mice per group. (J) Myocardial fibrosis (left ventricle) determined by hematoxylin and eosin staining in control wild-type and Fgr−/− mice subjected to permanent ischemia and analyzed after 28 days. The fibrosis area is represented at right; n are mice per group. All data from (A, B) was analyzed by one-way ANOVA with Tukey’s multiple comparison test; (G–I) was analyzed by two-way ANOVA with Tukey’s multiple comparisons test. All other panels are compared by two-tailed unpaired-t test (C-F and J). The number of replicates (n) per group is indicated in each panel.
Extended Data Fig. 10 Protection from nephrotoxic injury by targeting Fgr.
(A) Schematic of the nephrotoxic injury model (top) and setup of conditions combining endotoxin (LPS) with increasing amounts of nephrotoxic serum (NTS), resulting in gradual increase in markers of kidney damage in serum and urine; n are number of mice per dose. (B) Transmission electron micrograph of kidney venules showing an example of intravascular occlusion in the NTS-treated mice, from 2 mice and 25–30 images analyzed. (C) Levels of the indicated metabolites in plasma of control and Fgr−/− mice before and after induction of glomerulonephritis with LPS plus NTS. The control group was treated with LPS only; n are mice analyzed per group. (D–E) Mice reconstituted with marrow from wild-type or Fgr−/− donor mice were infected with C. albicans (D) or S. aureus (E) and infection progression was measured by weight loss, and in the case of C. albicans infection by scoring the fungal load in kidneys (CFU); n are mice analyzed per group and data in (A) was analyzed by unpaired t-test. For (D–E) groups were compared by two-way ANOVA analyses for weight loss, and unpaired t-test for CFUs. Data in (C) was analyzed by one-way ANOVA with Tukey’s multiple comparisons test. (F) Scheme modeling neutrophil states, transitions, and delivery of inflammatory signals to the host tissues from B3 cells. Each transition is proposed to be caused by different signals, e.g. delivered by platelets and PSGL-1 for initial transition from B1 to B2, and via Fgr for transitions from B2 to B3. The number of replicates (n) per group is indicated in each panel.
Supplementary information
Supplementary Information
This file contains Supplementary Tables 1–5 and captions for Supplementary Videos 1–6.
Supplementary Video 1
Multiphoton imaging of leukocytes in influenza-infected trachea.
Supplementary Video 2
Multiphoton imaging of leukocytes during ischaemia–reperfusion injury.
Supplementary Video 3
Multiphoton imaging of DCs and neutrophils during laser-burn injury.
Supplementary Video 4
Spinning-disk imaging of inflamed venules for behavioural analysis of intravascular neutrophils.
Supplementary Video 5
Three behavioural states for intravascular neutrophils within inflamed venules.
Supplementary Video 6
Filtering out subcellular fragments.
Source data
Rights and permissions
About this article
Cite this article
Crainiciuc, G., Palomino-Segura, M., Molina-Moreno, M. et al. Behavioural immune landscapes of inflammation. Nature 601, 415–421 (2022). https://doi.org/10.1038/s41586-021-04263-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-021-04263-y
This article is cited by
-
Cecelia: a multifunctional image analysis toolbox for decoding spatial cellular interactions and behaviour
Nature Communications (2025)
-
BEHAV3D: a 3D live imaging platform for comprehensive analysis of engineered T cell behavior and tumor response
Nature Protocols (2024)
-
Immune-associated biomarkers identification for diagnosing carotid plaque progression with uremia through systematical bioinformatics and machine learning analysis
European Journal of Medical Research (2023)
-
BEHAV3D: an imaging and transcriptomics platform that unravels T cell antitumor activity
Nature Biotechnology (2023)
-
Bridging live-cell imaging and next-generation cancer treatment
Nature Reviews Cancer (2023)
