Super-Resolution Imaging Through the Human Skull

Boiling Histotripsy in Ex Vivo Human Brain: Proof-of-concept

OBJECTIVE

Non-invasive surgical approaches, including boiling histotripsy (BH), are currently being developed for the treatment of brain disorders aiming to avoid craniotomy and exposure of intervening tissues, and, thus, minimize associated complications. This work aimed to demonstrate the feasibility of BH for mechanical fractionation of human brain tissues ex vivo under B-mode guidance, with preliminary measurements of tissue stiffness via shear wave elastography.

METHODS

Young's moduli of 25 human autopsy brain samples obtained from de-identified patients of 51-91 y old (median 77 y old) were measured via shear wave elastography prior to BH sonication. Seventeen volumetric BH lesions (1-4 layers of 5 × 5 points with a 1- mm step) were produced near brain surface (n = 10), in white matter (n = 3), in thalamus (n = 2), and globus pallidus (n = 2) using 12 element 1.5 MHz sector transducer under B-mode guidance with 10 ms or 2 ms pulses delivered 10 or 15 times per sonication point with 1% duty cycle. After treatment, the lesions were evaluated grossly through bisection, histologically with hematoxylin and eosin staining, and ultrastructurally via scanning and transmission electron microscopy.

RESULTS

Young's moduli of autopsy brain samples were lower in older patients (from 32.9 ± 6.6 kPa in 51 y olds to 10 ± 2 kPa in 91 y olds) and at higher temperature (6%-50% lower at 37°С vs 23°С), and were within the range observed clinically. All tested BH treatments performed near the brain surface (i.e., mostly in gray matter) resulted in formation of well-demarcated rectangular lesions with homogenized content and sharp boundaries, with majority of residual fragments below 100 microns. The use of shorter pulses (2 ms vs 10 ms) accelerated the treatment at least threefold, and the highest liquefaction rate was 568 mm3/min. White matter was more resistant to BH vs gray matter: at least 15 pulses of 2 ms duration were required per each sonication point, and the liquefaction rate was three times lower. The ability of BH to produce lesions in thalamus and globus pallidus was also confirmed.

CONCLUSION

This work presents the first demonstration of BH proof-of-concept in human brain tissues ex vivo under B-mode guidance with clinically relevant treatment rates.

Transcranial ultrafast ultrasound Doppler imaging: A phantom study

Ultrafast ultrasound Doppler imaging facilitates the assessment of cerebral hemodynamics with high spatio-temporal resolution. However, the significant acoustic impedance mismatch between the skull and soft tissue results in phase aberrations, which can compromise the quality of transcranial imaging and introduce biases in velocity and direction quantification of blood flow. This paper proposed an aberration correction method that combines deep learning-based skull sound speed modelling with ray theory to realize transcranial plane-wave imaging and ultrafast Doppler imaging. The method was validated through phantom experiments using a linear array with a center frequency of 6.25 MHz, 128 elements, and a pitch of 0.3 mm. The results demonstrated an improvement in the imaging quality of intracranial targets when using the proposed method. After aberration correction, the average locating deviation decreased from 1.40 mm to 0.27 mm in the axial direction, from 0.50 mm to 0.20 mm in the lateral direction, and the average full-width-at-half-maximum (FWHM) decreased from 1.37 mm to 0.97 mm for point scatterers. For circular inclusions, the average contrast-to-noise ratio (CNR) improved from 8.1 dB to 11.0 dB, and the average eccentricity decreased from 0.36 to 0.26. Furthermore, the proposed method was applied to transcranial ultrafast Doppler flow imaging. The results showed a significant improvement in accuracy and quality for blood Doppler flow imaging. The results in the absence of the skull were considered as the reference, and the average normalized root-mean-square errors of the axial velocity component on the five selected axial profiles were reduced from 17.67% to 8.02% after the correction.

A transcranial multiple waves suppression method for plane wave imaging based on Radon transform

Transcranial ultrasound imaging presents a significant challenge due to the intricate interplay between ultrasound waves and the heterogeneous human skull. The skull's presence induces distortion, refraction, multiple scattering, and reflection of ultrasound signals, thereby complicating the acquisition of high-quality images. Extracting reflections from the entire waveform is crucial yet exceedingly challenging, as intracranial reflections are often obscured by strong amplitude direct waves and multiple scattering. In this paper, a multiple wave suppression method for ultrasound plane wave imaging is proposed to mitigate the impact of skull interference. Drawing upon prior research, we developed an enhanced high-resolution linear Radon transform using the maximum entropy principle and Bayesian method, facilitating wavefield separation. We detailed the process of wave field separation in the Radon domain through simulation of a model with a high velocity layer. When plane waves emitted at any steering angles, both multiple waves and first arrival waves manifested as distinct energy points. In the brain simulation, we contrasted the characteristic differences between skull reflection and brain-internal signal in Radon domain, and demonstrated that multiples suppression method reduces side and grating lobe levels by approximately 30 dB. Finally, we executed in vitro experiments using a monkey skull to separate weak intracranial reflection signals from strong skull reflections, enhancing the contrast-to-noise ratio by 85 % compared to conventional method using full waveform. This study deeply explores the effect of multiples on effective signal separation, addresses the complexity of wavefield separation, and verifies its efficacy through imaging, thereby significantly advancing ultrasound transcranial imaging techniques.

Element Position Calibration for Matrix Array Transducers with Multiple Disjoint Piezoelectric Panels

