BIO-INSPIRED COMPUTATIONAL TECHNIQUES APPLIED TO THE CLUSTERING AND VISUALIZATION OF SPATIO-TEMPORAL GEOSPATIAL DATA

The coverage and volume of geo-referenced datasets are extensive and incessantly¦growing. The systematic capture of geo-referenced information generates large volumes¦of spatio-temporal data to be analyzed. Clustering and visualization play a key¦role in the exploratory data analysis and the extraction of knowledge embedded in¦these data. However, new challenges in visualization and clustering are posed when¦dealing with the special characteristics of this data. For instance, its complex structures,¦large quantity of samples, variables involved in a temporal context, high dimensionality¦and large variability in cluster shapes.¦The central aim of my thesis is to propose new algorithms and methodologies for¦clustering and visualization, in order to assist the knowledge extraction from spatiotemporal¦geo-referenced data, thus improving making decision processes.¦I present two original algorithms, one for clustering: the Fuzzy Growing Hierarchical¦Self-Organizing Networks (FGHSON), and the second for exploratory visual data analysis:¦the Tree-structured Self-organizing Maps Component Planes. In addition, I present¦methodologies that combined with FGHSON and the Tree-structured SOM Component¦Planes allow the integration of space and time seamlessly and simultaneously in¦order to extract knowledge embedded in a temporal context.¦The originality of the FGHSON lies in its capability to reflect the underlying structure¦of a dataset in a hierarchical fuzzy way. A hierarchical fuzzy representation of¦clusters is crucial when data include complex structures with large variability of cluster¦shapes, variances, densities and number of clusters. The most important characteristics¦of the FGHSON include: (1) It does not require an a-priori setup of the number¦of clusters. (2) The algorithm executes several self-organizing processes in parallel.¦Hence, when dealing with large datasets the processes can be distributed reducing the¦computational cost. (3) Only three parameters are necessary to set up the algorithm.¦In the case of the Tree-structured SOM Component Planes, the novelty of this algorithm¦lies in its ability to create a structure that allows the visual exploratory data analysis¦of large high-dimensional datasets. This algorithm creates a hierarchical structure¦of Self-Organizing Map Component Planes, arranging similar variables' projections in¦the same branches of the tree. Hence, similarities on variables' behavior can be easily¦detected (e.g. local correlations, maximal and minimal values and outliers).¦Both FGHSON and the Tree-structured SOM Component Planes were applied in¦several agroecological problems proving to be very efficient in the exploratory analysis¦and clustering of spatio-temporal datasets.¦In this thesis I also tested three soft competitive learning algorithms. Two of them¦well-known non supervised soft competitive algorithms, namely the Self-Organizing¦Maps (SOMs) and the Growing Hierarchical Self-Organizing Maps (GHSOMs); and the¦third was our original contribution, the FGHSON. Although the algorithms presented¦here have been used in several areas, to my knowledge there is not any work applying¦and comparing the performance of those techniques when dealing with spatiotemporal¦geospatial data, as it is presented in this thesis.¦I propose original methodologies to explore spatio-temporal geo-referenced datasets¦through time. Our approach uses time windows to capture temporal similarities and¦variations by using the FGHSON clustering algorithm. The developed methodologies¦are used in two case studies. In the first, the objective was to find similar agroecozones¦through time and in the second one it was to find similar environmental patterns¦shifted in time.¦Several results presented in this thesis have led to new contributions to agroecological¦knowledge, for instance, in sugar cane, and blackberry production.¦Finally, in the framework of this thesis we developed several software tools: (1)¦a Matlab toolbox that implements the FGHSON algorithm, and (2) a program called¦BIS (Bio-inspired Identification of Similar agroecozones) an interactive graphical user¦interface tool which integrates the FGHSON algorithm with Google Earth in order to¦show zones with similar agroecological characteristics.

Spatial, temporal and subcellular localization of islet-brain 1 (IB1), a homologue of JIP-1, in mouse brain.

