Whole-brain MR-registered cryo-imaging of a porcine-human glioma model to compare contrast agent biodistributions.

As rapidly accelerating technology, fluorescence guided surgery (FGS) has the potential to place molecular information directly into the surgeon's field of view by imaging administered fluorescent contrast agents in real time, circumnavigating pre-operative MR registration challenges with brain deformation. The most successful implementation of FGS is 5-ALA-PpIX guided glioma resection which has been linked to improved patient outcomes. While FGS may offer direct in-field guidance, fluorescent contrast agent distributions are not as familiar to the surgical community as Gd-MRI uptake, and may provide discordant information from previous Gd-MRI guidance. Thus, a method to assess and validate consistency between fluorescence-labeled tumor regions and Gd-enhanced tumor regions could aid in understanding the correlation between optical agent fluorescence and Gd-enhancement. Herein, we present an approach for comparing whole-brain fluorescence biodistributions with Gd-enhancement patterns on a voxel-by-voxel basis using co-registered fluorescent cryo-volumes and Gd-MRI volumes. In this initial study, a porcine-human glioma xenograft model was administered 5-ALA-PpIX, imaged with MRI, and euthanized 22 hours following 5-ALA administration. Following euthanization, the extracted brain was imaged with the cryo-macrotome system. After image processing steps and non-rigid, point-based registration, the fluorescence cryo-volume and Gd-MRI volume were compared for similarity metrics including: image similarity, tumor shape similarity, and classification similarity. This study serves as a proof-of-principle in validating our screening approach for quantitatively comparing 3D biodistributions between optical agents and Gd-based agents.

Considerations for NIR-I and short-wave infrared (SWIR) fluorescence imaging within a clinical operating room.

Short-wave infrared (SWIR/NIR-II) fluorescence imaging has received increased attention for use in fluorescence-guided surgery (FGS) due to the potential for higher resolution imaging of subsurface structures and reduced autofluorescence compared to conventional NIR-I imaging. As with any fluorescence imaging modality introduced in the operating room, an appropriate accounting of contaminating background signal from other light sources in the operating room is an important step. Herein, we report the background signals in the SWIR and NIR-I emitted from commonly-used equipment in the OR, such as ambient and operating lights, LCD screens and surgical guidance systems. These results can guide implementation of protocols to reduce background signal.

Small dense LDL cholesterol in human subjects with different chronic inflammatory diseases.

BACKGROUND AND AIMS: Chronic inflammatory diseases (CID) are associated with a profound increase in cardiovascular (CV) risk resulting in reduced life expectancy. However, LDL-cholesterol is reported to be low in CID patients which is referred to as the "LDL paradoxon". The aim of the present study was to investigate whether LDL-particles in CID exhibit an increased content of the highly atherogenic small-dense LDL subfraction (sdLDL). METHODS AND RESULTS: In this prospective, single center, observational study we enrolled 141 patients with CID (RA n = 59, inflammatory bowel disease (IBD) n = 35, ankylosing spondylitis (SpA) n = 25, Psoriasis n = 22) in 2011 through 2013 to evaluate sdLDL levels before as well as 6 and 26 weeks after initiation of different anti-cytokine therapies (anti-TNFα, anti-IL-6R antibodies). sdLDL levels were compared to 141 healthy individuals in a case control design. Compared to healthy controls, all CID patients displayed a significantly higher sdLDL content within the LDL cholesterol fraction: RA 35.0 ± 9.2% (p < 0.001), SpA 42.5 ± 10.5% (p < 0.001), IBD 37.5 ± 7.1% (p < 0.001), Psoriasis 33.6 ± 4.6% (p < 0.01). Furthermore, the sdLDL/LDL ratio was significantly higher in male compared to female RA subjects (p < 0.05). Neither anti-TNFα nor anti-IL6R medication altered sdLDL levels despite a significant improvement of disease activity. CONCLUSION: In several different chronic inflammatory disease entities, LDL-cholesterol is shifted toward a pro-atherogenic phenotype due to an increased sdLDL content which might in part explain the LDL paradoxon. Since premature CV disease is a major burden of affected patients, specifically targeting lipid metabolism should be considered routinely in clinical patient care. CLINICAL TRIALS: Registration at German Clinical Trial Register (DRKS): DRKS00005285.

