Comparative fractal analysis of cultured glia derived from optic nerve and brain demonstrate different rates of morphological differentiation

Microglial morphometric analysis: so many options, so little consistency

Quantification of microglial activation through morphometric analysis has long been a staple of the neuroimmunologist's toolkit. Microglial morphological phenomics can be conducted through either manual classification or constructing a digital skeleton and extracting morphometric data from it. Multiple open-access and paid software packages are available to generate these skeletons via semi-automated and/or fully automated methods with varying degrees of accuracy. Despite advancements in methods to generate morphometrics (quantitative measures of cellular morphology), there has been limited development of tools to analyze the datasets they generate, in particular those containing parameters from tens of thousands of cells analyzed by fully automated pipelines. In this review, we compare and critique the approaches using cluster analysis and machine learning driven predictive algorithms that have been developed to tackle these large datasets, and propose improvements for these methods. In particular, we highlight the need for a commitment to open science from groups developing these classifiers. Furthermore, we call attention to a need for communication between those with a strong software engineering/computer science background and neuroimmunologists to produce effective analytical tools with simplified operability if we are to see their wide-spread adoption by the glia biology community.

Effects of Kinesiology Tape on Non-linear Center of Mass Dispersion During the Y Balance Test

Static taping of the ankle or knee joint is a common method of reducing risk of injury by providing mechanical stability. An alternative taping technique employs kinesiology tape, which has the additional benefit of improving functionality by stimulating proprioception. There is substantial disagreement whether kinesiology tape shows significant differences in proprioception and postural stability as compared to rigid/static tape when applied at the lower limb. The current study investigated the effects of kinesiology tape and static tape during a Y Balance Test on center of mass as an indicator for postural stability. Forty-one individuals, free of injury, performed the Y Balance Test under the three conditions; no tape, kinesiology tape, and static tape applied at the lower limb to the quadriceps, triceps surae and ankle joint. All participants completed the Y Balance Test to determine whether any significant differences could be observed using center of mass movement as a surrogate measure for balance and proprioception. The Minkowski-Bouligand and box-counting fractal dimension analyses were used as measures of the dynamic changes in the center of mass whilst undertaking the Y Balance Test. Statistical analyses included the Kruskal Wallis test to allow for non-normally distributed data and a Bonferroni corrected pairwise T-test as a post hoc test to ascertain pairwise differences between the three taping conditions. Significance was set at 0.05. The fractal analyses of the dynamic changes in center of mass showed significant differences between the control and both the static tape and kinesiology tape groups (p = 0.021 and 0.009, respectively). The current study developed a novel measure of dynamic changes in the center of mass during a set movement that indicated real-time processing effects during a balance task associated with the type of taping used to enhance postural stability.

FGF2 deficit during development leads to specific neuronal cell loss in the enteric nervous system

The largest part of the peripheral nervous system is the enteric nervous system (ENS). It consists of an intricate network of several enteric neuronal subclasses with distinct phenotypes and functions within the gut wall. The generation of these enteric phenotypes is dependent upon appropriate neurotrophic support during development. Glial cell line-derived neurotrophic factor (GDNF) and fibroblast growth factor-2 (FGF2) play an important role in the differentiation and function of the ENS. A lack of GDNF or its receptor (Ret) causes intestinal aganglionosis in mice, while fibroblast growth factor receptor signaling antagonist is identified as regulating proteins in the GDNF/Ret signaling in the developing ENS. Primary myenteric plexus cultures and wholemount preparations of wild type (WT) and FGF2-knockout mice were used to analyze distinct enteric subpopulations. Fractal dimension (D) as a measure of self-similarity is an excellent tool to analyze complex geometric shape and was applied to classify the subclasses of enteric neurons concerning their individual morphology. As a consequence of a detailed analysis of subpopulation variations, wholemount preparations were stained for the calcium binding proteins calbindin and calretinin. The fractal analysis showed a reliable consistence of subgroups with different fractal dimensions (D) in each culture investigated. Seven different neuronal subtypes could be differentiated according to a rising D. Within the same D, the neurite length revealed significant differences between wild type and FGF2-knockout cultures, while the subclass distribution was also altered. Depending on the morphological characteristics, the reduced subgroup was supposed to be a secretomotor neuronal type, which could be confirmed by calbindin and calretinin staining of the wholemount preparations. These revealed a reduction up to 40 % of calbindin-positive neurons in the FGF2-knockout mouse. We therefore consider FGF2 playing a more important role in the fine-tuning of the ENS during development as previously assumed.

