Modelling tree population dynamics at the alpine and boreal tree-line ecotones in response to climate and land-use change

Summary Ecotones are sensitive to change because they contain high numbers of species living at the margin of their environmental tolerance. This is equally true of tree-lines, which are determined by attitudinal or latitudinal temperature gradients. In the current context of climate change, they are expected to undergo modifications in position, tree biomass and possibly species composition. Attitudinal and latitudinal tree-lines differ mainly in the steepness of the underlying temperature gradient: distances are larger at latitudinal tree-lines, which could have an impact on the ability of tree species to migrate in response to climate change. Aside from temperature, tree-lines are also affected on a more local level by pressure from human activities. These are also changing as a consequence of modifications in our societies and may interact with the effects of climate change. Forest dynamics models are often used for climate change simulations because of their mechanistic processes. The spatially-explicit model TreeMig was used as a base to develop a model specifically tuned for the northern European and Alpine tree-line ecotones. For the latter, a module for land-use change processes was also added. The temperature response parameters for the species in the model were first calibrated by means of tree-ring data from various species and sites at both tree-lines. This improved the growth response function in the model, but also lead to the conclusion that regeneration is probably more important than growth for controlling tree-line position and species' distributions. The second step was to implement the module for abandonment of agricultural land in the Alps, based on an existing spatial statistical model. The sensitivity of its most important variables was tested and the model's performance compared to other modelling approaches. The probability that agricultural land would be abandoned was strongly influenced by the distance from the nearest forest and the slope, bath of which are proxies for cultivation costs. When applied to a case study area, the resulting model, named TreeMig-LAb, gave the most realistic results. These were consistent with observed consequences of land-abandonment such as the expansion of the existing forest and closing up of gaps. This new model was then applied in two case study areas, one in the Swiss Alps and one in Finnish Lapland, under a variety of climate change scenarios. These were based on forecasts of temperature change over the next century by the IPCC and the HadCM3 climate model (ΔT: +1.3, +3.5 and +5.6 °C) and included a post-change stabilisation period of 300 years. The results showed radical disruptions at both tree-lines. With the most conservative climate change scenario, species' distributions simply shifted, but it took several centuries reach a new equilibrium. With the more extreme scenarios, some species disappeared from our study areas (e.g. Pinus cembra in the Alps) or dwindled to very low numbers, as they ran out of land into which they could migrate. The most striking result was the lag in the response of most species, independently from the climate change scenario or tree-line type considered. Finally, a statistical model of the effect of reindeer (Rangifer tarandus) browsing on the growth of Pinus sylvestris was developed, as a first step towards implementing human impacts at the boreal tree-line. The expected effect was an indirect one, as reindeer deplete the ground lichen cover, thought to protect the trees against adverse climate conditions. The model showed a small but significant effect of browsing, but as the link with the underlying climate variables was unclear and the model was not spatial, it was not usable as such. Developing the TreeMig-LAb model allowed to: a) establish a method for deriving species' parameters for the growth equation from tree-rings, b) highlight the importance of regeneration in determining tree-line position and species' distributions and c) improve the integration of social sciences into landscape modelling. Applying the model at the Alpine and northern European tree-lines under different climate change scenarios showed that with most forecasted levels of temperature increase, tree-lines would suffer major disruptions, with shifts in distributions and potential extinction of some tree-line species. However, these responses showed strong lags, so these effects would not become apparent before decades and could take centuries to stabilise. Résumé Les écotones son sensibles au changement en raison du nombre élevé d'espèces qui y vivent à la limite de leur tolérance environnementale. Ceci s'applique également aux limites des arbres définies par les gradients de température altitudinaux et latitudinaux. Dans le contexte actuel de changement climatique, on s'attend à ce qu'elles subissent des modifications de leur position, de la biomasse des arbres et éventuellement des essences qui les composent. Les limites altitudinales et latitudinales diffèrent essentiellement au niveau de la pente des gradients de température qui les sous-tendent les distance sont plus grandes pour les limites latitudinales, ce qui pourrait avoir un impact sur la capacité des espèces à migrer en réponse au changement climatique. En sus de la température, la limite des arbres est aussi influencée à un niveau plus local par les pressions dues aux activités humaines. Celles-ci sont aussi en mutation suite aux changements dans nos sociétés et peuvent interagir avec les effets du changement climatique. Les modèles de dynamique forestière sont souvent utilisés pour simuler les effets du changement climatique, car ils sont basés sur la modélisation de processus. Le modèle spatialement explicite TreeMig a été utilisé comme base pour développer un modèle spécialement adapté pour la limite des arbres en Europe du Nord et dans les Alpes. Pour cette dernière, un module servant à simuler des changements d'utilisation du sol a également été ajouté. Tout d'abord, les paramètres de la courbe de réponse à la température pour les espèces inclues dans le modèle ont été calibrées au moyen de données dendrochronologiques pour diverses espèces et divers sites des deux écotones. Ceci a permis d'améliorer la courbe de croissance du modèle, mais a également permis de conclure que la régénération est probablement plus déterminante que la croissance en ce qui concerne la position de la limite des arbres et la distribution des espèces. La seconde étape consistait à implémenter le module d'abandon du terrain agricole dans les Alpes, basé sur un modèle statistique spatial existant. La sensibilité des variables les plus importantes du modèle a été testée et la performance de ce dernier comparée à d'autres approches de modélisation. La probabilité qu'un terrain soit abandonné était fortement influencée par la distance à la forêt la plus proche et par la pente, qui sont tous deux des substituts pour les coûts liés à la mise en culture. Lors de l'application en situation réelle, le nouveau modèle, baptisé TreeMig-LAb, a donné les résultats les plus réalistes. Ceux-ci étaient comparables aux conséquences déjà observées de l'abandon de terrains agricoles, telles que l'expansion des forêts existantes et la fermeture des clairières. Ce nouveau modèle a ensuite été mis en application dans deux zones d'étude, l'une dans les Alpes suisses et l'autre en Laponie finlandaise, avec divers scénarios de changement climatique. Ces derniers étaient basés sur les prévisions de changement de température pour le siècle prochain établies par l'IPCC et le modèle climatique HadCM3 (ΔT: +1.3, +3.5 et +5.6 °C) et comprenaient une période de stabilisation post-changement climatique de 300 ans. Les résultats ont montré des perturbations majeures dans les deux types de limites de arbres. Avec le scénario de changement climatique le moins extrême, les distributions respectives des espèces ont subi un simple glissement, mais il a fallu plusieurs siècles pour qu'elles atteignent un nouvel équilibre. Avec les autres scénarios, certaines espèces ont disparu de la zone d'étude (p. ex. Pinus cembra dans les Alpes) ou ont vu leur population diminuer parce qu'il n'y avait plus assez de terrains disponibles dans lesquels elles puissent migrer. Le résultat le plus frappant a été le temps de latence dans la réponse de la plupart des espèces, indépendamment du scénario de changement climatique utilisé ou du type de limite des arbres. Finalement, un modèle statistique de l'effet de l'abroutissement par les rennes (Rangifer tarandus) sur la croissance de Pinus sylvestris a été développé, comme première étape en vue de l'implémentation des impacts humains sur la limite boréale des arbres. L'effet attendu était indirect, puisque les rennes réduisent la couverture de lichen sur le sol, dont on attend un effet protecteur contre les rigueurs climatiques. Le modèle a mis en évidence un effet modeste mais significatif, mais étant donné que le lien avec les variables climatiques sous jacentes était peu clair et que le modèle n'était pas appliqué dans l'espace, il n'était pas utilisable tel quel. Le développement du modèle TreeMig-LAb a permis : a) d'établir une méthode pour déduire les paramètres spécifiques de l'équation de croissance ä partir de données dendrochronologiques, b) de mettre en évidence l'importance de la régénération dans la position de la limite des arbres et la distribution des espèces et c) d'améliorer l'intégration des sciences sociales dans les modèles de paysage. L'application du modèle aux limites alpines et nord-européennes des arbres sous différents scénarios de changement climatique a montré qu'avec la plupart des niveaux d'augmentation de température prévus, la limite des arbres subirait des perturbations majeures, avec des glissements d'aires de répartition et l'extinction potentielle de certaines espèces. Cependant, ces réponses ont montré des temps de latence importants, si bien que ces effets ne seraient pas visibles avant des décennies et pourraient mettre plusieurs siècles à se stabiliser.

KCM, a general framework for codon-based model of evolution

With the advancement of high-throughput sequencing and dramatic increase of available genetic data, statistical modeling has become an essential part in the field of molecular evolution. Statistical modeling results in many interesting discoveries in the field, from detection of highly conserved or diverse regions in a genome to phylogenetic inference of species evolutionary history Among different types of genome sequences, protein coding regions are particularly interesting due to their impact on proteins. The building blocks of proteins, i.e. amino acids, are coded by triples of nucleotides, known as codons. Accordingly, studying the evolution of codons leads to fundamental understanding of how proteins function and evolve. The current codon models can be classified into three principal groups: mechanistic codon models, empirical codon models and hybrid ones. The mechanistic models grasp particular attention due to clarity of their underlying biological assumptions and parameters. However, they suffer from simplified assumptions that are required to overcome the burden of computational complexity. The main assumptions applied to the current mechanistic codon models are (a) double and triple substitutions of nucleotides within codons are negligible, (b) there is no mutation variation among nucleotides of a single codon and (c) assuming HKY nucleotide model is sufficient to capture essence of transition- transversion rates at nucleotide level. In this thesis, I develop a framework of mechanistic codon models, named KCM-based model family framework, based on holding or relaxing the mentioned assumptions. Accordingly, eight different models are proposed from eight combinations of holding or relaxing the assumptions from the simplest one that holds all the assumptions to the most general one that relaxes all of them. The models derived from the proposed framework allow me to investigate the biological plausibility of the three simplified assumptions on real data sets as well as finding the best model that is aligned with the underlying characteristics of the data sets. -- Avec l'avancement de séquençage à haut débit et l'augmentation dramatique des données géné¬tiques disponibles, la modélisation statistique est devenue un élément essentiel dans le domaine dé l'évolution moléculaire. Les résultats de la modélisation statistique dans de nombreuses découvertes intéressantes dans le domaine de la détection, de régions hautement conservées ou diverses dans un génome de l'inférence phylogénétique des espèces histoire évolutive. Parmi les différents types de séquences du génome, les régions codantes de protéines sont particulièrement intéressants en raison de leur impact sur les protéines. Les blocs de construction des protéines, à savoir les acides aminés, sont codés par des triplets de nucléotides, appelés codons. Par conséquent, l'étude de l'évolution des codons mène à la compréhension fondamentale de la façon dont les protéines fonctionnent et évoluent. Les modèles de codons actuels peuvent être classés en trois groupes principaux : les modèles de codons mécanistes, les modèles de codons empiriques et les hybrides. Les modèles mécanistes saisir une attention particulière en raison de la clarté de leurs hypothèses et les paramètres biologiques sous-jacents. Cependant, ils souffrent d'hypothèses simplificatrices qui permettent de surmonter le fardeau de la complexité des calculs. Les principales hypothèses retenues pour les modèles actuels de codons mécanistes sont : a) substitutions doubles et triples de nucleotides dans les codons sont négligeables, b) il n'y a pas de variation de la mutation chez les nucléotides d'un codon unique, et c) en supposant modèle nucléotidique HKY est suffisant pour capturer l'essence de taux de transition transversion au niveau nucléotidique. Dans cette thèse, je poursuis deux objectifs principaux. Le premier objectif est de développer un cadre de modèles de codons mécanistes, nommé cadre KCM-based model family, sur la base de la détention ou de l'assouplissement des hypothèses mentionnées. En conséquence, huit modèles différents sont proposés à partir de huit combinaisons de la détention ou l'assouplissement des hypothèses de la plus simple qui détient toutes les hypothèses à la plus générale qui détend tous. Les modèles dérivés du cadre proposé nous permettent d'enquêter sur la plausibilité biologique des trois hypothèses simplificatrices sur des données réelles ainsi que de trouver le meilleur modèle qui est aligné avec les caractéristiques sous-jacentes des jeux de données. Nos expériences montrent que, dans aucun des jeux de données réelles, tenant les trois hypothèses mentionnées est réaliste. Cela signifie en utilisant des modèles simples qui détiennent ces hypothèses peuvent être trompeuses et les résultats de l'estimation inexacte des paramètres. Le deuxième objectif est de développer un modèle mécaniste de codon généralisée qui détend les trois hypothèses simplificatrices, tandis que d'informatique efficace, en utilisant une opération de matrice appelée produit de Kronecker. Nos expériences montrent que sur un jeux de données choisis au hasard, le modèle proposé de codon mécaniste généralisée surpasse autre modèle de codon par rapport à AICc métrique dans environ la moitié des ensembles de données. En outre, je montre à travers plusieurs expériences que le modèle général proposé est biologiquement plausible.

