Magma generation at the easternmost section of the Hellenic arc: Hf, Nd, Pb and Sr isotope geochemistry of Nisyros and Yali volcanoes (Greece)

Geochemical and petrographical studies of lavas and ignimbrites from the Quaternary Nisyros-Yali volcanic system in the easternmost part of the Hellenic arc (Greece) reveal insight into magma generating processes. A compositional gap between 61 and 68 wt.% SiO2 is recognized that coincides with the stratigraphic distinction between pre-caldera and postcaldera volcanic units. Trace element systematics support the subdivision of Nisyros and Yali volcanic units into two distinct suites of rocks. The variation of Nd and Hf present day isotope data and the fact that they are distinct from the isotope compositions of MORB rule out an origin by pure differentiation and require assimilation of a crustal component. Lead isotope ratios of Nisyros and Yali volcanic rocks support mixing of mantle material with a lower crust equivalent. However, Sr-87/Sr-86 ratios of 0.7036-0.7048 are incompatible with a simple binary mixing scenario and give low depleted mantle extraction ages (< 0.1 Ga), in contrast with Pb model ages of 0.3 Ga and Hf and Nd model ages of ca. 0.8 Ga. The budget of fluid-mobile elements Sr and Pb is likely to be dominated by abundant hydrous fluids characterised by mantle-like Sr isotope ratios. Late stage fluids probably were enriched in CO2, needed to explain the high Th concentrations. The occurrence of hydrated minerals (e.g., amphibole) in the first post-caldera unit with the lowermost Sr-87/Sr-86 ratio of 0.7036 +/- 2 can be interpreted as the result of the increased water activity in the source. The presence of two different plagioclase phenocryst generations in the first lava subsequent to the caldera-causing event is indicative for a longer storage time of this magma at a shallower level. A model capable of explaining these observations involves three evolutionary stages. First stage, assimilation of lower crustal material by a primitive magma of mantle origin (as modelled by Nd-Hf isotope systematics). This stage ended by an interruption in replenishment that led to an increase of crystallization and, hence, an increase in viscosity, suppressing eruption. During this time gap, differentiation by fractional crystallization led to enrichment of incompatible species, especially aqueous fluids, to silica depolymerisation and to a decrease in viscosity, finally enabling eruption again in the third stage. (c) 2005 Elsevier B.V. All rights reserved.

Allostatic working memory processing in schizophrenia: Human and animal model coherence

The most evident symptoms of schizophrenia are severe impairment of cognitive functions like attention, abstract reasoning and working memory. The latter has been defined as the ability to maintain and manipulate on-line a limited amount of information. Whereas several studies show that working memory processes are impaired in schizophrenia, the specificity of this deficit is still unclear. Results obtained with a new paradigm, involving visuospatial, dynamic and static working memory processing, suggest that schizophrenic patients rely on a specific compensatory strategy. An animal model of schizophrenia with a transient deficit in glutathione during the development reveals similar substitutive processing, masking the impairment in working memory functions in specific test conditions only. Taken together, these results show coherence between working memory deficits in schizophrenic patients and in animal models. More generally, it is possible to consider that the pathological state may be interpreted as a reduced homeostatic reserve. However, this may be balanced in specific situations by efficient allostatic strategies. Thus, the pathological condition would remain latent in several situations, due to such allostatic regulations. However, to maintain a performance based on highly specific strategies requires in turn specific conditions, limitating adaptative resources in humans and in animals. In summary, we suggest that the psychological and physical load to maintain this rigid allostatic state is very high in patients and animal subjects.

Structure et évolution thermique de la croûte supérieure des chaînes de subduction : approches thermologique et modélisation numérique, dans les North Cascades (U.S.A) et de la zone de faille Motagua (Guatemala)

