Pixel-based analysis of pulmonary changes on CT lung images due to COVID-19 pneumonia

Objectives: Computed tomography (CT) plays a complementary role in the diagnosis of the pneumonia-burden of COVID-19 disease. However, the low contrast of areas of inflammation on CT images, areas of infection are difficult to identify. The purpose of this study is to develop a post-image-processing method for quantitative analysis of COVID-19 pneumonia-related changes in CT attenuation values using a pixel-based analysis rather than more commonly used clustered focal pneumonia volumes. The COVID-19 pneumonia burden is determined by experienced radiologists in the clinic. Previous AI software was developed for the measurement of COVID-19 lesions based on the extraction of local pneumonia features. In this respect, changes in the pixel levels beyond the clusters may be overlooked by deep learning algorithms. The proposed technique focuses on the quantitative measurement of COVID-19 related pneumonia over the entire lung in pixel-by-pixel fashion rather than only clustered focal pneumonia volumes.

Material and Methods: Fifty COVID-19 and 50 age-matched negative control patients were analyzed using the proposed technique and commercially available artificial intelligence (AI) software. The %pneumonia was calculated using the relative volume of parenchymal pixels within an empirically defined CT density range, excluding pulmonary airways, vessels, and fissures. One-way ANOVA analysis was used to investigate the statistical difference between lobar and whole lung %pneumonia in the negative control and COVID-19 cohorts.

Results: The threshold of high-and-low CT attenuation values related to pneumonia caused by COVID-19 were found to be between ₋642.4 HU and 143 HU. The %pneumonia of the whole lung, left upper, and lower lobes were 8.1 ± 4.4%, 6.1 ± 4.5, and 11.3 ± 7.3% for the COVID-19 cohort, respectively, and statistically different (P < 0.01). Additionally, the pixel-based methods correlate well with existing AI methods and are approximately four times more sensitive to pneumonia particularly at the upper lobes compared with commercial software in COVID-19 patients (P < 0.01).

Conclusion: Pixel-by-pixel analysis can accurately assess pneumonia in COVID-19 patients with CT. Pixel-based techniques produce more sensitive results than AI techniques. Using the proposed novel technique, %pneumonia could be quantitatively calculated not only in the clusters but also in the whole lung with an improved sensitivity by a factor of four compared to AI-based analysis. More significantly, pixel-by-pixel analysis was more sensitive to the upper lobe pneumonia, while AI-based analysis overlooked the upper lung pneumonia region. In the future, this technique can be used to investigate the efficiency of vaccines and drugs and post COVID-19 effects.

Impaired pulmonary ventilation beyond pneumonia in COVID-19: A preliminary observation

Background

Coronavirus disease 2019 (COVID-19) may severely impair pulmonary function and cause hypoxia. However, the association of COVID-19 pneumonia on CT with impaired ventilation remains unexplained. This pilot study aims to demonstrate the relationship between the radiological findings on COVID-19 CT images and ventilation abnormalities simulated in a computational model linked to the patients’ symptoms.

Methods

Twenty-five patients with COVID-19 and four test-negative healthy controls who underwent a baseline non-enhanced CT scan: 7 dyspneic patients, 9 symptomatic patients without dyspnea, and 9 asymptomatic patients were included. A 2D U-Net-based CT segmentation software was used to quantify radiological futures of COVID-19 pneumonia. The CT image-based full-scale airway network (FAN) flow model was employed to assess regional lung ventilation. Functional and radiological features were compared across groups and correlated with the clinical symptoms. Heterogeneity in ventilation distribution and ventilation defects associated with the pneumonia and the patients’ symptoms were assessed.

Results

Median percentage ventilation defects were 0.2% for healthy controls, 0.7% for asymptomatic patients, 1.2% for symptomatic patients without dyspnea, and 11.3% for dyspneic patients. The median of percentage pneumonia was 13.2% for dyspneic patients and 0% for the other groups. Ventilation defects preferentially affected the posterior lung and worsened with increasing pneumonia linearly (y = 0.91x + 0.99, R² = 0.73) except for one of the nine dyspneic patients who had disproportionally large ventilation defects (7.8% of the entire lung) despite mild pneumonia (1.2%). The symptomatic and dyspneic patients showed significantly right-skewed ventilation distributions (symptomatic without dyspnea: 0.86 ± 0.61, dyspnea 0.91 ± 0.79) compared to the patients without symptom (0.45 ± 0.35). The ventilation defect analysis with the FAN model provided a comparable diagnostic accuracy to the percentage pneumonia in identifying dyspneic patients (area under the receiver operating characteristic curve, 0.94 versus 0.96).

