Monitoring Dual VEGF Inhibition in Human Pancreatic Tumor Xenografts With Dynamic Contrast-Enhanced Ultrasound.

BACKGROUND: Association of drugs acting against different antiangiogenic mechanisms may increase therapeutic effect and reduce resistance. Noninvasive monitoring of changes in the antiangiogenic response of individual tumors could guide selection and administration of drug combinations. Noninvasive detection of early therapeutic response during dual, vertical targeting of the vascular endothelial growth factor pathway was investigated in an ectopic subcutaneous xenograft model for human pancreatic tumor. METHODS: Dynamic contrast-enhanced ultrasound 12 MHz was used to monitor tumor-bearing Naval Medical Research Institute mice beginning 15 days after tumor implantation. Mice received therapy from 15 to 29 days with sorafenib (N = 9), ziv-aflibercept (N = 11), combined antiangiogenic agents (N = 11), and placebo control (N = 14). Sorafenib (BAY 43-9006; Nexavar), a multikinase inhibitor acting on Raf kinase and receptor tyrosine kinases-including vascular endothelial growth factor receptors 2 and 3-was administered daily (60 mg/kg, per os). Ziv-aflibercept (ZALTRAP), a high-affinity ligand trap blocking the activity of vascular endothelial growth factor A, vascular endothelial growth factor B, and placental growth factor was administered twice per week (40 mg/kg, intraperitoneally). RESULTS: Functional evaluation with dynamic contrast-enhanced ultrasound indicated stable tumor vascularization for the control group while revealing significant and sustained reduction after 1 day of therapy in the combined group (P = .007). There was no survival benefit or penalty due to drug combination. The functional progression-free survival assessed with dynamic contrast-enhanced ultrasound was significantly higher for the 3 treated groups; whereas, the progression-free survival based on tumor size did not discriminate therapeutic effect. CONCLUSIONS: Dynamic contrast-enhanced ultrasound, therefore, presents strong potential to monitor microvascular modifications during antiangiogenic therapy, a key role to monitoring antiangiogenic combining therapy to adapt dose range drug.

Effects of Signal Saturation on QUS Parameter Estimates Based on High-Frequency-Ultrasound Signals Acquired From Isolated Cancerous Lymph Nodes.

Choosing an appropriate dynamic range (DR) for acquiring radio frequency (RF) data from a high-frequency-ultrasound (HFU) system is challenging because signals can vary greatly in amplitude as a result of focusing and attenuation effects. In addition, quantitative ultrasound (QUS) results are altered by saturated data. In this paper, the effects of saturation on QUS estimates of effective scatterer diameter (ESD) and effective acoustic concentration (EAC) were quantified using simulated and experimental RF data. Experimental data were acquired from 69 dissected human lymph nodes using a single-element transducer with a 26-MHz center frequency. Artificially saturated signals ( xc) were produced by thresholding the original unsaturated RF echo signals. Saturation severity was expressed using a quantity called saturate-signal-to-noise ratio (SSNR). Results indicated that saturation has little effect on ESD estimates. However, EAC estimates decreased significantly with decreasing SSNR. An EAC correction algorithm exploiting a linear relationship between EAC values over a range of SSNR values and l1 -norm of xc (i.e., the sum of absolute values of the true RF echo signal) is developed. The maximal errors in EAC estimates resulting from saturation were -8.05, -3.59, and -0.93 dB/mm3 with the RF echo signals thresholded to keep 5, 6, and 7-bit from the original 8-bit DR, respectively. The EAC correction algorithm reduced maximal errors to -3.71, -0.89, and -0.26 dB/mm3 when signals were thresholded at 5, 6, and 7-bit, respectively.

Local Transverse-Slice-Based Level-Set Method for Segmentation of 3-D High-Frequency Ultrasonic Backscatter From Dissected Human Lymph Nodes.

