Registration and fusion of contrast-enhanced MRI myocardial substrate maps and X-ray angiograms

The aim of this work is to provide the necessary methods to register and fuse the endo-epicardial signal intensity (SI) maps extracted from contrast-enhanced magnetic resonance imaging (ceMRI) with X-ray coronary ngiograms using an intrinsic registrationbased algorithm to help pre-planning and guidance of catheterization procedures. Fusion of angiograms with SI maps was treated as a 2D-3D pose estimation, where each image point is projected to a Plücker line, and the screw representation for rigid motions is minimized using a gradient descent method. The resultant transformation is applied to the SI map that is then projected and fused on each angiogram. The proposed method was tested in clinical datasets from 6 patients with prior myocardial infarction. The registration procedure is optionally combined with an iterative closest point algorithm (ICP) that aligns the ventricular contours segmented from two ventriculograms.

On-line MRI sequences for the evaluation of apple internal quality.

On-line dynamic MRI, which is oriented to industrial grading lines, requires high-speed sequences with motion correction artefacts. In this study two different types of motion correction sequences have been used and have been implemented in real-time (FLASH and UFLARE). They are based on T2* and T2 respectively and their selection depends on the expected contrast effect of the disorder: while watercore enhances bright areas due to higher fluid mobility, internal breakdown potentiates low signal due to texture degradation. For watercore study, five different apple cultivars were used (Normanda-18-, Fuji-35-, Helada-36-, Verde Doncella-54-, Esperiega-75-) along two seasons (2011 and 2012). In total 218 fruits were measured under both, static conditions (20 slices per fruit) and under dynamic conditions (3 repetitions without slice selection). For internal breakdown, Braeburn cultivar has been studied (in total 106 fruits) under both static (20 slices per fruit) and dynamic conditions (3 replicates with slice selection). Metrological aspects such as repeatability of dynamic images and subsequent histogram feature stability become of major interest for further industrial application. Segregation ability among varying degrees of disorder is also analyzed.

Comparison of linear and non-linear 2D+T registration methods for DE-MRI cardiac perfusion studies

A series of motion compensation algorithms is run on the challenge data including methods that optimize only a linear transformation, or a non-linear transformation, or both – first a linear and then a non-linear transformation. Methods that optimize a linear transformation run an initial segmentation of the area of interest around the left myocardium by means of an independent component analysis (ICA) (ICA-*). Methods that optimize non-linear transformations may run directly on the full images, or after linear registration. Non-linear motion compensation approaches applied include one method that only registers pairs of images in temporal succession (SERIAL), one method that registers all image to one common reference (AllToOne), one method that was designed to exploit quasi-periodicity in free breathing acquired image data and was adapted to also be usable to image data acquired with initial breath-hold (QUASI-P), a method that uses ICA to identify the motion and eliminate it (ICA-SP), and a method that relies on the estimation of a pseudo ground truth (PG) to guide the motion compensation.

Image processing methods for human brain connectivity analysis from in-vivo diffusion MRI

The structural connectivity of the brain is considered to encode species-wise and subject-wise patterns that will unlock large areas of understanding of the human brain. Currently, diffusion MRI of the living brain enables to map the microstructure of tissue, allowing to track the pathways of fiber bundles connecting the cortical regions across the brain. These bundles are summarized in a network representation called connectome that is analyzed using graph theory. The extraction of the connectome from diffusion MRI requires a large processing flow including image enhancement, reconstruction, segmentation, registration, diffusion tracking, etc. Although a concerted effort has been devoted to the definition of standard pipelines for the connectome extraction, it is still crucial to define quality assessment protocols of these workflows. The definition of quality control protocols is hindered by the complexity of the pipelines under test and the absolute lack of gold-standards for diffusion MRI data. Here we characterize the impact on structural connectivity workflows of the geometrical deformation typically shown by diffusion MRI data due to the inhomogeneity of magnetic susceptibility across the imaged object. We propose an evaluation framework to compare the existing methodologies to correct for these artifacts including whole-brain realistic phantoms. Additionally, we design and implement an image segmentation and registration method to avoid performing the correction task and to enable processing in the native space of diffusion data. We release PySDCev, an evaluation framework for the quality control of connectivity pipelines, specialized in the study of susceptibility-derived distortions. In this context, we propose Diffantom, a whole-brain phantom that provides a solution to the lack of gold-standard data. The three correction methodologies under comparison performed reasonably, and it is difficult to determine which method is more advisable. We demonstrate that susceptibility-derived correction is necessary to increase the sensitivity of connectivity pipelines, at the cost of specificity. Finally, with the registration and segmentation tool called regseg we demonstrate how the problem of susceptibility-derived distortion can be overcome allowing data to be used in their original coordinates. This is crucial to increase the sensitivity of the whole pipeline without any loss in specificity.

