Investigation of discrete imaging models and iterative image reconstruction in differential X-ray phase-contrast tomography.

Differential X-ray phase-contrast tomography (DPCT) refers to a class of promising methods for reconstructing the X-ray refractive index distribution of materials that present weak X-ray absorption contrast. The tomographic projection data in DPCT, from which an estimate of the refractive index distribution is reconstructed, correspond to one-dimensional (1D) derivatives of the two-dimensional (2D) Radon transform of the refractive index distribution. There is an important need for the development of iterative image reconstruction methods for DPCT that can yield useful images from few-view projection data, thereby mitigating the long data-acquisition times and large radiation doses associated with use of analytic reconstruction methods. In this work, we analyze the numerical and statistical properties of two classes of discrete imaging models that form the basis for iterative image reconstruction in DPCT. We also investigate the use of one of the models with a modern image reconstruction algorithm for performing few-view image reconstruction of a tissue specimen.

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.

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.

Analysis and spatial modeling of the indoor radon Swiss data

The present research deals with an important public health threat, which is the pollution created by radon gas accumulation inside dwellings. The spatial modeling of indoor radon in Switzerland is particularly complex and challenging because of many influencing factors that should be taken into account. Indoor radon data analysis must be addressed from both a statistical and a spatial point of view. As a multivariate process, it was important at first to define the influence of each factor. In particular, it was important to define the influence of geology as being closely associated to indoor radon. This association was indeed observed for the Swiss data but not probed to be the sole determinant for the spatial modeling. The statistical analysis of data, both at univariate and multivariate level, was followed by an exploratory spatial analysis. Many tools proposed in the literature were tested and adapted, including fractality, declustering and moving windows methods. The use of Quan-tité Morisita Index (QMI) as a procedure to evaluate data clustering in function of the radon level was proposed. The existing methods of declustering were revised and applied in an attempt to approach the global histogram parameters. The exploratory phase comes along with the definition of multiple scales of interest for indoor radon mapping in Switzerland. The analysis was done with a top-to-down resolution approach, from regional to local lev¬els in order to find the appropriate scales for modeling. In this sense, data partition was optimized in order to cope with stationary conditions of geostatistical models. Common methods of spatial modeling such as Κ Nearest Neighbors (KNN), variography and General Regression Neural Networks (GRNN) were proposed as exploratory tools. In the following section, different spatial interpolation methods were applied for a par-ticular dataset. A bottom to top method complexity approach was adopted and the results were analyzed together in order to find common definitions of continuity and neighborhood parameters. Additionally, a data filter based on cross-validation was tested with the purpose of reducing noise at local scale (the CVMF). At the end of the chapter, a series of test for data consistency and methods robustness were performed. This lead to conclude about the importance of data splitting and the limitation of generalization methods for reproducing statistical distributions. The last section was dedicated to modeling methods with probabilistic interpretations. Data transformation and simulations thus allowed the use of multigaussian models and helped take the indoor radon pollution data uncertainty into consideration. The catego-rization transform was presented as a solution for extreme values modeling through clas-sification. Simulation scenarios were proposed, including an alternative proposal for the reproduction of the global histogram based on the sampling domain. The sequential Gaussian simulation (SGS) was presented as the method giving the most complete information, while classification performed in a more robust way. An error measure was defined in relation to the decision function for data classification hardening. Within the classification methods, probabilistic neural networks (PNN) show to be better adapted for modeling of high threshold categorization and for automation. Support vector machines (SVM) on the contrary performed well under balanced category conditions. In general, it was concluded that a particular prediction or estimation method is not better under all conditions of scale and neighborhood definitions. Simulations should be the basis, while other methods can provide complementary information to accomplish an efficient indoor radon decision making.

Predictive analysis and mapping of indoor radon concentrations in Switzerland

