3D probability tomography: theoretical developments and applications to high-resolution geophysical prospecting

The theory of the 3D multipole probability tomography method (3D GPT) to image source poles, dipoles, quadrupoles and octopoles, of a geophysical vector or scalar field dataset is developed. A geophysical dataset is assumed to be the response of an aggregation of poles, dipoles, quadrupoles and octopoles. These physical sources are used to reconstruct without a priori assumptions the most probable position and shape of the true geophysical buried sources, by determining the location of their centres and critical points of their boundaries, as corners, wedges and vertices. This theory, then, is adapted to the geoelectrical, gravity and self potential methods. A few synthetic examples using simple geometries and three field examples are discussed in order to demonstrate the notably enhanced resolution power of the new approach. At first, the application to a field example related to a dipole–dipole geoelectrical survey carried out in the archaeological park of Pompei is presented. The survey was finalised to recognize remains of the ancient Roman urban network including roads, squares and buildings, which were buried under the thick pyroclastic cover fallen during the 79 AD Vesuvius eruption. The revealed anomaly structures are ascribed to wellpreserved remnants of some aligned walls of Roman edifices, buried and partially destroyed by the 79 AD Vesuvius pyroclastic fall. Then, a field example related to a gravity survey carried out in the volcanic area of Mount Etna (Sicily, Italy) is presented, aimed at imaging as accurately as possible the differential mass density structure within the first few km of depth inside the volcanic apparatus. An assemblage of vertical prismatic blocks appears to be the most probable gravity model of the Etna apparatus within the first 5 km of depth below sea level. Finally, an experimental SP dataset collected in the Mt. Somma-Vesuvius volcanic district (Naples, Italy) is elaborated in order to define location and shape of the sources of two SP anomalies of opposite sign detected in the northwestern sector of the surveyed area. The modelled sources are interpreted as the polarization state induced by an intense hydrothermal convective flow mechanism within the volcanic apparatus, from the free surface down to about 3 km of depth b.s.l..

Diffusion of laser polarized gases in MRI

The use of Magnetic Resonance Imaging (MRI) as a diagnostic tool is increasingly employing functional contrast agents to study or contrast entire mechanisms. Contrast agents in MRI can be classified in two categories. One type of contrast agents alters the NMR signal of the protons in its surrounding, e.g. lowers the T1 relaxation time. The other type enhances the Nuclear Magnetic Resonance (NMR) signal of specific nuclei. For hyperpolarized gases the NMR signal is improved up to several orders of magnitude. However, gases have a high diffusivity which strongly influences the NMR signal strength, hence the resolution and appearance of the images. The most interesting question in spatially resolved experiments is of course the achievable resolution and contrast by controlling the diffusivity of the gas. The influence of such diffusive processes scales with the diffusion coefficient, the strength of the magnetic field gradients and the timings used in the experiment. Diffusion may not only limit the MRI resolution, but also distort the line shape of MR images for samples, which contain boundaries or diffusion barriers within the sampled space. In addition, due to the large polarization in gaseous 3He and 129Xe, spin diffusion (different from particle diffusion) could play a role in MRI experiments. It is demonstrated that for low temperatures some corrections to the NMR measured diffusion coefficient have to be done, which depend on quantum exchange effects for indistinguishable particles. Physically, if these effects can not change the spin current, they can do it indirectly by modifying the velocity distribution of the different spin states separately, so that the subsequent collisions between atoms and therefore the diffusion coefficient can eventually be affected. A detailed study of the hyperpolarized gas diffusion coefficient is presented, demonstrating the absence of spin diffusion (different from particle diffusion) influence in MRI at clinical conditions. A novel procedure is proposed to control the diffusion coefficient of gases in MRI by admixture of inert buffer gases. The experimental measured diffusion agrees with theoretical simulations. Therefore, the molecular mass and concentration enter as additional parameters into the equations that describe structural contrast. This allows for setting a structural threshold up to which structures contribute to the image. For MRI of the lung this allows for images of very small structural elements (alveoli) only, or in the other extreme, all airways can be displayed with minimal signal loss due to diffusion.

Galactic fountain flows in high-resolution hydrodynamical simulations of galaxy disks

