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.

The feasibility of 350 μm spatial resolution coronary magnetic resonance angiography at 3 T in humans.

PURPOSE: The purposes of this study were to (1) develop a high-resolution 3-T magnetic resonance angiography (MRA) technique with an in-plane resolution approximate to that of multidetector coronary computed tomography (MDCT) and a voxel size of 0.35 × 0.35 × 1.5 mm³ and to (2) investigate the image quality of this technique in healthy participants and preliminarily in patients with known coronary artery disease (CAD). MATERIALS AND METHODS: A 3-T coronary MRA technique optimized for an image acquisition voxel as small as 0.35 × 0.35 × 1.5 mm³ (high-resolution coronary MRA [HRC]) was implemented and the coronary arteries of 22 participants were imaged. These included 11 healthy participants (average age, 28.5 years; 5 men) and 11 participants with CAD (average age, 52.9 years; 5 women) as identified on MDCT. In addition, the 11 healthy participants were imaged using a method with a more common spatial resolution of 0.7 × 1 × 3 mm³ (regular-resolution coronary MRA [RRC]). Qualitative and quantitative comparisons were made between the 2 MRA techniques. RESULTS: Normal vessels and CAD lesions were successfully depicted at 350 × 350 μm² in-plane resolution with adequate signal-to-noise ratio (SNR) and contrast-to-noise ratio. The CAD findings were consistent among MDCT and HRC. The HRC showed a 47% improvement in sharpness despite a reduction in SNR (by 72%) and in contrast-to-noise ratio (by 86%) compared with the regular-resolution coronary MRA. CONCLUSION: This study, as a first step toward substantial improvement in the resolution of coronary MRA, demonstrates the feasibility of obtaining at 3 T a spatial resolution that approximates that of MDCT. The acquisition in-plane pixel dimensions are as small as 350 × 350 μm² with a 1.5-mm slice thickness. Although SNR is lower, the images have improved sharpness, resulting in image quality that allows qualitative identification of disease sites on MRA consistent with MDCT.

Acid-sensing ion channels : modulation by serine-proteases and involvement in acid-induced action potential generation in hippocampal neurons