OBJECTIVE: To detect metastases in freshly excised human lymph nodes (LNs) using three-dimensional (3-D), high-frequency, quantitative ultrasound (QUS) methods, the LN parenchyma (LNP) must be segmented to preclude QUS analysis of data in regions outside the LNP and to compensate ultrasound attenuation effects due to overlying layers of LNP and residual perinodal fat (PNF). METHODS: After restoring the saturated radio-frequency signals from PNF using an approach based on smoothing cubic splines, the three regions, i.e., LNP, PNF, and normal saline (NS), in the LN envelope data are segmented using a new, automatic, 3-D, three-phase, statistical transverseslice-based level-set (STS-LS) method that amends Lankton's method. Due to ultrasound attenuation and focusing effects, the speckle statistics of the envelope data vary with imaged depth. Thus, to mitigate depth-related inhomogeneity effects, the STS-LS method employs gamma probabilitydensity functions to locally model the speckle statistics within consecutive transverse slices. RESULTS: Accurate results were obtained on simulated data. On a representative dataset of 54 LNs acquired from colorectal-cancer patients, the Dice similarity coefficient for LNP, PNF, and NS were 0.938 ± 0.025, 0.832 ± 0.086, and 0.968 ± 0.008, respectively, when compared to expert manual segmentation. CONCLUSION: The STS-LS outperforms the established methods based on global and local statistics in our datasets and is capable of accurately handling the depth-dependent effects due to attenuation and focusing. SIGNIFICANCE: This advance permits the automatic QUS-based cancer detection in the LNs. Furthermore, the STS-LS method could potentially be used in a wide range of ultrasound-imaging applications suffering from depth-dependent effects.

Automatic motion estimation using flow parameters for dynamic contrast-enhanced ultrasound.

Dynamic contrast-enhanced ultrasound (DCE-US) sequences are subject to motion which can disturb functional flow quantification. This can make estimated parameters more variable or unreliable. Methods that compensate for motion are therefore desirable. The most commonly used motion correction techniques in DCE-US register the images in the sequence with respect to a user-selected reference image. However, this image may not include all features that are representative of the whole sequence. Moreover, image-based registration neglects pertinent, functional-flow information contained in the DCE-US sequence. An operator-free method is proposed that combines the motion estimation and flow-parameter quantification (M/Q method) in a single mathematical framework. This method is based on a realistic multiplicative model of the DCE-US noise. By computing likelihood in this model, motion and flow parameters are both estimated iteratively. First, the maximization is accomplished by estimating functional and motion parameters. Then, a final registration based on a non-parametric temporal smoothing of the sequence is performed. This method is compared to a conventional (mutual information) registration method where all the images of the sequence are registered with respect to a reference image chosen by an expert. The two methods are evaluated on simulated sequences and DCE-US sequences acquired in patients (N = 15). The M/Q method demonstrates significantly (p < 0.05) lower Dice coefficients and Hausdorff distance than the conventional method on the simulated data sets. On the in vivo sequences analysed, the M/Q methods outperformed the conventional method in terms of mean Dice and Hausdorff distance on 80% of the sequences, and in terms of standard deviation of Dice and Hausdorff distance on 87% of the sequences.

Modeling the envelope statistics of three-dimensional high-frequency ultrasound echo signals from dissected human lymph nodes.

This work investigates the statistics of the envelope of three-dimensional (3D) high-frequency ultrasound (HFU) data acquired from dissected human lymph nodes (LNs). Nine distributions were employed, and their parameters were estimated using the method of moments. The Kolmogorov Smirnov (KS) metric was used to quantitatively compare the fit of each candidate distribution to the experimental envelope distribution. The study indicates that the generalized gamma distribution best models the statistics of the envelope data of the three media encountered: LN parenchyma, fat and phosphate-buffered saline (PBS). Furthermore, the envelope statistics of the LN parenchyma satisfy the pre-Rayleigh condition. In terms of high fitting accuracy and computationally efficient parameter estimation, the gamma distribution is the best choice to model the envelope statistics of LN parenchyma, while, the Weibull distribution is the best choice to model the envelope statistics of fat and PBS. These results will contribute to the development of more-accurate and automatic 3D segmentation of LNs for ultrasonic detection of clinically significant LN metastases.

Dual-mode registration of dynamic contrast-enhanced ultrasound combining tissue and contrast sequences.