Mapping strain exerted on blood vessel walls using deuterium double-quantum-filtered MRI

A technique is described for displaying distinct tissue layers of large blood vessel walls as well as measuring their mechanical strain. The technique is based on deuterium double-quantum-filtered (DQF) spectroscopic imaging. The effectiveness of the double-quantum filtration in suppressing the signal of bulk water is demonstrated on a phantom consisting of rat tail tendon fibers. Only intrafibrillar water is displayed, excluding all other signals of water molecules that reorient isotropically. One- and two-dimensional spectroscopic imaging of bovine aorta and coronary arteries show the characteristic DQF spectrum of each of the tissue layers. This property is used to obtain separate images of the outer layer, the tunica adventitia, or the intermediate layer, the tunica media, or both. To visualize the effect of elongation, the average residual quadrupole splitting <Δνq> is calculated for each pixel. Two-dimensional deuterium quadrupolar splitting images are obtained for a fully relaxed and a 55% elongated sample of bovine coronary artery. These images indicate that the strong effect of strain is associated with water molecules in the tunica adventitia whereas the DQF NMR signal of water in the tunica media is apparently strain-insensitive. After appropriate calibration, these average quadrupolar splitting images can be interpreted as strain maps.

Late-onset minor and major depression: early evidence for common neuroanatomical substrates detected by using MRI

The purpose of our study was to examine the neuroanatomical correlates of late-onset minor and major depression and to compare them with similar measures obtained from nondepressed controls. Our study groups were comprised of 18 patients with late-onset minor depression, 35 patients diagnosed with late-onset major depression, and 30 nondepressed controls. All subjects were scanned by using a 1.5-tesla MRI scanner. Absolute whole brain volume and normalized measures of prefrontal and temporal lobe volumes were obtained and used for comparison among groups. Our findings indicate that patients with minor depression present with specific neuroanatomical abnormalities that are comparable with the major depression group but significantly different from the controls. Normalized prefrontal lobe volumes show a significant linear trend with severity of depression, with volumes decreasing with illness severity. Whole brain volumes did not differ significantly among groups. These findings have broad implications for the biology of late-life depression and suggest that there may be common neurobiological substrates that underlie all clinically significant forms of late-onset mood disturbances.

Contributions of cell kill and posttreatment tumor growth rates to the repopulation of intracerebral 9L tumors after chemotherapy: An MRI study

The drought of progress in clinical brain tumor therapy provides an impetus for developing new treatments as well as methods for testing therapeutics in animal models. The inability of traditional assays to simultaneously measure tumor size, location, growth kinetics, and cell kill achieved by a treatment complicates the interpretation of therapy experiments in animal models. To address these issues, tumor volume measurements obtained from serial magnetic resonance images were used to noninvasively estimate cell kill values in individual rats with intracerebral 9L tumors after treatment with 0.5, 1, or 2 × LD10 doses of 1,3-bis(2-chloroethyl)-1-nitrosourea. The calculated cell kill values were consistently lower than those reported using traditional assays. A dose-dependent increase in 9L tumor doubling time after treatment was observed that significantly contributed to the time required for surviving cells to repopulate the tumor mass. This study reveals that increases in animal survival are not exclusively attributable to the fraction of tumor cells killed but rather are a function of the cell kill and repopulation kinetics, both of which vary with treatment dose.

Dynamic mapping at the laminar level of odor-elicited responses in rat olfactory bulb by functional MRI

