The usefulness of spatial pattern analysis in understanding the pathogenesis of neurodegenerative disorders, with particular reference to plaque formation in Alzheimer's disease

Discrete, microscopic lesions are developed in the brain in a number of neurodegenerative diseases. These lesions may not be randomly distributed in the tissue but exhibit a spatial pattern, i.e., a departure from randomness towards regularlity or clustering. The spatial pattern of a lesion may reflect its development in relation to other brain lesions or to neuroanatomical structures. Hence, a study of spatial pattern may help to elucidate the pathogenesis of a lesion. A number of statistical methods can be used to study the spatial patterns of brain lesions. They range from simple tests of whether the distribution of a lesion departs from random to more complex methods which can detect clustering and the size, distribution and spacing of clusters. This paper reviews the uses and limitations of these methods as applied to neurodegenerative disorders, and in particular to senile plaque formation in Alzheimer's disease.

The relationship between the spatial pattern of senile plaques and neurofibrillary tangles in Alzheimer's disease

The topographic pattern of senile plaques (SP) and neurofibrillary tangles (NFT) was studied in silver stained coronal sections of neocortex and hippocampus in ten cases of Alzheimer's disease (AD). Both lesions showed evidence of clustering in the tissue with many of the clusters being regularly spaced. The patterns of SP and NFT were compared 1) in the same cortical zone, 2) between upper and lower zones of the cortex and 3) in regions connected by either association fibres or the perforant path. Correlations between the lesions in the same cortical zone were found in 20% of the layers examined while correlations between upper and lower zones occurred in 64% of cortical regions examined. There was evidence that NFT in upper and lower cortex may be in register in some tissues. In addition, positive correlations were found between upper NFT and lower SP and negative correlations between upper SP and lower NFT in some tissues. Regular clustering of lesions was also observed in brain regions connected to one another suggesting that they develop on functinally related sets of neurons.

A spatial pattern analysis of beta-amyloid (Abeta) deposition in the temporal lobe in Alzheimer's disease

To determine the factors influencing the distribution of -amyloid (Abeta) deposits in Alzheimer's disease (AD), the spatial patterns of the diffuse, primitive, and classic A deposits were studied from the superior temporal gyrus (STG) to sector CA4 of the hippocampus in six sporadic cases of the disease. In cortical gyri and in the CA sectors of the hippocampus, the Abeta deposits were distributed either in clusters 200-6400 microm in diameter that were regularly distributed parallel to the tissue boundary or in larger clusters greater than 6400 microm in diameter. In some regions, smaller clusters of Abeta deposits were aggregated into larger 'superclusters'. In many cortical gyri, the density of Abeta deposits was positively correlated with distance below the gyral crest. In the majority of regions, clusters of the diffuse, primitive, and classic deposits were not spatially correlated with each other. In two cases, double immunolabelled to reveal the Abeta deposits and blood vessels, the classic Abeta deposits were clustered around the larger diameter vessels. These results suggest a complex pattern of Abeta deposition in the temporal lobe in sporadic AD. A regular distribution of Abeta deposit clusters may reflect the degeneration of specific cortico-cortical and cortico-hippocampal pathways and the influence of the cerebral blood vessels. Large-scale clustering may reflect the aggregation of deposits in the depths of the sulci and the coalescence of smaller clusters.

Fast algorithms for automatic mapping with space-limited covariance functions

