Evaluation et validation des équations de prédiction du risque de diabète dans la population lausannoise : étude CoLaus

Le diabète de type 2 est une maladie fréquente et en augmentation dans le monde entier. Malheureusement, elle est souvent diagnostiquée à un stade ou des complications sont déjà apparues. Depuis quelques années, des scores ont été développés pour identifier les sujets à risque de développer cette maladie. L'utilisation d'un tel score par le praticien pourrait amener ces patients à des mesures préventives, telles que le changement d'hygiène de vie, ou la prescription d'un traitement médicamenteux. Le but de notre étude est de comparer et de valider différents scores de risques de diabète de type 2 et de déterminer leur capacité à prédire la survenue de cette maladie dans la population de la cohorte CoLaus. Les premiers résultats, en étude transversale, ont tout d'abord montré de grandes différences quant à la population à risque d'un score à l'autre. En effet, le nombre de personnes à traiter varie considérablement selon la méthode utilisée. Ces différents scores ont donc nécessité une validation prospective. Ces résultats ont fait l'objet d'une publication (Schmid et col, Diabetes Care. 2011 Aug;34(8):1863-8). Au moyen des données du suivi à 5 ans, il est sorti qu'un score de risque utilisant des variables biologiques et cliniques, ainsi qu'un score utilisant des variables uniquement cliniques, obtenaient de très bon résultats quant à la prédiction du diabète de type 2. En effet, un des scores testés donne une valeur prédictive positive d'environ 20% à 5 ans, ce qui signifie qu'un patient « détecté » sur 5 pourrait bénéficier d'une intervention précoce. Toutefois, ces résultats concernent la population lausannoise et ne sont donc pas forcément applicables à l'ensemble de la population suisse. De plus, de plus amples études sont nécessaires évaluer l'efficacité d'un tel score dans la prévention du diabète en Suisse. Ces résultats ont fait l'objet d'une seconde publication (Schmid et col, Arch Intern Med. 2012 Jan 23;172(2):188-9). Dans un troisième volet de l'étude, l'impact de marqueurs génétiques a été évalué dans un sous- groupe de la population CoLaus. Les résultats n'ont toutefois montré qu'une très faible amélioration de la prédiction du risque en utilisant ces marqueurs. Ceci devrait nous encourager à intensifier les efforts de prévention sur le style de vie pour toute la population, plutôt qu'une approche ciblée sur les personnes génétiquement prédisposées. Ces résultats ont fait l'objet d'une troisième publication (Schmid et col, J Clin Endocrinol Metab. 2012 Apr 24. [Epub ahead of print]). La même démarche méthodologique a été utilisée pour évaluer l'importance pronostique de plusieurs marqueurs inflammatoires (interleukines 1 et 6, Τ Ν F-, protéine C-réactive) hépatiques (GT) ou adipocytaires (leptine et adiponectine) dans la survenue du diabète. Ces résultats sont actuellement soumis au Journal of Clinical Endocrinology and Metabolism).

CHARACTERIZATION OF EPSTEIN-BARR VIRUS-ENCODED LATENT MEMBRANE PROTEIN 1

Le virus d'Epstein-Barr (EBV), un virus de la famille des gammaherpesvirus, infecte plus de 95% de la population adulte mondiale. EBV est associé à plusieurs types de cancers dont le lymphome de Hodgkin, le lymphome de Burkitt et le carcinome nasopharyngé. La protéine membranaire de latence 1 (LMP1), l'oncogène principal d'EBV, est une protéine membranaire intégrale composée d'une petite extrémité N-terminale cytoplasmique, six segments transmembranaires (TMs) lié par de petites boucles et un long domaine C-terminale cytoplasmique. Le gène de LMP1, BNLF-1, est très polymorphe et plusieurs variants de la protéine LMP1 ont été décrits. Parmi les variants de LMP1 la majeure différence décrite est leur capacité à activer le facteur de transcription NF-κB. Nous avons défini des polymorphismes permettant aux variants d'avoir une activation accrue de NF-κB comparé au prototype B95-8 LMP1. Tous les polymorphismes cruciaux identifiés dans notre étude se trouvent dans les TMs 4 et 5 de LMP1. Nous avons étudié l'implication de chaque paire de TMs dans l'association à la membrane, l'auto-agrégation, la liaison aux partenaires cellulaires de LMP1 TRAF3 et β-TrCP, ainsi que pour NF-κB. De plus, nous avons décrit un nouveau rôle pour LMP1 consistant à inhiber l'activation contrôlée par MAVS de ISRE et du promoteur d'IFNβ. En résumé, nous avons observé que les différentes paires de TMs, ainsi que les deux boucles intracellulaires, ne sont pas équivalents. Dans l'ensemble, notre étude a montré que les TMs jouent un rôle clé dans les interactions protéine-protéine et la signalisation et qu'ils peuvent être considérés comme des régulateurs essentiels des activités de LMP1.

Le soutien aux autorités lors des votations fédérales de 1971 à 1999 dans les cantons suisses

Cet article vise à expliquer les variations inter-cantonales du succès des autorités ainsi que du soutien populaire qui leur est accordé en votation fédérale entre 1971 et 1999. A cette fin, l'auteur cherche à intégrer deux courants de recherches dans un cadre d'analyse unique en déterminant l'impact respectif des mots d'ordre des partis politiques et des facteurs socio-structurels sur le résultat de vote ainsi que leur effet conjugué. L'analyse empirique, qui est situé au niveau cantonal, est réalisée à l'aide d'une régression logistique et d'un modèle d'équations structurelles. Les résultats obtenus montrent que le caractère alémanique, catholique et agraire d'un canton se traduit par une réticence face à la politique des autorités, qui peut s'amplifier selon le type juridique de la votation et l'enjeu politique soumis au vote. Plus qu'ailleurs, les partis politiques de ces cantons s'alignent cependant fermement derrière la recommandation de vote des autorités et exercent, de ce fait, un effet atténuateur qui contribue au nivellement des disparités inter-cantonales. Ces résultats suggèrent ainsi que les partis politiques cantonaux réussissent, dans une certaine mesure, à contrôler le résultat de vote par le biais d'une concertation accrue durant la phase référendaire.

Utilisation d'EMLA® (lidocaïne/prilocaïne)de façon routinière lors de vaccinations chez l'enfant : anodin?

L'utilisation d'EMLA 5%®, une association de deux anesthésiques locaux, la lidocaïne et la prilocaïne, sous forme de crème ou de patch, s'est répandue de façon assez générale. Son utilisation lors de gestes invasifs, tels que ponctions veineuses ou injections intramusculaires, petite chirurgie comme ablation de verrues plantaires ou biopsies cutanées, s'est étendue à toutes sortes d'indications, telles qu'administration de vaccins ou, en dehors d'une indication médicale, lors de piercing ou de pose de boucles d'oreilles. S'il est vrai que l'accent a été mis ces dernières années sur l'importance d'une bonne antalgie, en particulier chez les jeunes enfants, la généralisation de l'utilisation d'anesthésiques topiques pour des gestes anodins doit faire évoquer les risques potentiels associés à de telles molécules.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Étude du mécanisme de transport des cations par la Na, K-ATPase et de son implication dans la migraine familiale hémiplégique de type 2