Distribution of Ixodes ricinus ticks and prevalence of tick-borne encephalitis virus among questing ticks in the Arctic Circle region of northern Norway.

This study investigated the geographical distribution of Ixodes ricinus and prevalence of the tick-borne encephalitis virus (TBEV) in northern Norway. Flagging for questing I. ricinus ticks was performed in areas ranging from Vikna in Nord-Trøndelag County, located 190km south of the Arctic Circle (66.3°N), to Steigen in Nordland County, located 155km north of the Arctic Circle. We found that ticks were abundant in both Vikna (64.5°N) and Brønnøy (65.1°N). Only a few ticks were found at locations ∼66°N, and no ticks were found at several locations up to 67.5°N. Real-time PCR (RT-PCR) analyses of the collected ticks (nymphs and adults) for the presence of TBEV revealed a low prevalence (0.1%) of TBEV among the nymphs collected in Vikna, while a prevalence of 0% to 3% was found among nymphs collected at five locations in Brønnøy. Adult ticks collected in Vikna and Brønnøy had higher rates of TBEV infection (8.6% and 0%-9.0%, respectively) than the nymphs. No evidence of TBEV was found in the few ticks collected further north of Brønnøy. This is the first report of TBEV being detected at locations up to 65.1°N. It remains to be verified whether viable populations of I. ricinus exist at locations north of 66°N. Future studies are warranted to increase our knowledge concerning tick distribution, tick abundance, and tick-borne pathogens in northern Norway.

Hyperspectral data processing improves PpIX contrast during fluorescence guided surgery of human brain tumors.

Fluorescence guided surgery (FGS) using aminolevulinic-acid (ALA) induced protoporphyrin IX (PpIX) provides intraoperative visual contrast between normal and malignant tissue during resection of high grade gliomas. However, maps of the PpIX biodistribution within the surgical field based on either visual perception or the raw fluorescence emissions can be masked by background signals or distorted by variations in tissue optical properties. This study evaluates the impact of algorithmic processing of hyperspectral imaging acquisitions on the sensitivity and contrast of PpIX maps. Measurements in tissue-simulating phantoms showed that (I) spectral fitting enhanced PpIX sensitivity compared with visible or integrated fluorescence, (II) confidence-filtering automatically determined the lower limit of detection based on the strength of the PpIX spectral signature in the collected emission spectrum (0.014-0.041 µg/ml in phantoms), and (III) optical-property corrected PpIX estimates were more highly correlated with independent probe measurements (r = 0.98) than with spectral fitting alone (r = 0.91) or integrated fluorescence (r = 0.82). Application to in vivo case examples from clinical neurosurgeries revealed changes to the localization and contrast of PpIX maps, making concentrations accessible that were not visually apparent. Adoption of these methods has the potential to maintain sensitive and accurate visualization of PpIX contrast over the course of surgery.

Viscoelastic power law parameters of in vivo human brain estimated by MR elastography.

The noninvasive imaging technique of magnetic resonance elastography (MRE) was used to estimate the power law behavior of the viscoelastic properties of the human brain in vivo. The mechanical properties for four volunteers are investigated using shear waves induced over a frequency range of 10-50Hz to produce a displacement field measured by magnetic resonance motion-encoding gradients. The average storage modulus (µR) reconstructed with non-linear inversion (NLI) increased from approximately 0.95 to 2.58kPa over the 10-50Hz span; the average loss modulus (µI) also increased from 0.29 to 1.25kPa over the range. These increases were modeled by independent power law (PL) relations for µR and µI returning whole brain, group mean exponent values of 0.88 and 1.07 respectively. Investigation of these exponents also showed distinctly different behavior in the region of the cerebral falx compared to other brain structures.

Suitability of poroelastic and viscoelastic mechanical models for high and low frequency MR elastography.