Conclusions

COVID-19 pneumonia segmentations from CT scans are accompanied by impaired pulmonary ventilation preferentially in dyspneic patients. Ventilation analysis with CT image-based computational modelling shows it is able to assess functional impairment in COVID-19 and potentially identify one of the aetiologies of hypoxia in patients with COVID-19 pneumonia.

Imaging Dynamic Expiration: Feasibility of MRI Spirometry Using Hyperpolarized Xenon Gas

Purpose

To examine the feasibility of imaging-based spirometry using high-temporal-resolution projection MRI and hyperpolarized xenon 129 (¹²⁹Xe) gas.

Materials and Methods

In this prospective exploratory study, five healthy participants (age range, 25-45 years; three men) underwent an MRI spirometry technique using inhaled hyperpolarized ¹²⁹Xe and rapid two-dimensional projection MRI. Participants inhaled ¹²⁹Xe, then performed a forced expiratory maneuver while in an MR imager. Images of the lungs during expiration were captured in time intervals as short as 250 msec. Volume-corrected images of the lungs at expiration commencement (0 second), 1 second after expiration, and 6 seconds after expiration were extracted to generate forced expiratory volume in 1 second (FEV₁), forced vital capacity (FVC), and FEV₁/FVC ratio pulmonary maps. For comparison, participants performed conventional spirometry in the sitting position using room air, in the supine position using room air, and in the supine position using a room air and ¹²⁹Xe mixture. Paired t tests with Bonferroni corrections for multiple comparisons were used for statistical analyses.

Results

The mean MRI-derived FEV₁/FVC value was lower in comparison with conventional spirometry (0.52 ± 0.03 vs 0.70 ± 0.05, P < .01), which may reflect selective ¹²⁹Xe retention. A secondary finding of this study was that 1 L of inhaled ¹²⁹Xe negatively impacted pulmonary function as measured by conventional spirometry (in supine position), which reduced measured FEV₁ (2.70 ± 0.90 vs 3.04 ± 0.85, P < .01) and FEV₁/FVC (0.70 ± 0.05 vs 0.79 ± 0.04, P < .01).

Conclusion

A forced expiratory maneuver was successfully imaged with hyperpolarized ¹²⁹Xe and high-temporal-resolution MRI. Derivation of regional lung spirometric maps was feasible.Keywords: MR-Imaging, MR-Dynamic Contrast Enhanced, MR-Functional Imaging, Pulmonary, Thorax, Diaphragm, Lung, Pleura, Physics Supplemental material is available for this article. © RSNA, 2021.

Lobar Ventilation in Patients with COPD Assessed with the Full-Scale Airway Network Flow Model and Xenon-enhanced Dual-Energy CT

Background The full-scale airway network (FAN) flow model shows excellent agreement with limited functional imaging data but requires further validation prior to clinical use. Purpose To validate the ventilation distributions computed with the FAN flow model with xenon ventilation from xenon-enhanced dual-energy (DE) CT in participants with chronic obstructive pulmonary disease (COPD). Materials and Methods In this prospective study, the FAN model extracted structural data from xenon-enhanced DE CT images of men with COPD scanned between June 2012 and July 2013 to compute gas ventilation dynamics. The ventilation distributions on the middle cross-section plane, percentage lobar ventilation, and ventilation heterogeneity quantified by the coefficient of variation (CV) were compared between xenon-enhanced DE CT imaging and the FAN model. The relationship between the ventilation parameters with the densitometry and pulmonary function test results was demonstrated. The agreements and correlations between the parameters were measured using the concordance correlation coefficient and the Pearson correlation coefficient. Results Twenty-two men with COPD (mean age, 67 years ± 7 [standard deviation]) were evaluated. The percentage lobar ventilation computed with FAN showed a strong positive correlation with xenon-enhanced DE CT data (r = 0.7, P < .001). Ninety-five percent of lobar ventilation CV differences lay within 95% confidence intervals. Correlations of the percentage lobar ventilation were negative for percentage emphysema (xenon-enhanced DE CT: r = -0.38, P < .001; FAN: r = -0.23, P = .02) but were positive for percentage normal tissue volume (xenon-enhanced DE CT: r = 0.78, P < .001; FAN: r = 0.45, P < .001). Lung CVs of FAN revealed negative correlations with the spirometry results (CV_FAN vs percentage predicted forced expiratory volume in 1 second: r = -0.75, P < .001; CV_FAN vs ratio of forced expiratory volume in 1 second to forced vital capacity: r = -0.67, P < .001). Conclusion The full-scale airway network modeled lobar ventilation in patients with chronic obstructive pulmonary disease correlated with the xenon-enhanced dual-energy CT imaging data. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Parraga and Eddy in this issue.