Two-dimensional ultrasound transducers enable the acquisition of fully volumetric data that have been demonstrated to provide greater diagnostic information in the clinical setting and are a critical tool for emerging ultrasound methods, such as super-resolution and functional imaging. This technology, however, is not without its limitations. Due to increased fabrication complexity, some matrix probes with disjoint piezoelectric panels may require initial calibration. In this manuscript, two methods for calibrating the element positions of the Vermon 1024-channel 8 MHz matrix transducer are detailed. This calibration is a necessary step for acquiring high resolution B-mode images while minimizing transducer-based image degradation. This calibration is also necessary for eliminating vessel-doubling artifacts in super-resolution images and increasing the overall signal-to-noise ratio (SNR) of the image. Here, we show that the shape of the point spread function (PSF) can be significantly improved and PSF-doubling artifacts can be reduced by up to 10 dB via this simple calibration procedure.

Proof of Concept of an Affordable, Compact and Transcranial Submillimeter Accurate Ultrasound-Based Tracking System

In neurosurgery, a current challenge is to provide localized therapy in deep and difficult-to-access brain areas with millimeter accuracy. In this prospect, new surgical devices such as microrobots are being developed, which require controlled inbrain navigation to ensure the safety and efficiency of the intervention. In this context, the device tracking technology has to answer a three-sided challenge: invasiveness, performance, and facility of use. Although ultrasound seems appropriate for transcranial tracking, the skull remains an obstacle because of its significant acoustic perturbations. A compact and affordable ultrasound-based tracking system that minimizes skull-related disturbances is proposed here. This system consists of three emitters fixed on the patient's head and a one-millimeter receiver embedded in the surgical device. The 3D position of the receiver is obtained by trilateration based on time of flight measurements. The system demonstrates a submillimeter tracking accuracy through an 8.9 mm thick skull plate phantom. This result opens multiple perspectives in terms of millimeter accurate navigation for a large number of neurobiomedical devices.

Controlled noninvasive modulation of deep brain regions in humans

Transcranial focused ultrasound provides noninvasive and reversible approaches for precise and personalized manipulations of brain circuits, with the potential to transform our understanding of brain function and treatments of brain dysfunction. However, effective applications in humans have been limited by the human head, which attenuates and distorts ultrasound severely and unpredictably. This has led to uncertain ultrasound intensities delivered into the brain. Here, we address this lingering barrier using a direct measurement approach that can be repeatedly applied to the human brain. The approach uses an ultrasonic scan of the head to measure and compensate for the attenuation of the ultrasound by all obstacles within the ultrasound path. No other imaging modality is required and the method is parameter-free and personalized to each subject. The approach accurately restores operators’ intended intensities inside ex-vivo human skulls. Moreover, the approach is critical for effective modulation of deep brain regions in humans. When applied, the approach modulates fMRI Blood Oxygen Level Dependent (BOLD) activity in disease-relevant deep brain regions. This tool unlocks the potential of emerging approaches based on low-intensity ultrasound for controlled manipulations of neural circuits in humans.

Transcranial focused ultrasound has had limited applications in humans due to the unpredictable distortions of ultrasound by the human head. Thomas Riis and colleagues report an approach which enables direct correction for the attenuation of ultrasound by the skull and hair, thus enabling controlled ultrasound therapies in humans.

Ray theory-based compounded plane wave ultrasound imaging for aberration corrected transcranial imaging: Phantom experiments and simulations

Compounded plane wave imaging (CPWI) allows high-frame-rate measurement and has been one of the most promising modalities for real-time brain imaging. However, ultrasonic brain imaging using the CPWI modality is usually performed with a worn thin or removal of the skull layer. Otherwise, the skull layer is expected to distort the ultrasonic wavefronts and significantly decrease intracranial imaging quality. The motivation of this study is to investigate a CPWI method for transcranial brain imaging with the skull layer. A coordinate transformation ray-tracing (CTRT) approach was proposed to track the distorted ultrasonic wavefronts and calculate the time delays for the ultrasound plane wave passing through the skull layer. With an accurate correction for the time delays in beamforming, the CTRT-based CPWI could achieve high-quality intracranial images with the presence of skulls. The proposed CTRT-based CPWI method was verified using a simplified three-layer transcranial model. The full-wave simulation demonstrated that CTRT could accurately (i.e., relative percentage error less than0.18%) track the distorted transmitting wavefront through skull. Compared with the CPWI without aberration correction, the CTRT-based CPWI provided high-quality intracranial imaging and could accurately localize intracranial point scatterers; specifically, positioning error decreases from 0.5 mm to 0.1 mm on average in the axial direction and from 0.7 mm to 0.1 mm on average in the lateral direction. As the compounded angles increased in the CTRT-based CPWI, the contrast improved by 16.2 dB on average for the region of interest, and the array performance indicator (representing resolution) decreased by 4.0 on average for the intracranial point scatterers. The CTRT is of low computational cost compared with full wave simulation. This study suggested that the proposed CTRT-based CPWI might have the potential for real-time and non-invasive transcranial aberration-corrected imaging.

Three-dimensional ultrasound matrix imaging

Matrix imaging paves the way towards a next revolution in wave physics. Based on the response matrix recorded between a set of sensors, it enables an optimized compensation of aberration phenomena and multiple scattering events that usually drastically hinder the focusing process in heterogeneous media. Although it gave rise to spectacular results in optical microscopy or seismic imaging, the success of matrix imaging has been so far relatively limited with ultrasonic waves because wave control is generally only performed with a linear array of transducers. In this paper, we extend ultrasound matrix imaging to a 3D geometry. Switching from a 1D to a 2D probe enables a much sharper estimation of the transmission matrix that links each transducer and each medium voxel. Here, we first present an experimental proof of concept on a tissue-mimicking phantom through ex-vivo tissues and then, show the potential of 3D matrix imaging for transcranial applications.

Aberration correction in abdominal histotripsy