PURPOSE: Descriptions of the structure of brain tissue as a porous cellular matrix support application of a poroelastic (PE) mechanical model which includes both solid and fluid phases. However, the majority of brain magnetic resonance elastography (MRE) studies use a single phase viscoelastic (VE) model to describe brain tissue behavior, in part due to availability of relatively simple direct inversion strategies for mechanical property estimation. A notable exception is low frequency intrinsic actuation MRE, where PE mechanical properties are imaged with a nonlinear inversion algorithm. METHODS: This paper investigates the effect of model choice at each end of the spectrum of in vivo human brain actuation frequencies. Repeat MRE examinations of the brains of healthy volunteers were used to compare image quality and repeatability for each inversion model for both 50 Hz externally produced motion and ≈1 Hz intrinsic motions. Additionally, realistic simulated MRE data were generated with both VE and PE finite element solvers to investigate the effect of inappropriate model choice for ideal VE and PE materials. RESULTS: In vivo, MRE data revealed that VE inversions appear more representative of anatomical structure and quantitatively repeatable for 50 Hz induced motions, whereas PE inversion produces better results at 1 Hz. Reasonable VE approximations of PE materials can be derived by equating the equivalent wave velocities for the two models, provided that the timescale of fluid equilibration is not similar to the period of actuation. An approximation of the equilibration time for human brain reveals that this condition is violated at 1 Hz but not at 50 Hz. Additionally, simulation experiments when using the "wrong" model for the inversion demonstrated reasonable shear modulus reconstructions at 50 Hz, whereas cross-model inversions at 1 Hz were poor quality. Attenuation parameters were sensitive to changes in the forward model at both frequencies, however, no spatial information was recovered because the mechanisms of VE and PE attenuation are different. CONCLUSIONS: VE inversions are simpler with fewer unknown properties and may be sufficient to capture the mechanical behavior of PE brain tissue at higher actuation frequencies. However, accurate modeling of the fluid phase is required to produce useful mechanical property images at the lower frequencies of intrinsic brain motions.

3D parallel-detection microwave tomography for clinical breast imaging.

A biomedical microwave tomography system with 3D-imaging capabilities has been constructed and translated to the clinic. Updates to the hardware and reconfiguration of the electronic-network layouts in a more compartmentalized construct have streamlined system packaging. Upgrades to the data acquisition and microwave components have increased data-acquisition speeds and improved system performance. By incorporating analog-to-digital boards that accommodate the linear amplification and dynamic-range coverage our system requires, a complete set of data (for a fixed array position at a single frequency) is now acquired in 5.8 s. Replacement of key components (e.g., switches and power dividers) by devices with improved operational bandwidths has enhanced system response over a wider frequency range. High-integrity, low-power signals are routinely measured down to -130 dBm for frequencies ranging from 500 to 2300 MHz. Adequate inter-channel isolation has been maintained, and a dynamic range >110 dB has been achieved for the full operating frequency range (500-2900 MHz). For our primary band of interest, the associated measurement deviations are less than 0.33% and 0.5° for signal amplitude and phase values, respectively. A modified monopole antenna array (composed of two interwoven eight-element sub-arrays), in conjunction with an updated motion-control system capable of independently moving the sub-arrays to various in-plane and cross-plane positions within the illumination chamber, has been configured in the new design for full volumetric data acquisition. Signal-to-noise ratios (SNRs) are more than adequate for all transmit/receive antenna pairs over the full frequency range and for the variety of in-plane and cross-plane configurations. For proximal receivers, in-plane SNRs greater than 80 dB are observed up to 2900 MHz, while cross-plane SNRs greater than 80 dB are seen for 6 cm sub-array spacing (for frequencies up to 1500 MHz). We demonstrate accurate recovery of 3D dielectric property distributions for breast-like phantoms with tumor inclusions utilizing both the in-plane and new cross-plane data.

Multiresolution MR elastography using nonlinear inversion.