Résumé La Na,K-ATPase est une protéine transmembranaire, présente dans toutes les cellules de mammifères et indispensable à la viabilité cellulaire. Elle permet le maintien des gradients sodiques et potassiques à l'origine du potentiel membranaire en transportant 3 Na+ en dehors de la cellule contre 2 K+, grâce à l'énergie fournie par l'hydrolyse d'une molécule d'ATP. Le potentiel membranaire est indispensable au maintien de l'excitabilité cellulaire et à la transmission de l'influx nerveux. Il semblerait que la Na,K-ATPase soit liée à l'hypertension et à certains troubles neurologiques comme la Migraine Familiale Hémiplégique (1VIFH). La MFH est une forme de migraine avec aura, qui se caractérise par une hémiparésie. Cette forme de migraine est très rare. Elle se transmet génétiquement sur un mode autosomique dominant. Plusieurs mutations localisées dans le gène de la Na,K-ATPase ont été identifiées durant ces 3 dernières années. C'est la première fois qu'une maladie génétique est associée au gène de la Na,K-ATPase. La compréhension du fonctionnement de cette protéine peut donner des informations sur les mécanismes conduisant à ces pathologies. On sait que la fonction d'une protéine est liée à sa structure. L'étude de sa fonction nécessite donc l'étude de sa structure. Alors que la structure de la SERCA a été déterminée à haute résolution, par cristallographie, celle de la Na,K-ATPase ne l'est toujours pas. Mais ces 2 ATPases présentent une telle homologie qu'un modèle de la Na,K-ATPase a pu être élaboré à partir de la structure de la SERCA. Les objectifs de cette étude sont d'une part, de comprendre le contrôle de l'accessibilité du K+ extracellulaire àses sites de liaison. Pour cela, nous avons ciblé cette étude sur la 2ìème et la 31eme boucle extracellulaire, qui relient respectivement les segments transmembranaires (STM) 3-4 et 5-6. Le choix s'est porté sur ces 2 boucles car elles bordent le canal des cations formés des 4ième' Sième et 6'ème hélices. D'autre part, nous avons également essayer de comprendre les effets des mutations, liées à la Migraine Familiale Hémiplégique de type 2 (MFH2), sur la fonctionnalité de la Na,K-ATPase. Alors que les STM et les domaines cytoplasmiques sont relativement proches entre la Na,KATPase et la SERCA, les boucles extracellulaires présentent des différences. Le modèle n'est donc pas une approche fiable pour déterminer la structure et la fonction des régions extracellulaires. Nous avons alors utilisé une approche fonctionnelle faisant appel à la mutation dirigée puis à l'étude de l'activité fonctionnelle de la Na,K ATPase par électrophysiologie sur des ovocytes de Xenopus. En conclusion, nous pouvons dire que la troisième boucle extracellulaire participerait à la structure de la voie d'entrée des cations et que la deuxième boucle extracellulaire semble impliquée dans le contrôle de l'accessibilité des ions K+àses sites de liaison. Concernant les mutations associées à la MFH2, nos résultats ont montré une forte diminution de l'activité fonctionnelle de la pompe Na,K, inférieure aux conditions physiologiques de fonctionnement, et pour une des mutations nous avons observés une diminution de l'affmité apparente au K+ externe. Nous poumons faire l'hypothèse que l'origine pathologique de la migraine est liée à une diminution de l'activité de la pompe à Na+. Summary The Na,K-ATPase is a transmembrane protein, present in all mammalian cells and is necessary for the viability of the cells. It maintains the gradients of Na+ and K+ involved in the membrane potential, by transporting 3Na+ out the cell, and 2K+ into the cell, using the energy providing from one ATP molecule hydrolysis. The membrane potential is necessary for the cell excitability and for the transmission of the nervous signal. Some evidence show that Na,K-ATPase is involved in hypertension and neurological disorders like the Familial Hemiplegic Migraine (FHM). La FHM is a rare form of migraine characterised by aura and hemiparesis and an autosomal dominant transmission. Several mutations linked to the Na,KATPase gene have been identified during these 3 last years. It's the first genetic disorder associated with the Na,K-ATPase gene. Understand the function of this protein is important to elucidate the mechanisms implicated in these pathologies. The function of a protein is linked with its structure. Thus, to know the function of a protein, we need to know its structure. While the Ca-ATPase (SERCA) has been crystallised with a high resolution, the structure of the Na,K-ATPase is not known. Because of the great homology between these 2 ATPases, a model of the Na,K-ATPase was realised by comparing with the structure of the SERCA. The aim of this study is on one side, understand the control of the extracellular K+ accessibility to their binding sites. Because of theirs closed proximity with the cation pathway, located between the 4th, 5th and 6th helices, we have targeted this study on the 2nd and the 3rd extracellular loops linking respectively the transmembrane segment (TMS) 3 and 4, and the TMS 5 and 6. And on the other side, we have tried to understand the functional effects of mutations linked with the Familial Hemiplegic Migraine Type 2 (FHM2). In contrast with the transmembrane segments and the cytoplasmic domains, the extracellular loops show lots of difference between Na,K-ATPase and SERCA, the model is not a good approach to know the structure and the function of the extracellular loops. Thus, we have used a functional approach consisting in directed mutagenesis and the study of the functional activity of the Na,K-ATPase by electrophysiological techniques with Xenopus oocytes. In conclusion, we have demonstrated that the third extracellular loop could participate in the structure of the entry of the cations pathway and that the second extracellular loop could control the K+ accessibility to their binding sites. Concerning the mutations associated with the FHM2, our results showed a strong decrease in the functional activity of the Na,K-pump under physiological conditions and for one of mutations, induce a decrease in the apparent external K+ affinity. We could make the hypothesis that the pathogenesis of migraine is related to the decrease in Na,K-pump activity. Résumé au large publique De la même manière que l'assemblage des mots forme des phrases et que l'assemblage des phrases forme des histoires, l'assemblage des cellules forme des organes et l'ensemble des organes constitue les êtres vivants. La fonction d'une cellule dans le corps humain peut se rapprocher de celle d'une usine hydroélectrique. La matière première apportée est l'eau, l'usine électrique va ensuite convertir l'eau en énergie hydraulique pour fournir de l'électricité. Le fonctionnement de base d'une cellule suit le même processus. La cellule a besoin de matières premières (oxygène, nutriments, eau...) pour produire une énergie sous forme chimique, l'ATP. Cette énergie est utilisée par exemple pour contracter les muscles et permet donc à l'individu de se déplacer. Morphologiquement la cellule est une sorte de petit sac rempli de liquide (milieu intracellulaire) baignant elle-même dans le liquide (milieu extracellulaire) composant le corps humain (un adulte est constitué environ de 65 % d'eau). La composition du milieu intracellulaire est différente de celle du milieu extracellulaire. Cette différence doit être maintenue pour que l'organisme fonctionne correctement. Une des différences majeures est la quantité de sodium. En effet il y a beaucoup plus de sodium à l'extérieur qu'à l'intérieur de la cellule. Bien que l'intérieur de la cellule soit isolé de l'extérieur par une membrane, le sodium arrive à passer à travers cette membrane, ce qui a tendance à augmenter la quantité de sodium dans la cellule et donc à diminuer sa différence de concentration entre le milieu extracellulaire et le milieu intracellulaire. Mais dans les membranes, il existe des pompes qui tournent et dont le rôle est de rejeter le sodium de la cellule. Ces pompes sont des protéines connues sous le nom de pompe à sodium ou Na,K-ATPase. On lui attribue le nom de Na,K-ATPase car en réalité elle rejette du sodium (Na) et en échange elle fait entrer dans la cellule du potassium (K), et pour fonctionner elle a besoin d'énergie (ATP). Lorsque les pompes à sodium ne fonctionnent pas bien, cela peut conduire à des maladies. En effet la Migraine Familiale Hémiplégique de type 2, est une migraine très rare qui se caractérise par l'apparition de la paralysie de la moitié d'un corps avant l'apparition du mal de tête. C'est une maladie génétique (altération qui modifie la fonction d'une protéine) qui touche la pompe à sodium située dans le cerveau. On a découvert que certaines altérations (mutations) empêchent les pompes à sodium de fonctionner correctement. On pense alors que le développement des migraines est en partie dû au fait que ces pompes fonctionnent moins bien. Il est important de bien connaître la fonction de ces pompes car cela permet de comprendre des mécanismes pouvant conduire à certaines maladies, comme les migraines. En biologie, la fonction d'une protéine est étudiée à travers sa structure. C'est pourquoi l'objectif de cette thèse a été d'étudier la structure de la Na,K-ATPase afin de mieux comprendre son mécanisme d'action.

