# Phenotypic clustering: a novel method for microglial morphology analysis

Franck Verdonk, Pascal Roux, Patricia Flamant, Laurence Fiette, Fernando A. Bozza, Sébastien Simard, Marc Lemaire, Benoit Plaud, Spencer L. Shorte, Tarek Sharshar, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Franck Verdonk, Pascal Roux, Patricia Flamant, Laurence Fiette, Fernando A. Bozza, et al., Phenotypic clustering: a novel method for microglial morphology analysis. Journal of Neuroinflammation, 2016, 13 (1), pp.153-. 10.1186/s12974-016-0614-7. pasteur-01389354

HAL Id: pasteur-01389354

https://pasteur.hal.science/pasteur-01389354

Submitted on 28 Oct 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
<figure><img src="image_2.png" /></figure>

Visual representation of Microglial cells in Frontal Cortex region, imaged on a confocal system (LSM700-Zeiss) with Plan-Apochromat 63X/NA 1.4 oil objective

Verdonk *et al.*
RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# Phenotypic clustering: a novel method for microglial morphology analysis

Franck Verdonk<sup>1,4,10,11</sup>, Pascal Roux<sup>2</sup>, Patricia Flamant<sup>1</sup>, Laurence Fiette<sup>1</sup>, Fernando A. Bozza<sup>3</sup>, Sébastien Simard<sup>2</sup>, Marc Lemaire<sup>4</sup>, Benoit Plaud<sup>5,9</sup>, Spencer L. Shorte<sup>2</sup>, Tarek Sharshar<sup>1,6,8,11</sup>, Fabrice Chrétien<sup>1,7,10,11*†</sup> and Anne Danckaert<sup>2*†</sup>

**Abstract**

**Background:** Microglial cells are tissue-resident macrophages of the central nervous system. They are extremely dynamic, sensitive to their microenvironment and present a characteristic complex and heterogeneous morphology and distribution within the brain tissue. Many experimental clues highlight a strong link between their morphology and their function in response to aggression. However, due to their complex "dendritic-like" aspect that constitutes the major pool of murine microglial cells and their dense network, precise and powerful morphological studies are not easy to realize and complicate correlation with molecular or clinical parameters.

**Methods:** Using the knock-in mouse model CX3CR1<sup>GFP/+</sup>, we developed a 3D automated confocal tissue imaging system coupled with morphological modelling of many thousands of microglial cells revealing precise and quantitative assessment of major cell features: cell density, cell body area, cytoplasm area and number of primary, secondary and tertiary processes. We determined two morphological criteria that are the complexity index (CI) and the covered environment area (CEA) allowing an innovative approach lying in (i) an accurate and objective study of morphological changes in healthy or pathological condition, (ii) an in situ mapping of the microglial distribution in different neuroanatomical regions and (iii) a study of the clustering of numerous cells, allowing us to discriminate different sub-populations.

**Results:** Our results on more than 20,000 cells by condition confirm at baseline a regional heterogeneity of the microglial distribution and phenotype that persists after induction of neuroinflammation by systemic injection of lipopolysaccharide (LPS). Using clustering analysis, we highlight that, at resting state, microglial cells are distributed in four microglial sub-populations defined by their CI and CEA with a regional pattern and a specific behaviour after challenge.

**Conclusions:** Our results counteract the classical view of a homogenous regional resting state of the microglial cells within the brain. Microglial cells are distributed in different defined sub-populations that present specific behaviour after pathological challenge, allowing postulating for a cellular and functional specialization. Moreover, this new experimental approach will provide a support not only to neuropathological diagnosis but also to study microglial function in various disease models while reducing the number of animals needed to approach the international ethical statements.

**Keywords:** Microglial cell morphology, Neuroinflammation, Automated high-content analysis, Clustering, Sub-population behaviour, Complexity index

<footnote>* Correspondence: fabrice.chretien@pasteur.fr; anne.danckaert@pasteur.fr</footnote>

<footnote>†Equal contributors</footnote>

<footnote n="1">¹Human Histopathology and Animal Models Unit, Infection and Epidemiology Department, Institut Pasteur, Paris, France</footnote>

<footnote n="2">²Imagopole - CITech, Institut Pasteur, Paris, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
## Background

Microglial cells are parenchymal tissue macrophages that account for 10 % of brain cells [1] involved in brain immune surveillance and homeostasis. In healthy conditions, they are involved in synaptic development and maintenance [2], neuronal survival [3] and phagocytosis [4] to maintain brain homeostasis. Microglia migrate to interact with other cell types (i.e. astrocyte and neurons) and produce a variety of factors required to induce neural progenitor differentiation [5] and neuronal apoptosis [6]. Their activation is the main component of the neuroinflammatory process, which can result either from direct brain insult or systemic inflammation leading to either neuroprotective [7, 8] or neurotoxic [9, 10] responses. Their activation is characterized by both immunological and morphological changes, including mainly a decrease in ramification up to the amoeboid form (i.e. large non-ramified cells). Hypotheses of an adaptation of their morphology to their specific function within the central nervous system (CNS) are formulated [11, 12] but neither could be actually substantiated due to two boundaries: (i) the difficulty-to-detect subtle morphological variations of microglia and also probably intermediate forms between the categories of activation that are historically described [13, 14] and (ii) the need to couple a very precise morphological single cell approach with a proteomic or transcriptomic study using specific techniques [15] taking into account the regionalization that is specific to the CNS. In fact, since the identification of microglial cells in 1932 by del Rio-Hortega [16], the assessment of microglial morphology allows characterizing their role or their level of activation [17, 18] but it remains challenging. The spectacular complexity of these cells, their morphological variability and their dense distribution within the tissue underlines the importance of a precise morphological characterization in a large number of cells distributed in the different functional CNS regions. The standard nowadays is the use of immunohistochemistry staining using CD45 or ionized calcium-binding adapter molecule 1 (Iba-1) localizing in the monocytic lineage and well expressed in microglia [19]. Those techniques, however, have some limitations. Indeed, since the expression of these markers depends on the intensity of microglial activation, immunostaining can be insufficient for accurately describing a “dendritic-like” or ramified phenotype [20, 21] that constitutes more than 90 % of the microglial cells in young mice [11]. Additionally, it required particular histological techniques such as paraffin embedding that may affect the precision of a morphological study. These limitations might be overcome by the use of transgenic mice, such as an Iba-1<sup>GFP/+</sup> [22] or CX3CR1<sup>GFP/+</sup> [23] mouse in which brain microglia express spontaneously green fluorescent protein (GFP) respectively under the control of the Iba-1 or the CX3CL1

(fractalkine) receptor locus. These transgenic models allow very precise visualization of the microglial ramifications requiring no immunostaining techniques. These models overcome the technical limits but ask the question of the quantification methods. Usually, quantification methods are based on semi-quantitative scoring and manual counting, making them time consuming, susceptible to inter- or intra-observer variability and imprecise. Furthermore, manual analyses are unable to assess either large numbers of cells or their network organization. Recently, innovative technical and/or mathematical methods have been developed allowing automated acquisition, fractal analysis [13] or segmentation of individual microglial cell shapes [24]. They have enabled to reliably assess the changes in microglial morphology albeit in limited numbers of cells (<70) [25–27]. Because of this limitation, the range of statistical analyses is also restricted while methods are available for analysing large numbers of cells and thus detecting subtle morphological phenotypes and changes therein. These methods are alternatives to conventional statistics and are able to exploit or highlight the major heterogeneity of cell populations in the same tissue, with an underlying organization that cannot be directly observed. Indeed, plasticity, reflected by slight morphological changes, is considered as a major functional property of microglial cells [28]. Using diverse criteria, it becomes now possible to discern precisely these subpopulations and structures by applying a clustering approach [29].

