Multi-scale support vector algorithms for hot spot detection and modelling

The algorithmic approach to data modelling has developed rapidly these last years, in particular methods based on data mining and machine learning have been used in a growing number of applications. These methods follow a data-driven methodology, aiming at providing the best possible generalization and predictive abilities instead of concentrating on the properties of the data model. One of the most successful groups of such methods is known as Support Vector algorithms. Following the fruitful developments in applying Support Vector algorithms to spatial data, this paper introduces a new extension of the traditional support vector regression (SVR) algorithm. This extension allows for the simultaneous modelling of environmental data at several spatial scales. The joint influence of environmental processes presenting different patterns at different scales is here learned automatically from data, providing the optimum mixture of short and large-scale models. The method is adaptive to the spatial scale of the data. With this advantage, it can provide efficient means to model local anomalies that may typically arise in situations at an early phase of an environmental emergency. However, the proposed approach still requires some prior knowledge on the possible existence of such short-scale patterns. This is a possible limitation of the method for its implementation in early warning systems. The purpose of this paper is to present the multi-scale SVR model and to illustrate its use with an application to the mapping of Cs137 activity given the measurements taken in the region of Briansk following the Chernobyl accident.

Evaluating sampling strategies and logistic regression methods for modelling complex land cover changes

The role of land cover change as a significant component of global change has become increasingly recognized in recent decades. Large databases measuring land cover change, and the data which can potentially be used to explain the observed changes, are also becoming more commonly available. When developing statistical models to investigate observed changes, it is important to be aware that the chosen sampling strategy and modelling techniques can influence results. We present a comparison of three sampling strategies and two forms of grouped logistic regression models (multinomial and ordinal) in the investigation of patterns of successional change after agricultural land abandonment in Switzerland. Results indicated that both ordinal and nominal transitional change occurs in the landscape and that the use of different sampling regimes and modelling techniques as investigative tools yield different results. Synthesis and applications. Our multimodel inference identified successfully a set of consistently selected indicators of land cover change, which can be used to predict further change, including annual average temperature, the number of already overgrown neighbouring areas of land and distance to historically destructive avalanche sites. This allows for more reliable decision making and planning with respect to landscape management. Although both model approaches gave similar results, ordinal regression yielded more parsimonious models that identified the important predictors of land cover change more efficiently. Thus, this approach is favourable where land cover change pattern can be interpreted as an ordinal process. Otherwise, multinomial logistic regression is a viable alternative.

Mass effects mediate coexistence in competing shrews.

Recent developments in metacommunity theory have raised awareness that processes occurring at regional scales might interfere with local dynamics and affect conditions for the local coexistence of competing species. Four main paradigms are recognized in this context (namely, neutral, patch-dynamics, species-sorting, and mass-effect), which differ according to the role assigned to ecological or life-history differences among competing species, as well as to the relative time scale of regional vs. local dynamics. We investigated the patterns of regional and local coexistence of two species of shrews (Crocidura russula and Sorex coronatus) sharing a similar diet (generalist insectivores) over four generations, in a spatially structured habitat at the altitudinal limit of their distributions. Local populations were small, and regional dynamics were strong, with high rates of extinction and recolonization. Niche analysis revealed significant habitat differentiation on a few important variables, including temperature and availability of winter resting sites. In sites suitable for both species, we found instances of local coexistence with no evidence of competitive exclusion. Patterns of temporal succession did not differ from random, with no suggestion of a colonization-competition trade-off. Altogether, our data provide support for the mass-effect paradigm, where regional coexistence is mediated by specialization on different habitat types, and local coexistence by rescue effects from source sites. The strong regional dynamics and demographic stochasticity, together with high dispersal rates, presumably contributed to mass effects by overriding local differences in specific competitive abilities.

Mutagenesis and modelling of the alpha(1b)-adrenergic receptor highlight the role of the helix 3/helix 6 interface in receptor activation.