Effects of fractal surface on C6 glioma cell morphogenesis and differentiation in vitro

Neurons and glial cells in the brain are surrounded by a fractal environment. A fractal alkylketene dimmer (AKD) surface was shown to provide such a biomimetic environment for glial cell culture. However, little is known about the effects of fractal surface on the complexity of cell morphology. In particular, whether fractal surface induces glial cell differentiation remains to be elucidated. The present work, thus determined the fractal dimension (FD) of cell complexity with a geometrically calculational parameter, the expressions of GFAP gene and protein in C6 glioma cells on fractal AKD, non-fractal AKD and PLL-coated surfaces. Fractal surface suppressed the proliferation of glioma cell, and significantly increased the length and number of cell process. Furthermore, the enhanced values of FD were accompanied with the expressions of GFAP gene and protein, especially that of gene. However, cells on non-fractal and PLL surface proliferated gradually along with the culture time, showing the fibroblast-like morphology, and accompanied with the consistent expressions of GFAP gene and protein. These results suggested that C6 glioma cell differentiation can be induced by fractal AKD surface.

A three-dimensional fractal analysis method for quantifying white matter structure in human brain

Fractal dimension (FD) is increasingly used to quantify complexity of brain structures. Previous research that analyzed FD of human brain mainly focused on two-dimensional measurements. In this study, we developed a three-dimensional (3D) box-counting method to measure FD of human brain white matter (WM) interior structure, WM surface and WM general structure simultaneously. This method, which firstly incorporates a shape descriptor (3D skeleton) representing interior structure and combines the three features, provides a more comprehensive characterization of WM structure. WM FD of different brain segments was computed to test robustness of the method. FDs of fractal phantoms were computed to test the accuracy of the method. The consistency of the computed and theoretical FD values suggests that our method is accurate in measuring FDs of fractals. Statistical analysis was performed to examine sensitivity of the method in detecting WM structure differences in a number of young and old subjects. FD values of the WM skeleton and surface were significantly greater in young than old individuals, indicating more complex WM structures in young people. These results suggest that our method is accurate in quantifying three-dimensional brain WM structures and sensitive in detecting age-related degeneration of the structures.

Fractal properties of the human genome

The fractal dimension of the human chromosomes and four other genomes were determined using the box counting method. Human chromosomes exhibited a fractal dimension (D) of about 0.8, while values for a bacteria, yeast, worm and plant were higher. Analysis of three human chromosomes over five orders of magnitude of scale (from 10(8) to 10(4) bp), showed D to be non-constant at the smaller scales, when introns were included as gaps. The relationship between D and gene density fit an empirical equation related to that expected from theory, and allowed for the calculation of the fractal initiator or self-similarity ratio. This value (0.57) was constant at all scales for human chromosomes, and was similar for other species, except for Arabidopsis.

Fractal dimension in human cerebellum measured by magnetic resonance imaging

Fractal dimension has been used to quantify the structures of a wide range of objects in biology and medicine. We measured fractal dimension of human cerebellum (CB) in magnetic resonance images of 24 healthy young subjects (12 men and 12 women). CB images were resampled to a series of image sets with different 3D resolutions. At each resolution, the skeleton of the CB white matter was obtained and the number of pixels belonging to the skeleton was determined. Fractal dimension of the CB skeleton was calculated using the box-counting method. The results indicated that the CB skeleton is a highly fractal structure, with a fractal dimension of 2.57 +/- 0.01. No significant difference in the CB fractal dimension was observed between men and women.

Characterization of Melanophore Morphology by Fractal Dimension Analysis