In efforts to develop a method allowing the assessment of the morphology of a large microglial population, and therefore for a better understanding of the microglial behaviour in different developmental, homeostatic and disease contexts, we herein propose an innovative strategy. Our approach is based on the automated acquisition of fluorescence and the measurement of morphological indexes in CX3CR1<sup>GFP/+</sup> transgenic mice that allows discriminating microglial sub-populations based on a clustering analysis. To validate this method, we compared clustering analysis to parametric and non-parametric statistical approaches assessing their capacity to detect inter-regional variability and post-stimulation changes in microglial morphology. In conclusion, our approach is extremely efficient, reproducible and accurate.

## Methods

### Animals

Male and female C57BL/6 JRj mice purchased from JANVIER LABS and in-house CX3CR1<sup>GFP/+</sup> mice aged from 9 to 11 weeks were used for these experiments. In the CX3CR1<sup>GFP/+</sup> model, the CX3CL1 receptor gene, the CX3CR1, was knocked-in with a GFP reporter gene [23]. This gene is constitutively turned on in microglial cells and thus allows us to image them selectively using the
GFP without any immunochemistry method. Mice were housed in cages in groups of seven, in a temperature-(22 ± 1.5 °C) and humidity-controlled environment, with a 12-h light/dark cycle. Mice were provided with food and water ad libitum according to international guidelines.

## Treatment conditions

Two experimental groups were considered. In the lipopolysaccharide (LPS) group (n = 6 by strain), the mice were injected intraperitoneally with 5 mg/kg of LPS from *Escherichia coli* serotype 055:B5 (Sigma-Aldrich) [30, 31] dissolved in 0.9 % saline. Twenty-four hours later, these mice were killed by cervical dislocation and the brain was collected. In parallel, mice belonging to the control group were not anaesthetized and were also killed by cervical dislocation before brain collection.

## Tissue preparation

After cervical dislocation, the brains were immediately removed and cut in a trans-sagittal plane in the interhemispheric fissure. Cerebral hemispheres were fixed during 24 h in 4 % buffered formalin (QPath, Labonord SAS, Templemars, France). Following fixation, tissue samples were sliced along a sagittal plane on a calibrated vibratome (VT1000 S, Leica, Germany) into 100-µm-thick free-floating slices. The most medial slices were used for analysis.

## Histological analysis and immunohistochemistry

Brain sections of the left hemisphere of C57BL/6 JRj mice were incubated with the rabbit antibody against Iba-1 (Wako Chemicals, Richmond, VA, 1:500) and revealed by the secondary antibody Dy488 (Jackson ImmunResearch Laboratories, Baltimore, PA). A classical protocol was used: rehydration, blocking with 20 % goat serum and 0.5 % Triton-X 100 for 2 h, incubation with primary antibody (Dako Diluent buffer, Glostrup, Denmark) overnight at 4 °C followed by incubation with secondary antibody 4 h at room temperature. The stained sections were mounted on slides and coverslipped

## Image acquisition and processing

Using a spinning disc confocal system (CellVoyager CV1000, Yokogawa, Japan) with a UPLSAPO 40×/NA 0.9 objective, sample areas were acquired as a square of 10 × 10 fields of view with a depth of 30 µm at 2-µm increments (16 focal depths) generating one volume in four regions of interest: striatum, frontal cortex, hippocampus and cerebellum. These regions were acquired sequentially allowing the coverage of approximately 3 mm² of tissue per region. Each field corresponds to a matrix of 920 × 920 pixels; the pixel size in X and Y dimensions is 0.19 µm according to the objective. The

488-nm laser was used to excite GFP or detect Iba-1 and thus to image the microglial cells.

Before the shape characterization analysis, focal stacks of each mosaic were reconstructed by combining images from the different focal depths. Each stack was subsequently divided into three 10-µm sub-volumes to allow a two-dimensional (2D) maximum intensity projection analysis (Fig. 1a), consistent with the average size of cell bodies. Mosaic, volume creation and maximum projection processing from confocal images were done using automated free plugins [32] of the ImageJ v1.50 software interface [33].

## Image analysis for characterization of microglial cell population

An automatic image analysis was performed consecutively on the three maximum projection mosaics described above, using a custom-designed script developed with the Acapella™ image analysis software (version 2.7, PerkinElmer Technologies, Waltham, USA). This script was subdivided into two subroutines: the first, for automated detection of processes (neurite detection module from Acapella™ [34]), generating morphological characteristics per cell, and the second, for the 2D in situ morphological cartography. The data workflow is illustrated step by step in Additional file 1.

## Microglial morphological criteria

Using the custom-designed script cited above, the following morphological criteria could be extracted for each microglial cell (Fig. 1b, c): a set of measured criteria as cell body and cytoplasm area, defined as the cell body area associated with the cytoplasmic area of the primary ramifications, expressed in µm²; branching characteristics such as the total number and length (µm) of ramifications and the number of primary, secondary and tertiary ramifications; and roundness (ratio between surface and perimeter squared of the cell body) and GFP intensity by whole cell.

A second set of calculated criteria extrapolated from the previous ones yielded the complexity index (CI) and the covered environment area (CEA). First, we defined the CI using two different criteria extracted from the Acapella™ script: the number of segments of each cell, a segment being defined as the length of process between two nodes, and the number of its primary ramifications. By dividing these two criteria, we obtain also a mean complexity by primary ramification for each microglial cell (Additional file 2). On the other hand, CEA represents the 2D total surface covered by its ramifications and defined as the area of the polygon formed by linking the extremities of its processes, expressed in µm². The areal density of microglial cells by region or by brain
<figure><img src="image_4.png" /><figcaption>**Fig. 1** The characterization of microglial cells by morphological criteria. a Confocal images, representing a sub-part of the analysed image in the frontal cortex region after maximum intensity projection. The individual microglia based on GFP fluorescence appears in white *outline*. The scale bar equals 50 µm. b Ramification detection based on GFP fluorescence with Acapella™ software. The segmented ramifications linked to an individual microglia are shown artificially in green, the unattributed ramifications in white. The scale bar equals 50 µm. c The morphological criteria to characterize a microglial cell. The cell body detection (blue) and cytoplasm area (pink) have performed as a starting point to characterize a microglial cell. The complexity index (green) and the covered environment area (CEA in orange) have been deduced from ramification detection. The scale bar equals 10 µm. d Two-dimensional cartography at a single cell resolution. Colours correspond to the range of complexity and CEA with a gradient from a low level of complexity and CEA (yellow) to a high level of complexity and CEA (red). The scale bar equals 50 µm. For illustration, the images are contrast adjusted to aid in visualizing the GFP expression</figcaption></figure>

was calculated by dividing the number of microglial cells selected by the scanned tissue area.

The CI revealed a completely distinct microglial phenotype, the *amoeboid cells*. The amoeboid or rod cells are characterized by a CI = 1 (no nodes) and are characteristic for activated cells, displaying engulfing, phagocytic properties [11]. Because of their particular role, we distinguish them from the other microglial cells.

## Elimination of outliers and redundant cells

We proceeded to the filtration of outliers by size (inferior to 10 µm² or superior to 500 µm²) and roundness (inferior to 0.7 considered as the limit of a noisy form) to eliminate artefacts due to tissue noise. Using the

property of our Acapella™ script to generate in situ 2D cartographies (Additional file 1), cell duplicates on two adjacent sub-volumes were rejected for analysis to avoid data duplication. This 2D localization in situ allowed us also to remove the cells located at the edges of the mosaic reasoning that they may be truncated.

## Manual and semi-automated method workflow