In transabdominal histotripsy, ultrasound pulses are focused on the body to noninvasively destroy soft tissues via cavitation. However, the ability to focus is limited by phase aberration, or decorrelation of the ultrasound pulses due to spatial variation in the speed of sound throughout heterogeneous tissue. Phase aberration shifts, broadens, and weakens the focus, thereby reducing the safety and efficacy of histotripsy therapy. This paper reviews and discusses aberration effects in histotripsy and in related therapeutic ultrasound techniques (e.g., high intensity focused ultrasound), with an emphasis on aberration by soft tissues. Methods for aberration correction are reviewed and can be classified into two groups: model-based methods, which use segmented images of the tissue as input to an acoustic propagation model to predict and compensate phase differences, and signal-based methods, which use a receive-capable therapy array to detect phase differences by sensing acoustic signals backpropagating from the focus. The relative advantages and disadvantages of both groups of methods are discussed. Importantly, model-based methods can correct focal shift, while signal-based methods can restore substantial focal pressure, suggesting that both methods should be combined in a 2-step approach. Aberration correction will be critical to improving histotripsy treatments and expanding the histotripsy treatment envelope to enable non-invasive, non-thermal histotripsy therapy for more patients.

The Emergence of Functional Ultrasound for Noninvasive Brain-Computer Interface

A noninvasive brain-computer interface is a central task in the comprehensive analysis and understanding of the brain and is an important challenge in international brain-science research. Current implanted brain-computer interfaces are cranial and invasive, which considerably limits their applications. The development of new noninvasive reading and writing technologies will advance substantial innovations and breakthroughs in the field of brain-computer interfaces. Here, we review the theory and development of the ultrasound brain functional imaging and its applications. Furthermore, we introduce latest advancements in ultrasound brain modulation and its applications in rodents, primates, and human; its mechanism and closed-loop ultrasound neuromodulation based on electroencephalograph are also presented. Finally, high-frequency acoustic noninvasive brain-computer interface is prospected based on ultrasound super-resolution imaging and acoustic tweezers.

Transcranial Ultrasound Imaging With Decomposition Descent Learning-Based Full Waveform Inversion

Noninvasive brain diagnosis is extremely important because of its efficiency, low cost, and painless nature in the prediction of stroke, cerebral hemorrhage, and other brain research. At present, achieving full 3-D quantitative ultrasonic imaging of the human brain is a cutting-edge challenge due to the complex structures of the human brain and the strong scattering caused by the skulls. In this article, we achieved quantitative ultrasonic imaging of inside-brain anomalies with our proposed method, the decomposition descent learning-based full waveform inversion (DDL-FWI). The proposed method adopts a linear residual decomposing technique to greatly alleviate the computation burden in fast inversion tomography (FIT) with enhanced convergence guaranteed by residual functions. Testing results in both simulation and laboratory experiments demonstrated that our method can achieve high-quality quantitative imaging of brain soft tissues and skulls even starting from homogeneous water background in 2-D, and this method is capable of reconstructing both complex brain tissues and clots in 2-D and 3-D cases using either clean or noisy signals, with a robust 3-D clot resolution as small as 18 mm and 2-D reconstruction speed in 11.20 s. Combined with advanced ultrasonic hardware, DDL-FWI can be easily trained and used for brain imaging efficiently that frees patients from harmful influences from traditional imaging techniques, e.g., ionization radiations from X-ray computed tomography (CT).

Ultrasound super-resolution imaging for differential diagnosis of breast masses

OBJECTIVE

Ultrasound imaging has been widely used in breast cancer screening. Recently, ultrasound super-resolution imaging (SRI) has shown the capability to break the diffraction limit to display microvasculature. However, the application of SRI on differential diagnosis of breast masses remains unknown. Therefore, this study aims to evaluate the feasibility and clinical value of SRI for visualizing microvasculature and differential diagnosis of breast masses.

METHODS

B mode, color-Doppler flow imaging (CDFI) and contrast-enhanced ultrasound (CEUS) images of 46 patients were collected respectively. SRI were generated by localizations of each possible contrast signals. Micro-vessel density (MVD) and microvascular flow rate (MFR) were calculated from SRI and time to peak (TTP), peak intensity (PI) and area under the curve (AUC) were obtained by quantitative analysis of CEUS images respectively. Pathological results were considered as the gold standard. Independent chi-square test and multivariate logistic regression analysis were performed using these parameters to examine the correlation.

RESULTS

The results showed that SRI technique could be successfully applied on breast masses and display microvasculature at a significantly higher resolution than the conventional CDFI and CEUS images. The results showed that the PI, AUC, MVD and MFR of malignant breast masses were significantly higher than those of benign breast masses, while TTP was significantly lower than that of benign breast masses. Among all five parameters, MVD showed the highest positive correlation with the malignancy of breast masses.

CONCLUSIONS

SRI is able to successfully display the microvasculature of breast masses. Compared with CDFI and CEUS, SRI can provide additional morphological and functional information for breast masses. MVD has a great potential in assisting the differential diagnosis of breast masses as an important imaging marker.

Ultrasound super-resolution imaging for the differential diagnosis of thyroid nodules: A pilot study

Objective

Ultrasound imaging provides a fast and safe examination of thyroid nodules. Recently, the introduction of super-resolution imaging technique shows the capability of breaking the Ultrasound diffraction limit in imaging the micro-vessels. The aim of this study was to evaluate its feasibility and value for the differentiation of thyroid nodules.

Methods

In this study, B-mode, contrast-enhanced ultrasound, and color Doppler flow imaging examinations were performed on thyroid nodules in 24 patients. Super-resolution imaging was performed to visualize the microvasculature with finer details. Microvascular flow rate (MFR) and micro-vessel density (MVD) within thyroid nodules were computed. The MFR and MVD were used to differentiate the benign and malignant thyroid nodules with pathological results as a gold standard.

Results