This study proposes a new method for automatic, iterative image registration in the context of dynamic contrast-enhanced ultrasound (DCE-US) imaging. By constructing a cost function of image registration using a combination of the tissue and contrast-microbubble responses, this new method, referred to as dual-mode registration, performs alignment based on both tissue and vascular structures. Data from five focal liver lesions (FLLs) were used for the evaluation. Automatic registration based on the dual-mode registration technique and tissue-mode registration obtained using the linear response image sequence alone were compared to manual alignment of the sequence by an expert. Comparison of the maximum distance between the transformations applied by the automatic registration techniques and those from expert manual registration reference showed that the dual-mode registration provided better precision than the tissue-mode registration for all cases. The reduction of maximum distance ranged from 0.25 to 9.3mm. Dual-mode registration is also significantly better than tissue-mode registration for the five sequences with p-values lower than 0.03. The improved sequence alignment is also demonstrated visually by comparison of images from the sequences and the video playbacks of the motion-corrected sequences. This new registration technique better maintains a selected region of interest (ROI) within a fixed position of the image plane throughout the DCE-US sequence. This should reduce motion-related variability of the echo-power estimations and, thus, contribute to more robust perfusion quantification with DCE-US.

A multiplicative model for improving microvascular flow estimation in dynamic contrast-enhanced ultrasound (DCE-US): theory and experimental validation.

Perfusion parameter estimation from dynamic contrast-enhanced ultrasound (DCE-US) data relies on fitting parametric models of flow to curves describing linear echo power as a function of time. The least squares criterion is generally used to fit these models to data. This criterion is optimal in the sense of maximum likelihood under the assumption of an additive white Gaussian noise. In the current work, it is demonstrated that this assumption is not held for DCEUS. A better-adapted maximum likelihood criterion based on a multiplicative model is proposed. It is tested on simulated bolus perfusion data and on 11 sequences acquired in vivo during bolus perfusion of contrast agent in the cortex of healthy murine kidney, an area where the perfusion is expected to be approximately homogeneous. Results on simulated data show a significant improvement (p < 0.05) of the precision and the accuracy for the estimations of perfusion parameters time to peak (TTP), wash-in rate (WiR), and mean transit time (MTT). On the 11 in vivo sequences, the new method leads to a significant reduction (p < 0.05) in the variation of parametric maps for 9 sequences for TTP and 10 sequences for WiR and MTT. The mean percent decreases of the coefficient of variation are 40%, 25%, and 59% for TTP, WiR, and MTT, respectively. This method should contribute to a more robust and accurate estimation of perfusion parameters and an improved resolution of parametric imaging.

Echo-power estimation from log-compressed video data in dynamic contrast-enhanced ultrasound imaging.

Ultrasound (US) scanners typically apply lossy, non-linear modifications to the US data for visualization purposes. The resulting images are then stored as compressed video data. Some system manufacturers provide dedicated software for quantification purposes to eliminate such processing distortions, at least partially. This is currently the recommended approach for quantitatively assessing changes in contrast-agent concentration from clinical data. However, the machine-specific access to US data and the limited set of analysis functionalities offered by each dedicated-software package make it difficult to perform comparable analyses with different US systems. The objective of this work was to establish if linearization of compressed video images obtained with an arbitrary US system can provide an alternative to dedicated-software analysis of machine-specific files for the estimation of echo-power. For this purpose, an Aplio 50 system (Toshiba Medical Systems, Tochigi, Japan), coupled with dedicated CHI-Q (Contrast Harmonic Imaging Quantification) software by Toshiba Medical Systems, was used. Results were compared with two approaches that apply algorithms to estimate relative echo-power from compressed video images: commercially available VueBox software by Bracco Suisse SA (Geneva, Switzerland) and in-laboratory software called PixPower. The echo-power estimated by CHI-Q analysis indicated a strong linear relationship versus agent concentration in vitro (R(2) ≥ 0.9996) for dynamic range (DR) settings of DR60 and DR80, with slopes between 9.22 and 9.57 dB/decade (p = 0.05). These values approach the theoretically predicted dependence of 10.0 dB/decade (equivalent to 3 dB for each concentration doubling). Echo-power estimations obtained from compressed video images with VueBox and PixPower also exhibited strong linear proportionality with concentration (R(2) ≥ 0.9996), with slopes between 9.30 and 9.68 dB/decade (p = 0.05). On an independent in vivo data set (N = 24), the difference in echo-power estimation between CHI-Q and each of the other two approaches was calculated after excluding regions that contain pixels affected by saturated or thresholded pixel values. The mean difference in estimates (expressed in decibels) was -0.25 dB between VueBox and CHI-Q (95% confidence interval: -0.75 to 0.26 dB) and -0.17 dB between PixPower and CHI-Q (95% confidence interval: -0.67 to 0.13 dB). To achieve linearization of data, one of the approaches (VueBox) requires calibration files provided by the software manufacturer for each machine type and setting. The other (PixPower) requires empirical correction of the imaging dynamic range based on ground truth data. These requirements could potentially be removed if US system manufacturers were willing to make relevant information on the applied processing publically available. Reliable echo-power estimation from linearized data would facilitate inclusion of different US systems in multicentric studies and more widespread implementation of emerging techniques for quantitative analysis of contrast ultrasound.

