Characterization of oscillations in the brain and cerebrospinal fluid using ultra-high field magnetic resonance imaging

Tiago Martins, Tales Santini, Bruno de Almeida, Minjie Wu, Kristine A. Wilckens, Davneet Minhas, James W. Ibinson, Howard J. Aizenstein, Tamer S. Ibrahim
doi: https://doi.org/10.1101/2023.12.05.23299452

Abstract

Development of innovative non-invasive neuroimaging methods and biomarkers are critical for studying brain disease. In this work, we have developed a methodology to characterize the frequency responses and spatial localization of oscillations and movements of cerebrospinal fluid (CSF) flow in the human brain. Using 7 Tesla human MRI and ultrafast echo-planar imaging (EPI), in-vivo images were obtained to capture CSF oscillations and movements. Physiological data was simultaneously collected and correlated with the 7T MR data. The primary components of CSF oscillations were identified using spectral analysis (with frequency bands identified around 0.3Hz, 1.2Hz and 2.4Hz) and were mapped spatially and temporally onto the MR image domain and temporally onto the physiological domain. The developed methodology shows a good consistency and repeatability (standard deviation of 0.052 and 0.078 for 0.3Hz and 1.2Hz bands respectively) in-vivo for potential brain dynamics and CSF flow and clearance studies.

1. Introduction

Clearance and exchange of brain fluids promote brain health by removing neurotoxic metabolic byproducts from the brain such as amyloid beta and tau (Nedergaard, 2013; Xie et al., 2013). Clearance of brain fluids is driven by pulsations of the perivascular space from autonomic nervous system (ANS) activity, which may vary as a function of brain states such as sleep and wakefulness (Hauglund et al., 2020; Herculano-Houzel, 2013; Xie et al., 2013) as well as brain diseases such as Alzheimer’s disease (Peng et al., 2016; Ramanathan et al., 2015) and major depressive disorder(Hock et al., 1998; Pomara et al., 2012). The lymphatic draining system of the brain tissue, also known as the glymphatic system, is described as convection of cerebrospinal fluid (CSF) between the peri-arterial and peri-venous spaces. This convective flow is partially due to the cardiac-induced blood flow pulsations along the arteries (Adolph et al., 1967; Iliff et al., 2013; Martin et al., 2012; Schroth & Klose, 1992). Water is propelled by the arterial pulsations through aquaporin channels and supports solute transport from extracellular interstitial spaces, through perivascular spaces, and into CSF spaces. CSF and waste products from the brain are then pushed from parenchyma to subarachnoid spaces and eventually cleared via arachnoid granulations and dural and nasal lymphatic vessels (Iliff et al., 2013; Kiviniemi et al., 2016; Leon et al., 2017; Rennels et al., 1990).

Development of quantitative CSF imaging methods is critical to understand factors that influence brain fluid clearance. T1-weighted magnetic resonance imaging (MRI) with intrathecal injection of a gadolinium (Gd)-based contrast agent has been used to characterize CSF flow in human participants with idiopathic normal pressure hydrocephalus (iNPH) and dementia (Eide & Ringstad, 2019; Ringstad et al., 2017, 2018). This technique has afforded fully quantitative, high-resolution imaging of CSF and interstitial fluid (ISF) flow throughout the whole head but is highly invasive due to necessary lumbar puncture. Gd may also deposit in the brain, inhibiting longitudinal study (Gulani et al., 2017).

Fast acquisition functional magnetic resonance imaging (fMRI) paradigms have also been used to characterize CSF dynamics in iNPH and Alzheimer’s disease (AD) patients and during sleep (Fultz et al., 2019; Shanks et al., 2019; Yamada et al., 2020; Yang et al., 2022). While non-invasive, these sequences have relatively poor signal-to-noise ratio and spatial resolution and have been limited to narrow fields of view encompassing only the 4th ventricular or cerebral aqueduct.

The bulk changes in blood volume at the capillary level could cause widespread fluctuations of measured signal intensity with the cardiac cycle. Furthermore, large vessel pulsatility may cause tissue movement and production of an influx of unsaturated blood into the slice of interest affecting the measured signal intensity in the areas adjacent to the vessels. This leads to a signal variation when using T2* echo-planar imaging (EPI) acquisitions (Dagli et al., 1999). Hence, fMRI and other techniques have been used to characterize different sources of pulsations in the brain (Biswal et al., 1995; Dagli et al., 1999; Kiviniemi et al., 2000; Poncelet et al., 1992; Purdon & Weisskoff, 1998). Thus, using MRI of CSF dynamics can highly inform the study of brain diseases and the role of sleep-wake states (Fultz et al., 2019; Xie et al., 2013).

Ultra-high field MRI (≥ 7 Tesla) provides a major advantage of increased signal-to-noise ratio (SNR). The enhanced SNR can be used either to increase the resolution of the images or to decrease the scanning time (with the use of higher acceleration factors). Other advantages of 7 Tesla (T) field strength are the higher sensitivity to blood-oxygen-level-dependent (BOLD) signal and better vasculature conspicuity (Moser et al., 2012; Santini, Wood, et al., 2021). The high signal-to-noise ratio (SNR) and fast acquisitions of 7T MRI scanners allow studies to perform analysis of blood and cerebrospinal fluid (CSF) flow within the brain (Scouten & Constable, 2008).