To assess the accuracy of the 2D reconstruction of the custom-designed script cited, we compared the criteria measured by the cell using our automated method with the analysis realized manually by three independent experimenters, as a technical benchmark: the cell body area, the number of segments, the number of primary ramifications and the measure of the CEA. The manual
method has been performed using Fiji environment (Additional file 3).

In parallel, to assess that the criteria analysis is not software specific, a free semi-automated method has been implemented in a Fiji macro using Skeletonize (2D/3D) and Sholl Analysis plugins [35, 36]. The cell body area, the number of segments and the number of primary ramifications have been extracted by this custom-designed macro.

These two methods were tested and compared with our automated method in two regions of interest, the hippocampus and the cerebellum, and in the two conditions, LPS and control.

## Data analysis and statistics

All the data extracted from Acapella™ after elimination of outliers and redundant cells were exploited following two levels: (i) by the brain and (ii) by specific region. We conducted this study using two different statistical methods. The first conventional method consisted in an approach considering only the median or the average for each morphological criterion by animal. The second approach, more powerful and original, presented results using *k*-means clustering by cell population.

To realize the conventional approach, Prism 6.0 (GraphPad Software Inc.®, USA) was used for statistical analysis by animal and regions of interest. Data were analysed via Mann-Whitney test or Student's *t* test after being assessed for normality of sample distribution. Inter-sample/inter-region variability was tested by ANOVA Kruskal-Wallis method. Qualitative traits (i.e. clustering phenotype distribution) were analysed with a chi-square ($\chi^2$) test. Spearman coefficient has been expressed to represent the correlation between two sets of data. Statistical significance is shown on the graphs ($*p < 0.05$; $**p < 0.01$; $***p < 0.001$; $****p < 0.0001$). Statistical tests used for each data set are indicated in the figure legends.

The required number of samples per group (*n*) has been evaluated with *pwr.t.test* R function [37], with $\alpha = 5\%$ and $1 - \beta = 90\%$.

## The clustering analysis

In a second time, Prism data were transferred into JMP® version 11.0 (Statistical Analysis System Institute Inc., USA) for a complete multivariate analysis by cell population. A principal component analysis (PCA) was performed to identify the correlation between the different analysed features. To detect and characterize the sub-population of microglial cells, a *k*-means clustering method, appropriated for a large set of data, has been applied (*k* = 4). The statistics of each cluster (mean and frequency) were used to characterize sub-populations and determine their phenotype, later named clusters 1 to 4. For each condition, the amount of microglial cells

analysed was about 810 by region, 2870 by brain or 20,000 by group.

## Data storage and annotations by in situ 2D cartography

All acquired and analysed mosaics have been imported into the OMERO (“OME Remote Objects technology”) image database [38] including visual results in 2D cartographies by phenotyped cell (Fig. 1d and Additional file 1). An open-source script (OMERO.csv) has been used to annotate automatically by textual information (i.e. sex, condition or clinical observation by sample or region appurtenance) our large set of imported images in the database.

## Results

Using data extracted from more than 20,000 cells per condition, we performed and compared two types of complementary statistical approaches at an inter-regional and an inter-group level. To validate the acquisition process, the GFP intensity of each microglial cell was measured. No difference in the microglial cell GFP intensity between the two conditions, whatever the brain region, was found (Table 1 and Additional file 4).

## Conventional statistical approach

For this conventional approach, two statistical analyses were performed: one to observe an inter-region variability (Figs. 2 and 3), the other to detect a difference inter-group by brain region (Tables 1 and 2).

## Microglial cell body area, cytoplasm area and density

Mean cell body and cytoplasm areas of microglial cells in the LPS group brains are significantly higher than in the control group brains (respectively, 149 vs. 74 µm² for the cytoplasm area in the hippocampus, *p* = 0.0012, Table 1). We found a high heterogeneity between the mean microglial cell body and cytoplasm areas in the different regions of the brain within a given group (*p* < 0.01 or *p* < 0.001). At the opposite, in the context of LPS-induced inflammation, the mean cell body area was homogeneous wherever the region (Fig. 2a). Moreover, in an

<table><caption>**Table 1** Morphological variability study for microglial cells between two groups</caption><thead><tr><th>Criteria</th><th>H</th><th>FC</th><th>S</th><th>C</th></tr></thead><tbody><tr><td>GFP intensity</td><td>0.2331</td><td>0.5245</td><td>0.5245</td><td>0.4394</td></tr><tr><td>Cell body area</td><td>0.0221</td><td>0.0221</td><td>0.0140</td><td>0.0012</td></tr><tr><td>Cytoplasm area</td><td>0.0012</td><td>0.0012</td><td>0.0012</td><td>0.0082</td></tr><tr><td>Complexity</td><td>0.1364</td><td>0.6154</td><td>0.6037</td><td>0.0012</td></tr><tr><td>CEA</td><td>0.9172</td><td>0.2925</td><td>0.6154</td><td>0.1783</td></tr><tr><td>Density</td><td>0.0734</td><td>0.9749</td><td>0.9172</td><td>0.1014</td></tr></tbody></table>

Values are expressed as Mann-Whitney exact *p* values; significant differences in italics
<figure><img src="image_5.png" /><figcaption>**Fig. 2** The inter-region variability by morphological criteria. Four regions have been explored: hippocampus (H), frontal cortex (FC), striatum (S) and cerebellum (C) in two different conditions, the control (**left column**) and the LPS (**right column**). **a** Historical parameters to characterize the microglial morphology: the cell body area and the cytoplasm area defined as the cell body area associated with the cytoplasmic area of the primary ramifications in µm². **b** Calculated criteria extrapolated from the Acapella™ script: the complexity index (CI) and the covered environment area (CEA), in µm². Data shown are means ± SD in the control and LPS groups (n = 7 and n = 6, respectively). The scale bars equal 10 µm. ANOVA Kruskal-Wallis test was used to compare the different regions. \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Characteristics of the amoeboid population and their inter-region variability. Bar charts represent the characteristics of the amoeboid cell morphology (characterized by a CI = 1, without nodes) and the distribution within the four explored regions: hippocampus (H), frontal cortex (FC), striatum (S) and cerebellum (C) in the two different conditions, the control (**left column**) and the LPS (**right column**). **a** Parameters to characterize the amoeboid cells morphology: the cell body area and the cytoplasm area defined as the cell body area associated with the cytoplasmic area of the primary ramifications in µm². **b** Parameters to characterize the amoeboid cell distribution: density calculated by dividing the number of microglial cells selected by the scanned tissue area (3.03 mm²) and frequency as the ratio between the number of amoeboid cells and the total number of microglial cells analysed. Data shown are means ± SD per condition (n = 7, n = 6 for control and LPS, respectively); we used ANOVA Kruskal-Wallis test. **p** < 0.01, ***p** < 0.001</figcaption></figure>
<table><caption>**Table 2** Morphological variability study for amoeboid between two groups</caption><thead><tr><th>Criteria</th><th>H</th><th>FC</th><th>S</th><th>C</th></tr></thead><tbody><tr><td>GFP intensity</td><td>0.7133</td><td>0.9021</td><td>0.7133</td><td>0.8135</td></tr><tr><td>Cell body area</td><td>0.0023</td><td>0.1014</td><td>0.8135</td><td>0.0082</td></tr><tr><td>Cytoplasm area</td><td>0.0221</td><td>0.0221</td><td>0.0012</td><td>0.1375</td></tr><tr><td>Density</td><td>0.1288</td><td>0.0221</td><td>0.1014</td><td>0.0734</td></tr><tr><td>Frequency</td><td>0.9021</td><td>0.2308</td><td>0.1276</td><td>0.5058</td></tr></tbody></table>