In this paper we discuss a fast Bayesian extension to kriging algorithms which has been used successfully for fast, automatic mapping in emergency conditions in the Spatial Interpolation Comparison 2004 (SIC2004) exercise. The application of kriging to automatic mapping raises several issues such as robustness, scalability, speed and parameter estimation. Various ad-hoc solutions have been proposed and used extensively but they lack a sound theoretical basis. In this paper we show how observations can be projected onto a representative subset of the data, without losing significant information. This allows the complexity of the algorithm to grow as O(n m 2), where n is the total number of observations and m is the size of the subset of the observations retained for prediction. The main contribution of this paper is to further extend this projective method through the application of space-limited covariance functions, which can be used as an alternative to the commonly used covariance models. In many real world applications the correlation between observations essentially vanishes beyond a certain separation distance. Thus it makes sense to use a covariance model that encompasses this belief since this leads to sparse covariance matrices for which optimised sparse matrix techniques can be used. In the presence of extreme values we show that space-limited covariance functions offer an additional benefit, they maintain the smoothness locally but at the same time lead to a more robust, and compact, global model. We show the performance of this technique coupled with the sparse extension to the kriging algorithm on synthetic data and outline a number of computational benefits such an approach brings. To test the relevance to automatic mapping we apply the method to the data used in a recent comparison of interpolation techniques (SIC2004) to map the levels of background ambient gamma radiation. © Springer-Verlag 2007.

A morphometric study of the spatial patterns of TDP-43 immunoreactive neuronal inclusions in frontotemporal lobar degeneration (FTLD) with progranulin (GRN) mutation

Mutations of the progranulin (GRN) gene are a major cause of familial frontotemporal lobar degeneration with transactive response (TAR) DNA-binding protein of 43 kDa (TDP-43) proteinopathy (FTLD-TDP). We studied the spatial patterns of TDP-43 immunoreactive neuronal cytoplasmic inclusions (NCI) and neuronal intranuclear inclusions (NII) in histological sections of the frontal and temporal lobe in eight cases of FTLD-TDP with GRN mutation using morphometric methods and spatial pattern analysis. In neocortical regions, the NCI were clustered and the clusters were regularly distributed parallel to the pia mater; 58% of regions analysed exhibiting this pattern. The NII were present in regularly distributed clusters in 35% of regions but also randomly distributed in many areas. In neocortical regions, the sizes of the regular clusters of NCI and NII were 400-800 µm, approximating to the size of the modular columns of the cortico-cortical projections, in 31% and 36% of regions respectively. The NCI and NII also exhibited regularly spaced clustering in sectors CA1/2 of the hippocampus and in the dentate gyrus. The clusters of NCI and NII were not spatially correlated. The data suggest degeneration of the cortico-cortical and cortico-hippocampal pathways in FTLD-TDP with GRN mutation, the NCI and NII affecting different clusters of neurons.

Parallel implementation of image segmentation algorithms on a network of transputers

Image segmentation is one of the most computationally intensive operations in image processing and computer vision. This is because a large volume of data is involved and many different features have to be extracted from the image data. This thesis is concerned with the investigation of practical issues related to the implementation of several classes of image segmentation algorithms on parallel architectures. The Transputer is used as the basic building block of hardware architectures and Occam is used as the programming language. The segmentation methods chosen for implementation are convolution, for edge-based segmentation; the Split and Merge algorithm for segmenting non-textured regions; and the Granlund method for segmentation of textured images. Three different convolution methods have been implemented. The direct method of convolution, carried out in the spatial domain, uses the array architecture. The other two methods, based on convolution in the frequency domain, require the use of the two-dimensional Fourier transform. Parallel implementations of two different Fast Fourier Transform algorithms have been developed, incorporating original solutions. For the Row-Column method the array architecture has been adopted, and for the Vector-Radix method, the pyramid architecture. The texture segmentation algorithm, for which a system-level design is given, demonstrates a further application of the Vector-Radix Fourier transform. A novel concurrent version of the quad-tree based Split and Merge algorithm has been implemented on the pyramid architecture. The performance of the developed parallel implementations is analysed. Many of the obtained speed-up and efficiency measures show values close to their respective theoretical maxima. Where appropriate comparisons are drawn between different implementations. The thesis concludes with comments on general issues related to the use of the Transputer system as a development tool for image processing applications; and on the issues related to the engineering of concurrent image processing applications.

The use of Meteosat satellite data for spatial rainfall estimations and hydrological simulations