Islet-brain 1 (IB1) was recently identified as a DNA-binding protein of the GLUT2 gene promoter. The mouse IB1 is the rat and human homologue of the Jun-interacting protein 1 (JIP-1) which has been recognized as a key player in the regulation of c-Jun amino-terminal kinase (JNK) mitogen-activated protein kinase (MAPK) pathways. JIP-1 is involved in the control of apoptosis and may play a role in brain development and aging. Here, IB1 was studied in adult and developing mouse brain tissue by in situ hybridization, Northern and Western blot analysis at cellular and subcellular levels, as well as by immunocytochemistry in brain sections and cell cultures. IB1 expression was localized in the synaptic regions of the olfactory bulb, retina, cerebral and cerebellar cortex and hippocampus in the adult mouse brain. IB1 was also detected in a restricted number of axons, as in the mossy fibres from dentate gyrus in the hippocampus, and was found in soma, dendrites and axons of cerebellar Purkinje cells. After birth, IB1 expression peaks at postnatal day 15. IB1 was located in axonal and dendritic growth cones in primary telencephalon cells. By biochemical and subcellular fractionation of neuronal cells, IB1 was detected both in the cytosolic and membrane fractions. Taken together with previous data, the restricted neuronal expression of IB1 in developing and adult brain and its prominent localization in synapses suggest that the protein may be critical for cell signalling in developing and mature nerve terminals.

Characterization of soil chemical properties of strawberry fields using principal component analysis

One of the largest strawberry-producing municipalities of Rio Grande do Sul (RS) is Turuçu, in the South of the State. The strawberry production system adopted by farmers is similar to that used in other regions in Brazil and in the world. The main difference is related to the soil management, which can change the soil chemical properties during the strawberry cycle. This study had the objective of assessing the spatial and temporal distribution of soil fertility parameters using principal component analysis (PCA). Soil sampling was based on topography, dividing the field in three thirds: upper, middle and lower. From each of these thirds, five soil samples were randomly collected in the 0-0.20 m layer, to form a composite sample for each third. Four samples were taken during the strawberry cycle and the following properties were determined: soil organic matter (OM), soil total nitrogen (N), available phosphorus (P) and potassium (K), exchangeable calcium (Ca) and magnesium (Mg), soil pH (pH), cation exchange capacity (CEC) at pH 7.0, soil base (V%) and soil aluminum saturation(m%). No spatial variation was observed for any of the studied soil fertility parameters in the strawberry fields and temporal variation was only detected for available K. Phosphorus and K contents were always high or very high from the beginning of the strawberry cycle, while pH values ranged from very low to very high. Principal component analysis allowed the clustering of all strawberry fields based on variables related to soil acidity and organic matter content.

Outcome evaluation of ankle osteoarthritis treatments: plantar pressure analysis during relatively long-distance walking.

Background : This study aimed to use plantar pressure analysis in relatively long-distance walking for objective outcome evaluation of ankle osteoarthritis treatments, i.e., ankle arthrodesis and total ankle replacement.Methods : Forty-seven subjects in four groups: three patient groups and controls, participated in the study. Each subject walked twice in 50-m trials. Plantar pressure under the pathological foot was measured using pressure insoles. Six parameters: initial contact time, terminal contact time, maximum force time, peak pressure time, maximum force and peak pressure were calculated and averaged over trials in ten regions of foot. The parameters in each region were compared between patient groups and controls and their effect size was estimated. Besides, the correlations between pressure parameters and clinical scales were calculated.Findings : We observed based on temporal parameters that patients postpone the heel-off event, when high force in forefoot and high ankle moment happens. Also based on maximum force and peak pressure, the patients apply smoothened maximum forces on the affected foot. In ten regions, some parameters showed improvements after total ankle replacement, some showed alteration of foot function after ankle arthrodesis and some others showed still abnormality after both surgical treatments. These parameters showed also significant correlation with clinical scales in at least two regions of foot.Interpretation : Plantar pressure parameters in relatively long-distance trials showed to be strong tools for outcome evaluation of ankle osteoarthritis treatments. (C) 2010 Elsevier Ltd. All rights reserved.

Spatial analysis of erythrocyte membrane fluctuations by digital holographic microscopy.

Red blood cell (RBC) membrane fluctuations provide important insights into cell states. We present a spatial analysis of red blood cell membrane fluctuations by using digital holographic microscopy (DHM). This interferometric and dye-free technique, possessing nanometric axial and microsecond temporal sensitivities enables to measure cell membrane fluctuations (CMF) on the whole cell surface. DHM acquisition is combined with a model which allows extracting the membrane fluctuation amplitude, while taking into account cell membrane topology. Uneven distribution of CMF amplitudes over the RBC surface is observed, showing maximal values in a ring corresponding to the highest points on the RBC torus as well as in some scattered areas in the inner region of the RBC. CMF amplitudes of 35.9+/-8.9 nm and 4.7+/-0.5 nm (averaged over the cell surface) were determined for normal and ethanol-fixed RBCs, respectively.