Values are expressed as Mann-Whitney exact *p* values; significant differences in italics

inter-regional comparison, the mean cytoplasm area was significantly greater in the cerebellum than in the other regions, in either the control or LPS groups (Additional file 4); it did not statistically differ between the hippocampus, frontal cortex and striatum, in both conditions. In an inter-group comparison (Table 1), we observed a significant increase of cell area (body and cytoplasm) in the LPS group; the microglial cell cytoplasm was twice greater after the LPS challenge regardless of the region. In each condition, the microglial density varied significantly among regions, with a higher mean density in the frontal cortex and a lower one in the cerebellum (Additional file 4). The microglial density did not differ between the control and LPS groups, in the total brain or in each region (Table 1).

### Complexity index (CI) and covered environment area (CEA)

Both microglial CI, defined as the mean complexity by primary ramification, and CEA, defined as the 2D total surface covered by microglial ramifications, showed no statistical difference between the two different groups. CI and CEA were significantly lower in the cerebellum than in the other regions, in both the control and the LPS group. We did not find any statistical difference between the hippocampus, frontal cortex and striatum regardless of the conditions (Fig. 2b). CI and CEA did not statistically differ between the LPS and control groups, in each region except the cerebellum (Table 1). The cerebellar microglial CI was significantly lower in the LPS group than in the control group (Additional file 4).

### Comparison with Iba-1 expression based on morphological criteria by sample

To strengthen our statements about the CX3CR1<sup>GFP/+</sup> mice model, we performed the same experiments based on Iba-1 expression in wildtype C57BL/6 mice and we were interested in the same morphological criteria (microglial cell body area, cytoplasm area, cytoplasm intensity, CI and CEA). In the hippocampus, we observed the same differences between the two groups compared to the GFP model with a significant increase of cell area (body and cytoplasm) in the LPS group whereas we observed no difference considering the CI and CEA

(Additional file 5). In the cerebellum, despite a trend comparable to what is seen in the GFP mouse model, there is no significant difference in the cell area. This isolated difference may be explained by the smaller number of cells analysed, about 10 times less, linked to the limits of the immunostaining on thick sections.

### Amoeboid cells

In each condition, the frequency of amoeboid cells, considered as a particular group of microglial cells with specific morphological characteristics (CI = 1) associated with a specific function, varied significantly among regions, with the lowest frequency in the frontal cortex than in the cerebellum (Fig. 3): control group: 5.64 ± 2.9 % and 14.73 ± 3.11 %, respectively, *p* = 0.0078; LPS group: 3.13 ± 1.27 % and 17.30 ± 5.53 %, respectively, *p* = 0.0018 (Additional file 4). The LPS challenge was associated with increased amoeboid cell body and cytoplasm areas among brain regions. There is no difference by region in terms of amoeboid frequencies between the two groups (Table 2 and Additional file 6).

In conclusion, using a conventional statistical approach, we found that (i) in comparison to other regions, cerebellar microglial cells presented a bigger cytoplasm, were less dense and complex, more frequently amoeboid and more responsive to LPS; (ii) the LPS challenge is associated with an increase, such as twice greater, in cell body and cytoplasm areas but not with a decrease in CI and CEA, indicating that there is no evidence for a "deramification" process.

### Comparison with benchmark methods

To assess the accuracy of our automated method, after random extraction of an analysed cell subset, three independent experimenters performed manual measures of the selected criteria as cell body area, CI and CEA. Statistical tests carried out by region showed the same trends between the control and LPS groups whatever the method considered (Additional file 7A). It is to be noted that the analysis time per cell is multiplied by 10 between the automated method and the manual method and the number of cells studied by condition is 400 times lower (Additional file 7C).

We also tested whether the results of the analysis were not software specific. Using a semi-automatic analysis with the Fiji software environment on a greater number of cells than in the manual comparison, we did not find any difference in the morphological criteria between the semi-automatic and automatic methods in each region. The results obtained showed a strong correlation between the two methods (Additional files 7B, C).
## Cell heterogeneity by condition

Although no difference was observed between the mean CI or CEA of the microglial cells in the two conditions, considering every cell of every brain in each condition, we found a high cell heterogeneity between the samples using Kruskal-Wallis test (p < 0.0001) illustrating the biological variability across individuals (Additional file 8), whatever the morphological criteria tested.

Such variability does not allow us to compare directly populations using classical statistical studies. Moreover, in order to distinguish sub-populations in our large datasets of cells, we pooled all the microglial cells from each group by region and also conducted a new statistical study using a clustering method.

## Statistical clustering approach

The principal component analysis (PCA) showed that the CI and CEA did not correlate, allowing proceeding to the cluster analysis based on these indexes. In the control group, based on the whole brain without region discrimination, the rates of the clusters 1, 2, 3 and 4 used to characterize sub-populations by the k-means clustering were 69, 18, 11 and 2%, respectively (Fig. 4a). The cutoff fixing the high or the low characteristic of one population was set as the average of each morphological criterion in the control group (Fig. 4b) and therefore allowed to discriminate four sub-populations (SP): SP1: low CEA and low CI (–/–); SP2: low CEA and high CI (–/+); SP3: high CEA and low CI (+/–); and SP4: high CEA and high CI (+/+). Considering the whole brain, the proportions of SP did not differ statistically between the control and LPS conditions (Fig. 4b). A contrario, the proportions of these sub-populations varied significantly among the regions in the control groups, as shown in Fig. 5. SP1 was most represented in the cerebellum (92% of the cells). SP2 represented from 28% of the cells in the striatum to 15% in the hippocampus. SP3 is present uniquely in the striatum at a 16% rate and SP4 from 3% in the striatum to 11% in the hippocampus. The regional distribution of the sub-populations varied significantly between the control and LPS conditions (p < 0.0001), with a significant increase in the SP4 in the striatum (from 3 to 46%) and at the opposite a significant decrease to disappearance of this population in the other areas. The SP3 increases largely in the frontal cortex, hippocampus and cerebellum. The SP1 remains stable except in the hippocampus with an increase up to 80% of the cell population (Fig. 5).

Cluster analysis based on CI and CEA revealed a regional pattern of microglial sub-population with particular responsiveness to LPS. Therefore, we obtained a new frequency of phenotyped cells and hereby defined sub-populations of microglial cells based on particular behaviour.

## Impact on the number of needed animals

The impact of a high-content automated approach is also important ethically. It is possible to assess the number of mice required to highlight the same differences as those observed using a classical statistical study and considering only the mean by animal. This evaluation of the number of subjects required (n) is based on the two-sample t test defined by the following formula:

$$n = 2 \times t \times \frac{\sigma^2}{(\bar{m}_A - \bar{m}_B)^2}$$

in which $\bar{m}_{A,B}$ are the mean of criteria for two sets of data and $\sigma$ is the common standard deviation of two samples with the hypothesis that $n_A = n_B$.

Considering a power of 90% and an alpha risk of 5%, the expected sample size showing exactly the same differences is from 1 to 27 times greater considering the CI depending on the brain area or from 6 to 100 times greater considering the CEA. This sample size is summarized in Table 3.

## Discussion

Microglial morphological analysis is the historical technique to describe the microglial cells and also the only way to study these cells within the complex environment of the central nervous system [39]. In the literature, four major microglial phenotypes are usually distinguished based on distinct morphological criteria [11]: ramified (presenting a small cell body and numerous branched ramifications) that constitute about 90% of the microglial pool between 1 and 3 months in the murine model and are considered to be the microglial “resting” state; primed (bigger cell body but unchanged ramification pattern compared to ramified phenotype); reactive (even bigger cell body, shorter, fewer and thicker ramifications); and amoeboid (two or less processes without any branch). Despite a probable link between morphology and function [40–42] and due to a large panel of slight microglial morphological changes [28], morphological analyses are subject to many hurdles, limiting their strength in terms of objectivity, precision and accuracy leading to subjective interpretations.

