Clustering and hot spot detection in socio-economic spatio-temporal data

Distribution of socio-economic features in urban space is an important source of information for land and transportation planning. The metropolization phenomenon has changed the distribution of types of professions in space and has given birth to different spatial patterns that the urban planner must know in order to plan a sustainable city. Such distributions can be discovered by statistical and learning algorithms through different methods. In this paper, an unsupervised classification method and a cluster detection method are discussed and applied to analyze the socio-economic structure of Switzerland. The unsupervised classification method, based on Ward's classification and self-organized maps, is used to classify the municipalities of the country and allows to reduce a highly-dimensional input information to interpret the socio-economic landscape. The cluster detection method, the spatial scan statistics, is used in a more specific manner in order to detect hot spots of certain types of service activities. The method is applied to the distribution services in the agglomeration of Lausanne. Results show the emergence of new centralities and can be analyzed in both transportation and social terms.

Decoupled post-glacial history in mutualistic plant-insect interactions: insights from the yellow loosestrife (Lysimachia vulgaris) and its associated oil-collecting bees (Macropis europaea and M-fulvipes)

AimWe take a comparative phylogeographical approach to assess whether three species involved in a specialized oil-rewarding pollination system (i.e. Lysimachia vulgaris and two oil-collecting bees within the genus Macropis) show congruent phylogeographical trajectories during post-glacial colonization processes. Our working hypothesis is that within specialized mutualistic interactions, where each species relies on the co-occurrence of the other for survival and/or reproduction, partners are expected to show congruent evolutionary trajectories, because they are likely to have followed parallel migration routes and to have shared glacial refugia. LocationWestern Palaearctic. MethodsOur analysis relies on the extensive sampling of 104 Western Palaearctic populations (totalling 434, 159 and 74 specimens of Lysimachiavulgaris, Macropiseuropaea and Macropisfulvipes, respectively), genotyped with amplified fragment length polymorphism. Based on this, we evaluated the regional genetic diversity (Shannon diversity and allele rarity index) and genetic structure (assessed using structure, population networks, isolation-by-distance and spatial autocorrelation metrics) of each species. Finally, we compared the general phylogeographical patterns obtained. ResultsContrary to our expectations, the analyses revealed phylogeographical signals suggesting that the investigated organisms demonstrate independent post-glacial trajectories as well as distinct contemporaneous demographic parameters, despite their mutualistic interaction. Main conclusionsThe mutualistic partners investigated here are likely to be experiencing distinct and independent evolutionary dynamics because of their contrasting life-history traits (e.g. dispersal abilities), as well as distinct hubs and migration routes. Such conditions would prevent and/or erase any signature of co-structuring of lineages in space and time. As a result, the lack of phylogeographical congruence driven by differences in life-history traits might have arisen irrespective of the three species having shared similar Pleistocene glacial refugia.

Connexins in lymphatic vessel physiology and disease.

Connexins are transmembrane proteins that form gap junction- and hemi-channels. Once inserted into the membrane, hemi-channels (connexons) allow for diffusion of ions and small molecules (<1kDa) between the extracellular space and the cytosol. Gap junction channels allow diffusion of similar molecules between the cytoplasms of adjacent cells. The expression and function of connexins in blood vessels has been intensely studied in the last few decades. In contrast, only a few studies paid attention to lymphatic vessels; convincing in vivo data with respect to expression patterns of lymphatic connexins and their functional roles have only recently begun to emerge. Interestingly, mutations in connexin genes have been linked to diseases of lymphatic vasculature, most notably primary and secondary lymphedema. This review summarizes the available data regarding lymphatic connexins. More specifically it addresses (i) early studies aimed at presence of gap junction-like structures in lymphatic vessels, (ii) more recent studies focusing on lymphatic connexins using genetically engineered mice, and (iii) results of clinical studies that have reported lymphedema-linked mutations in connexin genes.

Using environmental niche modeling to understand biological invasions in a changing world