Computer simulations on a new model of the alpha1b-adrenergic receptor based on the crystal structure of rhodopsin have been combined with experimental mutagenesis to investigate the role of residues in the cytosolic half of helix 6 in receptor activation. Our results support the hypothesis that a salt bridge between the highly conserved arginine (R143(3.50)) of the E/DRY motif of helix 3 and a conserved glutamate (E289(6.30)) on helix 6 constrains the alpha1b-AR in the inactive state. In fact, mutations of E289(6.30) that weakened the R143(3.50)-E289(6.30) interaction constitutively activated the receptor. The functional effect of mutating other amino acids on helix 6 (F286(6.27), A292(6.33), L296(6.37), V299(6.40,) V300(6.41), and F303(6.44)) correlates with the extent of their interaction with helix 3 and in particular with R143(3.50) of the E/DRY sequence.

New developments and applications in modelling occupational exposure to airborne contaminants

Résumé: L'évaluation de l'exposition aux nuisances professionnelles représente une étape importante dans l'analyse de poste de travail. Les mesures directes sont rarement utilisées sur les lieux même du travail et l'exposition est souvent estimée sur base de jugements d'experts. Il y a donc un besoin important de développer des outils simples et transparents, qui puissent aider les spécialistes en hygiène industrielle dans leur prise de décision quant aux niveaux d'exposition. L'objectif de cette recherche est de développer et d'améliorer les outils de modélisation destinés à prévoir l'exposition. Dans un premier temps, une enquête a été entreprise en Suisse parmi les hygiénistes du travail afin d'identifier les besoins (types des résultats, de modèles et de paramètres observables potentiels). Il a été constaté que les modèles d'exposition ne sont guère employés dans la pratique en Suisse, l'exposition étant principalement estimée sur la base de l'expérience de l'expert. De plus, l'émissions de polluants ainsi que leur dispersion autour de la source ont été considérés comme des paramètres fondamentaux. Pour tester la flexibilité et la précision des modèles d'exposition classiques, des expériences de modélisations ont été effectuées dans des situations concrètes. En particulier, des modèles prédictifs ont été utilisés pour évaluer l'exposition professionnelle au monoxyde de carbone et la comparer aux niveaux d'exposition répertoriés dans la littérature pour des situations similaires. De même, l'exposition aux sprays imperméabilisants a été appréciée dans le contexte d'une étude épidémiologique sur une cohorte suisse. Dans ce cas, certains expériences ont été entreprises pour caractériser le taux de d'émission des sprays imperméabilisants. Ensuite un modèle classique à deux-zone a été employé pour évaluer la dispersion d'aérosol dans le champ proche et lointain pendant l'activité de sprayage. D'autres expériences ont également été effectuées pour acquérir une meilleure compréhension des processus d'émission et de dispersion d'un traceur, en se concentrant sur la caractérisation de l'exposition du champ proche. Un design expérimental a été développé pour effectuer des mesures simultanées dans plusieurs points d'une cabine d'exposition, par des instruments à lecture directe. Il a été constaté que d'un point de vue statistique, la théorie basée sur les compartiments est sensée, bien que l'attribution à un compartiment donné ne pourrait pas se faire sur la base des simples considérations géométriques. Dans une étape suivante, des données expérimentales ont été collectées sur la base des observations faites dans environ 100 lieux de travail différents: des informations sur les déterminants observés ont été associées aux mesures d'exposition des informations sur les déterminants observés ont été associé. Ces différentes données ont été employées pour améliorer le modèle d'exposition à deux zones. Un outil a donc été développé pour inclure des déterminants spécifiques dans le choix du compartiment, renforçant ainsi la fiabilité des prévisions. Toutes ces investigations ont servi à améliorer notre compréhension des outils des modélisations ainsi que leurs limitations. L'intégration de déterminants mieux adaptés aux besoins des experts devrait les inciter à employer cet outil dans leur pratique. D'ailleurs, en augmentant la qualité des outils des modélisations, cette recherche permettra non seulement d'encourager leur utilisation systématique, mais elle pourra également améliorer l'évaluation de l'exposition basée sur les jugements d'experts et, par conséquent, la protection de la santé des travailleurs. Abstract Occupational exposure assessment is an important stage in the management of chemical exposures. Few direct measurements are carried out in workplaces, and exposures are often estimated based on expert judgements. There is therefore a major requirement for simple transparent tools to help occupational health specialists to define exposure levels. The aim of the present research is to develop and improve modelling tools in order to predict exposure levels. In a first step a survey was made among professionals to define their expectations about modelling tools (what types of results, models and potential observable parameters). It was found that models are rarely used in Switzerland and that exposures are mainly estimated from past experiences of the expert. Moreover chemical emissions and their dispersion near the source have also been considered as key parameters. Experimental and modelling studies were also performed in some specific cases in order to test the flexibility and drawbacks of existing tools. In particular, models were applied to assess professional exposure to CO for different situations and compared with the exposure levels found in the literature for similar situations. Further, exposure to waterproofing sprays was studied as part of an epidemiological study on a Swiss cohort. In this case, some laboratory investigation have been undertaken to characterize the waterproofing overspray emission rate. A classical two-zone model was used to assess the aerosol dispersion in the near and far field during spraying. Experiments were also carried out to better understand the processes of emission and dispersion for tracer compounds, focusing on the characterization of near field exposure. An experimental set-up has been developed to perform simultaneous measurements through direct reading instruments in several points. It was mainly found that from a statistical point of view, the compartmental theory makes sense but the attribution to a given compartment could ñó~be done by simple geometric consideration. In a further step the experimental data were completed by observations made in about 100 different workplaces, including exposure measurements and observation of predefined determinants. The various data obtained have been used to improve an existing twocompartment exposure model. A tool was developed to include specific determinants in the choice of the compartment, thus largely improving the reliability of the predictions. All these investigations helped improving our understanding of modelling tools and identify their limitations. The integration of more accessible determinants, which are in accordance with experts needs, may indeed enhance model application for field practice. Moreover, while increasing the quality of modelling tool, this research will not only encourage their systematic use, but might also improve the conditions in which the expert judgments take place, and therefore the workers `health protection.