The objective of the present study was to develop a new powerful tool and method for describing microglial morphology and for assessing its variability among brain regions in order to further characterize and understand its behaviour in the context of development, homeostasis and disease. We confirmed the added value of our method by studying the changes related to the LPS challenge, a well-characterized model to study inflammation. This method is based on a clustering analysis integrating new developed and automatically acquired morphological indexes, the complexity index (CI) and the
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 4** Approach by clustering to track sub-populations of microglial cells in the whole brain. **a** The scatter plots illustrate, at a single cell resolution, the CEA and CI characteristics and their frequency by cluster. The symbols "+" and "x" correspond to the centre of each cluster by the control and the LPS condition, respectively. The pie charts show the cluster frequencies by k-means clustering method (k = 4), and no significant difference has been observed between the two conditions using the chi-square test. **b** Four sub-populations have been defined by the cutoff (dotted lines) fixing the high (+) or the low (−) characteristic of one sub-population in the whole brain (WB). The cutoff was defined as the average of each morphological criterion (CI and CEA) in the control group. The centre of each cluster was plotted in the graph. The pie charts represent the proportions of sub-populations by condition. The same repartitions by sub-population as by cluster were observed</figcaption></figure>
<figure><img src="image_9.png" /></figure>

covered environment area (CEA), the area covered by microglial ramifications. These two criteria are critical since different experimental studies [43–45] highlighted certain kinetics not only in morphological changes after microglial stimulation including rapid microglial process growth, extension and reorientation towards injury but also in response to healthy central nervous system (CNS) environment [46–48].

For assessing the discriminatory power of our approach, we compared the clustering analysis to parametric and non-parametric statistical approaches. We found that the clustering analysis allowed us to identify different microglial phenotypes that are heterogeneously distributed among brain regions and presenting different

behaviours under the LPS challenge. The parametric statistical tests failed to identify these sub-populations.

As a first step, we compared available methods that are manual, qualitative or semi-quantitative to ours. Differences arise in (i) a number of cells analysed 30 times greater than reported by previous studies [25–27], (ii) an objective quantification of the morphological parameters providing a reliability in numerical and resolving power and (iii) providing a significant savings in time with an analysis that is 10 times faster. Second, while the CX3CR1<sup>GFP/+</sup> transgenic model is a standard tool to study microglial morphology in in vivo [49, 50] or ex vivo experiments [51, 52], we performed the same experiments based on Iba-1 expression in C57BL/6 mice