Three-dimensional quantitative ultrasound for detecting lymph node metastases.

PURPOSE: Detection of metastases in lymph nodes (LNs) is critical for cancer management. Conventional histological methods may miss metastatic foci. To date, no practical means of evaluating the entire LN volume exists. The aim of this study was to develop fast, reliable, operator-independent, high-frequency, quantitative ultrasound (QUS) methods for evaluating LNs over their entire volume to effectively detect LN metastases. METHODS: We scanned freshly excised LNs at 26 MHz and digitally acquired echo-signal data over the entire three-dimensional (3D) volume. A total of 146 LNs of colorectal, 26 LNs of gastric, and 118 LNs of breast cancer patients were enrolled. We step-sectioned LNs at 50-µm intervals and later compared them with 13 QUS estimates associated with tissue microstructure. Linear-discriminant analysis classified LNs as metastatic or nonmetastatic, and we computed areas (Az) under receiver-operator characteristic curves to assess classification performance. The QUS estimates and cancer probability values derived from discriminant analysis were depicted in 3D images for comparison with 3D histology. RESULTS: Of 146 LNs of colorectal cancer patients, 23 were metastatic; Az = 0.952 ± 0.021 (95% confidence interval [CI]: 0.911-0.993); sensitivity = 91.3% (specificity = 87.0%); and sensitivity = 100% (specificity = 67.5%). Of 26 LNs of gastric cancer patients, five were metastatic; Az = 0.962 ± 0.039 (95% CI: 0.807-1.000); sensitivity = 100% (specificity = 95.3%). A total of 17 of 118 LNs of breast cancer patients were metastatic; Az = 0.833 ± 0.047 (95% CI: 0.741-0.926); sensitivity = 88.2% (specificity = 62.5%); sensitivity = 100% (specificity = 50.5%). 3D cancer probability images showed good correlation with 3D histology. CONCLUSIONS: These results suggest that operator- and system-independent QUS methods allow reliable entire-volume LN evaluation for detecting metastases. 3D cancer probability images can help pathologists identify metastatic foci that could be missed using conventional methods.

Three-dimensional quantitative ultrasound to guide pathologists towards metastatic foci in lymph nodes.

The detection of metastases in freshly-excised lymph nodes from cancer patients during lymphadenectomy is critically important for cancer staging, treatment, and optimal patient management. Currently, conventional histologic methods suffer a high rate of false-negative determinations because pathologists cannot evaluate each excised lymph nodes in its entirety. Therefore, lymph nodes are undersampled and and small but clinically relevant metastatic regions can be missed. In this study, quantitative ultrasound (QUS) methods using high-frequency transducers (i.e., > 20 MHz) were developed and evaluated for their ability to detect and guide pathologists towards suspicious regions in lymph nodes. A custom laboratory scanning system was used to acquire radio-frequency (RF) data in 3D from excised lymph nodes using a 26-MHz center-frequency transducer. Overlapping 1-mm cylindrical regions-of-interest (ROIs) of the RF data were processed to yield 13 QUS estimates quantifying tissue microstructure and organization. These QUS methods were applied to more than 260 nodes from more than 160 colorectal-, gastric-, and breast-cancer patients. Cancer-detection performance was assessed for individual estimates and linear combinations of estimates. ROC results demonstrated excellent classification. For colorectal- and gastric-cancer nodes, the areas under the ROC curves (AUCs) were greater than 0.95. Slightly poorer results (AUC=0.85) were obtained for breast-cancer nodes. Images based on QUS parameters also permitted localization of cancer foci in some micrometastatic cases.