We have applied functional MRI (fMRI) based on blood oxygenation level-dependent (BOLD) image-contrast to map odor-elicited olfactory responses at the laminar level in the rat olfactory bulb (OB) elicited by iso-amyl acetate (10−2 dilution of saturated vapor) with spatial and temporal resolutions of 220×220×1,000 μm and 36 s. The laminar structure of the OB was clearly depicted by high-resolution in vivo anatomical MRI with spatial resolution of 110×110×1,000 μm. In repeated BOLD fMRI measurements, highly significant (P < 0.001) foci were located in the outer layers of both OBs. The occurrence of focal OB activity within a domain at the level of individual glomeruli or groups of glomeruli was corroborated on an intra- and inter-animal basis under anesthetized conditions with this noninvasive method. The dynamic studies demonstrated that the odor-elicited BOLD activations were highly reproducible on a time scale of minutes, whereas over tens of minutes the activations sometimes varied slowly. We found large BOLD signal (ΔS/S = 10–30%) arising from the olfactory nerve layer, which is devoid of synapses and composed of unmyelinated fibers and glial cells. Our results support previous studies with other methods showing that odors elicit activity within glomerular layer domains in the mammalian OB, and extend the analysis to shorter time periods at the level of individual glomeruli or groups of glomeruli. With further improvement, BOLD fMRI should be ideal for systematic analysis of the functional significance of individual glomeruli in olfactory information encoding and of spatiotemporal processing within the olfactory system.

PRIM: Proximity imaging of green fluorescent protein-tagged polypeptides

We report a serendipitous discovery that extends the impressive catalog of reporter functions performed by green fluorescent protein (GFP) or its derivatives. When two GFP molecules are brought into proximity, changes in the relative intensities of green fluorescence emitted upon excitation at 395 vs. 475 nm result. These spectral changes provide a sensitive ratiometric index of the extent of self-association that can be exploited to quantitatively image homo-oligomerization or clustering processes of GFP-tagged proteins in vivo. The method, which we term proximity imaging (PRIM), complements fluorescence resonance energy transfer between a blue fluorescent protein donor and a GFP acceptor, a powerful method for imaging proximity relationships between different proteins. However, unlike fluorescence resonance energy transfer (which is a spectral interaction), PRIM depends on direct contact between two GFP modules, which can lead to structural perturbations and concomitant spectral changes within a module. Moreover, the precise spatial arrangement of the GFP molecules within a given dimer determines the magnitude and direction of the spectral change. We have used PRIM to detect FK1012-induced dimerization of GFP fused to FK506-binding protein and clustering of glycosylphosphatidylinositol-anchored GFP at cell surfaces.

Imaging of H217O distribution in the brain of a live rat by using proton-detected 17O MRI

Imaging of H217O has a number of important applications. Mapping the distribution of H217O produced by oxidative metabolism of 17O-enriched oxygen gas may lead to a new method of metabolic functional imaging; regional cerebral blood flow also can be measured by measuring the H217O distribution after the injection of 17O-enriched physiological saline solution. Previous studies have proposed a method for indirect detection of 17O. The method is based on the shortening of the proton T2 in H217O solutions, caused by the residual 17O-1H scalar coupling and transferred to the bulk water via fast chemical exchange. It has been shown that the proton T2 of H217O solutions can be restored to that of H216O by irradiating the resonance frequency of the 17O nucleus. The indirect 17O image thus is obtained by taking the difference between two T2-weighted spin-echo images: one acquired after irradiation of the 17O resonance and one acquired without irradiation. It also has been established that, at relatively low concentrations of H217O, the indirect method yields an image that quantitatively reflects the H217O distribution in the sample. The method is referred to as PRIMO (proton imaging of oxygen). In this work, we show in vivo proton images of the H217O distribution in a rat brain after an i.v. injection of H217O-enriched physiological saline solution. Implementing the indirect detection method in an echo-planar imaging sequence enabled obtaining H217O images with good spatial and temporal resolution of few seconds.

In vivo NMR and MRI using injection delivery of laser-polarized xenon

Because xenon NMR is highly sensitive to the local environment, laser-polarized xenon could be a unique probe of living tissues. Realization of clinical and medical science applications beyond lung airspace imaging requires methods of efficient delivery of laser-polarized xenon to tissues, because of the short spin-lattice relaxation times and relatively low concentrations of xenon attainable in the body. Preliminary results from the application of a polarized xenon injection technique for in vivo 129Xe NMR/MRI are extrapolated along with a simple model of xenon transit to show that the peak local concentration of polarized xenon delivered to tissues by injection may exceed that delivered by respiration by severalfold.

Transposon Tc1-derived, sequence-tagged sites in Caenorhabditis elegans as markers for gene mapping