This work revolves around the creation of a method of acquisition and analysis that can be used as biomarker for study of central nervous system functioning and brain diseases. Using fast EPI and physiological acquisitions, we viewed and analyzed the CSF MR signal in real-time. We report CSF oscillation patterns through spectral analysis and apply the same methodology across different datasets to validate the observed results.

2. Methods

The overall design of this study is based on two main steps: in-vivo 7T image acquisition, concurrent with physiology measurements, and image processing with spectral analysis. The processing and analysis of the power-frequency spectrum and its corresponding spatial mappings were fine tuned for the processing parameters such as detection bandwidth and peak span, thresholding levels for masks and smoothing degrees for filtering.

The volunteers scanned for this work provided informed consent as part of an approved study by the University of Pittsburgh’s Institutional Review Board (identification number PRO17030036). For CSF flow data collection, five healthy volunteers (all females, age range 21-25 years old) were scanned to obtain EPI data from the whole brain, including the cerebellum. For comparison between physiological and MRI data, one additional healthy volunteer (male, age range 26-30 years old) was scanned and had physiological data collected simultaneously.

2.1 Image Acquisition

Images were acquired using a whole-body 7T MRI system (Siemens 7T MAGNETON) and with the human-connectome EPI multiband MR sequence (Moeller et al., 2010; Uğurbil et al., 2013). The sequence is capable of achieving fast acquisitions, high contrast to the CSF flow, and high sensitivity to BOLD signal, thus making it a good candidate for studies of sleep and neurodegenerative and psychological disorders. The imaging was done with an in-house developed and fabricated 16-channel Tic-Tac-Toe transmit array with a 32-channel receive head coil (Krishnamurthy et al., 2019; Santini et al., 2018; Santini, Wood, et al., 2021) that is load insensitive (Ibrahim et al., 2008; Kim et al., 2016; Santini, Wood, et al., 2021) and capable of whole brain homogenous imaging at 7T (Ibrahim et al., 2013). By using this coil design, we were able to acquire signal from the entire brain with minimal regions of significant excitation losses using the single transmit mode of the 7T scanner.

The acquired EPI images yield a real-time visualization of the CSF flow in the brain. The sequence was optimized to perform fast brain imaging. The main data acquisition was done with echo time (TE) of 20 ms (carefully chosen for future potential BOLD analysis), repetition time (TR) of 155ms, isotropic resolution of 2mm, and acceleration factor of 2. The field of view (FOV) was 192 mm x 192 mm x 6 mm. The acquisition was broken into 19 slabs of 3 axial slices each for a total of 57 slices, which is enough to have whole-brain coverage. A total of 600 volumes were sequentially acquired per slab in a single sequence run for an acquisition time of 1 minute and 36 seconds per slab.

During the development of the protocol, the EPI acquisition on the Volunteer 1 was done using TR of 152ms but only 15 slabs for a total of 45 axial slices. Another EPI acquisition on the same volunteer was also performed using TR of 51ms and a single slice.

Two spin-echo echo-planar images were also acquired for B0 field distortion correction with the same phase encoding (PE) direction of the EPI acquisition (PA) and with the opposite PE direction (AP). These acquisitions were performed for 57 slices, TE of 39.4ms, TR of 6000ms with matched parameters as the EPI sequence in terms of field of view, resolution, number of slices, echo spacing, and position.

A T1-weighted imaging (MPRAGE) sequence was used for proper localization of the EPI field of view and as a structural scan for the image processing. This acquisition was done using 0.75mm isotropic resolution, TR of 3000ms, TE of 2.17ms, and 256 slices for a coverage of 240 mm x 173 mm x 192 mm in total time of acquisition of ∼5 minutes.

2.2 Physiological Measurements

Electrocardiogram (ECG) and respiration signals were collected from a consented volunteer inside the MR scanner using MR compatible ECG leads and an expansion belt attached to the chest to track inflation and deflation of the chest during respiration. Acquisition was digitalized using BIOPAC system (“ECG,” n.d.; “Respiration Transducer for MRI | TSD221-MRI | Research | BIOPAC,” n.d.). The simultaneously collected data allows temporal signal analysis of both MR and physiologic signals. The imaging data acquired in conjunction with the physiological data used TR of 75ms, TE of 28ms, and slice thickness of 4 mm.

2.3 Image Processing

The processing pipeline was developed based on MATLAB (MATLAB - MathWorks, n.d.), ANTs (Avants et al., 2009), and FSL (Jenkinson et al., 2012) software. It consists of denoising, distortion correction, bias correction, and skull stripping of each dataset volume. The initial step was loading the slabs and merging them into a single dataset. Next, denoising was performed using a noise estimation tool with variance stabilization transformation (VST) for Rician-distributed noise (Foi, 2011). The Rician heteroscedastic noise is converted to a homoscedastic noise after the forward VST. The block-matching 4D (BM4D) denoising algorithm (Maggioni et al., 2013) can then be applied and the denoised image is obtained after the inverse VST. This tool has been used for other MRI applications(Santini, Koo, et al., 2021) and yields a good result when applied to EPI data. Distortion correction was done using the estimated B0 maps derived from the spin-echo sequence using topup tool (Andersson et al., 2003) from FSL software package. The generated map was used for correction of the EPI data. Then, the images were bias corrected using the N4(Tustison et al., 2010) tool from the ANTs software package with spline distance parameter of 200. The final skull stripping was performed using the FSL brain extraction tool (BET).