SUMMARY Acid-sensing ion channels (ASICs) are non-voltage gated sodium channels. They are activated by rapid extracellular acidification and generate an inactivating inward current. Four ASIC genes have been cloned: ASIC1, 2, 3 and 4, with variants a and b for ASIC1and AS1C2. ASICs are expressed in neurons of the central (CNS) and peripheral nervous system (PNS). In the CNS, ASICs have a role in learning, memory, as well as in neuronal death in ischemia. In the PNS, ASICs are involved in the perception of acid-induced pain, as well as in mechanoperception. In one part of my thesis project, we addressed the question of the mechanism of regulation of ASIC1 a by the serine protease trypsin at the molecular level. Trypsin modifies the function of ASIC1 a but not of ASIC1b. In order to identify the channel region responsible for this effect, we created chimeras between ASIC1 a and 1b. Subsequently, to identify the exact trypsin target(s), we mutated predicted trypsin sites in the region identified by the chimera. In the second part of a project, we investigated the role of ASICs at the cellular level, in neuronal signaling. Using the whole-cell patch clamp in hippocampal neuronal culture, we studied the potential involvement of ASICs in action potential (AP) generation. In the first part of the thesis work, we showed that trypsin modifies ASIC1a function: it shifts the pH activation and the steady-state inactivation curve towards more acidic values and accelerates the time course of the channel recovery from inactivation. We also showed that trypsin cleaves ASIC1a and that the functional effect and a channel cleavage correlate. In the inactivated state, channels cannot be modified by trypsin. Cleavage occurs in a channel region that is also important for inactivation of all ASICs; a part of this region is critical for the inhibition of ASIC1 a by the spider toxin Psalmotoxin1. In the second part of the thesis work, we showed that ASIC activity can modulate AP generation. ASIC activity by itself can induce trains of APs. In situations in which this activity by itself is not sufficient to induce APs, it can contribute to AP generation. During high neuronal activity, ASIC activity can block already existing trains of APs. In conclusion, depending on the activity of neuron in a particular moment, ASICs can differently modulate AP generation; they can induce, facilitate or inhibit APs. We also showed that trypsin changes the capability of ASICs to modulate AP generation by shifting the pH dependence to more acidic values, which adapts channel gating to pH conditions which may occur in pathological conditions such as ischemia. Our finding that trypsin modifies ASIC1 a function identifies a novel pharmacological tool, and proposes a mechanism of ASIC1a regulation that may have a physiological importance. The identification of the exact site of trypsin action gives insight to the molecular mechanisms of ASIC regulation. This work proposes a role in modulation of AP generation for ASICs in the CNS. RESUME Les canaux ASIC sont les canaux ioniques activés par l'acidification rapide extracellulaire. Activés, ils génèrent un courant entrant qui inactive en présence de stimulus acide. Quatre gènes ASIC ont été clonés, ASIC1, 2, 3 et 4, avec les variants a et b pour ASIC1 et 2. Les ASICs sont exprimés dans les neurones du système nerveux central (SNC) et périphérique (SNP). Dans le SNC, les ASIC ont un rôle dans le mémoire, apprentissage et la mort neuronale dans t'ischémie. Dans le SNP, ils ont un rôle dans la perception de la douleur et méchanosensation. Dans une partie de mon projet de thèse, nous avons étudié les mécanismes de la régulation d'ASIC1a par la sérine-protéase trypsine au niveau moléculaire. La trypsine modifie la fonction d'ASIC1a et pas ASIC1b. Nous avons créé les chimères entre ASIC1 a et 1 b, afin d'identifier la région du canal responsable pour l'effet. Pour identifier le(s) site(s) exactes de l'action de la trypsine, nous avons muté les sites potentiels de la trypsine dans la région identifiée par les chimères. Dans la deuxième partie du projet, nous avons étudié le rôle des ASICs au niveau cellulaire. En utilisant la technique du patch clamp dans les cultures des neurones de l'hippocampe, nous avons étudié l'implication des ASICs dans la génération des potentiels d'action (PA). Nous avons montré que la trypsine agit sur le canal ASIC1a ; elle décale l'activation et « steady-state » inactivation vers les valeurs plus acides, et elle raccourcit le temps du « recovery » du canal. La trypsine coupe ASIC1a sur le résidu K145 et l'effet fonctionnel et la coupure corrèlent. Nous avons identifié la région du canal responsable pour l'inactivation de tous les ASICs ; une partie de cette région est responsable pour ['inhibition d'ASIC1 a par la Psalmotoxinel . Nous avons montré que les ASICs peuvent moduler la génération des PAs. L'activité des ASICs peut induire les trains des PAs. Quand l'activité des ASICs n'est pas suffisante pour induire le PA, elle peut contribuer à sa génération. Pendant l'activité neuronale forte, l'activité des ASICs peut bloquer les trains des PAs qui existent déjà. En conclusion, dépendant de l'activité neuronale, les ASICs peuvent moduler la génération des PAs différemment ; ils peuvent induire, faciliter ou inhiber les PAs. La trypsine change la capacité des ASICs de moduler les PAs. Après l'action de la trypsine, les ASICs peuvent moduler la génération des PAs dans les conditions légèrement acides, suivies par les fluctuations du pH acide, qui peuvent exister dans l'ischémie. Le fait que la trypsine agit sur ASIC1a définit l'outil pharmacologique et propose le mécanisme de la régulation d'ASICI a qui pourrait avoir l'importance physiologique. L'identification du site de l'action de la trypsine éclaircit les mécanismes moléculaires de la régulation des ASICs. Cette étude propose un rôle des ASICs dans la modulation de la génération des PAs. Résumé pour le public large Les neurones sont les cellules de système nerveux dont la fonction est la signalisation. Comme toutes les autres cellules, les neurones ont une membrane qui sépare l'intérieur du milieu extérieur. Cette membrane est imperméable pour des particules chargées (ions). Dans cette membrane existent les protéines spécifiques, « canaux », qui permettent le transport des ions d'un côté de la membrane à l'autre, comme réponse aux stimuli différents. Ce transport des ions à travers la membrane génère un courant, qu'on peut mesurer. Ce courant est la base de la communication entre les neurones, ou, ce qu'on appelle la signalisation neuronale. Quand ce courant est suffisamment grand, il permet la génération du potentiel d'action, qui est le message principal de communication neuronale. Les canaux ASIC (acid-sensing ion channel), que nous étudions dans le laboratoire, sont activés par les acides. Les acides sont relâchés dans beaucoup de situations dans le système nerveux. Les ASIC ont été découverts récemment (en 1996), et nous ne connaissons pas encore très bien toutes les fonctions de ces canaux. Nous savons qu'ils ont un rôle dans le mémoire, apprentissage, la sensation de la douleur et l'infarctus cérébral. Dans la première partie de ce projet de thèse, nous avons voulu mieux comprendre comment fonctionnent ces canaux. Pour faire ça, nous avons étudié la régulation des ASICs par une protéine, trypsine, qui coupe le canal ASIC. Nous avons étudié ou exactement la trypsine coupe le canal et quels effets ça produit sur la fonction du canal. Dans la deuxième partie du projet de thèse, nous avons voulu mieux connaître comment le canal fonctionne au niveau de la cellule, comment il interagit avec les autres canaux et si il a un rôle dans la génération des potentiels d'action. Nous avons pu montrer que la trypsine change la fonction du canal, ce qui lui permet de fonctionner différemment. Nous avons aussi déterminé ou exactement ta trypsine coupe le canal. Au niveau de la cellule, nous avons montré que les ASIC peuvent moduler la génération des potentiels d'action, étant, dépendant de l'activité du neurone, soit activateurs, soit inhibiteurs. La trypsine est une molécule qui peut être libérée dans le système nerveux pendant certaines conditions, comme l'infarctus cérébral. A cause de ça, les connaissances que la trypsine agit sur le anal ASIC pourraient être important physiologiquement. La connaissance de l'endroit exacte ou la trypsine coupe le canal nous aide à mieux comprendre la relation structure-fonction du canal. La modulation de la génération des potentiels d'actions par les ASIC indique que ces canaux peuvent avoir un rôle important dans la signalisation neuronale.