Comparative genomics of the vertebrate insulin/TOR signal transduction pathway: A network-level analysis of selective pressures

Complexity of biological function relies on large networks of interacting molecules. However, the evolutionary properties of these networks are not fully understood. It has been shown that selective pressures depend on the position of genes in the network. We have previously shown that in the Drosophila insulin/target of rapamycin (TOR) signal transduction pathway there is a correlation between the pathway position and the strength of purifying selection, with the downstream genes being most constrained. In this study, we investigated the evolutionary dynamics of this well-characterized pathway in vertebrates. More specifically, we determined the impact of natural selection on the evolution of 72 genes of this pathway. We found that in vertebrates there is a similar gradient of selective constraint in the insulin/TOR pathway to that found in Drosophila. This feature is neither the result of a polarity in the impact of positive selection nor of a series of factors affecting selective constraint levels (gene expression level and breadth, codon bias, protein length, and connectivity). We also found that pathway genes encoding physically interacting proteins tend to evolve under similar selective constraints. The results indicate that the architecture of the vertebrate insulin/TOR pathway constrains the molecular evolution of its components. Therefore, the polarity detected in Drosophila is neither specific nor incidental of this genus. Hence, although the underlying biological mechanisms remain unclear, these may be similar in both vertebrates and Drosophila.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

Analysis of the mechanisms controlling the activity of flp1, a chromosomal passenger protein in the fission Schizosaccharomyces pombe

ABSTRACT In S. cerevisiae, the protein phosphatase Cdc14pwt is essential far mitotic exit through its contribution to reducing mitotic CDK activity. But Cdc14pwt also acts as a mare general temporal coordinator of mid and late mitotic events by controlling the partitioning of DNA, microtubule stability and cytokinesis. Cdc14pwt orthologs are well conserved from yeasts to humans, and sequence comparison revealed the presence of three domains, A, B and C, of which A and B form the catalytic domain. Cdc14pwt orthologs are regulated (in part) through cell cycle dependent changes in their localization. Some of them are thought to be kept inactive by sequestration in the nucleolus during interphase. This is the case for flp1pwt, the single identified Cdc14pwt ortholog in the fission yeast S. pombe. In early mitosis, flp1pwt leaves the nucleolus and localizes to the kinetochores, the contractile ring and the mitotic spindle, suggesting that it has multiple substrates and regulates many mitotic processes. flp1D cells show a high chromosome loss rate and septation defects, suggesting a role for flp1wt in the fidelity of chromosome transmission and cytokinesis. The aim of this study is to characterize the mechanisms underlying flp1pwt functions and the control of its activity. A structure-function analysis has revealed that the presence of both A and B domains is required for biological function and for proper flp1pwt mitotic localization. In contrast, the C domain of flp1pwt is responsible for its proper nucleolar localization in G2/interphase. My data suggest that dephosphorylation of substrates by flp1pwt is not necessary for any changes in localization of flp1pwt except that at the medial ring. In that particular case, the catalytic activity of flp1pwt is required for efficient localization, therefore revealing an additional level of regulation. All the functions of flp1pwt assayed to date require its catalytic activity, emphasizing the importance of further identification of its substrates. As described for other orthologs, the capability of selfinteraction and phosphorylation status might help to control flp1pwt activity. My data suggest that flp1pwt forms oligomers in vivo and that phosphorylation is not essential far localization changes of the protein. In addition, the hypophosphorylated form of flp1pwt might be specifically involved in the promotion of cytokinesis. The results of this study suggest that multiple modes of regulation including localization, selfassociation and phosphorylation allow a fine-tuning regulation of flp1pwt phosphatase activity, and more generally that of Cdc14pwt family of phosphatases. RESUME Chez la levure S. cerevisiae, la protéine phosphatase Cdc14pwt est essentielle pour la sortie de mitose du fait de sa contribution dans la réduction d'activité des CDK mitotiques. Comme elle contrôle également le partage de l'ADN, la stabilité des microtubules et la cytokinèse, Cdc14pwt est en fait considérée comme un coordinateur temporel général des évènements de milieu et de fin de mitose. Les orthologues de Cdc14pwt sont bien conservés, des levures jusqu'à l'espèce humaine. Des comparaisons de séquence ont révélé la présence de trois domaines A, B et C, les deux premiers constituant le domaine catalytique. Ils sont régulés (en partie) via des changements dans leur localisation, eux-mêmes dépendants du cycle cellulaire. Plusieurs de ces orthologues sont supposés inactivés par séquestration dans le nucléole en interphase, ce qui est le cas de flp1pwt le seul orthologue de Cdc14pwt identifié chez la levure fissipare S, pombe. En début de mitose, flp1pwt quitte le nucléole et localise au niveau des kinetochores, de l'anneau contractile d'actine et du fuseau mitotique, ce qui laisse supposer de multiples substrats et fonctions. Comme les cellules délétées pour le gène flp1wt présentent un taux élevé de perte de chromosome et des défauts de septation, flp1pwt semble jouer un rôle dans la fidélité de la transmission du matériel génétique et la cytokinèse. Le but de cette étude est de caractériser les mécanismes impliqués dans les fonctions assurées par flp1pwt d'une part, et dans le contrôle de son activité d'autre part. Une analyse structure-fonction a révélé que la présence simultanée des deux domaines A et B est requise pour la fonction biologique de flp1pwt et sa localisation correcte pendant la mitose. Par contre, le domaine C de flp1pwt confère une localisation nucléolaire adéquate en G2/interphase. Mes données suggèrent que la déphosphorylation de substrats par flp1pwt est dispensable pour sa localisation correcte excepté celle à l'anneau médian, qui requiert dans ce cas, l'activité catalytique de flp1pwt, révélant ainsi un niveau de régulation supplémentaire. Toutes les fonctions de flp1 pwt testées jusqu'à présent nécessitent également son activité catalytique, ce qui accentue l'importance de l'identification future de ses substrats. Comme cela a déjà été décrit pour d'autres orthologues, la capacité d'auto-intéraction et le niveau de phosphorylation pourraient contrôler l'activité de flp1pwt. En effet, mes données suggèrent que flp1pwt forme des oligomères in vivo et que la phosphorylation n'est pas essentielle pour les changements de localisation observés pour la protéine. De plus, la forme hypophosphorylée de flp1pwt pourrait être spécifiquement impliquée dans la promotion de la cytokinèse. De multiples modes de régulation incluant la localisation, l'auto-association et la phosphorylation semblent permettre un contrôle fin et subtil de l'activité de la phosphatase flp1pwt, et plus généralement celle des protéines de la famille de Cdc14pwt.

