Influence of the liquid ionic strength on the resonance frequency and shell parameters of lipid-coated microbubbles.

The correct measurement of the resonance frequency and shell properties of coated microbubbles (MBs) is essential in understanding and optimizing their response to ultrasound (US) exposure parameters. In diagnostic and therapeutic ultrasound, MBs are typically surrounded by blood; however, the influence of the medium charges on the MB resonance frequency has not been systematically studied using controlled measurements. This study aims to measure the medium charge interactions on MB behavior by measuring the frequency-dependent attenuation of the same size MBs in mediums with different charge densities. In-house lipid-coated MBs with C3F8 gas core were formulated. The MBs were isolated to a mean size of 2.35 µm using differential centrifugation. MBs were diluted to ≈8×105 MBs/mL in distilled water (DW), Phosphate-Buffered Saline solution (PBS1x) and PBS10x. The frequency-dependent attenuation of the MBs solutions was measured using an aligned pair of PVDF transducers with a center frequency of 10MHz and 100% bandwidth in the linear oscillation regime (7 kPa pressure amplitude). The MB shell properties were estimated by fitting the linear equation to experiments. Using a pendant drop tension meter, the surface tension at the equilibrium of ≈6 mm diameter size drops of the same MB shell was measured inside DW, PBS1x and PBS10x. The surface tension at the C3F8/solution interface was estimated by fitting the Young-Laplace equation from the recorded images. The frequency of the peak attenuation at different salinity levels was 13, 7.5 and 6.25 MHz in DW, PBS1x and PBS-10x, respectively. The attenuation peak increased by ≈140% with increasing ion density. MBs' estimated shell elasticity decreased by 64% between DW and PBS-1x and 36% between PBS-1x and PBS-10x. The drop surface tension reduced by 10.5% between DW and PBS-1x and by 5% between PBS-1x and PBS-10x, respectively. Reduction in the shell stiffness is consistent with the drop surface tension measurements. The shell viscosity was reduced by ≈40% between DW and PBS-1x and 42% between PBS-1x and PBS-10x. The reduction in the fitted stiffness and viscosity is possibly due to the formation of a densely charged layer around the shell, further reducing the effective surface tension on the MBs. The changes in the resonance frequency and estimated shell parameters were significant and may potentially help to better understand and explain bubble behavior in applications.

Probing the pressure dependence of sound speed and attenuation in bubbly media: Experimental observations, a theoretical model and numerical calculations.

The problem of attenuation and sound speed of bubbly media has remained partially unsolved. Comprehensive data regarding pressure-dependent changes of the attenuation and sound speed of a bubbly medium are not available. Our theoretical understanding of the problem is limited to linear or semi-linear theoretical models, which are not accurate in the regime of large amplitude bubble oscillations. Here, by controlling the size of the lipid coated bubbles (mean diameter of ≈5.4µm), we report the first time observation and characterization of the simultaneous pressure dependence of sound speed and attenuation in bubbly water below, at and above microbubbles resonance (frequency range between 1-3 MHz). With increasing acoustic pressure (between 12.5-100 kPa), the frequency of the peak attenuation and sound speed decreases while maximum and minimum amplitudes of the sound speed increase. We propose a nonlinear model for the estimation of the pressure dependent sound speed and attenuation with good agreement with the experiments. The model calculations are validated by comparing with the linear and semi-linear models predictions. One of the major challenges of the previously developed models is the significant overestimation of the attenuation at the bubble resonance at higher void fractions (e.g. 0.005). We addressed this problem by incorporating bubble-bubble interactions and comparing the results to experiments. Influence of the bubble-bubble interactions increases with increasing pressure. Within the examined exposure parameters, we numerically show that, even for low void fractions (e.g. 5.1×10-6) with increasing pressure the sound speed may become 4 times higher than the sound speed in the non-bubbly medium.

On the threshold of 1/2 order subharmonic emissions in the oscillations of ultrasonically excited bubbles.

The pressure threshold for 1/2 order subharmonic (SH) emissions and period doubling during the oscillations of ultrasonically excited bubbles is thought to be minimum when the bubble is sonicated with twice its resonance frequency (fr). This estimate is based on studies that simplified or neglected the effects of thermal damping. In this work, the nonlinear dynamics of ultrasonically excited bubbles is investigated accounting for the thermal dissipation. Results are visualized using bifurcation diagrams as a function of pressure. Here we show that, and depending on the gas, the pressure threshold for 1/2 order SHs can be minimum at a frequency between 0.5fr≤f≤0.6fr. In this frequency range, the generation of 1/2 order SHs are due to the occurrence of 5/2 order ultra-harmonic resonance. The stability of such oscillations is size dependent. For an air bubble immersed in water, only bubbles bigger than 1 µm in diameter are able to emit non-destructive SHs in these frequency ranges.