Modelling the distribution of bats in relation to landscape structure in a temperate mountain environment

1. Landscape modification is often considered the principal cause of population decline in many bat species. Thus, schemes for bat conservation rely heavily on knowledge about species-landscape relationships. So far, however, few studies have quantified the possible influence of landscape structure on large-scale spatial patterns in bat communities. 2. This study presents quantitative models that use landscape structure to predict (i) spatial patterns in overall community composition and (ii) individual species' distributions through canonical correspondence analysis and generalized linear models, respectively. A geographical information system (GIS) was then used to draw up maps of (i) overall community patterns and (ii) distribution of potential species' habitats. These models relied on field data from the Swiss Jura mountains. 3. Fight descriptors of landscape structure accounted for 30% of the variation in bat community composition. For some species, more than 60% of the variance in distribution could be explained by landscape structure. Elevation, forest or woodland cover, lakes and suburbs, were the most frequent predictors. 4. This study shows that community composition in bats is related to landscape structure through species-specific relationships to resources. Due to their nocturnal activities and the difficulties of remote identification, a comprehensive bat census is rarely possible, and we suggest that predictive modelling of the type described here provides an indispensable conservation tool.

High-resolution numerical modelling of flow-vegetation interactions

In this paper, we present and apply a new three-dimensional model for the prediction of canopy-flow and turbulence dynamics in open-channel flow. The approach uses a dynamic immersed boundary technique that is coupled in a sequentially staggered manner to a large eddy simulation. Two different biomechanical models are developed depending on whether the vegetation is dominated by bending or tensile forces. For bending plants, a model structured on the Euler-Bernoulli beam equation has been developed, whilst for tensile plants, an N-pendula model has been developed. Validation against flume data shows good agreement and demonstrates that for a given stem density, the models are able to simulate the extraction of energy from the mean flow at the stem-scale which leads to the drag discontinuity and associated mixing layer.

Evolution, maintenance and ecological consequences of genetic caste determination in " Pogonomyrmex " harvester ants