Cette thèse cible l'étude de la structure thermique de la croûte supérieure (<10km) dans les arcs magmatiques continentaux, et son influence sur l'enregistrement thermochronologique de leur exhumation et de leur évolution topographique. Nous portons notre regard sur deux chaînes de montagne appartenant aux Cordillères Américaines : Les Cascades Nord (USA) et la zone de faille Motagua (Guatemala). L'approche utilisée est axée sur la thermochronologie (U-Th-Sm)/He sur apatite et zircon, couplée avec la modélisation numérique de la structure thermique de la croûte. Nous mettons en évidence la variabilité à la fois spatiale et temporelle du gradient géothermique, et attirons l'attention du lecteur sur l'importance de prendre en compte la multitude des processus géologiques perturbant la structure thermique dans les chaînes de type cordillère, c'est à dire formées lors de la subduction océanique sous un continent.Une nouvelle approche est ainsi développée pour étudier et contraindre la perturbation thermique autour des chambres magmatiques. Deux profiles âge-elevation (U-Th-Sm)/He sur apatite et zircon, ont été collectées 7 km au sud du batholithe de Chilliwack, Cascades Nord. Les résultats montrent une variabilité spatiale et temporelle du gradient géothermique lors de l'emplacement magmatique qui peut être contrainte et séparé de l'exhumation. Durant l'emplacement de l'intrusion, la perturbation thermique y atteint un état d'équilibre (-80-100 °C/km) qui est fonction du flux de magma et de ia distance à la source du magma, puis rejoint 40 °C/km à la fin du processus d'emplacement magmatique.Quelques nouvelles données (U-Th)/He, replacées dans une compilation des données existantes dans les Cascades Nord, indiquent une vitesse d'exhumation constante (-100 m/Ma) dans le temps et l'espace entre 35 Ma et 2 Ma, associée à un soulèvement uniforme de la chaîne contrôlé par l'emplacement de magma dans la croûte durant toute l'activité de l'arc. Par contre, après ~2 Ma, le versant humide de la chaîne est affecté par une accélération des taux d'exhumation, jusqu'à 3 km de croûte y sont érodés. Les glaciations ont un triple effet sur l'érosion de cette chaîne: (1) augmentation des vitesses d'érosion, d'exhumation et de soulèvement la où les précipitations sont suffisantes, (2) limitation de l'altitude contrôlé par la position de Γ Ε LA, (3) élargissement du versant humide et contraction du versant aride de la chaîne.Les modifications des réseaux de drainage sont des processus de surface souvent sous-estimés au profil d'événements climatiques ou tectoniques. Nous proposons une nouvelle approche couplant une analyse géomorphologique, des données thermochronologiques de basse température ((U-Th-Sm)/He sur apatite et zircon), et l'utilisation de modélisation numérique thermo-cinématique pour les mettre en évidence et les dater; nous testons cette approche sur la gorge de la Skagit river dans les North Cascades.De nouvelles données (U-Th)/He sur zircons, complétant les données existantes, montrent que le déplacement horizontal le long de la faille transformante continentale Motagua, la limite des plaques Caraïbe/Amérique du Nord, a juxtaposé un bloc froid, le bloc Maya (s.s.), contre un bloque chaud, le bloc Chortis (s.s.) originellement en position d'arc. En plus de donner des gammes d'âges thermochronologiques très différents des deux côtés de la faille, le déplacement horizontal rapide (~2 cm/a) a produit un fort échange thermique latéral, résultant en un réchauffement du côté froid et un refroidissement du côté chaud de la zone de faille de Motagua.Enfin des données (U-Th-Sm)/He sur apatite témoignent d'un refroidissement Oligocène enregistré uniquement dans la croûte supérieure de la bordure nord de la zone de faille Motagua. Nous tenterons ultérieurement de reproduire ce découplage vertical de la structure thermique par la modélisation de la formation d'un bassin transtensif et de circulation de fluides le long de la faille de Motagua. - This thesis focuses on the influence of the dynamic thermal structure of the upper crust (<10km) on the thermochronologic record of the exhumational and topographic history of magmatic continental arcs. Two mountain belts from the American Cordillera are studied: the North Cascades (USA) and the Motagua fault zone (Guatemala). I use a combined approach coupling apatite and zircon (U-Th-Sm}/He thermochronology and thermo- kinematic numerical modelling. This study highlights the temporal and spatial variability of the geothermal gradient and the importance to take into account the different geological processes that perturb the thermal structure of Cordilleran-type mountain belts (i.e. mountain belts related to oceanic subduction underneath a continent}.We integrate apatite and zircon (U-Th)/He data with numerical thermo-kinematic models to study the relative effects of magmatic and surface processes on the thermal evolution of the crust and cooling patterns in the Cenozoic North Cascades arc (Washington State, USA). Two age-elevation profiles that are located 7 km south of the well-studied Chiliiwack intrusions shows that spatial and temporal variability in geothermal gradients linked to magma emplacement can be contrained and separated from exhumation processes. During Chiliiwack batholith emplacement at -35-20 Ma, the geothermal gradient of the country rocks increased to a very high steady-state value (80-100°C/km), which is likely a function of magma flux and the distance from the magma source area. Including temporally varying geothermal gradients in the analysis allows quantifying the thermal perturbation around magmatic intrusions and retrieving a relatively simple denudation history from the data.The synthesis of new and previously published (U-Th)/He data reveals that denudation of the Northern Cascades is spatially and temporally constant at -100 m/Ma between ~32 and ~2 Ma, which likely reflects uplift due to magmatic crustal thickening since the initiation of the Cenozoic stage of the continental magmatic arc. In contrast, the humid flank of the North Cascades is affected by a ten-fold acceleration in exhumation rate at ~2 Ma, which we interpret as forced by the initiation of glaciations; around 3 km of crust have been eroded since that time. Glaciations have three distinct effects on the dynamics of this mountain range: (1) they increase erosion, exhumation and uplift rates where precipitation rates are sufficient to drive efficient glacial erosion; (2) they efficiently limit the elevation of the range; (3) they lead to widening of the humid flank and contraction of the arid flank of the belt.Drainage reorganizations constitute an important agent of landscape evolution that is often underestimated to the benefit of tectonic or climatic events. We propose a new method that integrates geomorphology, low-temperature thermochronometry (apatite and zircon {U-Th-Sm)/He), and 3D numerical thermal-kinematic modelling to detect and date drainage instability producing recent gorge incision, and apply this approach to the Skagit River Gorge, North Cascades.Two zircon (U-Th)/He age-elevation profiles sampled on both sides of the Motagua Fault Zone (MFZ), the boundary between the North American and the Caribbean plates, combined with published thermochronological data show that strike-slip displacement has juxtaposed the cold Maya block (s.s.) against the hot, arc derived, Chortis block (s.s ), producing different age patterns on both sides of the fault and short-wavelength lateral thermal exchange, resulting in recent heating of the cool side and cooling of the hot side of the MFZ.Finally, an apatite (U-Th-Sm)/He age-elevation profile records rapid cooling at -35 Ma localized only in the upper crust along the northern side of the Motagua fault zone. We will try to reproduce these data by modeling the thermal perturbation resulting from the formation of a transtensional basin and of fluid flow activity along a crustal- scale strike-slip fault.