Three-dimensional high-frequency backscatter and envelope quantification of cancerous human lymph nodes.

Quantitative imaging methods using high-frequency ultrasound (HFU) offer a means of characterizing biological tissue at the microscopic level. Previously, high-frequency, 3-D quantitative ultrasound (QUS) methods were developed to characterize 46 freshly-dissected lymph nodes of colorectal-cancer patients. 3-D ultrasound radiofrequency data were acquired using a 25.6 MHz center-frequency transducer and each node was inked before tissue fixation to recover orientation after sectioning for 3-D histological evaluation. Backscattered echo signals were processed using 3-D cylindrical regions-of-interest (ROIs) to yield four QUS estimates associated with tissue microstructure (i.e., effective scatterer size, acoustic concentration, intercept and slope). These QUS estimates, obtained by parameterizing the backscatter spectrum, showed great potential for cancer detection. In the present study, these QUS methods were applied to 112 lymph nodes from 77 colorectal and gastric cancer patients. Novel QUS methods parameterizing the envelope statistics of the ROIs using Nakagami and homodyned-K distributions were also developed; they yielded four additional QUS estimates. The ability of these eight QUS estimates to classify lymph nodes and detect cancer was evaluated using receiver operating characteristics (ROC) curves. An area under the ROC curve of 0.996 with specificity and sensitivity of 95% were obtained by combining effective scatterer size and one envelope parameter based on the homodyned-K distribution. Therefore, these advanced 3-D QUS methods potentially can be valuable for detecting small metastatic foci in dissected lymph nodes.

Three-dimensional high-frequency characterization of cancerous lymph nodes.

High-frequency ultrasound (HFU) offers a means of investigating biologic tissue at the microscopic level. High-frequency, three-dimensional (3-D) quantitative-ultrasound (QUS) methods were developed to characterize freshly-dissected lymph nodes of cancer patients. Three-dimensional ultrasound data were acquired from lymph nodes using a 25.6-MHz center-frequency transducer. Each node was inked prior to tissue fixation to recover orientation after sectioning for 3-D histologic evaluation. Backscattered echo signals were processed using 3-D cylindrical regions-of-interest to yield four QUS estimates associated with tissue microstructure (i.e., effective scatterer size, acoustic concentration, intercept and slope). QUS estimates were computed following established methods using two scattering models. In this study, 46 lymph nodes acquired from 27 patients diagnosed with colon cancer were processed. Results revealed that fully-metastatic nodes could be perfectly differentiated from cancer-free nodes using slope or scatterer-size estimates. Specifically, results indicated that metastatic nodes had an average effective scatterer size (i.e., 37.1 +/- 1.7 microm) significantly larger (p < 0.05) than that in cancer-free nodes (i.e., 26 +/- 3.3 microm). Therefore, the 3-D QUS methods could provide a useful means of identifying small metastatic foci in dissected lymph nodes that might not be detectable using current standard pathology procedures.

Ultrasonic backscatter and attenuation (11-27 MHz) variation with collagen fiber distribution in ex vivo human dermis.