L'intersubjectivité en questions. Introduction

La notion d'intersubjectivité est au coeur des débats contemporains dans le champ de la psychologie. L'intersubjectivité traverse en effet les différentes approches de la psychologie, des neurosciences à la psychologie sociale, des perspectives développementales aux approches psychoaffectives de la personnalité, des contextes de l'apprentissage à ceux du soin. Cet ouvrage collectif, premier volume de la collection Actualités psychologiques, se propose de faire le point sur cette notion et d'apporter des pistes actuelles de compréhension aux enjeux de la place de l'autre et de l'altérité dans la construction du sujet. Des chercheurs de renommée internationale, de Suisse et de différents pays européens, apportent ici leur contribution, dans différents domaines de la psychologie, invitant au dialogue entre les modèles théoriques qui sous-tendent la notion d'intersubjectivité et entre les méthodologies qui permettent d'en appréhender la dynamique. Contribution exceptionnelle, cet ouvrage est à destination de tous, psychologues praticiens et professionnels de la relation, chercheurs et étudiants en psychologie ou dans des disciplines connexes. Il permet d'approcher l'intersubjectivité comme une notion tout à la fois familière et complexe, au service de la rencontre de l'autre.

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.

Reliability of methods used for the risk assessment of mixture of micropollutants in aquatic enviromments : from theoretical concepts to field observations

Ces dernières années, de nombreuses recherches ont mis en évidence les effets toxiques des micropolluants organiques pour les espèces de nos lacs et rivières. Cependant, la plupart de ces études se sont focalisées sur la toxicité des substances individuelles, alors que les organismes sont exposés tous les jours à des milliers de substances en mélange. Or les effets de ces cocktails ne sont pas négligeables. Cette thèse de doctorat s'est ainsi intéressée aux modèles permettant de prédire le risque environnemental de ces cocktails pour le milieu aquatique. Le principal objectif a été d'évaluer le risque écologique des mélanges de substances chimiques mesurées dans le Léman, mais aussi d'apporter un regard critique sur les méthodologies utilisées afin de proposer certaines adaptations pour une meilleure estimation du risque. Dans la première partie de ce travail, le risque des mélanges de pesticides et médicaments pour le Rhône et pour le Léman a été établi en utilisant des approches envisagées notamment dans la législation européenne. Il s'agit d'approches de « screening », c'est-à-dire permettant une évaluation générale du risque des mélanges. Une telle approche permet de mettre en évidence les substances les plus problématiques, c'est-à-dire contribuant le plus à la toxicité du mélange. Dans notre cas, il s'agit essentiellement de 4 pesticides. L'étude met également en évidence que toutes les substances, même en trace infime, contribuent à l'effet du mélange. Cette constatation a des implications en terme de gestion de l'environnement. En effet, ceci implique qu'il faut réduire toutes les sources de polluants, et pas seulement les plus problématiques. Mais l'approche proposée présente également un biais important au niveau conceptuel, ce qui rend son utilisation discutable, en dehors d'un screening, et nécessiterait une adaptation au niveau des facteurs de sécurité employés. Dans une deuxième partie, l'étude s'est portée sur l'utilisation des modèles de mélanges dans le calcul de risque environnemental. En effet, les modèles de mélanges ont été développés et validés espèce par espèce, et non pour une évaluation sur l'écosystème en entier. Leur utilisation devrait donc passer par un calcul par espèce, ce qui est rarement fait dû au manque de données écotoxicologiques à disposition. Le but a été donc de comparer, avec des valeurs générées aléatoirement, le calcul de risque effectué selon une méthode rigoureuse, espèce par espèce, avec celui effectué classiquement où les modèles sont appliqués sur l'ensemble de la communauté sans tenir compte des variations inter-espèces. Les résultats sont dans la majorité des cas similaires, ce qui valide l'approche utilisée traditionnellement. En revanche, ce travail a permis de déterminer certains cas où l'application classique peut conduire à une sous- ou sur-estimation du risque. Enfin, une dernière partie de cette thèse s'est intéressée à l'influence que les cocktails de micropolluants ont pu avoir sur les communautés in situ. Pour ce faire, une approche en deux temps a été adoptée. Tout d'abord la toxicité de quatorze herbicides détectés dans le Léman a été déterminée. Sur la période étudiée, de 2004 à 2009, cette toxicité due aux herbicides a diminué, passant de 4% d'espèces affectées à moins de 1%. Ensuite, la question était de savoir si cette diminution de toxicité avait un impact sur le développement de certaines espèces au sein de la communauté des algues. Pour ce faire, l'utilisation statistique a permis d'isoler d'autres facteurs pouvant avoir une influence sur la flore, comme la température de l'eau ou la présence de phosphates, et ainsi de constater quelles espèces se sont révélées avoir été influencées, positivement ou négativement, par la diminution de la toxicité dans le lac au fil du temps. Fait intéressant, une partie d'entre-elles avait déjà montré des comportements similaires dans des études en mésocosmes. En conclusion, ce travail montre qu'il existe des modèles robustes pour prédire le risque des mélanges de micropolluants sur les espèces aquatiques, et qu'ils peuvent être utilisés pour expliquer le rôle des substances dans le fonctionnement des écosystèmes. Toutefois, ces modèles ont bien sûr des limites et des hypothèses sous-jacentes qu'il est important de considérer lors de leur application. - Depuis plusieurs années, les risques que posent les micropolluants organiques pour le milieu aquatique préoccupent grandement les scientifiques ainsi que notre société. En effet, de nombreuses recherches ont mis en évidence les effets toxiques que peuvent avoir ces substances chimiques sur les espèces de nos lacs et rivières, quand elles se retrouvent exposées à des concentrations aiguës ou chroniques. Cependant, la plupart de ces études se sont focalisées sur la toxicité des substances individuelles, c'est à dire considérées séparément. Actuellement, il en est de même dans les procédures de régulation européennes, concernant la partie évaluation du risque pour l'environnement d'une substance. Or, les organismes sont exposés tous les jours à des milliers de substances en mélange, et les effets de ces "cocktails" ne sont pas négligeables. L'évaluation du risque écologique que pose ces mélanges de substances doit donc être abordé par de la manière la plus appropriée et la plus fiable possible. Dans la première partie de cette thèse, nous nous sommes intéressés aux méthodes actuellement envisagées à être intégrées dans les législations européennes pour l'évaluation du risque des mélanges pour le milieu aquatique. Ces méthodes sont basées sur le modèle d'addition des concentrations, avec l'utilisation des valeurs de concentrations des substances estimées sans effet dans le milieu (PNEC), ou à partir des valeurs des concentrations d'effet (CE50) sur certaines espèces d'un niveau trophique avec la prise en compte de facteurs de sécurité. Nous avons appliqué ces méthodes à deux cas spécifiques, le lac Léman et le Rhône situés en Suisse, et discuté les résultats de ces applications. Ces premières étapes d'évaluation ont montré que le risque des mélanges pour ces cas d'étude atteint rapidement une valeur au dessus d'un seuil critique. Cette valeur atteinte est généralement due à deux ou trois substances principales. Les procédures proposées permettent donc d'identifier les substances les plus problématiques pour lesquelles des mesures de gestion, telles que la réduction de leur entrée dans le milieu aquatique, devraient être envisagées. Cependant, nous avons également constaté que le niveau de risque associé à ces mélanges de substances n'est pas négligeable, même sans tenir compte de ces substances principales. En effet, l'accumulation des substances, même en traces infimes, atteint un seuil critique, ce qui devient plus difficile en terme de gestion du risque. En outre, nous avons souligné un manque de fiabilité dans ces procédures, qui peuvent conduire à des résultats contradictoires en terme de risque. Ceci est lié à l'incompatibilité des facteurs de sécurité utilisés dans les différentes méthodes. Dans la deuxième partie de la thèse, nous avons étudié la fiabilité de méthodes plus avancées dans la prédiction de l'effet des mélanges pour les communautés évoluant dans le système aquatique. Ces méthodes reposent sur le modèle d'addition des concentrations (CA) ou d'addition des réponses (RA) appliqués sur les courbes de distribution de la sensibilité des espèces (SSD) aux substances. En effet, les modèles de mélanges ont été développés et validés pour être appliqués espèce par espèce, et non pas sur plusieurs espèces agrégées simultanément dans les courbes SSD. Nous avons ainsi proposé une procédure plus rigoureuse, pour l'évaluation du risque d'un mélange, qui serait d'appliquer d'abord les modèles CA ou RA à chaque espèce séparément, et, dans une deuxième étape, combiner les résultats afin d'établir une courbe SSD du mélange. Malheureusement, cette méthode n'est pas applicable dans la plupart des cas, car elle nécessite trop de données généralement indisponibles. Par conséquent, nous avons comparé, avec des valeurs générées aléatoirement, le calcul de risque effectué selon cette méthode plus rigoureuse, avec celle effectuée traditionnellement, afin de caractériser la robustesse de cette approche qui consiste à appliquer les modèles de mélange sur les courbes SSD. Nos résultats ont montré que l'utilisation de CA directement sur les SSDs peut conduire à une sous-estimation de la concentration du mélange affectant 5 % ou 50% des espèces, en particulier lorsque les substances présentent un grand écart- type dans leur distribution de la sensibilité des espèces. L'application du modèle RA peut quant à lui conduire à une sur- ou sous-estimations, principalement en fonction de la pente des courbes dose- réponse de chaque espèce composant les SSDs. La sous-estimation avec RA devient potentiellement importante lorsque le rapport entre la EC50 et la EC10 de la courbe dose-réponse des espèces est plus petit que 100. Toutefois, la plupart des substances, selon des cas réels, présentent des données d' écotoxicité qui font que le risque du mélange calculé par la méthode des modèles appliqués directement sur les SSDs reste cohérent et surestimerait plutôt légèrement le risque. Ces résultats valident ainsi l'approche utilisée traditionnellement. Néanmoins, il faut garder à l'esprit cette source d'erreur lorsqu'on procède à une évaluation du risque d'un mélange avec cette méthode traditionnelle, en particulier quand les SSD présentent une distribution des données en dehors des limites déterminées dans cette étude. Enfin, dans la dernière partie de cette thèse, nous avons confronté des prédictions de l'effet de mélange avec des changements biologiques observés dans l'environnement. Dans cette étude, nous avons utilisé des données venant d'un suivi à long terme d'un grand lac européen, le lac Léman, ce qui offrait la possibilité d'évaluer dans quelle mesure la prédiction de la toxicité des mélanges d'herbicide expliquait les changements dans la composition de la communauté phytoplanctonique. Ceci à côté d'autres paramètres classiques de limnologie tels que les nutriments. Pour atteindre cet objectif, nous avons déterminé la toxicité des mélanges sur plusieurs années de 14 herbicides régulièrement détectés dans le lac, en utilisant les modèles CA et RA avec les courbes de distribution de la sensibilité des espèces. Un gradient temporel de toxicité décroissant a pu être constaté de 2004 à 2009. Une analyse de redondance et de redondance partielle, a montré que ce gradient explique une partie significative de la variation de la composition de la communauté phytoplanctonique, même après avoir enlevé l'effet de toutes les autres co-variables. De plus, certaines espèces révélées pour avoir été influencées, positivement ou négativement, par la diminution de la toxicité dans le lac au fil du temps, ont montré des comportements similaires dans des études en mésocosmes. On peut en conclure que la toxicité du mélange herbicide est l'un des paramètres clés pour expliquer les changements de phytoplancton dans le lac Léman. En conclusion, il existe diverses méthodes pour prédire le risque des mélanges de micropolluants sur les espèces aquatiques et celui-ci peut jouer un rôle dans le fonctionnement des écosystèmes. Toutefois, ces modèles ont bien sûr des limites et des hypothèses sous-jacentes qu'il est important de considérer lors de leur application, avant d'utiliser leurs résultats pour la gestion des risques environnementaux. - For several years now, the scientists as well as the society is concerned by the aquatic risk organic micropollutants may pose. Indeed, several researches have shown the toxic effects these substances may induce on organisms living in our lakes or rivers, especially when they are exposed to acute or chronic concentrations. However, most of the studies focused on the toxicity of single compounds, i.e. considered individually. The same also goes in the current European regulations concerning the risk assessment procedures for the environment of these substances. But aquatic organisms are typically exposed every day simultaneously to thousands of organic compounds. The toxic effects resulting of these "cocktails" cannot be neglected. The ecological risk assessment of mixtures of such compounds has therefore to be addressed by scientists in the most reliable and appropriate way. In the first part of this thesis, the procedures currently envisioned for the aquatic mixture risk assessment in European legislations are described. These methodologies are based on the mixture model of concentration addition and the use of the predicted no effect concentrations (PNEC) or effect concentrations (EC50) with assessment factors. These principal approaches were applied to two specific case studies, Lake Geneva and the River Rhône in Switzerland, including a discussion of the outcomes of such applications. These first level assessments showed that the mixture risks for these studied cases exceeded rapidly the critical value. This exceeding is generally due to two or three main substances. The proposed procedures allow therefore the identification of the most problematic substances for which management measures, such as a reduction of the entrance to the aquatic environment, should be envisioned. However, it was also showed that the risk levels associated with mixtures of compounds are not negligible, even without considering these main substances. Indeed, it is the sum of the substances that is problematic, which is more challenging in term of risk management. Moreover, a lack of reliability in the procedures was highlighted, which can lead to contradictory results in terms of risk. This result is linked to the inconsistency in the assessment factors applied in the different methods. In the second part of the thesis, the reliability of the more advanced procedures to predict the mixture effect to communities in the aquatic system were investigated. These established methodologies combine the model of concentration addition (CA) or response addition (RA) with species sensitivity distribution curves (SSD). Indeed, the mixture effect predictions were shown to be consistent only when the mixture models are applied on a single species, and not on several species simultaneously aggregated to SSDs. Hence, A more stringent procedure for mixture risk assessment is proposed, that would be to apply first the CA or RA models to each species separately and, in a second step, to combine the results to build an SSD for a mixture. Unfortunately, this methodology is not applicable in most cases, because it requires large data sets usually not available. Therefore, the differences between the two methodologies were studied with datasets created artificially to characterize the robustness of the traditional approach applying models on species sensitivity distribution. The results showed that the use of CA on SSD directly might lead to underestimations of the mixture concentration affecting 5% or 50% of species, especially when substances present a large standard deviation of the distribution from the sensitivity of the species. The application of RA can lead to over- or underestimates, depending mainly on the slope of the dose-response curves of the individual species. The potential underestimation with RA becomes important when the ratio between the EC50 and the EC10 for the dose-response curve of the species composing the SSD are smaller than 100. However, considering common real cases of ecotoxicity data for substances, the mixture risk calculated by the methodology applying mixture models directly on SSDs remains consistent and would rather slightly overestimate the risk. These results can be used as a theoretical validation of the currently applied methodology. Nevertheless, when assessing the risk of mixtures, one has to keep in mind this source of error with this classical methodology, especially when SSDs present a distribution of the data outside the range determined in this study Finally, in the last part of this thesis, we confronted the mixture effect predictions with biological changes observed in the environment. In this study, long-term monitoring of a European great lake, Lake Geneva, provides the opportunity to assess to what extent the predicted toxicity of herbicide mixtures explains the changes in the composition of the phytoplankton community next to other classical limnology parameters such as nutrients. To reach this goal, the gradient of the mixture toxicity of 14 herbicides regularly detected in the lake was calculated, using concentration addition and response addition models. A decreasing temporal gradient of toxicity was observed from 2004 to 2009. Redundancy analysis and partial redundancy analysis showed that this gradient explains a significant portion of the variation in phytoplankton community composition, even when having removed the effect of all other co-variables. Moreover, some species that were revealed to be influenced positively or negatively, by the decrease of toxicity in the lake over time, showed similar behaviors in mesocosms studies. It could be concluded that the herbicide mixture toxicity is one of the key parameters to explain phytoplankton changes in Lake Geneva. To conclude, different methods exist to predict the risk of mixture in the ecosystems. But their reliability varies depending on the underlying hypotheses. One should therefore carefully consider these hypotheses, as well as the limits of the approaches, before using the results for environmental risk management