Machine learning analysis and modeling of interest rate curves

The present research deals with the review of the analysis and modeling of Swiss franc interest rate curves (IRC) by using unsupervised (SOM, Gaussian Mixtures) and supervised machine (MLP) learning algorithms. IRC are considered as objects embedded into different feature spaces: maturities; maturity-date, parameters of Nelson-Siegel model (NSM). Analysis of NSM parameters and their temporal and clustering structures helps to understand the relevance of model and its potential use for the forecasting. Mapping of IRC in a maturity-date feature space is presented and analyzed for the visualization and forecasting purposes.

How Much Does It Cost? Optimization of Costs in Sequence Analysis of Social Science Data

One major methodological problem in analysis of sequence data is the determination of costs from which distances between sequences are derived. Although this problem is currently not optimally dealt with in the social sciences, it has some similarity with problems that have been solved in bioinformatics for three decades. In this article, the authors propose an optimization of substitution and deletion/insertion costs based on computational methods. The authors provide an empirical way of determining costs for cases, frequent in the social sciences, in which theory does not clearly promote one cost scheme over another. Using three distinct data sets, the authors tested the distances and cluster solutions produced by the new cost scheme in comparison with solutions based on cost schemes associated with other research strategies. The proposed method performs well compared with other cost-setting strategies, while it alleviates the justification problem of cost schemes.

Spatiotemporal analysis of multichannel EEG: CARTOOL.

This paper describes methods to analyze the brain's electric fields recorded with multichannel Electroencephalogram (EEG) and demonstrates their implementation in the software CARTOOL. It focuses on the analysis of the spatial properties of these fields and on quantitative assessment of changes of field topographies across time, experimental conditions, or populations. Topographic analyses are advantageous because they are reference independents and thus render statistically unambiguous results. Neurophysiologically, differences in topography directly indicate changes in the configuration of the active neuronal sources in the brain. We describe global measures of field strength and field similarities, temporal segmentation based on topographic variations, topographic analysis in the frequency domain, topographic statistical analysis, and source imaging based on distributed inverse solutions. All analysis methods are implemented in a freely available academic software package called CARTOOL. Besides providing these analysis tools, CARTOOL is particularly designed to visualize the data and the analysis results using 3-dimensional display routines that allow rapid manipulation and animation of 3D images. CARTOOL therefore is a helpful tool for researchers as well as for clinicians to interpret multichannel EEG and evoked potentials in a global, comprehensive, and unambiguous way.