Fractal or focal dimension (FD) analysis is a valuable tool to identify physiologic stimuli at the cellular and tissue levels that allows for quantification of cell perimeter complexity. The FD analysis was determined on fluorescence images of caffeine- or epinephrine-treated (or untreated control) killifish Fundulus heteroclitus (Linneaus) melanophores in culture. Cell perimeters were indicated by rhodamine-phalloidin labeling of cortical microfilaments using box-counting FD analysis. Caffeine-treated melanophores displayed dispersed melanosomes in cells with less serrated edges and reduced FD and complexity. Complexity in epinephrine-treated cells was significantly higher than the caffeine-treated cells or in the control. Cytoarchitectural variability of the cell perimeter is expected because cells change shape when cued with agents. Epinephrine-treated melanophores demonstrated aggregated melanosomes in cells with more serrated edges, significantly higher FD and thus complexity. Melanophores not treated with caffeine or epinephrine produced variable distributions of melanosomes and resulted in cells with variably serrated edges and intermediate FD with a larger SE of the regression and greater range of complexity. Dispersion of melanosomes occurs with rearrangements of the cytoskeleton to accommodate centrifugal distribution of melanosomes throughout the cell and to the periphery. The loading of melanosomes onto cortical microfilaments may provide a less complex cell contour, with the even distribution of the cytoskeleton and melanosomes. Aggregation of melanosomes occurs with rearrangements of the cytoskeleton to accommodate centripetal distribution of melanosomes. The aggregation of melanosomes may contribute to centripetal retraction of the cytoskeleton and plasma membrane. The FD analysis is, therefore, a convenient method to measure contrasting morphologic changes within stimulated cells.

Morphological development of microglia in the postnatal rat brain

Morphological transformation of lectin-positive microglia/macrophages in the developing rat cerebral hemisphere was analysed using quantitative methods. During the first postnatal month, the cells showed increases in their size and fractal dimension accompanied by a simultaneous decrease in their solidity. Regional variations in dynamics of the process indicated the existence of spatio-temporal developmental gradients within the cerebral hemisphere wall which might correspond with regional patterns of neuronal differentiation. Results of the present study prove that the quantitative methods can be the source of reliable data replacing subjective cell typologies.

Morphological transformations of cells immunopositive for GFAP, TrkA or p75 in the CA1 hippocampal area following transient global ischemia in the rat. A quantitative study

Transient global ischemia induces intensive neuronal degeneration in the hippocampal CA1 pyramidal layer, accompanied by reactive transformation of glial cells. Previously, we have shown using the double immunostaining method that the NGF receptors (NGFR) p75 and TrkA are expressed mainly on subpopulations of GFAP+ astrocytes, and this expression increases progressively after ischemia. In the presented study, we analyzed quantitatively the morphological transformations of cells immunopositive for GFAP or NGF receptors in the stratum radiatum of the CA1 hippocampal area in different survival periods after ischemia, evoked by 10-min cardiac arrest in adult rats. In control brains, NGF receptors were expressed only on small cells with poorly ramified processes. After ischemia, the NGFR+ cells increased in size and morphological complexity (measured using fractal analysis). However, even 2 weeks after ischemia these cells did not reach the size and value of the fractal dimension typical of the largest GFAP+ astrocytes. Moreover, the reaction of NGFR+ cells was significantly delayed in comparison with the total astrocyte population. The obtained results suggest that NGF receptors are expressed mainly by immature astrocytes and ischemia induces the maturation of these cells.

Development of rat CA1 neurones in acute versus organotypic slices: role of experience in synaptic morphology and activity

Despite their wide use, the physiological relevance of organotypic slices remains controversial. Such cultures are prepared at 5 days postnatal. Although some local circuitry remains intact, they develop subsequently in isolation from the animal and hence without plasticity due to experience. Development of synaptic connectivity and morphology might be expected to proceed differently under these conditions than in a behaving animal. To address these questions, patch-clamp techniques and confocal microscopy were used in the CA1 region of the rat hippocampus to compare acute slices from the third postnatal week with various stages of organotypic slices. Acute slices prepared at postnatal days (P) 14, 17 and 21 were found to be developmentally equivalent to organotypic slices cultured for 1, 2 and 3 weeks, respectively, in terms of development of synaptic transmission and dendritic morphology. The frequency of inhibitory and excitatory miniature synaptic currents increased in parallel. Development of dendritic length and primary branching as well as spine density and proportions of different spine types were also similar in both preparations,at these corresponding stages. The most notable difference between organotypic and acute slices was a four- to five-fold increase in the absolute frequency of glutamatergic (but not GABAergic)miniature postsynaptic currents in organotypic slices. This was probably related to an increase in complexity of higher order dendritic branching in organotypic slices, as measured by fractal analysis, resulting in an increased total synapse number. Both increased excitatory miniature synaptic current frequency and dendritic complexity were already established during the first week in culture. The level of complexity then stayed constant in both preparations over subsequent stages, with synaptic frequency increasing in parallel. Thus, although connectivity was greater in organotypic slices, once this was established, development continued in both preparations at are markably similar rate. We conclude that, for the parameters studied, changes seem to be preprogrammed by 5 days and their subsequent development is largely independent of environment.