2.4 Spectral Analysis

The frequency analysis was performed for each dataset individually and resulted in both a frequency power spectrum and a mask for brain localization of specific frequency bands. A frequency spectrum was calculated for selected points for validation of the findings across different brain regions.

After processing each individual slab of the EPI data, the frequency processing and analysis were performed in MATLAB and Python. The time series of each voxel was used to generate a frequency spectrum using fast Fourier transform (FFT). With the 600 volumes of 155ms TR acquisition, the frequency resolution of the frequency spectrum is 0.011Hz and the maximum frequency is 3.23Hz. The same frequency analysis for the 51ms TR data produces a much larger frequency spectrum of up to 9.8Hz. Therefore, frequency components higher than 3Hz can be observed and analyzed. The analysis was done using both the average of the 3D space and individual points manually selected.

Spatial analysis was done by creating image masks based on the localization of voxels with significant signal in each frequency band. Each frequency band was defined with a bandwidth of 0.3Hz and was centered at the local maxima of the frequency spectrum amplitude with a minimum peak distance 0.15Hz (0.5Hz was used for the dataset with 51ms TR). The power map of a given frequency band was determined voxel-wise by averaging respective power values within the frequency band. For better visualization, each power map was then binarized with a chosen threshold (75% of the peak amplitude of the corresponding frequency band) and spatially smoothed using a Gaussian filter (sigma of 1.6), generating the final masks for each frequency band. These masks were overlaid on the original EPI and T1 weighted acquisitions for anatomical reference. The T1 weighted image was registered with the average EPI image of the dataset using SPM12.

3. Results

A video was created based on the image series of the fast EPI data after acquisition and processing (Figure 1). The video visually indicates the presence of periodical signal from the CSF flow. To confirm the presence of physiological signals such as respiration and cardiac motion, CSF temporal data was aligned with measurements from the electrocardiogram and respiration belt for visual comparison of similarity between the physiological activities and the change in signal intensity from CSF regions (Figure 2a). The frequency spectrum of the datasets was also aligned following the same comparison as the time series data (Figure 2b). The two major signal bands were highlighted between the ECG and CSF data (around 1.1Hz) and the respiration belt and CSF data (around 0.3Hz).

Figure 1:
Figure 1:

Fast EPI acquisition (TR=100ms) showing signal changes due to CSF flow; axial slices with a spatial resolution of 1.53 x 1.53 x 3mm and a sagittal slice with spatial resolution of 1.5 x 1.5 x 4.4mm. The blue arrows point to regions of large variation in signal over time. A video showing these oscillations in real time is available at https://doi.org/10.6084/m9.figshare.24022932

Figure 2:
Figure 2:

Data acquisition along the cardiac cycle. Echo-planar imaging acquisition performed with concurrent physiological measurement of electrocardiogram in the 7T scanner. a) Time series of ECG, respiration belt, and EPI signals. The EPI signal was temporally aligned with the physiological data using an external trigger signal from the scanner. The red lines represent the R-peaks of the ECG. B) Frequency spectrum of the ECG, respiration belt, and EPI signals. The purple region highlights the common frequency band between the respiration and EPI signals, and the orange region highlights the common frequency band between the ECG and EPI signals.

To verify that various points of the brain contribute differently on the frequency spectrum, Figure 3 represents the frequency spectrum for 9 arbitrary points throughout the brain. The position of each point is described by the brain anatomy that it belongs to as shown on the top-right corner of each spectrum graph. Most of the points show frequencies around 1.2Hz. Depending on the position, the signal shows the 0.3Hz and/or the 2.4Hz bands.

Figure 3:
Figure 3:

Frequency spectrum for 9 selected points (a-i) throughout the brain for volunteer 1. Some points show higher intensity on the 1.2Hz and 2.4Hz bands (points a, b, e, f, g, h, and i) whereas other show more on the 0.3Hz band (points c, d, and g). The labels describe the brain anatomy where the data was obtained.

A frequency analysis was obtained from the EPI data. For each participant, the frequency spectrum was calculated for the average of all voxels in the dataset to validate the observations (Figure 4). The most significant frequency bands were highlighted after identification of the center frequency using a local maxima approach (findpeaks implementation in Python). For all volunteers, bands with similar center frequencies of approximately 0.3Hz, 1.2Hz, and 2.4Hz can be identified. Table 1 shows the center frequency for the bands calculated for each volunteer. The frequency bands with centers at 0.3Hz and 1.2Hz closely approximate the respiration and cardiac frequencies of a human adult (around 18 breaths per minute and 72 heart beats per minute, respectively). These bands can be identified as Band 1 and Band 3 on Table 1. The averages for all volunteers are 0.322Hz and 1.217Hz respectively with a standard deviation of 0.052 and 0.078 respectively.