Delayed ventilation assessment using fast dynamic hyperpolarised Xenon-129 magnetic resonance imaging

OBJECTIVES

To investigate the use of a fast dynamic hyperpolarised 129Xe ventilation magnetic resonance imaging (DXeV-MRI) method for detecting and quantifying delayed ventilation in patients with chronic obstructive pulmonary disease (COPD).

METHODS

Three male participants (age range 31-43) with healthy lungs and 15 patients (M/F = 12:3, age range = 48-73) with COPD (stages II-IV) underwent spirometry tests, quantitative chest computed tomography (QCT), and DXeV-MRI at 1.5-Tesla. Regional delayed ventilation was captured by measuring the temporal signal change in each lung region of interest (ROI) in comparison to that in the trachea. In addition to its qualitative assessment through visual inspection by a clinical radiologist, delayed ventilation was quantitatively captured by calculating a covariance measurement of the lung ROI and trachea signals, and quantified using both the time delay, and the difference between the integrated areas covered by the signal-time curves of the two signals.

RESULTS

Regional temporal ventilation, consistent with the expected physiological changes across a free breathing cycle, was demonstrated with DXeV-MRI in all patients. Delayed ventilation was observed in 13 of the 15 COPD patients and involved variable lung ROIs. This was in contrast to the control group, where no delayed ventilation was demonstrated (p = 0.0173).

CONCLUSIONS

DXeV-MRI offers a non-invasive way of detecting and quantifying delayed ventilation in patients with COPD, and provides physiological information on regional pulmonary function during a full breathing cycle.

KEY POINTS

• Dynamic xenon MRI allows for the non-invasive detection and measurement of delayed ventilation in COPD patients. • Dynamic xenon MRI during a free breathing cycle can provide unique information about pulmonary physiology and pulmonary disease pathophysiology. • With further validation, dynamic xenon MRI could offer a non-invasive way of measuring collateral ventilation which can then be used to guide lung volume reduction therapy (LVRT) for certain COPD patients.

CT-based Airway Flow Model to Assess Ventilation in Chronic Obstructive Pulmonary Disease: A Pilot Study

Background The lack of functional information in thoracic CT remains a limitation of its use in the clinical management of chronic obstructive pulmonary disease (COPD). Purpose To compare the distribution of pulmonary ventilation assessed by a CT-based full-scale airway network (FAN) flow model with hyperpolarized xenon 129 (¹²⁹Xe) MRI (hereafter, ¹²⁹Xe MRI) and technetium 99m-diethylenetriaminepentaacetic acid aerosol SPECT ventilation imaging (hereafter, V-SPECT) in participants with COPD. Materials and Methods In this prospective study performed between May and August 2017, pulmonary ventilation in participants with COPD was computed by using the FAN flow model. The modeled pulmonary ventilation was compared with functional imaging data from breath-hold time-series ¹²⁹Xe MRI and V-SPECT. FAN-derived ventilation images on the coronal plane and volumes of interest were compared with functional lung images. Percentage lobar ventilation estimated by the FAN model was compared with that measured at ¹²⁹Xe MRI and V-SPECT. The statistical significance of ventilation distribution between FAN and functional images was demonstrated with the Spearman correlation coefficient and χ² distance. Results For this study, nine participants (seven men [mean age, 65 years ± 5 {standard deviation}] and two women [mean age, 63 years ± 7]) with COPD that was Global Initiative for Chronic Obstructive Lung Disease stage II-IV were enrolled. FAN-modeled ventilation profile showed strong positive correlation with images from ¹²⁹Xe MRI (ρ = 0.67; P < .001) and V-SPECT (ρ = 0.65; P < .001). The χ² distances of the ventilation histograms in the volumes of interest between the FAN and ¹²⁹Xe MRI and FAN and V-SPECT were 0.16 ± 0.08 and 0.28 ± 0.14, respectively. The ratios of lobar ventilations in the models were linearly correlated to images from ¹²⁹Xe MRI (ρ = 0.67; P < .001) and V-SPECT (ρ = 0.59; P < .001). Conclusion A CT-based full-scale airway network flow model provided regional pulmonary ventilation information for chronic obstructive pulmonary disease and correlates with hyperpolarized xenon 129 MRI and technetium 99m-diethylenetriaminepentaacetic acid aerosol SPECT ventilation imaging. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Schiebler and Parraga in this issue.