PPAR gamma activators induce growth arrest and process extension in B12 oligodendrocyte-like cells and terminal differentiation of cultured oligodendrocytes

Peroxisome proliferator-activated receptors (PPARs) are key transcription factors in the control of lipid homeostasis and cell differentiation, but little is known about their function in oligodendrocytes, the major lipid-synthesizing cells in the central nervous system (CNS). Using the B12 oligodendrocyte-like cell line and rat spinal cord-derived oligodendrocytes, we evaluated the importance of PPARgamma in the maturation process of these cells. B12 cells express all PPAR isoforms (alpha, beta/delta, and gamma), as assessed by RT-PCR, Western-blot, and transactivation assays. B12 cells respond specifically to PPARgamma agonists by arresting cell proliferation and extending cell processes, events that are blocked by the PPARgamma antagonist GW9662. In addition, alkyl-dihydroxyacetone phosphate synthase (ADAPS), a key peroxisomal enzyme involved in the synthesis of myelin-rich lipid plasmalogens, is increased in PPARgamma agonist-treated B12 cells. In contrast with B12 cells, both immature and mature isolated spinal cord oligodendrocytes presented a high and similar expression level of ADAPS, as assessed by immunocytochemistry. However, as in B12 cells, isolated spinal cord oligodendrocytes were also found to respond specifically to PPARgamma agonists with a four-fold increase in the number of mature cells. Our data suggest a relevant role for PPARgamma in oligodendrocyte lipid metabolism and differentiation.

Identification of living oligodendrocyte developmental stages by fractal analysis of cell morphology

The Mandelbrot's fractal dimension (D), a measure of shape complexity, has been used to quantify the complex morphology of living cells. Previous studies on glial cells have shown that as cells increase in morphological complexity, their "D" value increases, suggesting that "D" could be used to estimate their stage of differentiation. In the present study the box-counting method was used to calculate the "D" values of rat cerebellar oligodendrocytes during their differentiation in primary culture. These values were correlated with the immunoreactivity of cells to antigenic markers commonly used for assessing their stages of differentiation: A2B5, O4 and anti-galactocerebroside (Gal-C). Our results show that changes of the fractal dimension during differentiation follow the well known pattern of markers expression by these cells. These results demonstrate that A2B5-, O4-, and Gal-C-expressing oligodendrocytes can be confidently estimated from their respective fractal dimension values. Based on this immunocytochemical calibration, the calculation of "D" allows an easy and fast determination of the developmental stage of living (unstained) oligodendrocytes before the study of their physiological characteristics. Using this method we precisely identified living oligodendrocyte progenitors and early pro-oligodendrocytes expressing voltage-activated sodium currents that is a common characteristic of these two immature developmental stages (Sontheimer et al. [1989b] Neuron 2:1135-1145).

The effect of type 1 astrocytes on neuronal complexity: a fractal analysis

Embryonic, ventral spinal cord neurons were grown on poly(d-lysine) (PDL) or on a monolayer of type 1 astrocytes. At various times from 6 h to 2 weeks postplating, cells were fluorescently labeled and fixed with 4% paraformaldehyde. The cell surface immunoreaction allowed visualization of neurons in their entirety, namely, cell bodies and various membranous extensions that included lamellipodia, growth cones, axons, and dendrites. Outlines were drawn for individual neurons and their fractal dimension (D) was calculated. Neurons on poly(d-lysine) reached a peak D at 3 days in vitro, 1 day later than neurons on astrocytes (2 days in vitro). The maximum D was greater for cells on poly(d-lysine) when compared with neurons on astrocytes. In a second experiment the maximum D was similar for neurons on both surfaces but neurons on PDL maintained a higher D for a much longer period than neurons on astrocytes. An examination of fluorescent images revealed that neurons on poly(d-lysine) exhibited lamellipodia and large growth cones for several days and these structures were likely responsible for the high D seen in these cells. These structures were rarely observed in neurons plated on astrocytes. Interestingly, D on both surfaces decreased to a similar value at between 1 and 2 weeks in vitro. The trend for D in these cultures, an initial increase to a peak value followed by a decrease to a stable value, is discussed in light of the chemical nature of the two surfaces and synapse formation and stabilization.