On the origin of the MR image phase contrast: an in vivo MR microscopy study of the rat brain at 14.1 T.

Recent studies at high magnetic fields using the phase of gradient-echo MR images have shown the ability to unveil cortical substructure in the human brain. To investigate the contrast mechanisms in phase imaging, this study extends, for the first time, phase imaging to the rodent brain. Using a 14.1 T horizontal bore animal MRI scanner for in vivo micro-imaging, images with an in-plane resolution of 33 microm were acquired. Phase images revealed, often more clearly than the corresponding magnitude images, hippocampal fields, cortical layers (e.g. layer 4), cerebellar layers (molecular and granule cell layers) and small white matter structures present in the striatum and septal nucleus. The contrast of the phase images depended in part on the orientation of anatomical structures relative to the magnetic field, consistent with bulk susceptibility variations between tissues. This was found not only for vessels, but also for white matter structures, such as the anterior commissure, and cortical layers in the cerebellum. Such susceptibility changes could result from variable blood volume. However, when the deoxyhemoglobin content was reduced by increasing cerebral blood flow (CBF) with a carbogen breathing challenge, contrast between white and gray matter and cortical layers was not affected, suggesting that tissue cerebral blood volume (and therefore deoxyhemoglobin) is not a major source of the tissue phase contrast. We conclude that phase variations in gradient-echo images are likely due to susceptibility shifts of non-vascular origin.

Three-dimensional magnetic resonance myocardial motion tracking from a single image plane.

Three-dimensional imaging for the quantification of myocardial motion is a key step in the evaluation of cardiac disease. A tagged magnetic resonance imaging method that automatically tracks myocardial displacement in three dimensions is presented. Unlike other techniques, this method tracks both in-plane and through-plane motion from a single image plane without affecting the duration of image acquisition. A small z-encoding gradient is subsequently added to the refocusing lobe of the slice-selection gradient pulse in a slice following CSPAMM acquisition. An opposite polarity z-encoding gradient is added to the orthogonal tag direction. The additional z-gradients encode the instantaneous through plane position of the slice. The vertical and horizontal tags are used to resolve in-plane motion, while the added z-gradients is used to resolve through-plane motion. Postprocessing automatically decodes the acquired data and tracks the three-dimensional displacement of every material point within the image plane for each cine frame. Experiments include both a phantom and in vivo human validation. These studies demonstrate that the simultaneous extraction of both in-plane and through-plane displacements and pathlines from tagged images is achievable. This capability should open up new avenues for the automatic quantification of cardiac motion and strain for scientific and clinical purposes.

