CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

Nava enhances ventilatory variability and diaphragmatic activity/tidal volume coupling

INTRODUCTION. Neurally Adjusted Ventilatory Assist (NAVA) is a new ventilatory mode in which ventilator settings are adjusted based on the electrical activity detected in the diaphragm (Eadi). This mode offers significant advantages in mechanical ventilation over standard pressure support (PS) modes, since ventilator input is determined directly from patient ventilatory demand. Therefore, it is expected that tidal volume (Vt) under NAVA would show better correlation with Eadi compared with PS, and exhibit greater variability due to the variability in the Eadi input to the ventilator. OBJECTIVES. To compare tidal volume variability in PS and NAVA ventilation modes, and its correlation with patient ventilatory demand (as characterized by maximum Eadi). METHODS. Acomparative study of patient-ventilator interaction was performed for 22 patients during standard PS with clinician determined ventilator settings; and NAVA, with NAVA gain set to ensure the same peak airway pressure as the total pressure obtained in PS. A 20 min continuous recording was performed in each ventilator mode. Respiratory rate, Vt, and Eadi were recorded. Tidal volume variance and Pearson correlation coefficient between Vt and Eadi were calculated for each patient. A periodogram was plotted for each ventilator mode and each patient, showing spectral power as a function of frequency to assess variability. RESULTS. Median, lower quartile and upper quartile values for Vt variance and Vt/Eadi correlation are shown in Table 1. The NAVA cohort exhibits substantially greater correlation and variance than the PS cohort. Power spectrums for Vt and Eadi are shown in Fig. 1 (PS and NAVA) for a typical patient. The enlarged section highlights how changes in Eadi are highly synchronized with NAVA ventilation, but less so for PS. CONCLUSIONS. There is greater variability in tidal volume and correlation between tidal volume and diaphragmatic electrical activity with NAVA compared to PS. These results are consistent with the improved patient-ventilator synchrony reported in the literature.

mitoKATP channel activation in the postanoxic developing heart protects E-C coupling via NO-, ROS-, and PKC-dependent pathways.

Whereas previous studies have shown that opening of the mitochondrial ATP-sensitive K(+) (mitoK(ATP)) channel protects the adult heart against ischemia-reperfusion injury, it remains to be established whether this mechanism also operates in the developing heart. Isolated spontaneously beating hearts from 4-day-old chick embryos were subjected to 30 min of anoxia followed by 60 min of reoxygenation. The chrono-, dromo-, and inotropic disturbances, as well as alterations of the electromechanical delay (EMD), reflecting excitation-contraction (E-C) coupling, were investigated. Production of reactive oxygen species (ROS) in the ventricle was determined using the intracellular fluorescent probe 2',7'-dichlorofluorescin (DCFH). Effects of the specific mitoK(ATP) channel opener diazoxide (Diazo, 50 microM) or the blocker 5-hydroxydecanoate (5-HD, 500 microM), the nitric oxide synthase (NOS) inhibitor N(G)-nitro-L-arginine methyl ester (L-NAME, 50 microM), the antioxidant N-(2-mercaptopropionyl)glycine (MPG, 1 mM), and the PKC inhibitor chelerythrine (Chel, 5 microM) on oxidative stress and postanoxic functional recovery were determined. Under normoxia, the baseline parameters were not altered by any of these pharmacological agents, alone or in combination. During the first 20 min of postanoxic reoxygenation, Diazo doubled the peak of ROS production and, interestingly, accelerated recovery of ventricular EMD and the PR interval. Diazo-induced ROS production was suppressed by 5-HD, MPG, or L-NAME, but not by Chel. Protection of ventricular EMD by Diazo was abolished by 5-HD, MPG, L-NAME, or Chel, whereas protection of the PR interval was abolished by L-NAME exclusively. Thus pharmacological opening of the mitoK(ATP) channel selectively improves postanoxic recovery of cell-to-cell communication and ventricular E-C coupling. Although the NO-, ROS-, and PKC-dependent pathways also seem to be involved in this cardioprotection, their interrelation in the developing heart can differ markedly from that in the adult myocardium.

Interhemispheric coupling between the posterior sylvian regions impacts successful auditory temporal order judgment

Introduction: Accurate registration of the relative timing between the occurrence of sensory events on a sub-second time scale is crucial for both sensory-motor and cognitive functions (Mauk and Buonomano, 2004; Habib, 2000). Support for this assumption comes notably from evidence that temporal processing impairments are implicated in a range of neurological and psychiatric conditions (e.g. Buhusi & Meck, 2005). For instance, deficits in fast auditory temporal integration have been regularly put forward as resulting in phonologic discrimination impairments at the basis of speech comprehension deficits characterizing e.g. dyslexia (Habib, 2000). At least two aspects of the brain mechanisms of temporal order judgment remain unknown. First, it is unknown when during the course of stimulus processing a temporal ,,stamp‟ is established to guide TOJ perception. Second, the extent of interplay between the cerebral hemispheres in engendering accurate TOJ performance is unresolved Methods: We investigated the spatiotemporal brain dynamics of auditory temporal order judgment (aTOJ) using electrical neuroimaging analyses of auditory evoked potentials (AEPs) recorded while participants completed a near-threshold task requiring spatial discrimination of left-right and right-left sound sequences. Results: AEPs to sound pairs modulated topographically as a function of aTOJ accuracy over the 39-77ms post-stimulus period, indicating the engagement of distinct configurations of brain networks during early auditory processing stages. Source estimations revealed that accurate and inaccurate performance were linked to bilateral posterior sylvian regions activity (PSR). However, activity within left, but not right, PSR predicted behavioral performance suggesting that left PSR activity during early encoding phases of pairs of auditory spatial stimuli appears critical for the perception of their order of occurrence. Correlation analyses of source estimations further revealed that activity between left and right PSR was significantly correlated in the inaccurate but not accurate condition, indicating that aTOJ accuracy depends on the functional de-coupling between homotopic PSR areas. Conclusions: These results support a model of temporal order processing wherein behaviorally relevant temporal information - i.e. a temporal 'stamp'- is extracted within the early stages of cortical processes within left PSR but critically modulated by inputs from right PSR. We discuss our results with regard to current models of temporal of temporal order processing, namely gating and latency mechanisms.