Figure 4:
Figure 4:

Frequency spectrum of the average signal intensity within the brain for 5 subjects. Regions of 0.3Hz around each peak frequency were highlighted to denote the bandwidth used for spatial mask creation.

View this table:
Table 1 Values of the calculated center frequency responses for each identified frequency band.

The masks created per frequency band (Figure 5) show a spatial localization for the frequency band centered at lower frequencies, e.g., 0.40Hz (Figure 5a) overlapping with brain regions with larger volume of CSF (the main ventricles and the brain periphery). Frequency bands centered at heart rate frequencies, e.g., 1.18Hz (Figure 5c) can be found on the regions with a stream of CSF (the main cerebral aqueduct). Similar patterns were observed for all volunteers as shown in Figure 6 as the mask for the heart rate band is demonstrated in each of the volunteer’s data.

Figure 5:
Figure 5:

Frequency spectrum of volunteer 3 with spatial localization of the signal from 5 separate frequency bands. For each band, a spatial mask was applied to the T1 weighted image. The acquired data can be obtained from the inferior region of the brain (cerebellum) up to the middle of the brain. The bandwidth for each band is 0.3Hz. The center frequencies are a) 0.4Hz, b) 0.76Hz, c) 1.18Hz, d) 1.96Hz, and e) 2.37Hz. The acquisition was done using an EPI sequence with TR=155ms with 19 slabs of 3 slices each for a total of 57 slices.

Figure 6:
Figure 6:

Visualization of the mask created for each volunteer at approximately the same position in the brain (bottom of the brain and top of the cerebellum and at approximately the same frequency band; a) for volunteer 1 at 1.17Hz; b) for volunteer 2 at 1.34Hz; c) for volunteer 3 at 1.18Hz; d) for volunteer 4 at 1.15Hz; and e) for volunteer 5 at 1.24Hz.

For the larger frequency spectrum (dataset with TR of 51ms), extra bands can be identified, and the center frequency of the most prominent band was calculated at 3.5Hz (Figure 7).

Figure 7:
Figure 7:

Frequency spectrum for volunteer 1 done using an EPI sequence single slice with TR=51ms. Maximum frequency of 9.8Hz and frequency bands of 0.3Hz were highlighted for better visualization of the main frequencies shown in the spectrum.

4. Discussion

We demonstrated a method to analyze the CSF motion in the human brain in-vivo using ultrafast 7T EPI acquisitions. The raw visualization of the real-time signal (Figure 1) shows in-vivo CSF motion. The flow of CSF within the ventricles and in the subarachnoid space can be clearly visualized with changes in signal intensity. The time series and the frequency spectrum comparison between the collected physiological data and the EPI data (Figure 2) shows a direct alignment between the two types of data where the cardiac and respiration cycles can be observed in the EPI MRI data. The frequency analysis also shows consistent results across multiple volunteers, with similar frequency spectrums are observed.

Compared to previous studies (Fultz et al., 2019; Shanks et al., 2019; Yamada et al., 2020; Yang et al., 2022), we were able to achieve a greater frequency range and improved frequency resolution. This was made possible by utilizing ultrafast acquisition times (ranging from 51ms to 150ms), which allowed for whole-brain spectral analysis up to 9.8Hz. Additionally, we employed high SNR and homogeneous images by using a 7T MRI with a customized RF coil system(Krishnamurthy et al., 2019; Santini et al., 2018; Santini, Wood, et al., 2021). To optimize sequence parameters, we tailored the flip-angle to maximize the signal of the CSF flow and adjusted other parameters to minimize susceptibility-related distortions. We also selected a TE that could potentially capture the BOLD signal if functional connectivity data are warranted. The creation of frequency masks allowed for an analysis of the spatial localization of each frequency band. The presence of lower frequencies responses (respiration band – 0.4Hz) in the brain periphery suggests a less turbulent fluid flow in regions with more space for fluid flow. Furthermore, the presence of heart rate frequencies (1.2Hz) in the ventricles further validates the analysis as the arterial pulse wave in the choroid plexus, for instance, is known to influence the CSF motion(Bilston et al., 2010; Iliff et al., 2013; Martin et al., 2012). Additionally, the presence of high frequencies (over 2Hz) responses can suggest a more turbulent flow that also aligns with regions of the main cerebral aqueduct.

This work provides a basis for potentially identifying new biomarkers for brain fluid dynamics. For example, the frequency spectrum can be analyzed for different brain diseases. The lower frequency bands (below 1Hz) contain physiological signals that corelate with the heart rate and breathing, so brain conditions that affect those variables can be analyzed directly from the MRI data. The magnitude of each band may also provide insights into the coupling between the heart and breathing rates with the CSF movement, which may indicate lower clearance rate. On the other hand, the higher frequency bands (above 1.8Hz), can be correlated with sleep cycles and potential sleep studies (Xie et al., 2013).

5. Conclusion

