Imaging the dynamics of cardiac fiber orientation in vivo using 3D Ultrasound Backscatter Tensor Imaging

Dynamics of multicellular swirling on micropatterned substrates

Chirality plays a crucial role in biology, as it is highly conserved and fundamentally important in the developmental process. To better understand the relationship between the chirality of individual cells and that of tissues and organisms, we develop a generalized mechanics model of chiral polarized particles to investigate the swirling dynamics of cell populations on substrates. Our analysis reveals that cells with the same chirality can form distinct chiral patterns on ring-shaped or rectangular substrates. Interestingly, our studies indicate that an excessively strong or weak individual cellular chirality hinders the formation of such chiral patterns. Our studies also indicate that there exists the influence distance of substrate boundaries in chiral patterns. Smaller influence distances are observed when cell-cell interactions are weaker. Conversely, when cell-cell interactions are too strong, multiple cells tend to be stacked together, preventing the formation of chiral patterns on substrates in our analysis. Additionally, we demonstrate that the interaction between cells and substrate boundaries effectively controls the chiral distribution of cellular orientations on ring-shaped substrates. This research highlights the significance of coordinating boundary features, individual cellular chirality, and cell-cell interactions in governing the chiral movement of cell populations and provides valuable mechanics insights into comprehending the intricate connection between the chirality of single cells and that of tissues and organisms.

Biopathologic Characterization and Grade Assessment of Breast Cancer With 3-D Multiparametric Ultrasound Combining Shear Wave Elastography and Backscatter Tensor Imaging

OBJECTIVE

Despite recent improvements in medical imaging, the final diagnosis and biopathologic characterization of breast cancers currently still requires biopsies. Ultrasound is commonly used for clinical examination of breast masses. B-Mode and shear wave elastography (SWE) are already widely used to detect suspicious masses and differentiate benign lesions from cancers. But additional ultrasound modalities such as backscatter tensor imaging (BTI) could provide relevant biomarkers related to tissue organization. Here we describe a 3-D multiparametric ultrasound approach applied to breast carcinomas in the aims of (i) validating the ability of BTI to reveal the underlying organization of collagen fibers and (ii) assessing the complementarity of SWE and BTI to reveal biopathologic features of diagnostic interest.

METHODS

Three-dimensional SWE and BTI were performed ex vivo on 64 human breast carcinoma samples using a linear ultrasound probe moved by a set of motors. Here we describe a 3-D multiparametric representation of the breast masses and quantitative measurements combining B-mode, SWE and BTI.

RESULTS

Our results reveal for the first time that BTI can capture the orientation of the collagen fibers around tumors. BTI was found to be a relevant marker for assessing cancer stages, revealing a more tangent tissue orientation for in situ carcinomas than for invasive cancers. In invasive cases, the combination of BTI and SWE parameters allowed for classification of invasive tumors with respect to their grade with an accuracy of 95.7%.

CONCLUSION

Our results highlight the potential of 3-D multiparametric ultrasound imaging for biopathologic characterization of breast tumors.

Biofabrication of Composite Bioink-Nanofiber Constructs: Effect of Rheological Properties of Bioinks on 3D (Bio)Printing and Cells Interaction with Aligned Touch Spun Nanofibers

This paper reports on a novel approach for the fabrication of composite multilayered bioink-nanofibers construct. This work achieves this by using a hands-free 3D (bio)printing integrated touch-spinning approach. Additionally, this work investigates the interaction of fibroblasts in different bioinks with the highly aligned touch-spun nanofibers. This work conducts a comprehensive characterization of the rheological properties of the inks, starting with low-strain oscillatory rheology to analyze the viscoelastic behavior, when the material structure remains intact. Moreover, this work performs amplitude sweeps to investigate the stability of the inks under large deformations, rotational rheology to examine the shear thinning profile, and a three-step creep experiment to study time-dependent rheological behavior. The obtained rheological results are correlated to visual observation of the flow behavior of inks. These behaviors span from an ink with zero-shear viscosity, very weak shear thinning, and no thixotropic behavior to inks exhibiting flow stress, pronounced shear thinning, and thixotropy. It is demonstrated that inks have an essential effect on cell behavior. While all bioinks allow a preferred directionality of the fibroblasts along the fiber direction, cells tend to form aggregates in bioinks with higher viscosity, and a considerable number of agglomerates are observed in the presence of laponite-RD.

Mimicking Transmural Helical Cardiomyofibre Orientation Using Bouligand-like Pore Structures in Ice-Templated Collagen Scaffolds

The helical arrangement of cardiac muscle fibres underpins the contractile properties of the heart chamber. Across the heart wall, the helical angle of the aligned fibres changes gradually across the range of 90-180°. It is essential to recreate this structural hierarchy in vitro for developing functional artificial tissue. Ice templating can achieve single-oriented pore alignment via unidirectional ice solidification with a flat base mould design. We hypothesise that the orientation of aligned pores can be controlled simply via base topography, and we propose a scalable base design to recapitulate the transmural fibre orientation. We have utilised finite element simulations for rapid testing of base designs, followed by experimental confirmation of the Bouligand-like orientation. X-ray microtomography of experimental samples showed a gradual shift of 106 ± 10°, with the flexibility to tailor pore size and spatial helical angle distribution for personalised medicine.

Collagenase treatment reduces the anisotropy of ultrasonic backscatter in rat myocardium by reducing collagen crosslinks

Dysregulation of collagen deposition, degradation, and crosslinking in the heart occur in response to increased physiological stress. Collagen content has been associated with ultrasonic backscatter (brightness), and we have shown that the anisotropy of backscatter can be used to measure myofiber alignment, that is, variation in the brightness of a left ventricular short-axis ultrasound. This study investigated collagen's role in anisotropy of ultrasonic backscatter; female Sprague-Dawley rat hearts were treated with a collagenase-containing solution, for either 10 or 30 min, or control solution for 30 min. Serial ultrasound images were acquired at 2.5-min intervals throughout collagenase treatment. Ultrasonic backscatter was assessed from anterior and posterior walls, where collagen fibrils are predominately aligned perpendicular to the angle of insonification, and the lateral and septal walls, where collagen is predominately aligned parallel to the angle of insonification. Collagenase digestion reduced backscatter anisotropy within the myocardium. Collagen remains present in the myocardium throughout collagenase treatment, but crosslinking is altered within 10 min. These data suggest that crosslinking of collagen modulates the anisotropy of ultrasonic backscatter. An Anisotropy Index, derived from differences in backscatter from parallel and perpendicularly aligned fibers, may provide a noninvasive index to monitor the progression and state of myocardial fibrosis.