Feedback from the most massive components of a young stellar cluster deeply affects the surrounding ISM driving an expanding over-pressured hot gas cavity in it. In spiral galaxies these structures may have sufficient energy to break the disk and eject large amount of material into the halo. The cycling of this gas, which eventually will fall back onto the disk, is known as galactic fountains. We aim at better understanding the dynamics of such fountain flow in a Galactic context, frame the problem in a more dynamic environment possibly learning about its connection and regulation to the local driving mechanism and understand its role as a metal diffusion channel. The interaction of the fountain with a hot corona is hereby analyzed, trying to understand the properties and evolution of the extraplanar material. We perform high resolution hydrodynamical simulations with the moving-mesh code AREPO to model the multi-phase ISM of a Milky Way type galaxy. A non-equilibrium chemical network is included to self consistently follow the evolution of the main coolants of the ISM. Spiral arm perturbations in the potential are considered so that large molecular gas structures are able to dynamically form here, self shielded from the interstellar radiation field. We model the effect of SN feedback from a new-born stellar cluster inside such a giant molecular cloud, as the driving force of the fountain. Passive Lagrangian tracer particles are used in conjunction to the SN energy deposition to model and study diffusion of freshly synthesized metals. We find that both interactions with hot coronal gas and local ISM properties and motions are equally important in shaping the fountain. We notice a bimodal morphology where most of the ejected gas is in a cold $10^4$ K clumpy state while the majority of the affected volume is occupied by a hot diffuse medium. While only about 20\% of the produced metals stay local, most of them quickly diffuse through this hot regime to great scales.

Diffusion-weighted MR imaging of native and transplanted kidneys

Applications of diffusion-weighted (DW) magnetic resonance (MR) imaging outside the brain have gained increasing importance in recent years. Owing to technical improvements in MR imaging units and faster sequences, the need for noninvasive imaging without contrast medium administration, mainly in patients with renal insufficiency, can be met successfully by applying this technique. DW MR imaging is quantified by the apparent diffusion coefficient (ADC), which provides information on diffusion and perfusion simultaneously. By using a biexponential fitting process of the DW MR imaging data, these two entities can be separated, because this type of fitting process can serve as an estimate of both the perfusion fraction and the true diffusion coefficient. DW MR imaging can be applied for functional evaluation of the kidneys in patients with acute or chronic renal failure. Impairment of renal function is accompanied by a decreased ADC. Acute ureteral obstruction leads to perfusion and diffusion changes in the affected kidney, and renal artery stenosis results in a decreased ADC. In patients with pyelonephritis, diffuse or focal changes in signal intensity are seen on the high-b-value images, with increased signal intensity corresponding to low signal intensity on the ADC map. The feasibility and reproducibility of DW MR imaging in patients with transplanted kidneys have already been demonstrated, and initial results seem to be promising for the assessment of allograft deterioration. Overall, performance of renal DW MR imaging, presuming that measurements are of high quality, will further boost this modality, particularly for early detection of diffuse renal conditions, as well as more accurate characterization of focal renal lesions.

A multi-center study: Intra-scan and inter-scan variability of diffusion spectrum imaging

The objective of this study was to investigate whether it is possible to pool together diffusion spectrum imaging data from four different scanners, located at three different sites. Two of the scanners had identical configuration whereas two did not. To measure the variability, we extracted three scalar maps (ADC, FA and GFA) from the DSI and utilized a region and a tract-based analysis. Additionally, a phantom study was performed to rule out some potential factors arising from the scanner performance in case some systematic bias occurred in the subject study. This work was split into three experiments: intra-scanner reproducibility, reproducibility with twin-scanner settings and reproducibility with other configurations. Overall for the intra-scanner and twin-scanner experiments, the region-based analysis coefficient of variation (CV) was in a range of 1%-4.2% and below 3% for almost every bundle for the tract-based analysis. The uncinate fasciculus showed the worst reproducibility, especially for FA and GFA values (CV 3.7-6%). For the GFA and FA maps, an ICC value of 0.7 and above is observed in almost all the regions/tracts. Looking at the last experiment, it was found that there is a very high similarity of the outcomes from the two scanners with identical setting. However, this was not the case for the two other imagers. Given the fact that the overall variation in our study is low for the imagers with identical settings, our findings support the feasibility of cross-site pooling of DSI data from identical scanners.

High interstitial fluid pressure is associated with low tumour penetration of diagnostic monoclonal antibodies applied for molecular imaging purposes