L'activité humaine affecte particulièrement la biodiversité, qui décline à une vitesse préoccupante. Parmi les facteurs réduisant la biodiversité, on trouve les espèces envahissantes. Symptomatiques d'un monde globalisé où l'échange se fait à l'échelle de la planète, certaines espèces, animales ou végétales, sont introduites, volontairement ou accidentellement par l'activité humaine (par exemple lors des échanges commerciaux ou par les voyageurs). Ainsi, ces espèces atteignent des régions qu'elles n'auraient jamais pu coloniser naturellement. Une fois introduites, l'absence de compétiteur peut les rendre particulièrement nuisibles. Ces nuisances sont plus ou moins directes, allant de problèmes sanitaires (p. ex. les piqûres très aigües des fourmis de feu, originaires d'Amérique du Sud et colonisant à une vitesse fulgurante les USA, l'Australie ou la Chine) à des nuisances sur la biodiversité (p. ex. les ravages de la perche du Nil sur la diversité unique des poissons Cichlidés du Lac Victoria). Il est donc important de pouvoir prévenir de telles introductions. De plus, pour le biologiste, ces espèces représentent une rare occasion de pouvoir comprendre les mécanismes évolutifs et écologiques qui expliquent le succès des envahissantes dans un monde où les équilibres sont bouleversés. Les modèles de niche environnementale sont un outil particulièrement utile dans le cadre de cette problématique. En reliant des observations d'espèces aux conditions environnementales où elles se trouvent, ils peuvent prédire la distribution potentielle des envahissantes, permettant d'anticiper et de mieux limiter leur impact. Toutefois, ils reposent sur des hypothèses pas évidentes à démontrer. L'une d'entre elle étant que la niche d'une espèce reste constante dans le temps, et dans l'espace. Le premier objectif de mon travail est de comparer si la niche d'une espèce envahissante diffère entre sa distribution d'origine native et celle d'origine introduite. En étudiant 50 espèces de plantes et 168 espèces de Mammifères, je démontre que c'est le cas et que par corolaire, il est possible de prédire leurs distributions. La deuxième partie de mon travail consiste à comprendre quelles seront les interactions entre le changement climatiques et les envahissantes, afin d'estimer leur impact sous un climat réchauffé. En étudiant la distribution de 49 espèces de plantes envahissantes, je démontre que les montagnes, régions relativement préservée par ce problème, deviendront bien plus exposées aux risques d'invasions biologiques. J'expose aussi comment les interactions entre l'activité humaine, le réchauffement climatique et les espèces envahissantes menacent la vigne sauvage en Europe et propose des zones géographiques particulièrement adaptée pour sa conservation. Enfin, à une échelle beaucoup plus locale, je montre qu'il est possible d'utiliser ces modèles de niches le long d'une rivière à une échelle extrêmement fine (1 mètre), potentiellement utile pour rationnaliser des mesures de conservations sur le terrain. - Biodiversity is significantly negatively affected by human activity. Invasive species are one of the most important factors causing biodiversity's decline. Intimately linked to the era of global trade, some plant or animal species can be accidentally or casually introduced with human activity (e.g. trade or travel). In this way, these species reach areas they could never reach through natural dispersal. Once naturalized, the lack of competitors can make these species highly noxious. Their effect is more or less direct, from sanitary problems (e.g. the harmful sting of Fire Ants, originating from South America and now spreading throughout USA, China and Australia) or can affect biodiversity (e.g. the Nile perch, devastating the one of the richest hotspot of Cichlid fishes diversity in Lake Victoria). It is thus important to prevent such harmful introductions. Moreover, invasive species represent for biologists one of the rare occasions to understand the evolutionary and ecological mechanisms behind the success of invaders in a world where natural equilibrium is already disturbed. Environmental niche models are particularly useful to tackle this problematic. By relating species observation to the environmental conditions where they occur, they can predict the potential distribution of invasive species, allowing a better anticipation and thus limiting their impact. However, they rely on strong assumption, one of the most important being that the modeled niche remains constant through space and time. The first aim of my thesis is to quantify the difference between the native and the invaded niche. By investigating 50 plant and 168 mammal species, I show that the niche is at least partially conserved, supporting for reliable predictions of invasive' s potential distributions. The second aim of my thesis is to understand the possible interactions between climate change and invasive species, such as to assess their impact under a warmer climate. By studying 49 invasive plant species, I show that mountain areas, which were relatively preserved, will become more suitable for biological invasions. Additionally, I show how interactions between human activity, global warming and invasive species are threatening the wild grapevine in Europe and propose geographical areas particularly adapted for conservation measures. Finally, at a much finer scale where conservation plannings ultimately take place, I show that it is possible to model the niche at very high resolution (1 meter) in an alluvial area allowing better prioritizations for conservation.