Nonlinear dynamics and bifurcation structure of ultrasonically excited lipid coated microbubbles.

In many applications, microbubbles (MBs) are encapsulated by a lipid coating to increase their stability. However, the complex behavior of the lipid coating including buckling and rupture sophisticates the dynamics of the MBs and as a result the dynamics of the lipid coated MBs (LCMBs) are not well understood. Here, we investigate the nonlinear behavior of the LCMBs by analyzing their bifurcation structure as a function of acoustic pressure. We show that, the LC can enhance the generation of period 2 (P2), P3, higher order subharmonics (SH), superharmonics and chaos at very low excitation pressures (e.g. 1 kPa). For LCMBs sonicated by their SH resonance frequency and in line with experimental observations with increasing pressure, P2 oscillations exhibit three stages: generation at low acoustic pressures, disappearance and re-generation. Within non-destructive oscillation regimes and by pressure amplitude increase, LCMBs can also exhibit two saddle node (SN) bifurcations resulting in possible abrupt enhancement of the scattered pressure. The first SN resembles the pressure dependent resonance phenomenon in uncoated MBs and the second SN resembles the pressure dependent SH resonance. Depending on the initial surface tension of the LCMBs, the nonlinear behavior may also be suppressed for a wide range of excitation pressures.

Nonlinear power loss in the oscillations of coated and uncoated bubbles: Role of thermal, radiation and encapsulating shell damping at various excitation pressures.

This study presents the fundamental equations governing the pressure dependent disipation mechanisms in the oscillations of coated bubbles. A simple generalized model (GM) for coated bubbles accounting for the effect of compressibility of the liquid is presented. The GM was then coupled with nonlinear ODEs that account for the thermal effects. Starting with mass and momentum conservation equations for a bubbly liquid and using the GM, nonlinear pressure dependent terms were derived for power dissipation due to thermal damping (Td), radiation damping (Rd) and dissipation due to the viscosity of liquid (Ld) and coating (Cd). The pressure dependence of the dissipation mechanisms of the coated bubble have been analyzed. The dissipated energies were solved for uncoated and coated 2-20 µm in bubbles over a frequency range of 0.25fr-2.5fr (fr is the bubble resonance) and for various acoustic pressures (1 kPa-300 kPa). Thermal effects were examined for air and C3F8 gas cores. In the case of air bubbles, as pressure increases, the linear thermal model looses accuracy and accurate modeling requires inclusion of the full thermal model. However, for coated C3F8 bubbles of diameter 1-8 µm, which are typically used in medical ultrasound, thermal effects maybe neglected even at higher pressures. For uncoated bubbles, when pressure increases, the contributions of Rd grow faster and become the dominant damping mechanism for pressure dependent resonance frequencies (e.g. fundamental and super harmonic resonances). For coated bubbles, Cd is the strongest damping mechanism. As pressure increases, Rd contributes more to damping compared to Ld and Td. For coated bubbles, the often neglected compressibility of the liquid has a strong effect on the oscillations and should be incorporated in models. We show that the scattering to damping ratio (STDR), a measure of the effectiveness of the bubble as contrast agent, is pressure dependent and can be maximized for specific frequency ranges and pressures.

Critical corrections to models of nonlinear power dissipation of ultrasonically excited bubbles.

Current models for calculating nonlinear power dissipation during the oscillations of acoustically excited bubbles generate non-physical values for the radiation damping (Drd) term for some frequency and pressure regions that include near resonance oscillations. Moreover, the ratio of the dissipated powers significantly deviate from the values that are calculated by the linear model at low amplitude oscillations (acoustic excitation pressure of PA=1 kPa and expansion ratio of <≊1.01). In high amplitude oscillation regimes (Pa⩾20 kPa), the dissipated power due to Drd deviates largely from the dissipated power as calculated by the widely accepted approach that uses the scattered power by the bubbles. We provide critical corrections to the present models. The validity of the results was examined in regimes of low amplitude oscillations and high amplitude oscillations. In the low amplitude regime, the ratio of the dissipated powers as calculated by the current and proposed model were compared with the linear model predictions. At higher amplitude oscillations, the dissipated power by radiation loss as calculated by the current and the proposed models were compared with the dissipated power calculated using the scattered power by the bubbles. We show that non-physical values are absent in the proposed model. Moreover, predictions of the proposed approach are identical to the predictions of the linear model and the dissipated power estimated using the scattered pressure by the bubble. We show that damping due to thermal effects, liquid viscosity and radiation heavily depend on the excitation pressure and that the linear model estimations are not valid even at pressures as low as 20 kPa.