Geochemistry and stable isotope composition of fresh alkaline porphyry copper tailings: Implications on sources and mobility of elements during transport and early stages of deposition

Prevention of acid mine drainage (AMD) in sulfide-containing tailings requires the identification of the geochemical processes and element pathways in the early stages of tailing deposition. However, analyses of recently deposited tailings in active tailings impoundments are scarce because mineralogical changes occur near the detection limits of many assays. This study shows that a detailed geochemical study which includes stable isotopes of water (delta H-2, delta O-18), dissolved sulfates (delta S-34, delta O-18) and hydrochernical parameter (pH, Eh, DOC, major and trace elements) from tailings samples taken at different depths in rainy and dry seasons allows the understanding of weathering (oxidation, dissolution, sorption, and desorption), water and element pathways, and mixing processes in active tailings impoundments. Fresh alkaline tailings (pH 9.2-10.2) from the Cu-Mo porphyry deposit in El Teniente, Chile had low carbonate (0.8-1.1 Wt-% CaCO3 equivalent) and sulfide concentrations (0.8-1.3 wt.%, mainly as pyrite). In the alkaline tailings water, Mo and Cu (up to 3.9 mg/L Mo and 0.016 mg/L Cu) were mobile as MoO42- and Cu (OH)(2)(0). During the flotation, tailings water reached equilibrium with gypsum (up to 738 mg/L Ca and 1765 mg/ L SO4). The delta S-34 VS. delta O-18 covariations of dissolved sulfate (2.3 to 4.5% delta S-34 and 4.1 to 6.0 % delta O-18) revealed the sulfate sources: the dissolution of primary sulfates (12.0 to 13.2%. delta S-34, 7.4 to 10.9%.delta O-18) and oxidation of primary sulfides (-6.7 to 1.7%. delta S-34). Sedimented tailings in the tailings impoundment can be divided into three layers with different water sources, element pathways, and geochemical processes. The deeper sediments (> 1 m depth) were infiltrated by catchment water, which partly replaced the original tailings water, especially during the winter season. This may have resulted in the change from alkaline to near-neutral pH and towards lower concentrations of most dissolved elements. The neutral pH and high DOC (up to 99.4 mg/L C) of the catchment water mobilized Cu (up to 0.25 mg/L) due to formation of organic Cu complexes; and Zn (up to 130 mg/L) due to dissolution of Zn oxides and desorption). At I m depth, tailings pore water obtained during the winter season was chemically and isotopically similar to fresh tailings water (pH 9.8-10.6, 26.7-35.5 mg/L Cl, 2.3-6.0 mg/L Mo). During the summer, a vadose zone evolved locally and temporarily up to 1.2 m depth. resulting in a higher concentration of dissolved solids in the pore water due to evaporation. During periodical new deposition of fresh tailings, the geochemistry of the surface layer was geochemically similar to fresh tailings. In periods without deposition, sulfide oxidation was suggested by decreasing pH (7.7-9.5), enrichment of MoO42- and SO42-, and changes in the isotopic composition of dissolved sulfates. Further enrichment for Na, K, Cl, SO4, Mg, Cu, and Mo (up to 23.8 mg/L Mo) resulted from capillary transport towards the surface followed by evaporation and the precipitation of highly soluble efflorescent salts (e.g., mirabilite, syngenite) at the tailing surface during summer. (C) 2008 Elsevier B.V. All rights reserved.

Characterization of PHO1 and PHO1;H1 gene function in Arabidopsis thaliana

