High frequency ultrasonic backscatter from erythrocyte suspension

Specification and Evaluation of Plasticizer Migration Simulants for Human Blood Products: A Delphi Study

Potentially toxic plasticizers are commonly added to polyvinyl chloride medical devices for transfusion in order to improve their flexibility and workability. As the plasticizers are not chemically bonded to the PVC, they can be released into labile blood products (LBPs) during storage. Ideally, LBPs would be used in laboratory studies of plasticizer migration from the medical device. However, short supply (i.e., limited stocks of human blood in collection centres) has prompted the development of specific simulants for each type of LBP in the evaluation of new transfusion devices. We performed a Delphi study with a multidisciplinary panel of 24 experts. In the first (qualitative) phase, the panel developed consensus definitions of the specification criteria to be met by each migration simulant. Next, we reviewed the literature on techniques for simulating the migration of plasticizers into LBPs. A questionnaire was elaborated and sent out to the experts, and the replies were synthesized in order to obtain a consensus. The qualitative study established specifications for each biological matrix (whole blood, red blood cell concentrate, plasma, and platelet concentrate) and defined the criteria required for a suitable LBP simulant. Ten criteria were suggested: physical and chemical characteristics, opacity, form, stability, composition, ability to mimic a particular clinical situation, ease and safety of use, a simulant–plastic interaction correlated with blood, and compatibility with analytical methods. The questionnaire data revealed a consensus on the use of natural products (such as pig’s blood) to mimic the four LBPs. Opinions diverged with regard to synthetic products. However, an isotonic solution and a rheological property modifier were considered to be of value in the design of synthetic simulants. Consensus reached by the Delphi group could be used as a database for the development of simulants used to assess the migration of plasticizers from PVC bags into LBPs.

High frequency ultrasound imaging and simulations of sea urchin oocytes

High frequency ultrasound backscatter signals from sea urchin oocytes were measured using a 40 MHz transducer and compared to numerical simulations. The Faran scattering model was used to calculate the ultrasound scattered from single oocytes in suspension. The urchin oocytes are non-nucleated with uniform size and biomechanical properties; the backscatter from each cell is similar and easy to simulate, unlike typical nucleated mammalian cells. The time domain signal measured from single oocytes in suspension showed two distinct peaks, and the power spectrum was periodic with minima spaced approximately 10 MHz apart. Good agreement to the Faran scattering model was observed. Measurements from tightly packed oocyte cell pellets showed similar periodic features in the power spectra, which was a result of the uniform size and consistent biomechanical properties of the cells. Numerical simulations that calculated the ultrasound scattered from individual oocytes within a three dimensional volume showed good agreement to the measured signals and B-scan images. A cepstral analysis of the signal was used to calculate the size of the cells, which was 78.7 μm (measured) and 81.4 μm (simulated). This work supports the single scattering approximation, where ultrasound is discretely scattered from single cells within a bulk homogeneous sample, and that multiple scattering has a negligible effect. This technique can be applied towards understanding the complex scattering behaviour from heterogeneous tissues.

Quantification of Very Low Concentrations of Leukocyte Suspensions In Vitro by High-Frequency Ultrasound