PURPOSE: Nonlinear inversion (NLI) in MR elastography requires discretization of the displacement field for a finite element (FE) solution of the "forward problem", and discretization of the unknown mechanical property field for the iterative solution of the "inverse problem". The resolution requirements for these two discretizations are different: the forward problem requires sufficient resolution of the displacement FE mesh to ensure convergence, whereas lowering the mechanical property resolution in the inverse problem stabilizes the mechanical property estimates in the presence of measurement noise. Previous NLI implementations use the same FE mesh to support the displacement and property fields, requiring a trade-off between the competing resolution requirements. METHODS: This work implements and evaluates multiresolution FE meshes for NLI elastography, allowing independent discretizations of the displacements and each mechanical property parameter to be estimated. The displacement resolution can then be selected to ensure mesh convergence, and the resolution of the property meshes can be independently manipulated to control the stability of the inversion. RESULTS: Phantom experiments indicate that eight nodes per wavelength (NPW) are sufficient for accurate mechanical property recovery, whereas mechanical property estimation from 50 Hz in vivo brain data stabilizes once the displacement resolution reaches 1.7 mm (approximately 19 NPW). Viscoelastic mechanical property estimates of in vivo brain tissue show that subsampling the loss modulus while holding the storage modulus resolution constant does not substantially alter the storage modulus images. Controlling the ratio of the number of measurements to unknown mechanical properties by subsampling the mechanical property distributions (relative to the data resolution) improves the repeatability of the property estimates, at a cost of modestly decreased spatial resolution. CONCLUSIONS: Multiresolution NLI elastography provides a more flexible framework for mechanical property estimation compared to previous single mesh implementations.

A spectrally constrained dual-band normalization technique for protoporphyrin IX quantification in fluorescence-guided surgery.

We report a dual-band normalization technique for in vivo quantification of the metabolic biomarker, protoporphyrin IX (PpIX), during brain tumor resection procedures. The accuracy of the approach was optimized in tissue simulating phantoms with varying absorption and scattering properties, validated with fluorimetric assessments on ex vivo brain tissue, and tested on human data acquired in vivo during fluorescence-guided surgery of brain tumors. The results demonstrate that the dual-band normalization technique allows PpIX concentrations to be accurately quantified by correction with reflectance data recorded and integrated within only two narrow wavelength intervals. The simplicity of the method lends itself to the enticing prospect that the method could be applicable to wide-field applications in quantitative fluorescence imaging and dosimetry in photodynamic therapy.

Analytic expression of fluorescence ratio detection correlates with depth in multi-spectral sub-surface imaging.

Here we derived analytical solutions to diffuse light transport in biological tissue based on spectral deformation of diffused near-infrared measurements. These solutions provide a closed-form mathematical expression which predicts that the depth of a fluorescent molecule distribution is linearly related to the logarithm of the ratio of fluorescence at two different wavelengths. The slope and intercept values of the equation depend on the intrinsic values of absorption and reduced scattering of tissue. This linear behavior occurs if the following two conditions are satisfied: the depth is beyond a few millimeters and the tissue is relatively homogeneous. We present experimental measurements acquired with a broad-beam non-contact multi-spectral fluorescence imaging system using a hemoglobin-containing diffusive phantom. Preliminary results confirm that a significant correlation exists between the predicted depth of a distribution of protoporphyrin IX molecules and the measured ratio of fluorescence at two different wavelengths. These results suggest that depth assessment of fluorescence contrast can be achieved in fluorescence-guided surgery to allow improved intra-operative delineation of tumor margins.

Subzone based magnetic resonance elastography using a Rayleigh damped material model.

PURPOSE: Recently, the attenuating behavior of soft tissue has been addressed in magnetic resonance elastography by the inclusion of a damping mechanism in the methods used to reconstruct the resulting mechanical property image. To date, this mechanism has been based on a viscoelastic model for material behavior. Rayleigh, or proportional, damping provides a more generalized model for elastic energy attenuation that uses two parameters to characterize contributions proportional to elastic and inertial forces. In the case of time-harmonic vibration, these two parameters lead to both the elastic modulus and the density being complex valued (as opposed to the case of pure viscoelasticity, where only the elastic modulus is complex valued). METHODS: This article presents a description of Rayleigh damping in the time-harmonic case, discussing the differences between this model and the viscoelastic damping models. In addition, the results from a subzone based Rayleigh damped elastography study of gelatin and tofu phantoms are discussed, along with preliminary results from in vivo breast data. RESULTS: Both the phantom and the tissue studies presented here indicate a change in the Rayleigh damping structure, described as Rayleigh composition, between different material types, with tofu and healthy tissue showing lower Rayleigh composition values than gelatin or cancerous tissue. CONCLUSIONS: It is possible that Rayleigh damping elastography and the concomitant Rayleigh composition images provide a mechanism for differentiating tissue structure in addition to measuring elastic stiffness and attenuation. Such information could be valuable in the use of Rayleigh damped magnetic resonance elastography as a diagnostic imaging tool.