Analysis of fractal dimension of O2A glial cells differentiating in vitro

Fractal dimension is a quantitative measure of morphological complexity. Glial cells of the oligodendrocyte-type 2 astrocyte (O2A) lineage exhibit increasing morphological complexity as they differentiate in vitro. Enriched populations of O2A progenitor cells isolated from neonatal rat cerebral hemispheres or optic nerves were allowed to differentiate in vitro, and their fractal dimensions were measured over time. The fractal dimensions of the maturing cells correlated with perceived complexity; cells with elaborate process branching had larger fractal dimensions than cells with a simpler morphology. An analysis of changes in fractal dimension revealed distinct rates of growth for both oligodendrocytes and type 2 astrocytes. The fractal dimension remained constant over a 10-fold range in optical magnification, demonstrating that cultured O2A glial cells exhibit self-similarity, a defining characteristic of fractal objects. These results illustrate that fractal dimension analysis of maturing cell populations is a useful method for quantitatively describing the process of cell differentiation.

Use of fractal theory in neuroscience: methods, advantages, and potential problems

Fractal analysis has already found widespread application in the field of neuroscience and is being used in many other areas. Applications are many and include ion channel kinetics of biological membranes and classification of neurons according to their branching characteristics. In this article we review some practical methods that are now available to allow the determination of the complexity and scaling relationships of anatomical and physiological patterns. The problems of describing fractal dimensions are discussed and the concept of fractal dimensionality is introduced. Several related methodological considerations, such as preparation of the image and estimation of the fractal dimensions from the data points, as well as the advantages and problems of fractal geometric analysis, are discussed.

Morphology of reactive microglia in the injured cerebral cortex. Fractal analysis and complementary quantitative methods

The present study focuses on application of quantitative methods measuring differences between particular morphological types of microglial cells as well as between their proliferating and non-proliferating examples. On the basis of subjective classification, microglial cells of three morphological types (ramified, hypertrophied and bushy) were selected from the neocortex of injured rat brain. Thereafter, the morphological complexity of each cell was assessed by calculation its fractal dimension as well as its form factor, convexity, ramification factor and solidity. The fractal dimension seemed a good parameter for detecting small changes in the space-filing capacity of cells, for example, it shows differences between ramified cells from control and injured brains. This measure seemed insensitive to some aspects of cell morphology. To obtain precise quantification of observed changes other morphological parameters had to be applied. Proliferating and non-proliferating microglial cells displayed significant differences in their solidity and ramification factors, but not in fractal dimension and convexity. The results indicated that proliferating microglia were more massive and less-ramified but they did not reduce their spatial complexity.

Neurons and fractals: how reliable and useful are calculations of fractal dimensions?

In the past 15 years it has become possible to determine the fractal dimension (Df) of complex objects, including neurons, by automated image analysis methods. However, there are many unresolved issues that need to be addressed. In this paper we discuss how the Df calculated by different methods may vary and how fractal analysis may be of use for retinal ganglion cell characterization. The goal of this work was to acknowledge inherent sources of variation during measurement and evaluate current fractal analysis methods for describing structure. Our results show that different algorithms and even the same algorithm performed by different computer programs and/or experimenters may give different but consistent numerical values. All described methods demonstrated their suitability for classifying cat retinal ganglion cells into distinct groups. Our results reinforce the idea that comparison of measurements of different profiles using the same measurement method may be useful and valid even if an exact numeric value of the dimension is not realised in practice.

Fractal methods and results in cellular morphology--dimensions, lacunarity and multifractals