The human epithelial cell adhesion molecule (EpCAM) is highly expressed in a variety of clinical tumour entities. Although an antibody against EpCAM has successfully been used as an adjuvant therapy in colon cancer, this therapy has never gained wide-spread use. We have therefore investigated the possibilities and limitations for EpCAM as possible molecular imaging target using a panel of preclinical cancer models. Twelve human cancer cell lines representing six tumour entities were tested for their EpCAM expression by qPCR, flow cytometry analysis and immunocytochemistry. In addition, EpCAM expression was analyzed in vivo in xenograft models for tumours derived from these cells. Except for melanoma, all cell lines expressed EpCAM mRNA and protein when grown in vitro. Although they exhibited different mRNA levels, all cell lines showed similar EpCAM protein levels upon detection with monoclonal antibodies. When grown in vivo, the EpCAM expression was unaffected compared to in vitro except for the pancreatic carcinoma cell line 5072 which lost its EpCAM expression in vivo. Intravenously applied radio-labelled anti EpCAM MOC31 antibody was enriched in HT29 primary tumour xenografts indicating that EpCAM binding sites are accessible in vivo. However, bound antibody could only be immunohistochemically detected in the vicinity of perfused blood vessels. Investigation of the fine structure of the HT29 tumour blood vessels showed that they were immature and prone for higher fluid flux into the interstitial space. Consistent with this hypothesis, a higher interstitial fluid pressure of about 12 mbar was measured in the HT29 primary tumour via "wick-in-needle" technique which could explain the limited diffusion of the antibody into the tumour observed by immunohistochemistry.

High-resolution computer simulation of the dynamics of isoelectric focusing using carrier ampholytes: focusing with concurrent electrophoretic mobilization is an isotachophoretic process

Focusing of four hemoglobins with concurrent electrophoretic mobilization was studied by computer simulation. A dynamic electrophoresis simulator was first used to provide a detailed description of focusing in a 100-carrier component, pH 6-8 gradient using phosphoric acid as anolyte and NaOH as catholyte. These results are compared to an identical simulation except that the catholyte contained both NaOH and NaCl. A stationary, steady-state distribution of carrier components and hemoglobins is produced in the first configuration. In the second, the chloride ion migrates into and through the separation space. It is shown that even under these conditions of chloride ion flux a pH gradient forms. All amphoteric species acquire a slight positive charge upon focusing and the whole pattern is mobilized towards the cathode. The cathodic gradient end is stable whereas the anodic end is gradually degrading due to the continuous accumulation of chloride. The data illustrate that the mobilization is a cationic isotachophoretic process with the sodium ion being the leading cation. The peak height of the hemoglobin zones decreases somewhat upon mobilization, but the zones retain a relatively sharp profile, thus facilitating detection. The electropherograms that would be produced by whole column imaging and by a single detector placed at different locations along the focusing column are presented and show that focusing can be commenced with NaCl present in the catholyte at the beginning of the experiment. However, this may require detector placement on the cathodic side of the catholyte/sample mixture interface.

Application of 3D documentation and geometric reconstruction methods in traffic accident analysis: with high resolution surface scanning, radiological MSCT/MRI scanning and real data based animation

The examination of traffic accidents is daily routine in forensic medicine. An important question in the analysis of the victims of traffic accidents, for example in collisions between motor vehicles and pedestrians or cyclists, is the situation of the impact. Apart from forensic medical examinations (external examination and autopsy), three-dimensional technologies and methods are gaining importance in forensic investigations. Besides the post-mortem multi-slice computed tomography (MSCT) and magnetic resonance imaging (MRI) for the documentation and analysis of internal findings, highly precise 3D surface scanning is employed for the documentation of the external body findings and of injury-inflicting instruments. The correlation of injuries of the body to the injury-inflicting object and the accident mechanism are of great importance. The applied methods include documentation of the external and internal body and the involved vehicles and inflicting tools as well as the analysis of the acquired data. The body surface and the accident vehicles with their damages were digitized by 3D surface scanning. For the internal findings of the body, post-mortem MSCT and MRI were used. The analysis included the processing of the obtained data to 3D models, determination of the driving direction of the vehicle, correlation of injuries to the vehicle damages, geometric determination of the impact situation and evaluation of further findings of the accident. In the following article, the benefits of the 3D documentation and computer-assisted, drawn-to-scale 3D comparisons of the relevant injuries with the damages to the vehicle in the analysis of the course of accidents, especially with regard to the impact situation, are shown on two examined cases.