Optical breast shape capture and finite-element mesh generation for electrical impedance tomography.

X-ray mammography is the standard for breast cancer screening. The development of alternative imaging modalities is desirable because mammograms expose patients to ionizing radiation. Electrical impedance tomography (EIT) may be used to determine tissue conductivity, a property which is an indicator of cancer presence. EIT is also a low-cost imaging solution and does not involve ionizing radiation. In breast EIT, impedance measurements are made using electrodes placed on the surface of the patient's breast. The complex conductivity of the volume of the breast is estimated by a reconstruction algorithm. EIT reconstruction is a severely ill-posed inverse problem. As a result, noisy instrumentation and incorrect modelling of the electrodes and domain shape produce significant image artefacts. In this paper, we propose a method that has the potential to reduce these errors by accurately modelling the patient breast shape. A 3D hand-held optical scanner is used to acquire the breast geometry and electrode positions. We develop methods for processing the data from the scanner and producing volume meshes accurately matching the breast surface and electrode locations, which can be used for image reconstruction. We demonstrate this method for a plaster breast phantom and a human subject. Using this approach will allow patient-specific finite-element meshes to be generated which has the potential to improve the clinical value of EIT for breast cancer diagnosis.

An octahedral shear strain-based measure of SNR for 3D MR elastography.

A signal-to-noise ratio (SNR) measure based on the octahedral shear strain (the maximum shear strain in any plane for a 3D state of strain) is presented for magnetic resonance elastography (MRE), where motion-based SNR measures are commonly used. The shear strain, γ, is directly related to the shear modulus, µ, through the definition of shear stress, τ = µÎ³. Therefore, noise in the strain is the important factor in determining the quality of motion data, rather than the noise in the motion. Motion and strain SNR measures were found to be correlated for MRE of gelatin phantoms and the human breast. Analysis of the stiffness distributions of phantoms reconstructed from the measured motion data revealed a threshold for both strain and motion SNR where MRE stiffness estimates match independent mechanical testing. MRE of the feline brain showed significantly less correlation between the two SNR measures. The strain SNR measure had a threshold above which the reconstructed stiffness values were consistent between cases, whereas the motion SNR measure did not provide a useful threshold, primarily due to rigid body motion effects.

Effects of frequency- and direction-dependent elastic materials on linearly elastic MRE image reconstructions.

The mechanical model commonly used in magnetic resonance elastography (MRE) is linear elasticity. However, soft tissue may exhibit frequency- and direction-dependent (FDD) shear moduli in response to an induced excitation causing a purely linear elastic model to provide an inaccurate image reconstruction of its mechanical properties. The goal of this study was to characterize the effects of reconstructing FDD data using a linear elastic inversion (LEI) algorithm. Linear and FDD phantoms were manufactured and LEI images were obtained from time-harmonic MRE acquisitions with variations in frequency and driving signal amplitude. LEI responses to artificially imposed uniform phase shifts in the displacement data from both purely linear elastic and FDD phantoms were also evaluated. Of the variety of FDD phantoms considered, LEI appeared to tolerate viscoelastic data-model mismatch better than deviations caused by poroelastic and anisotropic mechanical properties in terms of visual image contrast. However, the estimated shear modulus values were substantially incorrect relative to independent mechanical measurements even in the successful viscoelastic cases and the variations in mean values with changes in experimental conditions associated with uniform phase shifts, driving signal frequency and amplitude were unpredictable. Overall, use of LEI to reconstruct data acquired in phantoms with FDD material properties provided biased results under the best conditions and significant artifacts in the worst cases. These findings suggest that the success with which LEI is applied to MRE data in tissue will depend on the underlying mechanical characteristics of the tissues and/or organs systems of clinical interest.