Shrinking microbubbles with microfluidics: mathematical modelling to control microbubble sizes.

Microbubbles have applications in industry and life-sciences. In medicine, small encapsulated bubbles (<10 µm) are desirable because of their utility in drug/oxygen delivery, sonoporation, and ultrasound diagnostics. While there are various techniques for generating microbubbles, microfluidic methods are distinguished due to their precise control and ease-of-fabrication. Nevertheless, sub-10 µm diameter bubble generation using microfluidics remains challenging, and typically requires expensive equipment and cumbersome setups. Recently, our group reported a microfluidic platform that shrinks microbubbles to sub-10 µm diameters. The microfluidic platform utilizes a simple microbubble-generating flow-focusing geometry, integrated with a vacuum shrinkage system, to achieve microbubble sizes that are desirable in medicine, and pave the way to eventual clinical uptake of microfluidically generated microbubbles. A theoretical framework is now needed to relate the size of the microbubbles produced and the system's input parameters. In this manuscript, we characterize microbubbles made with various lipid concentrations flowing in solutions that have different interfacial tensions, and monitor the changes in bubble size along the microfluidic channel under various vacuum pressures. We use the physics governing the shrinkage mechanism to develop a mathematical model that predicts the resulting bubble sizes and elucidates the dominant parameters controlling bubble sizes. The model shows a good agreement with the experimental data, predicting the resulting microbubble sizes under different experimental input conditions. We anticipate that the model will find utility in enabling users of the microfluidic platform to engineer bubbles of specific sizes.

Non-invasive monitoring of ultrasound-stimulated microbubble radiation enhancement using photoacoustic imaging.

Modulation of the tumour microvasculature has been demonstrated to affect the effectiveness of radiation, stimulating the search for anti-angiogenic and vascular-disrupting treatment modalities. Microbubbles stimulated by ultrasound have recently been demonstrated as a radiation enhancer when used with different cancer models including PC3. Here, photoacoustics imaging technique was used to assess this treatment's effects on haemoglobin levels and oxygen saturation. Correlations between this modality and power doppler assessments of blood flow, and histology measurements of vascular integrity and cell death were also investigated. Xenograft prostate tumours in SCID mice were treated with 0, 2, or 8 Gy radiation combined with microbubbles exposed to 500 kHz ultrasound at a peak negative pressure of 0, 570, and 750 kPa. Tumours were assessed and levels of total haemoglobin, oxygen saturation were measured using photoacoustics before and 24 hours after treatment along with power doppler measured blood flow. Mice were then sacrificed and tumours were assessed for cell death and vascular composition using immunohistochemistry. Treatments using 8 Gy and microbubbles resulted in oxygen saturation decreasing by 28 ± 10% at 570 kPa and 25 ± 29% at 750 kPa, which corresponded to 44 ± 9% and 40 ± 14% respective decreases in blood flow as measured with power doppler. Corresponding histology indicated 31 ± 5% at 570 kPa and 37 ± 5% at 750 kPa in terms of cell death. There were drops in intact vasculature of 15 ± 2% and 20 ± 2%, for treatments at 570 kPa and 750 kPa. In summary, photoacoustic measures of total haemoglobin and oxygen saturation paralleled changes in power doppler indicators of blood flow. Destruction of tumour microvasculature with microbubble-enhanced radiation also led to decreases in blood flow and was associated with increases in cell death and decreases in intact vasculature as detected with CD31 labeling.

Monitoring structural changes in cells with high-frequency ultrasound signal statistics.

We investigate the use of signal envelope statistics to monitor and quantify structural changes during cell death using an in vitro cell model. Using a f/2.35 transducer (center frequency 20 MHz), ultrasound backscatter data were obtained from pellets of acute myeloid leukemia cells treated with a DNA-intercolating chemotherapy drug, as well as from pellets formed with mixtures of treated and untreated cells. Simulations of signals from pellets of mixtures of cells were generated as a summation of point scatterers. The signal envelope statistics were examined by fitting the Rayleigh and generalized gamma distributions. The fit parameters of the generalized gamma distribution showed sensitivity to structural changes in the cells. The scale parameter showed a 200% increase (p<0.05) between untreated and cells treated for 24 h. The shape parameter showed a 50% increase (p<0.05) over 24 h. Experimental results showed reasonable agreement with simulations. The results indicate that high-frequency ultrasound signal statistics can be used to monitor structural changes within a very low percentage of treated cells in a population, raising the possibility of using this technique in vivo.