Accurate measurement of very low cerebrospinal fluid (CSF) white blood cell (WBC) concentration is key to the diagnosis of bacterial meningitis, lethal if not promptly treated. Here we show that high frequency ultrasound (HFUS) can detect CSF WBC in vitro in concentrations relevant to meningitis diagnosis with a much finer precision than gold standard manual counting in a Fuchs-Rosenthal chamber. WBC concentrations in a mock CSF model, in the range 0-50 WBC/μL, have been tested and compared to gold standard ground truth. In this range, excellent agreement (Cohen's kappa [κ] = 0.78-90) (Cohen 1960) was observed between HFUS and the gold standard method. The presented experimental set-up allowed us to detect WBC concentrations as low as 2 cells/μL. HFUS shows promise as a low-cost, reliable and automated technology to measure very low CSF WBC concentrations for the diagnosis of early meningitis.

Ultrasonic scattering measurements of a live single cell at 86 MHz

Cell separation and sorting techniques have been employed biomedical applications such as cancer diagnosis and cell gene expression analysis. The capability to accurately measure ultrasonic scattering properties from cells is crucial in making an ultrasonic cell sorter a reality if ultrasound scattering is to be used as the sensing mechanism as well. To assess the performance of sensing and identifying live single cells with high-frequency ultrasound, an 86-MHz lithium niobate press-focused single-element acoustic transducer was used in a high-frequency ultrasound scattering measurement system that was custom designed and developed for minimizing noise and allowing better mobility. Peak-to-peak echo amplitude, integrated backscatter (IB) coefficient, spectral parameters including spectral slope and intercept, and midband fit from spectral analysis of the backscattered echoes were measured and calculated from a live single cell of two different types on an agar surface: leukemia cells (K562 cells) and red blood cells (RBCs). The amplitudes of echo signals from K562 cells and RBCs were 48.25 ± 11.98 mV(pp) and 56.97 ± 7.53 mV(pp), respectively. The IB coefficient was -89.39 ± 2.44 dB for K562 cells and -89.00 ± 1.19 dB for RBCs. The spectral slope and intercept were 0.30 ± 0.19 dB/MHz and -56.07 ± 17.17 dB, respectively, for K562 cells and 0.78 ± 0.092 dB/MHz and -98.18 ± 8.80 dB, respectively, for RBCs. Midband fits of K562 cells and RBCs were -31.02 ± 3.04 dB and -33.51 ± 1.55 dB, respectively. Acoustic cellular discrimination via these parameters was tested by Student's t-test. Their values, except for the IB value, showed statistically significant difference (p < 0.001). This paper reports for the first time that ultrasonic scattering measurements can be made on a live single cell with a highly focused high-frequency ultrasound microbeam at 86 MHz. These results also suggest the feasibility of ultrasonic scattering as a sensing mechanism in the development of ultrasonic cell sorters.

High-resolution tissue Doppler imaging of the zebrafish heart during its regeneration

The human heart cannot regenerate after injury, whereas the adult zebrafish can fully regenerate its heart even after 20% of the ventricle is amputated. Many studies have begun to reveal the cellular and molecular mechanisms underlying this regenerative process, which have exciting implications for human cardiac diseases. However, the dynamic functions of the zebrafish heart during regeneration are not yet understood. This study established a high-resolution echocardiography for tissue Doppler imaging (TDI) of the zebrafish heart to explore the cardiac functions during different regeneration phases. Experiments were performed on AB-line adult zebrafish (n=40) in which 15% of the ventricle was surgically removed. An 80-MHz ultrasound TDI based on color M-mode imaging technology was employed. The cardiac flow velocities and patterns from both the ventricular chamber and myocardium were measured at different regeneration phases relative to the day of amputation. The peak velocities of early diastolic inflow, early diastolic myocardial motion, late diastolic myocardial motion, early diastolic deceleration slope, and heart rate were increased at 3 days after the myocardium amputation, but these parameters gradually returned to close to their baseline values for the normal heart at 7 days after amputation. The peak velocities of late diastolic inflow, ventricular systolic outflow, and systolic myocardial motion did not significantly differ during the heart regeneration.

Backscattering measurement from a single microdroplet

Backscattering measurements for acoustically trapped lipid droplets were undertaken by employing a P[VDF-TrFE] broadband transducer of f-number = 1, with a bandwidth of 112%. The wide bandwidth allowed the transmission of the 45 MHz trapping signal and the 15 MHz sensing signal using the same transducer. Tone bursts at 45 MHz were first transmitted by the transducer to hold a single droplet at the focus (or the center of the trap) and separate it from its neighboring droplets by translating the transducer perpendicularly to the beam axis. Subsequently, 15 MHz probing pulses were sent to the trapped droplet and the backscattered RF echo signal received by the same transducer. The measured beam width at 15 MHz was measured to be 120 μ m. The integrated backscatter (IB) coefficient of an individual droplet was determined within the 6-dB bandwidth of the transmit pulse by normalizing the power spectrum of the RF signal to the reference spectrum obtained from a flat reflector. The mean IB coefficient for droplets with a 64 μ m average diameter (denoted as cluster A) was -107 dB, whereas it was -93 dB for 90-μm droplets (cluster B). The standard deviation was 0.9 dB for each cluster. The experimental values were then compared with those computed with the T-matrix method and a good agreement was found: the difference was as small as 1 dB for both clusters. These results suggest that this approach might be useful as a means for measuring ultrasonic backscattering from a single microparticle, and illustrate the potential of acoustic sensing for cell sorting.

High-frequency attenuation and backscatter measurements of rat blood between 30 and 60 MHz

There has recently been a great deal of interest in noninvasive high-frequency ultrasound imaging of small animals such as rats due to their being the preferred animal model for gene therapy and cancer research. Improving the interpretation of the obtained images and furthering the development of the imaging devices require a detailed knowledge of the ultrasound attenuation and backscattering of biological tissue (e.g. blood) at high frequencies. In the present study, the attenuation and backscattering coefficients of the rat red blood cell (RBC) suspensions and whole blood with hematocrits ranging from 6% to 40% were measured between 30 and 60 MHz using a modified substitution approach. The acoustic parameters of porcine blood under the same conditions were also measured in order to compare differences in the blood properties between these two animals. For porcine blood, both whole blood and RBC suspension were stirred at a rotation speed of 200 rpm. Three different rotation speeds of 100, 200 and 300 rpm were carried out for rat blood experiments. The attenuation coefficients of both rat and porcine blood were found to increase linearly with frequency and hematocrit (the values of coefficients of determination (r(2)) are around 0.82-0.97 for all cases). The average attenuation coefficient of rat whole blood with a hematocrit of 40% increased from 0.26 Nepers mm(-1) at 30 MHz to 0.47 Nepers mm(-1) at 60 MHz. The maximum backscattering coefficients of both rat and porcine RBC suspensions were between 10% and 15% hematocrits at all frequencies. The fourth-power dependence of backscatter on frequency was approximately valid for rat RBC suspensions with hematocrits between 6% and 40%. However, the frequency dependence of the backscatter estimate deviates from a fourth-power law for porcine RBC suspension with hematocrit higher than 20%. The backscattering coefficient plateaued for hematocrits higher than 15% in porcine blood, but for rat blood it was maximal around a hematocrit of 20% at the same rotation speed, and shifted to a hematocrit of 10% at a higher speed. The backscattering properties of rat RBCs in plasma are similar to those of RBCs in saline at a higher rotation speed. The differences in attenuation and backscattering between rat and porcine blood may be attributed to RBCs' being smaller and the RBC aggregation level being lower for rat blood than for porcine blood.

High-frequency ultrasound backscattering by blood: analytical and semianalytical models of the erythrocyte cross section

This paper proposes analytical and semianalytical models of the ultrasonic backscattering cross section (BCS) of various geometrical shapes mimicking a red blood cell (RBC) for frequencies varying from 0 to 90 MHz. By assuming the first-order Born approximation and by modeling the shape of a RBC by a realistic biconcave volume, different scattering behaviors were identified for increasing frequencies. For frequencies below 18 MHz, a RBC can be considered a Rayleigh scatterer. For frequencies less than 39 MHz, the general concept of acoustic inertia tensor is introduced to describe the variation of the BCS with the frequency and the incidence direction. For frequencies below 90 MHz, ultrasound backscattering by a RBC is equivalent to backscattering by a cylinder of height 2 microm and diameter 7.8 microm. These results lay the basis of ultrasonic characterization of RBC aggregation by proposing a method that distinguishes the contribution of the individual RBC acoustical characteristics from collective effects, on the global blood backscattering coefficient. A new method of data reduction that models the frequency dependence of the ultrasonic BCS of micron-sized weak scatterers is also proposed. Applications of this method are in tissue characterization as well as in hematology.

Echogenicity variations from porcine blood II: the "bright ring" under oscillatory flow

Echogenicity variations from porcine blood were observed in a mock flow loop under pulsatile flow in a series of experiments (Paeng et al. 2004). In this paper, oscillatory flow was generated to further investigate the cyclic and radial variation of blood echogenicity and its origin and mechanisms by several parameters, including stroke volume, stroke rate, mean steady flow and transducer angle, using a GE LOGIQ 700 Expert system. The echogenicity at the center of the tube was enhanced during acceleration and lower during deceleration, and the expansion and collapse of the "bright ring" was observed twice per cycle. The "black hole," a central echo-poor zone surrounded by a hyperechoic zone, was barely observable under oscillatory flow, and these patterns differed from those under pulsatile flow. The cyclic and radial variation of echogenicity under oscillatory flow was affected by such hemodynamic parameters as stroke volume, stroke rate and mean steady flow. It was suggested that rouleaux might be aligned at an angle of about 25 degrees relative to the tube axis during the acceleration phase, based on the experimental results reaching a maximum of the echogenicity variation at a transducer angle of 25 degrees. Radial distribution of rouleaux alignments was proposed to be another important factor to blood echogenicity variation, in addition to combined effects of shear rate and flow acceleration on erythrocyte aggregation and blood echogenicity. The weak cyclic variation of echogenicity was also observed from the porcine erythrocyte suspensions under pure oscillatory flow, but not under pulsatile flow. It is postulated that the echogenicity variations from erythrocyte suspensions are from red cell deformation.

Cyclic and radial variation of the Doppler power from porcine whole blood

The Doppler power from porcine blood was observed in a mock flow loop to have cyclic and radial variation during a pulsatile cycle. It was found to decrease with shear rate under steady flow, except near the center of the tube at which other mechanisms such as the effects of radial distribution on the rouleaux might be involved. Under pulsatile flow, the timing of the peak of the Doppler power measured at the center of the tube became closer to the peak systole from 20 to 60 beats/minute (BPM), and the power and velocity peaks coincided at 60 BPM. The overall radial variation of the Doppler power during a whole pulsatile cycle was prominent due to the increase of shear rate from the center to 4.5 mm radial position within a tube of 6.35 mm radius. The cyclic variation of the Doppler power varied with the radial position, being relatively large at the center, reaching a minimum at an intermediate radial position, and increasing again near the wall. The peak of the Doppler power occurred at early systole near the tube wall and lagged the flow closer to the center. The "black hole" phenomenon was observed only over portions of the flow cycle. All these complex variations of the Doppler power across the tube over a cycle are thought to be the result of red cell aggregation, which can be affected by shear rate and acceleration.

Modeling the frequency dependence (5-120 MHz) of ultrasound backscattering by red cell aggregates in shear flow at a normal hematocrit

The frequency dependence of the ultrasound signal backscattered by blood in shear flow was studied using a simulation model. The ultrasound backscattered signal was computed with a linear model that considers the characteristics of the ultrasound system and tissue acoustic properties. The tissue scattering properties were related to the position and shape of the red blood cells (RBCs). A 2D microrheological model simulated the RBC dynamics in a Couette shear flow system. This iterative model, described earlier [Biophys. J. 82, 1696-1710 (2002)], integrates the hydrodynamic effect of the flow, as well as adhesive and repulsive forces between RBCs. RBC aggregation was simulated at 40% hematocrit and shear rates of 0.05-2 s(-1). The RBC aggregate sizes ranged, on average, from 3.3 RBCs at 2 s(-1) to 33.5 cells at 0.05 s(-1). The ultrasound backscattered power was studied at frequencies between 5-120 MHz and insonification angles between 0-180 degrees. At frequencies below approximately 30 MHz, the ultrasound backscattered power increased as the shear rate was decreased and the size of the aggregates was raised. A totally different scattering behavior was noted above 30 MHz. Typical spectral slopes of the backscattered power (log-log scale) between 5-25 MHz equaled 3.8, whereas slopes down to 0.6 were measured at 0.05 s(-1), between 40-60 MHz. The ultrasound backscattered power was shown to be angle dependent at low frequencies (5-25 MHz). The anisotropy persisted at high frequencies (>25 MHz) for small aggregates (at 2 s(-1)). In conclusion, this study sheds some light on the blood backscattering behavior with an emphasis on the non-Rayleigh regime. Additional experimental studies may be necessary to validate the simulation results, and to fully understand the relation between the ultrasound backscattered power, level of RBC aggregation, shear rate, frequency, and insonification angle.

A new high resolution color flow system using an eigendecomposition-based adaptive filter for clutter rejection

We present a new signal processing strategy for high frequency color flow mapping in moving tissue environments. A new application of an eigendecomposition-based clutter rejection filter is presented with modifications to deal with high blood-to-clutter ratios (BCR). Additionally, a new method for correcting blood velocity estimates with an estimated tissue motion profile is detailed. The performance of the clutter filter and velocity estimation strategies is quantified using a new swept-scan signal model. In vivo color flow images are presented to illustrate the potential of the system for mapping blood flow in the microcirculation with external tissue motion.

A new high resolution color flow system using an eigendecomposition-based adaptive filter for clutter rejection

We present a new signal processing strategy for high frequency color flow mapping in moving tissue environments. A new application of an eigendecomposition-based clutter rejection filter is presented with modifications to deal with high blood-to-clutter ratios (BCR). Additionally, a new method for correcting blood velocity estimates with an estimated tissue motion profile is detailed. The performance of the clutter filter and velocity estimation strategies is quantified using a new swept-scan signal model. In vivo color flow images are presented to illustrate the potential of the system for mapping blood flow in the microcirculation with external tissue motion.

Spatial velocity profile in mouse embryonic aorta and Doppler-derived volumetric flow: a preliminary model

Characterizing embryonic circulatory physiology requires accurate cardiac output and flow data. Despite recent applications of high-frequency ultrasound Doppler to the study of embryonic circulation, current Doppler analysis of volumetric flow is relatively crude. To improve Doppler derivation of volumetric flow, we sought a preliminary model of the spatial velocity profile in the mouse embryonic dorsal aorta using ultrasound biomicroscopy (UBM)-Doppler data. Embryonic hematocrit is 0.05-0.10 so rheologic properties must be insignificant. Low Reynolds numbers (<500) and Womersley parameters (<0.76) suggest laminar flow. UBM demonstrated a circular dorsal aortic cross section with no significant tapering. Low Dean numbers (<100) suggest the presence of minimal skewing of the spatial velocity profile. The inlet length allows for fully developed flow. There is no apparent aortic wall pulsatility. Extrapolation of prior studies to these vessel diameters (300-350 microm) and flow velocities (~50-200 mm/s) suggests parabolic spatial velocity profiles. Therefore, mouse embryonic dorsal aortic blood flow may correspond to Poiseuille flow in a straight rigid tube with parabolic spatial velocity profiles. As a first approximation, these results are an important step toward precise in utero ultrasound characterization of blood flow within the developing mammalian circulation.

The significance of shape and orientation in single-particle weak-scatterer models

Particles that have a density and compressibility comparable to those of the surrounding medium (weak scatterers) are often modeled as spherical scatterers of equivalent volume. The simplicity and symmetry of the spherical model does not however account for the effect of the angle of incidence of the incident field onto the particle. By performing a comparative study of scattering models for the case of a fluid sphere, it is shown that models derived using the Born approximation closely match the exact solution for weak scatterers. This approximation is therefore applied to produce a model for the fluid disk, which is used to investigate the significance of particle shape and orientation. Scattering by a red blood cell of given volume modeled as a sphere is compared to the result obtained by approximating its shape as a disk of varying aspect ratio and orientation. The spherical model is shown to provide a good description for frequencies up to 20 MHz, beyond which particle shape becomes significant. This effect could go undetected if the scattered field is only observed at 90 degrees relative to the direction of the incident field. Nevertheless, the significance of particle shape and orientation might form the basis of a novel detection technique.

High-frequency backscatter and attenuation measurements of porcine erythrocyte suspensions between 30-90 MHz

There are now diagnostic ultrasonic imaging devices that operate at very high frequencies (VHF) of 20 MHz and beyond for clinical applications in ophthalmology, dermatology and vascular surgery. To be able to better interpret these images and to further the development of these devices, knowledge of ultrasonic attenuation and scattering of biologic tissues, such as blood, in the high-frequency range is crucial. VHF attenuation and backscatter experiments were made on porcine red blood cell (RBC) suspensions, for which much data on attenuation and backscatter can be found in the literature in the lower frequency range. Attenuation and backscatter at hematocrits of 6%, 10%, 15%, 20%, 25% and 30% from 30 to 90 MHz were measured using a modified substitution method that allows the utilization of focused transducers. The results show that the attenuation coefficient from all suspensions increased linearly with frequency and the backscatter coefficient for low hematocrit suspensions was found to have a maximum between 10% and 15%. At higher hematocrits, a decrease in the frequency-dependence was observed, possibly indicating that Rayleigh scattering is no longer valid because the wavelength in the VHF range is comparable to the size of a porcine RBC.

A system-based approach to modeling the ultrasound signal backscattered by red blood cells

A system-based model is proposed to describe and simulate the ultrasound signal backscattered by red blood cells (RBCs). The model is that of a space-invariant linear system that takes into consideration important biological tissue stochastic scattering properties as well as the characteristics of the ultrasound system. The formation of the ultrasound signal is described by a convolution integral involving a transducer transfer function, a scatterer prototype function, and a function representing the spatial arrangement of the scatterers. The RBCs are modeled as nonaggregating spherical scatterers, and the spatial distribution of the RBCs is determined using the Percus-Yevick packing factor. Computer simulations of the model are used to study the power backscattered by RBCs as a function of the hematocrit, the volume of the scatterers, and the frequency of the incident wave (2-500 MHz). Good agreement is obtained between the simulations and theoretical and experimental data for both Rayleigh and non-Rayleigh scattering conditions. In addition to these results, the renewal process theory is proposed to model the spatial arrangement of the scatterers. The study demonstrates that the system-based model is capable of accurately predicting important characteristics of the ultrasound signal backscattered by blood. The model is simple and flexible, and it appears to be superior to previous one- and two-dimensional simulation studies.

Ultrasonic backscattering from porcine whole blood of varying hematocrit and shear rate under pulsatile flow

It was shown previously that ultrasonic scattering from whole blood varies during a flow cycle under pulsatile flow both in vitro and in vivo. It has been postulated that this cyclic variation may be associated with the dynamics of red cell aggregation because the shearing force acting on the red cell aggregates across the lumen is a function of time during a flow cycle. In all studies, the local shear rate variation as a function of time is unknown. The effect of shear rate on the red cell aggregation and, thus, on ultrasonic scattering from blood can only be merely speculated. One solution to this problem is to estimate the shear rate in a flow conduit by finite element analysis (FEA). An FEA computational fluid dynamics (CFD) tool was used to calculate local shear rate in a series of experiments in which ultrasonic backscattering from porcine whole blood under pulsatile flow was measured as a function of hematocrit and shear rate intravascularly with a 10-MHz catheter-mounted transducer in a mock flow loop. The results show that, at 20 beats per min (BPM), the magnitudes of the cyclic variation for hematocrits at 30, 40, and 50% were approximately 4 dB. However, at 60 BPM, the magnitude of cyclic variation was found to be minimal. The results also confirm previous findings that the amplitude and the timing of the peak of ultrasonic backscattering from porcine whole blood under pulsatile flow during a flow cycle are dependent upon the shear rate and hematocrit in a complicated way.

An approach for measuring ultrasonic backscattering from biological tissues with focused transducers

When the standard substitution method is used with a focused transducer to measure the backscattering coefficient from biological tissues including blood, it yields erroneous results. Extending the backscattering measurements to frequencies beyond 15 MHz necessitates the use of focused transducers because of the worsened signal-to-noise ratio--caused by the increased attenuation and the smaller transducer aperture size--in order to make the measurements close to the transducer. An approach which allows the use of focused transducers in backscattering measurements has been developed. It has been used to measure the backscattering coefficient of red cell suspensions of hematocrit ranging from a few percent to 30% in the frequency range from 5 MHz to 30 MHz. The results at hematocrits below 20% agree well with those obtained with the standard substitution method, although they differ as the hematocrit is increased beyond 20%. The experimental results also show that the fourth-power dependence of backscatter on frequency is in general approximately valid for suspended erythrocytes of hematocrit between 6% and 30%.