Summary Division of labor between reproducers (queens) and helpers (workers) is the main characteristic of social insect societies and at the root of their ecological success. Kin selection models predict that phenotypic differences between queens and workers should result from environmental rather than from genetic differences. However, genetic effects on queen and worker differentiation were found in two populations-of Pogonomyrmex harvester ants. Each of the two populations is composed of two genetically distinct lineages. Queens (which can be of either lineage) generally mate with males of their own and of the alternate lineage and produce two types of female offspring, those fertilized by males of the queens' lineage which develop into queens and those fertilized by males of the alternate lineage which develop into workers. All four lineages were further suggested to be themselves of hybrid origin between-the species P: barbatus and P. rugosus, in which queens and workers do not differ genetically. In a first set of experiments, we tested if female caste determination (the differentiation into queens and workers) in the lineages was genetically hardwired and if it was associated with costs in terms of the ability to optimally allocate resources to the production of queens and workers. To this end we first mated queens of-two lineages to a single male. Queens mated to a male of the alternate lineage successfully raised worker offspring whereas queens mated to a male of their own lineage almost always failed to produce workers. This reveals that pure-lineage individuals have lost the ability to develop into workers. Second, we analyzed offspring produced by naturally mated queens. During the stage of colony founding when only workers are produced, naturally mated queens laid a high proportion of pure-lineage eggs but the large majority of these eggs failed to develop. As a consequence, the number of offspring produced by incipient colonies decreased linearly with the proportion of pure-lineage eggs laid by queens. Moreover, queens of the lineage most commonly represented in a given population produced more pure-lineage eggs, in line with the view that they mate randomly with the two types of males and indiscriminately use their sperm. Altogether these results predict frequency-dependent founding success for pairs of lineages because queens of the more common lineage will produce more pure-lineage eggs and their colonies be less successful during the stage of colony founding. To describe the distribution of populations characterized with genetic caste determination relative to the populations with environmental caste determination we genotyped queens and workers collected during a large survey of -additional populations. Genetic caste determination associated with pairs of interbreeding lineages was frequent and widespread in the studied range and we identified four additional lineages displaying genetic caste determination. Overall, there were thus eight highly differentiated lineages with genetic caste determination. These lineages always co-occurred in the same complementary lineage pairs. Three of the four lineage pairs appeared to have a common origin, while their relationship with the forth could not be resolved. The genetic survey also revealed that, in addition to being genetically isolated from one another, all eight lineages were genetically distinct from P. rugosus and P. barbatus, even when colonies of interbreeding lineages co-occurred with colonies of either putative parent at the same site. This raised the question of the mechanisms involved in the reproductive isolation between the lineages and the parental species and between the two lineages of a lineage pair. At a site where one lineage pair co-occurred with P. rugosus, we identified two pre-zygotic mechanisms (differences in timing for mating flights between P. rugosus and the lineage pair and assortative mating) and one post-zygotic mechanism (high levels of hybrid unviablility) which in combination may largely account for the reproductive isolation between the lineages and their parental species. The mechanisms accounting for the reproductive isolation between the two lineages of a lineage pair varied across lineage pairs. In one lineage pair, inter-lineage individuals exclusively occurred in the sterile worker caste, raising the possibility that inter-lineage eggs have completely lost the ability to develop into queens in this lineage pair and that there is thus no opportunity for gene flow. In each of the three remaining lineage pairs, inter-lineage queens were produced by a minority of colonies. In these lineage pairs, colonies headed by inter-lineage queens failed to grow sufficiently to produce reproductive individuals which may account for the reproductive isolation between co-occurring lineages in three lineage pairs. In conclusion, the results of this thesis show that genetic caste determination is costly but widespread in Pogonomyrmex harvester ants. Reproductive isolation among the lineages and between the lineages and the parental species as well as frequency-dependent founding success for co-occurring lineages may contribute to the persistence of this extraordinary system. Résumé La division du travail entre individus reproducteurs (les reines) et individus non-reproducteurs (ouvrières) représente la caractéristique principale des sociétés d'insectes et est à la base de leur succès écologique. Des modèles de sélection de parentèle prédisent que les différences phénotypiques entre reines et ouvrières devraient provenir d'effets environnementaux plutôt que de différences génétiques. Malgré ce fait, des effets génétiques sur la différentiation entre reines et ouvrières ont été montrés dans deux populations de fourmis moissonneuses du genre Pogonomyrmex. Chacune des deux populations est composée de deux lignées génétiquement distinctes. Les reines de chaque lignée s'accouplent en général avec des mâles de leur propre lignée ainsi qu'avec des mâles de l'autre lignée et produisent deux types d'oeufs, ceux qui sont fécondés par les mâles de leur propre lignée qui se développent en nouvelles reines et ceux qui sont fécondés par les mâles de l'autre lignée qui se développent en ouvrières. Il a été suggéré que les lignées sont elles-mêmes des hybrides entre les deux espèces P. barbatus et P. rugosus. Dans ces deux espèces, les reines et ouvrières ne sont pas génétiquement distinctes. Dans une première série d'expériences, nous avons testé si la détermination de la caste femelle (le développement en reine ou en ouvrière) est génétiquement rigide et si elle est associée à des coûts en terme de capacité à allouer de façon optimale les ressources pour la production de reines et d'ouvrières. Pour cela nous avons accouplé des reines de deux lignées avec un seul mâle. Les reines accouplées avec un mâle de l'autre lignée ont élevé de nouvelles ouvrières avec succès alors que les reines accouplées avec un mâle de leur propre lignée ont presque toujours échoué à produire des ouvrières. Ceci montre que les individus de lignée pure ont perdu la capacité de se développer en ouvrière. Deuxièmement, nous avons analysé la descendance de reines qui se sont accouplées naturellement. Durant le stade de fondation de la colonie, où seules des ouvrières sont élevées, les reines accouplées naturellement ont pondu une grande proportion d'oeufs de lignée pure mais la majorité de ces derniers ne se sont pas développés. En conséquence, le nombre de descendants produits par des colonies fondatrices diminuait linéairement avec la proportion des oeufs de lignée pure pondus par la reine en accord avec l'hypothèse que les reines s'accouplent au hasard avec les deux types de mâles et utilisent leur sperme aléatoirement. Dans l'ensemble; ces résultats prédisent un succès de fondation fréquence-dépendant pour les deux lignées, car les reines de la lignée la plus fréquente produiront .plus d'oeufs de lignée pure et leurs colonies auront moins de succès lors de la fondation de colonies par rapport aux colonies de la lignée la moins fréquente. Pour décrire la distribution des-populations caractérisées par une détermination génétique des castes par rapport aux populations caractérisées par une détermination environnementale des castes, nous avons génotypé des reines et des ouvrières qui ont été collectées lors d'une analyse de populations supplémentaires. La détermination génétique des castes associée à des croisements entre lignées est fréquente et largement répartie dans l'aire étudiée. Nous avons identifié quatre lignées supplémentaires, ayant une détermination génétique des castes, pour un total de huit lignées. Ces huit lignées forment quatre paires de lignées et on ne trouve jamais deux lignées de paires différentes, dans une population. Trois des quatre paires de lignées s'avèrent avoir une origine commune alors que leur relation avec la quatrième paire de lignées n'a pas pu être résolue. L'analyse génétique de populations supplémentaires a également révélé qu'en plus d'être génétiquement isolées les unes des autres, les huit lignées sont génétiquement distinctes de P. rugosus et P. barbatus même si les colonies d'une paire de lignées se trouvent en sympatrie avec l'une ou l'autre des espèces parentales. Ceci relève la question des mécanismes impliqués dans l'isolation reproductive entre les lignées et les espèces parentales ainsi qu'entre les deux lignées d'une paire. En étudiant un site où une paire de lignées se trouve en sympatrie avec P. rugosus, nous avons identifié deux mécanismes pré-zygotiques (des différences dans le timing du vol nuptial entre P. rugosus et les lignées et des accouplements assortis) ainsi qu'un mécanisme post-zygotique (un niveau élevé de non-viabilité des hybrides). En combinaison, ces mécanismes peuvent largement expliquer l'isolement reproductif entre les lignées et leurs espèces parentales. Les mécanismes contribuant à l'isolement reproductif entre les deux lignées d'une paire variaient entre paires de lignées. Dans une paire, les individus de génotype inter-lignée se trouvent uniquement dans la caste stérile des ouvrières, suggérant qu'il n'y a pas d'opportunité pour avoir du flux de gènes entre les deux lignées ce cette paire. Dans chacune des trois autres paires de lignées des nouvelles reines de génotype inter-lignée sont produites par une minorité de colonies. Par contre, les colonies avec une reine mère de génotype inter-lignée ne se développent pas suffisamment pour produire des individus reproducteurs. Ceci peut donc expliquer pourquoi il n'y a pas de flux de gènes entre les deux lignées de trois paires. En conclusion, les résultats de cette thèse montrent que la détermination génétique de la caste est coûteuse mais très répandue chez les fourmis. moissonneuses du genre Pogonomyrmex. L'isolement reproductif des lignées entre elles et avec les espèces parentales, ainsi qu'un succès de fondation fréquence-dépendant contribuent à la persistance de ce système extraordinaire.