Summary Phosphorus is one of the major macronutrients required for plant growth and development. Plant roots acquire phosphorus as inorganic phosphate (Pi), which is further distributed to the shoot, via the transpiration stream and root pressure, where Pi is imported again into cells. PHO1 in Arabidopsis has been identified as a protein involved in the loading of Pi into the root xylem. PHO1 does not have any homology to described Pi transporters including the Pht1 family of H+/ Pi cotransporters. PHO1 bears two domains, SPX and EXS domains, previously identified in Saccharomyces cerevisiae proteins involved in Pi transport and/or sensing, or in sorting proteins to endomembranes. Phylogenetic analysis of the PHO1 gene family revealed the presence of three clusters, with PHO1 and PHO1;H1 forming one cluster. The biological significance behind this cluster was demonstrated by the complementation of the pho1 mutant with only PHO1 and PHO1;H1, of all the PHO1 family members, when expressed under the PHO1 promoter. PHO1 has been shown to be expressed mostly in the root vascular cylinder and at low level in the shoot. PHO1;H1 had a different expression pattern, being expressed in both root and shoot vascular cylinder to the same level, with the levels in leaves increasing with the leaf maturity, suggesting additional role of PHO1;H1 in the Pi mobilization in leaves. In order to further explore the role of PHO1, Pi dynamics was studied on plants expressing PHO1 at different levels compared to the wild type: PHO1 overexpressors, PHO1 underexpressors and the pho1 mutant. Overexpression of the PHO1 protein in the shoot vascular tissue was shown to lead to increased Pi efflux out of the leaf cells and Pi accumulation in the shoot xylem apoplast compared to wild type, confirming the hypothesized role of PHO1 in xylem loading with Pi. The overexpression of PHO1 in the shoot was responsible far both changed Pi dynamic and stunted growth of PHO1 overexpressors, as shown by grafting experiments between wild type and PHO1 overexpressor. We found a ca. 2 fold decrease of shoot phosphorus and a 5-10 fold decrease in vacuolar Pi content in the PHO1 underexpressors and the pho1 null mutant compared to wild type, consistent with the role of PHO1 in the transfer of Pi from the root to the shoot. Shoot Pi deficiency results in a poor growth of the pho1 mutant. Grafting experiments between pho1 and wild type confirmed that both Pi deficiency and stunt growth of the pho1 mutant were dependent on the pho1 root, further supporting the importance of PHO1 in the root xylem loading with Pi. The pho1 mutant and the PHO1 underexpressors accumulated 8-15 fold more Pi in the root relative to wild type. In contrast to the pho1 mutant, the growth of PHO1 underexpressors was not impaired by the low shoat Pi content. This finding suggests that either PHO1 protein or root Pi concentration is important in Pi signaling and development of Pi deficiency symptoms leading to reduced growth. Résumé Le phosphore est l'un des nutriments essentiels à la croissance et au développement des plantes. Les racines absorbent le phosphore sous forme de phosphate inorganique (Pi) qui est dirigé, par la transpiration et la pression de la racine, vers les feuilles où le phosphate est acquis par les cellules. La protéine PHO1 a été démontrée indispensable au chargement du Pi dans le xylème des racines d'Arabidopsis. PHO1 ne démontre pas d'homologie aux transporteurs de Pi connus, incluant la famille Pht1 de cotransporteurs H+/Pi qui ont comme fonction le transport du phosphate à l'intérieur de la cellule. PHO1 contient deux domaines, SPX et EXS, aussi présents dans des protéines de Saccharomyces cerevisiae impliquées dans le transport ou la perception du phosphate, ou dans la localisation des protéines vers différentes membranes. Le génome d'Arabidopsis contient onze gènes homologues à PHO1. Neuf de ces homologues sont répartis en trois groupes. PHO1 et PHO1;H1 forment un de ces groupes. Nos travaux ont démontré que seuls PHO1;H1 et PHO1, sous contrôle du promoteur PHO1, peuvent complémenter le mutant pho1. PHO1 est exprimé principalement dans le cylindre vasculaire de la racine et faiblement dans la partie aérienne. Le degré d'expression de PHO1;H1 est similaire dans le cylindre vasculaire de la racine et des feuilles. Ceci suggère que PHO1;H1 est aussi impliqué dans la mobilisation du Pi dans les feuilles, en plus de son rôle dans le transfert du Pi dans le xylème des racines. Afin de mieux explorer le rôle de PHO1, la dynamique du phosphate a été observée dans trois lignées de plantes transgéniques: un sur-expresseur de PHO1, un sous-expresseur de PHO1 et le mutant pho1. La sur-expression de PHO1 dans le tissue vasculaire des feuilles a provoqué l'efflux du Pi vers l'espace apoplastic du xylème, ce qui confirme le rôle de PHO1 dans le chargement du Pi dans le xylème. La sur-expressìon de PHO1 dans la rosette est responsable d'un changement de la dynamique du Pi et de la diminution de la croissance, ce qui fut démontré par une expérience de greffe de la rosette du sur-expresseur de PHO1 sur les racines du sauvage. On a observé pour le sous-expresseur de PHO1 et le mutant pho1 une diminution du phosphore d'environ 2 fais au niveau des feuilles, et une diminution de 5-10 fois du Pi dans les vacuoles des feuilles, par rapport au sauvage. Ceci confirme le rôle proposé de PHO1 dans le transfert du Pi des racines aux feuilles. La carence de Pi chez pho1 implique une diminution de la taille de la rosette. Pour expliquer ce phénotype une autre expérience de greffe démontra que la cause de ce changement provenait des racines. Ceci renforce l'hypothèse de l'importance du rôle de PHO1 dans le xylème de la racine pour le chargement du Pi. Le mutant phot et le sous-expresseur de PHO1 accumulent 8-15 fois plus de Pi dans leurs racines comparé au sauvage. Cependant, contrairement au phot mutant, le sous-expresseur de PHO1 avait une croissance comparable au sauvage malgré le niveau bas du Pi dans les feuilles. Ceci suggère que la taille de la rosette lors d'une carence en Pi chez Arabidopsis serait la conséquence d'un changement de concentration de Pi dans les racines ou d'une influence de la protéine PHO1.