We present an approach to map large numbers of Tc1 transposon insertions in the genome of Caenorhabditis elegans. Strains have been described that contain up to 500 polymorphic Tc1 insertions. From these we have cloned and shotgun sequenced over 2000 Tc1 flanks, resulting in an estimated set of 400 or more distinct Tc1 insertion alleles. Alignment of these sequences revealed a weak Tc1 insertion site consensus sequence that was symmetric around the invariant TA target site and reads CAYATATRTG. The Tc1 flanking sequences were compared with 40 Mbp of a C. elegans genome sequence. We found 151 insertions within the sequenced area, a density of ≈1 Tc1 insertion in every 265 kb. As the rest of the C. elegans genome sequence is obtained, remaining Tc1 alleles will fall into place. These mapped Tc1 insertions can serve two functions: (i) insertions in or near genes can be used to isolate deletion derivatives that have that gene mutated; and (ii) they represent a dense collection of polymorphic sequence-tagged sites. We demonstrate a strategy to use these Tc1 sequence-tagged sites in fine-mapping mutations.

Assessment and discrimination of odor stimuli in rat olfactory bulb by dynamic functional MRI

Dynamic blood oxygenation level-dependent functional MRI was applied at 7 T in the rat olfactory bulb (OB) with pulsed delivery of iso-amyl acetate (IAA) and limonene. Acquisition times for single-slice and whole OB data were 8 and 32 s, respectively, with spatial resolution of 220 × 220 × 250 μm. On an intrasubject basis, short IAA exposures of 0.6 min separated by 3.5-min intervals induced reproducible spatial activity patterns (SAPs) in the olfactory nerve layer, glomerular layer, and external plexiform layer. During long exposures (≈10 min), the initially dominant dorsal SAPs declined in intensity and area, whereas in some OB regions, the initially weak ventral/lateral SAPs increased first and then decreased. The SAPs of different concentrations were topologically similar, which implies that whereas an odor at various concentrations activates the same subsets of receptor cells, different concentrations are assessed and discriminated by variable magnitudes of laminarspecific activations. IAA and limonene reproducibly activated different subsets of receptor cells with some overlaps. Whereas qualitative topographical agreement was observed with results from other methods, the current dynamic blood oxygenation level-dependent functional MRI results can provide quantitative SAPs of the entire OB.

An approach to probe some neural systems interaction by functional MRI at neural time scale down to milliseconds

In this paper, we demonstrate an approach by which some evoked neuronal events can be probed by functional MRI (fMRI) signal with temporal resolution at the time scale of tens of milliseconds. The approach is based on the close relationship between neuronal electrical events and fMRI signal that is experimentally demonstrated in concurrent fMRI and electroencephalographic (EEG) studies conducted in a rat model with forepaw electrical stimulation. We observed a refractory period of neuronal origin in a two-stimuli paradigm: the first stimulation pulse suppressed the evoked activity in both EEG and fMRI signal responding to the subsequent stimulus for a period of several hundred milliseconds. When there was an apparent site–site interaction detected in the evoked EEG signal induced by two stimuli that were primarily targeted to activate two different sites in the brain, fMRI also displayed signal amplitude modulation because of the interactive event. With visual stimulation using two short pulses in the human brain, a similar refractory phenomenon was observed in activated fMRI signals in the primary visual cortex. In addition, for interstimulus intervals shorter than the known latency time of the evoked potential induced by the first stimulus (≈100 ms) in the primary visual cortex of the human brain, the suppression was not present. Thus, by controlling the temporal relation of input tasks, it is possible to study temporal evolution of certain neural events at the time scale of their evoked electrical activity by noninvasive fMRI methodology.

Mental chronometry using latency-resolved functional MRI

Vascular responses to neural activity are exploited as the basis of a number of brain imaging techniques. The vascular response is thought to be too slow to resolve the temporal sequence of events involved in cognitive tasks, and hence, imaging studies of mental chronometry have relied on techniques such as the evoked potential. Using rapid functional MRI (fMRI) of single trials of two simple behavioral tasks, we demonstrate that while the microvascular response to the onset of neural activity is delayed consistently by several seconds, the relative timing between the onset of the fMRI responses in different brain areas appears preserved. We examined a number of parameters that characterize the fMRI response and determined that its onset time is best defined by the inflection point from the resting baseline. We have found that fMRI onset latencies determined in this manner correlate well with independently measurable parameters of the tasks such as reaction time or stimulus presentation time and can be used to determine the origin of processing delays during cognitive or perceptual tasks with a temporal accuracy of tens of milliseconds and spatial resolution of millimeters.