Super-resolution imaging (SRI) technique can be successfully applied on human thyroid nodules to visualize the microvasculature with finer details and obtain the useful clinical information MVD and MFR to help differential diagnosis. The results suggested that the mean value of the MFR within benign thyroid nodule was 16.76 ± 6.82 mm/s whereas that within malignant thyroid was 9.86 ± 4.54 mm/s. The mean value of the MVD within benign thyroid was 0.78 while the value for malignant thyroid region was 0.59. MFR and MVD within the benign thyroid nodules were significantly higher than those within the malignant thyroid nodules respectively (p < 0.01).

Conclusions

This study demonstrates the feasibility of ultrasound super-resolution imaging to show micro-vessels of human thyroid nodules via a clinical ultrasound platform. The important imaging markers, such as MVD and MFR, can be derived from SRI to provide more useful clinical information. It has the potential to be a new tool for aiding differential diagnosis of thyroid nodules.

3D transcranial ultrasound localization microscopy for discrimination between ischemic and hemorrhagic stroke in early phase

Early diagnosis is a critical part of the emergency care of cerebral hemorrhages and ischemia. A rapid and accurate diagnosis of strokes reduces the delays to appropriate treatments and a better functional recovery. Currently, CTscan and MRI are the gold standards with constraints of accessibility, availability, and possibly some contraindications. The development of Ultrasound Localization Microscopy (ULM) has enabled new perspectives to conventional transcranial ultrasound imaging with increased sensitivity, penetration depth, and resolution. The possibility of volumetric imaging has increased the field-of-view and provided a more precise description of the microvascularisation. In this study, rats (n = 9) were subjected to thromboembolic ischemic stroke or intracerebral hemorrhages prior to volumetric ULM at the early phases after onsets. Although the volumetric ULM performed in the early phase of ischemic stroke revealed a large hypoperfused area in the cortical area of the occluded artery, it showed a more diffused hypoperfusion in the hemorrhagic model. Respective computations of a Microvascular Diffusion Index highlighted different patterns of perfusion loss during the first 24 h of these two strokes’ subtypes. Our study provides the first proof that this methodology should allow early discrimination between ischemic and hemorrhagic stroke with a potential toward diagnosis and monitoring in clinic.

A theranostic 3D ultrasound imaging system for high resolution image-guided therapy

Microbubble contrast agents are a diagnostic tool with broad clinical impact and an increasing number of indications. Many therapeutic applications have also been identified. Yet, technologies for ultrasound guidance of microbubble-mediated therapy are limited. In particular, arrays that are capable of implementing and imaging microbubble-based therapy in three dimensions in real-time are lacking. We propose a system to perform and monitor microbubble-based therapy, capable of volumetric imaging over a large field-of-view. To propel the promise of the theranostic treatment strategies forward, we have designed and tested a unique array and system for 3D ultrasound guidance of microbubble-based therapeutic protocols based on the frequency, temporal and spatial requirements. Methods: Four 256-channel plane wave scanners (Verasonics, Inc, WA, USA) were combined to control a 1024-element planar array with 1.3 and 2.5 MHz therapeutic and imaging transmissions, respectively. A transducer aperture of ~40×15 mm was selected and Field II was applied to evaluate the point spread function. In vitro experiments were performed on commercial and custom phantoms to assess the spatial resolution, image contrast and microbubble-enhanced imaging capabilities. Results: We found that a 2D array configuration with 64 elements separated by λ-pitch in azimuth and 16 elements separated by 1.5λ-pitch in elevation ensured the required flexibility. This design, of 41.6 mm × 16 mm, thus provided both an extended field-of-view, up to 11 cm x 6 cm at 10 cm depth and steering of ±18° in azimuth and ±12° in elevation. At a depth of 16 cm, we achieved a volume imaging rate of 60 Hz, with a contrast ratio and resolution, respectively, of 19 dB, 0.8 mm at 3 cm and 20 dB and 2.1 mm at 12.5 cm. Conclusion: A single 2D array for both imaging and therapeutics, integrated with a 1024 channel scanner can guide microbubble-based therapy in volumetric regions of interest.

Intensity distribution segmentation in ultrafast Doppler combined with scanning laser confocal microscopy for assessing vascular changes associated with ageing in murine hippocampi

The hippocampus plays an important role in learning and memory, requiring high-neuronal oxygenation. Understanding the relationship between blood flow and vascular structure—and how it changes with ageing—is physiologically and anatomically relevant. Ultrafast Doppler (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μDoppler) and scanning laser confocal microscopy (SLCM) are powerful imaging modalities that can measure in vivo cerebral blood volume (CBV) and post mortem vascular structure, respectively. Here, we apply both imaging modalities to a cross-sectional and longitudinal study of hippocampi vasculature in wild-type mice brains. We introduce a segmentation of CBV distribution obtained from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μDoppler and show that this mice-independent and mesoscopic measurement is correlated with vessel volume fraction (VVF) distribution obtained from SLCM—e.g., high CBV relates to specific vessel locations with large VVF. Moreover, we find significant changes in CBV distribution and vasculature due to ageing (5 vs. 21 month-old mice), highlighting the sensitivity of our approach. Overall, we are able to associate CBV with vascular structure—and track its longitudinal changes—at the artery-vein, venules, arteriole, and capillary levels. We believe that this combined approach can be a powerful tool for studying other acute (e.g., brain injuries), progressive (e.g., neurodegeneration) or induced pathological changes.

Refraction-Corrected Transcranial Ultrasound Imaging Through the Human Temporal Window Using a Single Probe

