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.

Structure and energetics of mixed 4He-3He drops

Using a finite-range density functional, we have investigated the energetics and structural features of mixed helium clusters. The possibility of doping the cluster with a molecule of sulfur hexafluoride is also considered. It is seen that the repulsion introduced by the impurity strongly modifies the properties of the smallest drops. Although only a qualitative comparison is possible, the gross features displayed by our calculations are in agreement with recent experimental findings.

From nonwetting to prewetting: The asymptotic behavior of 4He drops on alkali substrates

We investigate the spreading of 4He droplets on alkali-metal surfaces at zero temperature, within the frame of finite range density-functional theory. The equilibrium configurations of several 4HeN clusters and their asymptotic trend with increasing particle number N, which can be traced to the wetting behavior of the quantum fluid, are examined for nanoscopic droplets. We discuss the size effects inferring that the asymptotic properties of large droplets correspond to those of the prewetting film.

MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

Effect of the nanoparticles on the structure and crystallization of amorphous silicon thin films produced by rf glow discharge

Thin films of nanostructured silicon (ns-Si:H) were deposited by plasma-enhanced chemical vapor deposition in the presence of silicon nanoparticles at 100 C substrate temperature using silane and hydrogen gas mixture under continuous wave (cw) plasma conditions. The nanostructure of the films has been demonstrated by diverse ways: transmission electron microscopy, Raman spectroscopy and x-ray diffraction, which have shown the presence of ordered silicon clusters (1!=2 nm) embedded in an amorphous silicon matrix. Due to the presence of these ordered domains, the films crystallize faster than standard hydrogenated amorphous silicon samples, as evidenced by electrical measurements during the thermal annealing.

Microwave absorption and magnetization tunneling in Mn12-acetate molecular clusters

In this paper we study the effect of microwave absorption on the quantum relaxation rate of Mn12 molecular clusters. We have determined first the resonant frequencies of a microwave resonator containing a single crystal of Mn12-Acetate and measured initial isothermal magnetization curves while microwave power was put into the resonator. We have found that the tunneling rate changes one order of magnitude for certain frequencies. This suggests that the microwave shaking of the nuclear spin and molecular vibrational degrees of freedom is responsible for the huge increasing of the tunneling rate.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Synovial sarcomas usually metastasize after >5 years: a multicenter retrospective analysis with minimum follow-up of 10 years for survivors.

Background: Synovial sarcoma (SS) is a malignant soft tissue sarcoma with a poor prognosis because of late local recurrence and distant metastases. To our knowledge, no studies have minimum follow-up of 10 years that evaluate long-term outcomes for survivors. Patients and methods: Data on 62 patients who had been treated for SS from 1968 to 1999 were studied retrospectively in a multicenter study. Mean follow-up of living patients was 17.2 years and of dead patients 7.7 years. Results: Mean age at diagnosis was 35.4 years (range 6-82 years). Overall survival was 38.7%. The 5-year survival was 74.2%; 10-year survival was 61.2%; and 15-year survival was 46.5%. Fifteen patients (24%) died of disease after 10 years of follow-up. Local recurrence occurred after a mean of 3.6 years (range 0.5-14.9 years) and metastases at a mean of 5.7 years (range 0.5-16.3 years). Only four patients were treated technically correctly with a planned biopsy followed by a wide resection or amputation. Factors associated with significantly worse prognosis included larger tumor size, metastases at the time of diagnosis, high-grade histology, trunk-related disease, and lack of wide resection as primary surgical treatment. Conclusions: In SS, metastases develop late with high mortality. Patients with SS should be followed for >10 years.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Basement lithostratigraphy of the Adula nappe: implications for Palaeozoic evolution and Alpine kinematics

The Adula nappe belongs to the Lower Penni- nic domain of the Central Swiss Alps. It consists mostly of pre-Triassic basement lithologies occurring as strongly folded and sheared gneisses of various types with mafic boudins. We propose a new lithostratigraphy for the northern Adula nappe basement that is supported by detailed field investigations, U-Pb zircon geochronology, and whole-rock geochemistry. The following units have been identified: Cambrian clastic metasediments with abundant carbonate lenses and minor bimodal magmatism (Salahorn Formation); Ordovician metapelites associated with amphibolite boudins with abundant eclogite relicts representing oceanic metabasalts (Trescolmen Formation); Ordovician peraluminous metagranites of calc-alkaline affinity ascribed to subduction-related magmatism (Ga- renstock Augengneiss); Ordovician metamorphic volcanic- sedimentary deposits (Heinisch Stafel Formation); Early Permian post-collisional granites recording only Alpine orogenic events (Zervreila orthogneiss). All basement lithologies except the Permian granites record a Vari- scan ? Alpine polyorogenic metamorphic history. They document a complex Paleozoic geotectonic evolution consistent with the broader picture given by the pre- Mesozoic basement framework in the Alps. The internal consistency of the Adula basement lithologies and the stratigraphic coherence of the overlying Triassic sediments suggest that most tectonic contacts within the Adula nappe are pre-Alpine in age. Consequently, me ́lange models for the Tertiary emplacement of the Adula nappe are not consistent and must be rejected. The present-day structural complexity of the Adula nappe is the result of the intense Alpine ductile deformation of a pre-structured entity.