FUNCTIONAL PRINCIPAL COMPONENTS MODEL FOR HIGH-DIMENSIONAL BRAIN IMAGING

We establish a fundamental equivalence between singular value decomposition (SVD) and functional principal components analysis (FPCA) models. The constructive relationship allows to deploy the numerical efficiency of SVD to fully estimate the components of FPCA, even for extremely high-dimensional functional objects, such as brain images. As an example, a functional mixed effect model is fitted to high-resolution morphometric (RAVENS) images. The main directions of morphometric variation in brain volumes are identified and discussed.

Cerebral connectivity and psychotic personality traits : A diffusion tensor imaging study

This study aims to investigate the relationship between regional connectivity in the brain white matter and the presence of psychotic personality traits, in healthy subjects with psychotic traits. Thirteen healthy controls were administered the MMPI-2, to assess psychotic traits and, according to MMPI results, a dichotomization into a group of "high-psychotic" and "low-psychotic" was performed. Diffusion tensor imaging (DTI) was used as a non-invasive measure, in order to obtain information about the fractional anisotropy (FA), an intravoxel index of local connectivity and, by means of a voxelwise approach, the between-group differences of the FA values were calculated. The "high-psychotic" group showed higher FA in the left arcuate fasciculus. Subjects with low scores for psychotic traits had significantly higher FA in the corpus callosum, right arcuate fasciculus, and fronto-parietal fibers. In line with previous brain imaging studies of schizophrenia spectrum disorders, our results suggest that psychotic personality traits are related to altered connectivity and brain asymmetry.

The amygdala in schizophrenia: a trimodal magnetic resonance imaging study

In schizophrenic psychoses, structural and functional alterations of the amygdala have been demonstrated by several neuroimaging studies. However, postmortem examinations on the brains of schizophrenics did not confirm the volume changes reported by volumetric magnetic resonance imaging (MRI) studies. In order to address these contradictory findings and to further elucidate the possibly underlying pathophysiological process of the amygdala, we employed a trimodal MRI design including high-resolution volumetry, diffusion tensor imaging (DTI), and quantitative magnetization transfer imaging (qMTI) in a sample of 14 schizophrenic patients and 14 matched controls. Three-dimensional MRI volumetry revealed a significant reduction of amygdala raw volumes in the patient group, while amygdala volumes normalized for intracranial volume did not differ between the two groups. The regional diffusional anisotropy of the amygdala, expressed as inter-voxel coherence (COH), showed a marked and significant reduction in schizophrenics. Assessment of qMTI parameters yielded significant group differences for the T2 time of the bound proton pool and the T1 time of the free proton pool, while the semi-quantitative magnetization transfer ratio (MTR) did not differ between the groups. The application of multimodal MRI protocols is diagnostically relevant for the differentiation between schizophrenic patients and controls and provides a new strategy for the detection and characterization of subtle structural alterations in defined regions of the living brain.

Characterization of white matter alterations in phenylketonuria by magnetic resonance relaxometry and diffusion tensor imaging

A multimodal MR study including relaxometry, diffusion tensor imaging (DTI), and MR spectroscopy was performed on patients with classical phenylketonuria (PKU) and matched controls, to improve our understanding of white matter (WM) lesions. Relaxometry yields information on myelin loss or malformation and may substantiate results from DTI attributed to myelin changes. Relaxometry was used to determine four brain compartments in normal-appearing brain tissue (NABT) and in lesions: water in myelin bilayers (myelin water, MW), water in gray matter (GM), water in WM, and water with long relaxation times (cerebrospinal fluid [CSF]-like signals). DTI yielded apparent diffusion coefficients (ADCs) and fractional anisotropies. MW and WM content were reduced in NABT and in lesions of PKU patients, while CSF-like signals were significantly increased. ADC values were reduced in PKU lesions, but also in the corpus callosum. Diffusion anisotropy was reduced in lesions because of a stronger decrease in the longitudinal than in the transverse diffusion. WM content and CSF-like components in lesions correlated with anisotropy and ADC. ADC values in lesions and in the corpus callosum correlated negatively with blood and brain phenylalanine (Phe) concentrations. Intramyelinic edema combined with vacuolization is a likely cause of the WM alterations. Correlations between diffusivity and Phe concentrations confirm vulnerability of WM to high Phe concentrations.