Régulation du canal à sodium épithélial par le sodim extracellulaire et développement d'un microsystème intégré pour l'étude électrophysiologique sur ovocytes de "Xenopus laevis"

Résumé : La première partie de ce travail de thèse est consacrée au canal à sodium épithélial (ENaC), l'élément clé du transport transépithélial de Na+ dans le néphron distal, le colon et les voies aériennes. Ce canal est impliqué dans certaines formes génétiques d'hypo- et d'hypertension (PHA I, syndrome de Liddle), mais aussi, indirectement, dans la mucoviscidose. La réabsorption transépithéliale de Na+ est principalement régulée par des hormones (aldostérone, vasopressine), mais aussi directement par le Na+, via deux phénomènes distincts, la « feedback inhibition » et la « self-inhibition » (SI). Ce second phénomène est dépendant de la concentration de Na+ extracellulaire, et montre une cinétique rapide (constante de temps d'environ 3 s). Son rôle physiologique serait d'assurer l'homogénéité de la réabsorption de Na+ et d'empêcher que celle-ci soit excessive lorsque les concentrations de Na+ sont élevées. Différents éléments appuient l'hypothèse de la présence d'un site de détection de la concentration du Na+ extracellulaire sur ENaC, gouvernant la SI. L'objectif de ce premier projet est de démontrer l'existence du site de détection impliqué dans la SI et de déterminer ses propriétés physiologiques et sa localisation. Nous avons montré que les caractéristiques de la SI (en termes de sélectivité et affinité ionique) sont différentes des propriétés de conduction du canal. Ainsi, nos résultats confirment l'hypothèse de l'existence d'un site de détection du Na+ (responsable de la transmission de l'information au mécanisme de contrôle de l'ouverture du canal), différent du site de conduction. Par ailleurs, ce site présente une affinité basse et indépendante du voltage pour le Na+ et le Li+ extracellulaires. Le site semble donc être localisé dans le domaine extracellulaire, plutôt que transmembranaire, de la protéine. L'étape suivante consiste alors à localiser précisément le site sur le canal. Des études précédentes, ainsi que des résultats préliminaires récemment obtenus, mettent en avant le rôle dans la self-inhibition du premiers tiers des boucles extracellulaires des sous-unités α et γ du canal. Le second projet tire son origine des limitations de la méthode classique pour l'étude des canaux ioniques, après expression dans les ovocytes de Xenopus laevis, par la méthode du voltage-clamp à deux électrodes, en particulier les limitations dues à la lenteur des échanges de solutions. En outre, cette méthode souffre de nombreux désavantages (manipulations délicates et peu rapides, grands volumes de solution requis). Plusieurs systèmes améliorés ont été élaborés, mais aucun ne corrige tous les désavantages de la méthode classique Ainsi, l'objectif ici est le développement d'un système, pour l'étude électrophysiologique sur ovocytes, présentant les caractéristiques suivantes : manipulation des cellules facilitée et réduite, volumes de solution de perfusion faibles et vitesse rapide d'échange de la perfusion. Un microsystème intégré sur une puce a été élaboré. Ces capacités de mesure ont été testées en utilisant des ovocytes exprimant ENaC. Des résultats similaires (courbes IV, courbes dose-réponse au benzamil) à ceux obtenus avec le système traditionnel ont été enregistrés avec le microsystème. Le temps d'échange de solution a été estimé à ~20 ms et des temps effectifs de changement ont été déterminés comme étant 8 fois plus court avec le nouveau système comparé au classique. Finalement, la SI a été étudiée et il apparaît que sa cinétique est 3 fois plus rapide que ce qui a été estimé précédemment avec le système traditionnel et son amplitude de 10 à 20 % plus importante. Le nouveau microsystème intégré apparaît donc comme adapté à la mesure électrophysiologique sur ovocytes de Xenopus, et possèdent des caractéristiques appropriées à l'étude de phénomènes à cinétique rapide, mais aussi à des applications de type « high throughput screening ». Summary : The first part of the thesis is related to the Epithelial Sodium Channel (ENaC), which is a key component of the transepithelial Na+ transport in the distal nephron, colon and airways. This channel is involved in hypo- and hypertensive syndrome (PHA I, Liddle syndrome), but also indirectly in cystic fibrosis. The transepithelial reabsorption of Na+ is mainly regulated by hormones (aldosterone, vasopressin), but also directly by Na+ itself, via two distinct phenomena, feedback inhibition and self-inhibition. This latter phenomenon is dependant on the extracellular Na+ concentration and has rapid kinetics (time constant of about 3 s). Its physiological role would be to prevent excessive Na+ reabsorption and ensure this reabsorption is homogenous. Several pieces of evidence enable to propose the hypothesis of an extracellular Na+ sensing site on ENaC, governing self-inhibition. The aim of this first project is to demonstrate the existence of the sensing site involved in self-inhibition and to determine its physiological properties and localization. We show self-inhibition characteristics (ionic selectivity and affinity) are different from the conducting properties of the channel. Our results support thus the hypothesis that the Na+ sensing site (responsible of the transmission of the information about the extracellular Na+ concentration to the channel gating mechanism), is different from the channel conduction site. Furthermore, the site has a low and voltage-insensitive affinity for extracellular Na+ or Li+. This site appears to be located in the extracellular domain rather than in the transmembrane part of the channel protein. The next step is then to precisely localize the site on the channel. Some previous studies and preliminary results we recently obtained highlight the role of the first third of the extracellular loop of the α and γ subunits of the channel in self-inhibition. The second project originates in the limitation of the classical two-electrode voltageclamp system classically used to study ion channels expressed in Xenopus /aevis oocytes, in particular limitations related to the slow solution exchange time. In addition, this technique undergoes several drawbacks (delicate manipulations, time consumption volumes). Several improved systems have been built up, but none corrected all these detriments. The aim of this second study is thus to develop a system for electrophysiological study on oocytes featuring an easy and reduced cell handling, small necessary perfusion volumes and fast fluidic exchange. This last feature establishes the link with the first project, as it should enable to improve the kinetics analysis of self-inhibition. A PDMS chip-based microsystem has been elaborated. Its electrophysiological measurement abilities have been tested using oocytes expressing ENaC. Similar measurements (IV curves of benzamil-sensitive currents, benzamil dose-response curves) have been obtained with this system, compared to the traditional one. The solution exchange time has been estimated at N20 ms and effective exchange times (on inward currents) have been determined as 8 times faster with the novel system compared to the classical one. Finally, self-inhibition has been studied and it appears its kinetics is 3 times faster and its amplitude 10 to 20 % higher than what has been previously estimated with the traditional system. The novel integrated microsystem appears therefore to be convenient for electrophysiological measurement on Xenopus oocytes, and displays features suitable for the study of fast kinetics phenomenon, but also high throughput screening applications. Résumé destiné large public : Le corps humain est composé d'organes, eux-mêmes constitués d'un très grand nombre de cellules. Chaque cellule possède une paroi appelée membrane cellulaire qui sépare l'intérieur de cette cellule (milieu intracellulaire) du liquide (milieu extracellulaire) dans lequel elle baigne. Le maintien de la composition stable de ce milieu extracellulaire est essentiel pour la survie des cellules et donc de l'organisme. Le sodium est un des composants majeurs du milieu extracellulaire, sa quantité dans celui-ci doit être particulièrement contrôlée. Le sodium joue en effet un rôle important : il conditionne le volume de ce liquide extracellulaire, donc, par la même, du sang. Ainsi, une grande quantité de sodium présente dans ce milieu va de paire avec une augmentation du volume sanguin, ce qui conduit l'organisme à souffrir d'hypertension. On se rend donc compte qu'il est très important de contrôler la quantité de sodium présente dans les différents liquides de l'organisme. Les apports de sodium dans l'organisme se font par l'alimentation, mais la quantité de sodium présente dans le liquide extracellulaire est contrôlée de manière très précise par le rein. Au niveau de cet organe, on appelle urine primaire le liquide résultant de la filtration du sang. Elle contient de nombreuses substances, des petites molécules, dont l'organisme a besoin (sodium, glucose...), qui sont ensuite récupérées dans l'organe. A la sortie du rein, l'urine finale ne contient plus que l'excédent de ces substances, ainsi que des déchets à éliminer. La récupération du sodium est plus ou moins importante, en fonction des ajustements à apporter à la quantité présente dans le liquide extracellulaire. Elle a lieu grâce à la présence de protéines, dans les membranes des cellules du rein, capables de le transporter et de le faire transiter de l'urine primaire vers le liquide extracellulaire, qui assurera ensuite sa distribution dans l'ensemble de l'organisme. Parmi ces protéines « transporteurs de sodium », nous nous intéressons à une protéine en particulier, appelée ENaC. Il a été montré qu'elle jouait un rôle important dans cette récupération de sodium, elle est en effet impliquée dans des maladies génétiques conduisant à l'hypo- ou à l'hypertension. De précédents travaux ont montré que lorsque le sodium est présent en faible quantité dans l'urine primaire, cette protéine permet d'en récupérer une très grande partie. A l'inverse, lorsque cette quantité de sodium dans l'urine primaire est importante, sa récupération par le biais d'ENaC est réduite. On parle alors d'autorégulation : la protéine elle-même est capable d'adapter son activité de transport en fonction des conditions. Ce phénomène d'autorégulation constitue a priori un mécanisme préventif visant à éviter une trop grande récupération de sodium, limitant ainsi les risques d'hypertension. La première partie de ce travail de thèse a ainsi consisté à clarifier le mécanisme d'autorégulation de la protéine ENaC. Ce phénomène se caractérise en particulier par sa grande vitesse, ce qui le rend difficile à étudier par les méthodes traditionnelles. Nous avons donc, dans une deuxième partie, développé un nouveau système permettant de mieux décrire et analyser cette « autorégulation » d'ENaC. Ce second projet a été mené en collaboration avec l'équipe de Martin Gijs de l'EPFL.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Dissection of hormone-responsive plasma membrane to nucleus signaling of Bravis Radix : its role in vascular differentiation and root meristem growth