Transcranial ultrasound imaging (TUI) is a diagnostic modality with numerous applications, but unfortunately, it is hindered by phase aberration caused by the skull. In this article, we propose to reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. First, the compressional sound velocity of the aberrator (i.e., the skull) is estimated using the bidirectional headwave technique. The medium is described with four layers (i.e., lens, water, skull, and water), and a fast marching method calculates the travel times between individual array elements and image pixels. Finally, a delay-and-sum algorithm is used for image reconstruction with coherent compounding. The point spread function (PSF) in a wire phantom image and reconstructed with the conventional technique (using a constant sound speed throughout the medium), and the proposed method was quantified with numerical synthetic data and experiments with a bone-mimicking plate and a human skull, compared with the PSF achieved in a ground truth image of the medium without the aberrator (i.e., the bone plate or skull). A phased-array transducer (P4-1, ATL/Philips, 2.5 MHz, 96 elements, pitch = 0.295 mm) was used for the experiments. The results with the synthetic signals, the bone-mimicking plate, and the skull indicated that the proposed method reconstructs the scatterers with an average lateral/axial localization error of 0.06/0.14 mm, 0.11/0.13 mm, and 1.0/0.32 mm, respectively. With the human skull, an average contrast ratio (CR) and full-width-half-maximum (FWHM) of 37.1 dB and 1.75 mm were obtained with the proposed approach, respectively. This corresponds to an improvement of CR and FWHM by 7.1 dB and 36% compared with the conventional method, respectively. These numbers were 12.7 dB and 41% with the bone-mimicking plate.

Acoustic properties across the human skull

Transcranial ultrasound is emerging as a noninvasive tool for targeted treatments of brain disorders. Transcranial ultrasound has been used for remotely mediated surgeries, transient opening of the blood-brain barrier, local drug delivery, and neuromodulation. However, all applications have been limited by the severe attenuation and phase distortion of ultrasound by the skull. Here, we characterized the dependence of the aberrations on specific anatomical segments of the skull. In particular, we measured ultrasound propagation properties throughout the perimeter of intact human skulls at 500 kHz. We found that the parietal bone provides substantially higher transmission (average pressure transmission 31 ± 7%) and smaller phase distortion (242 ± 44 degrees) than frontal (13 ± 2%, 425 ± 47 degrees) and occipital bone regions (16 ± 4%, 416 ± 35 degrees). In addition, we found that across skull regions, transmission strongly anti-correlated (R=-0.79) and phase distortion correlated (R=0.85) with skull thickness. This information guides the design, positioning, and skull correction functionality of next-generation devices for effective, safe, and reproducible transcranial focused ultrasound therapies.

Aging-related cerebral microvascular changes visualized using ultrasound localization microscopy in the living mouse

Aging-related cognitive decline is an emerging health crisis; however, no established unifying mechanism has been identified for the cognitive impairments seen in an aging population. A vascular hypothesis of cognitive decline has been proposed but is difficult to test given the requirement of high-fidelity microvascular imaging resolution with a broad and deep brain imaging field of view, which is restricted by the fundamental trade-off of imaging penetration depth and resolution. Super-resolution ultrasound localization microscopy (ULM) offers a potential solution by exploiting circulating microbubbles to achieve a vascular resolution approaching the capillary scale without sacrificing imaging depth. In this report, we apply ULM imaging to a mouse model of aging and quantify differences in cerebral vascularity, blood velocity, and vessel tortuosity across several brain regions. We found significant decreases in blood velocity, and significant increases in vascular tortuosity, across all brain regions in the aged cohort, and significant decreases in blood volume in the cerebral cortex. These data provide the first-ever ULM measurements of subcortical microvascular dynamics in vivo within the context of the aging brain and reveal that aging has a major impact on these measurements.

Transcranial Neuromodulation Array With Imaging Aperture for Simultaneous Multifocus Stimulation in Nonhuman Primates

Even simple behaviors arise from the simultaneous activation of multiple regions in the brain. Thus, the ability to simultaneously stimulate multiple regions within a brain circuit should allow for better modulation of function. However, performing simultaneous multifocus ultrasound neuromodulation introduces challenges to transducer design. Using 3-D Fullwave simulations, we have designed an ultrasound neuromodulation array for nonhuman primates that: 1) can simultaneously focus on multiple targets and 2) include an imaging aperture for additional functional imaging. This design is based on a spherical array, with 128 15-mm elements distributed in a spherical helix pattern. It is shown that clustering the elements tightly around the 65-mm imaging aperture located at the top of the array improves targeting at shallow depths, near the skull surface. Spherical arrays have good focusing capabilities through the skull at the center of the array, but focusing on off-center locations is more challenging due to the natural geometric configuration and the angle of incidence with the skull. In order to mitigate this, the 64 elements closest to the aperture were rotated toward and focusing on a shallow target, and the 64 elements farthest from the aperture were rotated toward and focusing on a deeper target. Data illustrated that this array produced focusing on the somatosensory cortex with a gain of 4.38 and to the thalamus with a gain of 3.82. To improve upon this, the array placement was optimized based on phase aberration simulations, allowing for the elements with the largest impact on the gain at each focal point to be found. This optimization resulted in an array design that can focus on the somatosensory cortex with a gain of 5.19 and the thalamus with a gain of 4.45. Simulations were also performed to evaluate the ability of the array to focus on 28 additional brain regions, showing that off-center target regions can be stimulated, but those closer to the skull will require corrective steps to deliver the same amount of energy to those locations. This simulation and design process can be adapted to an individual monkey or human skull morphologies and specific target locations within individuals by using orientable 3-D printing of the transducer case and by electronic phase aberration correction.

Characterization of the Ultrasound Localization Microscopy Resolution Limit in the Presence of Image Degradation