Electron microscopy of high pressure frozen samples: bridging the gap between cellular ultrastructure and atomic resolution

Transmission electron microscopy has provided most of what is known about the ultrastructural organization of tissues, cells, and organelles. Due to tremendous advances in crystallography and magnetic resonance imaging, almost any protein can now be modeled at atomic resolution. To fully understand the workings of biological "nanomachines" it is necessary to obtain images of intact macromolecular assemblies in situ. Although the resolution power of electron microscopes is on the atomic scale, in biological samples artifacts introduced by aldehyde fixation, dehydration and staining, but also section thickness reduces it to some nanometers. Cryofixation by high pressure freezing circumvents many of the artifacts since it allows vitrifying biological samples of about 200 mum in thickness and immobilizes complex macromolecular assemblies in their native state in situ. To exploit the perfect structural preservation of frozen hydrated sections, sophisticated instruments are needed, e.g., high voltage electron microscopes equipped with precise goniometers that work at low temperature and digital cameras of high sensitivity and pixel number. With them, it is possible to generate high resolution tomograms, i.e., 3D views of subcellular structures. This review describes theory and applications of the high pressure cryofixation methodology and compares its results with those of conventional procedures. Moreover, recent findings will be discussed showing that molecular models of proteins can be fitted into depicted organellar ultrastructure of images of frozen hydrated sections. High pressure freezing of tissue is the base which may lead to precise models of macromolecular assemblies in situ, and thus to a better understanding of the function of complex cellular structures.

High-resolution computer simulation of electrophoretic mobilization in isoelectric focusing

Cationic and anionic electrophoretic mobilization for focusing of hemoglobins (Hb's) in the presence of 100 carrier ampholytes covering a pI range of 6.00-7.98 was studied by computer simulation at a constant current density of 300 A/m(2). Electropherograms that would be produced by whole column imaging and by single detectors placed at different locations along the focusing column are presented. Upon mobilization, peak heights of the Hb zones decrease, but the zones retain a relatively sharp constant profile and are migrating at a constant velocity. A further peak decrease occurs during readjustment at the locations of the original buffer/column interfaces, indicating that detection sensitivity is the lowest at these locations. An anionic carrier ampholyte mobility smaller than that of its cationic species produces a cathodic drift which is smaller than the transport rate used for electrophoretic mobilization. Compared to the case with equal mobilities of carrier ampholyte species, a small increase (decrease) is predicted for the cationic (anionic) mobilization rate within the focusing column. Simulation data suggest that electrophoretic mobilization after focusing and focusing with concurrent electrophoretic mobilization are comparable isotachophoretic processes that occur when there is an uninterrupted flux of an ion through the focusing column. Cathodic drift caused by unequal mobilities of the species of carrier ampholytes, electrophoretic mobilization, and decomposition occurring at the pH gradient edges are related electrophoretic processes.

High-resolution computer simulation of the dynamics of isoelectric focusing: in quest of more realistic input parameters for carrier ampholytes

The impact of the systematic variation of either DeltapK(a) or mobility of 140 biprotic carrier ampholytes on the conductivity profile of a pH 3-10 gradient was studied by dynamic computer simulation. A configuration with the greatest DeltapK(a) in the pH 6-7 range and uniform mobilities produced a conductivity profile consistent with that which is experimentally observed. A similar result was observed when the neutral (pI = 7) ampholyte is assigned the lowest mobility and mobilities of the other carriers are systematically increased as their pI's recede from 7. When equal DeltapK(a) values and mobilities are assigned to all ampholytes a conductivity plateau in the pH 5-9 region is produced which does not reflect what is seen experimentally. The variation in DeltapK(a) values is considered to most accurately reflect the electrochemical parameters of commercially available mixtures of carrier ampholytes. Simulations with unequal mobilities of the cationic and anionic species of the carrier ampholytes show either cathodic (greater mobility of the cationic species) or anodic (greater mobility of the anionic species) drifts of the pH gradient. The simulated cationic drifts compare well to those observed experimentally in a capillary in which the focusing of three dyes was followed by whole column optical imaging. The cathodic drift flattens the acidic portion of the gradient and steepens the basic part. This phenomenon is an additional argument against the notion that focused zones of carrier ampholytes have no electrophoretic flux.