<table><caption>**Table 3** Expected sample size to obtain significant differences between the two groups using a conventional approach</caption><thead><tr><th></th><th>CTRL (m&#773;<sub>LPS</sub>)</th><th>σ</th><th>n</th><th>CTRL (m&#773;<sub>LPS</sub>)</th><th>σ</th><th>n</th></tr></thead><tbody><tr><td></td><td>H</td><td></td><td></td><td>FC</td><td></td><td></td></tr><tr><td>CI</td><td>5.9 (4.9)</td><td>1.2</td><td>36</td><td>5.7 (6.1)</td><td>1.2</td><td>189</td></tr><tr><td>CEA (µm<sup>2</sup>)</td><td>953 (908)</td><td>253</td><td>643</td><td>950 (1184)</td><td>349</td><td>47</td></tr><tr><td>Cytoplasm area (µm<sup>2</sup>)</td><td>74.4 (149.4)</td><td>42.2</td><td>7</td><td>62.4 (123.0)</td><td>33.0</td><td>7</td></tr><tr><td></td><td>S</td><td></td><td></td><td>C</td><td></td><td></td></tr><tr><td>CI</td><td>7.9 (7.3)</td><td>1.6</td><td>114</td><td>4.5 (3.6)</td><td>0.5</td><td>7</td></tr><tr><td>CEA (µm<sup>2</sup>)</td><td>1378 (1499)</td><td>435</td><td>272</td><td>450 (360)</td><td>130</td><td>45</td></tr><tr><td>Cytoplasm area (µm<sup>2</sup>)</td><td>66.3 (129.2)</td><td>34.1</td><td>7</td><td>123.4 (177.0)</td><td>37.9</td><td>11</td></tr></tbody></table>

Where *m&#773;* is the mean of criteria for the control (LPS) set of data, *σ* is the common standard deviation of the two groups and *n* is the expected number of samples to obtain a significant difference
and observed the same trends than those observed with the GFP model considering the main morphological criteria. Nonetheless, this approach allows to analyse 10 times less cells than the CX3CR1<sup>GFP/+</sup> thanks to the immunohistochemistry techniques, confirming the interest of the use of a knock-in model.

The relevance of CI and CEA indexes, easily assessed with automated methods, is confirmed by their ability to identify morphological sub-types with a particular regional distribution and sensitivity to LPS. The lack of collinearity between these two derived indexes indicates that they provide specific morphological information that may reflect particular functional status. One may argue that for a given covered area, hypo- and hyperamified microglial cells have different roles.

The last methodological input concerns the analysis by *k*-means clustering, considered to be a simple but efficient algorithm to define sub-populations of identical phenotyped cells. It divides the data into *k* clusters, minimizing the squared distance between each data point to the centre of its cluster. The main interests of this algorithm are its speed and ease of interpretation, and it is particularly adapted for identifying sub-groups in a large dataset by cellular heterogeneity recognition [53]. The principal limitations are that it requires an a priori specification of the number of cluster centres and has a strong sensitivity to outliers and noise. Together with Kongsui and colleagues [54] in a recent study interested in the structural alterations of the microglial cells within the prefrontal cortex in rats following LPS injection, we do not find any statistical difference in the mean values of CI or CEA. Kongsui and colleagues raise the hypothesis (i) that microglial process alteration is a later phenomenon or (ii) that a substantially larger group of cells studied associated with improved analytical approaches may reveal differences. The results of our cluster analysis based on these indexes that identify major changes in the global and regional distribution of morphological sub-types support the second hypothesis. This discrepancy clearly illustrates the discriminating power of cluster analysis, when compared to parametric tests.

The regional heterogeneity of microglial morphology and the microglial effect of the LPS challenge have long been studied [1, 55]. In our study, we confirm this regionalization of the microglial distribution with a decrease in density, CI and CEA indexes from the frontal cortex to the cerebellum. Injection of LPS was associated with an increase in the cytoplasmic area and in the proportion of amoeboid cells. Kozlowski and Weimer found the same trend into the cortex in their study in 2012, correlated with an overexpression of Iba-1 and CD68 [24]. In a more original way, we observed a regional susceptibility to LPS thanks to the microglial morphology, as it is also observed at a protein level by various

experimental studies [56, 57]. This regionalization of pathophysiological processes at cellular and proteic levels supports the clinical and behavioural specific responses to neurological challenges [58, 59]. A method like ours may contribute to assess the nature of the microenvironmental factors involved in microglial shape and reactivity [60, 61].

Since the cluster analysis enables the assessment of a large population of cells per animal, it dramatically reduces the number of animals needed for testing a hypothesis. For instance, it would have been necessary to sacrifice 6 to 100 times more mice to observe a statistical difference in CI (or CEA) between the LPS- and non-LPS-treated groups, using a parametric test. This is a major ethical advantage and in compliance with the European and American requests.¹ Thanks to the large number of cells analysed in the same animal with our automated method, it is also possible to use an accurate statistical approach and consequently to dramatically improve the ethical considerations of experimental works.

Although we have found and confirmed the heterogeneity of microglial cells both in the resting and the inflammatory brain, a functional assay to correlate the morphology with function is still required. Many markers have been described to characterize different microglial activation stages, but the use of these markers only would not have been sufficiently precise, and more parameters are thus needed to further link morphology and microglial function. However, our results further confirm data from the literature as Kozlowski and Weimer showed it [24].

One of the strengths of this approach is that other functional or morphological parameters could be included in the future according to the focus of the experimenters to highlight different behaviour or improve the precision of the study. It is also possible that increasing the number of parameters could lead to the discovery of novel microglial cell states before, during and after pathological conditions.

Moreover, our technique benefits from the ease of accessibility and will allow all labs to use these parameters as a tool to characterize the microglial behaviour and further understand in a standardized manner its role in healthy and diseased condition.

Two main fields could benefit from this approach. First, through the use of optical sectioning microscopy, it is possible and easy with our method to work on a three-dimensional network. The appearance of clusters of complex cells or of similar activation states may confirm what is already highlighted in 2D or, on the contrary, may reveal new behaviours. Second, in subsequent studies, it would be useful to correlate the precise morphology of a cell that we are able to identify with its
function, possibly using single cell techniques in situ
after morphological analysis within tissue.

Conclusions

Our study presents an automated approach of morphological analysis coupled with a high-content statistical study that allows highlighting sub-populations of microglial cells, even in a healthy condition, counteracting the classical view of a homogenous resting state. These sub-populations present specific behaviour after neuroinflammatory challenge induced by LPS, allowing postulating a cellular specialization specified by its morphology. Moreover, our study confirms the need to work region by region considering this type of cells. This statistical approach may become the cornerstone of any study involving dendritic-shaped cells while seeking the reduction of the number of animals in accordance with the international guidelines for animal welfare.

Endnote

¹Directive 2010/63/EU of the European Parliament and of the Council.

Additional files

**Additional file 1:** High-content analysis workflow overview, from acquisition to statistics for the seizure of a large number of brain regions of mice, their classification and associated statistical analysis. Acquisition pipeline: four regions of interest were scanned using a spinning disc confocal system (CV1000-Yokagama): striatum (S), frontal cortex (FC), hippocampus (H) and cerebellum (C) representing approximately 13 mm² of the entire brain surface with a depth of 30 µm. The voxel size is equal to 0.19 × 0.19 × 2.0 µm³, respectively, for the X, Y and Z dimensions by stack. Analysis pipeline: using a macro implemented in ImageJ free software, each data stack was divided into three sub-volumes followed by a maximum projection for two-dimensional (2D) analysis. Each sub-volume (v1 to v3) was analysed by Acapella™ software to extract morphological criteria for each microglial cell. The original image data and generated in situ 2D cartographies were stored in Image Database for visualization, sharing and clinical annotations. At the end, statistics were done with the data extracted from this analysis pipeline. (TIF 1480 kb)

**Additional file 2:** The complexity index as new morphological criteria. The complexity index (CI) for each microglial cell was defined by dividing two different criteria: the number of segments of each cell and the number of its primary ramifications. (A) A typical schematic microglial cell presenting a circular cell body area (in red) and some ramified processes. Each of them is composed of one primary ramification (white arrow) and several sub-ramifications separated by nodes (blue arrows). One segment is defined as the length of process between two nodes. Each segment is visible in a single colour. (B) In the left column, individual microglia based on GFP fluorescence appears in white outline. In the right column, schematic representation of the microglial cell characterized by its CI (white text). The scale bars equal 10 µm. (TIF 676 kb)

**Additional file 3:** Data analysis workflow with manual and semi-automated methods. The manual method has been performed using Fiji environment and selection drawing tools to measure the cell body area, the number of roots, ramifications and the CEA. The semi-automated method has been implemented in a Fiji macro using successively Analyze Particles, Skeletonize (2D/3D) and Sholl Analysis plugins. The scale bars equal 10 µm. (TIF 1987 kb)

**Additional file 4:** Descriptive statistics for microglial cells by condition. (PDF 13 kb)

**Additional file 5:** Characterization of microglial cells by morphological criteria based on Iba-1 expression. Two regions have been explored: hippocampus (H) and cerebellum (C) in the control or LPS conditions. The scatter plots illustrate, by analysed cell, the body area (in blue), cytoplasm area (in pink) and intensity (in grey), CI (in green) or CEA (in orange) characteristics for each animal in both groups. Data shown are means ± SD. The scale bars equal 10 µm. The Mann-Whitney test was used to compare the control and LPS groups (respectively, *n* = 5 and *n* = 6).*p* < 0.05. (TIF 262 kb)

**Additional file 6:** Descriptive statistics for amoeboid cells by condition. (PDF 53 kb)

**Additional file 7:** The comparison of different quantitative analysis methods for microglial morphological criteria by cell. (A) The scatter plots illustrate, at a single cell resolution, the body area, CI or CEA characteristics for the automatic or manual method in the control and LPS groups. The Mann-Whitney test was used to compare the control and LPS groups (respectively, *n* = 5 and *n* = 6).*p* < 0.05. (B) The correlation plots of two morphological criteria between semi-automatic and automatic analysis with cell body area (in blue) and CI (in green) at a single cell resolution. Values indicate the Spearman correlation coefficient. The lines represent the linear regression. (C) The table indicates the number of cells (Nb cells), the cell body area, the CI, the CEA and the mean analysis time per cell for each condition (control and LPS groups) and each region (hippocampus and cerebellum). (TIF 424 kb)

**Additional file 8:** Quantitative analysis of microglial cell morphology in the whole brain considering CI and CEA. The scatter plots illustrate, at a single cell resolution, the CI or CEA characteristics for each animal in both groups. ANOVA Kruskal-Wallis was used to test the inter-sample heterogeneity. Data shown are means ± SD. **** *p* < 0.0001. (TIF 438 kb)

Abbreviations

CEA, Covered environment area; CI, Complexity index; CNS, Central nervous system; GFP, Green fluorescent protein; Iba-1, Ionized calcium binding adapter molecule 1; LPS, Lipopolysaccharide; OMERO, OME Remote Objects technology; PCA, Principal component analysis; SP, Sub-population

Acknowledgements

The authors thank Pierre Rocheteau and Daniel Fiole for advice and comments on the manuscript.

Funding

This work was supported by the Fondation des Gueules Cassées, the Comité d'Interface SFAR (Société Française d'Anesthésie-Réanimation), SRLF (Société de Réanimation de Langue Française), INSERM and the Fondation Pierre Deniker. Franck Verdonk received financial support from ANRT (Association Nationale de la Recherche et de la Technologie) through an Air Liquide-CIFRE contract 2012/1315. The Imagopole is part of the France Biolmaging infrastructure supported by the French National Research Agency (ANR-10-INSB-04-01, "Investments for the future") and is grateful for support from the Conseil de la Région Ile-de-France (programme Sesame 2007, project Imagopole, S.L. Shorte). SS's work was supported by the Wellcome Trust (095931/Z/11/Z; awardee Jason Swedlow, Dundee University and SLS local partner).

Availability of data and materials

All relevant data are within the paper and its additional files. The raw data studied are stored in the OMERO image database not yet accessible from an external account. A set of associated analysis scripts and ImageJ macro can be sent on request to the corresponding authors.

Authors' contributions

FV and FC conceived and designed the study. FC and LF coordinated the study. PF performed technical support for animal preparation. FV and PR carried out and supervised the image acquisitions. AD conceived, designed and coordinated the image and statistical analysis. SS designed the automatic annotation tools. FV, FC, AD and TS wrote the manuscript. SLS, FAB, BP and ML contributed in the critical revision of the manuscript. All authors read and approved the final manuscript.
**Competing interests**