It is estimated that around 230 people die each year due to radon (222Rn) exposure in Switzerland. 222Rn occurs mainly in closed environments like buildings and originates primarily from the subjacent ground. Therefore it depends strongly on geology and shows substantial regional variations. Correct identification of these regional variations would lead to substantial reduction of 222Rn exposure of the population based on appropriate construction of new and mitigation of already existing buildings. Prediction of indoor 222Rn concentrations (IRC) and identification of 222Rn prone areas is however difficult since IRC depend on a variety of different variables like building characteristics, meteorology, geology and anthropogenic factors. The present work aims at the development of predictive models and the understanding of IRC in Switzerland, taking into account a maximum of information in order to minimize the prediction uncertainty. The predictive maps will be used as a decision-support tool for 222Rn risk management. The construction of these models is based on different data-driven statistical methods, in combination with geographical information systems (GIS). In a first phase we performed univariate analysis of IRC for different variables, namely the detector type, building category, foundation, year of construction, the average outdoor temperature during measurement, altitude and lithology. All variables showed significant associations to IRC. Buildings constructed after 1900 showed significantly lower IRC compared to earlier constructions. We observed a further drop of IRC after 1970. In addition to that, we found an association of IRC with altitude. With regard to lithology, we observed the lowest IRC in sedimentary rocks (excluding carbonates) and sediments and the highest IRC in the Jura carbonates and igneous rock. The IRC data was systematically analyzed for potential bias due to spatially unbalanced sampling of measurements. In order to facilitate the modeling and the interpretation of the influence of geology on IRC, we developed an algorithm based on k-medoids clustering which permits to define coherent geological classes in terms of IRC. We performed a soil gas 222Rn concentration (SRC) measurement campaign in order to determine the predictive power of SRC with respect to IRC. We found that the use of SRC is limited for IRC prediction. The second part of the project was dedicated to predictive mapping of IRC using models which take into account the multidimensionality of the process of 222Rn entry into buildings. We used kernel regression and ensemble regression tree for this purpose. We could explain up to 33% of the variance of the log transformed IRC all over Switzerland. This is a good performance compared to former attempts of IRC modeling in Switzerland. As predictor variables we considered geographical coordinates, altitude, outdoor temperature, building type, foundation, year of construction and detector type. Ensemble regression trees like random forests allow to determine the role of each IRC predictor in a multidimensional setting. We found spatial information like geology, altitude and coordinates to have stronger influences on IRC than building related variables like foundation type, building type and year of construction. Based on kernel estimation we developed an approach to determine the local probability of IRC to exceed 300 Bq/m3. In addition to that we developed a confidence index in order to provide an estimate of uncertainty of the map. All methods allow an easy creation of tailor-made maps for different building characteristics. Our work is an essential step towards a 222Rn risk assessment which accounts at the same time for different architectural situations as well as geological and geographical conditions. For the communication of 222Rn hazard to the population we recommend to make use of the probability map based on kernel estimation. The communication of 222Rn hazard could for example be implemented via a web interface where the users specify the characteristics and coordinates of their home in order to obtain the probability to be above a given IRC with a corresponding index of confidence. Taking into account the health effects of 222Rn, our results have the potential to substantially improve the estimation of the effective dose from 222Rn delivered to the Swiss population.

Real-time MR imaging of myocardial regional function using strain-encoding (SENC) with tissue through-plane motion tracking.

PURPOSE: To implement real-time myocardial strain-encoding (SENC) imaging in combination with tracking the tissue displacement in the through-plane direction. MATERIALS AND METHODS: SENC imaging was combined with the slice-following technique by implementing three-dimensional (3D) selective excitation. Certain adjustments were implemented to reduce scan time to one heartbeat. A total of 10 volunteers and five pigs were scanned on a 3T MRI scanner. Spatial modulation of magnetization (SPAMM)-tagged images were acquired on planes orthogonal to the SENC planes for comparison. Myocardial infarction (MI) was induced in two pigs and the resulting SENC images were compared to standard delayed-enhancement (DE) images. RESULTS: The strain values computed from SENC imaging with slice-following showed significant difference from those acquired without slice-following, especially during systole (P < 0.01). The strain curves computed from the SENC images with and without slice-following were similar to those computed from the orthogonal SPAMM images, with and without, respectively, tracking the tag line displacement in the strain direction. The resulting SENC images showed good agreement with the DE images in identifying MI in infarcted pigs. CONCLUSION: Correction of through-plane motion in real-time cardiac functional imaging is feasible using slice-following. The strain measurements are more accurate than conventional SENC measurements in humans and animals, as validated with conventional MRI tagging.