Debating the centre : the theory of principal qualificand (mukhyavisesya) in classical Indian philosophy

The theory of language has occupied a special place in the history of Indian thought. Indian philosophers give particular attention to the analysis of the cognition obtained from language, known under the generic name of śābdabodha. This term is used to denote, among other things, the cognition episode of the hearer, the content of which is described in the form of a paraphrase of a sentence represented as a hierarchical structure. Philosophers submit the meaning of the component items of a sentence and their relationship to a thorough examination, and represent the content of the resulting cognition as a paraphrase centred on a meaning element, that is taken as principal qualificand (mukhyaviśesya) which is qualified by the other meaning elements. This analysis is the object of continuous debate over a period of more than a thousand years between the philosophers of the schools of Mimāmsā, Nyāya (mainly in its Navya form) and Vyākarana. While these philosophers are in complete agreement on the idea that the cognition of sentence meaning has a hierarchical structure and share the concept of a single principal qualificand (qualified by other meaning elements), they strongly disagree on the question which meaning element has this role and by which morphological item it is expressed. This disagreement is the central point of their debate and gives rise to competing versions of this theory. The Mïmāmsakas argue that the principal qualificand is what they call bhāvanā ̒bringing into being̒, ̒efficient force̒ or ̒productive operation̒, expressed by the verbal affix, and distinct from the specific procedures signified by the verbal root; the Naiyāyikas generally take it to be the meaning of the word with the first case ending, while the Vaiyākaranas take it to be the operation expressed by the verbal root. All the participants rely on the Pāninian grammar, insofar as the Mimāmsakas and Naiyāyikas do not compose a new grammar of Sanskrit, but use different interpretive strategies in order to justify their views, that are often in overt contradiction with the interpretation of the Pāninian rules accepted by the Vaiyākaranas. In each of the three positions, weakness in one area is compensated by strength in another, and the cumulative force of the total argumentation shows that no position can be declared as correct or overall superior to the others. This book is an attempt to understand this debate, and to show that, to make full sense of the irreconcilable positions of the three schools, one must go beyond linguistic factors and consider the very beginnings of each school's concern with the issue under scrutiny. The texts, and particularly the late texts of each school present very complex versions of the theory, yet the key to understanding why these positions remain irreconcilable seems to lie elsewhere, this in spite of extensive argumentation involving a great deal of linguistic and logical technicalities. Historically, this theory arises in Mimāmsā (with Sabara and Kumārila), then in Nyāya (with Udayana), in a doctrinal and theological context, as a byproduct of the debate over Vedic authority. The Navya-Vaiyākaranas enter this debate last (with Bhattoji Dïksita and Kaunda Bhatta), with the declared aim of refuting the arguments of the Mïmāmsakas and Naiyāyikas by bringing to light the shortcomings in their understanding of Pāninian grammar. The central argument has focused on the capacity of the initial contexts, with the network of issues to which the principal qualificand theory is connected, to render intelligible the presuppositions and aims behind the complex linguistic justification of the classical and late stages of this debate. Reading the debate in this light not only reveals the rationality and internal coherence of each position beyond the linguistic arguments, but makes it possible to understand why the thinkers of the three schools have continued to hold on to three mutually exclusive positions. They are defending not only their version of the principal qualificand theory, but (though not openly acknowledged) the entire network of arguments, linguistic and/or extra-linguistic, to which this theory is connected, as well as the presuppositions and aims underlying these arguments.

Investigation of classical epidemiological links between patients harbouring identical, non-predominant meticillin-resistant Staphylococcus aureus genotypes and lessons for epidemiological tracking.