Membrane polarization at the excitation threshold induced by external electric fields in cardiomyocytes of rats at different developmental stages

External electric fields (E), used for cardiac pacing and defibrillation/cardioversion, induce a spatially variable change in cardiomyocyte transmembrane potential (ΔVm) that depends on cell geometry and E orientation. This study investigates E-induced ΔVm in cardiomyocytes from rats at different ages, which show marked size/geometry variation. Using a tridimensional numerical electromagnetic model recently proposed (NM3D), it was possible: (a) to evaluate the suitability of the simpler, prolate spheroid analytical model (PSAM) to calculate amplitude and location of ΔVm maximum (ΔVmax) for E = 1 V.cm-1; and (b) to estimate the ΔVmax required for excitation (ΔVT) from experimentally determined threshold E values (ET). Ventricular myocytes were isolated from neonatal, weaning, adult, and aging Wistar rats. NM3D was constructed as the extruded 2D microscopy cell image, while measured minor and major cell dimensions were used for PSAM. Acceptable ΔVm estimates can be obtained with PSAM from paralelepidal cells for small θ. ET, but not ΔVT, was higher for neonate cells. ΔVT was significantly greater in the cell from older animals, which indicate lower responsiveness to E associated with aging, rather than with altered cell geometry/dimensions. ΔVT might be used as a non-invasive indicator of cell excitability as it is little affected by cell geometry/size.

Graft-host coupling changes can lead to engraftment arrhythmia: a computational study

After myocardial infarction (MI), a significant portion of heart muscle is replaced with scar tissue, progressively leading to heart failure. Human pluripotent stem cell-derived cardiomyocytes (hPSC-CM) offer a promising option for improving cardiac function after MI. However, hPSC-CM transplantation can lead to engraftment arrhythmia (EA). EA is a transient phenomenon arising shortly after transplantation then spontaneously resolving after a few weeks. The underlying mechanism of EA is unknown. We hypothesize that EA may be explained partially by time-varying, spatially heterogeneous, graft-host electrical coupling. Here, we created computational slice models derived from histological images that reflect different configuration of grafts in the infarcted ventricle. We ran simulations with varying degrees of connection imposed upon the graft-host perimeter to assess how heterogeneous electrical coupling affected EA with non-conductive scar, slow-conducting scar and scar replaced by host myocardium. We also quantified the effect of variation in intrinsic graft conductivity. Susceptibility to EA initially increased and subsequently decreased with increasing graft-host coupling, suggesting the waxing and waning of EA is regulated by progressive increases in graft-host coupling. Different spatial distributions of graft, host and scar yielded markedly different susceptibility curves. Computationally replacing non-conductive scar with host myocardium or slow-conducting scar, and increasing intrinsic graft conductivity both demonstrated potential means to blunt EA vulnerability. These data show how graft location, especially relative to scar, along with its dynamic electrical coupling to host, can influence EA burden; moreover, they offer a rational base for further studies aimed to define the optimal delivery of hPSC-CM injection. KEY POINTS: Human pluripotent stem cell-derived cardiomyocytes (hPSC-CM) hold great cardiac regenerative potential but can also cause engraftment arrhythmias (EA). Spatiotemporal evolution in the pattern of electrical coupling between injected hPSC-CMs and surrounding host myocardium may explain the dynamics of EA observed in large animal models. We conducted simulations in histology-derived 2D slice computational models to assess the effects of heterogeneous graft-host electrical coupling on EA propensity, with or without scar tissue. Our findings suggest spatiotemporally heterogeneous graft-host coupling can create an electrophysiological milieu that favours graft-initiated host excitation, a surrogate metric of EA susceptibility. Removing scar from our models reduced but did not abolish the propensity for this phenomenon. Conversely, reduced intra-graft electrical connectedness increased the incidence of graft-initiated host excitation. The computational framework created for this study can be used to generate new hypotheses, targeted delivery of hPSC-CMs.

Backscatter tensor imaging and 3D speckle tracking for simultaneous ex vivo structure and deformation measurement of myocardium

OBJECTIVE

Biaxial mechanical testing is a common method for elucidation of mechanical properties of excised ventricular myocardium, especially in the context of structural remodeling that accompanies heart disease. Current imaging strategies in biaxial testing are based on optical camera imaging of the tissue surface, thus providing no information about the tissue microstructure and limiting strain measurements to two dimensions. Here, these limitations are overcome by replacing the camera with ultrasound imaging in order to measure both transmural fiber orientation and 3D tissue deformation during biaxial testing.

METHODS

Quasi-static biaxial mechanical testing is applied to four samples of excised porcine ventricular myocardium (two left- and two right-ventricular tissues). During testing, a rotational scan of an ultrasound linear array provides data for both backscatter tensor imaging and 3D speckle tracking, from which transmural fiber orientation and tissue deformation are computed, respectively. Ultrasound-derived fiber orientation and tissue strain are validated against histology and camera surface imaging, respectively.

DISCUSSION

Ultrasound-derived fiber angle and tissue strain exhibit good accuracy, with root-mean-square errors of 9.9° and 1.2% strain, respectively. Further investigation into the optimization of backscatter tensor imaging is warranted. Replacing the rotational scan of a linear array with volume imaging with a matrix array will improve the technique.

CONCLUSION

Ultrasound imaging can replace the optical camera measurement during biaxial mechanical testing of ventricular myocardium in order to accurately provide measurements of transmural fiber orientation and tissue strain. In situ knowledge of transmural fiber structure and tissue deformation can enhance the inverse problem used to determine tissue mechanical properties from biaxial testing.

Impact of Ventricular Geometric Characteristics on Myocardial Stiffness Assessment Using Shear-Wave Velocity in Healthy Children and Young Adults

BACKGROUND