The development of non-invasive neuroimaging methods and biomarkers is essential for studying brain diseases. This work presented a novel methodology to characterize the frequency spectrum and spatial localization of CSF oscillations and movements in the human brain. The use of ultrafast EPI in conjunction with 7T human MRI and simultaneous collection of physiological data enabled the identification of primary components of CSF oscillations and their mapping spatially and temporally onto the MR image and physiological domains. The methodology showed good consistency and repeatability in-vivo, making it a promising tool for potential brain dynamics and CSF flow/clearance studies. These findings may have significant implications in the diagnosis and treatment of brain diseases, and further research is necessary to explore the potential of this methodology in clinical studies.

Data and Code Availability

The codes and dataset used in this work are available upon request.

Author Contributions

Conceptualization: TM, TS, MW, KAW, DM, JWI, HJA, TSI; Methodology: TM, TS, MW, JWI, TSI; Software: TM, TS, MW; Formal analysis: TM, TS; Investigation: TM, TS, JWI, TSI; Resources: JWI, TSI; Writing - Original Draft: TM, TS, KAW, DM, JWI, HJA, TSI; Writing - Review & Editing: TM, TS, KAW, DM, JWI, HJA, TSI; Supervision: TSI; Funding acquisition HJA, TSI.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

Declaration of Competing Interests

The authors declare no competing interest.

Acknowledgments and Funding

This work was supported by NIH R01AG063525 and R01MH111265. The author Tiago Martins was partially supported by the CAPES Foundation, Ministry of Education of Brazil, 88881.128222/2016-01. This research was supported in part by the University of Pittsburgh Center for Research Computing, RRID:SCR_022735, through the resources provided. Specifically, this work used the H2P cluster, which is supported by NSF award number OAC-2117681.