Forensic intelligence framework. Part II: study of the main generic building blocks and challenges through the examples of illicit drugs and false identity documents monitoring

The development of forensic intelligence relies on the expression of suitable models that better represent the contribution of forensic intelligence in relation to the criminal justice system, policing and security. Such models assist in comparing and evaluating methods and new technologies, provide transparency and foster the development of new applications. Interestingly, strong similarities between two separate projects focusing on specific forensic science areas were recently observed. These observations have led to the induction of a general model (Part I) that could guide the use of any forensic science case data in an intelligence perspective. The present article builds upon this general approach by focusing on decisional and organisational issues. The article investigates the comparison process and evaluation system that lay at the heart of the forensic intelligence framework, advocating scientific decision criteria and a structured but flexible and dynamic architecture. These building blocks are crucial and clearly lay within the expertise of forensic scientists. However, it is only part of the problem. Forensic intelligence includes other blocks with their respective interactions, decision points and tensions (e.g. regarding how to guide detection and how to integrate forensic information with other information). Formalising these blocks identifies many questions and potential answers. Addressing these questions is essential for the progress of the discipline. Such a process requires clarifying the role and place of the forensic scientist within the whole process and their relationship to other stakeholders.

Spatio-temporal differentiation and sociality in spiders.

Species that differ in their social system, and thus in traits such as group size and dispersal timing, may differ in their use of resources along spatial, temporal, or dietary dimensions. The role of sociality in creating differences in habitat use is best explored by studying closely related species or socially polymorphic species that differ in their social system, but share a common environment. Here we investigate whether five sympatric Anelosimus spider species that range from nearly solitary to highly social differ in their use of space and in their phenology as a function of their social system. By studying these species in Serra do Japi, Brazil, we find that the more social species, which form larger, longer-lived colonies, tend to live inside the forest, where sturdier, longer lasting vegetation is likely to offer better support for their nests. The less social species, which form single-family groups, in contrast, tend to occur on the forest edge where the vegetation is less robust. Within these two microhabitats, species with longer-lived colonies tend to occupy the potentially more stable positions closer to the core of the plants, while those with smaller and shorter-lived colonies build their nests towards the branch tips. The species further separate in their use of common habitat due to differences in the timing of their reproductive season. These patterns of habitat use suggest that the degree of sociality can enable otherwise similar species to differ from one another in ways that may facilitate their co-occurrence in a shared environment, a possibility that deserves further consideration.

The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling.

Predicting which species will occur together in the future, and where, remains one of the greatest challenges in ecology, and requires a sound understanding of how the abiotic and biotic environments interact with dispersal processes and history across scales. Biotic interactions and their dynamics influence species' relationships to climate, and this also has important implications for predicting future distributions of species. It is already well accepted that biotic interactions shape species' spatial distributions at local spatial extents, but the role of these interactions beyond local extents (e.g. 10 km(2) to global extents) are usually dismissed as unimportant. In this review we consolidate evidence for how biotic interactions shape species distributions beyond local extents and review methods for integrating biotic interactions into species distribution modelling tools. Drawing upon evidence from contemporary and palaeoecological studies of individual species ranges, functional groups, and species richness patterns, we show that biotic interactions have clearly left their mark on species distributions and realised assemblages of species across all spatial extents. We demonstrate this with examples from within and across trophic groups. A range of species distribution modelling tools is available to quantify species environmental relationships and predict species occurrence, such as: (i) integrating pairwise dependencies, (ii) using integrative predictors, and (iii) hybridising species distribution models (SDMs) with dynamic models. These methods have typically only been applied to interacting pairs of species at a single time, require a priori ecological knowledge about which species interact, and due to data paucity must assume that biotic interactions are constant in space and time. To better inform the future development of these models across spatial scales, we call for accelerated collection of spatially and temporally explicit species data. Ideally, these data should be sampled to reflect variation in the underlying environment across large spatial extents, and at fine spatial resolution. Simplified ecosystems where there are relatively few interacting species and sometimes a wealth of existing ecosystem monitoring data (e.g. arctic, alpine or island habitats) offer settings where the development of modelling tools that account for biotic interactions may be less difficult than elsewhere.