Ultrasound localization microscopy (ULM) has been able to overcome the diffraction limit of ultrasound imaging. The resolution limit of ULM has been previously modeled using the Cramér-Rao lower bound (CRLB). While this model has been validated in a homogeneous medium, it estimates a resolution limit, which has not yet been achieved in vivo. In this work, we investigated the effects of three sources of image degradation on the resolution limit of ULM. The Fullwave simulation tool was used to simulate acquisitions of transabdominal contrast-enhanced data at depth. The effects of reverberation clutter, trailing clutter, and phase aberration were studied. The resolution limit, in the presence of reverberation clutter alone, was empirically measured to be up to 39 times worse in the axial dimension and up to 2.1 times worse in the lateral dimension than the limit predicted by the CRLB. While reverberation clutter had an isotropic impact on the resolution, trailing clutter had a constant impact on both dimensions across all signal-to-trailing-clutter ratios (STCR). Phase aberration had a significant impact on the resolution limit over the studied analysis ranges. Phase aberration alone degraded the resolution limit up to 70 and 160 [Formula: see text] in the lateral and axial dimensions, respectively. These results illustrate the importance of phase aberration correction and clutter filtering in ULM postprocessing. The analysis results were demonstrated through the simulation of the ULM process applied to a cross-tube model that was degraded by each of the three aforementioned sources of degradation.

3D Transcranial Ultrasound Localization Microscopy in the Rat Brain with a Multiplexed Matrix Probe

OBJECTIVE

Ultrasound Localization Microscopy (ULM) provides images of the microcirculation in-depth in living tissue. However, its implementation in two-dimension is limited by the elevation projection and tedious plane-by-plane acquisition. Volumetric ULM alleviates these issues and can map the vasculature of entire organs in one acquisition with isotropic resolution. However, its optimal implementation requires many independent acquisition channels, leading to complex custom hardware.

METHODS

In this article, we implemented volumetric ultrasound imaging with a multiplexed 32 x 32 probe driven by a single commercial ultrasound scanner. We propose and compare three different sub-aperture multiplexing combinations for localization microscopy in silico and in vitro with a flow of microbubbles in a canal. Finally, we evaluate the approach for micro-angiography of the rat brain.The "light" combination allows a higher maximal volume rate than the "full" combination while maintaining the field of view and resolution.

RESULTS

In the rat brain, 100,000 volumes were acquired within 7 min with a dedicated ultrasound sequence and revealed vessels down to 31 m in diameter with flows from 4.3 mm/s to 28.4 mm/s.

CONCLUSION

This work demonstrates the ability to perform a complete angiography with unprecedented resolution in the living rats brain with a simple and light setup through the intact skull.

SIGNIFICANCE

We foresee that it might contribute to democratize 3D ULM for both preclinical and clinical studies.

Transcranial Theranostic Ultrasound for Pre-Planning and Blood-Brain Barrier Opening: A Feasibility Study Using an Imaging Phased Array In Vitro and In Vivo

Focused ultrasound (FUS) for blood-brain barrier (BBB) opening is a safe, reversible and non-invasive strategy for targeted drug delivery to the brain, however extensive pre-planning strategies are necessary for successful FUS-mediated BBB opening through the structurally complex primate skull.

OBJECTIVE

This work aims to demonstrate a pre-planning pipeline consisting of transcranial simulations and in vitro experimentation used to inform synchronous BBB opening and power cavitation imaging (PCI) with a single theranostic ultrasound (TUS) phased array.

METHODS

Acoustic wave propagation simulation readouts of pressure attenuation and focal shift through clinical-CT and micro-CT-based primate skull models were compared, while the latter were used to determine the impact of beam steering angle on focal shift and pressure attenuation. In vitro experimentation with a channel phantom enabled characterization of skull-induced receive focal shift (RFS), while in vivo BBB opening and PCI using in silico and in vitro pre-planning information was conducted using a custom Verasonics/MATLAB script.

RESULTS

Simulations confirmed steering angle dependent transcranial focal shift and pressure attenuation, while in vitro experiments revealed minimal (0.30-1.50 mm) skull-induced RFS. In vivo rodent experiments with overlaid primate skull fragments demonstrated successful TUS-mediated BBB opening and spatially correlated power cavitation images (PCI) with regions of BBB opening on T1-weighted magnetic resonance images (MRI).

CONCLUSION

Herein, we demonstrate the feasibility for TUS-mediated BBB opening in vivo using in silico and in vitro pre-planning information.

SIGNIFICANCE

TUS as an ultrasound-guided modality for BBB opening could be a promising alternative to current FUS-mediated BBB opening configurations in the clinic.

Effect of Skull Porous Trabecular Structure on Transcranial Ultrasound Imaging in the Presence of Elastic Wave Mode Conversion at Varying Incidence Angle

With the advancement of aberration correction techniques, transcranial ultrasound imaging has exhibited great potential in applications such as imaging neurological function and guiding therapeutic ultrasound. However, the feasibility of transcranial imaging varies among individuals because of the differences in skull acoustic properties. To better understand the fundamental mechanisms underlying the variation in imaging performance, the effect of the structure of the porous trabecular bone on transcranial imaging performance (i.e., target localization errors and resolution) was investigated for the first time through the use of elastic wave simulations and experiments. Simulation studies using high-resolution computed tomography data from ex vivo skull samples revealed that imaging at large incidence angles reduced the target localization error for skulls having low porosity; however, as skull porosity increased, large angles of incidence resulted in degradation of resolution and increased target localization errors. Experimental results indicate that imaging at normal incidence introduced a localization error of 1.85 ± 0.10 mm, while imaging at a large incidence angle (40°) resulted in an increased localization error of 6.54 ± 1.33 mm and caused a single point target to no longer appear as a single, coherent target in the resulting image, which is consistent with simulation results. This first investigation of the effects of skull microstructure on transcranial ultrasound imaging indicates that imaging performance is highly dependent on the porosity of the skull, particularly at non-normal angles of incidence.

Transcranial ultrafast ultrasound localization microscopy of brain vasculature in patients