Diastolic myocardial stiffness (MS) can serve as a key diagnostic parameter for congenital or acquired heart diseases. Using shear modulus and shear-wave velocity (SWV), shear-wave elastography (SWE) is an emerging ultrasound-based technique that can allow noninvasive assessment of MS. However, MS extrinsic parameters such as left ventricular geometric characteristics could affect shear-wave propagation. The aims of this study were to determine a range of normal values of MS using SWE in age groups of healthy children and young adults and to explore the impact of left ventricular geometric characteristics on SWE.

METHODS

Sixty healthy volunteers were recruited in the study and divided into 2 groups: neonates (0-1 months old, n = 15) and >1 month old (1 month to 45 years of age, n = 45). SWE was performed using the Verasonics Vantage systems with a phased-array ultrasound probe. The anteroseptal basal segment was assessed in two views. SWE was electrocardiographically triggered during the end-diastolic phase. Conventional echocardiography was performed to assess ventricular function and anatomy. Results are presented as stiffness values along with mean velocity measurements and SDs. Simple and multivariate linear regression analyses were performed.

RESULTS

For neonates, mean MS was 1.87 ± 0.79 kPa (range, 0.59-2.91 kPa; mean SWV, 1.37 ± 0.57 m/sec), with high variability and no correlation with age (P = .239). For this age group, no statistically significant correlation was found between MS and any demographic or echocardiographic parameters (P > .05). For the >1 month old group, a mean MS value of 1.67 ± 0.53 kPa was observed (range, 0.6-3 kPa; mean SWV, 1.29 ± 0.49 m/sec) for healthy volunteers. When paired for age, no sex-related difference was observed (P = .55). In univariate linear regression analysis, age (r = 0.83, P < .01), diastolic interventricular septal thickness (r = 0.72, P < .01), and left ventricular end-diastolic diameter (r = 0.67, P < .01) were the parameters with the highest correlation coefficients with MS. In a multiple linear regression analysis incorporating these three parameters as cofounding factors, age was the only statistically significant parameters (r = 0.81, P = .02).

CONCLUSION

Diastolic MS increases linearly in children and young adults. Diastolic MS correlates more robustly with age than with myocardial and left ventricular geometric characteristics. However, the geometry affects SWV, implying the need to determine well-established boundaries in future studies for the clinical application of SWE.

Assessing engineered tissues and biomaterials using ultrasound imaging: In vitro and in vivo applications

Quantitative assessment of the structural, functional, and mechanical properties of engineered tissues and biomaterials is fundamental to their development for regenerative medicine applications. Ultrasound (US) imaging is a non-invasive, non-destructive, and cost-effective technique capable of longitudinal and quantitative monitoring of tissue structure and function across centimeter to sub-micron length scales. Here we present the fundamentals of US to contextualize its application for the assessment of biomaterials and engineered tissues, both in vivo and in vitro. We review key studies that demonstrate the versatility and broad capabilities of US for clinical and pre-clinical biomaterials research. Finally, we highlight emerging techniques that further extend the applications of US, including for ultrafast imaging of biomaterials and engineered tissues in vivo and functional monitoring of stem cells, organoids, and organ-on-a-chip systems in vitro.

A tomographic microscopy-compatible Langendorff system for the dynamic structural characterization of the cardiac cycle

Introduction

Cardiac architecture has been extensively investigated ex vivo using a broad spectrum of imaging techniques. Nevertheless, the heart is a dynamic system and the structural mechanisms governing the cardiac cycle can only be unveiled when investigating it as such.

Methods

This work presents the customization of an isolated, perfused heart system compatible with synchrotron-based X-ray phase contrast imaging (X-PCI).

Results

Thanks to the capabilities of the developed setup, it was possible to visualize a beating isolated, perfused rat heart for the very first time in 4D at an unprecedented 2.75 μm pixel size (10.6 μm spatial resolution), and 1 ms temporal resolution.

Discussion

The customized setup allows high-spatial resolution studies of heart architecture along the cardiac cycle and has thus the potential to serve as a tool for the characterization of the structural dynamics of the heart, including the effects of drugs and other substances able to modify the cardiac cycle.

Ultrasound-Based Estimation of Fibre-Directional Strain: A Simulation Study

Left ventricular (LV) strains are typically represented with respect to the imaging axes. Contraction within the myocardium occurs along myofibres, which vary in orientation. Therefore, a mismatch exists between the direction in which strain is calculated and that in which contraction occurs. In this study, ultrasound-based fibre orientation and 3-D strain estimation were combined to calculate the fibre-directional strain. Three-dimensional ultrasound volumes were created by simulating simple geometrical phantoms and a phantom based on a finite-element (FE) model of LV mechanics. Fibre-like structures were embedded within tissue-mimicking scatterers. Strains were applied to the numerical phantom, whereas the FE phantom was deformed based on the LV model. Fibre orientation was accurately estimated for both phantoms. There was poor agreement in axial and elevational strains (root mean square error = 29.9% and 12.3%), but good agreement in lateral and fibre-directional strains (root mean square error = 6.4% and 5.9% respectively), which aligned in the midwall. Simplifications to reduce computational complexity caused poor axial and elevational strain estimation. However, calculation of fibre-directional strain from single-modality ultrasound volumes was successful. Further studies, in ex vivo setups because of the fundamental limitations of currently available transducers, are needed to verify real-world performance of the method.

Refraction-corrected backscatter tensor imaging of excised porcine ventricular myocardium

Backscatter tensor imaging (BTI) is performed on excised porcine right- and left-ventricular myocardium to estimate the transmural myofiber orientation. Calculation of the backscatter spatial coherence employs measured sound speeds of the myocardium and the fluid that separates the tissue from the imaging array to account for effects of refraction during the delay-and-sum beamforming calculation. Compared to the assumption of a homogeneous sound speed in the imaging region, accounting for refraction yields significantly increased average spatial coherence as well as contrast of spatial coherence between the along- and across-fiber directions, thus improving sensitivity of BTI for myofiber orientation estimation.

Proof of Concept of 3-D Backscatter Tensor Imaging Tomography for Non-invasive Assessment of Human Breast Cancer Collagen Organization