Real-time motion correction in navigator-gated free-breathing double-oblique submillimeter 3D right coronary artery magnetic resonance angiography.

RATIONALE AND OBJECTIVES: The purpose of this study was the investigation of the impact of real-time adaptive motion correction on image quality in navigator-gated, free-breathing, double-oblique three-dimensional (3D) submillimeter right coronary magnetic resonance angiography (MRA). MATERIALS AND METHODS: Free-breathing 3D right coronary MRA with real-time navigator technology was performed in 10 healthy adult subjects with an in-plane spatial resolution of 700 x 700 microm. Identical double-oblique coronary MR-angiograms were performed with navigator gating alone and combined navigator gating and real-time adaptive motion correction. Quantitative objective parameters of contrast-to-noise ratio (CNR) and vessel sharpness and subjective image quality scores were compared. RESULTS: Superior vessel sharpness, increased CNR, and superior image quality scores were found with combined navigator gating and real-time adaptive motion correction (vs. navigator gating alone; P < 0.01 for all comparisons). CONCLUSION: Real-time adaptive motion correction objectively and subjectively improves image quality in 3D navigator-gated free-breathing double-oblique submillimeter right coronary MRA.

ZHARP: three-dimensional motion tracking from a single image plane.

Three-dimensional imaging and quantification of myocardial function are essential steps in the evaluation of cardiac disease. We propose a tagged magnetic resonance imaging methodology called zHARP that encodes and automatically tracks myocardial displacement in three dimensions. Unlike other motion encoding techniques, zHARP encodes both in-plane and through-plane motion in a single image plane without affecting the acquisition speed. Postprocessing unravels this encoding in order to directly track the 3-D displacement of every point within the image plane throughout an entire image sequence. Experimental results include a phantom validation experiment, which compares zHARP to phase contrast imaging, and an in vivo study of a normal human volunteer. Results demonstrate that the simultaneous extraction of in-plane and through-plane displacements from tagged images is feasible.

Brain Connectivity Analysis in Children. Effects of Prematurity. and Prenatal Growth Restriction