This paper discusses the concepts of fractal geometry in a cellular biological context. It defines the concept of the fractal dimension. D, as a measure of complexity and illustrates the two different general ways of quantitatively measuring D by length-related and mass-related methods. Then, these several Ds are compared and contrasted. A goal of the paper is to find methods other than length-related measures that can distinguish between two objects that have the same D but are structurally different. The mass-related D is shown potentially to be such a measure. The concept of lacunarity, L, is defined and methods of measuring L are illustrated. L is also shown to be a potentially distinguishing measure. Finally, the notion of multifracticality is defined and illustrated to exist in certain individual nerve and glial cells.

Fractal analysis of cytoskeleton rearrangement in cardiac muscle during head-down tilt

Head-down tilt by tail suspension of the rat produces a volume, but not pressure, load on the heart. One response of the heart is cytoskeleton rearrangement, a phenomenon commonly referred to as disruption. In these experiments, we used fractal analysis as a means to measure complexity of the microtubule structures at 8 and 18 h after imposition of head-down tilt. Microtubules in whole tissue cardiac myocytes were stained with fluorescein colchicine and were visualized by confocal microscopy. The fractal dimensions (D) of the structures were calculated by the dilation method, which involves successively dilating the outline perimeter of the microtubule structures and measuring the area enclosed. The head-down tilt resulted in a progressive decrease in D (decreased complexity) when measured at small dilations of the perimeter, but the maximum D (maximum complexity) of the microtubule structures did not change with treatment. Analysis of the fold change in complexity as a function of the dilation indicates an almost twofold decrease in microtubule complexity at small kernel dilations. This decrease in complexity is associated with a more Gaussian distribution of microtubule diameters, indicating a less structured microtubule cytoskeleton. We interpret these data as a microtubule rearrangement, rather than erosion, because total tubulin flourescence was not different between groups. This conclusion is supported by F-actin fluorescence data indicating a dispersed structure without loss of actin.

Complexity analysis of oligodendroglial processes expressing myelin-associated glycoprotein

Oligodendroglia synthesize myelin in the mammalian central nervous system. Mature oligodendroglia have been identified in culture by two criteria; the expression of molecules characteristic of myelin, such as galactocerebroside (galC) and myelin-associated glycoprotein (MAG), and the elaboration of complex processes. Myelin gene expression can be documented by the binding of specific antibodies and antisera to the myelin-specific molecules; process complexity can be described by the fractal dimension, D. In this study, anti-MAG antisera was used to document MAG expression in the processes of oligodendroglia. Eighty percent of the galC+ oligodendroglia bound anti-MAG antiserum. With time in culture, MAG immunoreactivity seemed to extend from the cell soma into the oligodendroglial processes. To quantify this observation, fractal dimensions were calculated using either galC or MAG immunoreactivity to visualize oligodendroglial processes. A fractal dimension of 1.5 was calculated for O1+ processes by day 4 of culture; this value for D remained constant over the course of 1 month in culture. The fractal dimension calculated for MAG+ processes increased from 1.2 to 1.5 over the course of 28 days in culture. This change in fractal dimension confirms our visual impression that galC-containing processes acquire MAG slowly over the course of several weeks in culture.

Retinoic acid regulates the development of oligodendrocyte precursor cells in vitro

Cultures of oligodendrocyte precursor cells can be grown from brain hemispheres of newborn rats. These cells, also called O-2A progenitor cells, can differentiate in vitro into oligodendrocytes or type 2 astrocytes. Basic FGF and PDGF are known to stimulate their proliferation and delay their differentiation. Lack or excess of retinoic acid (RA) has been known for a long time to alter brain development suggesting that this compound is involved in normal brain development. Here we report that RA partially inhibits both the proliferation and the differentiation of oligodendrocyte precursor cells. It also down-regulates the mitogenic effect of bFGF on these cells while keeping them in an immature stage. RA is more effective than bFGF in inhibiting myelin basic protein mRNA expression in these cells, and like bFGF, it preserves their bipotential character. RA nuclear receptors RAR-alpha and their transcripts are expressed in oligodendrocyte precursor cells as seen by Western blot, Northern blot and in situ hybridization. The expression of RAR-alpha transcripts is stimulated transiently by RA alone or associated to bFGF. The expression of RAR-beta transcripts is not constitutive and is induced by RA alone or associated to bFGF and to a lesser extent by bFGF alone. These results suggest that retinoids participate in the control of the development of glial cells of the oligodendrocyte lineage.