Magnetic resonance imaging correlates of first-episode psychosis in young adult male patients: combined analysis of grey and white matter.

Background: Several patterns of grey and white matter changes have been separately described in young adults with first-episode psychosis. Concomitant investigation of grey and white matter densities in patients with first-episode psychosis without other psychiatric comorbidities that include all relevant imaging markers could provide clues to the neurodevelopmental hypothesis in schizophrenia. Methods: We recruited patients with first-episode psychosis diagnosed according to the DSM-IV-TR and matched controls. All participants underwent magnetic resonance imaging (MRI). Voxel-based morphometry (VBM) analysis and mean diffusivity voxel-based analysis (VBA) were used for grey matter data. Fractional anisotropy and axial, radial and mean diffusivity were analyzed using tract-based spatial statistics (TBSS) for white matter data. Results: We included 15 patients and 16 controls. The mean diffusivity VBA showed significantly greater mean diffusivity in the first-episode psychosis than in the control group in the lingual gyrus bilaterally, the occipital fusiform gyrus bilaterally, the right lateral occipital gyrus and the right inferior temporal gyrus. Moreover, the TBSS analysis revealed a lower fractional anisotropy in the first-episode psychosis than in the control group in the genu of the corpus callosum, minor forceps, corticospinal tract, right superior longitudinal fasciculus, left middle cerebellar peduncle, left inferior longitudinal fasciculus and the posterior part of the fronto-occipital fasciculus. This analysis also revealed greater radial diffusivity in the first-episode psychosis than in the control group in the right corticospinal tract, right superior longitudinal fasciculus and left middle cerebellar peduncle. Limitations: The modest sample size and the absence of women in our series could limit the impact of our results. Conclusion: Our results highlight the structural vulnerability of grey matter in posterior areas of the brain among young adult male patients with first-episode psychosis. Moreover, the concomitant greater radial diffusivity within several regions already revealed by the fractional anisotropy analysis supports the idea of a late myelination in patients with first-episode psychosis.

Individual prediction of cognitive decline in mild cognitive impairment using support vector machine-based analysis of diffusion tensor imaging data.

Although cross-sectional diffusion tensor imaging (DTI) studies revealed significant white matter changes in mild cognitive impairment (MCI), the utility of this technique in predicting further cognitive decline is debated. Thirty-five healthy controls (HC) and 67 MCI subjects with DTI baseline data were neuropsychologically assessed at one year. Among them, there were 40 stable (sMCI; 9 single domain amnestic, 7 single domain frontal, 24 multiple domain) and 27 were progressive (pMCI; 7 single domain amnestic, 4 single domain frontal, 16 multiple domain). Fractional anisotropy (FA) and longitudinal, radial, and mean diffusivity were measured using Tract-Based Spatial Statistics. Statistics included group comparisons and individual classification of MCI cases using support vector machines (SVM). FA was significantly higher in HC compared to MCI in a distributed network including the ventral part of the corpus callosum, right temporal and frontal pathways. There were no significant group-level differences between sMCI versus pMCI or between MCI subtypes after correction for multiple comparisons. However, SVM analysis allowed for an individual classification with accuracies up to 91.4% (HC versus MCI) and 98.4% (sMCI versus pMCI). When considering the MCI subgroups separately, the minimum SVM classification accuracy for stable versus progressive cognitive decline was 97.5% in the multiple domain MCI group. SVM analysis of DTI data provided highly accurate individual classification of stable versus progressive MCI regardless of MCI subtype, indicating that this method may become an easily applicable tool for early individual detection of MCI subjects evolving to dementia.

Nonlinear analysis of human physical activity patterns in health and disease.

The reliable and objective assessment of chronic disease state has been and still is a very significant challenge in clinical medicine. An essential feature of human behavior related to the health status, the functional capacity, and the quality of life is the physical activity during daily life. A common way to assess physical activity is to measure the quantity of body movement. Since human activity is controlled by various factors both extrinsic and intrinsic to the body, quantitative parameters only provide a partial assessment and do not allow for a clear distinction between normal and abnormal activity. In this paper, we propose a methodology for the analysis of human activity pattern based on the definition of different physical activity time series with the appropriate analysis methods. The temporal pattern of postures, movements, and transitions between postures was quantified using fractal analysis and symbolic dynamics statistics. The derived nonlinear metrics were able to discriminate patterns of daily activity generated from healthy and chronic pain states.