Introduction: Survival of children born prematurely or with very low birth weight has increased dramatically, but the long term developmental outcome remains unknown. Many children have deficits in cognitive capacities, in particular involving executive domains and those disabilities are likely to involve a central nervous system deficit. To understand their neurostructural origin, we use DTI. Structurally segregated and functionally regions of the cerebral cortex are interconnected by a dense network of axonal pathways. We noninvasively map these pathways across cortical hemispheres and construct normalized structural connection matrices derived from DTI MR tractography. Group comparisons of brain connectivity reveal significant changes in fiber density in case of children with poor intrauterine grown and extremely premature children (gestational age<28 weeks at birth) compared to control subjects. This changes suggest a link between cortico-axonal pathways and the central nervous system deficit. Methods: Sixty premature born infants (5-6 years old) were scanned on clinical 3T scanner (Magnetom Trio, Siemens Medical Solutions, Erlangen, Germany) at two hospitals (HUG, Geneva and CHUV, Lausanne). For each subject, T1-weighted MPRAGE images (TR/TE=2500/2.91,TI=1100, resolution=1x1x1mm, matrix=256x154) and DTI images (30 directions, TR/TE=10200/107, in-plane resolution=1.8x1.8x2mm, 64 axial, matrix=112x112) were acquired. Parent(s) provided written consent on prior ethical board approval. The extraction of the Whole Brain Structural Connectivity Matrix was performed following (Cammoun, 2009 and Hagmann, 2008). The MPARGE images were registered using an affine registration to the non-weighted-DTI and WM-GM segmentation performed on it. In order to have equal anatomical localization among subjects, 66 cortical regions with anatomical landmarks were created using the curvature information, i.e. sulcus and gyrus (Cammoun et al, 2007; Fischl et al, 2004; Desikan et al, 2006) with freesurfer software (http://surfer.nmr.mgh.harvard.edu/). Tractography was performed in WM using an algorithm especially designed for DTI/DSI data (Hagmann et al., 2007) and both information were then combined in a matrix. Each row and column of the matrix corresponds to a particular ROI. Each cell of index (i,j) represents the fiber density of the bundle connecting the ROIs i and j. Subdividing each cortical region, we obtained 4 Connectivity Matrices of different resolution (33, 66, 125 and 250 ROI/hemisphere) for each subject . Subjects were sorted in 3 different groups, namely (1) control, (2) Intrauterine Growth Restriction (IUGR), (3) Extreme Prematurity (EP), depending on their gestational age, weight and percentile-weight score at birth. Group-to-group comparisons were performed between groups (1)-(2) and (1)-(3). The mean age at examination of the three groups were similar. Results: Quantitative analysis were performed between groups to determine fibers density differences. For each group, a mean connectivity matrix with 33ROI/hemisphere resolution was computed. On the other hand, for all matrix resolutions (33,66,125,250 ROI/hemisphere), the number of bundles were computed and averaged. As seen in figure 1, EP and IUGR subjects present an overall reduction of fibers density in both interhemispherical and intrahemispherical connections. This is given quantitatively in table 1. IUGR subjects presents a higher percentage of missing fiber bundles than EP when compared to control subjects (~16% against 11%). When comparing both groups to control subjects, for the EP subjects, the occipito-parietal regions seem less interhemispherically connected whilst the intrahemispherical networks present lack of fiber density in the lymbic system. Children born with IUGR, have similar reductions in interhemispherical connections than the EP. However, the cuneus and precuneus connections with the precentral and paracentral lobe are even lower than in the case of the EP. For the intrahemispherical connections the IUGR group preset a loss of fiber density between the deep gray matter structures (striatum) and the frontal and middlefrontal poles, connections typically involved in the control of executive functions. For the qualitative analysis, a t-test comparing number of bundles (p-value<0.05) gave some preliminary significant results (figure 2). Again, even if both IUGR and EP appear to have significantly less connections comparing to the control subjects, the IUGR cohort seems to present a higher lack of fiber density specially relying the cuneus, precuneus and parietal areas. In terms of fiber density, preliminary Wilcoxon tests seem to validate the hypothesis set by the previous analysis. Conclusions: The goal of this study was to determine the effect of extreme prematurity and poor intrauterine growth on neurostructural development at the age of 6 years-old. This data indicates that differences in connectivity may well be the basis for the neurostructural and neuropsychological deficit described in these populations in the absence of overt brain lesions (Inder TE, 2005; Borradori-Tolsa, 2004; Dubois, 2008). Indeed, we suggest that IUGR and prematurity leads to alteration of connectivity between brain structures, especially in occipito-parietal and frontal lobes for EP and frontal and middletemporal poles for IUGR. Overall, IUGR children have a higher loss of connectivity in the overall connectivity matrix than EP children. In both cases, the localized alteration of connectivity suggests a direct link between cortico-axonal pathways and the central nervous system deficit. Our next step is to link these connectivity alterations to the performance in executive function tests.

Cerebellar cortical layers: in vivo visualization with structural high-field-strength MR imaging.

Purpose: To perform in vivo imaging of the cerebellum with an in-plane resolution of 120 mm to observe its cortical granular and molecular layers by taking advantage of the high signal-to-noise ratio and the increased magnetic susceptibility-related contrast available at high magnetic field strength such as 7 T. Materials and Methods: The study was approved by the institutional review board, and all patients provided written consent. Three healthy persons (two men, one woman; mean age, 30 years; age range, 28-31 years) underwent MR imaging with a 7-T system. Gradient-echo images (repetition time msec/echo time msec, 1000/25) of the human cerebellum were acquired with a nominal in-plane resolution of approximately 120 mum and a section thickness of 1 mm. Results: Structures with dimensions as small as 240 mum, such as the granular and molecular layers in the cerebellar cortex, were detected in vivo. The detection of these structures was confirmed by comparing the contrast obtained on T2*-weighted and phase images with that obtained on images of rat cerebellum acquired at 14 T with 30 mum in-plane resolution. Conclusion: In vivo cerebellar imaging at near-microscopic resolution is feasible at 7 T. Such detailed observation of an anatomic area that can be affected by a number of neurologic and psychiatric diseases, such as stroke, tumors, autism, and schizophrenia, could potentially provide newer markers for diagnosis and follow-up in patients with such pathologic conditions. (c) RSNA, 2010.

Pizgrischite, (Cu,Fe)Cu14PbBi17S35, a new sulfosalt from the Swiss Alps: Description, crystal structure and occurrence

Pizgrischite, (Cu,Fe)Cu14PbBi17S35, is a new mineral species named after the type locality, Piz Grisch Mountain, Val Ferrera, Graubunden, Switzerland. This sulfosalt occurs as thin, striated, metallic lead-grey blades measuring up to I cm in length, embedded in quartz and associated with tetrahedrite, chalcopyrite, pyrite, sphalerite, emplectite and derivatives of the aikinite-bismuthinite series. In plane-polarized light, the new species is brownish grey with no perceptible pleochroism; under crossed nicols in oil immersion, it presents a weak anisotropy with dark brown tints. Minimum and maximum reflectance values (in %) in air are: 40.7-42.15 (470 nm), 41.2-43.1 (546 nm), 41.2-43.35 (589 nm) and 40.7-43.3 (650 nm). Cleavage is perfect along 001 I and well developed on {010}. Abundant polysynthetic twinning is observed on (010). The mean micro-indentation hardness is 190 kg/mm(2) (Mohs hardness 3.3), and the calculated density is 6.58 g/cm(3). Electron-microprobe analyses yield (wt%; mean result of seven analyses): Cu 16.48, Pb 2.10, Fe 0.77, Bi 60.70, Sb 0.35, S 19.16, Se 0.04, total 99.60. The resulting empirical chemical formula is (Cu15.24Fe0.80Pb0.60)(Sigma 16.64)(Bi17.07Sb0.17)(Sigma 17.24)(S35.09Se0.03)(Sigma 35.12), in accordance with the formula derived from the single-crystal refinement of the structure, (Cu,Fe)Cu14PbBi17S35. Pizgrischite is monoclinic, space group C2/m, with the following unit-cell parameters: a 35.054(2), b3.91123(I), c43.192(2) angstrom, beta 96.713(4)degrees, V5881.24 angstrom(3), Z=4. The strongest seven X-ray powder-diffraction lines [d in angstrom (I)(hkl)] are: 5.364(40)((6) over bar 04), 4.080(50)((8) over bar 05), 3.120(40)(118), 3.104(68)((3) over bar 18), 2.759(53) ((9) over bar 11),2.752(44)(910) and 1.956(100)(020). The crystal structure is an expanded monoclinic derivative of kupcikite. Pizgrischite belongs to the cuprobismutite series of bismuth sulfosalts but, sensu stricto, it is not a homologue of cuprobismutite. At the type locality. pizarischite is the result of the Alpine metamorphism under greenschist-facies conditions of pre-Tertiary hydrothermal Cu-Bi mineralization.

Central and Cortical Gray Mater Segmentation of Magnetic Resonance Images of the Fetal Brain

Motivation. The study of human brain development in itsearly stage is today possible thanks to in vivo fetalmagnetic resonance imaging (MRI) techniques. Aquantitative analysis of fetal cortical surfacerepresents a new approach which can be used as a markerof the cerebral maturation (as gyration) and also forstudying central nervous system pathologies [1]. However,this quantitative approach is a major challenge forseveral reasons. First, movement of the fetus inside theamniotic cavity requires very fast MRI sequences tominimize motion artifacts, resulting in a poor spatialresolution and/or lower SNR. Second, due to the ongoingmyelination and cortical maturation, the appearance ofthe developing brain differs very much from thehomogenous tissue types found in adults. Third, due tolow resolution, fetal MR images considerably suffer ofpartial volume (PV) effect, sometimes in large areas.Today extensive efforts are made to deal with thereconstruction of high resolution 3D fetal volumes[2,3,4] to cope with intra-volume motion and low SNR.However, few studies exist related to the automatedsegmentation of MR fetal imaging. [5] and [6] work on thesegmentation of specific areas of the fetal brain such asposterior fossa, brainstem or germinal matrix. Firstattempt for automated brain tissue segmentation has beenpresented in [7] and in our previous work [8]. Bothmethods apply the Expectation-Maximization Markov RandomField (EM-MRF) framework but contrary to [7] we do notneed from any anatomical atlas prior. Data set &Methods. Prenatal MR imaging was performed with a 1-Tsystem (GE Medical Systems, Milwaukee) using single shotfast spin echo (ssFSE) sequences (TR 7000 ms, TE 180 ms,FOV 40 x 40 cm, slice thickness 5.4mm, in plane spatialresolution 1.09mm). Each fetus has 6 axial volumes(around 15 slices per volume), each of them acquired inabout 1 min. Each volume is shifted by 1 mm with respectto the previous one. Gestational age (GA) ranges from 29to 32 weeks. Mother is under sedation. Each volume ismanually segmented to extract fetal brain fromsurrounding maternal tissues. Then, in-homogeneityintensity correction is performed using [9] and linearintensity normalization is performed to have intensityvalues that range from 0 to 255. Note that due tointra-tissue variability of developing brain someintensity variability still remains. For each fetus, ahigh spatial resolution image of isotropic voxel size of1.09 mm is created applying [2] and using B-splines forthe scattered data interpolation [10] (see Fig. 1). Then,basal ganglia (BS) segmentation is performed on thissuper reconstructed volume. Active contour framework witha Level Set (LS) implementation is used. Our LS follows aslightly different formulation from well-known Chan-Vese[11] formulation. In our case, the LS evolves forcing themean of the inside of the curve to be the mean intensityof basal ganglia. Moreover, we add local spatial priorthrough a probabilistic map created by fitting anellipsoid onto the basal ganglia region. Some userinteraction is needed to set the mean intensity of BG(green dots in Fig. 2) and the initial fitting points forthe probabilistic prior map (blue points in Fig. 2). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed as described in [8]. Results.The case study presented here has 29 weeks of GA. Thehigh resolution reconstructed volume is presented in Fig.1. The steps of BG segmentation are shown in Fig. 2.Overlap in comparison with manual segmentation isquantified by the Dice similarity index (DSI) equal to0.829 (values above 0.7 are considered a very goodagreement). Such BG segmentation has been applied on 3other subjects ranging for 29 to 32 GA and the DSI hasbeen of 0.856, 0.794 and 0.785. Our segmentation of theinner (red and blue contours) and outer cortical surface(green contour) is presented in Fig. 3. Finally, torefine the results we include our WM segmentation in theFreesurfer software [12] and some manual corrections toobtain Fig.4. Discussion. Precise cortical surfaceextraction of fetal brain is needed for quantitativestudies of early human brain development. Our workcombines the well known statistical classificationframework with the active contour segmentation forcentral gray mater extraction. A main advantage of thepresented procedure for fetal brain surface extraction isthat we do not include any spatial prior coming fromanatomical atlases. The results presented here arepreliminary but promising. Our efforts are now in testingsuch approach on a wider range of gestational ages thatwe will include in the final version of this work andstudying as well its generalization to different scannersand different type of MRI sequences. References. [1]Guibaud, Prenatal Diagnosis 29(4) (2009). [2] Rousseau,Acad. Rad. 13(9), 2006, [3] Jiang, IEEE TMI 2007. [4]Warfield IADB, MICCAI 2009. [5] Claude, IEEE Trans. Bio.Eng. 51(4) (2004). [6] Habas, MICCAI (Pt. 1) 2008. [7]Bertelsen, ISMRM 2009 [8] Bach Cuadra, IADB, MICCAI 2009.[9] Styner, IEEE TMI 19(39 (2000). [10] Lee, IEEE Trans.Visual. And Comp. Graph. 3(3), 1997, [11] Chan, IEEETrans. Img. Proc, 10(2), 2001 [12] Freesurfer,http://surfer.nmr.mgh.harvard.edu.

Three-dimensional high-resolution fast spin-echo coronary magnetic resonance angiography.

Due to SNR constraints, current "bright-blood" 3D coronary MRA approaches still suffer from limited spatial resolution when compared to conventional x-ray coronary angiography. Recent 2D fast spin-echo black-blood techniques maximize signal for coronary MRA at no loss in image spatial resolution. This suggests that the extension of black-blood coronary MRA with a 3D imaging technique would allow for a further signal increase, which may be traded for an improved spatial resolution. Therefore, a dual-inversion 3D fast spin-echo imaging sequence and real-time navigator technology were combined for high-resolution free-breathing black-blood coronary MRA. In-plane image resolution below 400 microm was obtained. Magn Reson Med 45:206-211, 2001.

Magnetization transfer-based 3D visualization of foot peripheral nerves.

PURPOSE: To investigate magnetization transfer (MT) effects as a new source of contrast for imaging and tracking of peripheral foot nerves. MATERIALS AND METHODS: Two sets of 3D spoiled gradient-echo images acquired with and without a saturation pulse were used to generate MT ratio (MTR) maps of 260 μm in-plane resolution for eight volunteers at 3T. Scan parameters were adjusted to minimize signal loss due to T2 dephasing, and a dedicated coil was used to improve the inherently low signal-to-noise ratio of small voxels. Resulting MTR values in foot nerves were compared with those in surrounding muscle tissue. RESULTS: Average MTR values for muscle (45.5 ± 1.4%) and nerve (21.4 ± 3.1%) were significantly different (P < 0.0001). In general, the difference in MTR values was sufficiently large to allow for intensity-based segmentation and tracking of foot nerves in individual subjects. This procedure was termed MT-based 3D visualization. CONCLUSION: The MTR serves as a new source of contrast for imaging of peripheral foot nerves and provides a means for high spatial resolution tracking of these structures. The proposed methodology is directly applicable on standard clinical MR scanners and could be applied to systemic pathologies, such as diabetes.

Contrast agent-enhanced, free-breathing, three-dimensional coronary magnetic resonance angiography.

For free-breathing, high-resolution, three-dimensional coronary magnetic resonance angiography (MRA), the use of intravascular contrast agents may be helpful for contrast enhancement between coronary blood and myocardium. In six patients, 0.1 mmol/kg of the intravascular contrast agent MS-325/AngioMARK was given intravenously followed by double-oblique, free-breathing, three-dimensional inversion-recovery coronary MRA with real-time navigator gating and motion correction. Contrast-enhanced, three-dimensional coronary MRA images were compared with images obtained with a T2 prepulse (T2Prep) without exogenous contrast. The contrast-enhanced images demonstrated a 69% improvement in the contrast-to-noise ratio (6.6 +/- 1.1 vs. 11.1 +/- 2.5; P < 0.01) compared with the T2Prep approach. By using the intravascular agent, extensive portions (> 80 mm) of the native left and right coronary system could be displayed consistently with sub-millimeter in-plane resolution. The intravascular contrast agent, MS-325/AngioMARK, leads to a considerable enhancement of the blood/muscle contrast for coronary MRA compared with T2Prep techniques. The clinical value of the agent remains to be defined in a larger patient series. J. Magn. Reson. Imaging 1999;10:790-799.

Breathhold three-dimensional coronary magnetic resonance angiography using real-time navigator technology.

The acquisition duration of most three-dimensional (3D) coronary magnetic resonance angiography (MRA) techniques is considerably prolonged, thereby precluding breathholding as a mechanism to suppress respiratory motion artifacts. Splitting the acquired 3D volume into multiple subvolumes or slabs serves to shorten individual breathhold duration. Still, problems associated with misregistration due to inconsistent depths of expiration and diaphragmatic drift during sustained respiration remain to be resolved. We propose the combination of an ultrafast 3D coronary MRA imaging sequence with prospective real-time navigator technology, which allows correction of the measured volume position. 3D volume splitting using prospective real-time navigator technology, was successfully applied for 3D coronary MRA in five healthy individuals. An ultrafast 3D interleaved hybrid gradient-echoplanar imaging sequence, including T2Prep for contrast enhancement, was used with the navigator localized at the basal anterior wall of the left ventricle. A 9-cm-thick volume, with in-plane spatial resolution of 1.1 x 2.2 mm, was acquired during five breathholds of 15-sec duration each. Consistently, no evidence of misregistration was observed in the images. Extensive contiguous segments of the left anterior descending coronary artery (48 +/- 18 mm) and the right coronary artery (75 +/- 5 mm) could be visualized. This technique has the potential for screening for anomalous coronary arteries, making it well suited as part of a larger clinical MR examination. In addition, this technique may also be applied as a scout scan, which allows an accurate definition of imaging planes for subsequent high-resolution coronary MRA.