Molecular and Functional Characterization of Neuroblastoma-Initiating Cells

Neuroblastoma (NB) is the most common extracranial malignant tumor in young children and arises at any site of the sympathetic nervous system. The disease exhibits a remarkable phenotypic diversity ranging from spontaneous regression to fatal disease. Poor outcome results from a rapidly progressive, metastatic and drug-resistant disease. Recent studies have suggested that solid tumors may arise from a minor population of cancer stem cells (CSCs) with stem cell markers and typical properties such as self-renewal ability, asymmetric division and drug resistance. In this model, CSCs possess the exclusive ability to initiate and maintain the tumor, and to produce distant metastases. Tumor cell subpopulations with stem-like phenotypes have indeed been identified in several cancer including leukemia, breast, brain and colon cancers. CSC hypothesis still needs to be validated in the other cancers including NB.NB originates from neural crest-derived malignant sympatho-adrenal cells. We have identified rare cells that express markers in conformity with neural crest stem cells and their derived lineages within primary NB tissue and cell lines, leading us to postulate the existence of CSCs in NB tumors.In the absence of specific markers to isolate CSCs, we adapted to NB tumor cells the sphere functional assay, based on the ability of stem cells to grow as spheres in non-adherent conditions. By serial passages of spheres from bone marrow NB metastases, a subset of cells was gradually selected and its specific gene expression profile identified by micro-array time-course analysis. The differentially expressed genes in spheres are enriched in genes implicated in development including CD133, ABC-transporters, WNT and NOTCH genes, identified in others solid cancers as CSCs markers, and other new markers, all referred by us as the Neurosphere Expression Profile (NEP). We confirmed the presence of a cell subpopulation expressing a combination of the NEP markers within a few primary NB samples.The tumorigenic potential of NB spheres was assayed by in vivo tumor growth analyses using orthotopic (adrenal glands) implantations of tumor cells into immune-compromised mice. Tumors derived from the sphere cells were significantly more frequent and were detected earlier compared to whole tumor cells. However, NB cells expressing the neurosphere-associated genes and isolated from the bulk tumors did not recapitulate the CSC-like phenotype in the orthotopic model. In addition, the NB sphere cells lost their higher tumorigenic potential when implanted in a subcutaneous heterotopic in vivo model.These results highlighted the complex behavior of CSC functions and led us to consider the stem-like NB cells as a dynamic and heterogeneous cell population influenced by microenvironment signals.Our approach identified for the first time candidate genes that may be associated with NB self-renewal and tumorigenicity and therefore would establish specific functional targets for more effective therapies in aggressive NB.

The tumor microenvironment and its contribution to tumor evolution toward metastasis.

Cancer cells acquire cell-autonomous capacities to undergo limitless proliferation and survival through the activation of oncogenes and inactivation of tumor suppressor genes. Nevertheless, the formation of a clinically relevant tumor requires support from the surrounding normal stroma, also referred to as the tumor microenvironment. Carcinoma-associated fibroblasts, leukocytes, bone marrow-derived cells, blood and lymphatic vascular endothelial cells present within the tumor microenvironment contribute to tumor progression. Recent evidence indicates that the microenvironment provides essential cues to the maintenance of cancer stem cells/cancer initiating cells and to promote the seeding of cancer cells at metastatic sites. Furthermore, inflammatory cells and immunomodulatory mediators present in the tumor microenvironment polarize host immune response toward specific phenotypes impacting tumor progression. A growing number of studies demonstrate a positive correlation between angiogenesis, carcinoma-associated fibroblasts, and inflammatory infiltrating cells and poor outcome, thereby emphasizing the clinical relevance of the tumor microenvironment to aggressive tumor progression. Thus, the dynamic and reciprocal interactions between tumor cells and cells of the tumor microenvironment orchestrate events critical to tumor evolution toward metastasis, and many cellular and molecular elements of the microenvironment are emerging as attractive targets for therapeutic strategies.

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.

Qualitative analysis of barriers and facilitators encountered by HIV patients in an ART adherence programme.