Tumor growth, similarly to several other pathologies, tends to change the structural orientation of soft tissue fibers, which can become relevant markers for diagnosis. Current diagnosis protocols may require a biopsy for histological analysis, which is an invasive, painful and stressful procedure with a minimum turnaround time of 2 d. Otherwise, diagnosis may involve the use of complex methods with limited availability such as diffusion tensor imaging (magnetic resonance diffusion tensor imaging), which is not widely used in medical practice. Conversely, advanced methodologies in ultrasound imaging such as backscatter tensor imaging (BTI) might become a routine procedure in clinical practice at a limited cost. This method evaluates the local organization of soft tissues based on the spatial coherence of their backscattered ultrasonic echoes. Previous work has proven that BTI applied with matrix probes enables measurement of the orientation of soft tissue fibers, especially in the myocardium. The aims of the study described here were (i) to present for the first time a methodology for performing BTI in a volume on ex vivo human breast tumors using a linear probe and (ii) to display a first proof of concept of the link between BTI measurements and the orientation of collagen fibers.

Myocardial mesostructure and mesofunction

The complex and highly organized structural arrangement of some five billion cardiomyocytes directs the coordinated electrical activity and mechanical contraction of the human heart. The characteristic transmural change in cardiomyocyte orientation underlies base-to-apex shortening, circumferential shortening, and left ventricular torsion during contraction. Individual cardiomyocytes shorten ∼15% and increase in diameter ∼8%. Remarkably, however, the left ventricular wall thickens by up to 30-40%. To accommodate this, the myocardium must undergo significant structural rearrangement during contraction. At the mesoscale, collections of cardiomyocytes are organized into sheetlets, and sheetlet shear is the fundamental mechanism of rearrangement that produces wall thickening. Herein, we review the histological and physiological studies of myocardial mesostructure that have established the sheetlet shear model of wall thickening. Recent developments in tissue clearing techniques allow for imaging of whole hearts at the cellular scale, whereas magnetic resonance imaging (MRI) and computed tomography (CT) can image the myocardium at the mesoscale (100 µm to 1 mm) to resolve cardiomyocyte orientation and organization. Through histology, cardiac diffusion tensor imaging (DTI), and other modalities, mesostructural sheetlets have been confirmed in both animal and human hearts. Recent in vivo cardiac DTI methods have measured reorientation of sheetlets during the cardiac cycle. We also examine the role of pathological cardiac remodeling on sheetlet organization and reorientation, and the impact this has on ventricular function and dysfunction. We also review the unresolved mesostructural questions and challenges that may direct future work in the field.

Spatial Coherence in Medical Ultrasound: A Review

Traditional pulse-echo ultrasound imaging heavily relies on the discernment of signals based on their relative magnitudes but is limited in its ability to mitigate sources of image degradation, the most prevalent of which is acoustic clutter. Advances in computing power and data storage have made it possible for echo data to be alternatively analyzed through the lens of spatial coherence, a measure of the similarity of these signals received across an array. Spatial coherence is not currently explicitly calculated on diagnostic ultrasound scanners but a large number of studies indicate that it can be employed to describe image quality, to adaptively select system parameters and to improve imaging and target detection. With the additional insights provided by spatial coherence, it is poised to play a significant role in the future of medical ultrasound. This review details the theory of spatial coherence in pulse-echo ultrasound and key advances made over the last few decades since its introduction in the 1980s.

In vivo whole brain microvascular imaging in mice using transcranial 3D Ultrasound Localization Microscopy

Background

Non-invasive high-resolution imaging of the cerebral vascular anatomy and function is key for the study of intracranial aneurysms, stenosis, arteriovenous malformations, and stroke, but also neurological pathologies, such as degenerative diseases. Direct visualization of the microvascular networks in the whole brain remains however challenging in vivo.

Methods

In this work, we performed 3D ultrafast ultrasound localization microscopy (ULM) using a 2D ultrasound matrix array and mapped the whole-brain microvasculature and flow at microscopic resolution in C57Bl6 mice in vivo.

Findings

We demonstrated that the mouse brain vasculature can be imaged directly through the intact skull at a spatial resolution of 20 µm and over the whole brain depth and at high temporal resolution (750 volumes.s⁻¹). Individual microbubbles were tracked to estimate the flow velocities that ranged from 2 mm.s⁻¹ in arterioles and venules up to 100 mm.s⁻¹ in large vessels. The vascular maps were registered automatically with the Allen atlas in order to extract quantitative vascular parameters such as local flow rates and velocities in regions of interest.

Interpretation

We show the potential of 3D ULM to provide new insights into whole-brain vascular flow in mice models at unprecedented vascular scale for an in vivo technique. This technology is highly translational and has the potential to become a major tool for the clinical investigation of the cerebral microcirculation.

Funding

This study was supported by the European Research Council under the European Union's Seventh Framework Program (FP/2007-2013) / ERC Grant Agreement n° 311025 and by the Fondation Bettencourt-Schueller under the program “Physics for Medicine”. We acknowledge the ART (Technological Research Accelerator) biomedical ultrasound program of INSERM.

Impact of Intraventricular Septal Fiber Orientation on Cardiac Electromechanical Function

Cardiac fiber direction is an important factor determining the propagation of electrical activity, as well as the development of mechanical force. In this article, we imaged the ventricles of several species with special attention to the intraventricular septum to determine the functional consequences of septal fiber organization. First, we identified a dual-layer organization of the fiber orientation in the intraventricular septum of ex vivo sheep hearts using diffusion tensor imaging at high field MRI. To expand the scope of the results, we investigated the presence of a similar fiber organization in five mammalian species (rat, canine, pig, sheep, and human) and highlighted the continuity of the layer with the moderator band in large mammalian species. We implemented the measured septal fiber fields in three-dimensional electromechanical computer models to assess the impact of the fiber orientation. The downward fibers produced a diamond activation pattern superficially in the right ventricle. Electromechanically, there was very little change in pressure volume loops although the stress distribution was altered. In conclusion, we clarified that the right ventricular septum has a downwardly directed superficial layer in larger mammalian species, which can have modest effects on stress distribution.

NEW & NOTEWORTHY A dual-layer organization of the fiber orientation in the intraventricular septum was identified in ex vivo hearts of large mammals. The RV septum has a downwardly directed superficial layer that is continuous with the moderator band. Electrically, it produced a diamond activation pattern. Electromechanically, little change in pressure volume loops were noticed but stress distribution was altered. Fiber distribution derived from diffusion tensor imaging should be considered for an accurate strain and stress analysis.

Experimental Study of Aperiodic Plane Wave Imaging for Ultrafast 3-D Ultrasound Imaging