This ex vivo study explores the relationship of ultrasonic attenuation and backscatter to dermal microarchitecture by comparing ultrasonic measurements of these parameters (11-27 MHz) to a microscopic analysis of three parameters describing the collagen distribution (mean thickness and spacing of collagen bundles along the insonification direction and the percent area occupied by collagen). Skin samples (N= 31) were obtained from patients undergoing breast or abdominal reduction surgery. Radio-frequency (rf) signals were acquired in a B-scan format using an ultrasound system developed for skin imaging (Ultrasons Technologies, Tours, France). Ultrasonic data were analyzed to calculate average integrated backscatter (IBS in dB) and frequency dependence of backscatter (n, dimensionless) of each specimen at depths centered approximately 370,620 and 880 microm beneath the skin surface. Average integrated attenuation coefficient (IA in dB.cm(-1)) and frequency dependence of attenuation coefficient (beta in dB.cm(-1).MHz(-1)) were estimated across the depth between 240 and 1,000 microm. The three collagen distribution parameters were estimated using digitized microcopic fields from matched regions of histological sections stained with hematoxylin-eosin-saffron. No significant correlation was identified between collagen distribution parameters and IA or beta. For the most superficial depth studied in abdominal skin, n was inversely correlated to collagen bundle thickness (r = -0.67,p = 0.002) and percent area (r = -0.65,p = 0.003). At the same depth, IBS was inversely correlated to percent area of collagen (r = -0.51,p = 0.03). The rather high collagen packing (48 to 82% area) measured in histological sections and the inverse relationship observed between IBS and percent area of collagen suggest that a packing factor should be included in models relating skin collagen distribution to ultrasound spectral parameters. A better understanding of the relationship between ultrasound parameters and the microarchitecture of the dermis should help to interpret changes in ultrasonic parameters observed during in vivo ultrasonic skin examinations.

Singular spectrum analysis applied to backscattered ultrasound signals from in vitro human cancellous bone specimens.

Mean scatterer spacing (MSS) holds particular promise for the detection of changes in quasiperiodic tissue microstructures such as may occur during development of disease in the liver, spleen, or bones. Many techniques that may be applied for MSS estimation (temporal and spectral autocorrelation, power spectrum and cepstrum, higher order statistics, and quadratic transformation) characterize signals that contain a mixture of periodic and nonperiodic contributions. In contrast, singular spectrum analysis (SSA), a method usually applied in nonlinear dynamics, first identifies components of signals corresponding to periodic structures and, second, identifies dominant periodicity. Thus, SSA may better separate periodic structures from nonperiodic structures and noise. Using an ultrasound echo simulation model, we previously demonstrated SSA's potential to identify MSS of structures in quasiperiodic scattering media. The current work aims to observe the behavior of MSS estimation by SSA using ultrasound measurements in phantom materials (two parallel, nylon-line phantoms and four foam phantoms of different densities). The SSA was able to estimate not only the nylon-line distances but also nylon-line thickness. The method also was sensitive to the average pore-size differences of the four sponges. The algorithms then were applied to characterize human cancellous bone microarchitectures. Using 1-MHz center-frequency, radio-frequency ultrasound signals, MSS was measured in 24 in vitro bone samples and ranged from 1.0 to 1.7 mm. The SSA MSS estimates correlate significantly to MSS measured independently from synchrotron microtomography, r2 = 0.68. Thus, application of SSA to backscattered ultrasound signals seems to be useful for providing information linked to tissue microarchitecture that is not evident from clinical images.

Optimization of attenuation estimation in reflection for in vivo human dermis characterization at 20 MHz.

In vivo skin attenuation estimators must be applicable to backscattered radio frequency signals obtained in a pulse-echo configuration. This work compares three such estimators: short-time Fourier multinarrowband (MNB), short-time Fourier centroid shift (FC), and autoregressive centroid shift (ARC). All provide estimations of the attenuation slope (beta, dB x cm(-1) x MHz(-1)); MNB also provides an independent estimation of the mean attenuation level (IA, dB x cm(-1)). Practical approaches are proposed for data windowing, spectral variance characterization, and bandwidth selection. Then, based on simulated data, FC and ARC were selected as the best (compromise between bias and variance) attenuation slope estimators. The FC, ARC, and MNB were applied to in vivo human skin data acquired at 20 MHz to estimate betaFC, betaARC, and IA(MNB), respectively (without diffraction correction, between 11 and 27 MHz). Lateral heterogeneity had less effect and day-to-day reproducibility was smaller for IA than for beta. The IA and betaARC were dependent on pressure applied to skin during acquisition and IA on room and skin-surface temperatures. Negative values of IA imply that IA and beta may be influenced not only by skin's attenuation but also by structural heterogeneity across dermal depth. Even so, IA was correlated to subject age and IA, betaFC, and betaARC were dependent on subject gender. Thus, in vivo attenuation measurements reveal interesting variations with subject age and gender and thus appeared promising to detect skin structure modifications.