Improved predictive mapping of indoor radon concentrations using ensemble regression trees based on automatic clustering of geological units.

PURPOSE: According to estimations around 230 people die as a result of radon exposure in Switzerland. This public health concern makes reliable indoor radon prediction and mapping methods necessary in order to improve risk communication to the public. The aim of this study was to develop an automated method to classify lithological units according to their radon characteristics and to develop mapping and predictive tools in order to improve local radon prediction. METHOD: About 240 000 indoor radon concentration (IRC) measurements in about 150 000 buildings were available for our analysis. The automated classification of lithological units was based on k-medoids clustering via pair-wise Kolmogorov distances between IRC distributions of lithological units. For IRC mapping and prediction we used random forests and Bayesian additive regression trees (BART). RESULTS: The automated classification groups lithological units well in terms of their IRC characteristics. Especially the IRC differences in metamorphic rocks like gneiss are well revealed by this method. The maps produced by random forests soundly represent the regional difference of IRCs in Switzerland and improve the spatial detail compared to existing approaches. We could explain 33% of the variations in IRC data with random forests. Additionally, the influence of a variable evaluated by random forests shows that building characteristics are less important predictors for IRCs than spatial/geological influences. BART could explain 29% of IRC variability and produced maps that indicate the prediction uncertainty. CONCLUSION: Ensemble regression trees are a powerful tool to model and understand the multidimensional influences on IRCs. Automatic clustering of lithological units complements this method by facilitating the interpretation of radon properties of rock types. This study provides an important element for radon risk communication. Future approaches should consider taking into account further variables like soil gas radon measurements as well as more detailed geological information.

The fourth extracellular loop of the alpha subunit of Na,K-ATPase. Functional evidence for close proximity with the second extracellular loop.

Na,K-ATPase is the main active transport system that maintains the large gradients of Na(+) and K(+) across the plasma membrane of animal cells. The crystal structure of a K(+)-occluding conformation of this protein has been recently published, but the movements of its different domains allowing for the cation pumping mechanism are not yet known. The structure of many more conformations is known for the related calcium ATPase SERCA, but the reliability of homology modeling is poor for several domains with low sequence identity, in particular the extracellular loops. To better define the structure of the large fourth extracellular loop between the seventh and eighth transmembrane segments of the alpha subunit, we have studied the formation of a disulfide bond between pairs of cysteine residues introduced by site-directed mutagenesis in the second and the fourth extracellular loop. We found a specific pair of cysteine positions (Y308C and D884C) for which extracellular treatment with an oxidizing agent inhibited the Na,K pump function, which could be rapidly restored by a reducing agent. The formation of the disulfide bond occurred preferentially under the E2-P conformation of Na,K-ATPase, in the absence of extracellular cations. Using recently published crystal structure and a distance constraint reproducing the existence of disulfide bond, we performed an extensive conformational space search using simulated annealing and showed that the Tyr(308) and Asp(884) residues can be in close proximity, and simultaneously, the SYGQ motif of the fourth extracellular loop, known to interact with the extracellular domain of the beta subunit, can be exposed to the exterior of the protein and can easily interact with the beta subunit.

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.

Garnet growth in the Nufenen Pass area (Swiss Alps)