Genomic copy number aberrations in glioblastoma, response to therapy and molecular machanisms

Résumé : Le glioblastome (GBM, WHO grade IV) est la tumeur cérébrale primaire la plus fréquente et la plus maligne, son pronostic reste très réservé et sa réponse aux différents traitements limitée. Récemment, une étude clinique randomisée (EORTC 26981/NCIC CE.3) a démontré que le traitement combiné de temozolomide et radiothérapie (RT/TMZ) est le meilleur dans les cas de GBM nouvellement diagnostiqués [1]. Cependant, seul un sous-groupe de patients bénéficie du traitement RT/TMZ et même parmi eux, leur survie reste très limitée. Pour tenter de mieux comprendre les réponses au traitement RT/TMZ, la biologie du GBM, identifier d'autres facteurs de résistance et découvrir de nouvelles cibles aux traitements, nous avons conduit une analyse moléculaire étendue à 73 patients inclus dans cette étude clinique. Nous avons complété les résultats moléculaires déjà obtenus par un profil génomique du nombre de copies par Array Comparative Genomic Hybridization. Afin d'atteindre nos objectifs, nous avons analysé en parallèle les données cliniques des patients et leurs profils moléculaires. Nos résultats confirment des analyses connues dans le domaine des aberrations du nombre de copies (CNA) et de profils du glioblastome. Nous avons observé une bonne corrélation entre le CNA génomique et l'expression de l'ARN messager dans le glioblastome et identifié un nouveau modèle de CNA du chromosome 7 pouvant présenter un intérêt clinique. Nous avons aussi observé par l'analyse du CNA que moins de 10% des glioblastomes conservent leurs mécanismes de suppression de tumeurs p53 et Rb1. Nous avons aussi observé que l'amplification du CDK4 peut constituer un facteur supplémentaire de résistance au traitement RT/TMZ, cette observation nécessite confirmation sur un plus grand nombre d'analyses. Nous avons montré que dans notre analyse des profils moléculaires et cliniques, il n'est pas possible de différencier le GBM à composante oligodendrogliale (GBM-O) du glioblastome. En superposant les profils moléculaires et les modèles expérimentaux in vitro, nous avons identifié WIF-1 comme un gène suppresseur de tumeur probable et une activation du signal WNT dans la pathologie du glioblastome. Ces observations pourraient servir à une meilleure compréhension de cette maladie dans le futur. Abstract : Glioblastoma, (GBM, WHO grade IV) is the most malignant and most frequent primary brain tumor with a very poor prognosis and response to therapy. A recent randomized clinical trial (EORTC26981/NCIC CE.3) established RT/TMZ as the 1St effective chemo-radiation therapy in newly diagnosed GBM [1]. However only a genetic subgroup of patients benefit from RT/TMZ and even in this subgroup overall survival remains very dismal. To explain the observed response to RT/TMZ, have a better understanding of GBM biology, identify other resistance factors and discover new drugable targets a comprehensive molecular analysis was performed in 73 of these GBM trial cohort. We complemented the available molecular data with a genomic copy number profiling by Array Comparative Genomic Hybridization. We proceeded to align the molecular profiles and the Clinical data, to meet our project objectives. Our data confirm known GBM Copy Number Aberrations and profiles. We observed a good correlation of genomic CN and mRNA expression in GBM, and identified new interesting CNA pattern for chromosome 7 with a potential clinical value. We also observed that by copy number aberration data alone, less than 10% of GBM have an intact p53 and Rb1 tumor .suppressor pathways. We equally observed that CDK4 amplification might constitute an additional RT/TMZ resistant factor, an observation that will need confirmation in a larger data set. We show that the molecular and clinical profiles in our data set, does not support the identification of GBM-O as a new entity in GBM. By combining the molecular profiles and in vitro model experiments we identify WIF1 as a potential GBM TSG and an activated WNT signaling as a pathologic event in GBM worth incorporation in attempts to better understand and impact outcome in this disease.