High-frequency ultrasound scattering from microspheres and single cells.

Assessing the proportion of biological cells in a volume of interest undergoing structural changes, such as cell death, using high-frequency ultrasound (20-100 MHz), requires the development of a theoretical model of scattering by any arbitrary cell ensemble. A prerequisite to building such a model is to know the scattering by a single cell in different states. In this paper, a simple model for the high-frequency acoustic scattering by one cell is proposed. A method for deducing the backscatter transfer function from a single, subresolution scatterer is also devised. Using this method, experimental measurements of backscatter from homogeneous, subresolution polystyrene microspheres and single, viable eukaryotic cells, acquired across a broad, continuous range of frequencies were compared with elastic scattering theory and the proposed cell scattering model, respectively. The resonant features observed in the backscatter transfer function of microspheres were found to correspond accurately to theoretical predictions. Using the spacing of the major spectral peaks in the transfer functions obtained experimentally, it is possible to predict microsphere diameters with less than 4% error. Such good agreement was not seen between the cell model and the measured backscatter from cells. Possible reasons for this discrepancy are discussed.

Attenuation mapping for monitoring thermal therapy using ultrasound transmission imaging.

The use of an ultrasound (US) transmission imaging system to monitor attenuation changes during tissue heating was investigated. This work presents preliminary results of images obtained from an acoustic camera before, during and after heating tissue phantoms using a heated needle. Two types of tissue-mimicking phantoms were used, agar and polyacrylamide-based. Regions of interests were chosen in images obtained from the real-time imaging system, and the pixel intensity values before, during and after heating were compared. In both phantoms, a decrease in image intensities was observed during heating, indicating an increase in tissue attenuation. Additionally, an irreversible change in image intensity was observed in regions close to the heat source. The reversibility of the intensity change was shown to be a function of the distance from the heating needle to the selected region. Initial results indicate that US transmission imaging can be used to monitor thermal therapy.

Ultrasonic spectral parameter characterization of apoptosis.

Ultrasound (US) spectral analysis methods are used to analyze the radiofrequency (RF) data collected from cell pellets exposed to chemotherapeutics that induce apoptosis and other chemicals that induce nuclear transformations. Calibrated backscatter spectra from regions-of-interest (ROI) were analyzed using linear regression techniques to calculate the spectral slope and midband fit. Two f/2 transducers, with operating frequencies of 30 and 34 MHz (relative bandwidths of 93% and 78%, respectively) were used with a custom-made imaging system that enabled the collection of the raw RF data. For apoptotic cells, the spectral slope increased from 0.37 dB/MHz before drug exposure to 0.57 dB/MHz 24 h after, corresponding to a change in effective scatterer radius from 8.7 to 3.2 microm. The midband fit increased in a time-dependent fashion, peaking at 13dB 24 h after exposure. The statistical deviation of the spectral parameters was in close agreement with theoretical predictions. The results provide a framework for using spectral parameter methods to monitor apoptosis in in vitro and in in vivo systems and are being used to guide the design of system and signal analysis parameters.

Comparison of thermal damage calculated using magnetic resonance thermometry, with magnetic resonance imaging post-treatment and histology, after interstitial microwave thermal therapy of rabbit brain.