Time-series hyperpolarized xenon-129 MRI of lobar lung ventilation of COPD in comparison to V/Q-SPECT/CT and CT

Purpose

To derive lobar ventilation in patients with chronic obstructive pulmonary disease (COPD) using a rapid time-series hyperpolarized xenon-129 (HPX) magnetic resonance imaging (MRI) technique and compare this to ventilation/perfusion single-photon emission computed tomography (V/Q-SPECT), correlating the results with high-resolution computed tomography (CT) and pulmonary function tests (PFTs).

Materials and methods

Twelve COPD subjects (GOLD stages I–IV) participated in this study and underwent HPX-MRI, V/Q-SPECT/CT, high-resolution CT, and PFTs. HPX-MRI was performed using a novel time-series spiral k-space sampling approach. Relative percentage ventilations were calculated for individual lobe for comparison to the relative SPECT lobar ventilation and perfusion. The absolute HPX-MRI percentage ventilation in each lobe was compared to the absolute CT percentage emphysema score calculated using a signal threshold method. Pearson’s correlation and linear regression tests were performed to compare each imaging modality.

Results

Strong correlations were found between the relative lobar percentage ventilation with HPX-MRI and percentage ventilation SPECT (r = 0.644; p < 0.001) and percentage perfusion SPECT (r = 0.767; p < 0.001). The absolute CT percentage emphysema and HPX percentage ventilation correlation was also statistically significant (r = 0.695, p < 0.001). The whole lung HPX percentage ventilation correlated with the PFT measurements (FEV₁ with r = − 0.886, p < 0.001*, and FEV₁/FVC with r = − 0.861, p < 0.001*) better than the whole lung CT percentage emphysema score (FEV₁ with r = − 0.635, p = 0.027; and FEV₁/FVC with r = − 0.652, p = 0.021).

Conclusion

Lobar ventilation with HPX-MRI showed a strong correlation with lobar ventilation and perfusion measurements derived from SPECT/CT, and is better than the emphysema score obtained with high-resolution CT.

Key Points

• The ventilation hyperpolarized xenon-129 MRI correlates well with ventilation and perfusion with SPECT/CT with the advantage of higher temporal and spatial resolution.

• The hyperpolarized xenon-129 MRI correlates with the PFT measurements better than the high-resolution CT with the advantage of avoiding the use of ionizing radiation.

Electronic supplementary material

The online version of this article (10.1007/s00330-018-5888-y) contains supplementary material, which is available to authorized users.

Magnetic resonance imaging of the time course of hyperpolarized 129Xe gas exchange in the human lungs and heart

Purpose

To perform magnetic resonance imaging (MRI), human lung imaging, and quantification of the gas-transfer dynamics of hyperpolarized xenon-129 (HPX) from the alveoli into the blood plasma.

Materials and methods

HPX MRI with iterative decomposition of water and fat with echo asymmetry and least-square estimation (IDEAL) approach were used with multi-interleaved spiral k-space sampling to obtain HPX gas and dissolved phase images. IDEAL time-series images were then obtained from ten subjects including six normal subjects and four patients with pulmonary emphysema to test the feasibility of the proposed technique for capturing xenon-129 gas-transfer dynamics (XGTD). The dynamics of xenon gas diffusion over the entire lung was also investigated by measuring the signal intensity variations between three regions of interest, including the left and right lungs and the heart using Welch’s t test.

