New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Brain palpation from physiological vibrations using MRI
Edited by David A. Weitz, Harvard University, Cambridge, MA, and approved September 1, 2015 (received for review May 21, 2015)

Significance
It is commonly supposed that noise obscures but does not contain useful information. However, in wave physics and especially, seismology, scientists developed some tools known as “noise correlation” to extract useful information and construct images from the random vibrations of a medium. Living tissues are full of unexploited vibrations as well. In this manuscript, we show that noise correlation techniques in the brain using MRI can conduct to a tomography related to the stiffness that physicians estimate during a palpation examination. This MRI brain seismology in rupture with other techniques is experimentally shown in vitro and in vivo.
Abstract
We present a magnetic resonance elastography approach for tissue characterization that is inspired by seismic noise correlation and time reversal. The idea consists of extracting the elasticity from the natural shear waves in living tissues that are caused by cardiac motion, blood pulsatility, and any muscle activity. In contrast to other magnetic resonance elastography techniques, this noise-based approach is, thus, passive and broadband and does not need any synchronization with sources. The experimental demonstration is conducted in a calibrated phantom and in vivo in the brain of two healthy volunteers. Potential applications of this “brain palpation” approach for characterizing brain anomalies and diseases are foreseen.
The complexity of wave fields can sometimes be an advantage for imaging. Such is the case in multiple scattering or reverberating media, where wave fields contain information about their sources and the medium itself. Turning this wave noise into useful measurements through correlation techniques has provided a breakthrough in a wide variety of domains, which range from seismology (1) to acoustics (2, 3) and electromagnetism (4). Living tissue is also full of unexploited vibrations. Their detection with ultrafast ultrasound scanners that can reach thousands of frames per second (5⇓–7) has recently opened up the medical field to correlation techniques and therefore, passive elastography (8). However, ultrasound is not suitable for brain imaging. MRI can image the brain, but its relatively low acquisition rate of a few frames per second is an issue. Synchronization with the shear-wave source is, thus, necessary (9, 10), which penalizes its potential implementation based on natural shear waves.
We describe a magnetic resonance elastography (MRE) -based method that is free from the need for synchronization and any controlled source. This approach extracts information related to the mechanical properties of the soft tissue from hundreds of snapshots of randomly fluctuating shear-wave fields. The key to decrypting the complex field is correlation or similarly but from a physical point of view, time reversal (11, 12). Not only does this wide-band approach maximize the signal-to-noise ratio, as any matched filter would, but it also avoids the Nyquist–Shannon problem that is inherent to slow imaging devices. Indeed, although time information is definitely lost, the spatial information is still present and allows shear-wavelength tomography to be conducted. This wavelength is closely related to the shear elasticity and thus, the intuitive estimation of the stiffness felt by physicians during palpation examination. To start, this concept is shown using MRI in a calibrated elastography phantom under randomly sampled vibrations. Arterial pulsation can produce motion in the brain as high as 1 mm (13); the resulting natural shear-wave field is analyzed through correlation algorithms, and passive brain palpation reconstructions are presented.
Results
Phantom Experiments.
As shown in the experimental setup (Fig. 1A), vibrations are generated synchronously from 12 contact points of a device mounted on a piezoelectric source and working with a 80.3-s modulated sinusoid excitation. The emitted frequencies range from 20 to 200 Hz. These vibrations propagate mainly as shear waves within a phantom designed for elastography tests (Computerized Imaging Reference Systems, Inc.) that contains four calibrated inclusions. The emitted waves reverberate inside this soft, solid cavity, and the resulting complex field is measured every 1.5 s using an MRI scanner. It is important to point out first that the vibrations driving are continuously repeated during the whole acquisition time and second, that the modulated sinusoid duration of 80.3 s divided by the imaging frame rate of 1.5 s is not an integer to avoid any cyclic sampling of the field. Each frame contains 128 × 128 pixels. The acquisition of 320 frames for the out-of-plane component of the displacement field was carried out at ∼0.7 Hz on six slices that were equally spaced by 3 mm, with 3 × 3 × 3-mm3 isotropic voxels, over 8 min. According to the Nyquist–Shannon limit, the maximum shear-wave frequency was, thus, ∼0.35 Hz, which is far below the central frequency of shear waves in the experiment. The first six consecutive snapshots of the displacement field of a total of 320 snapshots are shown in Fig. 1B. These snapshots clearly show the rectangular phantom surrounded by the random-phase noise. The elastic field shows complex spatial patterns with characteristic dimensions driven by the wavelengths associated with the 20- to 200-Hz frequency content of the sweep source signal. Small wavelengths are visible at t = 1.5 s, whereas larger wavelengths dominate at t = 3 s. A correlation coefficient smaller than 5% between consecutive snapshots of wave pattern confirms the random property of the diffuse elastic field. The extraction of the elastic parameters from 320 snapshots was then carried out using noise correlation.
(A) Experimental setup. The complex shear-wave field is generated from 12 point sources located at the surface and mounted on an MRI-compatible piezoelectric vibrator. The excitation signal is an 80.3-s modulated sinusoid, with the frequency range of 20 to 200 Hz. (B) Phase images that represent the out-of-plane displacement field uy in the zx plane of the inclusions. These first six snapshots were extracted from a movie of 320 images that was acquired using a gradient echo MRE sequence every 1.5 s.
Theory.
Let
Phantom correlation results. (A) Time-reversal focal spots for three different virtual source locations S in (Top) the phantom background (25 kPa), (Middle) the hard inclusion (80 kPa), and (Bottom) the soft inclusion (8 kPa). The virtual source S can be set arbitrarily to any location. (B) T2-weighted image of the phantom. (C) Shear-wavelength tomography. The four inclusions are clearly apparent, and their local wavelength estimations are proportional to the square root of their local elasticity.
The dimensions of the focal spot vary according to the local stiffness of the medium at the location
Time-reversal focal spots on three different virtual source locations S in the inclusions presenting different Young moduli: hard (80 kPa), medium hard (45 kPa), medium soft (14 kPa), and soft (8 kPa). An accurate estimation of the focus spot size resulting from correlation follows its theoretical dependence as the square root of the Young modulus.
Shear-wavelength tomography. The four 2 cm-sized inclusions are clearly apparent in the four different slices equally spaced by 3 mm. The reconstruction in Lower Left is presented in the article (Fig. 2C).
Brain Experiments.
The aim of passive elastography is to use natural shear-wave propagation in the human body. In the brain, the natural motion is caused by arterial pulsatility and cerebrospinal fluid exchange (23, 24). The MRE sequence presented in this section is the same as the one described for the phantom, except for the three following parameters: the voxel size of six slices is 4 × 4 × 4 mm3 to improve the signal to noise ratio (SNR), the total time of the in vivo experiments is 11 min, and three components of the displacements are acquired. The out-of-plane component only is investigated here. The first six 64 × 64-pixel snapshots shown in Fig. 3A show that the MRE sequence described in Materials and Methods can reveal the natural motion in the brain. The amplitudes are relatively small at only 5 µm, but the signal is clearly apparent. Because correlation computation is equivalent to an adaptive filter, the focal spot clearly emerges from the noise (Fig. 3B). The existence of these time-reversal focus spots is a convincing clue that natural motions in the brain are transported by shear waves. They are responsible for correlation decreasing down to negative values. The black −6-dB isolevel boundary lines more clearly illustrate the different sizes of the focal spots. With a mean of 10 cm, the wavelength estimations are within the expected range. Indeed, assuming that the motion is proportional to the intracranial pressure (25) and given the frequency-dependent sensitivity of the 50-Hz motion-encoding gradient used in our in vivo experiment, the central frequency of displacement (as encoded in the phase images) is expected to be ∼1 Hz (13); 10% of energy is nonetheless present at a frequency of 15 Hz. It is expected that the wavelength estimation is dominated by the higher frequencies. Therefore, the shear-wave speed in the brain is around 1.5 m/s (26, 27). Therefore, the expected smaller wavelength ranges from 8 to 15 cm. In addition, longer wavelengths are observed on the periphery of the frontal and parietal lobes than in the center of the brain (Fig. 3B), which is clearly apparent on the shear-wavelength tomography reconstructions (Fig. 4). The anatomic details of T2 images are apparent on the shear-wavelength tomography, especially along the longitudinal fissure in the axial view (Fig. 4A). Finally, the longest wavelengths are visible in the posterior region of the brain (Fig. 4).
In vivo brain results. (A) Phase representation of the displacement field. Six snapshots were extracted from a movie of 144 images. These in vivo measurements of brain motion were acquired with a gradient echo MRE sequence every 1.5 s. (B) Time-reversal focal spots in the brain for four different virtual source locations S. The isolevel black boundary lines stress the shear-wavelength variations.
In vivo brain passive MRE. (A) Axial view of (Left) the T2-weighted image and (Right) its corresponding shear-wavelength tomography. (B) Sagittal view of (Left) the T2-weighted image and (Right) its corresponding shear-wavelength tomography.
Discussion
The general results extracted from the first brain correlation reconstructions were verified on the same volunteer at a 1-mo interval as well as a second healthy volunteer (Figs. S3 and S4). Healthy volunteers were on their backs. The gravity-induced compression might enhance the local wavelength in the posterior brain. Additional studies of passive MRE are necessary to claim with confidence that these red posterior regions of the wavelength tomography reconstructions correspond to stiffer brain tissue (Fig. 4 and Fig. S4). Similarly, the apparent stiffer periphery region is not confirmed by MRE observations of the literature. Because it has to be a wave property effect, it is conjectured to be either a nonuniform frequency content of the shear-wave field or a change in the nature of the elastic wave with internal structure. At last, given that MRE includes three components of the field, that two orthogonal wavelength estimations are possible, and that the acquisition involves six different planes, we, therefore, end up with 36 possible tomography reconstructions. It is beyond the scope of this paper to offer detailed descriptions of each of these, although they are partly considered in Figs. S1, S2, S3, and S4. A promising possibility resides in the wave vector estimation that should partly reveal the brain anisotropy. In conclusion, the whole set of these experimental data brings us to the following statement: although natural tissue vibrations are often considered as disturbing noise, they nonetheless offer promising opportunities for medical imaging when correlation methods can be applied.
Time-reversal focal spots on four different virtual source locations S. They are measured from a sagittal slice in the brain. Isolevel black curves stress the shear-wavelength variations. The corresponding shear-wavelength tomography is shown in the article (Fig. 4B).
Comparison between two volunteers of (Left) the axial view of the T2-weighted image and (Right) its corresponding shear-wavelength tomography. The tomography on volunteer 2 was acquired 1 mo after the one presented in the manuscript (Fig. 4A). Although in different regions of the brain, the similarity between reconstructions is striking.
Materials and Methods
Instantaneous MRI snapshots are captured using a gradient echo single-shot echo planar imaging sequence (28) developed at the Mayo Clinic and installed on a 1.5-T system (General Electric). The displacement is encoded in the phase of the magnetic resonance signal using a 20-mT/m bipolar motion-encoding gradient applied over 10 ms (i.e., +20 mT/m over 5 ms immediately followed by −20 mT/m over 5 ms). The repetition and echo times are 1.5 s and 45 ms, respectively. In the phantom experiment, the field of view is 128 × 128 × 3 mm, whereas it is 64 × 64 × 4 mm in the brain experiment. The phantom is designed for elastography experiments by Computerized Imaging Reference Systems, Inc. The piezoelectric vibrator used in the phantom experiment is a PPA20M Cedrat Technologies. Institutional review board approval from Inserm was given and written informed consent was obtained from both volunteers.
Acknowledgments
We acknowledge R. Ehman and K. Glaser from the Mayo Clinic for making the magnetic resonance elastography (MRE) sequences available.
Footnotes
- ↵1To whom correspondence should be addressed. Email: stefan.catheline{at}inserm.fr.
Author contributions: R.S., J.-Y.C., and S.C. designed research; A.Z., R.S., A.-H.D., S.L., and S.C. performed research; J.-M.M. and O.R. analyzed data; and A.Z. and S.C. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1509895112/-/DCSupplemental.
References
- ↵
- ↵
- ↵
- ↵.
- Lerosey G,
- de Rosny J,
- Tourin A,
- Fink M
- ↵.
- Sandrin L,
- Catheline S,
- Tanter M,
- Hennequin X,
- Fink M
- ↵
- ↵
- ↵
- ↵.
- Muthupillai R, et al.
- ↵
- ↵
- ↵.
- Fink M
- ↵
- ↵.
- Campillo M,
- Paul A
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Applied Physical Sciences
- Biological Sciences
- Medical Sciences