Satellite information, in combination with conventional point source measurements, can be a valuable source of information. This thesis is devoted to the spatial estimation of areal rainfall over a region using both the measurements from a dense and sparse network of rain-gauges and images from the meteorological satellites. A primary concern is to study the effects of such satellite assisted rainfall estimates on the performance of rainfall-runoff models. Low-cost image processing systems and peripherals are used to process and manipulate the data. Both secondary as well as primary satellite images were used for analysis. The secondary data was obtained from the in-house satellite receiver and the primary data was obtained from an outside source. Ground truth data was obtained from the local Water Authority. A number of algorithms are presented that combine the satellite and conventional data sources to produce areal rainfall estimates and the results are compared with some of the more traditional methodologies. The results indicate that the satellite cloud information is valuable in the assessment of the spatial distribution of areal rainfall, for both half-hourly as well as daily estimates of rainfall. It is also demonstrated how the performance of the simple multiple regression rainfall-runoff model is improved when satellite cloud information is used as a separate input in addition to rainfall estimates from conventional means. The use of low-cost equipment, from image processing systems to satellite imagery, makes it possible for developing countries to introduce such systems in areas where the benefits are greatest.

Density and spatial pattern of β-amyloid (Aβ) deposits in corticobasal degeneration

Corticobasal degeneration (CBD) is a rare, progressive movement disorder characterized neuropathologically by widespread neuronal and glial pathology including tau-immunoreactive neuronal cytoplasmic inclusions (NCI), oligodendroglial inclusions (GI), and astrocytic plaques (AP). However, ß -amyloid (A ß) deposits have been observed in the cerebral cortex and/or hippocampus in some cases of CBD. To clarify the role of Aß deposition in CBD, the densities and spatial patterns of the Aß deposits were studied in three cases. In two cases, expressing apolipoprotein E (APOE) genotypes 2/3 or 3/3, the densities of the Aß deposits were similar to those in normal elderly brain. In the remaining case, expressing APOE genotype 3/4, Aß deposition was observed throughout the cerebral cortex, sectors CA1 and CA2 of the hippocampus, and the molecular layer of the dentate gyrus. The densities of the Aß deposits in this case were typical of those observed in Alzheimer's disease (AD). In the three cases, clustering of Aß deposits, with clusters ranging in size from 200 to >6400 µm in diameter, was evident in 25/27 (93%) of analyses. In addition, the clusters of Aß deposits were regularly distributed parallel to the tissue boundary in 52% of analyses, a spatial pattern similar to that observed in AD. These results suggest: (1) in some CBD cases, Aß pathology is age-related, (2) more extensive Aß deposition is observed in some cases, the density and spatial patterns of the Aß deposits being similar to AD, and (3) extensive deposition of Aß in CBD may be associated with APOE allele e4.

Exploring ant-based algorithms for gene expression data analysis

Objective: Recently, much research has been proposed using nature inspired algorithms to perform complex machine learning tasks. Ant colony optimization (ACO) is one such algorithm based on swarm intelligence and is derived from a model inspired by the collective foraging behavior of ants. Taking advantage of the ACO in traits such as self-organization and robustness, this paper investigates ant-based algorithms for gene expression data clustering and associative classification. Methods and material: An ant-based clustering (Ant-C) and an ant-based association rule mining (Ant-ARM) algorithms are proposed for gene expression data analysis. The proposed algorithms make use of the natural behavior of ants such as cooperation and adaptation to allow for a flexible robust search for a good candidate solution. Results: Ant-C has been tested on the three datasets selected from the Stanford Genomic Resource Database and achieved relatively high accuracy compared to other classical clustering methods. Ant-ARM has been tested on the acute lymphoblastic leukemia (ALL)/acute myeloid leukemia (AML) dataset and generated about 30 classification rules with high accuracy. Conclusions: Ant-C can generate optimal number of clusters without incorporating any other algorithms such as K-means or agglomerative hierarchical clustering. For associative classification, while a few of the well-known algorithms such as Apriori, FP-growth and Magnum Opus are unable to mine any association rules from the ALL/AML dataset within a reasonable period of time, Ant-ARM is able to extract associative classification rules.