Results

The technique enabled the acquisition of HPX gas and dissolved phase compartment images in a single breath-hold interval of 8 s. The y-intersect of the XGTD curves were also found to be statistically lower in the patients with lung emphysema than in the healthy group (p < 0.05).

Conclusion

This time-series IDEAL technique enables the visualization and quantification of inhaled xenon from the alveoli to the left ventricle with a clinical gradient strength magnet during a single breath-hold, in healthy and diseased lungs.

Key Points

• The proposed hyperpolarized xenon-129 gas and dissolved magnetic resonance imaging technique can provide regional and temporal measurements of xenon-129 gas-transfer dynamics.

• Quantitative measurement of xenon-129 gas-transfer dynamics from the alveolar to the heart was demonstrated in normal subjects and pulmonary emphysema.

• Comparison of gas-transfer dynamics in normal subjects and pulmonary emphysema showed that the proposed technique appears sensitive to changes affecting the alveoli, pulmonary interstitium, and capillaries.

Electronic supplementary material

The online version of this article (10.1007/s00330-018-5853-9) contains supplementary material, which is available to authorized users.

Fast dynamic ventilation MRI of hyperpolarized 129 Xe using spiral imaging

Purpose

To develop and optimize a rapid dynamic hyperpolarized ¹²⁹Xe ventilation (DXeV) MRI protocol and investigate the feasibility of capturing pulmonary signal‐time curves in human lungs.

Theory and Methods

Spiral k‐space trajectories were designed with the number of interleaves N_int = 1, 2, 4, and 8 corresponding to voxel sizes of 8 mm, 5 mm, 4 mm, and 2.5 mm, respectively, for field of view = 15 cm. DXeV images were acquired from a gas‐flow phantom to investigate the ability of N_int = 1, 2, 4, and 8 to capture signal‐time curves. A finite element model was constructed to investigate gas‐flow dynamics corroborating the experimental signal‐time curves. DXeV images were also carried out in six subjects (three healthy and three chronic obstructive pulmonary disease subjects).

Results

DXeV images and numerical modelling of signal‐time curves permitted the quantification of temporal and spatial resolutions for different numbers of spiral interleaves. The two‐interleaved spiral (N_int = 2) was found to be the most time‐efficient to obtain DXeV images and signal‐time curves of whole lungs with a temporal resolution of 624 ms for 13 slices. Signal‐time curves were well matched in three healthy volunteers. The Spearman's correlations of chronic obstructive pulmonary disease subjects were statistically different from three healthy subjects (P < 0.05).

Conclusion

The N_int = 2 spiral demonstrates the successful acquisition of DXeV images and signal‐time curves in healthy subjects and chronic obstructive pulmonary disease patients. Magn Reson Med 79:2597–2606, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Quantification of regional early stage gas exchange changes using hyperpolarized (129)Xe MRI in a rat model of radiation-induced lung injury

PURPOSE

To assess the feasibility of hyperpolarized (HP) (129)Xe MRI for detection of early stage radiation-induced lung injury (RILI) in a rat model involving unilateral irradiation by assessing differences in gas exchange dynamics between irradiated and unirradiated lungs.

METHODS

The dynamics of gas exchange between alveolar air space and pulmonary tissue (PT), PT and red blood cells (RBCs) was measured using single-shot spiral iterative decomposition of water and fat with echo asymmetry and least-squares estimation images of the right and left lungs of two age-matched cohorts of Sprague Dawley rats. The first cohort (n = 5) received 18 Gy irradiation to the right lung using a (60)Co source and the second cohort (n = 5) was not irradiated and served as the healthy control. Both groups were imaged two weeks following irradiation when radiation pneumonitis (RP) was expected to be present. The gas exchange data were fit to a theoretical gas exchange model to extract measurements of pulmonary tissue thickness (LPT) and relative blood volume (VRBC) from each of the right and left lungs of both cohorts. Following imaging, lung specimens were retrieved and percent tissue area (PTA) was assessed histologically to confirm RP and correlate with MRI measurements.

RESULTS