Abstract : Textural division of a mineral in pyramids, with their apices located at the centre of the mineral and their bases corresponding to the mineral faces is called textural sector zoning. Textural sector zoning is observed in many metamorphic minerals like andalousite and garnet. Garnets found in the graphite rich black shales of the Mesozoic cover of the Gotthard Massif display textural sector zoning. The morphology of this sector zoning is not the same in different types of black shales observed in the Nufenen pass area. Garnets in foliated black shales display a well developed sector zoning while garnets found in cm-scale layered black shales display well developed sectors in the direction of the schistosity plane. This sector zoning is always associated with up to 30μm sized birefringent lamellae emanating radial from the sector boundaries. They alternate with isotrope lamellae. The garnet forming reaction was determined using singular value decomposition approach and results compared to thermodynamic calculations. It is of the form chl + mu + cc + cld = bt + fds + ank + gt + czo and is similar in both layered and foliated black shales. The calculated X(O) is close to 0.36 and does not significantly vary during the metamorphic history of the rock. This corresponds to X CO2, X CH4, and X H2O BSE imaging of garnets on oriented-cuts revealed that the orientation of the lamellae found within the sectors is controlled by crystallography. BSE imaging and electron microprobe analysis revealed that these lamellae are calcium rich compared to the isotropic lamellae. The addition of Ca to an almandine rich garnet causes a small distortion of the X site and potentially, ordering. Ordered and disordered garnet might have very similar free energies for this composition. Hence, two garnets with different composition can be precipitated with minor overstepping of the reaction. It is enough that continued nucleation of a new garnet layer slightly prefers the same structure to assure a fiber-like growth of both garnet compositions side by side. This hypothesis is in agreement with the thermodynamic properties of the garnet solid solution described in the literature and could explain the textures observed in garnets with these compositions. To understand the differences in sector zoning morphology, and crystal growth kinetics, crystal size distribution were determined in several samples using 2D spatial analysis of slab surfaces. The same nucleation rate law was chosen for all cases. Different growth rate law for non-layered black shales and layered black shales were used. Garnet in layered black shales grew according to a growth rate law of the form R=kt ½. The transport of nutrient is the limiting factor. Transport will occur preferentially on the schistosity planes. The shapes of the garnets in such rocks are therefore ovoid with the longest axis parallel to the schistosity planes. Sector zoning is less developed with sectors present only parallel to the schistosity planes. Garnet in non-layered blackshales grew according to a growth rate law of the form R=kt. The limiting factor is the attachment at the surface of the garnet. Garnets in these rocks will display a well developed sector zoning in all directions. The growth rate law is thus influenced by the texture of the rock. It favours or hinders the transport of nutrient to the mineral surface. Résumé : La zonation sectorielle texturale consiste en la division d'un cristal en pyramides dont les sommets sont localisés au centre du minéral. La base de ces pyramides correspond aux faces du minéral. Ce type de zonation est fréquemment observé dans les minéraux métamorphiques tels que l'andalousite ou le grenat. Les grenats présents dans les marnes riches en graphites de la couverture Mésozoïque du Massif du Gotthard présent une zonation sectorielle texturale. La morphologie de cette zonation n'est pas la même dans les marnes litées et dans les marnes foliées. Les grenats des marnes foliées montrent des secteurs bien développés dans 3 directions. Les grenats des marnes litées montrent des secteurs développés uniquement dans la direction des plans de schistosité. Cette zonation sectorielle est toujours associée à des lamelles biréfringentes de quelques microns de large qui partent de la limite des secteurs et qui sont perpendiculaires aux faces du grenat. Ces lamelles alternent avec des lamelles isotropes. La réaction de formation du grenat a été déterminée par calcul matriciel et thermodynamique. La réaction est de la forme chl + mu + cc + cld= bt + fds + ank + gt + czo. Elle est similaire dans les roches litées et dans les roches foliées. L'évaluation des conditions fluides montrent que le X(O) est proche de 0.36 et ne change pas de façon significative durant l'histoire métamorphique de la roche. Des images BSE sur des coupes orientées ont révélé que l'orientation de lamelles biréfringentes est contrôlée parla crystallographie. La comparaison des analyses à la microsonde électronique et des images BSE révèle également que les lamelles biréfringentes sont plus riches en calcium que les lamelles isotropes. L'addition de calcium va déformer légèrement le site X et ainsi créer un ordre sur ce site. L'énergie interne d'un grenat ordré et d'un grenat désordonné sont suffisamment proches pour qu'un léger dépassement de l'énergie de la réaction de formation permette la coexistence des 2 types de grenat dans le même minéral. La formation de lamelles est expliquée par le fait qu'un grenat préférera la même structure. Ces observations sont en accord avec la thermodynamique des solutions solides du grenat et permet d'expliquer les structures similaires observées dans des grenats provenant de lithologies différentes. Une étude de la distribution des tailles des grenats et une modélisation de la croissance a permis de mettre en évidence 2 mécanismes de croissance différents suivant la texture de la roche. Dans les 2 cas, la loi de nucléation est la même. Dans les roches litées, la loi de croissance est de forme R=kt½. Le transport des nutriments est le facteur limitant. Ce transport a lieu préférentiellement dans la direction des niveaux de schistosité. Les grenats ont une forme légèrement allongée car la croissance des secteurs est facilitée sur les niveaux de schistosité. La croissance des grenats dans les roches foliées suit une loi de croissance de la forme R=kt. Les seuls facteurs limitant la croissance sont les processus d'attachement à la surface du grenat. La loi de croissance de ces grenats est donc contrainte par la texture de la roche. Cela se marque par des différences dans la morphologie de la zonation sectorielle.