AbstractPlants continuously grow during their complete life span and understanding the mechanisms that qualitatively regulate their traits remains a challenging topic in biology. The hormone auxin has been identified as a crucial molecule for shaping plant growth, as it has a role in most developmental processes. In the root, the directional, so-called polar transport of auxin generates a peak of concentration that specifies and maintains the stem cell niche and a subsequent gradient of decreasing concentration that also regulates cell proliferation and differentiation. For these reasons, auxin is considered the main morphogen of the root, as it is fundamental for its organization and maintenance. Recently, in Arabidopsis thaliana, a natural variation screen allowed the discovery of BREVIS RADIX (BRX) gene as a limiting factor for auxin responsive gene expression and thus for root growth.In this study, we discovered that BRX is a direct target of auxin that positively feeds back on auxin signaling, as a transcriptional co-regulator, through interaction with the Auxin Response Factor (ARF) MONOPTEROS (MP), modulating the auxin gene response magnitude during the transition between division and differentiation in the root meristem. Moreover, we provide evidence that BRX is activated at the plasma membrane level as an associated protein before moving into the nucleus to modulate cellular growth.To investigate the discrepancy between the auxin concentration and the expression pattern of its downstream targets, we combined experimental and computational approaches. Expression profiles deviating from the auxin gradient could only be modeled after intersection of auxin activity with the observed differential endocytosis pattern and with positive auto- regulatory feedback through plasma- membrane-to-nucleus transfer of BRX. Because BRX is required for expression of certain auxin response factor targets, our data suggest a cell-type-specific endocytosis-dependent input into transcriptional auxin perception. This input sustains expression of a subset of auxin-responsive genes across the root meristem's division and transition zones and is essential for meristem growth. Thus, the endocytosis pattern provides specific positional information to modulate auxin response. RésuméLes plantes croissent continuellement tout au long de leur cycle de vie. Comprendre et expliquer les mécanismes impliqués dans ce phénomène reste à l'heure actuelle, un défi. L'hormone auxine a été identifiée comme une molécule essentielle à la régulation de la croissance des plantes, car impliquée dans la plupart des processus développementaux. Dans la racine, le transport polaire de l'auxine, par la génération d'un pic de concentration, spécifie et maintient la niche de cellules souches, et par la génération d'un gradient de concentration, contrôle la prolifération et la différentiation cellulaire. Puisque l'auxine est essentielle pour l'organisation et la maintenance du système racinaire, il est considéré comme son principal morphogène. Récemment, dans la plante modèle, Arabidopsis thalinana, un criblage des variations génétique a permis d'identifier le gène Brevis radix (BRX) comme facteur limitant l'expression des gènes de réponse à l'auxine et par là même, la croissance de la racine.Dans ce travail, nous avons découvert que BRX est une cible direct de l'auxine qui rétroactive positivement le signalement de l'hormone, agissant ainsi comme un régulateur transcriptionnel à travers l'interaction avec la protéine Monopteros (MP) de la famille des facteurs de réponse à l'auxine (Auxin Responsive Factor, ARF), et modulant ainsi la magnitude de la réponse des gènes reliés à l'auxine durant la division et la différentiation cellulaire dans le méristème de la racine. De plus, nous fournissons des preuves que BRX est activées au niveau de la membrane plasmique, tel une protéine associée se déplaçant à l'intérieur du noyau et modulant la croissance cellulaire.Pour mener à bien l'investigation des divergences entre la concentration de l'auxine et les schémas d'expression de ses propres gènes cibles, nous avons combiné les approches expérimentales et computationnelles. Les profiles d'expressions déviant du gradient d'auxine pourraient seulement être modéliser après intersection de l'activité de l'auxine avec les schémas différentiels d'endocytose observés et les boucles de rétroaction positives et autorégulatrices par le transfert de BRX de la membrane plasmique au noyau. Puisque BRX est requis pour l'expression de certains gènes cibles des facteurs de réponse à l'auxine, nos données suggèrent une contribution dépendante d'une endocytose spécifique au type de cellule dans la perception transcriptionnelle à l'auxine Cette contribution soutient l'expression d'un sous-set de gène de réponse à l'auxine dans la division du méristème racinaire et la zone de transition, et par conséquent, est essentielle pour la croissance méristematique. Ainsi, le schéma d'endocytose fournit des informations positionnelles spécifiques à la modulation de la réponse à l'auxine.