Statistically significant differences in LPT and VRBC were observed between the irradiated and non-irradiated cohorts. In particular, LPT of the right and left lungs was increased approximately 8.2% and 5.0% respectively in the irradiated cohort. Additionally, VRBC of the right and left lungs was decreased approximately 36.1% and 11.7% respectively for the irradiated cohort compared to the non-irradiated cohort. PTA measurements in both right and left lungs were increased in the irradiated group compared to the non-irradiated cohort for both the left (P < 0.05) and right lungs (P < 0.01) confirming the presence of RP. PTA measurements also correlated with the MRI measurements for both the non-irradiated (r = 0.79, P < 0.01) and irradiated groups (r = 0.91, P < 0.01).

CONCLUSIONS

Regional RILI can be detected two weeks post-irradiation using HP (129)Xe MRI and analysis of gas exchange curves. This approach correlates well with histology and can potentially be used clinically to assess radiation pneumonitis associated with early RILI to improve radiation therapy outcomes.

Hyperpolarized (129) Xe imaging of the rat lung using spiral IDEAL

PURPOSE

To implement and optimize a single-shot spiral encoding strategy for rapid 2D IDEAL projection imaging of hyperpolarized (Hp) (129) Xe in the gas phase, and in the pulmonary tissue (PT) and red blood cells (RBCs) compartments of the rat lung, respectively.

THEORY AND METHODS

A theoretical and experimental point spread function analysis was used to optimize the spiral k-space read-out time in a phantom. Hp (129) Xe IDEAL images from five healthy rats were used to: (i) optimize flip angles by a Bloch equation analysis using measured kinetics of gas exchange and (ii) investigate the feasibility of the approach to characterize the exchange of Hp (129) Xe.

RESULTS

A read-out time equal to approximately 1.8 × T2* was found to provide the best trade-off between spatial resolution and signal-to-noise ratio (SNR). Spiral IDEAL approaches that use the entire dissolved phase magnetization should give an SNR improvement of a factor of approximately three compared with Cartesian approaches with similar spatial resolution. The IDEAL strategy allowed imaging of gas, PT, and RBC compartments with sufficient SNR and temporal resolution to permit regional gas exchange measurements in healthy rats.

CONCLUSION

Single-shot spiral IDEAL imaging of gas, PT and RBC compartments and gas exchange is feasible in rat lung using Hp (129) Xe. Magn Reson Med 76:566-576, 2016. © 2015 Wiley Periodicals, Inc.

Imaging the electro-kinetic response of biological tissues with optical coherence tomography

We demonstrate the feasibility of using optical coherence tomography (OCT) to detect and image an electro-kinetic response: electric-field induced optical changes (EIOC) in soft biological tissues. A low-frequency electric field was applied to ex vivo samples of porcine heart tissues, while OCT signals were acquired continuously. Experimental results show that the amplitude of the OCT signal change is proportional to the amplitude and inversely proportional to the frequency of the applied electric field. We show that the nonconductive component of the sample was eliminated in the normalized EIOC image. To the best our knowledge, this is the first time a two-dimensional image related to the electro-kinetic response of soft tissues is obtained with depth resolution. Since electro-kinetic properties can change during cancerogenesis, EIOC imaging can potentially be used for cancer detection.

Reversibility of electric-field-induced mechanical changes in soft tissues

Recently, ultrasound has been used to study the physiological-level electric-field-induced mechanical changes (EIMC) in general soft biological tissues. Here, we present the experimental results on the dependence of EIMC on the polarity of the applied electric field. We applied an ac voltage source to heart tissues and monitored the amplitude changes and time shifting of the ultrasound echoes. The shifting of the echoes was decomposed into a trend component and a fluctuation (feature) component. The changes in amplitude and the fluctuation component of the time shift, but not the trend component of the time shift, can be mostly reversed by reversing the polarity of the applied voltage. The polarity-dependence study reveals two different mechanisms underlying EIMC.

Electric-field induced strain in biological tissues

This paper reports a new effect whereby a physiological-level direct-current electrical field (at 1.4 V/cm) can induce time-varying mechanical strain in various types of biological tissues and gel phantoms. This effect cannot be explained by the piezoelectric effect, tissue contraction, temperature changes, and electrorestriction. The induced strain in tissues was analyzed by processing ultrasound echo signals. The sample expanded perpendicularly to the applied electric field. The expansion rate depended on the history of the applied electric field. The speed of sound changed little compared with the expansion. The new effect might be related to electrokinetic effects.