Clinical application of high-temperature thermal therapy as a treatment for solid tumours requires an accurate and close to real-time method for assessing tissue damage. Imaging methods that detect structural changes during heating may underestimate the extent of thermal damage. This is due to the occurrence of delayed damage manifested at tissue locations exposed to temperatures lower than those required to cause immediate structural changes. An alternative approach is to measure temperature and then calculate the expected damage based on the temperature history at each tissue location. Magnetic resonance (MR) imaging methods now allow temperature maps of the target and surrounding tissues to be generated in almost real-time. The aim of this work was to evaluate whether thermal damage zones calculated on the basis of MR thermometry maps measured during heating correspond to actual tissue damage as measured after treatment by histological methods and MR imaging. Four male rabbits were treated with high-temperature thermal therapy delivered in the brain by a single microwave antenna operating at 915 MHz. MR scanning was performed before, during and after treatment in a 1.5 T whole-body scanner. Temperature maps were produced using the proton resonance frequency (PRF) shift method of MR thermometry. In addition, conventional T1-weighted and T2-weighted spin-echo images were acquired after treatment. Thermal damage zones corresponding to cell death, microvascular blood flow stasis and protein coagulation were calculated using an Arrhenius analysis of the MR temperature/time course data. The calculated zones were compared with the lesions seen on histopathological examination of the brains which were removed within 6-8 h of treatment. The results showed that calculated damage zones based on MR thermometry agreed well with areas of damage as assessed using histology after heating was completed. The data suggest that real-time calculations of final expected thermal damage based on an Arrhenius analysis of MR temperature data may provide a useful method of real-time monitoring of thermal therapy when combined with conventional T2-weighted images taken after treatment.

The effects of dynamic optical properties during interstitial laser photocoagulation.

A nonlinear mathematical model was developed and experimentally validated to investigate the effects of changes in optical properties during interstitial laser photocoagulation (ILP). The effects of dynamic optical properties were calculated using the Arrhenius damage model, resulting in a nonlinear optothermal response. This response was experimentally validated by measuring the temperature rise in albumen and polyacrylamide phantoms. A theoretical study of ILP in liver was conducted constraining the peak temperatures below the vaporization threshold. The temperature predictions varied considerably between the static and dynamic scenarios, and were confirmed experimentally in phantoms. This suggests that the Arrhenius model can be used to predict dynamic changes in optical and thermal fields. An increase in temperature rise due to a decrease in light penetration within the coagulated region during ILP of the liver was also demonstrated. The kinetics of ILP are complex and nonlinear due to coagulation, which changes the tissue properties during treatment. These complex effects can be adequately modelled using an Arrhenius damage formulation.

Ultrasound imaging of apoptosis: high-resolution non-invasive monitoring of programmed cell death in vitro, in situ and in vivo.

A new non-invasive method for monitoring apoptosis has been developed using high frequency (40 MHz) ultrasound imaging. Conventional ultrasound backscatter imaging techniques were used to observe apoptosis occurring in response to anticancer agents in cells in vitro, in tissues ex vivo and in live animals. The mechanism behind this ultrasonic detection was identified experimentally to be the subcellular nuclear changes, condensation followed by fragmentation, that cells undergo during apoptosis. These changes dramatically increase the high frequency ultrasound scattering efficiency of apoptotic cells over normal cells (25- to 50-fold change in intensity). The result is that areas of tissue undergoing apoptosis become much brighter in comparison to surrounding viable tissues. The results provide a framework for the possibility of using high frequency ultrasound imaging in the future to non-invasively monitor the effects of chemotherapeutic agents and other anticancer treatments in experimental animal systems and in patients.

An investigation of the flow dependence of temperature gradients near large vessels during steady state and transient tissue heating.

Temperature distributions measured during thermal therapy are a major prognostic factor of the efficacy and success of the procedure. Thermal models are used to predict the temperature elevation of tissues during heating. Theoretical work has shown that blood flow through large blood vessels plays an important role in determining temperature profiles of heated tissues. In this paper, an experimental investigation of the effects of large vessels on the temperature distribution of heated tissue is performed. The blood flow dependence of steady state and transient temperature profiles created by a cylindrical conductive heat source and an ultrasound transducer were examined using a fixed porcine kidney as a flow model. In the transient experiments, a 20 s pulse of hot water, 30 degrees C above ambient, heated the tissues. Temperatures were measured at selected locations in steps of 0.1 mm. It was observed that vessels could either heat or cool tissues depending on the orientation of the vascular geometry with respect to the heat source and that these effects are a function of flow rate through the vessels. Temperature gradients of 6 degrees C mm(-1) close to large vessels were routinely measured. Furthermore, it was observed that the temperature gradients caused by large vessels depended on whether the heating source was highly localized (i.e. a hot needle) or more distributed (i.e. external ultrasound). The gradients measured near large vessels during localized heating were between two and three times greater than the gradients measured during ultrasound heating at the same location, for comparable flows. Moreover, these gradients were more sensitive to flow variations for the localized needle heating. X-ray computed tomography data of the kidney vasculature were in good spatial agreement with the locations of all of the temperature variations measured. The three dimensional vessel path observed could account for the complex features of the temperature profiles. The flow dependences of the transient temperature profiles near large vessels during the pulsed experiments were consistent with the temperature distributions measured in the steady state experiments and provided unique insights into the process of convective heat transfer in tissues. Finally, it was shown that even for very short treatment times (3-20 s), large vessels had significant effects on the tissue temperature distributions.