Role of CLOCK and circadian rhythms in calcium homeostasis

Le maintien d'une concentration sanguine constante de calcium est d'une importance cruciale et trois organes participent à la balance calcique normale : les reins, les intestins et les os. La concentration plasmatique de calcium est strictement régulée par l'hormone parathyroïdienne (PTH) et par la vitamine D. Des variations circadiennes de la PTH, de la vitamine D ainsi que du calcium plasmatique ont été décrites précédemment chez l'humain ainsi que chez le rat. Ces rythmes de PTH dans le sérum sont importants pour la régulation du remodelage de l'os. En effet, il a été montré chez les souris C57BL/6J que des injections de PTH une fois par jour mènent à une augmentation de la densité minérale de l'os alors que l'infusion en continu de PTH est associée à une diminution de cette densité. La vitamine D joue également un rôle fondamental dans la physiologie osseuse, car un déficit en vitamine D peut conduire à une ostéomalacie. Cependant la fonction des oscillations de vitamine D au niveau de l'homéostasie osseuse reste inconnue. L'horloge circadienne est un système interne de contrôle biologique du temps générant des rythmes de 24 heures dans l'expression des gènes, ainsi que dans la physiologie et le comportement. Ce contrôle s'opère par des boucles rétroactives positives et négatives de l'expression de gènes circadiens tels que CLOCK, BMAL1, CRY1 et 2 ou PERI et 2. Dans ce travail, nous avons émis l'hypothèse que l'homéostasie calcique est sous le contrôle de l'horloge circadienne. Dans un premier temps, nous avons montré chez les souris C57BL/6J des variations journalières des concentrations de calcium, de PTH et de vitamine D dans le sang, ainsi que de calcium dans les urines. Nous avons également démontré des changements au niveau de l'expression rénale des gènes importants dans l'homéostasie du calcium, tant au niveau de l'ARN messager que des protéines. Ensuite, pour analyser le rôle du système de l'horloge circadienne dans l'homéostasie du calcium, nous avons étudié des souris dans lesquelles a été supprimé le gène CLOCK crucial pour la fonction de l'horloge et nous avons comparé ces souris à des souris de type sauvage de même portée. Les souris CLOCK-I- étaient hypercalciuriques à chaque moment de la journée. Cependant le rythme circadien de l'excrétion de calcium était préservé. Le taux de calcium plasmatique ne différait pas entre les génotypes, mais les souris CLOCK -/- ne montraient pas de variations journalières de ce paramètre. Une perte du rythme journalier était également observée pour les niveaux de vitamine D, perte qui pourrait être une cause de l'altération de la micro-architecture osseuse révélée chez les souris CLOCK-/-. En effet, ces souris montrent une diminution du nombre de trabécules, de leur volume ainsi que de leur surface, ce qui suggère la présence d'ostéoporose. Nous avons également trouvé que le rythme de l'expression de l'ARN messager de CYP27B1 était aboli dans les reins des souris CLOCK -/-, ce qui peut expliquer l'altération du rythme de la vitamine D. Les taux sanguins de PTH étaient comparables entre les souris CLOCK -/- et de type sauvage. Dans les reins, une augmentation de l'expression de l'ARN messager de TRPV5 et NCX1 a été constatée, ce qui suggérerait une augmentation de la réabsorption de calcium dans le tubule convoluté distal et dans le tubule connecteur. Dans les intestins, la réabsorption calcique était diminuée, chez les souris CLOCK-I-, fait confirmé par une diminution des niveaux d'ARN messager de TRPV6 et PMCAL. En résumé, la suppression du gène CLOCK chez les souris a conduit à une hypercalciurie, une altération du rythme des taux plasmatiques de calcium et de vitamine D et à une détérioration de l'architecture osseuse. Pour conclure, ces résultats montrent que l'horloge circadienne est essentielle à l'homéostasie calcique ainsi qu'à la physiologie des os. - L'ostéoporose affecte environ 22 millions de femmes et 5.5 millions d'hommes en Europe, réduisant significativement leur qualité de vie et a causé 3.5 millions de nouvelles fractures en 2010. Les dépenses totales liées à ces fractures ont atteint 37 milliards d'euro et ce coût devrait augmenter de 25% d'ici à 2025. Le nombre de nouvelles fractures dues à l'ostéoporose à travers le monde est estimé à environ 1000 par heure. Parmi les causes de l'ostéoporose, le déficit én calcium et/ou en vitamine D joue un rôle important, mais il existe également des causes génétiques ou liées à des facteurs comme les hormones sexuelles (estrogènes, testostérone), l'âge, le tabac, le poids corporel, certains médicaments,... La vie est rythmique : ceci est dû à l'alternance naturelle du jour et de la nuit et de ses effets sur le corps. La prise alimentaire, par exemple, est un processus qui a lieu pendant la phase active, qui est prévisible (il se produit toujours au même moment) et qui peut être anticipé par le corps. Pour cela, une horloge interne est présente dans chaque cellule du corps et est synchronisée par la lumière du jour, entre autres stimuli. Cette horloge indique la phase du jour et régule l'expression de gènes impliqués dans les différents processus qui nécessitent une anticipation. Pendant mon travail de thèse, je me suis demandé si des îythmes circadiens (c'est-à-dire d'une durée d'environ 24 heures et indépendants des stimuli externes) étaient observables'pour les gènes régulant les flux de calcium dans le corps et si l'interruption de ces rythmes pouvait mener à des altérations de la qualité de l'os. J'ai d'abord travaillé avec des souris normales et j'ai pu montrer la présence de rythmes au niveau du calcium sanguin et urinaire, mais également au niveau des hormones et gènes qui contrôlent le métabolisme du calcium dans le corps, comme la vitamine D et l'hormone parathyroidienne. De manière intéressante, j'ai observé que la plupart de ces gènes ont un rythme synchronisé. J'ai ensuite utilisé un modèle de souris dans lequel l'horloge interne a été génétiquement invalidée et j'ai montré que ces souris présentent une augmentation de leur excrétion urinaire de calcium et un rythme circadien altéré de la vitamine D dans le sang. Ces souris absorbent aussi moins bien le calcium intestinal et présentent une ostéoporose marquée. Ce travail montre donc que l'horloge interne est nécessaire pour établir un rythme circadiens de certains facteurs influant les flux de calcium dans l'organisme, comme la vitamine D, et que la perturbation de ces rythmes mène à une dérégulation du métabolisme osseux. Ainsi, la perturbation de l'horloge interne peut causer une ostéoporose et une hypercalciurie qui pourraient aboutir à la formation de fractures et de calculs rénaux. L'extrapolation de ces observations chez l'homme ou à des changements plus subtiles des rythmes circadiens, comme le décalage horaire, restent à montrer. Cette recherche a démontré que les rythmes circadiens des mécanismes de régulation des flux de calcium dans l'organisme sont essentiels au maintien d'un squelette normal et suggère que les perturbations des rythmes circadiens pourraient être une nouvelle cause de l'ostéoporose. - Maintaining constant calcium concentration in the plasma is of a crucial importance and three organs participate in normal calcium balance - kidney, gut and bone. Plasma calcium concentration is strictly regulated by parathyroid hormone (PTH) and vitamin D. Circadian variations of PTH, vitamin D and plasma calcium were previously described in humans, as well as in rats. Rhythms in serum PTH are important for balanced bone remodelling. Indeed in C57BL/6J mice, PTH injection once per day leads to an increase in bone mineral density (BMD), whilst continuous infusion is associated with decreased BMD. Vitamin D also plays a crucial role in bone physiology, since the deficiency in vitamin D can lead to rickets/osteomalacia. However, the role of vitamin D rhythms in bone homeostasis remains unknown. The circadian clock is an. internal time-keeping system generating rhythms in gene expression with 24h periodicity, in physiology and in behaviour. It is operated by positive- and negative-feedback loops of circadian genes, such as CLOCK, BMAL1, CRY1 and 2 or PERI and 2. In this work, we hypothesized, that calcium homeostasis is under the control of the circadian clock. First, we showed daily variations in urinary calcium and serum calcium, PTH and l,25(OH)2 vitamin D, together with renal mRNA and protein levels of genes involved in calcium homeostasis in C57BL/6J mice. Second, and to investigate the role of the circadian clock system in calcium handling, we studied mice lacking the gene CLOCK crucial for fonction of the clock system and compared them to the WT littermates. CLOCK-/- mice were hypercalciuric at all timepoints of the day. However, the circadian rhythm of calcium excretion was preserved. Serum calcium levels did not differ between the genotypes, but CLOCK-/- mice did not exhibit daily variation for this parameter. Loss of rhythm was observed also for serum l,25(OH)2 vitamin D levels, which may be one of the causes of altered bone microarchitecture that was revealed in CLOCK-/- mice. They displayed increased trabecular separation and decreased trabecular number, trabecular bone volume and trabecular bone surface, suggestive of osteoporosis. We found that the rhythm of the mRNA expression of CYP27B1 was abolished in the kidney of CLOCK-/- mice, which could induce the altered rhythm of l,25(OH)2 vitamin. Serum PTH levels were comparable between CLOCK-/- and WT mice. In the kidney, increased mRNA expression of TRPV5 and NCX1 suggests increased calcium reabsorption in the distal convoluted and connecting tubule. In the gut, intestinal calcium absorption was decreased in CLOCK¬/- mice, confirmed by decreased mRNA levels of TRPV6 and PMCA1. In summary, deletion of the CLOCK gene in mice conducts to hypercalciuria, alteration of the rhythm in serum calcium and l,25(OH)2D levels, and impainnent of their bone microarchitecture. In conclusion, these data show that the circadian clock system is essential in calcium homeostasis and bone physiology.