According to molecular epidemiology theory, two isolates belong to the same chain of transmission if they are similar according to a highly discriminatory molecular typing method. This has been demonstrated in outbreaks, but is rarely studied in endemic situations. Person-to-person transmission cannot be established when isolates of meticillin-resistant Staphylococcus aureus (MRSA) belong to endemically predominant genotypes. By contrast, isolates of infrequent genotypes might be more suitable for epidemiological tracking. The objective of the present study was to determine, in newly identified patients harbouring non-predominant MRSA genotypes, whether putative epidemiological links inferred from molecular typing could replace classical epidemiology in the context of a regional surveillance programme. MRSA genotypes were defined using double-locus sequence typing (DLST) combining clfB and spa genes. A total of 1,268 non-repetitive MRSA isolates recovered between 2005 and 2006 in Western Switzerland were typed: 897 isolates (71%) belonged to four predominant genotypes, 231 (18%) to 55 non-predominant genotypes, and 140 (11%) were unique. Obvious epidemiological links were found in only 106/231 (46%) patients carrying isolates with non-predominant genotypes suggesting that molecular surveillance identified twice as many clusters as those that may have been suspected with classical epidemiological links. However, not all of these molecular clusters represented person-to-person transmission. Thus, molecular typing cannot replace classical epidemiology but is complementary. A prospective surveillance of MRSA genotypes could help to target epidemiological tracking in order to recognise new risk factors in hospital and community settings, or emergence of new epidemic clones.

Color of berry and juice of 'Isabel' grape treated with abscisic acid in different ripening stages

Abstract : The objective of this work was to evaluate the effect of (S)-cis-abscisic acid (S-ABA) application at different ripening stages, in increasing phenolic compounds and color of berry and juice of 'Isabel' grape (Vitis labrusca). The evaluated treatments were: control, without S-ABA application; 400 mg L-1S-ABA applied 7 days before veraison (DBV) + 400 mg L-1S-ABA at 35 days after first application (DAFA); 400 mg L-1S-ABA applied at veraison (V) + 400 mg L-1S-ABA at 35 DAFA; and 400 mg L-1S-ABA applied 7 days after veraison (DAV) + 400 mg L-1S-ABA at 35 DAFA. There was no difference among treatments regarding the physical characteristics of berries and clusters, as well as total polyphenols in berry and juice. However, there was an increase in total anthocyanins in berry and juice with S-ABA application. Colorimetric variables indicated the increase in color of berry treated with S-ABA. Juices produced from grapes treated with S-ABA were more appreciated by tasters. The treatments with 400 mg L-1S-ABA applied 7 days before, during, or 7 days after veraison, combined with an additional application 35 days after the first one, increment total anthocyanin concentration and color of berry and juice of 'Isabel' grape, with better juice acceptance, without affecting total polyphenol concentration.

Identification of in vitro HSC fate regulators by differential lipid raft clustering.

Most hematopoietic stem cells (HSC) in the bone marrow reside in a quiescent state and occasionally enter the cell cycle upon cytokine-induced activation. Although the mechanisms regulating HSC quiescence and activation remain poorly defined, recent studies have revealed a role of lipid raft clustering (LRC) in HSC activation. Here, we tested the hypothesis that changes in lipid raft distribution could serve as an indicator of the quiescent and activated state of HSCs in response to putative niche signals. A semi-automated image analysis tool was developed to map the presence or absence of lipid raft clusters in live HSCs cultured for just one hour in serum-free medium supplemented with stem cell factor (SCF). By screening the ability of 19 protein candidates to alter lipid raft dynamics, we identified six factors that induced either a marked decrease (Wnt5a, Wnt3a and Osteopontin) or increase (IL3, IL6 and VEGF) in LRC. Cell cycle kinetics of single HSCs exposed to these factors revealed a correlation of LRC dynamics and proliferation kinetics: factors that decreased LRC slowed down cell cycle kinetics, while factors that increased LRC led to faster and more synchronous cycling. The possibility of identifying, by LRC analysis at very early time points, whether a stem cell is activated and possibly committed upon exposure to a signaling cue of interest could open up new avenues for large-scale screening efforts.

Grx5 is a mitochondrial glutaredoxin required for the activity of iron/sulfur enzymes

Yeast cells contain a family of three monothiol glutaredoxins: Grx3, 4, and 5. Absence of Grx5 leads to constitutive oxidative damage, exacerbating that caused by external oxidants. Phenotypic defects associated with the absence of Grx5 are suppressed by overexpression ofSSQ1 and ISA2, two genes involved in the synthesis and assembly of iron/sulfur clusters into proteins. Grx5 localizes at the mitochondrial matrix, like other proteins involved in the synthesis of these clusters, and the mature form lacks the first 29 amino acids of the translation product. Absence of Grx5 causes: 1) iron accumulation in the cell, which in turn could promote oxidative damage, and 2) inactivation of enzymes requiring iron/sulfur clusters for their activity. Reduction of iron levels in grx5 null mutants does not restore the activity of iron/sulfur enzymes, and cell growth defects are not suppressed in anaerobiosis or in the presence of disulfide reductants. Hence, Grx5 forms part of the mitochondrial machinery involved in the synthesis and assembly of iron/sulfur centers.