Spatial patterns of the tau pathology in progressive supranuclear palsy

Progressive supranuclear palsy (PSP) is characterized neuropathologically by neuronal loss, gliosis, and the presence of tau-immunoreactive neuronal and glial cell inclusions affecting subcortical and some cortical regions. The objectives of this study were to determine (1) the spatial patterns of the tau-immunoreactive pathology, viz., neurofibrillary tangles (NFT), oligodendroglial inclusions (GI), tufted astrocytes (TA), and Alzheimer's disease-type neuritic plaques (NP) in PSP and (2) to investigate the spatial correlations between the histological features. Post-mortem material of cortical and subcortical regions of eight PSP cases was studied. Spatial pattern analysis was applied to the NFT, GI, TA, NP, abnormally enlarged neurons (EN), surviving neurons, and glial cells. NFT, GI, and TA were distributed either at random or in regularly distributed clusters. The EN and NP were mainly randomly distributed. Clustering of NFT and EN was more frequent in the cortex and subcortical regions, respectively. Variations in NFT density were not spatially correlated with the densities of either GI or TA, but were positively correlated with the densities of EN and surviving neurons in some regions. (1) NFT were the most widespread tau-immunoreactive pathology in PSP being distributed randomly in subcortical regions and in regular clusters in cortical regions, (2) GI and TA were more localized and exhibited a regular pattern of clustering in subcortical regions, and (3) neuronal and glial cell pathologies were not spatially correlated. © 2012 Springer-Verlag.

A multiple pheromone ant clustering algorithm

Ant colony optimisation algorithms model the way ants use pheromones for marking paths to important locations in their environment. Pheromone traces are picked up, followed, and reinforced by other ants but also evaporate over time. Optimal paths attract more pheromone and less useful paths fade away. The main innovation of the proposed Multiple Pheromone Ant Clustering Algorithm (MPACA) is to mark objects using many pheromones, one for each value of each attribute describing the objects in multidimensional space. Every object has one or more ants assigned to each attribute value and the ants then try to find other objects with matching values, depositing pheromone traces that link them. Encounters between ants are used to determine when ants should combine their features to look for conjunctions and whether they should belong to the same colony. This paper explains the algorithm and explores its potential effectiveness for cluster analysis. © 2014 Springer International Publishing Switzerland.

The multiple pheromone Ant clustering algorithm

Ant Colony Optimisation algorithms mimic the way ants use pheromones for marking paths to important locations. Pheromone traces are followed and reinforced by other ants, but also evaporate over time. As a consequence, optimal paths attract more pheromone, whilst the less useful paths fade away. In the Multiple Pheromone Ant Clustering Algorithm (MPACA), ants detect features of objects represented as nodes within graph space. Each node has one or more ants assigned to each feature. Ants attempt to locate nodes with matching feature values, depositing pheromone traces on the way. This use of multiple pheromone values is a key innovation. Ants record other ant encounters, keeping a record of the features and colony membership of ants. The recorded values determine when ants should combine their features to look for conjunctions and whether they should merge into colonies. This ability to detect and deposit pheromone representative of feature combinations, and the resulting colony formation, renders the algorithm a powerful clustering tool. The MPACA operates as follows: (i) initially each node has ants assigned to each feature; (ii) ants roam the graph space searching for nodes with matching features; (iii) when departing matching nodes, ants deposit pheromones to inform other ants that the path goes to a node with the associated feature values; (iv) ant feature encounters are counted each time an ant arrives at a node; (v) if the feature encounters exceed a threshold value, feature combination occurs; (vi) a similar mechanism is used for colony merging. The model varies from traditional ACO in that: (i) a modified pheromone-driven movement mechanism is used; (ii) ants learn feature combinations and deposit multiple pheromone scents accordingly; (iii) ants merge into colonies, the basis of cluster formation. The MPACA is evaluated over synthetic and real-world datasets and its performance compares favourably with alternative approaches.