Time-harmonic magnetic resonance elastography of the normal feline brain.

Imaging of the mechanical properties of in vivo brain tissue could eventually lead to non-invasive diagnosis of hydrocephalus, Alzheimer's disease and other pathologies known to alter the intracranial environment. The purpose of this work is to (1) use time-harmonic magnetic resonance elastography (MRE) to estimate the mechanical property distribution of cerebral tissue in the normal feline brain and (2) compare the recovered properties of grey and white matter. Various in vivo and ex vivo brain tissue property measurement strategies have led to the highly variable results that have been reported in the literature. MR elastography is an imaging technique that can estimate mechanical properties of tissue non-invasively and in vivo. Data was acquired in 14 felines and elastic parameters were estimated using a globo-regional nonlinear image reconstruction algorithm. Results fell within the range of values reported in the literature and showed a mean shear modulus across the subject group of 7-8 kPa with all but one animal falling within 5-15 kPa. White matter was statistically stiffer (p<0.01) than grey matter by about 1 kPa on a per subject basis. To the best of our knowledge, the results reported represent the most extensive set of estimates in the in vivo brain which have been based on MRE acquisition of the three-dimensional displacement field coupled to volumetric shear modulus image reconstruction achieved through nonlinear parameter estimation. However, the inter-subject variation in mean shear modulus indicates the need for further study, including the possibility of applying more advanced models to estimate the relevant tissue mechanical properties from the data.

Sensitivity study of an ultrasound coupled transrectal electrical impedance tomography system for prostate imaging.

In 2009, prostate cancer ranked as the most common cancer and the second most fatal cancer in men in the United States. Unfortunately, the current clinical diagnostic methods (e.g. prostate-specific antigen (PSA), digital rectal examination, endorectal MRI, transrectal ultrasound, biopsy) used for detecting and staging prostate cancer are limited. It has been shown that cancerous prostate tissue has significantly different electrical properties when compared to benign tissues. Based on these electrical property findings, a transrectal electrical impedance tomography (TREIT) system is proposed as a novel prostate imaging modality. The TREIT system comprises an array of electrodes interfaced with a clinical transrectal ultrasound (TRUS) probe. We evaluate this imaging system through a series of phantom imaging experiments to assess the system's ability to image high and low contrast objects at various positions. We found that the TREIT system can easily discern high contrast inclusions of 1 cm in diameter at distances centered at two times the radius of the TREIT probe away from the probe surface. Furthermore, this technology's ability to detect low contrast inclusions suggests that it has the potential to successfully detect prostate cancer.

Electrical impedance tomography reconstruction for three-dimensional imaging of the prostate.

Transrectal electrical impedance tomography (TREIT) has been proposed as an adjunct modality for enhancing standard clinical ultrasound (US) imaging of the prostate. The proposed TREIT probe has an array of electrodes adhered to the surface of a cylindrical US probe that is introduced inside of the imaging volume. Reconstructing TREIT images in the open-domain geometry established with this technique poses additional challenges to those encountered with closed-domain geometries, present in more conventional EIT systems, because of the rapidly decaying current densities at increasing distances from the probe surface. We developed a finite element method (FEM)-based dual-mesh reconstruction algorithm which employs an interpolation scheme for linking a fine forward mesh with a coarse grid of pixels, used for conductivity estimation. Simulation studies using the developed algorithm demonstrate the feasibility of imaging moderately contrasting inclusions at distances of three times the probe radius from the probe surface and at multiple angles about the probe's axis. The large, dense FEM meshes used here require significant computational effort. We have optimized our reconstruction algorithm with multi-core processing hardware and efficient parallelized computational software packages to achieve a speedup of 9.3 times when compared to a more traditional Matlab-based, single CPU solution. The simulation findings and computational optimization provide a state-of-the-art reconstruction platform for use in further evaluating transrectal electrical impedance tomography.

Monitoring of hemodynamic changes induced in the healthy breast through inspired gas stimuli with MR-guided diffuse optical imaging.