The contribution of patch topology and demographic parameters to PVA predictions: the case of the European tree frog

Population viability analyses (PVA) are increasingly used in metapopulation conservation plans. Two major types of models are commonly used to assess vulnerability and to rank management options: population-based stochastic simulation models (PSM such as RAMAS or VORTEX) and stochastic patch occupancy models (SPOM). While the first set of models relies on explicit intrapatch dynamics and interpatch dispersal to predict population levels in space and time, the latter is based on spatially explicit metapopulation theory where the probability of patch occupation is predicted given the patch area and isolation (patch topology). We applied both approaches to a European tree frog (Hyla arborea) metapopulation in western Switzerland in order to evaluate the concordances of both models and their applications to conservation. Although some quantitative discrepancies appeared in terms of network occupancy and equilibrium population size, the two approaches were largely concordant regarding the ranking of patch values and sensitivities to parameters, which is encouraging given the differences in the underlying paradigms and input data.

Contrasting magma types and timing of intrusion in the Permian layered mafic complex of Mont Collon (Western Alps, Valais, Switzerland): evidence from U/Pb zircon and 40Ar/39Ar amphibole dating

We have selected and dated three contrasting rock-types representative of the magmatic activity within the Permian layered mafic complex of Mont Collon, Austroalpine Dent Blanche nappe, Western Alps. A pegmatitic gabbro associated to the main cumulus sequence yields a concordant U/Pb zircon age of 284.2 +/- 0.6 Ma, whereas a pegmatitic granite dike crosscutting the latter yields a concordant age of 282.9 +/- 0.6 Ma. A Fe-Ti-rich ultrabasic lamprophyre, crosscutting all other lithologies of the complex, yields an 40Ar/39Ar plateau age of 260.2 +/- 0.7 Ma on a kaersutite concentrate. All ages are interpreted as magmatic. Sub-contemporaneous felsic dikes within the Mont Collon complex are ascribed to anatectic back-veining from the country-rock, related to the emplacement of the main gabbroic body in the continental crust, which is in accordance with new isotopic data. The lamprophyres have isotopic compositions typical of a depleted mantle, in contrast to those of the cumulate gabbros, close to values of the Bulk Silicate Earth. This indicates either contrasting sources for the two magma pulses - the subcontinental lithospheric mantle for the gabbros and the underlying asthenosphere for the lamprophyres - or a single depleted lithospheric source with variable degrees of crustal contamination of the gabbroic melts during their emplacement in the continental crust. The Mont Collon complex belongs to a series of Early Permian mafic massifs, which emplaced in a short time span about 285-280 Ma ago, in a limited sector of the post-Variscan continental crust now corresponding to the Austroalpine/ Southern Alpine domains and Corsica. This magmatic activity was controlled in space and time by crustal-scale transtensional shear zones.

Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

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 to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

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.

Cheminformatic studies of molecular properties and their influence on bioactivity in the context of drug discovery