Separate neuronal and glial Na+,K+-ATPase isoforms regulate glucose utilization in response to membrane depolarization and elevated extracellular potassium.

The role of cell type-specific Na+,K+-ATPase isozymes in function-related glucose metabolism was studied using differentiated rat brain cell aggregate cultures. In mixed neuron-glia cultures, glucose utilization, determined by measuring the rate of radiolabeled 2-deoxyglucose accumulation, was markedly stimulated by the voltage-dependent sodium channel agonist veratridine (0.75 micromol/L), as well as by glutamate (100 micromol/L) and the ionotropic glutamate receptor agonist N-methyl-D-aspartate (NMDA) (10 micromol/L). Significant stimulation also was elicited by elevated extracellular potassium (12 mmol/L KCl), which was even more pronounced at 30 mmol/L KCl. In neuron-enriched cultures, a similar stimulation of glucose utilization was obtained with veratridine, specific ionotropic glutamate receptor agonists, and 30 mmol/L but not 12 mmol/L KCl. The effects of veratridine, glutamate, and NMDA were blocked by specific antagonists (tetrodotoxin, CNQX, or MK801, respectively). Low concentrations of ouabain (10(-6) mol/L) prevented stimulation by the depolarizing agents but reduced only partially the response to 12 mmol/L KCl. Together with previous data showing cell type-specific expression of Na+,K+-ATPase subunit isoforms in these cultures, the current results support the view that distinct isoforms of Na+,K+-ATPase regulate glucose utilization in neurons in response to membrane depolarization, and in glial cells in response to elevated extracellular potassium.

The fourth transmembrane segment of the Na,K-ATPase alpha subunit: a systematic mutagenesis study.

The Na,K-ATPase is a major ion-motive ATPase of the P-type family responsible for many aspects of cellular homeostasis. To determine the structure of the pathway for cations across the transmembrane portion of the Na,K-ATPase, we mutated 24 residues of the fourth transmembrane segment into cysteine and studied their function and accessibility by exposure to the sulfhydryl reagent 2-aminoethyl-methanethiosulfonate. Accessibility was also examined after treatment with palytoxin, which transforms the Na,K-pump into a cation channel. Of the 24 tested cysteine mutants, seven had no or a much reduced transport function. In particular cysteine mutants of the highly conserved "PEG" motif had a strongly reduced activity. However, most of the non-functional mutants could still be transformed by palytoxin as well as all of the functional mutants. Accessibility, determined as a 2-aminoethyl-methanethiosulfonate-induced reduction of the transport activity or as inhibition of the membrane conductance after palytoxin treatment, was observed for the following positions: Phe(323), Ile(322), Gly(326), Ala(330), Pro(333), Glu(334), and Gly(335). In accordance with a structural model of the Na,K-ATPase obtained by homology modeling with the two published structures of sarcoplasmic and endoplasmic reticulum calcium ATPase (Protein Data Bank codes 1EUL and 1IWO), the results suggest the presence of a cation pathway along the side of the fourth transmembrane segment that faces the space between transmembrane segments 5 and 6. The phenylalanine residue in position 323 has a critical position at the outer mouth of the cation pathway. The residues thought to form the cation binding site II ((333)PEGL) are also part of the accessible wall of the cation pathway opened by palytoxin through the Na,K-pump.