Simulation of stromatolite growth using diffusion-limited aggregation and a discrete model of biogeochemical exchanges in microbial mat

Resume : Mieux comprendre les stromatolithes et les tapis microbiens est un sujet important en biogéosciences puisque cela aide à l'étude des premières formes de vie sur Terre, a mieux cerner l'écologie des communautés microbiennes et la contribution des microorganismes a la biominéralisation, et même à poser certains fondements dans les recherches en exobiologie. D'autre part, la modélisation est un outil puissant utilisé dans les sciences naturelles pour appréhender différents phénomènes de façon théorique. Les modèles sont généralement construits sur un système d'équations différentielles et les résultats sont obtenus en résolvant ce système. Les logiciels disponibles pour implémenter les modèles incluent les logiciels mathématiques et les logiciels généraux de simulation. L'objectif principal de cette thèse est de développer des modèles et des logiciels pour aider a comprendre, via la simulation, le fonctionnement des stromatolithes et des tapis microbiens. Ces logiciels ont été développés en C++ en ne partant d'aucun pré-requis de façon a privilégier performance et flexibilité maximales. Cette démarche permet de construire des modèles bien plus spécifiques et plus appropriés aux phénomènes a modéliser. Premièrement, nous avons étudié la croissance et la morphologie des stromatolithes. Nous avons construit un modèle tridimensionnel fondé sur l'agrégation par diffusion limitée. Le modèle a été implémenté en deux applications C++: un moteur de simulation capable d'exécuter un batch de simulations et de produire des fichiers de résultats, et un outil de visualisation qui permet d'analyser les résultats en trois dimensions. Après avoir vérifié que ce modèle peut en effet reproduire la croissance et la morphologie de plusieurs types de stromatolithes, nous avons introduit un processus de sédimentation comme facteur externe. Ceci nous a mené a des résultats intéressants, et permis de soutenir l'hypothèse que la morphologie des stromatolithes pourrait être le résultat de facteurs externes autant que de facteurs internes. Ceci est important car la classification des stromatolithes est généralement fondée sur leur morphologie, imposant que la forme d'un stromatolithe est dépendante de facteurs internes uniquement (c'est-à-dire les tapis microbiens). Les résultats avancés dans ce mémoire contredisent donc ces assertions communément admises. Ensuite, nous avons décidé de mener des recherches plus en profondeur sur les aspects fonctionnels des tapis microbiens. Nous avons construit un modèle bidimensionnel de réaction-diffusion fondé sur la simulation discrète. Ce modèle a été implémenté dans une application C++ qui permet de paramétrer et exécuter des simulations. Nous avons ensuite pu comparer les résultats de simulation avec des données du monde réel et vérifier que le modèle peut en effet imiter le comportement de certains tapis microbiens. Ainsi, nous avons pu émettre et vérifier des hypothèses sur le fonctionnement de certains tapis microbiens pour nous aider à mieux en comprendre certains aspects, comme la dynamique des éléments, en particulier le soufre et l'oxygène. En conclusion, ce travail a abouti à l'écriture de logiciels dédiés à la simulation de tapis microbiens d'un point de vue tant morphologique que fonctionnel, suivant deux approches différentes, l'une holistique, l'autre plus analytique. Ces logiciels sont gratuits et diffusés sous licence GPL (General Public License). Abstract : Better understanding of stromatolites and microbial mats is an important topic in biogeosciences as it helps studying the early forms of life on Earth, provides clues re- garding the ecology of microbial ecosystems and their contribution to biomineralization, and gives basis to a new science, exobiology. On the other hand, modelling is a powerful tool used in natural sciences for the theoretical approach of various phenomena. Models are usually built on a system of differential equations and results are obtained by solving that system. Available software to implement models includes mathematical solvers and general simulation software. The main objective of this thesis is to develop models and software able to help to understand the functioning of stromatolites and microbial mats. Software was developed in C++ from scratch for maximum performance and flexibility. This allows to build models much more specific to a phenomenon rather than general software. First, we studied stromatolite growth and morphology. We built a three-dimensional model based on diffusion-limited aggregation. The model was implemented in two C++ applications: a simulator engine, which can run a batch of simulations and produce result files, and a Visualization tool, which allows results to be analysed in three dimensions. After verifying that our model can indeed reproduce the growth and morphology of several types of stromatolites, we introduced a sedimentation process as an external factor. This lead to interesting results, and allowed to emit the hypothesis that stromatolite morphology may be the result of external factors as much as internal factors. This is important as stromatolite classification is usually based on their morphology, imposing that a stromatolite shape is dependant on internal factors only (i.e. the microbial mat). This statement is contradicted by our findings, Second, we decided to investigate deeper the functioning of microbial mats, We built a two-dimensional reaction-diffusion model based on discrete simulation, The model was implemented in a C++ application that allows setting and running simulations. We could then compare simulation results with real world data and verify that our model can indeed mimic the behaviour of some microbial mats. Thus, we have proposed and verified hypotheses regarding microbial mats functioning in order to help to better understand them, e.g. the cycle of some elements such as oxygen or sulfur. ln conclusion, this PhD provides a simulation software, dealing with two different approaches. This software is free and available under a GPL licence.

Estimation of glomerular filtration rate in hospitalized patients : are we overestimating renal function?