Plate tectonics of the Altaids

Abstract: The Altaids consist in a huge accretionary-type belt extending from Siberia through Mon-golia, northern China, Kyrgyzstan and Kazakhstan. They were formed from the Vendian through the Jurassic by the accretion of numerous displaced and exotic terranes (e.g. island arc, ribbon microcontinent, seamount, basaltic plateau, back-arc basin). The number, nature and origin of the terranes differ according to the palaeotectonic models of the different authors. Thanks to a geo- dynamic study (i.e. definition of tectonic settings and elaboration of geodynamic scenarios) and plate tectonics modelling, this work aims to present an alternative model explaining the Palaeozoic palaeotectonic evolution of the Altaids. Based on a large set of compiled geological data related to palaeogeography and geodyna¬mic (e.g. sedimentology, stratigraphy, palaeobiogeography, palaeomagnetism, magmatism, me- tamorphism, tectonic...), a partly new classification of the terranes and sutures implicated in the formation of the Altaids is proposed. In the aim to elaborate plate tectonics reconstructions, it is necessary to fragment the present arrangement of continents into consistent geological units. To avoid confusion with existing terminology (e.g. tectonic units, tectono-stratigraphic units, micro- continents, terranes, blocks...), the new concept of "Geodynamic Units (GDU)" was introduced. A terrane may be formed by a set of GDUs. It consists of a continental and/or oceanic fragment which has its own kinematic and geodynamic evolution for a given period. With the same ap-proach, the life span and type of the disappeared oceans is inferred thanks to the study of the mate-rial contained in suture zones. The interpretation of the tectonic settings within the GDUs comple-ted by the restoration of oceans leads to the elaboration of geodynamic scenarios. Since the Wilson cycle was presented in 1967, numerous works demonstrated that the continental growth is more complex and results from diverse geodynamic scenarios. The identification of these scenarios and their exploitation enable to elaborate plate tectonics models. The models are self-constraining (i.e. space and time constraints) and contest or confirm in turn the geodynamic scenarios which were initially proposed. The Altaids can be divided into three domains: (1) the Peri-Siberian, (2) the Kazakhstan, and (3) the Tarim-North China domains. The Peri-Siberian Domain consists of displaced (i.e. Sayan Terrane Tuva-Mongolian, Lake-Khamsara Terrane) and exotic terranes (i.e. Altai-Mongolian and Khangai-Argunsky Terrane) accreted to Siberia from the Vendian through the Ordovician. Fol-lowing the accretion of these terranes, the newly formed Siberia active margin remained active un-til its part collision with the Kazakhstan Superterrane in the Carboniferous. The eastern part of the active margin (i.e. East Mongolia) continued to act until the Permian when the North-China Tarim Superterrane collided with it. The geodynamic evolution of the eastern part of the Peri-Siberian Domain (i.e. Eastern Mongolia and Siberia) is complicated by the opening of the Mongol-Okhotsk Ocean in the Silurian. The Kazakhstan Domain is composed of several continental terranes of East Gondwana origin amalgamated together during the Ordovician-Silurian time. After these different orogenic events, the Kazakhstan Superterrane evolved as a single superterrane until its collision with a Tarim-North China related-terrane (i.e. Tianshan-Hanshan Terrane) and Siberian Continent during the Devonian. This new organisation of the continents imply a continued active margin from Siberia, to North China through the Kazakhstan Superterrane and the closure of the Junggar- Balkash Ocean which implied the oroclinal bending of the Kazakhstan Superterrane during the entire Carboniferous. The formation history of the Tarim-North China Domain is less complex. The Cambrian northern passive margin became active in the Ordovician. In the Silurian, the South Tianshan back-arc Ocean was open and led to the formation of the Tianshan-Hanshan Terrane which collided with the Kazakhstan Superterrane during the Devonian. The collision between Siberia and the eastern part of the Tarim-North China continents (i.e. Inner Mongolia), implied by the closure of the Solonker Ocean, took place in the Permian. Since this time, the major part of the Altaids was formed, the Mongol-Okhotsk Ocean only was still open and closed during the Jurassic. Résumé: La chaîne des Altaïdes est une importante chaîne d'accrétion qui s'étend en Sibérie, Mon-golie, Chine du Nord, Kirghizstan et Kazakhstan. Elle s'est formée durant la période du Vendian au Jurassique par l'accrétion de nombreux terranes déplacés ou exotiques (par exemple arc océa-nique, microcontinent, guyot, plateau basaltique, basin d'arrière-arc...). Le nombre, la nature ou encore l'origine diffèrent selon les modèles paléo-tectoniques proposés par les différents auteurs. Grâce à une étude géodynamique (c'est-à-dire définition des environnements tectoniques et éla-boration de scénarios géodynamiques) et à la modélisation de la tectonique des plaques, ce travail propose un modèle alternatif expliquant l'évolution paléo-tectonique des Altaïdes. Basé sur une large compilation de données géologiques pertinentes en termes de paléo-géographie et de géodynamique (par exemple sédimentologie, stratigraphie, paléo-biogéographie, paléomagnétisme, magmatisme, métamorphisme, tectonique...), une nouvelle classification des terranes et des sutures impliqués dans la formation des Altaïdes est proposée. Dans le but d'élabo¬rer des reconstructions de plaques tectoniques, il est nécessaire de fragmenter l'arrangement actuel des continents en unités tectoniques cohérentes. Afin d'éviter les confusions avec la terminolo¬gie existante (par exemple unité tectonique, unité tectono-stratigraphique, microcontinent, block, terrane...), le nouveau concept d' "Unité Géodynamique (UGD)" a été introduit. Un terrane est formé d'une ou plusieurs UGD et représente un fragment océanique ou continental défini pas sa propre cinétique et évolution géodynamique pour une période donnée. Parallèlement, la durée de vie et le type des océans disparus (c'est-à-dire principal ou secondaire) est déduite grâce à l'étude du matériel contenu dans les zones de sutures. L'interprétation des environnements tectoniques des UGD associés à la restauration des océans mène à l'élaboration de scénarios géodynamiques. Depuis que le Cycle de Wilson a été présenté en 1967, de nombreux travaux ont démontré que la croissance continentale peut résulter de divers scénarios géodynamiques. L'identification et l'ex-ploitation de ces scénarios permet finalement l'élaboration de modèles de tectonique des plaques. Les modèles sont auto-contraignants (c'est-à-dire contraintes spatiales et temporelles) et peuvent soit contester ou confirmer les scénarios géodynamiques initialement proposés. Les Altaïdes peuvent être divisées en trois domaines : (1) le Domaine Péri-Sibérien, (2) le Domaine Kazakh, et (3) le Domaine Tarim-Nord Chinois. Le Domaine Péri-Sibérien est composé de terranes déplacés (c'est-à-dire Terrane du Sayan, Tuva-Mongol et Lake-Khamsara) et exotiques (c'est-à-dire Terrane Altai-Mongol et Khangai-Argunsky) qui ont été accrétés au craton Sibérien durant la période du Vendien à l'Ordovicien. Suite à l'accrétion de ces terranes, la marge sud-est de la Sibérie nouvellement formée reste active jusqu'à sa collision partielle avec le Superterrane Ka-zakh au Carbonifère. La partie est de la marge active (c'est-à-dire Mongolie de l'est) continue son activité jusqu'au Permien lors de sa collision avec le Superterrane Tarim-Nord Chinois. L'évolu¬tion géodynamique de la partie est du Domaine Sibérien est compliquée par l'ouverture Silurienne de l'Océan Mongol-Okhotsk qui disparaîtra seulement au Jurassique. Le Domaine Kazakh est composé de plusieurs terranes d'origine est-Gondwanienne accrétés les uns avec les autres avant ou pendant le Silurien inférieur et leurs evolution successive sous la forme d'un seul superterrane. Le Superterrane Kazakh collisione avec un terrane Tarim-Nord Chinois (c'est-à-dire Terrane du Tianshan-Hanshan) durant le Dévonien et le continent Sibérien au Dévonien supérieur. Ce nouvel agencement des plaques induit une marge active continue le long des continents Sibérien, Kazakh et Nord Chinois et la fermeture de l'Océan Junggar-Balkash qui provoque le plissement oroclinal du Superterrane Kazakh durant le Carbonifère. L'histoire de la formation du Domaine Tarim-Nord Chinois est moins complexe. La marge passive nord Cambrienne devient active à l'Ordovicien et l'ouverture Silurienne du bassin d'arrière-arc du Tianshan sud mène à la formation du terrane du Tianshan-Hanshan. La collision Dévonienne entre ce dernier et le Superterrane Kazakh provoque la fermerture de l'Océan Tianshan sud. Finalement, la collision entre la Sibérie et la partie est du continent Tarim-Nord Chinois (c'est-à-dire Mongolie Intérieure) prend place durant le Permien suite à la fermeture de l'Océan Solonker. La majeure partie des Altaïdes est alors formée, seul l'Océan Mongol-Okhotsk est encore ouvert. Ce dernier se fermera seulement au Jurassique.

Niche, distribution and global change : modeling insights into biogeography and conservation biology Olivier Broennimann

RESUME Les évidences montrant que les changements globaux affectent la biodiversité s'accumulent. Les facteurs les plus influant dans ce processus sont les changements et destructions d'habitat, l'expansion des espèces envahissantes et l'impact des changements climatiques. Une évaluation pertinente de la réponse des espèces face à ces changements est essentielle pour proposer des mesures permettant de réduire le déclin actuel de la biodiversité. La modélisation de la répartition d'espèces basée sur la niche (NBM) est l'un des rares outils permettant cette évaluation. Néanmoins, leur application dans le contexte des changements globaux repose sur des hypothèses restrictives et demande une interprétation critique. Ce travail présente une série d'études de cas investiguant les possibilités et limitations de cette approche pour prédire l'impact des changements globaux. Deux études traitant des menaces sur les espèces rares et en danger d'extinction sont présentées. Les caractéristiques éco-géographiques de 118 plantes avec un haut degré de priorité de conservation sont revues. La prévalence des types de rareté sont analysées en relation avec leur risque d'extinction UICN. La revue souligne l'importance de la conservation à l'échelle régionale. Une évaluation de la rareté à échelle globale peut être trompeuse pour certaine espèces car elle ne tient pas en compte des différents degrés de rareté que présente une espèce à différentes échelles spatiales. La deuxième étude test une approche pour améliorer l'échantillonnage d'espèces rares en incluant des phases itératives de modélisation et d'échantillonnage sur le terrain. L'application de l'approche en biologie de la conservation (illustrée ici par le cas du chardon bleu, Eryngium alpinum), permettrait de réduire le temps et les coûts d'échantillonnage. Deux études sur l'impact des changements climatiques sur la faune et la flore africaine sont présentées. La première étude évalue la sensibilité de 227 mammifères africains face aux climatiques d'ici 2050. Elle montre qu'un nombre important d'espèces pourrait être bientôt en danger d'extinction et que les parcs nationaux africains (principalement ceux situé en milieux xériques) pourraient ne pas remplir leur mandat de protection de la biodiversité dans le futur. La seconde étude modélise l'aire de répartition en 2050 de 975 espèces de plantes endémiques du sud de l'Afrique. L'étude propose l'inclusion de méthodes améliorant la prédiction des risques liés aux changements climatiques. Elle propose également une méthode pour estimer a priori la sensibilité d'une espèce aux changements climatiques à partir de ses propriétés écologiques et des caractéristiques de son aire de répartition. Trois études illustrent l'utilisation des modèles dans l'étude des invasions biologiques. Une première étude relate l'expansion de la laitue sáuvage (Lactuca serriola) vers le nord de l'Europe en lien avec les changements du climat depuis 250 ans. La deuxième étude analyse le potentiel d'invasion de la centaurée tachetée (Centaures maculosa), une mauvaise herbe importée en Amérique du nord vers 1890. L'étude apporte la preuve qu'une espèce envahissante peut occuper une niche climatique différente après introduction sur un autre continent. Les modèles basés sur l'aire native prédisent de manière incorrecte l'entier de l'aire envahie mais permettent de prévoir les aires d'introductions potentielles. Une méthode alternative, incluant la calibration du modèle à partir des deux aires où l'espèce est présente, est proposée pour améliorer les prédictions de l'invasion en Amérique du nord. Je présente finalement une revue de la littérature sur la dynamique de la niche écologique dans le temps et l'espace. Elle synthétise les récents développements théoriques concernant le conservatisme de la niche et propose des solutions pour améliorer la pertinence des prédictions d'impact des changements climatiques et des invasions biologiques. SUMMARY Evidences are accumulating that biodiversity is facing the effects of global change. The most influential drivers of change in ecosystems are land-use change, alien species invasions and climate change impacts. Accurate projections of species' responses to these changes are needed to propose mitigation measures to slow down the on-going erosion of biodiversity. Niche-based models (NBM) currently represent one of the only tools for such projections. However, their application in the context of global changes relies on restrictive assumptions, calling for cautious interpretations. In this thesis I aim to assess the effectiveness and shortcomings of niche-based models for the study of global change impacts on biodiversity through the investigation of specific, unsolved limitations and suggestion of new approaches. Two studies investigating threats to rare and endangered plants are presented. I review the ecogeographic characteristic of 118 endangered plants with high conservation priority in Switzerland. The prevalence of rarity types among plant species is analyzed in relation to IUCN extinction risks. The review underlines the importance of regional vs. global conservation and shows that a global assessment of rarity might be misleading for some species because it can fail to account for different degrees of rarity at a variety of spatial scales. The second study tests a modeling framework including iterative steps of modeling and field surveys to improve the sampling of rare species. The approach is illustrated with a rare alpine plant, Eryngium alpinum and shows promise for complementing conservation practices and reducing sampling costs. Two studies illustrate the impacts of climate change on African taxa. The first one assesses the sensitivity of 277 mammals at African scale to climate change by 2050 in terms of species richness and turnover. It shows that a substantial number of species could be critically endangered in the future. National parks situated in xeric ecosystems are not expected to meet their mandate of protecting current species diversity in the future. The second study model the distribution in 2050 of 975 endemic plant species in southern Africa. The study proposes the inclusion of new methodological insights improving the accuracy and ecological realism of predictions of global changes studies. It also investigates the possibility to estimate a priori the sensitivity of a species to climate change from the geographical distribution and ecological proprieties of the species. Three studies illustrate the application of NBM in the study of biological invasions. The first one investigates the Northwards expansion of Lactuca serriola L. in Europe during the last 250 years in relation with climate changes. In the last two decades, the species could not track climate change due to non climatic influences. A second study analyses the potential invasion extent of spotted knapweed, a European weed first introduced into North America in the 1890s. The study provides one of the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. Models fail to predict the current full extent of the invasion, but correctly predict areas of introduction. An alternative approach, involving the calibration of models with pooled data from both ranges, is proposed to improve predictions of the extent of invasion on models based solely on the native range. I finally present a review on the dynamic nature of ecological niches in space and time. It synthesizes the recent theoretical developments to the niche conservatism issues and proposes solutions to improve confidence in NBM predictions of the impacts of climate change and species invasions on species distributions.

Radiolarian biochronology and Pacific origin of Jurassic-Cretaceous oceanic terranes in Costa Rica and Puerto Rico

AbstractAs demonstrated during several recent geological conferences, there is still a large debate concerning the origins of the Mesozoic oceanic remnants on the Caribbean Plate. The geodynamic models describing the Mesozoic history of the Caribbean realm can be divided into two main categories based on the origin of the Caribbean Plate: 1) An in situ origin between the Americas; 2) A Pacific origin and an eastward transport relative to the Americas. The study of the ribbon-bedded radiolarite is a key in determining the origins of associated Mesozoic oceanic terranes and may help to achieve a general agreement regarding the basic principles on the evolution of the Caribbean Plate. The Early Jurassic to early Late Cretaceous Bermeja Complex of Puerto Rico, witch contains serpentinized peridotite, altered basalt, amphibolite, and chert (Mariquita Chert Formation), and the contemporaneous Santa Rosa Accretionary Complex, which crops out in several half-windows along the south shores of the Santa Elena Peninsula in northwestern Costa Rica, are two of these little-known and crucial ophiolitic mélanges. The Manzanillo and Matambú fore-arc Terranes of the Nicoya Peninsula in the northwestern Costa Rica, which contain Late Cretaceous to Early Paleogene radiolarian-bearing siliceous mudstones and cherts associated with arc-derived mafic to intermediate volcaniclastics, bring important information on the history of the western active margin of the Caribbean Plate. A systematic radiolarian study of these three regions is presented herein in three different articles.The radiolarian biochronology of the Mariquita Chert Formation of the Bermeja Complex presented in this work indicate an early Middle Jurassic to early Late Cretaceous (late Bajocian-early Callovian to middle Albian-middle Cenomanian) age for the Mariquita Chert Formation. The illustrated assemblages contain 150 species, of which 3 are new (Pantanellium karinae, Loopus bermejaense, and L. boricus), and belonging to 59 genera. A review of the previous radiolarian published works on this formation and the results of this study suggest that the Bermeja Complex ranges in age from Middle Jurassic to early Late Cretaceous (late Aalenian to middle Cenomanian) and also reveal a possible feature of the complex, which is the youngling of radiolarian cherts from north to south, evoking a polarity of accretion. On the basis of a currently exhaustive inventory of the ribbonbedded radiolaritic facies on the Caribbean Plate, a re-examination of the distribution of Middle Jurassic sediments associated with oceanic crust from the Caribbean realm, and a paleoceanographical argumentation on the water currents, we come to the conclusion that the radiolarite and associated Mesozoic oceanic terranes of the Caribbean Plate are of Pacific origin. The previous argument for a Pacific origin of the Bermeja Complex presented by Montgomery et al. (1994a), based on their radiolarian age and their estimation of the oldest Proto-Caribbean oceanic crust, is nowadays seriously questionable, owing to the recent progresses in radiolarian biostratigraphy and new discoveries on the age of the first oceanic crust spreading between the Americas. Furthermore, we interpret the radiolarian Parvicingulidae-rich assemblages in the low-latitude Caribbean context as potential indicators of upwelling or land nutrients inputs, instead of indicators of paleolatitudes,as firstly stated by Pessagno and Blome (1986). Eventually, a discussion on the origin of the cherts of the Mariquita Formation illustrated by Middle Jurassic to middle Cretaceous geodynamic models of the Pacific and Caribbean realms bring up the possibility that the rocks of the Bermeja Complex are remnants of two different oceans.The Santa Rosa Accretionary Complex contains various oceanic assemblages of alkaline basalt, radiolarite and polymictic breccias. The radiolarian biochronology (19 illustrated assemblages, 232 species belonging to 63 genera) presented in this work indicate an Early Jurassic to early Late Cretaceous (early Pliensbachian to earliest Turonian) age for the sediments associated with oceanic basalts or recovered from blocks in breccias or megabreccias from the Santa Rosa Accretionary Complex. This study brings to light the Early Jurassic age of a sequence of ribbon-bedded radiolarite, which was previously thought to be of Cretaceous age, intruded by alkaline basalts sills. The presence of Early Jurassic large reworked blocks of radiolarite in a polymictic megabreccia, firstly reported by De Wever et al. (1985) is confirmed. Therefore, the alkaline basalt associated with these radiolarites could be of Jurassic age. In the Carrizal tectonic window, Middle Jurassic radiolarian chert blocks and Early Cretaceous brick-red ribbon-bedded radiolarites overlying pillow basalts are interpreted as fragments of a Middle Jurassic oceanic basement accreted to an Early Cretaceous oceanic plate, in an intra-oceanic subduction context. Whereas, knobby radiolarites and black shale at Playa Carrizal are indicative of a shallower middle Cretaceous paleoenvironment. Other younger oceanic remnants documented the rapid approach of the site of sedimentation to a subduction trench during the late Early Cretaceous (AlbianCenomanian), maybe early Late Cretaceous (Turonian).In total, 60 species belonging to 34 genera were present in relatively well-preserved radiolarian faunas from volcaniclastics and associated pelagic and hemipelagic rocks of the Matambú and Manzanillo terranes, ranging in age from Late Cretaceous to Early Paleogene (middle Turonian-Santonian to late Thanetian-Ypresian). This study shows that radiolarians can provide significant biostratigraphic control in the Nicoya Peninsula where very similar lithologies of different ages are present. Two radiolarian samples directly date the Berrugate Formation for the first time (middle Turonian-Santonian and Coniacian-Santonian). These ages allow to determine a volcanic arc activity on the western edge of the future Caribbean Plate at least since the Santonian that could have lasted through the middle Turonian-early Campanian interval by stratigraphic superposition. Moreover on the basis of these radiolarian ages, the Loma Chumico Formation of Albian age, and the Berrugate Formation of middle Turonian-early Maastrichtian age, can now be clearly differentiated. Two samples from the Sabana Grande Formation give a Coniacian-Santonian age and a Coniacian-Campanian age and indicate that there is a stratigraphic gap of ~10 million years between this formation and the underlying Albian Loma Chumico Formation.RésuméComme cela a pu se vérifier à plusieurs reprises lors de conférences géologiques récentes, le débat sur l'origine des terrains océaniques mésozoïques de la Plaque Caraïbes est toujours d'actualité. Les modèles géodynamiques décrivant l'histoire de la région caraïbes peuvent être classés en deux catégories basées sur l'origine de la Plaque Caraïbes : 1) Une origine in situ entre les Amériques ; 2) Une origine Pacifique et un transport vers l'est, par rapport aux Amériques. L'étude des radiolarites rubanées est capitale pour la détermination de l'origine des terrains océaniques allochtones du Mésozoïque et peut être utile pour parvenir à un compromis général concernant les principes basiques de l'évolution de la Plaque Caraïbes. Le complexe de Bermeja à Porto Rico qui est constitué de péridotites serpentinisées, de basaltes altérés, d'amphibolites et de cherts (Formation des Cherts de Mariquita), et le Complexe d'Accrétion de Santa Rosa qui affleure dans plusieurs demi-fenêtres tectoniques au sud de la Péninsule de Santa Elena au nord-ouest du Costa Rica sont deux de ces mélanges ophiolitiques peu décrits et déterminants. Les terrains de fore-arc de Manzanillo et de Matambu dans la Péninsule de Nicoya au nord-ouest du Costa Rica qui sont composés de calcaires siliceux et de cherts riches en radiolaires associés à du matériel volcanique d'arc mafique à intermédiaire, apportent d'importantes informations sur l'histoire de la marge active occidentale de la Plaque Caraïbe. Une étude systématique des radiolaires de ces trois régions est présentée dans ce travail sous forme de trois articles.La biochronologie des radiolaires de la Formation des Cherts de Mariquita du Complexe d'Accrétion de Santa Rosa présentée dans ce travail indique un âge Jurassique Moyen inférieur à Crétacé Supérieur inférieur (Bajocien supérieur-Callovien inférieur à Albien moyen-Cénomanien moyen) pour la Formation des Cherts de Mariquita. Les assemblages illustrés contiennent 150 espèces, parmis lesquelles 3 sont nouvelles (Pantanellium karinae, Loopus bermejaense et L. boricus), et appartenant à 59 genres différents. Une révision des travaux publiés précédemment sur les radiolaires de cette formation, ainsi que les résultats de cette étude suggèrent que le Complexe de Bermeja a un âge allant du Jurassique moyen au Crétacé Supérieur inférieur (Aalénien supérieur à Cénomanien moyen) et révèle aussi une caractéristique éventuelle du complexe qui est le rajeunissement des radiolarites du nord au sud, évoquant une polarité d'accrétion. Sur la base d'un inventaire actuellement exhaustif du facies radiolaritique rubané sur la Plaque Caraïbes, d'un nouvel examen de la distribution globale des sédiments du Jurassique Moyen associés à de la croûte océanique et d'une argumentation paléocéanographique sur les courants, nous arrivons à la conclusion que les radiolarites et les unités tectoniques océaniques du Mésozoïque associées de la Plaque Caraïbes sont d'origine pacifique. L'argument antérieur pour une origine pacifique du Complexe de Bermeja présenté par Montgomery et al. (1994a), basé sur leur âge à radiolaire et leur estimation de l'âge de la plus vieille croûte océanique des Proto-Caraïbes, est sérieusement remis en question aujourd'hui, en raison des progrès récents de la biostratigraphie des radiolaires et des nouvelles découvertes concernant l'âge du début de l'océanisation entre les Amériques. En outre, dans le contexte de basses latitudes des Caraïbes, nous interprétons les assemblages à radiolaires riches en Parvicingulidae comme étant des indicateurs potentiels d'apports en nutriments des zones d'uppwelling ou des terres, plutôt que des indicateurs de paléolatitudes, comme exposer pour la première fois par Pessagno et Blome (1986). Finalement, une discussion sur l'origine des cherts de la Formation de Mariquita illustrée par des modèles géodynamiques du Jurassique Moyen au Crétacé moyen des régions pacifique et caraïbes, fait poindre la possibilité que les roches du Complexe de Bermeja proviennent de deux océans différents.Le Complexe d'Accrétion de Santa Rosa contient plusieurs assemblages océaniques différents de basaltes alcalins, radiolarites et brèches polymictes. La biochronologie des radiolaires (19 assemblages illustrés, 232 espèces appartenant à 63 genres) présentée dans ce second travail indique un âge Jurassique Inférieur à Crétacé Supérieur inférieur (Pliensbachien inférieur à Turonien initial) pour les sédiments associés aux basaltes océaniques ou provenant de blocs dans des brèches ou des mégabrèches du Complexe d'Accrétion de Santa Rosa. Cette étude met en évidence l'âge Jurassique Inférieur d'une séquence de radiolarites rubanées entrecoupée de sills de basaltes alcalins, dont l'âge estimé était précédemment le Crétacé.La présence de blocs plurimétriques de radiolarites d'âge Jurassique Inférieur remaniés dans une mégabrèche polymicte, dont la présence avait été signalée par De Wever et al. (1985), est confirmée. Par conséquent, les basaltes alcalins associés à ces radiolarites pourraient aussi être d'âge Jurassique. Dans la fenêtre tectonique de Carrizal, des blocs de radiolarites d'âge Jurassique Moyen et des radiolarites du Crétacé Inférieur recouvrant des basaltes en coussins sont interprétés comme des fragments d'une croûte océanique d'âge Jurassique Moyen accrétés à une plaque océanique d'âge Crétacé Inférieur, dans un contexte de subduction intra-océanique. Alors que dans la même zone, les radiolarites « noueuses » et les argiles noires associées sont interprétées comme des indicateurs d'un milieu peu profond au Crétacé. D'autres fragments océaniques plus jeunes documentent une approche rapide du lieu de sédimentation vers une fosse de subduction pendant le Crétacé Inférieur supérieur (Albien-Cénomanien), peut-être Crétacé Supérieur (Turonien).Au total, 60 espèces appartenant à 34 genres ont été déterminées à partir de faunes à radiolaires relativement bien préservées, extraites de roches volcanoclastiques et pélagiques à hémipélagiques associées, provenant des terrains de Matambu et Manzanillo et ayant des âges compris entre le Crétacé Supérieur et le Paléogène Inférieur (Turonien moyen-Santonien à Thanétien supérieur-Yprésien). Cette étude montre que les radiolaires peuvent fournir un contrôle stratigraphique significatif dans la Péninsule de Nicoya, où des lithologies similaires, mais d'âges différents sont présentes. Deux échantillons à radiolaires permettent de dater la Formation de Berrugate pour la première fois (Turonien moyen-Santonien et Coniacien-Santonien). Ces âges permettent d'établir une activité volcanique d'arc le long de la marge occidentale de la futur Plaque Caraïbes au moins depuis le Santonien et qui pourrait avoir durée jusqu'au Turonien moyen-Campanien inférieur. De plus, sur la base de ces âges à radiolaires, la Formation de Loma Chumico d'âge Albien, et la Formation de Berrugate d'âge Turonien moyen-Maastrichtien inférieur, peuvent maintenant être différenciées. Deux échantillons de la Formation de Sabana Grande donnent des âges Coniacien-Santonien et Coniacien-Campanien et indiquent qu'il existe une lacune stratigraphique d'environ 10 millions d'années entre cette formation et la Formation de Loma Chumico sous-jacente d'âge Albien.