OBJECTIVE

Although plane wave imaging (PWI) with multiple plane waves (PWs) steered at different angles enables ultrafast three-dimensional (3-D) ultrasonic imaging, there is still a challenging tradeoff between image quality and frame rate. To address this challenge, we recently proposed the aperiodic PWI (APWI) with mathematical analysis and simulation study. In this paper, we demonstrate the feasibility of APWI and evaluate the performance with phantom and in vivo experiments.

METHODS

APWI with a concentric ring angle pattern (APWI-C) and APWI with a sunflower pattern (APWI-S) are evaluated. For experimental verification of the methods, the experimental results are compared with simulation results in terms of the mainlobe-to-sidelobe ratio. In addition, the performance of APWI is compared with that of conventional PWI by using a commercial phantom. To examine the potential for clinical use of APWI, a gallstone-mimicking phantom study and an in vivo carotid artery experiment are also conducted.

RESULTS

In the phantom study, the APWI methods provide a contrast ratio approximately 23 dB higher than that of PWI. In a gallstone mimicking experiment, the proposed methods yield 3-D rendered stone images more similar to the real stones than PWI. In the in vivo carotid artery images, APWI reduces the clutter artifacts inside the artery.

CONCLUSION

Phantom and in vivo studies show that the APWI enhances the contrast without compromising the spatial resolution and frame rate.

SIGNIFICANCE

This study experimentally demonstrates the feasibility and advantage of APWI for ultrafast 3-D ultrasonic imaging.

Ex Vivo Evaluation of Mechanical Anisotropic Tissues with High-Frequency Ultrasound Shear Wave Elastography

The use of imaging devices to assess directional mechanics of tissues is highly desirable. This is because the directional mechanics depend on fiber orientation, and altered directional mechanics are closely related to the pathological status of tissues. However, measuring directional mechanics in tissues with high-stiffness is challenging due to the difficulty of generating localized displacement in these tissues using acoustic radiation force, a general method for generating displacement in ultrasound-based elastography. In addition, common ultrasound probes do not provide rotational function, which makes the measurement of directional mechanics inaccurate and unreliable. Therefore, we developed a high-frequency ultrasound mechanical wave elastography system that can accommodate a wide range of tissue stiffness and is also equipped with a motorized rotation stage for precise imaging of directional mechanics. A mechanical shaker was applied to the elastography system to measure tissues with high-stiffness. Phantom and ex vivo experiments were performed. In the phantom experiments, the lateral and axial resolution of the system were determined to be 144 μm and 168 μm, respectively. In the ex vivo experiments, we used swine heart and cartilage, both of which are considered stiff. The elastography system allows us to acquire the directional mechanics with high angular resolution in the heart and cartilage. The results demonstrate that the developed elastography system is capable of imaging a wide range of tissues and has high angular resolution. Therefore, this system might be useful for the diagnostics of mechanically anisotropic tissues via ex vivo tests.

Optical Ultrastructure of Large Mammalian Hearts Recovers Discordant Alternans by In Silico Data Assimilation

Understanding and predicting the mechanisms promoting the onset and sustainability of cardiac arrhythmias represent a primary concern in the scientific and medical communities still today. Despite the long-lasting effort in clinical and physico-mathematical research, a critical aspect to be fully characterized and unveiled is represented by spatiotemporal alternans patterns of cardiac excitation. The identification of discordant alternans and higher-order alternating rhythms by advanced data analyses as well as their prediction by reliable mathematical models represents a major avenue of research for a broad and multidisciplinary scientific community. Current limitations concern two primary aspects: 1) robust and general-purpose feature extraction techniques and 2) in silico data assimilation within reliable and predictive mathematical models. Here, we address both aspects. At first, we extend our previous works on Fourier transformation imaging (FFI), applying the technique to whole-ventricle fluorescence optical mapping. Overall, we identify complex spatial patterns of voltage alternans and characterize higher-order rhythms by a frequency-series analysis. Then, we integrate the optical ultrastructure obtained by FFI analysis within a fine-tuned electrophysiological mathematical model of the cardiac action potential. We build up a novel data assimilation procedure demonstrating its reliability in reproducing complex alternans patterns in two-dimensional computational domains. Finally, we prove that the FFI approach applied to both experimental and simulated signals recovers the same information, thus closing the loop between the experiment, data analysis, and numerical simulations.

Architected fibrous scaffolds for engineering anisotropic tissues

Mimicking the native three-dimensional microenvironment is of crucial importance when biofabricating a new healthcare material. One aspect of the native tissue that is often omitted when designing a suitable scaffold is its anisotropy. Not only is matching native mechanical properties important when designing implantable scaffolds or healthcare materials, but matching physiological structure is also important as many cell populations respond differently to fiber orientation. Therefore, novel aligned electrospun scaffolds with varying fiber angles and spacing of bundles were created and mechanically characterized. Through controlling the angle between the fibers in each layer of the scaffold, a range of different physiological anisotropic mechanical properties were achieved that encompasses values found in native tissues. Extrapolation of this mechanical data allowed for any native tissue's anisotropic Young's modulus to be mimicked by electrospinning fibers at a particular angle. These electrospun scaffolds were then incorporated with cell-laden hydrogels to create hybrid structures that contain the benefits of both scaffolding techniques with the ability to encapsulate cells in the hydrogel. To conclude, this study develops a novel bundled fiber scaffold that was architected to yield anisotropic properties matching native tissues.

Cardiac mechanostructure: Using mechanics and anisotropy as inspiration for developing epicardial therapies in treating myocardial infarction

The mechanical environment and anisotropic structure of the heart modulate cardiac function at the cellular, tissue and organ levels. During myocardial infarction (MI) and subsequent healing, however, this landscape changes significantly. In order to engineer cardiac biomaterials with the appropriate properties to enhance function after MI, the changes in the myocardium induced by MI must be clearly identified. In this review, we focus on the mechanical and structural properties of the healthy and infarcted myocardium in order to gain insight about the environment in which biomaterial-based cardiac therapies are expected to perform and the functional deficiencies caused by MI that the therapy must address. From this understanding, we discuss epicardial therapies for MI inspired by the mechanics and anisotropy of the heart focusing on passive devices, which feature a biomaterials approach, and active devices, which feature robotic and cellular components. Through this review, a detailed analysis is provided in order to inspire further development and translation of epicardial therapies for MI.