Estimer la filtration glomérulaire chez les personnes âgées, tout en tenant compte de la difficulté supplémentaire d'évaluer leur masse musculaire, est difficile et particulièrement important pour la prescription de médicaments. Le taux plasmatique de la creatinine dépend à la fois de la fraction d'élimination rénale et extra-rénale et de la masse musculaire. Actuellement, pour estimer là filtration glomérulaire différentes formules sont utilisées, qui se fondent principalement sur la valeur de la créatinine. Néanmoins, en raison de la fraction éliminée par les voies tubulaires et intestinales la clairance de la créatinine surestime généralement le taux de filtration glomérulaire (GFR). Le but de cette étude est de vérifier la fiabilité de certains marqueurs et algorithmes de la fonction rénale actuellement utilisés et d'évaluer l'avantage additionnel de prendre en considération la masse musculaire mesurée par la bio-impédance dans une population âgée (> 70 ans) et avec une fonction rénale chronique compromise basée sur MDRD eGFR (CKD stades lll-IV). Dans cette étude, nous comparons 5 équations développées pour estimer la fonction rénale et basées respectivement sur la créatinine sérique (Cockcroft et MDRD), la cystatine C (Larsson), la créatinine combinée à la bêta-trace protéine (White), et la créatinine ajustée à la masse musculaire obtenue par analyse de la bio-impédance (MacDonald). La bio-impédance est une méthode couramment utilisée pour estimer la composition corporelle basée sur l'étude des propriétés électriques passives et de la géométrie des tissus biologiques. Cela permet d'estimer les volumes relatifs des différents tissus ou des fluides dans le corps, comme par exemple l'eau corporelle totale, la masse musculaire (=masse maigre) et la masse grasse corporelle. Nous avons évalué, dans une population âgée d'un service interne, et en utilisant la clairance de l'inuline (single shot) comme le « gold standard », les algorithmes de Cockcroft (GFR CKC), MDRD, Larsson (cystatine C, GFR CYS), White (beta trace protein, GFR BTP) et Macdonald (GFR = ALM, la masse musculaire par bio-impédance. Les résultats ont montré que le GFR (mean ± SD) mesurée avec l'inuline et calculée avec les algorithmes étaient respectivement de : 34.9±20 ml/min pour l'inuline, 46.7±18.5 ml/min pour CKC, 47.2±23 ml/min pour CYS, 54.4±18.2ml/min pour BTP, 49±15.9 ml/min pour MDRD et 32.9±27.2ml/min pour ALM. Les courbes ROC comparant la sensibilité et la spécificité, l'aire sous la courbe (AUC) et l'intervalle de confiance 95% étaient respectivement de : CKC 0 68 (055-0 81) MDRD 0.76 (0.64-0.87), Cystatin C 0.82 (0.72-0.92), BTP 0.75 (0.63-0.87), ALM 0.65 (0.52-0.78). ' En conclusion, les algorithmes comparés dans cette étude surestiment la GFR dans la population agee et hospitalisée, avec des polymorbidités et une classe CKD lll-IV. L'utilisation de l'impédance bioelectrique pour réduire l'erreur de l'estimation du GFR basé sur la créatinine n'a fourni aucune contribution significative, au contraire, elle a montré de moins bons résultats en comparaison aux autres equations. En fait dans cette étude 75% des patients ont changé leur classification CKD avec MacDonald (créatinine et masse musculaire), contre 49% avec CYS (cystatine C), 56% avec MDRD,52% avec Cockcroft et 65% avec BTP. Les meilleurs résultats ont été obtenus avec Larsson (CYS C) et la formule de Cockcroft.

The Mont Collon mafic complex (Austroalpine Dent Blanche nappe) : permian evolution of the Western European mantle

Résumé: Le complexe du Mont Collon (nappe de la Dent Blanche, Austroalpin) est l'un des exemples les mieux préservés du magmatisme mafique permien des Alpes occidentales. Il est composé d'affleurements discontinus et d'une stratification magmatique en son centre (Dents de Bertol) et est composé à 95% de roches mafiques cumulatives (gabbros à olivine et/ou cpx, anorthositiques, troctolites, wehrlites et wehrlites à plagioclase) et localement de quelques gabbros pegmatitiques. Ces faciès sont recoupés par de nombreux filons acides (aphtes, pegmatites quartziques, microgranodiorites et filons anorthositiques) et mafiques tardifs (dikes mélanocrates riches en Fe et Ti). Les calculs thermométriques (équilibre olivine-augite) montrent des températures de 1070-1120 ± 6°C, tandis que le thermomètre amphibole-plagioclase indique une température de 740 ± 40°C à 0.5 GPa pour les amphiboles magmatiques tardives. La geobarométrie sur pyroxène donne des pressions moyennes de 0.3-0.6 GPa, indiquant un emplacement dans la croûte moyenne. De plus, les températures obtenues sur des amphiboles coronitiques indiquent des températures de l'ordre de 700 ± 40°C confirmant que les réactions coronitiques apparaissent dans des conditions subsolidus. Les âges concordants U/Pb sur zircons de 284.2 ± 0.6 et 282.9 ± 0.6 Ma obtenus sur un gabbro pegmatitique et une pegmatitique quartzique, sont interprétés comme des âges de cristallisation. Les datations 40Ar/39Ar sur amphiboles des filons mélanocrates donnent un âge plateau de 260.2 ± 0.7 Ma, qui est probablement très proche de l'âge de cristallisation. Ainsi, cet age 40Ar/39Ar indique un second évènement magmatique au sein du complexe. Les compositions des roches totales en éléments majeurs et traces montrent peu de variations, ainsi que le Mg# (75-80). Les éléments traces enregistrent le caractère cumulatif des roches (anomalie positive en Eu) et révèlent des anomalies négatives systématiques en Nb, Ta, Zr, Hf et Ti dans les faciès basiques. Le manque de corrélation entre éléments majeurs et traces est caractéristique d'un processus de cristallisation in situ impliquant une quantité variable de liquide interstitiel (L) entre les phases cumulus. Les distributions des éléments traces dans les minéraux sont homogènes, indiquant une rééquilibration .subsolidus entre cristaux et liquide interstitiel. Un modèle quantitatif basé sur les équations de cristallisation in situ de Langmuir reproduisent correctement les concentrations en terres rares légères des minéraux cumulatifs montrant la présence de 0 à 35% de liquide interstitiel L pour des degrés de différenciation F de 0 à 45%, par rapport au faciès les moins évolués du complexe. En outre, les valeurs de L sont bien corrélées avec les proportions modales d'amphibole interstitielle et les concentrations en éléments incompatibles des roches (Zr, Nb). Le liquide parental calculé des cumulats du Mont Collon est caractérisé par un enrichissement relatif en terres rares légères et Th, un appauvrissement en terres rares lourdes typique d'une affinité transitionnelle (T-MORB) et une forte anomalie négative en Nb-Ta. Les roches cumulatives montrent des compositions isotopiques en Nd-Sr proches de la terre globale silicatée (BSE), soit 0.6<εNdi<+3.2, 0.7045<87Sr/86Sri<0.7056. Les rapports initiaux en Pb indiquent une source dans le manteau enrichi subcontinental lithosphérique, préalablement contaminé par des sédiments océaniques. Les dikes mélanocrates Fe-Ti sont représentatifs de liquides et ont des spectres de terres rares enrichis, une anomalie positive en Nb-Ta et des εNdi de +7, des 87Sr/86Sri de 0.703 et des rapports initiaux en Pb, similaires à ceux des basaltes d'île océanique, indiquant une source asthénosphérique modérément appauvrie. Ainsi, la fusion partielle du manteau lithosphérique subcontinental est induite par l'amincissement post-orogénique et la remontée de l'asthénosphère. Les filons mélanocrates proviennent, après délamination du manteau lithosphérique, de la fusion de l'asthénosphère. Abstract The early Permian Mont Collon mafic complex (Dent Blanche nappe, Austroalpine nappe system) is one of the best preserved examples of the Permian mafic magmatism in the Western Alps. It is composed of discontinuous exposures and a well-preserved magmatic layering (the Dents de Bertol cliff) crops out in the center part of the complex. It mainly consists of cumulative mafic rocks, which represent 95 vol-% of the mafic complex (ol- and cpx-bearing gabbros and rare anorthositic layers, troctolites, wehrlites and plagioclase-wehrlites) and locally pegmatitic gabbros. All these facies are crosscut by widespread acidic (aplites, quartz-rich pegmatites, microgranodiorites) and late mafic Fe-Ti melanocratic dikes. Olivine-augite thermometric calculations yield a range of 1070-1120 ± 6°C, while amphibole-plagioclase thermometer yields a temperature of 740 ± 40°C at 0.5 GPa. Pyroxene geobarometry points to a pressure of 0.3-0.6 GPa, indicating a middle crustal level of emplacement. Moreover, temperature calculations on the Mont Conon coronitic amphiboles indicate temperatures of 700 ± 40°C, close to those calculated for magmatic amphiboles. These temperatures confirm that coronitic reactions occurred at subsolidus conditions. ID-TIMS U/Pb zircon ages of 284.2 ± 0.6 and 282.9 ± 0.6 Ma obtained on a pegmatitic gabbro and a quartz-pegmatitic dike, respectively, were interpreted as the crystallization ages of these rocks. 40Ar/39Ar dating on amphiboles from Fe-Ti melanocratic dikes yields a plateau age of 260.2 ± 0.7 Ma, which is probably very close to the crystallization age. Consequently, this 40Ar/P39Ar age indicates a second magmatic event. Whole-rock major- and trace-element compositions show little variation across the whole intrusion and Mg-number stays within a narrow range (75-80). Trace-element concentrations record the cumulative nature of the rocks (e.g. positive Eu anomaly) and reveal systematic Nb, Ta, Zr, Hf and Ti negative anomalies for all basic facies. The lack of correlation between major and trace elements is characteristic of an in situ crystallization process involving variable amounts of interstitial liquid (L) trapped between the cumulus mineral phases. LA-ICPMS measurements show that trace-element distributions in minerals are homogeneous, pointing to subsolidus re-equilibration between crystals and interstitial melts. A quantitative modeling based on Langmuir's in situ crystallization equation successfully reproduced the Rare Earth Element (REE) concentrations in cumulitic minerals. The calculated amounts of interstitial liquid L vary between 0 and 35% for degrees of differentiation F of 0 to 45%, relative to the least evolved facies of the intrusion. Furthermore, L values are well correlated with the modal proportions of interstitial amphibole and whole-rock incompatible trace-element concentrations (e.g. Zr, Nb) of the tested samples. The calculated parental melt of the Mont Collon cumulates is characterized by a relative enrichment in Light REE and Th, a depletion in Heavy REE, typical of a transitional affinity (T-MORB), and strong negative Nb-Ta anomaly. Cumulative rocks display Nd-Sr isotopic compositions close to the BSE (-0.6 < εNdi < +3.2, 0.7045 < 87Sr/86Sri < 0.7056). Initial Pb ratios point to an origin from the melting of an enriched subcontinental lithospheric mantle source, previously contaminated at the source by oceanic sediments. The contrasted alkaline Fe-Ti melanocratic dikes are representative of liquids. They display enriched fractionated REE patterns, a positive Nb-Ta anomaly and εNdi of +7, 87Sr/86Sri of 0.703 and initial Pb ratios, all reminiscent of Ocean Island Basalt-type rocks, pointing to a moderately