Phosphorus and carbon burial during cretaceous oceanic anoxic events, links with climatic changes

According to Jenkyns (2010), oceanic anoxic events (OAE) record profound changes in the climatic and paleoceanographic state of the planet and represent major disturbances in the global carbon cycle. One of the most studied OAEs on a worldwide scale is the Cenomanian-Turonian OAE 2, which is characterized by a pronounced positive excursion in carbon-isotope records and the important accumulation of organic-rich sediments. The section at Gongzha (Tibet) and the sections at Barranca and Axaxacualco (Mexico) are located in remote parts of the Tethys, and show δ13C records, which are well correlated with those of classical Tethyan sections. Both sections, however, do not exhibit the presence of organic-rich sediments. Phosphorus Mass Accumulation Rates (PMAR) in Tibet show a pattern similar to that observed in the Tethys by Mort et al. (2007), which suggests enhanced Ρ regeneration during the OAE 2 time interval, though there is no evidence for anoxic conditions in Tibet. Ρ appears here to have been mainly driven by detrital influx and sea-level fluctuations. The sections at Barranca and Axaxacualco show that the Mexican carbonate platform persisted during this anoxic event, which allowed the evolution of platform fauna otherwise not present in Tethyan sections. The persistence of this carbonate platform close to the Caribbean Igneous Plateau, which is thought to have released bio-limiting metals, is explained by local uplift which delayed the drowning of the platform and a specific oceanic circulation that permitted the preservation of oligotrophic conditions in the area. The Coniacian-Santonian OAE (OAE3) appears to have been more dependent on local conditions than OAE2. The presence of black shales associated with OAE3 appear to have been restricted to shallow-water settings and epicontinental seas in areas located around the Atlantic Ocean. The sections at Olazagutia (Spain), and Ten Mile - Arbor Park (USA), two potential Global Boundary Stratotype Sections and Points (GSSP) sites, are devoid of organic-rich sediments and lack a δ13C positive excursion around the C-S boundary. The Gabal Ekma section (Sinai, Egypt) exhibits accumulations of organic-rich sediments, in addition to phosphorite bone beds layers, which may have been linked to an epicontinental upwelling zone and/or storm inputs. Our data suggest that OAE 3 is rarely expressed by truly anoxic conditions and seems to have been linked to local conditions rather than global paleoenvironmental change. The evidence for detrital-P being the likely cause of Ρ fluctuations during the OAEs studied here does not negate the idea that anoxia was the principal driver of these fluctuations in the western Tethys. However, an explanation is required as to why the Ρ accumulation signatures are mirrored in both oxic and anoxic sedimentary successions. 'Eustatic/climatic' and 'productivity/anoxic' models may have both operated simultaneously in different parts of the world depending on local conditions, both producing similar trends in Ρ accumulation. - Selon Jenkyns (2010), les événements anoxiques océaniques enregistrent de profonds changements dans le climat et la paléoceanographie de la planète et représente des perturbations majeures du cycle du carbone. L'un des plus étudiés à l'échelle mondiale est l'ΟΑΕ2 du Cénomanien-Turonien, qui est caractérisé par une très forte excursion positive des isotopes du carbone et une importante accumulation de sédiments riche en matière organique. La section de Gongzha (Tibet) et les sections de Barranca et Axaxcualco (Mexique) sont situées aux confins de la Téthys, et enregistrent une courbe isotopique en δ13C parfaitement corrélable avec les sections téthysiennes, mais ne montre pas d'accumulation de black shales. Le taux de phosphore en accumulation de masses (PMAR) au Tibet montre un pattern similaire observé également par Mort et al. (2007) dans la Téthys, suggérant un model de régénération du Ρ durant l'anoxie, cependant aucune conditions anoxiques régnent dans la région du Tibet. Ρ apparaît donc principalement guidé par le détritisme et les fluctuations du niveau marin. Les sections de Barranca et d'Axaxacualco montrent que la plateforme carbonatée mexicaine persiste durant cet événement anoxique, et permet le développement d'une faune de plateforme qui n'est pas présente dans les sections téthysiennes. La persistance de cette plateforme carbonatée si proche du plateau Caribéen, qui est connu pour le relâchement de métaux bio-limitant, peut être expliqué par un soulèvement tectonique local qui inhibe l'ennoiement de la plateforme et une circulation océanique spécifique qui permet la préservation de conditions oligotrophiques dans cette région. L'événement anoxique océanique du Coniacien-Santonien apparaît plus dépendant des conditions locales que pour l'ΟΑΕ2. Les black shales associés à POAE3 sont restreints aux zones situées autour de l'océan Atlantique et plus particulièrement aux eaux peu profondes et épicontinentales. Les sections d'Olazagutia (Espagne), Ten Mile Creek et Arbor Park (USA), qui sont deux potentielles sections GSSP (Sections de stratotype de limite globaux et de points), ne montre pas d'accumulation de black shales et pas de forte excursion positive en δ13C autour de la limite C-S. La section de Gabal Ekma (Sinai, Egypte) montre des accumulations de black shales, en plus des couches de phosphorites et d'accumulation d'os (« bone beds »), vraisemblablement lié à des zones active d'upwelling épicontinentale et/ou d'apport de tempêtes. Nos données suggèrent que l'OAE 3 est rarement exprimé par de vraies conditions anoxiques et semble être plus lié à des conditions plus locales que des changements paléo-environnementaux globaux, comme observés pour le Cénomanien- Turonien. Les arguments pour un modèle lié au phosphore détritique qui serait la cause des fluctuations du phosphore total durant les OAEs, n'écartent pas l'idée que l'anoxie est la principale cause de ces fluctuations dans les sections riches en matière organique de l'Ouest téthysien. Cependant une explication est nécessaire pour comprendre pourquoi la signature de l'accumulation du phosphore est semblable dans les successions sédimentaires déposées dans des conditions oxygénées et anoxiques. Les modèles « Eustatisme/Climat » et « Productivité/anoxie » ont simultanément opéré dans les différentes parties du monde dépendant de conditions locales, et ont produit des tendances similaires en accumulation de phosphore.