Architecture for an Ultrasound Advanced Open Platform With an Arbitrary Number of Independent Channels

Ultrasound open platforms are programmable and flexible tools for the development and test of novel methods. In most cases, they embed the electronics for the independent control of (maximum) 256 probe elements. However, a higher number of channels is needed for the control of 2-D array probes. This paper presents a system architecture that, through the hardware and software synchronization of multiple ULA-OP 256 scanners, may implement advanced open platforms with an arbitrary number of channels. The proposed solution needs a single personal computer, maintains real-time features, and preserves portability. A prototype demonstrator, composed of two ULA-OP 256 scanners connected to 512 elements of a matrix array, was implemented and tested according to different channel configurations. Experiments performed under MATLAB control confirmed that by doubling the number of elements (from 256 to 512) the signal-to-noise and contrast ratios improve by 9 dB and 3 dB, respectively. Furthermore, as a full 512-channel scanner, the demonstrator can produce real-time B-mode images at 18 Hz, high enough for probe positioning during acquisitions. Also, the demonstrator permitted the implementation of a new high frame rate, bi-plane, triplex modality. All probe elements are excited to simultaneously produce two planar, perpendicular diverging waves. Each scanner independently processes the echoes received by the 256 connected elements to beamform 1300 frames per second. For each insonified plane, good quality morphological (B-mode), qualitative (color flow-), and quantitative (spectral-) Doppler images are finally shown in real-time by a dedicated interface.

GPU-Based Simulation of Echocardiography Volumes Using Quantitative Fiber-Angle-to-Backscatter Measurements

The intensity of backscattered ultrasound signal from heart muscle is known to be related to the angle between cardiac fibers and the insonification direction. In this work, a GPU-based method of simulating three-dimensional (3D) echocardiographic images from an empirically derived angle-to-backscatter relationship is developed and validated. Images of a rotating fiber phantom are simulated, and it is validated that the angle-to-backscatter relationship is accurately reflected by the simulated envelope data. In a second experiment, echocardiography images are simulated from a diffusion tensor magnetic resonance imaging (DT-MRI) volume of a canine heart to demonstrate that the method produces view-dependent speckle. 3D volumes of a parametrically generated ideal left ventricle phantom are also simulated and processed into fiber orientation maps using the underlying quantitative parameters. Images are simulated based on the characteristics of a 35-by-32 two-dimensional (2D) matrix array probe and a clinical one-dimensional (1D) phased array probe. The processed fiber volumes exhibit good agreement with the virtual phantom's ground truth, having an average acute angle error (AAE) of less than 10 degrees for both probes. The simulation method is fast and opens a new approach for ultrasound fiber imaging.

Local scattering ultrasound imaging

Ultrasonic imaging is a widely used tool for detection, localisation and characterisation of material inhomogeneities with important applications in many fields. This task is particularly challenging when imaging in a complex medium, where the ultrasonic wave is scattered by the material microstructure, preventing detection and characterisation of weak targets. Fundamentally, the maximum information that can be experimentally obtained from each material region consists of a set of reflected signals for different incident waves. However, these data are not directly accessible from the raw measurements, which represent a superposition of reflections from all scatterers in the medium. Here we show, that a complete set of transmitter–receiver data encodes sufficient information in order to achieve full spatio–temporal separation of transmitter–receiver data, corresponding to different local scattering areas. We show that access to the local scattering data can provide valuable benefits for many applications. More importantly, this technique enables fundamentally new approaches, exploiting the angular distribution of the scattering amplitude and phase of each local scattering region. Here we demonstrate how the local scattering directivity can be used to build the local scattering image, releasing the full potential and richness of the transmit–receive data. As a proof of concept, we demonstrate the detection of small inclusions in various highly scattering materials using numerical and experimental examples. The described principles are very general and can be applied to any research field where the phased array technology is employed.

Mechanical stimuli for left ventricular growth during pressure overload

BACKGROUND

The mechanical stimulus (i.e. stress or stretch) for growth occurring in the pressure-overloaded left ventricle (LV) is not exactly known.

OBJECTIVE

To address this issue, we investigate the correlation between local ventricular growth (indexed by local wall thickness) and the local acute changes in mechanical stimuli after aortic banding.

METHODS

LV geometric data were extracted from 3D echo measurements at baseline and 2 weeks in the aortic banding swine model (n = 4). We developed and calibrated animal-specific finite element (FE) model of LV mechanics against pressure and volume waveforms measured at baseline. After the simulation of the acute effects of pressure-overload, the local changes of maximum, mean and minimum myocardial stretches and stresses in three orthogonal material directions (i.e., fiber, sheet and sheet-normal) over a cardiac cycle were quantified. Correlation between mechanical quantities and the corresponding measured local changes in wall thickness was quantified using the Pearson correlation number (PCN) and Spearman rank correlation number (SCN).

RESULTS

At 2 weeks after banding, the average septum thickness decreased from 10.6 ± 2.92mm to 9.49 ± 2.02mm, whereas the LV free-wall thickness increased from 8.69 ± 1.64mm to 9.4 ± 1.22mm. The FE results show strong correlation of growth with the changes in maximum fiber stress (PCN = 0.5471, SCN = 0.5111) and changes in the mean sheet-normal stress (PCN= 0.5266, SCN = 0.5256). Myocardial stretches, however, do not have good correlation with growth.

CONCLUSION

These results suggest that fiber stress is the mechanical stimuli for LV growth in pressure-overload.

4D Ultrafast Ultrasound Imaging of Naturally Occurring Shear Waves in the Human Heart