Changes in cerebral blood flow are associated with stroke, aneurysms, vascular cognitive impairment, neurodegenerative diseases and other pathologies. Brain angiograms, typically performed via computed tomography or magnetic resonance imaging, are limited to millimetre-scale resolution and are insensitive to blood-flow dynamics. Here we show that ultrafast ultrasound localization microscopy of intravenously injected microbubbles enables transcranial imaging of deep vasculature in the adult human brain at microscopic resolution and the quantification of haemodynamic parameters. Adaptive speckle tracking to correct for micrometric brain-motion artefacts and ultrasonic-wave aberrations induced during transcranial propagation allowed us to map the vascular network of tangled arteries to functionally characterize blood-flow dynamics at a resolution of up to 25 μm and to detect blood vortices in a small deep-seated aneurysm in a patient. Ultrafast ultrasound localization microscopy may facilitate the understanding of brain haemodynamics and of how vascular abnormalities in the brain are related to neurological pathologies.

Modeling Ultrasound Propagation in the Moving Brain: Applications to Shear Shock Waves and Traumatic Brain Injury

Traumatic brain injury (TBI) studies on the living human brain are experimentally infeasible due to ethical reasons and the elastic properties of the brain degrade rapidly postmortem. We present a simulation approach that models ultrasound propagation in the human brain, while it is moving due to the complex shear shock wave deformation from a traumatic impact. Finite difference simulations can model ultrasound propagation in complex media such as human tissue. Recently, we have shown that the fullwave finite difference approach can also be used to represent displacements that are much smaller than the grid size, such as the motion encountered in shear wave propagation from ultrasound elastography. However, this subresolution displacement model, called impedance flow, was only implemented and validated for acoustical media composed of randomly distributed scatterers. Herein, we propose a generalization of the impedance flow method that describes the continuous subresolution motion of structured acoustical maps, and in particular of acoustical maps of the human brain. It is shown that the average error in simulating subresolution displacements using impedance flow is small when compared to the acoustical wavelength ( λ /1702). The method is then applied to acoustical maps of the human brain with a motion that is imposed by the propagation of a shear shock wave. This motion is determined numerically with a custom piecewise parabolic method that is calibrated to ex vivo observations of shear shocks in the porcine brain. Then the fullwave simulation tool is used to model transmit-receive imaging sequences based on an L7-4 imaging transducer. The simulated radio frequency data are beamformed using a conventional delay-and-sum method and a normalized cross-correlation method designed for shock wave tracking is used to determine the tissue motion. This overall process is an in silico reproduction of the experiments that were previously performed to observe shear shock waves in fresh porcine brain. It is shown that the proposed generalized impedance flow method accurately captures the shear wave motion in terms of the wave profile, shock front characteristics, odd harmonic spectrum generation, and acceleration at the shear shock front. We expect that this approach will lead to improvements in image sequence design that takes into account the aberration and multiple reflections from the brain and in the design of tracking algorithms that can more accurately capture the complex brain motion that occurs during a traumatic impact. These methods of modeling ultrasound propagation in moving media can also be applied to other displacements, such as those generated by shear wave elastography or blood flow.

Ultrasound for the Brain: A Review of Physical and Engineering Principles, and Clinical Applications

The emergence of new ultrasound technologies has improved our understanding of the brain functions and offered new opportunities for the treatment of brain diseases. Ultrasound has become a valuable tool in preclinical animal and clinical studies as it not only provides information about the structure and function of brain tissues but can also be used as a therapy alternative for brain diseases. High-resolution cerebral flow images with high sensitivity can be acquired using novel functional ultrasound and super-resolution ultrasound imaging techniques. The noninvasive treatment of essential tremors has been clinically approved and it has been demonstrated that the ultrasound technology can revolutionize the currently existing treatment methods. Microbubble-mediated ultrasound can remotely open the blood-brain barrier enabling targeted drug delivery in the brain. More recently, ultrasound neuromodulation received a great amount of attention due to its noninvasive and deep penetration features and potential therapeutic benefits. This review provides a thorough introduction to the current state-of-the-art research on brain ultrasound and also introduces basic knowledge of brain ultrasound including the acoustic properties of the brain/skull and engineering techniques for ultrasound. Ultrasound is expected to play an increasingly important role in the diagnosis and therapy of brain diseases.

Imaging With Therapeutic Acoustic Wavelets-Short Pulses Enable Acoustic Localization When Time of Arrival is Combined With Delay and Sum

Passive acoustic mapping (PAM) is an algorithm that reconstructs the location of acoustic sources using an array of receivers. This technique can monitor therapeutic ultrasound procedures to confirm the spatial distribution and amount of microbubble activity induced. Current PAM algorithms have an excellent lateral resolution but have a poor axial resolution, making it difficult to distinguish acoustic sources within the ultrasound beams. With recent studies demonstrating that short-length and low-pressure pulses-acoustic wavelets-have the therapeutic function, we hypothesized that the axial resolution could be improved with a quasi-pulse-echo approach and that the resolution improvement would depend on the wavelet's pulse length. This article describes an algorithm that resolves acoustic sources axially using time of flight and laterally using delay-and-sum beamforming, which we named axial temporal position PAM (ATP-PAM). The algorithm accommodates a rapid short pulse (RaSP) sequence that can safely deliver drugs across the blood-brain barrier. We developed our algorithm with simulations (k-wave) and in vitro experiments for one-, two-, and five-cycle pulses, comparing our resolution against that of two current PAM algorithms. We then tested ATP-PAM in vivo and evaluated whether the reconstructed acoustic sources mapped to drug delivery within the brain. In simulations and in vitro, ATP-PAM had an improved resolution for all pulse lengths tested. In vivo, experiments in mice indicated that ATP-PAM could be used to target and monitor drug delivery into the brain. With acoustic wavelets and time of flight, ATP-PAM can locate acoustic sources with a vastly improved spatial resolution.

Ultrasound Technologies for Imaging and Modulating Neural Activity