Limited sampling strategies for monitoring tacrolimus in paediatric liver transplant recipients

Ce travail de recherche a été réalisé dans le laboratoire de pharmacologie clinique, au Centre Hospitalier Universitaire Sainte-Justine, à Montréal. C'est une étude rétrospective basée sur le suivi thérapeutique du Tacrolimus prescrit chez les enfants après transplantation hépatique. Ce suivi est nécessaire car le Tacrolimus possède une importante variabilité pharmacocinétique inter et intra-individuelle ainsi qu'un index thérapeutique très étroit. Actuellement, l'individualisation des doses prescrites est basée sur la mesure de la concentration de base - du médicament dans le sang (C0), mais des études récentes montrent que cette mesure ne reflète pas précisément l'exposition du Tacrolimus dans l'organisme chez les enfants. Le meilleur reflet de cette exposition est la mesure de l'aire sous la courbe (AUC). Cependant, cette dernière implique la mesure de multiples concentrations tout au long de l'intervalle entre 2 doses de médicament (Tacrolimus: 12 heures) ce qui est long, cher et impraticable en ambulatoire. De nouvelles méthodes utilisant un nombre limité de prélèvements ont donc été développées pour prédire au mieux cette AUC. Ce sont les "Limited sampling strategies" ou LSS. La plupart de ces LSS pour le Tacrolimus ont été développées et validées chez des patients transplantés adultes et leur application directe chez les transplantés pédiatriques n'est pas possible en raison de différences importantes au niveau des paramètres pharmacocinétiques du médicament entre ces deux populations. Aussi, le but de ce travail était de développer et valider, pour la première fois, des LSS chez les enfants transplantés hépatiques. Pour cela, une analyse de 36 profils pharmacocinétiques de 28 patients transplantés hépatiques âgés de 0.4- 18.5 ans a été effectuée. Tous les profils ont été réalisés au Centre Hospitalier Universitaire Sainte-Justine entre janvier 2007 et janvier 2009. Les LSS comportant au maximum 4 mesures de concentration ont été développées en utilisant une analyse de régression multiple. Parmi tous les modèles obtenus, cinq ont été sélectionnés sur la base de critères précis puis validés selon la méthode décrite par Sheiner et Beal.¦Les résultats montrent que ces cinq modèles peuvent prédire l'AUC du Tacrolimus avec une précision cliniquement acceptable de ± 15% alors que la C0 présente la plus faible corrélation avec l'AUC.¦En conclusion, cette étude confirme que la C0 ne permet pas de prédire de manière efficace l'exposition du Tacrolimus dans l'organisme dans notre population de patients pédiatriques contrairement aux LSS analysées qui offrent une méthode pratique et fiable. Par ailleurs, en permettant d'obtenir une estimation précise et simplifiée de l'AUC complète du Tacrolimus chez les patients, ces LSS ouvrent la porte à de futures études prospectives visant à mieux définir l'AUC cible du médicament et à déterminer si le suivi basé sur la mesure de l'AUC est plus efficace et plus sûr que celui basé sur la mesure de la C0.