References

  1. ↵
    Adolph, R., Fukusumi, H., & Fowler, N. (1967). Origin of cerebrospinal fluid pulsations. American Journal of Physiology-Legacy Content, 212(4), 840–846. doi:10.1152/ajplegacy.1967.212.4.840
    OpenUrlCrossRef
  2. ↵
    Andersson, J. L. R., Skare, S., & Ashburner, J. (2003). How to correct susceptibility distortions in spin-echo echo-planar images: Application to diffusion tensor imaging. NeuroImage, 20(2), 870–888. doi:10.1016/S1053-8119(03)00336-7
    OpenUrlCrossRef
  3. ↵
    Avants, B. B., Tustison, N., & Song, G. (2009). Advanced normalization tools (ANTS). Insight j, 2(365), 1–35.
    OpenUrl
  4. ↵
    Bilston, L. E., Stoodley, M. A., & Fletcher, D. F. (2010). The influence of the relative timing of arterial and subarachnoid space pulse waves on spinal perivascular cerebrospinal fluid flow as a possible factor in syrinx development: Laboratory investigation. Journal of Neurosurgery, 112(4), 808–813. doi:10.3171/2009.5.JNS08945
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    Biswal, B., Yetkin, F. Z., Haughton, V. M., & Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar mri. Magnetic Resonance in Medicine, 34(4), 537–541. doi:10.1002/mrm.1910340409
    OpenUrlCrossRefPubMedWeb of Science
  6. ↵
    Dagli, M. S., Ingeholm, J. E., & Haxby, J. V. (1999). Localization of Cardiac-Induced Signal Change in fMRI. NeuroImage, 9(4), 407–415. doi:10.1006/nimg.1998.0424
    OpenUrlCrossRefPubMedWeb of Science
  7. ↵
    ECG: Cardiology | Research | BIOPAC. (n.d.). BIOPAC Systems, Inc. Retrieved April 20, 2021, from https://www.biopac.com/application/ecg-cardiology/
  8. ↵
    Eide, P. K., & Ringstad, G. (2019). Delayed clearance of cerebrospinal fluid tracer from entorhinal cortex in idiopathic normal pressure hydrocephalus: A glymphatic magnetic resonance imaging study. Journal of Cerebral Blood Flow & Metabolism, 39(7), 1355–1368. doi:10.1177/0271678X18760974
    OpenUrlCrossRefPubMed
  9. ↵
    Foi, A. (2011). Noise estimation and removal in MR imaging: The variance-stabilization approach. 2011 IEEE International Symposium on Biomedical Imaging: From Nano to Macro, 1809–1814. doi:10.1109/ISBI.2011.5872758
    OpenUrlCrossRef
  10. ↵
    Fultz, N. E., Bonmassar, G., Setsompop, K., Stickgold, R. A., Rosen, B. R., Polimeni, J. R., & Lewis, L. D. (2019). Coupled electrophysiological, hemodynamic, and cerebrospinal fluid oscillations in human sleep. Science, 366(6465), 628–631. doi:10.1126/science.aax5440
    OpenUrlAbstract/FREE Full Text
  11. ↵
    Gulani, V., Calamante, F., Shellock, F. G., Kanal, E., & Reeder, S. B. (2017). Gadolinium deposition in the brain: Summary of evidence and recommendations. The Lancet Neurology, 16(7), 564–570. doi:10.1016/S1474-4422(17)30158-8
    OpenUrlCrossRefPubMed
  12. ↵
    Hauglund, N. L., Pavan, C., & Nedergaard, M. (2020). Cleaning the sleeping brain – the potential restorative function of the glymphatic system. Current Opinion in Physiology, 15, 1–6. doi:10.1016/j.cophys.2019.10.020
    OpenUrlCrossRef
  13. ↵
    Herculano-Houzel, S. (2013). Sleep It Out. Science, 342(6156), 316–317. doi:10.1126/science.1245798
    OpenUrlAbstract/FREE Full Text
  14. ↵
    Hock, C., Golombowski, S., Müller-Spahn, F., Naser, W., Beyreuther, K., Mönning, U., Schenk, D., Vigo-Pelfrey, C., Bush, A. M., Moir, R., Tanzi, R. E., Growdon, J. H., & Nitsch, R. M. (1998). Cerebrospinal Fluid Levels of Amyloid Precursor Protein and Amyloid β-Peptide in Alzheimer’s Disease and Major Depression – Inverse Correlation with Dementia Severity. European Neurology, 39(2), 111–118. doi:10.1159/000007917
    OpenUrlCrossRefPubMedWeb of Science
  15. ↵
    Ibrahim, T. S., Hue, Y.-K., Boada, F. E., & Gilbert, R. (2008). Tic Tac Toe: Highly-Coupled, Load Insensitive Tx/Rx Array and a Quadrature Coil Without Lumped Capacitors. 16th, 438.
  16. ↵
    Ibrahim, T. S., Krishnamurthy, N., Wood, S., Raval, S., & Kim, H. (2013). 20-To-8 channel Tx array with 32-channel adjustable receive-only insert for 7T head imaging. The 21st International Society of Magnetic Resonance in Medicine Annual Meeting.
  17. ↵
    Iliff, J. J., Wang, M., Zeppenfeld, D. M., Venkataraman, A., Plog, B. A., Liao, Y., Deane, R., & Nedergaard, M. (2013). Cerebral Arterial Pulsation Drives Paravascular CSF–Interstitial Fluid Exchange in the Murine Brain. Journal of Neuroscience, 33(46), 18190–18199. doi:10.1523/JNEUROSCI.1592-13.2013
    OpenUrlAbstract/FREE Full Text
  18. ↵
    Jenkinson, M., Beckmann, C. F., Behrens, T. E. J., Woolrich, M. W., & Smith, S. M. (2012). FSL. NeuroImage, 62(2), 782–790. doi:10.1016/j.neuroimage.2011.09.015
    OpenUrlCrossRefPubMedWeb of Science
  19. ↵
    Kim, J., Krishnamurthy, N., Santini, T., Zhao, Y., Zhao, T., Bae, K. T., & Ibrahim, T. S. (2016). Experimental and numerical analysis of B1+ field and SAR with a new transmit array design for 7T breast MRI. Journal of Magnetic Resonance, 269, 55–64. doi:10.1016/j.jmr.2016.04.012
    OpenUrlCrossRef
  20. ↵
    Kiviniemi, V., Jauhiainen, J., Tervonen, O., Pääkkö, E., Oikarinen, J., Vainionpää, V., Rantala, H., & Biswal, B. (2000). Slow vasomotor fluctuation in fMRI of anesthetized child brain. Magnetic Resonance in Medicine, 44(3), 373–378. doi:10.1002/1522-2594(200009)44:3<373::AID-MRM5>3.0.CO;2-P
    OpenUrlCrossRefPubMedWeb of Science
  21. ↵
    Kiviniemi, V., Wang, X., Korhonen, V., Keinänen, T., Tuovinen, T., Autio, J., LeVan, P., Keilholz, S., Zang, Y.-F., Hennig, J., & Nedergaard, M. (2016). Ultra-fast magnetic resonance encephalography of physiological brain activity – Glymphatic pulsation mechanisms? Journal of Cerebral Blood Flow & Metabolism, 36(6), 1033–1045. doi:10.1177/0271678X15622047
    OpenUrlCrossRefPubMed
  22. ↵
    Krishnamurthy, N., Santini, T., Wood, S., Kim, J., Zhao, T., Aizenstein, H. J., & Ibrahim, T. S. (2019). Computational and experimental evaluation of the Tic-Tac-Toe RF coil for 7 Tesla MRI. PLOS ONE, 14(1), e0209663. doi:10.1371/journal.pone.0209663
    OpenUrlCrossRefPubMed
  23. ↵
    Leon, M. J. de, Li, Y., Okamura, N., Tsui, W. H., Saint-Louis, L. A., Glodzik, L., Osorio, R. S., Fortea, J., Butler, T., Pirraglia, E., Fossati, S., Kim, H.-J., Carare, R. O., Nedergaard, M., Benveniste, H., & Rusinek, H. (2017). Cerebrospinal Fluid Clearance in Alzheimer Disease Measured with Dynamic PET. Journal of Nuclear Medicine, 58(9), 1471–1476. doi:10.2967/jnumed.116.187211
    OpenUrlAbstract/FREE Full Text
  24. ↵
    Maggioni, M., Katkovnik, V., Egiazarian, K., & Foi, A. (2013). Nonlocal Transform-Domain Filter for Volumetric Data Denoising and Reconstruction. IEEE Transactions on Image Processing, 22(1), 119–133. doi:10.1109/TIP.2012.2210725
    OpenUrlCrossRefPubMed
  25. ↵
    Martin, B. A., Reymond, P., Novy, J., Balédent, O., & Stergiopulos, N. (2012). A coupled hydrodynamic model of the cardiovascular and cerebrospinal fluid system. American Journal of Physiology-Heart and Circulatory Physiology, 302(7), H1492–H1509. doi:10.1152/ajpheart.00658.2011
    OpenUrlCrossRefPubMed
  26. ↵
    MATLAB - MathWorks. (n.d.). Retrieved April 20, 2021, from https://www.mathworks.com/products/matlab.html
  27. ↵
    Moeller, S., Yacoub, E., Olman, C. A., Auerbach, E., Strupp, J., Harel, N., & Uğurbil, K. (2010). Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI. Magnetic Resonance in Medicine, 63(5), 1144–1153. doi:10.1002/mrm.22361
    OpenUrlCrossRefPubMed
  28. ↵
    Moser, E., Stahlberg, F., Ladd, M. E., & Trattnig, S. (2012). 7-T MR—from research to clinical applications? NMR in Biomedicine, 25(5), 695–716. doi:10.1002/nbm.1794
    OpenUrlCrossRef
  29. ↵
    Nedergaard, M. (2013). Garbage Truck of the Brain. Science, 340(6140), 1529–1530. doi:10.1126/science.1240514
    OpenUrlAbstract/FREE Full Text
  30. ↵
    Peng, W., Achariyar, T. M., Li, B., Liao, Y., Mestre, H., Hitomi, E., Regan, S., Kasper, T., Peng, S., Ding, F., Benveniste, H., Nedergaard, M., & Deane, R. (2016). Suppression of glymphatic fluid transport in a mouse model of Alzheimer’s disease. Neurobiology of Disease, 93, 215–225. doi:10.1016/j.nbd.2016.05.015
    OpenUrlCrossRefPubMed
  31. ↵
    Pomara, N., Bruno, D., Sarreal, A. S., Hernando, R. T., Nierenberg, J., Petkova, E., Sidtis, J. J., Wisniewski, T. M., Mehta, P. D., Pratico, D., Zetterberg, H., & Blennow, K. (2012). Lower CSF Amyloid Beta Peptides and Higher F2-Isoprostanes in Cognitively Intact Elderly Individuals With Major Depressive Disorder. American Journal of Psychiatry, 169(5), 523–530. doi:10.1176/appi.ajp.2011.11081153
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    Poncelet, B. P., Wedeen, V. J., Weisskoff, R. M., & Cohen, M. S. (1992). Brain parenchyma motion: Measurement with cine echo-planar MR imaging. Radiology, 185(3), 645–651. doi:10.1148/radiology.185.3.1438740
    OpenUrlCrossRefPubMedWeb of Science
  33. ↵
    Purdon, P. L., & Weisskoff, R. M. (1998). Effect of temporal autocorrelation due to physiological noise and stimulus paradigm on voxel-level false-positive rates in fMRI. Human Brain Mapping, 6(4), 239–249. doi:10.1002/(SICI)1097-0193(1998)6:4<239::AID-HBM4>3.0.CO;2-4
    OpenUrlCrossRefPubMedWeb of Science
  34. ↵
    Ramanathan, A., Nelson, A., Sagare, A., & Zlokovic, B. (2015). Impaired vascular-mediated clearance of brain amyloid beta in Alzheimer’s disease: The role, regulation and restoration of LRP1. Frontiers in Aging Neuroscience, 7. https://www.frontiersin.org/article/10.3389/fnagi.2015.00136
  35. ↵
    Rennels, M. L., Blaumanis, O. R., & Grady, P. A. (1990). Rapid solute transport throughout the brain via paravascular fluid pathways. Advances in Neurology, 52, 431–439.
    OpenUrlPubMed
  36. ↵
    Respiration Transducer for MRI | TSD221-MRI | Research | BIOPAC. (n.d.). BIOPAC Systems, Inc. Retrieved April 20, 2021, from https://www.biopac.com/product/respiration-transducer-for-mri/
  37. ↵
    Ringstad, G., Valnes, L. M., Dale, A. M., Pripp, A. H., Vatnehol, S.-A. S., Emblem, K. E., Mardal, K.-A., & Eide, P. K. (2018). Brain-wide glymphatic enhancement and clearance in humans assessed with MRI. JCI Insight, 3(13), e121537. doi:10.1172/jci.insight.121537
    OpenUrlCrossRefPubMed
  38. ↵
    Ringstad, G., Vatnehol, S. A. S., & Eide, P. K. (2017). Glymphatic MRI in idiopathic normal pressure hydrocephalus. Brain, 140(10), 2691–2705. doi:10.1093/brain/awx191
    OpenUrlCrossRefPubMed
  39. ↵
    Santini, T., Koo, M., Farhat, N., Campos, V. P., Alkhateeb, S., Vieira, M. A. C., Butters, M. A., Rosano, C., Aizenstein, H. J., Mettenburg, J., Novelli, E. M., & Ibrahim, T. S. (2021). Analysis of hippocampal subfields in sickle cell disease using ultrahigh field MRI. NeuroImage: Clinical, 30, 102655. doi:10.1016/j.nicl.2021.102655
    OpenUrlCrossRef
  40. ↵
    Santini, T., Wood, S., Krishnamurthy, N., Martins, T., Aizenstein, H. J., & Ibrahim, T. S. (2021). Improved 7 Tesla transmit field homogeneity with reduced electromagnetic power deposition using coupled Tic Tac Toe antennas. Scientific Reports, 11(1), Article 1. doi:10.1038/s41598-020-79807-9
    OpenUrlCrossRef
  41. ↵
    Santini, T., Zhao, Y., Wood, S., Krishnamurthy, N., Kim, J., Farhat, N., Alkhateeb, S., Martins, T., Koo, M., Zhao, T., Aizenstein, H. J., & Ibrahim, T. S. (2018). In-vivo and numerical analysis of the eigenmodes produced by a multi-level Tic-Tac-Toe head transmit array for 7 Tesla MRI. PLOS ONE, 13(11), e0206127. doi:10.1371/journal.pone.0206127
    OpenUrlCrossRefPubMed
  42. ↵
    Schroth, G., & Klose, U. (1992). Cerebrospinal fluid flow. Neuroradiology, 35(1), 10–15. doi:10.1007/BF00588271
    OpenUrlCrossRefPubMedWeb of Science
  43. ↵
    Scouten, A., & Constable, R. T. (2008). VASO-based calculations of CBV change: Accounting for the dynamic CSF volume. Magnetic Resonance in Medicine, 59(2), 308–315. doi:10.1002/mrm.21427
    OpenUrlCrossRefPubMed
  44. ↵
    Shanks, J., Bloch, K. M., Laurell, K., Cesarini, K. G., Fahlström, M., Larsson, E.-M., & Virhammar, J. (2019). Aqueductal CSF Stroke Volume Is Increased in Patients with Idiopathic Normal Pressure Hydrocephalus and Decreases after Shunt Surgery. American Journal of Neuroradiology, 40(3), 453–459. doi:10.3174/ajnr.A5972
    OpenUrlAbstract/FREE Full Text
  45. ↵
    Tustison, N. J., Avants, B. B., Cook, P. A., Zheng, Y., Egan, A., Yushkevich, P. A., & Gee, J. C. (2010). N4ITK: Improved N3 bias correction. IEEE Transactions on Medical Imaging, 29(6), 1310–1320. doi:10.1109/TMI.2010.2046908
    OpenUrlCrossRefPubMedWeb of Science
  46. ↵
    Uğurbil, K., Xu, J., Auerbach, E. J., Moeller, S., Vu, A. T., Duarte-Carvajalino, J. M., Lenglet, C., Wu, X., Schmitter, S., Van de Moortele, P. F., Strupp, J., Sapiro, G., De Martino, F., Wang, D., Harel, N., Garwood, M., Chen, L., Feinberg, D. A., Smith, S. M., … Yacoub, E. (2013). Pushing spatial and temporal resolution for functional and diffusion MRI in the Human Connectome Project. NeuroImage, 80, 80–104. doi:10.1016/j.neuroimage.2013.05.012
    OpenUrlCrossRefPubMedWeb of Science
  47. ↵
    Xie, L., Kang, H., Xu, Q., Chen, M. J., Liao, Y., Thiyagarajan, M., O’Donnell, J., Christensen, D. J., Nicholson, C., Iliff, J. J., Takano, T., Deane, R., & Nedergaard, M. (2013). Sleep Drives Metabolite Clearance from the Adult Brain. Science, 342(6156), 373–377. doi:10.1126/science.1241224
    OpenUrlAbstract/FREE Full Text
  48. ↵
    Yamada, S., Ishikawa, M., Ito, H., Yamamoto, K., Yamaguchi, M., Oshima, M., & Nozaki, K. (2020). Cerebrospinal fluid dynamics in idiopathic normal pressure hydrocephalus on four-dimensional flow imaging. European Radiology, 30(8), 4454–4465. doi:10.1007/s00330-020-06825-6
    OpenUrlCrossRefPubMed
  49. ↵
    Yang, H.-C. (Shawn), Inglis, B., Talavage, T. M., Nair, V. V., Yao, J. (Fiona), Fitzgerald, B., Schwichtenberg, A. J., & Tong, Y. (2022). Coupling between cerebrovascular oscillations and CSF flow fluctuations during wakefulness: An fMRI study. Journal of Cerebral Blood Flow & Metabolism, 0271678X221074639. doi:10.1177/0271678X221074639
    OpenUrlCrossRef
Back to top
Posted December 06, 2023.
Download PDF
Data/Code
Email
Characterization of oscillations in the brain and cerebrospinal fluid using ultra-high field magnetic resonance imaging
Tiago Martins, Tales Santini, Bruno de Almeida, Minjie Wu, Kristine A. Wilckens, Davneet Minhas, James W. Ibinson, Howard J. Aizenstein, Tamer S. Ibrahim
medRxiv 2023.12.05.23299452; doi: https://doi.org/10.1101/2023.12.05.23299452
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area