CD8beta endows CD8 with efficient coreceptor function by coupling T cell receptor/CD3 to raft-associated CD8/p56(lck) complexes.

The extraordinary sensitivity of CD8+ T cells to recognize antigen impinges to a large extent on the coreceptor CD8. While several studies have shown that the CD8beta chain endows CD8 with efficient coreceptor function, the molecular basis for this is enigmatic. Here we report that cell-associated CD8alphabeta, but not CD8alphaalpha or soluble CD8alphabeta, substantially increases the avidity of T cell receptor (TCR)-ligand binding. To elucidate how the cytoplasmic and transmembrane portions of CD8beta endow CD8 with efficient coreceptor function, we examined T1.4 T cell hybridomas transfected with various CD8beta constructs. T1.4 hybridomas recognize a photoreactive Plasmodium berghei circumsporozoite (PbCS) peptide derivative (PbCS (4-azidobezoic acid [ABA])) in the context of H-2K(d), and permit assessment of TCR-ligand binding by TCR photoaffinity labeling. We find that the cytoplasmic portion of CD8beta, mainly due to its palmitoylation, mediates partitioning of CD8 in lipid rafts, where it efficiently associates with p56(lck). In addition, the cytoplasmic portion of CD8beta mediates constitutive association of CD8 with TCR/CD3. The resulting TCR-CD8 adducts exhibit high affinity for major histocompatibility complex (MHC)-peptide. Importantly, because CD8alphabeta partitions in rafts, its interaction with TCR/CD3 promotes raft association of TCR/CD3. Engagement of these TCR/CD3-CD8/lck adducts by multimeric MHC-peptide induces activation of p56(lck) in rafts, which in turn phosphorylates CD3 and initiates T cell activation.

Vasopressin-dependent coupling between sodium transport and water flow in a mouse cortical collecting duct cell line.

Water balance is achieved through the ability of the kidney to control water reabsorption in the connecting tubule and the collecting duct. In a mouse cortical collecting duct cell line (mCCD(c11)), physiological concentrations of arginine vasopressin increased both electrogenic, amiloride-sensitive, epithelial sodium channel (ENaC)-mediated sodium transport measured by the short-circuit current (Isc) method and water flow (Jv apical to basal) measured by gravimetry with similar activation coefficient K(1/2) (6 and 12 pM, respectively). Jv increased linearly according to the osmotic gradient across the monolayer. A small but highly significant Jv was also measured under isoosmotic conditions. To test the coupling between sodium reabsorption and water flow, mCCD(c11) cells were treated for 24 h under isoosmotic condition with either diluent, amiloride, vasopressin or vasopressin and amiloride. Isc, Jv, and net chemical sodium fluxes were measured across the same monolayers. Around 30% of baseline and 50% of vasopressin-induced water flow is coupled to an amiloride-sensitive, ENaC-mediated, electrogenic sodium transport, whereas the remaining flow is coupled to an amiloride-insensitive, nonelectrogenic sodium transport mediated by an unknown electroneutral transporter. The mCCD(c11) cell line is a first example of a mammalian tight epithelium allowing quantitative study of the coupling between sodium and water transport. Our data are consistent with the 'near isoosmotic' fluid transport model.

Echo-time independent signal modulations for strongly coupled systems in triple echo localization schemes: an extension of S-PRESS editing.

The double spin-echo point resolved spectroscopy sequence (PRESS) is a widely used method and standard in clinical MR spectroscopy. Existence of important J-modulations at constant echo times, depending on the temporal delays between the rf-pulses, have been demonstrated recently for strongly coupled spin systems and were exploited for difference editing, removing singlets from the spectrum (strong-coupling PRESS, S-PRESS). A drawback of this method for in vivo applications is that large signal modulations needed for difference editing occur only at relatively long echo times. In this work we demonstrate that, by simply adding a third refocusing pulse (3S-PRESS), difference editing becomes possible at substantially shorter echo times while, as applied to citrate, more favorable lineshapes can be obtained. For the example of an AB system an analytical description of the MR signal, obtained with this triple refocusing sequence (3S-PRESS), is provided.

Feasibility and electromagnetic compatibility study of the ClearPEM front-end electronics for simultaneous PET-MR imaging

In this work we present a first feasibility study of the ClearPEM technology for simultaneous PET-MR imaging. The mutual electromagnetic interference (EMI) effects between both systems were evaluated on a 7 T magnet by characterizing the response behavior of the ClearPEM detectors and front-end electronics to pulsed RF power and switched magnetic field gradients; and by analyzing the MR system performance degradation from noise pickup into the RF receiver chain, and from magnetic susceptibility artifacts caused by PET front-end materials.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

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

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

Parvalbumin interneurons mediate neuronal circuitry-neurogenesis coupling in the adult hippocampus.

Using immunohistology, electron microscopy, electrophysiology and optogenetics, we found that proliferating adult mouse hippocampal neural precursors received immature GABAergic synaptic inputs from parvalbumin-expressing interneurons. Recently shown to suppress adult quiescent neural stem cell activation, parvalbumin interneuron activation promoted newborn neuronal progeny survival and development. Our results suggest a niche mechanism involving parvalbumin interneurons that couples local circuit activity to the diametric regulation of two critical early phases of adult hippocampal neurogenesis.