Vers un naturalisme social; À la croisée des sciences sociales et des sciences cognitives

Pour les sociologues, l'explication de l'action ne peut reposer sur la nature première de l'esprit (les mécanismes universels de détection et de traitement de l'information dont les sciences cognitives prétendent rendre compte) mais sur sa nature seconde (les raisons d'agir plus ou moins collectives qui le meublent). Même les paradigmes sociologiques les plus opposés s'accordent pour traiter le social et le culturel comme le résultat d'un patient travail de dénaturalisation qui creuse un fossé incompressible entre l'ordre réducteur des causes naturelles et l'ordre respectable, de facture mentale ou sociale, des raisons, entre les instincts de base de l'organisme et la logique de haut-niveau des significations partagées et des normes culturelles. Or, l'objectif du naturalisme est précisément le dépassement de la rupture interne entre l'humain et l'animal, entre l'histoire culturelle et l'histoire biologique. Il vise à harmoniser, dans le sens de "rendre compatible", les hypothèses et les résultats des sciences sociales avec ceux des sciences naturelles, entreprise qui incite la sociologie à expliciter et à réviser les modèles cognitifs et anthropologiques auxquels elle fait implicitement appel. Un des apports potentiels des sciences cognitives et, plus généralement, du naturalisme est donc d'améliorer, de remplacer ou de falsifier, au sens poppérien du terme, les modèles cognitifs et les conceptions anthropologiques sur lesquels se basent les sciences de la société.

L'intersubjectivité en questions : agrégat ou principe organisateur en psychologie ?

La notion d'intersubjectivité est au coeur des débats contemporains dans le champ de la psychologie. L'intersubjectivité traverse en effet les différentes approches de la psychologie, des neurosciences à la psychologie sociale, des perspectives développementales aux approches psychoaffectives de la personnalité, des contextes de l'apprentissage à ceux du soin. Cet ouvrage collectif, premier volume de la collection Actualités psychologiques, se propose de faire le point sur cette notion et d'apporter des pistes actuelles de compréhension aux enjeux de la place de l'autre et de l'altérité dans la construction du sujet. Des chercheurs de renommée internationale, de Suisse et de différents pays européens, apportent ici leur contribution, dans différents domaines de la psychologie, invitant au dialogue entre les modèles théoriques qui sous-tendent la notion d'intersubjectivité et entre les méthodologies qui permettent d'en appréhender la dynamique. Contribution exceptionnelle, cet ouvrage est à destination de tous, psychologues praticiens et professionnels de la relation, chercheurs et étudiants en psychologie ou dans des disciplines connexes. Il permet d'approcher l'intersubjectivité comme une notion tout à la fois familière et complexe, au service de la rencontre de l'autre.

Regulation of frizzled receptor activity at the plasma membrane : an interplay of the receptor, the trimeric G protein Go, and RGS protein and a scaffolding protein Kermit

G-protein-signaling pathways convey extracellular signals inside the cells and regulate distinct physiological responses. This type of signaling pathways consists of three major components: G-protein-coupled receptors (GPCRs), heterotrimeric G proteins (G-proteins) and downstream effectors. Upon ligand binding, GPCRs activate heterotrimeric G proteins to initiate the signaling cascade. Dysfunction of GPCR signaling correlates with numerous diseases such as diabetes, nervous and immune system deficiency, and cancer. As the signaling switcher, G-proteins (Gs, Gq/11, G12/13, and Gi/o) have been an appealing topic of research for decades. A heterotrimeric G-protein is composed of three subunits, the guanine nucleotide associated a-subunit, ß and y subunits. In general, the duration of signaling is determined by the lifetime of activated (GTP bound) Ga subunits. Identification of novel communication partners of Ga subunits appears to be an attractive way to understand the machinery of GPCR signaling. In our lab, we mainly focus on Gao, which is abundantly expressed in the nervous system. Here we present two novel interacting partners of Drosophila Gao: Dhit and Kermit, identified through yeast two-hybrid screening and genetic screening respectively. Dhit is characterized by a small size with a conserved RGS domain and an N-terminal cysteine rich motif. The RGS domain possesses the GAP (GTPase activating protein) activity towards G proteins. However, we found that Dhit exerts not only the GAP activity but also the GDI (guanine nucleotide dissociation inhibitor) activity towards Gao. The unexpected GDI activity is preserved in GAIP/RGS19 - a mammalian homologue of Dhit. Further experiments confirmed the GDI activity of Dhit and GAIP/RGS19 in Drosophila and mammalian cell models. Therefore, we propose that Dhit and its mammalian homologues modulate GPCR signaling by a double suppression of Ga subunits - suppression of their nucleotide exchange with GTP and acceleration of their hydrolysis of GTP. Kermit/GEPC was first identified as a binding partner of GAIP/RGS19 in a yeast two- hybrid screen. Instead of interacting with the Drosophila homologue of GAIP/RGS19 (Dhit), Kermit binds to Gao in vivo and in vitro. The functional consequence of Kermit/Gao interaction is the regulation of localization of Vang (one of the planar cell polarity core components) at the apical membrane. Overall, my work elaborated the action of Gao with its two interaction partners in Gao- mediated signaling pathway. Conceivably, the understanding of GPCR signaling including Gao and its regulators or effectors will ultimately shed light on future pharmaceutical research. - Les voies de signalisation médiées par les protéines G transmettent des signaux extracellulaires à l'intérieur des cellules pour réguler des réponses physiologiques distinctes. Cette voie de signalisation consiste en trois composants majeurs : les récepteurs couplés aux protéines G (GPCRs), les protéines G hétérotrimériques (G-proteins) et les effecteurs en aval. Suite à la liaison du ligand, les GPCRs activent les protéines G hétérotrimériques qui initient la cascade de signalisation. Des dysfonctions dans la signalisation médiée par les GPCRs sont corrélées avec de nombreuses maladies comme le diabète, des déficiences immunes et nerveuses, ainsi que le cancer. Puisque la voie de signalisation s'active et se désactive, les protéines G (Gs, Gq/11, G12/13 et Gi/o) ont été un sujet de recherche attrayant pendant des décennies. Une protéine G hétérotrimérique est composée de trois sous-unités, la sous-unité a associée au nucléotide guanine, ainsi que les sous-unités ß et y. En général, la durée du signal est déterminée par le temps de demi-vie des sous-unités Ga activées (Ga liées au GTP). Identifier de nouveaux partenaires de communication des sous-unités Ga se révèle être un moyen attractif de comprendre la machinerie de la signalisation par les GPCRs. Dans notre laboratoire nous nous sommes concentrés principalement sur Gao qui est exprimée de manière abondante dans le système nerveux. Nous présentons ici deux nouveaux partenaires qui interagissent avec Gao chez la drosophile: Dhit et Kermit, qui ont été identifiés respectivement par la méthode du yeast two-hybrid et par criblage génétique. Dhit est caractérisé par une petite taille, avec un domaine RGS conservé et un motif N- terminal riche en cystéines. Le domaine RGS contient une activité GAP (GTPase activating protein) pour les protéines G. Toutefois, nous avons découvert que Dhit exerce non seulement une activité GAP mais aussi une activité GDI (guanine nucleotide dissociation inhibitor) à l'égard de Gao. Cette activité GDI inattendue est préservée dans RGS19 - un homologue de Dhit chez les mammifères. Des expériences supplémentaires ont confirmé l'activité GDI de Dhit et de RGS19 chez Drosophila melanogaster et les modèles cellulaires mammifères. Par conséquent, nous proposons que Dhit et ses homologues mammifères modulent la signalisation GPCR par une double suppression des sous-unités Ga - suppression de leur nucléotide d'échange avec le GTP et une accélération dans leur hydrolyse du GTP. Kermit/GIPC a été premièrement identifié comme un partenaire de liaison de RGS19 dans le criblage par yeast two-hybrid. Au lieu d'interagir avec l'homologue chez la drosophile de RGS19 (Dhit), Kermit se lie à Gao in vivo et in vitro. La conséquence fonctionnelle de l'interaction Kermit/Gao est la régulation de la localisation de Vang, un des composants essentiel de la polarité planaire cellulaire, à la membrane apicale. Globalement, mon travail a démontré l'action de Gao avec ses deux partenaires d'interaction dans la voie de signalisation médiée par Gao. La compréhension de la signalisation par les GPCRs incluant Gao et ses régulateurs ou effecteurs aboutira à mettre en lumière de futurs axes dans la recherche pharmacologique.

Late Cretaceous to Eocene geology of the South Central American forearc area (southern Costa Rica and western Panama) : initiation and evolution of an intra-oceanic convergent margin

Résumé : L'arc volcanique du sud de l'Amérique Centrale se situe sur la marge SW de la Plaque Caraïbe, au-dessus des plaques subduites de Cocos et Nazca. Il s'agit de l'un des arcs intra-océaniques les plus étudiés au monde, qui est généralement considéré comme s'étant développé à la fin du Crétacé le long d'un plateau océanique (le Plateau Caraïbe ou CLIP) et se trouvant actuellement dans un régime de subduction érosive. Au cours des dernières décennies, des efforts particuliers ont été faits pour comprendre les processus liés à la subduction sur la base d'études géophysiques et géochimiques. Au sud du Costa Rica et à l'ouest du Panama, des complexes d'accrétions et structures à la base de l'arc volcanique ont été exposés grâce à la subduction de rides asismiques et de failles transformantes. Des affleurements, situés jusqu'à seulement 15 km de la fosse, offrent une possibilité unique de mieux comprendre quelques uns des processus ayant lieu le long de la zone de subduction. Nous présentons de nouvelles contraintes sur l'origine de ces affleurements en alliant une étude de terrain poussée, de nouvelles données géochimiques, sédimentaires et paléontologiques, ainsi que des observations structurales effectuées en télédétection. Une nouvelle stratigraphie tectonique entre le Campanien et l'Éocène est définie pour la région d'avant-arc située entre la Péninsule d'Osa (Costa Rica) et la Péninsule d'Azuero (Panama). Nos résultats montrent que la partie externe de la marge est composée d'un arrangement complexe de roches ignées et de séquences sédimentaires de recouvrement qui comprennent principalement le socle de l'arc, des roches d'arc primitif, des fragments de monts sous-marins accrétés et des mélanges d'accrétion. Des preuves sont données pour le développement de l'arc volcanique du sud de l'Amérique Centrale sur un plateau océanique. Le début de la subduction le long de la marge SW de la Plaque Caraïbe a eu lieu au Campanien et a généré des roches d'arc primitif caractérisées par des affinités géochimiques particulières, globalement intermédiaires entre des affinités de plateau et d'arc insulaire. L'arc était mature au Maastrichtien et formait un isthme essentiellement continu entre l'Amérique du Nord et l'Amérique du Sud. Ceci a permis la migration de faunes terrestres entre les Amériques et pourrait avoir contribué à la crise fin Crétacé -Tertiaire en réduisant les courants océaniques subéquatoriaux entre le Pacifique et l'Atlantique. Plusieurs unités composées de fragments de monts sous-marins accrétés sont définies. La nature et l'arrangement structural de ces unités définissent de nouvelles contraintes sur les modes d'accrétion des monts sous-marins/îles océaniques et sur l'évolution de la marge depuis la formation de la zone de subduction. Entre la fin du Crétacé et l'Éocène moyen, la marge a enregistré plusieurs épisodes ponctuels d'accrétion de monts sous-marins alternant avec de la subduction érosive. A l'Éocène moyen, un événement tectonique régional pourrait avoir causé un fort couplage entre les plaques supérieure et inférieure, menant à des taux plus important d'accrétion de monts sous-marins. Durant cette période, la situation le long de la marge était très semblable à la situation actuelle et caractérisée par la présence de monts sous-marins subductants et l'absence d'accrétion de sédiments. L'enregistrement géologique montre qu'il n'est pas possible d'attribuer une nature érosive ou accrétionnaire à la marge dans le passé ou -par analogie- aujourd'hui, parce que (1) les processus d'accrétion et érosifs varient fortement spatialement et temporellement et (2) il est impossible d'évaluer la quantité exacte de matériel tectoniquement enlevé à la marge depuis le début de la subduction. Au sud du Costa Rica, certains fragments de monts sous-marins accrétés sont représentatifs d'une interaction entre une ride et un point chaud dans le Pacifique au Crétacé terminal/Paléocène. L'existence de ces fragments de monts sous-marins et la morphologie du fond de l'Océan Pacifique indiquent que la formation de la ride de Cocos-Nazca s'est formée au moins ~40 Ma avant l'âge proposé par les modèles tectoniques actuels. Au Panama, nous avons identifié une île océanique d'âge début Éocène qui a été accrétée à l'Éocène moyen. L'accrétion a eu lieu à très faible profondeur par détachement de l'île dans la fosse, et a mené à une exceptionnelle préservation des structures volcaniques. Des affleurement comprenant aussi bien des parties basses et hautes de l'édifice volcanique on été étudiées, depuis la phase sous-marine bouclier jusqu'à la phase subaérienne post-bouclier. La stratigraphie nous a permis de différencier les laves de la phase sous-marine de celles de la phase subaérienne. La composition des laves indique une diminution progressive de l'intensité de la fusion partielle de la source et une diminution de la température des laves produites durant les derniers stades de l'activité volcanique. Nous interprétons ces changements comme étant liés à l'éloignement progressif de l'île océanique de la zone de fusion ou point chaud. Abstract The southern Central American volcanic front lies on the SW edge of the Caribbean Plate, inboard of the subducting Cocos and Nazca Plates. It is one of the most studied intra-oceanic convergent margins around the world, which is generally interpreted to have developed in the late Cretaceous along an oceanic plateau (the Caribbean Large Igneous Province or CLIP) and to be currently undergoing a regime of subduction erosion. In the last decades a particular effort has been made to understand subduction-related processes on the basis of geophysical and geochemical studies. In southern Costa Rica and western Panama accretionary complexes and structures at the base of the volcanic front have been exposed in response to subduction of aseismic ridges and transforms. Onland exposures are located as close as to 15 km from the trench and provide a unique opportunity to better understand some of the processes occurring along the subduction zone. We provide new constraints on the origins of these exposures by integrating a comprehensive field work, new geochemical, sedimentary and paleontological data, as well as structural observations based on remote imaging. A new Campanian to Eocene tectonostratigraphy is defined for the forearc area located between the Osa Peninsula (Costa Rica) and the Azuero Peninsula (Panama). Our results show that the outer margin is composed of a complicated arrangement of igneous complexes and overlapping sedimentary sequences that essentially comprise an arc basement, primitive island-arc rocks, accreted seamount fragments and accretionary mélanges. Evidences are provided for the development of the southern Central American arc on the top an oceanic plateau. The subduction initiation along the SW edge of the Caribbean Plate occurred in the Campanian and led to formation of primitive island-arc rocks characterized by unusual geochemical affinities broadly intermediate between plateau and arc affinities. The arc was mature in the Maastrichtian and was forming a predominantly continuous landbridge between the North and South Americas. This allowed migration of terrestrial fauna between the Americas and may have contributed to the Cretaceous-Tertiary crisis by limiting trans-equatorial oceanic currents between the Pacific and the Atlantic. Several units composed of accreted seamount fragments are defined. The nature of the units and their structural arrangement provide new constraints on the modes of accretion of seamounts/oceanic islands and on the evolution of the margin since subduction initiation. Between the late Cretaceous and the middle Eocene, the margin recorded several local episodes of seamount accretion alternating with tectonic erosion. In the middle Eocene a regional tectonic event may have triggered strong coupling between the overriding and subducting plates, leading to higher rates of seamount accretion. During this period the situation along the margin was very similar to the present and characterized by subducting seamounts and absence of sediment accretion. The geological record shows that it is not possible to ascribe an overall erosive or accretionary nature to the margin in the past and, by analogy, today, because (1) accretionary and erosive processes exhibit significant lateral and temporal variations and (2) it is impossible to estimate the exact amount of material tectonically eroded from the margin since subduction initiation. In southern Costa Rica, accreted seamount fragments point toward a plume-ridge interaction in the Pacific in the late Cretaceous/Paleocene. This occurrence of accreted seamount fragments and morphology of the Pacific Ocean floor is indicative of the formation of the Cocos-Nazca spreading system at least ~40 Ma prior to the age proposed in current tectonic models. In Panama, we identified a remarkably-well preserved early Eocene oceanic island that accreted in the middle Eocene. The accretion probably occurred at very shallow depth by detachment of the island in the trench and led to an exceptional preservation of the volcanic structures. Exposures of both deep and superficial parts of the volcanic edifice have been studied, from the submarine-shield to subaerial-postshield stages. The stratigraphy allowed us to distinguish lavas produced during the submarine and subaerial stages. The lava compositions likely define a progressive diminution of source melting and a decrease in the temperature of erupted melts in the latest stages of volcanic activity. We interpret these changes to primarily reflect the progressive migration of the oceanic island out of the melting region or hotspot.