Experimental evaluation of two simple thermal models using transient temperature analysis.

Thermal models are used to predict temperature distributions of heated tissues during thermal therapies. Recent interest in short duration high temperature therapeutic procedures necessitates the accurate modelling of transient temperature profiles in heated tissues. Blood flow plays an important role in tissue heat transfer and the resultant temperature distribution. This work examines the transient predictions of two simple mathematical models of heat transfer by blood flow (the bioheat transfer equation model and the effective thermal conductivity equation model) and compares their predictions to measured transient temperature data. Large differences between the two models are predicted in the tissue temperature distribution as a function of blood flow for a short heat pulse. In the experiments a hot water needle, approximately 30 degrees C above ambient, delivered a 20 s heating pulse to an excised fixed porcine kidney that was used as a flow model. Temperature profiles of a thermocouple that primarily traversed the kidney cortex were examined. Kidney locations with large vessels were avoided in the temperature profile analysis by examination of the vessel geometry using high resolution computed tomography angiography and the detection of the characteristic large vessel localized cooling or heating patterns in steady-state temperature profiles. It was found that for regions without large vessels, predictions of the Pennes bioheat transfer equation were in much better agreement with the experimental data when compared to predictions of the scalar effective thermal conductivity equation model. For example, at a location r approximately 2 mm away from the source, the measured delay time was 10.6 +/- 0.5 s compared to predictions of 9.4 s and 5.4 s of the BHTE and ETCE models, respectively. However, for the majority of measured locations, localized cooling and heating effects were detected close to large vessels when the kidney was perfused. Finally, it is shown that increasing flow in regions without large vessels minimally perturbs temperature profiles for short exposure times; regions with large vessels still have a significant effect.

A theoretical comparison of energy sources--microwave, ultrasound and laser--for interstitial thermal therapy.

A number of heating sources are available for minimally invasive thermal therapy of tumours. The purpose of this work was to compare, theoretically, the heating characteristics of interstitial microwave, laser and ultrasound sources in three tissue sites: breast, brain and liver. Using a numerical method, the heating patterns, temperature profiles and expected volumes of thermal damage were calculated during standard treatment times with the condition that tissue temperatures were not permitted to rise above 100 degrees C (to ensure tissue vaporization did not occur). Ideal spherical and cylindrical applicators (200 microm and 800 microm radii respectively) were modelled for each energy source to demonstrate the relative importance of geometry and energy attenuation in determining heating and thermal damage profiles. The theoretical model included the effects of the collapse of perfusion due to heating. Heating patterns were less dependent on the energy source when small spherical applicators were modelled than for larger cylindrical applicators due to the very rapid geometrical decrease in energy with distance for the spherical applicators. For larger cylindrical applicators, the energy source was of greater importance. In this case, the energy source with the lowest attenuation coefficient was predicted to produce the largest volume of thermally coagulated tissue, in each tissue site.

Magnetic resonance imaging of temperature changes during interstitial microwave heating: a phantom study.

Changes in magnetic resonance (MR) signals during interstitial microwave heating are reported, and correlated with simultaneously acquired temperature readings from three fiber-optic probes implanted in a polyacrylamide gel phantom. The heating by a MR-compatible microwave antenna did not interfere with simultaneous MR image data acquisition. MR phase-difference images were obtained using a fast two-dimensional-gradient echo sequence. From these images the temperature-sensitive resonant frequency of the 1H nuclei was found to decrease approximately by 0.008 ppm/ degree C. The method and results presented here demonstrate that noninvasive MR-temperature imaging can be performed simultaneously with interstitial microwave thermal treatment.

Ultrasonic biomicroscopy of viable, dead and apoptotic cells.

Ultrasonic imaging is frequently used in medical diagnosis to differentiate normal and tumour tissues. Here we investigate if distinct types of cell death can be discriminated through the use of ultrasound biomicroscopy. By using a well-controlled system in vitro, we demonstrate that this imaging modality can be used to differentiate living cells, dead cells and cells that have died by programmed cell death or apoptosis. The results indicate a greater than twofold ultrasound backscatter signal from apoptotic cells in comparison to viable cells, whereas heat-killed cells exhibit an intermediate level of ultrasound backscatter. The results have potential implications in the study of disease-related biological processes involving apoptosis.