Molecular shape has long been known to be an important property for the process of molecular recognition. Previous studies postulated the existence of a drug-like shape space that could be used to artificially bias the composition of screening libraries, with the aim to increase the chance of success in Hit Identification. In this work, it was analysed to which extend this assumption holds true. Normalized Principal Moments of Inertia Ratios (NPRs) have been used to describe the molecular shape of small molecules. It was investigated, whether active molecules of diverse targets are located in preferred subspaces of the NPR shape space. Results illustrated a significantly stronger clustering than could be expected by chance, with parts of the space unlikely to be occupied by active compounds. Furthermore, a strong enrichment of elongated, rather flat shapes could be observed, while globular compounds were highly underrepresented. This was confirmed for a wide range of small molecule datasets from different origins. Active compounds exhibited a high overlap in their shape distributions across different targets, making a purely shape­ based discrimination very difficult. An additional perspective was provided by comparing the shapes of protein binding pockets with those of their respective ligands. Although more globular than their ligands, it was observed that binding sites shapes exhibited a similarly skewed distribution in shape space: spherical shapes were highly underrepresented. This was different for unoccupied binding pockets of smaller size. These were on the contrary identified to possess a more globular shape. The relation between shape complementarity and exhibited bioactivity was analysed; a moderate correlation between bioactivity and parameters including pocket coverage, distance in shape space, and others could be identified, which reflects the importance of shape complementarity. However, this also suggests that other aspects are of relevance for molecular recognition. A subsequent analysis assessed if and how shape and volume information retrieved from pocket or respective reference ligands could be used as a pre-filter in a virtual screening approach. ln Lead Optimization compounds need to get optimized with respect to a variety of pararneters. Here, the availability of past success stories is very valuable, as they can guide medicinal chemists during their analogue synthesis plans. However, although of tremendous interest for the public domain, so far only large corporations had the ability to mine historical knowledge in their proprietary databases. With the aim to provide such information, the SwissBioisostere database was developed and released during this thesis. This database contains information on 21,293,355 performed substructural exchanges, corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1,948 molecular targets representing 30 target classes, and on their impact on bioactivity . A user-friendly interface was developed that provides facile access to these data and is accessible at http//www.swissbioisostere.ch. The ChEMBL database was used as primary data source of bioactivity information. Matched molecular pairs have been identified in the extracted and cleaned data. Success-based scores were developed and integrated into the database to allow re-ranking of proposed replacements by their past outcomes. It was analysed to which degree these scores correlate with chemical similarity of the underlying fragments. An unexpectedly weak relationship was detected and further investigated. Use cases of this database were envisioned, and functionalities implemented accordingly: replacement outcomes are aggregatable at the assay level, and it was shawn that an aggregation at the target or target class level could also be performed, but should be accompanied by a careful case-by-case assessment. It was furthermore observed that replacement success depends on the activity of the starting compound A within a matched molecular pair A-B. With increasing potency the probability to lose bioactivity through any substructural exchange was significantly higher than in low affine binders. A potential existence of a publication bias could be refuted. Furthermore, often performed medicinal chemistry strategies for structure-activity-relationship exploration were analysed using the acquired data. Finally, data originating from pharmaceutical companies were compared with those reported in the literature. It could be seen that industrial medicinal chemistry can access replacement information not available in the public domain. In contrast, a large amount of often-performed replacements within companies could also be identified in literature data. Preferences for particular replacements differed between these two sources. The value of combining different endpoints in an evaluation of molecular replacements was investigated. The performed studies highlighted furthermore that there seem to exist no universal substructural replacement that always retains bioactivity irrespective of the biological environment. A generalization of bioisosteric replacements seems therefore not possible. - La forme tridimensionnelle des molécules a depuis longtemps été reconnue comme une propriété importante pour le processus de reconnaissance moléculaire. Des études antérieures ont postulé que les médicaments occupent préférentiellement un sous-ensemble de l'espace des formes des molécules. Ce sous-ensemble pourrait être utilisé pour biaiser la composition de chimiothèques à cribler, dans le but d'augmenter les chances d'identifier des Hits. L'analyse et la validation de cette assertion fait l'objet de cette première partie. Les Ratios de Moments Principaux d'Inertie Normalisés (RPN) ont été utilisés pour décrire la forme tridimensionnelle de petites molécules de type médicament. Il a été étudié si les molécules actives sur des cibles différentes se co-localisaient dans des sous-espaces privilégiés de l'espace des formes. Les résultats montrent des regroupements de molécules incompatibles avec une répartition aléatoire, avec certaines parties de l'espace peu susceptibles d'être occupées par des composés actifs. Par ailleurs, un fort enrichissement en formes allongées et plutôt plates a pu être observé, tandis que les composés globulaires étaient fortement sous-représentés. Cela a été confirmé pour un large ensemble de compilations de molécules d'origines différentes. Les distributions de forme des molécules actives sur des cibles différentes se recoupent largement, rendant une discrimination fondée uniquement sur la forme très difficile. Une perspective supplémentaire a été ajoutée par la comparaison des formes des ligands avec celles de leurs sites de liaison (poches) dans leurs protéines respectives. Bien que plus globulaires que leurs ligands, il a été observé que les formes des poches présentent une distribution dans l'espace des formes avec le même type d'asymétrie que celle observée pour les ligands: les formes sphériques sont fortement sous­ représentées. Un résultat différent a été obtenu pour les poches de plus petite taille et cristallisées sans ligand: elles possédaient une forme plus globulaire. La relation entre complémentarité de forme et bioactivité a été également analysée; une corrélation modérée entre bioactivité et des paramètres tels que remplissage de poche, distance dans l'espace des formes, ainsi que d'autres, a pu être identifiée. Ceci reflète l'importance de la complémentarité des formes, mais aussi l'implication d'autres facteurs. Une analyse ultérieure a évalué si et comment la forme et le volume d'une poche ou de ses ligands de référence pouvaient être utilisés comme un pré-filtre dans une approche de criblage virtuel. Durant l'optimisation d'un Lead, de nombreux paramètres doivent être optimisés simultanément. Dans ce contexte, la disponibilité d'exemples d'optimisations réussies est précieuse, car ils peuvent orienter les chimistes médicinaux dans leurs plans de synthèse par analogie. Cependant, bien que d'un extrême intérêt pour les chercheurs dans le domaine public, seules les grandes sociétés pharmaceutiques avaient jusqu'à présent la capacité d'exploiter de telles connaissances au sein de leurs bases de données internes. Dans le but de remédier à cette limitation, la base de données SwissBioisostere a été élaborée et publiée dans le domaine public au cours de cette thèse. Cette base de données contient des informations sur 21 293 355 échanges sous-structuraux observés, correspondant à 5 586 462 remplacements uniques mesurés dans 35 039 tests contre 1948 cibles représentant 30 familles, ainsi que sur leur impact sur la bioactivité. Une interface a été développée pour permettre un accès facile à ces données, accessible à http:/ /www.swissbioisostere.ch. La base de données ChEMBL a été utilisée comme source de données de bioactivité. Une version modifiée de l'algorithme de Hussain et Rea a été implémentée pour identifier les Matched Molecular Pairs (MMP) dans les données préparées au préalable. Des scores de succès ont été développés et intégrés dans la base de données pour permettre un reclassement des remplacements proposés selon leurs résultats précédemment observés. La corrélation entre ces scores et la similarité chimique des fragments correspondants a été étudiée. Une corrélation plus faible qu'attendue a été détectée et analysée. Différents cas d'utilisation de cette base de données ont été envisagés, et les fonctionnalités correspondantes implémentées: l'agrégation des résultats de remplacement est effectuée au niveau de chaque test, et il a été montré qu'elle pourrait également être effectuée au niveau de la cible ou de la classe de cible, sous réserve d'une analyse au cas par cas. Il a en outre été constaté que le succès d'un remplacement dépend de l'activité du composé A au sein d'une paire A-B. Il a été montré que la probabilité de perdre la bioactivité à la suite d'un remplacement moléculaire quelconque est plus importante au sein des molécules les plus actives que chez les molécules de plus faible activité. L'existence potentielle d'un biais lié au processus de publication par articles a pu être réfutée. En outre, les stratégies fréquentes de chimie médicinale pour l'exploration des relations structure-activité ont été analysées à l'aide des données acquises. Enfin, les données provenant des compagnies pharmaceutiques ont été comparées à celles reportées dans la littérature. Il a pu être constaté que les chimistes médicinaux dans l'industrie peuvent accéder à des remplacements qui ne sont pas disponibles dans le domaine public. Par contre, un grand nombre de remplacements fréquemment observés dans les données de l'industrie ont également pu être identifiés dans les données de la littérature. Les préférences pour certains remplacements particuliers diffèrent entre ces deux sources. L'intérêt d'évaluer les remplacements moléculaires simultanément selon plusieurs paramètres (bioactivité et stabilité métabolique par ex.) a aussi été étudié. Les études réalisées ont souligné qu'il semble n'exister aucun remplacement sous-structural universel qui conserve toujours la bioactivité quel que soit le contexte biologique. Une généralisation des remplacements bioisostériques ne semble donc pas possible.