The authors declare that they have no competing interests.

**Consent for publication**

Not applicable.

**Ethics approval and consent to participate**

All protocols were reviewed by the Institut Pasteur, the competent authority, for compliance with the French and European regulations on Animal Welfare and with Public Health Service recommendations. This project has been reviewed and approved (# 2013-0044) by the Institut Pasteur Ethics Committee (C2EA 89 - CETEA).

**Author details**

¹Human Histopathology and Animal Models Unit, Infection and Epidemiology Department, Institut Pasteur, Paris, France. ²Imagopole - CI Tech, Institut Pasteur, Paris, France. ³ICU, Instituto de Pesquisa Clínica Evandro Chagas, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil. ⁴Air Liquide Santé International, World Business Line Healthcare, Medical R&D, Paris-Saclay Research Center, 1 chemin de la Porte des Loges, Jouy-en-Josas, France. ⁵Department of Anaesthesiology and Surgical Intensive Care, Saint-Louis University Hospital of Paris, Paris, France. ⁶Department of Intensive Care, Raymond Poincare University Hospital, Garches, France. ⁷Laboratoire hospitalo-universitaire de Neuropathologie, Centre Hospitalier Sainte Anne, Paris, France. ⁸Versailles Saint Quentin University, Versailles, France. ⁹Paris Diderot University, Paris, France. ¹⁰Paris Descartes University, Sorbonne Paris Cité, Paris, France. ¹¹TRIGGERSEP, F-CRIN Network, Toulouse, France.

Received: 14 January 2016 Accepted: 6 June 2016

Published online: 17 June 2016

**References**

1. Lawson LJ, Perry VH, Dri P, Gordon S. Heterogeneity in the distribution and morphology of microglia in the normal adult mouse brain. Neuroscience. 1990;39(1):151-70.

2. Tremblay ME, Lowery RL, Majewska AK. Microglial interactions with synapses are modulated by visual experience. PLoS Biol. 2010;8(11):e1000527.

3. Ueno M, Fujita Y, Tanaka T, Nakamura Y, Kikuta J, Ishii M, et al. Layer V cortical neurons require microglial support for survival during postnatal development. Nat Neurosci. 2013;16(5):543-51.

4. Marin-Teva JL, Dusart I, Colin C, Gervais A, van Rooijen N, Mallat M. Microglia promote the death of developing Purkinje cells. Neuron. 2004;41(4):535-47.

5. Hsieh J, Aimone JB, Kaspar BK, Kuwabara T, Nakashima K, Gage FH. IGF-I instructs multipotent adult neural progenitor cells to become oligodendrocytes. J Cell Biol. 2004;164(1):111-22.

6. Paolicelli RC, Bolasco G, Pagani F, Maggi L, Scianni M, Panzanelli P, et al. Synaptic pruning by microglia is necessary for normal brain development. Science. 2011;333(6048):1456-8.

7. Block ML, Zecca L, Hong JS. Microglia-mediated neurotoxicity: uncovering the molecular mechanisms. Nat Rev Neurosci. 2007;8(1):57-69.

8. Hanisch UK, Kettenmann H. Microglia: active sensor and versatile effector cells in the normal and pathologic brain. Nat Neurosci. 2007;10(11):1387-94.

9. Doens D, Fernandez PL. Microglia receptors and their implications in the response to amyloid beta for Alzheimer's disease pathogenesis. J Neuroinflammation. 2014;11:48.

10. More SV, Kumar H, Kim JS, Song SY, Choi DK. Cellular and molecular mediators of neuroinflammation in the pathogenesis of Parkinson's disease. Mediators Inflamm. 2013;2013:952375.

11. Torres-Platas SG, Comeau S, Rachalski A, Bo GD, Cruceanu C, Turecki G, et al. Morphometric characterization of microglial phenotypes in human cerebral cortex. J Neuroinflammation. 2014;11:12.

12. Vela D, Dalmau I, Gonzalez B, Castellano B. Morphology and distribution of microglial cells in the young and adult mouse cerebellum. J Comp Neurol. 1995;361(4):602-16.

13. Karperien A, Ahammer H, Jelinek HF. Quantitating the subtleties of microglial morphology with fractal analysis. Front Cell Neurosci. 2013;7:3.

14. Perego C, Fumagalli S, De Simoni MG. Temporal pattern of expression and colocalization of microglia/macrophage phenotype markers following brain ischemic injury in mice. J Neuroinflammation. 2011;8:174.

15. Lovatt D, Ruble BK, Lee J, Dueck H, Kim TK, Fisher S, et al. Transcriptome in vivo analysis (TIVA) of spatially defined single cells in live tissue. Nat Methods. 2014;11(2):190-6.

16. del Rio-Hortega P. Microglia. Cytology and Cellular Pathology of the Nervous System, vol. 2. New York: W. Penfield, Hoeber; 1932. p. 482-534.

17. Aguzzi A, Barres BA, Bennett ML. Microglia: scapegoat, saboteur, or something else? Science. 2013;339(6116):156-61.

18. Stence N, Waite M, Dailey ME. Dynamics of microglial activation: a confocal time-lapse analysis in hippocampal slices. Glia. 2001;33(3):256-66.

19. Ohyu J, Marumo G, Ozawa H, Takashima S, Nakajima K, Kohsaka S, et al. Early axonal and glial pathology in fetal sheep brains with leukomalacia induced by repeated umbilical cord occlusion. Brain Dev. 1999;21(4):248-52.

20. Ito D, Imai Y, Ohsawa K, Nakajima K, Fukuuchi Y, Kohsaka S. Microglia-specific localisation of a novel calcium binding protein, Iba1. Brain Res Mol Brain Res. 1998;57(1):1-9.

21. Ito D, Tanaka K, Suzuki S, Dembo T, Fukuuchi Y. Enhanced expression of Iba1, ionized calcium-binding adapter molecule 1, after transient focal cerebral ischemia in rat brain. Stroke. 2001;32(5):1208-15.

22. Hirasawa T, Ohsawa K, Imai Y, Ondo Y, Akazawa C, Uchino S, et al. Visualization of microglia in living tissues using Iba1-EGFP transgenic mice. J Neurosci Res. 2005;81(3):357-62.

23. Jung S, Aliberti J, Graemmel P, Sunshine MJ, Kreutzberg GW, Sher A, et al. Analysis of fractalkine receptor CX(3)CR1 function by targeted deletion and green fluorescent protein reporter gene insertion. Mol Cell Biol. 2000;20(11):4106-14.

24. Kozlowski C, Weimer RM. An automated method to quantify microglia morphology and application to monitor activation state longitudinally in vivo. PLoS One. 2012;7(2):e31814.

25. Baron R, Babcock AA, Nemirovsky A, Finsen B, Monsonego A. Accelerated microglial pathology is associated with Abeta plaques in mouse models of Alzheimer's disease. Aging Cell. 2014;13(4):584-95.

26. Rangroo Thrane V, Thrane AS, Chang J, Alleluia V, Nagelhus EA, Nedergaard M. Real-time analysis of microglial activation and motility in hepatic and hyperammonemic encephalopathy. Neuroscience. 2012;220:247-55.

27. Valous NA, Lahrmann B, Zhou W, Veltkamp R, Grabe N. Multistage histopathological image segmentation of Iba1-stained murine microglias in a focal ischemia model: methodological workflow and expert validation. J Neurosci Methods. 2013;213(2):250-62.

28. Perry VH, Cunningham C, Holmes C. Systemic infections and inflammation affect chronic neurodegeneration. Nat Rev Immunol. 2007;7(2):161-7.

29. Ingalhalikar M, Smith AR, Bloy L, Gur R, Roberts TP, Verma R. Identifying sub-populations via unsupervised cluster analysis on multi-edge similarity graphs. Med Image Comput Comput Assist Interv. 2012;15(Pt 2):254-61.

30. Fan K, Li D, Zhang Y, Han C, Liang J, Hou C, et al. The induction of neuronal death by up-regulated microglial cathepsin H in LPS-induced neuroinflammation. J Neuroinflammation. 2015;12:54.

31. Qin L, Wu X, Block ML, Liu Y, Breese GR, Hong JS, et al. Systemic LPS causes chronic neuroinflammation and progressive neurodegeneration. Glia. 2007;55(5):453-62.

32. Preibisch S, Saalfeld S, Tomancak P. Globally optimal stitching of tiled 3D microscopic image acquisitions. Bioinformatics. 2009;25(11):1463-5.

33. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9(7):676-82.

34. Vallotton P, Lagerstrom R, Sun C, Buckley M Wang D, De Silva M, et al. Automated analysis of neurite branching in cultured cortical neurons using HCA-Vision. Cytometry A. 2007;71(10):889-95.

35. Arganda-Carreras I, Fernandez-Gonzalez R, Munoz-Barrutia A, Ortiz-De-Solorzano C. 3D reconstruction of histological sections: application to mammary gland tissue. Microsc Res Tech. 2010;73(11):1019-29.

36. Kroner A, Greenhalgh AD, Zarruk JG, Passos Dos Santos R, Gastel M, David S. TNF and increased intracellular iron alter macrophage polarization to a detrimental M1 phenotype in the injured spinal cord. Neuron. 2014;83(5):1098-116.

37. R Data Analysis Examples: University of California; [cited 2015 12 June 2005]. Available from: http://www.ats.ucla.edu/stat/r/dae/t_test_power2.htm.

38. Allan C, Burel JM, Moore J, Blackburn C, Linkert M, Loynton S, et al. OMERO: flexible, model-driven data management for experimental biology. Nat Methods. 2012;9(3):245-53.
39. Berry M, Butt AM, Wilkin G, Perry VH, Structure and function of glia in the central nervous system. In: Graham D, Lantos P, editors. Greenfield's neuropathology, 7th ed. London: Arnold; 2002. p. 75-122.

40. Graeber MB, Streit WJ. Microglia: biology and pathology. Acta Neuropathol. 2010;119(1):89-105.

41. Kettenmann H, Hanisch UK, Noda M, Verkratsky A. Physiology of microglia. Physiol Rev. 2011;91(2):461-553.

42. Soltys Z, Ziaja M, Pawlinski R, Setkowicz Z, Janeczko K. Morphology of reactive microglia in the injured cerebral cortex. Fractal analysis and complementary quantitative methods. J Neurosci Res. 2001;63(1):90-7.

43. Beynon SB, Walker FR. Microglial activation in the injured and healthy brain: what are we really talking about? Practical and theoretical issues associated with the measurement of changes in microglial morphology. Neuroscience. 2012;225:162-71.

44. Davalos D, Grutzendler J, Yang G, Kim JV, Zuo Y, Jung S, et al. ATP mediates rapid microglial response to local brain injury in vivo. Nat Neurosci. 2005; 8(6):752-8.

45. Nimmerjahn A, Kirchhoff F, Helmchen F. Resting microglial cells are highly dynamic surveillants of brain parenchyma in vivo. Science. 2005;308(5726):1314-8.

46. Dissing-Olesen L, LeDue JM, Rungta RL, Hefendehl JK, Choi HB, MacVicar BA. Activation of neuronal NMDA receptors triggers transient ATP-mediated microglial process outgrowth. J Neurosci. 2014;34(32):10511-27.

47. Hinwood M, Tynan RJ, Charnley JL, Beynon SB, Day TA, Walker FR. Chronic stress induced remodeling of the prefrontal cortex: structural re-organization of microglia and the inhibitory effect of minocycline. Cereb Cortex. 2013;23(8):1784-97.

48. Walker FR, Nilsson M, Jones K. Acute and chronic stress-induced disturbances of microglial plasticity, phenotype and function. Curr Drug Targets. 2013;14(11):1262-76.

49. Fourgeaud L, Traves PG, Tufail Y, Leal-Bailey H, Lew ED, Burrola PG, et al. TAM receptors regulate multiple features of microglial physiology. Nature. 2016;532(7598):240-4.

50. Avignone E, Lepleux M, Angibaud J, Nagerl UV. Altered morphological dynamics of activated microglia after induction of status epilepticus. J Neuroinflammation. 2015;12:202.

51. Umekawa T, Osman AM, Han W, Ikeda T, Blomgren K. Resident microglia, rather than blood-derived macrophages, contribute to the earlier and more pronounced inflammatory reaction in the immature compared with the adult hippocampus after hypoxia-ischemia. Glia. 2015;63(12):2220-30.

52. Kim I, Mlsna LM, Yoon S, Le B, Yu S, Xu D, et al. A postnatal peak in microglial development in the mouse hippocampus is correlated with heightened sensitivity to seizure triggers. Brain Behav. 2015;5(12):e00403.

53. Altschuler SJ, Wu LF. Cellular heterogeneity: do differences make a difference? Cell. 2010;141(4):559-63.

54. Kongsui R, Johnson SJ, Graham BA, Nilsson M, Walker FR. A combined cumulative threshold spectra and digital reconstruction analysis reveal structural alterations of microglia within the prefrontal cortex following low-dose LPS administration. Neuroscience. 2015;310:629-40.

55. Kongsui R, Beynon SB, Johnson SJ, Walker FR. Quantitative assessment of microglial morphology and density reveals remarkable consistency in the distribution and morphology of cells within the healthy prefrontal cortex of the rat. J Neuroinflammation. 2014;11:182.

56. Pintado C, Revilla E, Vizuete ML, Jimenez S, Garcia-Cuervo L, Vitorica J, et al. Regional difference in inflammatory response to LPS-injection in the brain: role of microglia cell density. J Neuroimmunol. 2011;238(1-2):44-51.

57. Kim WG, Mohney RP, Wilson B, Jeohn GH, Liu B, Hong JS. Regional difference in susceptibility to lipopolysaccharide-induced neurotoxicity in the rat brain: role of microglia. J Neurosci. 2000;20(16):6309-16.

58. Bonow RH, Aid S, Zhang Y, Becker KG, Bosetti F. The brain expression of genes involved in inflammatory response, the ribosome, and learning and memory is altered by centrally injected lipopolysaccharide in mice. Pharmacogenomics J. 2009;9(2):116-26.

59. Sun J, Zhang S, Zhang X, Zhang X, Dong H, Qian Y. IL-17A is implicated in lipopolysaccharide-induced neuroinflammation and cognitive impairment in aged rats via microglial activation. J Neuroinflammation. 2015;12:165.

60. De Lucia C, Rinchon A, Olmos-Alonso A, Riecken K, Fehse B, Boche D, et al. Microglia regulate hippocampal neurogenesis during chronic neurodegeneration. Brain Behav Immun. 2015.

61. Scheffel J, Regen T, Van Rossum D, Seifert S, Ribes S, Nau R, et al. Toll-like receptor activation reveals developmental reorganization and unmasks responder subsets of microglia. Glia. 2012;60(12):1930-43.

Submit your next manuscript to BioMed Central and we will help you at every step:

• We accept pre-submission inquiries

• Our selector tool helps you to find the most relevant journal

• We provide round the clock customer support

• Convenient online submission

• Thorough peer review

• Inclusion in PubMed and all major indexing services

• Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central