The objectives were to develop a novel three-dimensional technology for imaging naturally occurring shear wave (SW) propagation, demonstrate feasibility on human volunteers and quantify SW velocity in different propagation directions. Imaging of natural SWs generated by valve closures has emerged to obtain a direct measurement of cardiac stiffness. Recently, natural SW velocity was assessed in two dimensions on parasternal long axis view under the assumption of a propagation direction along the septum. However, in this approach the source localization and the complex three-dimensional propagation wave path was neglected making the speed estimation unreliable. High volume rate transthoracic acquisitions of the human left ventricle (1100 volume/s) was performed with a 4D ultrafast echocardiographic scanner. Four-dimensional tissue velocity cineloops enabled visualization of aortic and mitral valve closure waves. Energy and time of flight mapping allowed propagation path visualization and source localization, respectively. Velocities were quantified along different directions. Aortic and mitral valve closure SW velocities were assessed for the three volunteers with low standard deviation. Anisotropic propagation was also found suggesting the necessity of using a three-dimensional imaging approach. Different velocities were estimated for the three directions for the aortic (3.4± 0.1 m/s, 3.5± 0.3 m/s, 5.4± 0.7 m/s) and the mitral (2.8± 0.5 m/s, 2.9± 0.3 m/s, 4.6± 0.7 m/s) valve SWs. 4D ultrafast ultrasound alleviates the limitations of 2D ultrafast ultrasound for cardiac SW imaging based on natural SW propagations and enables a comprehensive measurement of cardiac stiffness. This technique could provide stiffness mapping of the left ventricle.

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

Cardiac multi-scale investigation of the right and left ventricle ex vivo: a review

The heart is a complex multi-scale system composed of components integrated at the subcellular, cellular, tissue and organ levels. The myocytes, the contractile elements of the heart, form a complex three-dimensional (3D) network which enables propagation of the electrical signal that triggers the contraction to efficiently pump blood towards the whole body. Cardiovascular diseases (CVDs), a major cause of mortality in developed countries, often lead to cardiovascular remodeling affecting cardiac structure and function at all scales, from myocytes and their surrounding collagen matrix to the 3D organization of the whole heart. As yet, there is no consensus as to how the myocytes are arranged and packed within their connective tissue matrix, nor how best to image them at multiple scales. Cardiovascular imaging is routinely used to investigate cardiac structure and function as well as for the evaluation of cardiac remodeling in CVDs. For a complete understanding of the relationship between structural remodeling and cardiac dysfunction in CVDs, multi-scale imaging approaches are necessary to achieve a detailed description of ventricular architecture along with cardiac function. In this context, ventricular architecture has been extensively studied using a wide variety of imaging techniques: ultrasound (US), optical coherence tomography (OCT), microscopy (confocal, episcopic, light sheet, polarized light), magnetic resonance imaging (MRI), micro-computed tomography (micro-CT) and, more recently, synchrotron X-ray phase contrast imaging (SR X-PCI). Each of these techniques have their own set of strengths and weaknesses, relating to sample size, preparation, resolution, 2D/3D capabilities, use of contrast agents and possibility of performing together with in vivo studies. Therefore, the combination of different imaging techniques to investigate the same sample, thus taking advantage of the strengths of each method, could help us to extract the maximum information about ventricular architecture and function. In this review, we provide an overview of available and emerging cardiovascular imaging techniques for assessing myocardial architecture ex vivo and discuss their utility in being able to quantify cardiac remodeling, in CVDs, from myocyte to whole organ.

Association between Sarcopenia/Lower Muscle Mass and Short-Term Regression of Deep Vein Thrombosis Using Direct Oral Anticoagulants

Advanced age, obesity, and muscle weakness are independent factors in the onset of deep vein thrombosis (DVT). Recently, an association between sarcopenia and DVT has been reported. We hypothesized that sarcopenia related factors, observed by ultrasonography, are associated with the regression effect on the thrombus following anticoagulation therapy. The present study focused on gastrocnemius muscle (GCM) thickness and the GCM's internal echogenic brightness. We examined the association with DVT regression following direct oral anticoagulants (DOACs) treatment.The prospective cohort study period was between October 2017 and August 2018. We enrolled 46 patients diagnosed with DVT by ultrasonography, who were aged >60 years old and treated with DOACs. Sarcopenia was evaluated using the Asian Working Group for Sarcopenia flowchart. The average DOACs treatment period was 94 days, and 29 patients exhibited thrombus regression. On univariate logistic regression analysis, sarcopenia, average GCM diameter index, and gastrocnemius integrated backscatter index were significantly associated with thrombus regression. In a multivariate model, only the average GCM diameter index correlated with thrombus regression.The average GCM diameter index is associated with DVT regression treated with DOACs. Considering the GCM diameter during DVT treatment can be a marker to make a decision for the treatment of DVT.

High Frame Rate Ultrasound for Electromechanical Wave Imaging to Differentiate Endocardial From Epicardial Myocardial Activation

Differentiation between epicardial and endocardial ventricular activation remains a challenge despite the latest technologies available. The aim of the present study was to develop a new tool method, based on electromechanical wave imaging (EWI), to improve arrhythmogenic substrate activation analysis. Experiments were conducted on left ventricles (LVs) of four isolated working mode swine hearts. The protocol aimed at demonstrating that different patterns of mechanical activation could be observed whether the ventricle was in sinus rhythm, paced from the epicardium or from the endocardium. A total of 72 EWI acquisitions were recorded on the anterior, lateral and posterior segments of the LV. A total of 54 loop records were blindly assigned to two readers. EWI sequences interpretations were correct in 89% of cases. The overall agreement rate between the two readers was 83%. When in a paced ventricle, the origin of the wave front was focal and originated from the endocardium or the epicardium. In sinus rhythm, wave front was global and activated within the entire endocardium toward the epicardium at a speed of 1.7 ± 0.28 m·s^-1. Wave front speeds were respectively measured when the endocardium or the epicardium were paced at a speed of 1.1 ± 0.35 m·s^-1 versus 1.3 ± 0.34 m·s^-1 (p = NS). EWI activation mapping allows activation localization within the LV wall and calculation of the wave front propagation speed through the muscle. In the future, this technology could help localize activation within the LV thickness during complex ablation procedures.

Myocardial Fiber Mapping of Rat Hearts, Using Apparent Backscatter, with Histologic Validation