Spatial variation in the temporal change of male and female melanic ornamentation in the barn owl.

Because the magnitude of selection can vary between sexes and in space and time, sexually antagonistic selection is difficult to demonstrate. In a Swiss population of barn owls (Tyto alba), a heritable eumelanic colour trait (size of black spots on ventral feathers) was positively selected with respect to yearling survival only in females. It remains unclear whether the absence of negative selection in males is typical in this species. To tackle this issue indirectly, we measured the size of black spots in 1733 skin specimens collected by museums from 1816 to 2001 in seven European countries and in the Middle-East. The temporal change in spot size was sex- and country-specific. In males, spots became smaller particularly in three countries (Middle-East, Italy and Switzerland). In females, the size of spots increased significantly in two countries (UK and Spain) and decreased in two others (Germany and Switzerland). Because migration and phenotypic plasticity cannot explain these results, selection is the most likely cause. The weaker temporal change in spot size in females than males may be because of the combined effect of strong genetic correlation between the sexes and stronger negative selection in males than positive selection in females. We thus suggest that in the barn owl, spot size (or genetically correlated traits) is sexually antagonistically selected and that its pattern of selection may account for the maintenance of its variation and sexual dimorphism.

Data-driven clustering : new methods and applications

Abstract : This work is concerned with the development and application of novel unsupervised learning methods, having in mind two target applications: the analysis of forensic case data and the classification of remote sensing images. First, a method based on a symbolic optimization of the inter-sample distance measure is proposed to improve the flexibility of spectral clustering algorithms, and applied to the problem of forensic case data. This distance is optimized using a loss function related to the preservation of neighborhood structure between the input space and the space of principal components, and solutions are found using genetic programming. Results are compared to a variety of state-of--the-art clustering algorithms. Subsequently, a new large-scale clustering method based on a joint optimization of feature extraction and classification is proposed and applied to various databases, including two hyperspectral remote sensing images. The algorithm makes uses of a functional model (e.g., a neural network) for clustering which is trained by stochastic gradient descent. Results indicate that such a technique can easily scale to huge databases, can avoid the so-called out-of-sample problem, and can compete with or even outperform existing clustering algorithms on both artificial data and real remote sensing images. This is verified on small databases as well as very large problems. Résumé : Ce travail de recherche porte sur le développement et l'application de méthodes d'apprentissage dites non supervisées. Les applications visées par ces méthodes sont l'analyse de données forensiques et la classification d'images hyperspectrales en télédétection. Dans un premier temps, une méthodologie de classification non supervisée fondée sur l'optimisation symbolique d'une mesure de distance inter-échantillons est proposée. Cette mesure est obtenue en optimisant une fonction de coût reliée à la préservation de la structure de voisinage d'un point entre l'espace des variables initiales et l'espace des composantes principales. Cette méthode est appliquée à l'analyse de données forensiques et comparée à un éventail de méthodes déjà existantes. En second lieu, une méthode fondée sur une optimisation conjointe des tâches de sélection de variables et de classification est implémentée dans un réseau de neurones et appliquée à diverses bases de données, dont deux images hyperspectrales. Le réseau de neurones est entraîné à l'aide d'un algorithme de gradient stochastique, ce qui rend cette technique applicable à des images de très haute résolution. Les résultats de l'application de cette dernière montrent que l'utilisation d'une telle technique permet de classifier de très grandes bases de données sans difficulté et donne des résultats avantageusement comparables aux méthodes existantes.