PURPOSE: The modulation of tissue hemodynamics has important clinical value in medicine for both tumor diagnosis and therapy. As an oncological tool, increasing tissue oxygenation via modulation of inspired gas has been proposed as a method to improve cancer therapy and determine radiation sensitivity. As a radiological tool, inducing changes in tissue total hemoglobin may provide a means to detect and characterize malignant tumors by providing information about tissue vascular function. The ability to change and measure tissue hemoglobin and oxygenation concentrations in the healthy breast during administration of three different types of modulated gas stimuli (oxygen/ carbogen, air/carbogen, and air/oxygen) was investigated. METHODS: Subjects breathed combinations of gases which were modulated in time. MR-guided diffuse optical tomography measured total hemoglobin and oxygen saturation in the breast every 30 s during the 16 min breathing stimulus. Metrics of maximum correlation and phase lag were calculated by cross correlating the measured hemodynamics with the stimulus. These results were compared to an air/air control to determine the hemodynamic changes compared to the baseline physiology. RESULTS: This study demonstrated that a gas stimulus consisting of alternating oxygen/carbogen induced the largest and most robust hemodynamic response in healthy breast parenchyma relative to the changes that occurred during the breathing of room air. This stimulus caused increases in total hemoglobin and oxygen saturation during the carbogen phase of gas inhalation, and decreases during the oxygen phase. These findings are consistent with the theory that oxygen acts as a vasoconstrictor, while carbogen acts as a vasodilator. However, difficulties in inducing a consistent change in tissue hemoglobin and oxygenation were observed because of variability in intersubject physiology, especially during the air/oxygen or air/carbogen modulated breathing protocols. CONCLUSIONS: MR-guided diffuse optical imaging is a unique tool that can measure tissue hemodynamics in the breast during modulated breathing. This technique may have utility in determining the therapeutic potential of pretreatment tissue oxygenation or in investigating vascular function. Future gas modulation studies in the breast should use a combination of oxygen and carbogen as the functional stimulus. Additionally, control measures of subject physiology during air breathing are critical for robust measurements.

Sensitivity study and optimization of a 3D electric impedance tomography prostate probe.

In current clinical practice, the primary diagnostic method for testing for prostate cancer is ultrasound-guided biopsy. In this paper, we consider using a sonolucent array of electrodes, printed on a thin Kapton layer and positioned on the imaging window of a transrectal ultrasound probe, as a method for providing coregistered electrical and ultrasound imaging of the prostate. As the electrical properties of malignant tissues have been shown to differ significantly from benign tissues, the estimation of the electrical properties is expected to be helpful in distinguishing certain beginning pathologies from cancer and in improving the detection rate that current biopsy methods provide. One of the main difficulties in estimating electrical properties of tissues with this electrode configuration is the rapid decay of the sensitivity with distance from the sensing array. In order to partially overcome this difficulty, we propose to use prior information from the ultrasound (US). Specifically we intend to delineate the boundaries of the prostate from the US, to subdivide the organ into a small number of voxels and to estimate the conductivity as constant on each of these subvolumes. We use a 3D forward model based on the finite element method for studying the sensitivity of a simulated segmented prostate for three different electrode array designs. The three designs present different electrode areas and inter-electrode gaps. Larger electrodes are desirable as they present a better contact, but we show that as they result in smaller inter-electrode gaps, shunting currents can be significant and the sensitivity is reduced. Because our clinical measurement system employs a single current source, we consider tetrapolar measurement patterns for evaluating these electrode configurations. Optimal measurement patterns are well defined for adaptive systems, where multiple currents are injected at the same time. For the electrode array designs we consider, which are three dimensional, there are no established systematic methods for forming sets of linearly independent tetrapolar measurement patterns. We develop a novel method for automatically computing a full set of independent tetrapolar measurement patterns that maximizes the sensitivity in a region of interest (ROI). We use these patterns in the forward modeling and sensitivity studies. In addition to the electrode arrays on the probe, we study the use of a further configuration, where a distal electrode is positioned on the exterior of the body and used for current injection.