Etude de caractères manuscrits : de la caractérisation morphologique à l'individualisation du scripteur

PLAN DE LA RECHERCHE La présentation de ces travaux de recherche est divisée en neuf chapitres : Chapitre 1. Les éléments qui concourent d'un point de vue essentiellement théorique à l'individualité de l'écriture manuscrite sont décrits. Dans un premier temps sont abordés les aspects anatomo-fonctionnels de l'écriture, permettant de saisir quelles sont les contributions musculaires et neuronales de l'individualité de l'écriture. Ensuite sont présentées les caractéristiques de l'écriture telles qu'elles sont étudiées et classifiées par les praticiens, ainsi que la notion de leur variation naturelle. Enfin, les "lois" dites fondamentales de l'écriture sont mentionnées et discutées. Chapitre 2. Les problèmes soulevés lors de débats contradictoires menés au sein de la justice nord-américaine, qui ont porté sur l'admissibilité des expertises d'écritures manuscrites, sont exposés. Ces débats ont mis en évidence la nécessité d'entreprendre des recherches scientifiques approfondies afin d'obtenir des données sur la variabilité des écritures manuscrites. Chapitre 3. Après une présentation de la nécessité de la récolte et de l'analyse de données quantitatives caractérisant l'écriture manuscrite, diverses études ayant cherché à réduire la part de subjectivité au sein du processus d'analyse des écritures sont résumées. Les études se basant sur des mesures ou des classifications effectuées visuellement sont traitées d'abord, suivies des études visant à estimer la fréquence d'apparition de certaines caractéristiques dans certaines populations de scripteurs, puis des études plus récentes taisant intervenir des méthodes de vérification et d'identification de scripteurs. Chapitre 4. Le choix de soumettre la forme des boucles des caractères manuscrits à une étude statistique est justifié, tant par l'importance présumée de son polymorphisme que de sa fréquente utilisation par les praticiens du domaine d'intérêt. Chapitre 5. Les descripteurs de Fourier, utilisés pour caractériser la forme du contour des boucles manuscrites, sont introduits. La mise en évidence de la diversité de leurs domaines d'application précède leur description mathématique, en mettant notamment en avant les avantages de cette technique parmi d'autres méthodes de caractérisation de la Tonne de contours. Chapitre 6. La première partie de la recherche expérimentale est consacrée à la validation de la méthodologie. L'objectif de cette étape consiste à montrer si la procédure d'analyse d'images et les descripteurs de Fourier permettent de démontrer des différences entre des groupes de boucles dont la Tonne est visuellement distincte. Chapitre 7. La variabilité des paramètres de la forme de boucles mesurés au sein d'une population plus importante de scripteurs est étudiée. Des tendances générales de forme permettant de grouper .des scripteurs sont soulignées, puis les caractéristiques particulières pouvant permettre d'individualiser un scripteur sont signalées. Chapitre 8. L'influence de la taille sur la forme des boucles de caractères manuscrits est évaluée, afin de déterminer si des boucles de la même lettre et d'un même scripteur mais de taille différente présentent des aspects de forme distincts. Des tendances générales de modification de la forme ainsi que certaines singularités sont mises en exergue. Chapitre 9. Des modèles d'évaluation permettant de quantifier la valeur probante de l'indice que représente la forme des boucles manuscrites sont exposés. Leur application, qui fait intervenir le rapport de vraisemblance, confirme le potentiel d'individualisation du scripteur mis en évidence précédemment. D'autre part, l'utilisation du rapport de vraisemblance comme métrique confirme les résultats obtenus par diverses méthodes statistiques dans les chapitres précédents. Enfin, la discussion générale résume les apports principaux de cette recherche et en dégage les conséquences pour le praticien, puis propose des pistes possibles pour la poursuite de la recherche.