Endocytosis, autophagy and JNK inhibition in neuroprotection against excitotoxicity and neonatal cerebral ischemia

Abstract : Neonatal stroke occurs in 1 out of 4000 live births and usually leads to serious motor and cognitive disabilities. Ischemic brain injury results from a complex of pathophysiological events that evolve over space and time making it difficult to devise successful therapy. To date, there are no effective treatments for perinatal brain damage. Most clinical trials of neuroprotectaot drugs have failed because of their side-effects. For this reason it is important to find ways to target drugs specifically into the stressed cells. In this study we plan to contribute to the development of an efficient neuroprotective strategy against excitotoxic cell death in the neonate. In order to achieve this goal, several strategies were followed. A recently described phenomenon of induced endocytosis associated with excitotoxicity was more deeply investigated. As a simplified model we used dissociated cortical neurons exposed to an excitotoxic dose of NMDA, and we showed that this phenomenon depends on clathrin and dynamin. Using a model of neonatal focal cerebral ischemia, we demonstrated that the excitotoxicity-related endocytosis targets molecules such as TAT peptides into stressed neurons. These appear to be viable, raising the possibility of using this phenomenon as a doorway for neuroprotection. One part of the project was devoted to the study of the TAT-conjugated JNK inhibitory peptide, D-JNKI1. Adose-response study showed strong neuroprotection over a wide dose-range in the case of delayed administration (either intravenous or intraperitoneal). Since D-JNKI1 is aTAT-linked peptide, we investigated the role of its own NMDA-induced endocytosis in its neuroprotective efficacy. Furthermore, we showed that this endocytosis is JNK dependent, and that D-JNKI1 regulates its own uptake. We additionally studied the different types of cell death involved in a model of neonatal focal cerebral ischemia. Necrosis occurred rapidly in the center of the lesion whereas apoptosis and autophagic cell death occurred late at the lesion border. Inhibiting apoptosis was not protective, but use of autophagy inhibitor 3methyladenine provided a strong neuroprotection. Finally, combining two neuroprotectants that target different intracellular pathways was neuroprotective in a severe model of cerebral ischemia where neither of the drugs was efficient when administered individually. Résumé : L'ischémie néonatale connaît une incidence de 1 naissance sur 4000, entraînant généralement de sérieux dysfonctionnements moteurs et cognitifs. L'ischémie cérébrale résulte d'évènements physiopathologiques complexes qui évoluent dans l'espace et le temps rendant difficile la conception de thérapies efficaces. A l'heure actuelle, aucun traitement n'existe pour lutter contre les accidents vasculaires cérébraux qui se produisent autour de la naissance. La plupart des essais cliniques concernant des molécules neuroprotectrices ont échoué du fait de leurs effets secondaires néfastes. Pour cette raison, il est important de trouver des moyens de cibler les drogues dans les cellules stressées spécifiquement. Dans cette étude nous visons à participer au développement d'une stratégie neuroprotectrice efficace contre l'ischémie cérébrale chez le nouveau-né. Dans ce but, plusieurs stratégies ont été poursuivies. Un nouveau phénomène d'endocytose induite par un stimulus excitotoxique a été récemment décrit. Une partie de cette étude va consister à mieux comprendre ce phénomène. Pour céla, nous avons utilisé comme modèle d'étude simplifié des cultures dissociées de neurones corticaux exposées à une dose excitotoxique de NMDA. Nous avons ainsi montré que cette endocytose associée à l'excitotoxicité dépend de la clathrine et de la dynamine. A l'aide d'un modèle d'ischémie cérébrale focale chez le raton de 12 jours, nous avons démontré que cette endocytose induite par l'excitotoxicité permet de cibler des molécules diverses et en particulier les peptides TAT dans les neurones stressés. Ces neurones fortement endocytiques apparaissent comme étant encore viables, ouvrant la possibilité d'utiliser cette endocytose comme moyen d'entrée pour des molécules thérapeutiques. Une partie du projet a été consacrée à l'étude d'un inhibiteur de la voie JNK, couplé au TAT, appelé D-JNKI1. Des études de dose réponse du D-JNKI1 ont été réalisées chez l'animal, testant les effets d'une administration retardée en injection intraveineuse ou intra péritonéale. Ces études démontrent qu'une large gamme de dose permet d'obCenir une réduction de la taille de la lésion. Comme D-JNK11 est couplé au peptide TAT, nous avons étudié la contribution que sa propre endocytose lors de l'excitotoxicité apporte à ses effets protecteurs. Par ailleurs, nous avons montré que cette endocytose induite par l'excitotoxicité dépend de la voie de signalisation JNK et que D-JNK11 est donc capable de réguler sa propre entrée. Nous avons en parallèle étudié les différents types de mort cellulaires impliqués dans le développement de la lésion dans un modèle sévère d'ischémie cérébrale chez le raton nouveau-né. La mort cellulaire par nécrose se développe rapidement dans le centre de la lésion alors que les morts cellulaires par apoptose et autophagique vont apparaître plus tard et au bord de la lésion. Inhiber l'apoptose n'a pas permis de réduire la taille de la lésion alors que l'utilisation d'un inhibiteur d'autophagie, la 3-méthyladénine, procure une forte neuroprotection. Finalement, la combinaison de deux peptides qui ciblent différentes voies de signalisation intracellulaire permet d'obtenir une bonne protection dans le modèle d'ischémie sévère dans lequel aucun des deux peptides administré séparément n'a donné d'effets bénéfiques.