Myocardial fiber architecture is a physiologically important regulator of ejection fraction, strain and pressure development. Apparent ultrasonic backscatter has been shown to be a useful method for recreating the myocardial fiber architecture in human-sized sheep hearts because of the dependence of its amplitude on the relative orientation of a myofiber to the angle of ultrasonic insonification. Thus, the anisotropy of the backscatter signal is linked to and provides information about the fiber orientation. In this study, we sought to determine whether apparent backscatter could be used to measure myofiber orientation in rodent hearts. Fixed adult-rat hearts were imaged intact, and both a transmural cylindrical core and transmural wedge of the left ventricular free wall were imaged. Cylindrical core samples confirmed that backscatter anisotropy could be measured in rat hearts. Ultrasound and histologic analysis of transmural myocardial wedge samples confirmed that the apparent backscatter could be reproducibly mapped to fiber orientation (angle of the fiber relative to the direction of insonification). These data provided a quantitative relationship between the apparent backscatter and fiber angle that was applied to whole-heart images. Myocardial fiber architecture was successfully measured in rat hearts. Quantifying myocardial fiber architecture, using apparent backscatter, provides a number of advantages, including its scalable use from rodents to man, its rapid low-cost acquisition and minimal contraindications. The method outlined in this study provides a method for investigators to begin detailed assessments of how the myocardial fiber architecture changes in preclinical disease models, which can be immediately translated into the clinic.

4D simultaneous tissue and blood flow Doppler imaging: revisiting cardiac Doppler index with single heart beat 4D ultrafast echocardiography

The goal of this study was to demonstrate the feasibility of semi-automatic evaluation of cardiac Doppler indices in a single heartbeat in human hearts by performing 4D ultrafast echocardiography with a dedicated sequence of 4D simultaneous tissue and blood flow Doppler imaging. 4D echocardiography has the potential to improve the quantification of major cardiac indices by providing more reproducible and less user dependent measurements such as the quantification of left ventricle (LV) volume. The evaluation of Doppler indices, however, did not benefit yet from 4D echocardiography because of limited volume rates achieved in conventional volumetric color Doppler imaging but also because spectral Doppler estimation is still restricted to a single location. High volume rate (5200 volume s^-1) transthoracic simultaneous tissue and blood flow Doppler acquisitions of three human LV were performed using a 4D ultrafast echocardiography scanner prototype during a single heartbeat. 4D color flow, 4D tissue Doppler cineloops and spectral Doppler at each voxel were computed. LV outflow tract, mitral inflow and basal inferoseptal locations were automatically detected. Doppler indices were derived at these locations and were compared against clinical 2D echocardiography. Blood flow Doppler indices E (early filling), A (atrial filling), E/A ratio, S (systolic ejection) and cardiac output were assessed on the three volunteers. Simultaneous tissue Doppler indices e' (mitral annular velocity peak), a' (late velocity peak), e'/a' ratio, s' (systolic annular velocity peak), E/e' ratio were also estimated. Standard deviations on three independent acquisitions were averaged over the indices and was found to be inferior to 4% and 8.5% for Doppler flow and tissue Doppler indices, respectively. Comparison against clinical 2D echocardiography gave a p  value larger than 0.05 in average indicating no significant differences. 4D ultrafast echocardiography can quantify the major cardiac Doppler indices in a single heart beat acquisition.

Structural and Biomechanical Adaptations of Right Ventricular Remodeling - in Pulmonary Arterial Hypertension - Reduces Left Ventricular Rotation During Contraction: A Computational Study

Pulmonary hypertension (PH) is a degenerative disease characterized by progressively increased right ventricular (RV) afterload that leads to ultimate functional decline [1]. Recent observational studies have documented a decrease in left ventricular (LV) torsion during ejection, with preserved LV ejection fraction (EF) in pediatric and adult PH patients [2-4]. The objective of this study was to develop a computational model of the bi-ventricular heart and use it to evaluate changes in LV torsion mechanics in response to mechanical, structural, and hemodynamic changes in the RV free-wall. The heart model revealed that LV apex rotation and torsion were decreased when increasing RV mechanical rigidity and during re-orientation of RV myocardial fibers. Furthermore, structural changes to the RV appear to have a notable impact on RV EF, but little influence on LV EF. Finally, RV pressure overload exponentially increased LV myocardial stress. The computational results found in this study are consistent with clinical observations in adult and pediatric PH patients, which reveal a decrease in LV torsion with preserved LV EF [3, 4]. Furthermore, discovered causes of decreased LV torsion are consistent with RV structural adaptations seen in PH rodent studies [5], which might also explain suspected stress-induced changes in LV myocardial gene/protein expression.

Experimental 3-D Ultrasound Imaging with 2-D Sparse Arrays using Focused and Diverging Waves

Three dimensional ultrasound (3-D US) imaging methods based on 2-D array probes are increasingly investigated. However, the experimental test of new 3-D US approaches is contrasted by the need of controlling very large numbers of probe elements. Although this problem may be overcome by the use of 2-D sparse arrays, just a few experimental results have so far corroborated the validity of this approach. In this paper, we experimentally compare the performance of a fully wired 1024-element (32 × 32) array, assumed as reference, to that of a 256-element random and of an “optimized” 2-D sparse array, in both focused and compounded diverging wave (DW) transmission modes. The experimental results in 3-D focused mode show that the resolution and contrast produced by the optimized sparse array are close to those of the full array while using 25% of elements. Furthermore, the experimental results in 3-D DW mode and 3-D focused mode are also compared for the first time and they show that both the contrast and the resolution performance are higher when using the 3-D DW at volume rates up to 90/second which represent a 36x speed up factor compared to the focused mode.

Effect of microstructural elongation on backscattered field: Intensity measurement and multiple scattering estimation with a linear transducer array

The effect of microstructural elongation on ultrasonic backscattered fields was studied. Two methods for determining the elongation direction of macrozones in titanium alloys, using the anisotropic spatial coherence of the backscattered field, are presented. Both methods use a phased array attached on a rotative holder that records the array response matrix at several angles. Two titanium alloys were investigated: TA6V and Ti17. TA6V exhibited a strong macrozone elongation, whereas Ti17 macrozones were found equiaxial. The first method is based on the measurement of backscattered intensity in function of the probe angle relative to the macrozones elongation direction. An angular dependence of backscattered intensity is observed in presence of elongated scatterers, and their elongation direction is collinear with the probe direction corresponding to a minimal intensity. This variability is linked to both piezoelectric shape and the backscattered field spatial properties. The second method is based on the measurement of the relative proportion of single to multiple scattering in a diffusive media, using a simplified version of the single scattering filter developed in Aubry and Derode (2009). It allows the measurement of the level of multiple scattering: both titanium alloys exhibited strong multiple scattering. The elongation direction was determined as the direction of minimal multiple scattering. Furthermore, these results were confirmed by the measurement of the coherent backscattering cone on both samples.