Spatial ecological complexity measures in GRASS GIS

Good estimates of ecosystem complexity are essential for a number of ecological tasks: from biodiversity estimation, to forest structure variable retrieval, to feature extraction by edge detection and generation of multifractal surface as neutral models for e.g. feature change assessment. Hence, measuring ecological complexity over space becomes crucial in macroecology and geography. Many geospatial tools have been advocated in spatial ecology to estimate ecosystem complexity and its changes over space and time. Among these tools, free and open source options especially offer opportunities to guarantee the robustness of algorithms and reproducibility. In this paper we will summarize the most straightforward measures of spatial complexity available in the Free and Open Source Software GRASS GIS, relating them to key ecological patterns and processes.

A spatial pattern analysis of β-amyloid (Aβ) deposition in the temporal lobe in Alzheimer's disease

To determine the factors influencing the distribution of β-amyloid (Aβ) deposits in Alzheimer's disease (AD), the spatial patterns of the diffuse, primitive, and classic Aβ deposits were studied from the superior temporal gyrus (STG) to sector CA4 of the hippocampus in six sporadic cases of the disease. In cortical gyri and in the CA sectors of the hippocampus, the Aβ deposits were distributed either in clusters 200-6400 μm in diameter that were regularly distributed parallel to the tissue boundary or in larger clusters greater than 6400 μm in diameter. In some regions, smaller clusters of Aβ deposits were aggregated into larger 'superclusters'. In many cortical gyri, the density of Aβ deposits was positively correlated with distance below the gyral crest. In the majority of regions, clusters of the diffuse, primitive, and classic deposits were not spatially correlated with each other. In two cases, double immunolabelled to reveal the Aβ deposits and blood vessels, the classic Aβ deposits were clustered around the larger diameter vessels. These results suggest a complex pattern of Aβ deposition in the temporal lobe in sporadic AD. A regular distribution of Aβ deposit clusters may reflect the degeneration of specific cortico-cortical and cortico-hippocampal pathways and the influence of the cerebral blood vessels. Large-scale clustering may reflect the aggregation of deposits in the depths of the sulci and the coalescence of smaller clusters.

Clustering of tau-immunoreactive pathology in chronic traumatic encephalopathy

Chronic traumatic encephalopathy (CTE) is a neurodegenerative disorder which may result from repetitive brain injury. A variety of tau-immunoreactive pathologies are present, including neurofibrillary tangles (NFT), neuropil threads (NT), dot-like grains (DLG), astrocytic tangles (AT), and occasional neuritic plaques (NP). In tauopathies, cellular inclusions in the cortex are clustered within specific laminae, the clusters being regularly distributed parallel to the pia mater. To determine whether a similar spatial pattern is present in CTE, clustering of the tau-immunoreactive pathology was studied in the cortex, hippocampus, and dentate gyrus in 11 cases of CTE and 7 cases of Alzheimer’s disease neuropathologic change (ADNC) without CTE. In CTE: (1) all aspects of tau-immunoreactive pathology were clustered and the clusters were frequently regularly distributed parallel to the tissue boundary, (2) clustering was similar in two CTE cases with minimal co-pathology compared with cases with associated ADNC or TDP-43 proteinopathy, (3) in a proportion of cortical gyri, estimated cluster size was similar to that of cell columns of the cortico-cortical pathways, and (4) clusters of the tau-immunoreactive pathology were infrequently spatially correlated with blood vessels. The NFT and NP in ADNC without CTE were less frequently randomly or uniformly distributed and more frequently in defined clusters than in CTE. Hence, the spatial pattern of the tau-immunoreactive pathology observed in CTE is typical of the tauopathies but with some distinct differences compared to ADNC alone. The spread of pathogenic tau along anatomical pathways could be a factor in the pathogenesis of the disease.