Background Medication adherence is a complex, dynamic and changing behaviour that is affected by a variety of factors, including the patient's beliefs and life circumstances. Studies have highlighted barriers to medication adherence (e.g., unmanaged side effects or a lack of social support), as well as facilitators of medication adherence (e.g., technical simplicity of treatment and psychological acceptance of the disease). Since August 2004, in Lausanne (Switzerland), physicians have referred patients who are either experiencing or are at risk of experiencing problems with their HIV antiretroviral treatment (ART) to a routine interdisciplinary ART adherence programme. This programme consists of multifactorial intervention including electronic drug monitoring (MEMS(TM)). Objective This study's objective was to identify the barriers and facilitators encountered by HIV patients with suboptimal medication adherence (≤90 % adherence over the study period). Setting The community pharmacy of the Department of Ambulatory Care and Community Medicine in Lausanne (Switzerland). Method The study consisted of a retrospective, qualitative, thematic content analysis of pharmacists' notes that were taken during semi-structured interviews with patients and conducted as part of the ART adherence programme between August 2004 and May 2008. Main outcome measure Barriers and facilitators encountered by HIV patients. Results Barriers to and facilitators of adherence were identified for the 17 included patients. These factors fell into three main categories: (1) cognitive, emotional and motivational; (2) environmental, organisational and social; and (3) treatment and disease. Conclusion The pharmacists' notes revealed that diverse barriers and facilitators were discussed during medication adherence interviews. Indeed, the results showed that the 17 non-adherent patients encountered barriers and benefited from facilitators. Therefore, pharmacists should inquire about all factors, regardless of whether they have a negative or a positive impact on medication adherence, and should consider all dimensions of patient adherence. The simultaneous strengthening of facilitators and better management of barriers may allow healthcare providers to tailor care to a patient's specific needs and support each individual patient in improving his medication-related behaviour.

Evolution and functional characterisation of the IR chemosensory receptors in the Drosophila larval gustatory system

Chemosensory receptor gene families encode divergent proteins capable of detecting a huge diversity of environmental stimuli that are constantly changing over evolutionary time as organisms adapt to distinct ecological niches. While olfaction is dedicated to the detection of volatile compounds, taste is key to assess food quality for nutritional value and presence of toxic substances. The sense of taste also provides initial signals to mediate endocrine regulation of appetite and food metabolism and plays a role in kin recognition. The fruit fly Drosophila melanogaster is a very good model for studying smell and taste because these senses are very important in insects and because a broad variety of genetic tools are available in Drosophila. Recently, a family of 66 chemosensory receptors, the Ionotropic Receptors (IRs) was described in fruit flies. IRs are distantly related to ionotropic glutamate receptors (iGluRs), but their evolutionary origin from these synaptic receptors is unclear. While 16 IRs are expressed in the olfactory system, nothing is known about the other members of this repertoire. In this thesis, I describe bioinformatic, expression and functional analyses of the IRs aimed at understanding how these receptors have evolved, and at characterising the role of the non-olfactory IRs. I show that these have emerged at the basis of the protostome lineage and probably have acquired their sensory function very early. Moreover, although several IRs are conserved across insects, there are rapid and dramatic changes in the size and divergence of IR repertoires across species. I then performed a comprehensive analysis of IR expression in the larva of Drosophila melanogaster, which is a good model to study taste and feeding mechanisms as it spends most of its time eating or foraging. I found that most of the divergent members of the IR repertoire are expressed in both peripheral and internal gustatory neurons, suggesting that these are involved in taste perception. Finally, through the establishment of a new neurophysiological assay in larvae, I identified for the first time subsets of IR neurons that preferentially detect sugars and amino acids, indicating that IRs might be involved in sensing these compounds. Together, my results indicate that IRs are an evolutionarily dynamic and functionally versatile family of receptors. In contrast to the olfactory IRs that are well-conserved, gustatory IRs are rapidly evolving species-specific receptors that are likely to be involved in detecting a wide variety of tastants. - La plupart des animaux possèdent de grandes familles de récepteurs chimiosensoriels dont la fonction est de détecter l'immense diversité de composés chimiques présents dans l'environnement. Ces récepteurs évoluent en même temps que les organismes s'adaptent à leur écosystème. Il existe deux manières de percevoir ces signaux chimiques : l'olfaction et le goût. Alors que le système olfactif perçoit les composés volatiles, le sens du goût permet d'évaluer, par contact, la qualité de la nourriture, de détecter des substances toxiques et de réguler l'appétit et le métabolisme. L'un des organismes modèles les plus pertinents pour étudier le sens du goût est le stade larvaire de la mouche du vinaigre Drosophila melanogaster. En effet, la principale fonction du stade larvaire est de trouver de la nourriture et de manger. De plus, il est possible d'utiliser tous les outils génétiques développés chez la drosophile. Récemment, une nouvelle famille de 66 récepteurs chimiosensoriels appelés Récepteurs Ionotropiques (IRs) a été découverte chez la drosophile. Bien que leur orogine soit peu claire, ces récepteurs sont similaires aux récepteurs ionotropiques glutamatergiques impliqués dans la transmission synaptique. 16 IRs sont exprimés dans le système olfactif de la mouche adulte, mais pour l'instant on ne connaît rien des autres membres de cette famille. Durant ma thèse, j'ai effectué des recherches sur l'évolution de ces récepteurs ainsi que sur l'expression et la fonction des IRs non olfactifs. Je démontre que les IRs sont apparus chez l'ancêtre commun des protostomiens et ont probablement acquis leur fonction sensorielle très rapidement. De plus, bien qu'un certain nombre d'IRs olfactifs soient conservés chez les insectes, d'importantes variations dans la taille et la divergence des répertoires d'IRs entre les espèces ont été constatées. J'ai également découvert qu'un grand nombre d'IRs non olfactifs sont exprimés dans différents organes gustatifs, ce qui leur confère probablement une fonction dans la perception des goûts. Finalement, pour la première fois, des neurones exprimant des IRs ont été identifiés pour leur fonction dans la perception de sucres et d'acides aminés chez la larve. Mes résultats présentent les IRs comme une famille très dynamique, aux fonctions très variées, qui joue un rôle tant dans l'odorat que dans le goût, et dont la fonction est restée importante tout au long de l'évolution. De plus, l'identification de neurones spécialisés dans la perception de certains composés permettra l'étude des circuits neuronaux impliqués dans le traitement de ces informations.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