Population pharmacokinetic modelling and evaluation of different dosage regimens for darunavir and ritonavir in HIV-infected individuals.

OBJECTIVES: Darunavir is a protease inhibitor that is administered with low-dose ritonavir to enhance its bioavailability. It is prescribed at standard dosage regimens of 600/100 mg twice daily in treatment-experienced patients and 800/100 mg once daily in naive patients. A population pharmacokinetic approach was used to characterize the pharmacokinetics of both drugs and their interaction in a cohort of unselected patients and to compare darunavir exposure expected under alternative dosage regimens. METHODS: The study population included 105 HIV-infected individuals who provided darunavir and ritonavir plasma concentrations. Firstly, a population pharmacokinetic analysis for darunavir and ritonavir was conducted, with inclusion of patients' demographic, clinical and genetic characteristics as potential covariates (NONMEM(®)). Then, the interaction between darunavir and ritonavir was studied while incorporating levels of both drugs into different inhibitory models. Finally, model-based simulations were performed to compare trough concentrations (Cmin) between the recommended dosage regimen and alternative combinations of darunavir and ritonavir. RESULTS: A one-compartment model with first-order absorption adequately characterized darunavir and ritonavir pharmacokinetics. The between-subject variability in both compounds was important [coefficient of variation (CV%) 34% and 47% for darunavir and ritonavir clearance, respectively]. Lopinavir and ritonavir exposure (AUC) affected darunavir clearance, while body weight and darunavir AUC influenced ritonavir elimination. None of the tested genetic variants showed any influence on darunavir or ritonavir pharmacokinetics. The simulations predicted darunavir Cmin much higher than the IC50 thresholds for wild-type and protease inhibitor-resistant HIV-1 strains (55 and 550 ng/mL, respectively) under standard dosing in >98% of experienced and naive patients. Alternative regimens of darunavir/ritonavir 1200/100 or 1200/200 mg once daily also had predicted adequate Cmin (>550 ng/mL) in 84% and 93% of patients, respectively. Reduction of darunavir/ritonavir dosage to 600/50 mg twice daily led to a 23% reduction in average Cmin, still with only 3.8% of patients having concentrations below the IC50 for resistant strains. CONCLUSIONS: The important variability in darunavir and ritonavir pharmacokinetics is poorly explained by clinical covariates and genetic influences. In experienced patients, treatment simplification strategies guided by drug level measurements and adherence monitoring could be proposed.