Visualizing and perturbing neural activity on a brain-wide scale in model animals and humans is a major goal of neuroscience technology development. Established electrical and optical techniques typically break down at this scale due to inherent physical limitations. In contrast, ultrasound readily permeates the brain, and in some cases the skull, and interacts with tissue with a fundamental resolution on the order of 100 μm and 1 ms. This basic ability has motivated major efforts to harness ultrasound as a modality for large-scale brain imaging and modulation. These efforts have resulted in already-useful neuroscience tools, including high-resolution hemodynamic functional imaging, focused ultrasound neuromodulation, and local drug delivery. Furthermore, recent breakthroughs promise to connect ultrasound to neurons at the genetic level for biomolecular imaging and sonogenetic control. In this article, we review the state of the art and ongoing developments in ultrasonic neurotechnology, building from fundamental principles to current utility, open questions, and future potential.

Transcranial activation and imaging of low boiling point phase-change contrast agents through the temporal bone using an ultrafast interframe activation ultrasound sequence

PURPOSE

As a cavitation enhancer, low boiling point phase-change contrast agents (PCCA) offer potential for ultrasound-mediated drug delivery in applications including intracerebral hemorrhage or brain tumors. In addition to introducing cavitation, ultrasound imaging also has the ability to provide guidance and monitoring of the therapeutic process by localizing delivery events. However, the highly attenuating skull poses a challenge for achieving an image with useful contrast. In this study, the feasibility of transcranial activation and imaging of low boiling point PCCAs through the human temporal bone was investigated by using a low frequency ultrafast interframe activation ultrasound (UIAU) imaging sequence with singular value decomposition-based denoising.

METHODS

Lipid-shelled PCCAs filled with decafluorobutane were activated and imaged at 37°C in tissue-mimicking phantoms both without and with an ex vivo human skull using the new UIAU sequence and a low frequency diagnostic transducer array at frequencies from 1.5 to 3.5 MHz. A singular value decomposition-based denoising filter was developed and used to further enhance transcranial image contrast. The contrast-to-tissue ratio (CTR) and contrast enhancement (CE) of UIAU was quantitatively evaluated and compared with the amplitude modulation pulse inversion (AMPI) and vaporization detection imaging (VDI) approaches.

RESULTS

Image results demonstrate enhanced contrast in the phantom channel with suppressed background when the low boiling point PCCA was activated both without and with an ex vivo human skull using the UIAU sequence. Quantitative results show that without the skull, low frequency UIAU imaging provided significantly higher image contrast (CTR ≥ 18.56 dB and CE ≥ 18.66 dB) than that of AMPI and VDI (P < 0.05). Transcranial imaging results indicated that the CE of UIAU (≥18.80 dB) was significantly higher than AMPI for free-field activation pressures of 5 and 6 MPa. The CE of UIAU is also significantly higher than that of VDI when PCCAs were activated at 2.5 MHz and 3 MHz (P < 0.05). The CTR (23.30 [20.07-25.56] dB) of denoised UIAU increased by 12.58 dB relative to the non-denoised case and was significantly higher than that of AMPI at an activation pressure of 4 MPa (P < 0.05).

CONCLUSIONS

Results indicate that low boiling point PCCAs can be activated and imaged at low frequencies including imaging through the temporal bone using the UIAU sequence. The UIAU imaging approach provides higher contrast than AMPI and VDI, especially at lower activation pressures with additional removal of electronic noise.

Early Ultrafast Ultrasound Imaging of Cerebral Perfusion correlates with Ischemic Stroke outcomes and responses to treatment in Mice

In the field of ischemic cerebral injury, precise characterization of neurovascular hemodynamic is required to select candidates for reperfusion treatments. It is thus admitted that advanced imaging-based approaches would be able to better diagnose and prognose those patients and would contribute to better clinical care. Current imaging modalities like MRI allow a precise diagnostic of cerebral injury but suffer from limited availability and transportability. The recently developed ultrafast ultrasound could be a powerful tool to perform emergency imaging and long term follow-up of cerebral perfusion, which could, in combination with MRI, improve imaging solutions for neuroradiologists. Methods: In this study, in a model of in situ thromboembolic stroke in mice, we compared a control group of non-treated mice (N=10) with a group receiving the gold standard pharmacological stroke therapy (N=9). We combined the established tool of magnetic resonance imaging (7T MRI) with two innovative ultrafast ultrasound methods, ultrafast Doppler and Ultrasound Localization Microscopy, to image the cerebral blood volumes at early and late times after stroke onset and compare with the formation of ischemic lesions. Results: Our study shows that ultrafast ultrasound can be used through the mouse skull to monitor cerebral perfusion during ischemic stroke. In our data, the monitoring of the reperfusion following thrombolytic within the first 2 h post stroke onset matches ischemic lesions measured 24 h. Moreover, similar results can be made with Ultrasound Localization Microscopy which could make it applicable to human patients in the future. Conclusion: We thus provide the proof of concept that in a mouse model of thromboembolic stroke with an intact skull, early ultrafast ultrasound can be indicative of responses to treatment and cerebral tissue fates following stroke. It brings new tools to study ischemic stroke in preclinical models and is the first step prior translation to the clinical settings.

Super-resolution Ultrasound Imaging

The majority of exchanges of oxygen and nutrients are performed around vessels smaller than 100 μm, allowing cells to thrive everywhere in the body. Pathologies such as cancer, diabetes and arteriosclerosis can profoundly alter the microvasculature. Unfortunately, medical imaging modalities only provide indirect observation at this scale. Inspired by optical microscopy, ultrasound localization microscopy has bypassed the classic compromise between penetration and resolution in ultrasonic imaging. By localization of individual injected microbubbles and tracking of their displacement with a subwavelength resolution, vascular and velocity maps can be produced at the scale of the micrometer. Super-resolution ultrasound has also been performed through signal fluctuations with the same type of contrast agents, or through switching on and off nano-sized phase-change contrast agents. These techniques are now being applied pre-clinically and clinically for imaging of the microvasculature of the brain, kidney, skin, tumors and lymph nodes.