Risk and global change : developing scientific methods for advocacy and awareness raising

Les catastrophes sont souvent perçues comme des événements rapides et aléatoires. Si les déclencheurs peuvent être soudains, les catastrophes, elles, sont le résultat d'une accumulation des conséquences d'actions et de décisions inappropriées ainsi que du changement global. Pour modifier cette perception du risque, des outils de sensibilisation sont nécessaires. Des méthodes quantitatives ont été développées et ont permis d'identifier la distribution et les facteurs sous- jacents du risque.¦Le risque de catastrophes résulte de l'intersection entre aléas, exposition et vulnérabilité. La fréquence et l'intensité des aléas peuvent être influencées par le changement climatique ou le déclin des écosystèmes, la croissance démographique augmente l'exposition, alors que l'évolution du niveau de développement affecte la vulnérabilité. Chacune de ses composantes pouvant changer, le risque est dynamique et doit être réévalué périodiquement par les gouvernements, les assurances ou les agences de développement. Au niveau global, ces analyses sont souvent effectuées à l'aide de base de données sur les pertes enregistrées. Nos résultats montrent que celles-ci sont susceptibles d'être biaisées notamment par l'amélioration de l'accès à l'information. Elles ne sont pas exhaustives et ne donnent pas d'information sur l'exposition, l'intensité ou la vulnérabilité. Une nouvelle approche, indépendante des pertes reportées, est donc nécessaire.¦Les recherches présentées ici ont été mandatées par les Nations Unies et par des agences oeuvrant dans le développement et l'environnement (PNUD, l'UNISDR, la GTZ, le PNUE ou l'UICN). Ces organismes avaient besoin d'une évaluation quantitative sur les facteurs sous-jacents du risque, afin de sensibiliser les décideurs et pour la priorisation des projets de réduction des risques de désastres.¦La méthode est basée sur les systèmes d'information géographique, la télédétection, les bases de données et l'analyse statistique. Une importante quantité de données (1,7 Tb) et plusieurs milliers d'heures de calculs ont été nécessaires. Un modèle de risque global a été élaboré pour révéler la distribution des aléas, de l'exposition et des risques, ainsi que pour l'identification des facteurs de risque sous- jacent de plusieurs aléas (inondations, cyclones tropicaux, séismes et glissements de terrain). Deux indexes de risque multiples ont été générés pour comparer les pays. Les résultats incluent une évaluation du rôle de l'intensité de l'aléa, de l'exposition, de la pauvreté, de la gouvernance dans la configuration et les tendances du risque. Il apparaît que les facteurs de vulnérabilité changent en fonction du type d'aléa, et contrairement à l'exposition, leur poids décroît quand l'intensité augmente.¦Au niveau local, la méthode a été testée pour mettre en évidence l'influence du changement climatique et du déclin des écosystèmes sur l'aléa. Dans le nord du Pakistan, la déforestation induit une augmentation de la susceptibilité des glissements de terrain. Les recherches menées au Pérou (à base d'imagerie satellitaire et de collecte de données au sol) révèlent un retrait glaciaire rapide et donnent une évaluation du volume de glace restante ainsi que des scénarios sur l'évolution possible.¦Ces résultats ont été présentés à des publics différents, notamment en face de 160 gouvernements. Les résultats et les données générées sont accessibles en ligne (http://preview.grid.unep.ch). La méthode est flexible et facilement transposable à des échelles et problématiques différentes, offrant de bonnes perspectives pour l'adaptation à d'autres domaines de recherche.¦La caractérisation du risque au niveau global et l'identification du rôle des écosystèmes dans le risque de catastrophe est en plein développement. Ces recherches ont révélés de nombreux défis, certains ont été résolus, d'autres sont restés des limitations. Cependant, il apparaît clairement que le niveau de développement configure line grande partie des risques de catastrophes. La dynamique du risque est gouvernée principalement par le changement global.¦Disasters are often perceived as fast and random events. If the triggers may be sudden, disasters are the result of an accumulation of actions, consequences from inappropriate decisions and from global change. To modify this perception of risk, advocacy tools are needed. Quantitative methods have been developed to identify the distribution and the underlying factors of risk.¦Disaster risk is resulting from the intersection of hazards, exposure and vulnerability. The frequency and intensity of hazards can be influenced by climate change or by the decline of ecosystems. Population growth increases the exposure, while changes in the level of development affect the vulnerability. Given that each of its components may change, the risk is dynamic and should be reviewed periodically by governments, insurance companies or development agencies. At the global level, these analyses are often performed using databases on reported losses. Our results show that these are likely to be biased in particular by improvements in access to information. International losses databases are not exhaustive and do not give information on exposure, the intensity or vulnerability. A new approach, independent of reported losses, is necessary.¦The researches presented here have been mandated by the United Nations and agencies working in the development and the environment (UNDP, UNISDR, GTZ, UNEP and IUCN). These organizations needed a quantitative assessment of the underlying factors of risk, to raise awareness amongst policymakers and to prioritize disaster risk reduction projects.¦The method is based on geographic information systems, remote sensing, databases and statistical analysis. It required a large amount of data (1.7 Tb of data on both the physical environment and socio-economic parameters) and several thousand hours of processing were necessary. A comprehensive risk model was developed to reveal the distribution of hazards, exposure and risk, and to identify underlying risk factors. These were performed for several hazards (e.g. floods, tropical cyclones, earthquakes and landslides). Two different multiple risk indexes were generated to compare countries. The results include an evaluation of the role of the intensity of the hazard, exposure, poverty, governance in the pattern and trends of risk. It appears that the vulnerability factors change depending on the type of hazard, and contrary to the exposure, their weight decreases as the intensity increases.¦Locally, the method was tested to highlight the influence of climate change and the ecosystems decline on the hazard. In northern Pakistan, deforestation exacerbates the susceptibility of landslides. Researches in Peru (based on satellite imagery and ground data collection) revealed a rapid glacier retreat and give an assessment of the remaining ice volume as well as scenarios of possible evolution.¦These results were presented to different audiences, including in front of 160 governments. The results and data generated are made available online through an open source SDI (http://preview.grid.unep.ch). The method is flexible and easily transferable to different scales and issues, with good prospects for adaptation to other research areas. The risk characterization at a global level and identifying the role of ecosystems in disaster risk is booming. These researches have revealed many challenges, some were resolved, while others remained limitations. However, it is clear that the level of development, and more over, unsustainable development, configures a large part of disaster risk and that the dynamics of risk is primarily governed by global change.

Unexplored aspects of bureaucratic autonomy: a state of the field and ways forward

This article first provides a selective overview of the literature on bureaucratic autonomy and identifies different approaches to this topic. The second section discusses three major sets of open questions, which will be tackled in the contributions to this special issue: the subjective, dynamic and relational nature of autonomy; the complex linkages between tasks, organizational forms, and national path dependencies on the one hand and autonomy and performance on the other hand; and the interplay between autonomy, accountability and democratic legitimacy.

Stochasticity and bistability in horizontal transfer control of a genomic island in Pseudomonas.

Genomic islands (GEI) comprise a recently recognized large family of potentially mobile DNA elements and play an important role in the rapid differentiation and adaptation of bacteria. Most importantly, GEIs have been implicated in the acquisition of virulence factors, antibiotic resistances or toxic compound metabolism. Despite detailed information on coding capacities of GEIs, little is known about the regulatory decisions in individual cells controlling GEI transfer. Here, we show how self-transfer of ICEclc, a GEI in Pseudomonas knackmussii B13 is controlled by a series of stochastic processes, the result of which is that only a few percent of cells in a population will excise ICEclc and launch transfer. Stochastic processes have been implicated before in producing bistable phenotypic transitions, such as sporulation and competence development, but never before in horizontal gene transfer (HGT). Bistability is instigated during stationary phase at the level of expression of an activator protein InrR that lays encoded on ICEclc, and then faithfully propagated to a bistable expression of the IntB13 integrase, the enzyme responsible for excision and integration of the ICEclc. Our results demonstrate how GEI of a very widespread family are likely to control their transfer rates. Furthermore, they help to explain why HGT is typically confined to few members within a population of cells. The finding that, despite apparent stochasticity, HGT rates can be modulated by external environmental conditions provides an explanation as to why selective conditions can promote DNA exchange.

Proximity of excitatory synapses and astroglial gap junctions in layer IV of the mouse barrel cortex.

Neurons and astrocytes, the two major cell populations in the adult brain, are characterized by their own mode of intercellular communication--the synapses and the gap junctions (GJ), respectively. In addition, there is increasing evidence for dynamic and metabolic neuroglial interactions resulting in the modulation of synaptic transmission at the so-called "tripartite synapse". Based on this, we have investigated at the ultrastructural level how excitatory synapses (ES) and astroglial GJ are spatially distributed in layer IV of the barrel cortex of the adult mouse. We used specific antibodies for connexin (Cx) 30 and 43 to identify astroglial GJ, these two proteins are known to be present in the majority of astroglial GJ in the cerebral cortex. In electron-microscopic images, we measured the distance between two ES, between two GJ and between a GJ and its nearest ES. We found a ratio of two GJ per three ES in the hollow and septal areas. Taking into account the size of an astrocyte domain, the high density of GJ suggests the occurrence of reflexive type, i.e. GJ between processes of the same astrocyte. Interestingly, the distance between an ES and an astroglial GJ was found to be significantly lower than that between either two synapses or between two GJ. These observations indicate that the two modes of cell-to-cell communication are not randomly distributed in layer IV of the barrel cortex. Consequently, this feature may provide the morphological support for the recently reported functional interactions between neuronal circuits and astroglial networks.