Phylogenetic conservatism in plant phenology

Phenological events - defined points in the life cycle of a plant or animal - have been regarded as highly plastic traits, reflecting flexible responses to various environmental cues. The ability of a species to track, via shifts in phenological events, the abiotic environment through time might dictate its vulnerability to future climate change. Understanding the predictors and drivers of phenological change is therefore critical. Here, we evaluated evidence for phylogenetic conservatism - the tendency for closely related species to share similar ecological and biological attributes - in phenological traits across flowering plants. We aggregated published and unpublished data on timing of first flower and first leaf, encompassing 4000 species at 23 sites across the Northern Hemisphere. We reconstructed the phylogeny for the set of included species, first, using the software program Phylomatic, and second, from DNA data. We then quantified phylogenetic conservatism in plant phenology within and across sites. We show that more closely related species tend to flower and leaf at similar times. By contrasting mean flowering times within and across sites, however, we illustrate that it is not the time of year that is conserved, but rather the phenological responses to a common set of abiotic cues. Our findings suggest that species cannot be treated as statistically independent when modelling phenological responses.Synthesis. Closely related species tend to resemble each other in the timing of their life-history events, a likely product of evolutionarily conserved responses to environmental cues. The search for the underlying drivers of phenology must therefore account for species' shared evolutionary histories.