1
|
Goodman SS, Haysley S, Jennings SG. Human Olivocochlear Effects: A Statistical Detection Approach Applied to the Cochlear Microphonic Evoked by Swept Tones. J Assoc Res Otolaryngol 2024:10.1007/s10162-024-00956-z. [PMID: 38954166 DOI: 10.1007/s10162-024-00956-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Accepted: 06/12/2024] [Indexed: 07/04/2024] Open
Abstract
The human medial olivocochlear (MOC) reflex was assessed by observing the effects of contralateral acoustic stimulation (CAS) on the cochlear microphonic (CM) across a range of probe frequencies. A frequency-swept probe tone (125-4757 Hz, 90 dB SPL) was presented in two directions (up sweep and down sweep) to normal-hearing young adults. This study assessed MOC effects on the CM in individual participants using a statistical approach that calculated minimum detectable changes in magnitude and phase based on CM signal-to-noise ratio (SNR). Significant increases in CM magnitude, typically 1-2 dB in size, were observed for most participants from 354 to 1414 Hz, where the size and consistency of these effects depended on participant, probe frequency, sweep direction, and SNR. CAS-related phase lags were also observed, consistent with CM-based MOC studies in laboratory animals. Observed effects on CM magnitude and phase were in the opposite directions of reported effects on otoacoustic emissions (OAEs). OAEs are sensitive to changes in the motility of outer hair cells located near the peak region of the traveling wave, while the effects of CAS on the CM likely originate from MOC-related changes in the conductance of outer hair cells located in the basal tail of the traveling wave. Thus, MOC effects on the CM are complementary to those observed for OAEs.
Collapse
Affiliation(s)
- Shawn S Goodman
- Department of Communication Sciences and Disorders, University of Iowa, Iowa City, IA, USA
| | - Sarah Haysley
- Department of Communication Sciences and Disorders, University of Utah, Salt Lake City, UT, USA
| | - Skyler G Jennings
- Department of Communication Sciences and Disorders, University of Utah, Salt Lake City, UT, USA.
| |
Collapse
|
2
|
Li S, Wang F, Gao S. New non-local mean methods for MRI denoising based on global self-similarity between values. Comput Biol Med 2024; 174:108450. [PMID: 38608325 DOI: 10.1016/j.compbiomed.2024.108450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Revised: 03/20/2024] [Accepted: 04/07/2024] [Indexed: 04/14/2024]
Abstract
Magnetic resonance imaging (MRI) is a non-invasive medical imaging technique that provides high-resolution 3D images and valuable insights into human tissue conditions. Even at present, the refinement of denoising methods for MRI remains a crucial concern for improving the quality of the images. This study aims to improve the prefiltered rotationally invariant non-local principal component analysis (PRI-NL-PCA) algorithm. We relaxed the original restrictions using particle swarm optimization to determine optimal parameters for the PCA part of the original algorithm. In addition, we adjusted the prefiltered rotationally invariant non-local mean (PRI-NLM) part by traversing the signal intensities of voxels instead of their spatial positions to reduce duplicate calculations and expand the search volume to the whole image when estimating voxels' signal intensities. The new method demonstrated superior denoising performance compared to the original approach. Moreover, in most cases, the new algorithm ran faster. Furthermore, our proposed method can also be applied to process Gaussian noise in natural images and has the potential to enhance other NLM-based denoising algorithms.
Collapse
Affiliation(s)
- Shiao Li
- Institute of Medical Technology, Peking University Health Science Center, Haidian District College Road No. 38, 100191, Beijing, China.
| | - Fei Wang
- Key Laboratory of Carcinogenesis and Translational Research, Department of Radiation Oncology, Beijing Cancer Hospital, Haidian District Fucheng Road No. 52, 100142, Beijing, China.
| | - Song Gao
- Institute of Medical Technology, Peking University Health Science Center, Haidian District College Road No. 38, 100191, Beijing, China.
| |
Collapse
|
3
|
Pineda Guzman RA, Naughton N, Majumdar S, Damon B, Kersh ME. Assessment of Mechanically Induced Changes in Helical Fiber Microstructure Using Diffusion Tensor Imaging. Ann Biomed Eng 2024; 52:832-844. [PMID: 38151645 DOI: 10.1007/s10439-023-03420-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 12/04/2023] [Indexed: 12/29/2023]
Abstract
Noninvasive methods to detect microstructural changes in collagen-based fibrous tissues are necessary to differentiate healthy from damaged tissues in vivo but are sparse. Diffusion Tensor Imaging (DTI) is a noninvasive imaging technique used to quantitatively infer tissue microstructure with previous work primarily focused in neuroimaging applications. Yet, it is still unclear how DTI metrics relate to fiber microstructure and function in musculoskeletal tissues such as ligament and tendon, in part because of the high heterogeneity inherent to such tissues. To address this limitation, we assessed the ability of DTI to detect microstructural changes caused by mechanical loading in tissue-mimicking helical fiber constructs of known structure. Using high-resolution optical and micro-computed tomography imaging, we found that static and fatigue loading resulted in decreased sample diameter and a re-alignment of the macro-scale fiber twist angle similar with the direction of loading. However, DTI and micro-computed tomography measurements suggest microstructural differences in the effect of static versus fatigue loading that were not apparent at the bulk level. Specifically, static load resulted in an increase in diffusion anisotropy and a decrease in radial diffusivity suggesting radially uniform fiber compaction. In contrast, fatigue loads resulted in increased diffusivity in all directions and a change in the alignment of the principal diffusion direction away from the constructs' main axis suggesting fiber compaction and microstructural disruptions in fiber architecture. These results provide quantitative evidence of the ability of DTI to detect mechanically induced changes in tissue microstructure that are not apparent at the bulk level, thus confirming its potential as a noninvasive measure of microstructure in helically architected collagen-based tissues, such as ligaments and tendons.
Collapse
Affiliation(s)
| | - Noel Naughton
- Beckman Institute for Advanced Science & Technology, University of Illinois Urbana-Champaign, Urbana, IL, USA
| | - Shreyan Majumdar
- Beckman Institute for Advanced Science & Technology, University of Illinois Urbana-Champaign, Urbana, IL, USA
| | - Bruce Damon
- Beckman Institute for Advanced Science & Technology, University of Illinois Urbana-Champaign, Urbana, IL, USA
- Carle Clinical Imaging Research Program, Stephens Family Clinical Research Institute, Carle Health, Urbana, IL, USA
- Department of Bioengineering, University of Illinois Urbana-Champaign, Urbana, IL, USA
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Radiology and Radiological Science, Vanderbilt University, Nashville, TN, USA
- Carle Illinois College of Medicine, University of Illinois Urbana-Champaign, Urbana, IL, USA
| | - Mariana E Kersh
- Department of Mechanical Science & Engineering, University of Illinois Urbana-Champaign, Urbana, IL, USA.
- Beckman Institute for Advanced Science & Technology, University of Illinois Urbana-Champaign, Urbana, IL, USA.
- Carle Illinois College of Medicine, University of Illinois Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
4
|
Lee HH, Tian Q, Sheft M, Coronado-Leija R, Ramos-Llorden G, Abdollahzadeh A, Fieremans E, Novikov DS, Huang SY. The effects of axonal beading and undulation on axonal diameter estimation from diffusion MRI: Insights from simulations in human axons segmented from three-dimensional electron microscopy. NMR IN BIOMEDICINE 2024; 37:e5087. [PMID: 38168082 PMCID: PMC10942763 DOI: 10.1002/nbm.5087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 11/19/2023] [Accepted: 11/21/2023] [Indexed: 01/05/2024]
Abstract
The increasing availability of high-performance gradient systems in human MRI scanners has generated great interest in diffusion microstructural imaging applications such as axonal diameter mapping. Practically, sensitivity to axon diameter in diffusion MRI is attained at strong diffusion weightings b , where the deviation from the expected 1 / b scaling in white matter yields a finite transverse diffusivity, which is then translated into an axon diameter estimate. While axons are usually modeled as perfectly straight, impermeable cylinders, local variations in diameter (caliber variation or beading) and direction (undulation) are known to influence axonal diameter estimates and have been observed in microscopy data of human axons. In this study, we performed Monte Carlo simulations of diffusion in axons reconstructed from three-dimensional electron microscopy of a human temporal lobe specimen using simulated sequence parameters matched to the maximal gradient strength of the next-generation Connectome 2.0 human MRI scanner ( ≲ 500 mT/m). We show that axon diameter estimation is accurate for nonbeaded, nonundulating fibers; however, in fibers with caliber variations and undulations, the axon diameter is heavily underestimated due to caliber variations, and this effect overshadows the known overestimation of the axon diameter due to undulations. This unexpected underestimation may originate from variations in the coarse-grained axial diffusivity due to caliber variations. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | - Maxina Sheft
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard–MIT Health Sciences and Technology, Cambridge, Massachusetts, USA
| | - Ricardo Coronado-Leija
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Gabriel Ramos-Llorden
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| | - Ali Abdollahzadeh
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Dmitry S. Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, New York, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, New York, USA
| | - Susie Y. Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, USA
- Harvard Medical School, Boston, Massachusetts, USA
| |
Collapse
|
5
|
Fouto AR, Henriques RN, Golub M, Freitas AC, Ruiz-Tagle A, Esteves I, Gil-Gouveia R, Silva NA, Vilela P, Figueiredo P, Nunes RG. Impact of truncating diffusion MRI scans on diffusional kurtosis imaging. MAGMA (NEW YORK, N.Y.) 2024:10.1007/s10334-024-01153-y. [PMID: 38393541 DOI: 10.1007/s10334-024-01153-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Revised: 01/09/2024] [Accepted: 02/07/2024] [Indexed: 02/25/2024]
Abstract
OBJECTIVE Diffusional kurtosis imaging (DKI) extends diffusion tensor imaging (DTI), characterizing non-Gaussian diffusion effects but requires longer acquisition times. To ensure the robustness of DKI parameters, data acquisition ordering should be optimized allowing for scan interruptions or shortening. Three methodologies were used to examine how reduced diffusion MRI scans impact DKI histogram-metrics: 1) the electrostatic repulsion model (OptEEM); 2) spherical codes (OptSC); 3) random (RandomTRUNC). MATERIALS AND METHODS Pre-acquired diffusion multi-shell data from 14 female healthy volunteers (29±5 years) were used to generate reordered data. For each strategy, subsets containing different amounts of the full dataset were generated. The subsampling effects were assessed on histogram-based DKI metrics from tract-based spatial statistics (TBSS) skeletonized maps. To evaluate each subsampling method on simulated data at different SNRs and the influence of subsampling on in vivo data, we used a 3-way and 2-way repeated measures ANOVA, respectively. RESULTS Simulations showed that subsampling had different effects depending on DKI parameter, with fractional anisotropy the most stable (up to 5% error) and radial kurtosis the least stable (up to 26% error). RandomTRUNC performed the worst while the others showed comparable results. Furthermore, the impact of subsampling varied across distinct histogram characteristics, the peak value the least affected (OptEEM: up to 5% error; OptSC: up to 7% error) and peak height (OptEEM: up to 8% error; OptSC: up to 11% error) the most affected. CONCLUSION The impact of truncation depends on specific histogram-based DKI metrics. The use of a strategy for optimizing the acquisition order is advisable to improve DKI robustness to exam interruptions.
Collapse
Affiliation(s)
- Ana R Fouto
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal.
| | | | - Marc Golub
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Andreia C Freitas
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Amparo Ruiz-Tagle
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Inês Esteves
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Raquel Gil-Gouveia
- Neurology Department, Hospital da Luz, Lisbon, Portugal
- Center for Interdisciplinary Research in Health, Universidade Católica Portuguesa, Lisbon, Portugal
| | - Nuno A Silva
- Learning Health, Hospital da Luz, Lisbon, Portugal
| | - Pedro Vilela
- Imaging Department, Hospital da Luz, Lisbon, Portugal
| | - Patrícia Figueiredo
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| | - Rita G Nunes
- Institute for Systems and Robotics-Lisboa and Department of Bioengineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal
| |
Collapse
|
6
|
Cerdán Cerdá A, Toschi N, Treaba CA, Barletta V, Herranz E, Mehndiratta A, Gomez-Sanchez JA, Mainero C, De Santis S. A translational MRI approach to validate acute axonal damage detection as an early event in multiple sclerosis. eLife 2024; 13:e79169. [PMID: 38192199 PMCID: PMC10776086 DOI: 10.7554/elife.79169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Accepted: 12/05/2023] [Indexed: 01/10/2024] Open
Abstract
Axonal degeneration is a central pathological feature of multiple sclerosis and is closely associated with irreversible clinical disability. Current noninvasive methods to detect axonal damage in vivo are limited in their specificity and clinical applicability, and by the lack of proper validation. We aimed to validate an MRI framework based on multicompartment modeling of the diffusion signal (AxCaliber) in rats in the presence of axonal pathology, achieved through injection of a neurotoxin damaging the neuronal terminal of axons. We then applied the same MRI protocol to map axonal integrity in the brain of multiple sclerosis relapsing-remitting patients and age-matched healthy controls. AxCaliber is sensitive to acute axonal damage in rats, as demonstrated by a significant increase in the mean axonal caliber along the targeted tract, which correlated with neurofilament staining. Electron microscopy confirmed that increased mean axonal diameter is associated with acute axonal pathology. In humans with multiple sclerosis, we uncovered a diffuse increase in mean axonal caliber in most areas of the normal-appearing white matter, preferentially affecting patients with short disease duration. Our results demonstrate that MRI-based axonal diameter mapping is a sensitive and specific imaging biomarker that links noninvasive imaging contrasts with the underlying biological substrate, uncovering generalized axonal damage in multiple sclerosis as an early event.
Collapse
Affiliation(s)
| | - Nicola Toschi
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical SchoolBostonUnited States
- Department of Biomedicine and Prevention, University of Rome Tor VergataRomeItaly
| | - Constantina A Treaba
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical SchoolBostonUnited States
| | - Valeria Barletta
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical SchoolBostonUnited States
| | - Elena Herranz
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical SchoolBostonUnited States
| | - Ambica Mehndiratta
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical SchoolBostonUnited States
| | - Jose A Gomez-Sanchez
- Instituto de Neurociencias de Alicante, CSIC-UMHSan Juan de AlicanteSpain
- Instituto de Investigación Sanitaria y Biomédica de Alicante (ISABIAL)AlicanteSpain
- Millennium Nucleus for the Study of Pain (MiNuSPain)SantiagoChile
| | - Caterina Mainero
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Harvard Medical SchoolBostonUnited States
| | - Silvia De Santis
- Instituto de Neurociencias de Alicante, CSIC-UMHSan Juan de AlicanteSpain
| |
Collapse
|
7
|
Manzano-Patron JP, Moeller S, Andersson JLR, Ugurbil K, Yacoub E, Sotiropoulos SN. DENOISING DIFFUSION MRI: CONSIDERATIONS AND IMPLICATIONS FOR ANALYSIS. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.24.550348. [PMID: 37546835 PMCID: PMC10402048 DOI: 10.1101/2023.07.24.550348] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Development of diffusion MRI (dMRI) denoising approaches has experienced considerable growth over the last years. As noise can inherently reduce accuracy and precision in measurements, its effects have been well characterised both in terms of uncertainty increase in dMRI-derived features and in terms of biases caused by the noise floor, the smallest measurable signal given the noise level. However, gaps in our knowledge still exist in objectively characterising dMRI denoising approaches in terms of both of these effects and assessing their efficacy. In this work, we reconsider what a denoising method should and should not do and we accordingly define criteria to characterise the performance. We propose a comprehensive set of evaluations, including i) benefits in improving signal quality and reducing noise variance, ii) gains in reducing biases and the noise floor and improving, iii) preservation of spatial resolution, iv) agreement of denoised data against a gold standard, v) gains in downstream parameter estimation (precision and accuracy), vi) efficacy in enabling noise-prone applications, such as ultra-high-resolution imaging. We further provide newly acquired complex datasets (magnitude and phase) with multiple repeats that sample different SNR regimes to highlight performance differences under different scenarios. Without loss of generality, we subsequently apply a number of exemplar patch-based denoising algorithms to these datasets, including Non-Local Means, Marchenko-Pastur PCA (MPPCA) in the magnitude and complex domain and NORDIC, and compare them with respect to the above criteria and against a gold standard complex average of multiple repeats. We demonstrate that all tested denoising approaches reduce noise-related variance, but not always biases from the elevated noise floor. They all induce a spatial resolution penalty, but its extent can vary depending on the method and the implementation. Some denoising approaches agree with the gold standard more than others and we demonstrate challenges in even defining such a standard. Overall, we show that dMRI denoising performed in the complex domain is advantageous to magnitude domain denoising with respect to all the above criteria.
Collapse
Affiliation(s)
| | - Steen Moeller
- Center for Magnetic Resonance Research, University of Minnesota, USA
| | | | - Kamil Ugurbil
- Center for Magnetic Resonance Research, University of Minnesota, USA
| | - Essa Yacoub
- Center for Magnetic Resonance Research, University of Minnesota, USA
| | - Stamatios N Sotiropoulos
- Sir Peter Mansfield Imaging Centre, School of Medicine, University of Nottingham, UK
- Wellcome Centre for Integrative Neuroimaging, University of Oxford, UK
- Nottingham Biomedical Research Centre, Queen's Medical Centre, University of Nottingham, UK
| |
Collapse
|
8
|
Dai E, Zhu A, Yang GK, Quah K, Tan ET, Fiveland E, Foo TKF, McNab JA. Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient. Neuroimage 2023; 279:120328. [PMID: 37586445 PMCID: PMC10529993 DOI: 10.1016/j.neuroimage.2023.120328] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 07/17/2023] [Accepted: 08/12/2023] [Indexed: 08/18/2023] Open
Abstract
Measuring the time/frequency dependence of diffusion MRI is a promising approach to distinguish between the effects of different tissue microenvironments, such as membrane restriction, tissue heterogeneity, and compartmental water exchange. In this study, we measure the frequency dependence of diffusivity (D) and kurtosis (K) with oscillating gradient diffusion encoding waveforms and a diffusion kurtosis imaging (DKI) model in human brains using a high-performance, head-only MAGNUS gradient system, with a combination of b-values, oscillating frequencies (f), and echo time that has not been achieved in human studies before. Frequency dependence of diffusivity and kurtosis are observed in both global and local white matter (WM) and gray matter (GM) regions and characterized with a power-law model ∼Λ*fθ. The frequency dependences of diffusivity and kurtosis (including changes between fmin and fmax, Λ, and θ) vary over different WM and GM regions, indicating potential microstructural differences between regions. A trend of decreasing kurtosis over frequency in the short-time limit is successfully captured for in vivo human brains. The effects of gradient nonlinearity (GNL) on frequency-dependent diffusivity and kurtosis measurements are investigated and corrected. Our results show that the GNL has prominent scaling effects on the measured diffusivity values (3.5∼5.5% difference in the global WM and 6∼8% difference in the global cortex) and subsequently affects the corresponding power-law parameters (Λ, θ) while having a marginal influence on the measured kurtosis values (<0.05% difference) and power-law parameters (Λ, θ). This study expands previous OGSE studies and further demonstrates the translatability of frequency-dependent diffusivity and kurtosis measurements to human brains, which may provide new opportunities to probe human brain microstructure in health and disease.
Collapse
Affiliation(s)
- Erpeng Dai
- Department of Radiology, Stanford University, Stanford, CA, USA.
| | | | - Grant K Yang
- Department of Radiology, Stanford University, Stanford, CA, USA; Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Kristin Quah
- Department of Radiology, Stanford University, Stanford, CA, USA; Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Ek T Tan
- Department of Radiology and Imaging, Hospital for Special Surgery, New York, NY, USA
| | | | | | | |
Collapse
|
9
|
Schiavi S, Palombo M, Zacà D, Tazza F, Lapucci C, Castellan L, Costagli M, Inglese M. Mapping tissue microstructure across the human brain on a clinical scanner with soma and neurite density image metrics. Hum Brain Mapp 2023; 44:4792-4811. [PMID: 37461286 PMCID: PMC10400787 DOI: 10.1002/hbm.26416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Revised: 05/02/2023] [Accepted: 06/23/2023] [Indexed: 08/05/2023] Open
Abstract
Soma and neurite density image (SANDI) is an advanced diffusion magnetic resonance imaging biophysical signal model devised to probe in vivo microstructural information in the gray matter (GM). This model requires acquisitions that include b values that are at least six times higher than those used in clinical practice. Such high b values are required to disentangle the signal contribution of water diffusing in soma from that diffusing in neurites and extracellular space, while keeping the diffusion time as short as possible to minimize potential bias due to water exchange. These requirements have limited the use of SANDI only to preclinical or cutting-edge human scanners. Here, we investigate the potential impact of neglecting water exchange in the SANDI model and present a 10-min acquisition protocol that enables to characterize both GM and white matter (WM) on 3 T scanners. We implemented analytical simulations to (i) evaluate the stability of the fitting of SANDI parameters when diminishing the number of shells; (ii) estimate the bias due to potential exchange between neurites and extracellular space in such reduced acquisition scheme, comparing it with the bias due to experimental noise. Then, we demonstrated the feasibility and assessed the repeatability and reproducibility of our approach by computing microstructural metrics of SANDI with AMICO toolbox and other state-of-the-art models on five healthy subjects. Finally, we applied our protocol to five multiple sclerosis patients. Results suggest that SANDI is a practical method to characterize WM and GM tissues in vivo on performant clinical scanners.
Collapse
Affiliation(s)
- Simona Schiavi
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI)University of GenoaGenoaItaly
| | - Marco Palombo
- CUBRIC, School of PsychologyCardiff UniversityCardiffUK
- School of Computer Science and InformaticsCardiff UniversityCardiffUK
| | | | - Francesco Tazza
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI)University of GenoaGenoaItaly
| | - Caterina Lapucci
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI)University of GenoaGenoaItaly
- HNSR, IRRCS Ospedale Policlinico San MartinoGenoaItaly
| | - Lucio Castellan
- Department of NeuroradiologyIRCCS Ospedale Policlinico San MartinoGenoaItaly
| | - Mauro Costagli
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI)University of GenoaGenoaItaly
- Laboratory of Medical Physics and Magnetic ResonanceIRCCS Stella MarisPisaItaly
| | - Matilde Inglese
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI)University of GenoaGenoaItaly
- IRCCS Ospedale Policlinico San MartinoGenoaItaly
| |
Collapse
|
10
|
Lee HH, Tian Q, Sheft M, Coronado-Leija R, Ramos-Llorden G, Abdollahzadeh A, Fieremans E, Novikov DS, Huang SY. The influence of axonal beading and undulation on axonal diameter mapping. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.19.537494. [PMID: 37131702 PMCID: PMC10153226 DOI: 10.1101/2023.04.19.537494] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
We consider the effect of non-cylindrical axonal shape on axonal diameter mapping with diffusion MRI. Practical sensitivity to axon diameter is attained at strong diffusion weightings b , where the deviation from the 1 / b scaling yields the finite transverse diffusivity, which is then translated into axon diameter. While axons are usually modeled as perfectly straight, impermeable cylinders, the local variations in diameter (caliber variation or beading) and direction (undulation) have been observed in microscopy data of human axons. Here we quantify the influence of cellular-level features such as caliber variation and undulation on axon diameter estimation. For that, we simulate the diffusion MRI signal in realistic axons segmented from 3-dimensional electron microscopy of a human brain sample. We then create artificial fibers with the same features and tune the amplitude of their caliber variations and undulations. Numerical simulations of diffusion in fibers with such tunable features show that caliber variations and undulations result in under- and over-estimation of axon diameters, correspondingly; this bias can be as large as 100%. Given that increased axonal beading and undulations have been observed in pathological tissues, such as traumatic brain injury and ischemia, the interpretation of axon diameter alterations in pathology may be significantly confounded.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| | - Maxina Sheft
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard-MIT Health Sciences and Technology, Cambridge, MA 02139, USA
| | - Ricardo Coronado-Leija
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Gabriel Ramos-Llorden
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| | - Ali Abdollahzadeh
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Dmitry S. Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY 10016, USA
- Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY 10016, USA
| | - Susie Y. Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA 02129,USA
- Harvard Medical School, Boston, MA 02115, USA
| |
Collapse
|
11
|
Huang H, Liu B, Xu Y, Zhou W. Synthetic-to-real domain adaptation with deep learning for fitting the intravoxel incoherent motion model of diffusion-weighted imaging. Med Phys 2023; 50:1614-1622. [PMID: 36308503 DOI: 10.1002/mp.16031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 10/03/2022] [Accepted: 10/03/2022] [Indexed: 01/15/2023] Open
Abstract
BACKGROUND Intravoxel incoherent motion (IVIM) is a type of diffusion-weighted imaging (DWI), and IVIM model parameters (water molecule diffusion rate Dt , pseudo-diffusion coefficient Dp , and tissue perfusion fraction Fp ) have been widely used in the diagnosis and characterization of malignant lesions. PURPOSE This study proposes a deep-learning model with synthetic-to-real domain adaptation to fit the IVIM model parameters of DWI. METHODS Ninety-eight consecutive patients diagnosed with hepatocellular carcinoma between January 2017 and September 2020 were included in the study, and routine IVIM-DWI serial examinations were performed using a 3.0 T magnetic resonance imaging system in preoperative MR imaging. The proposed method is mainly composed of two modules: a convolutional neural network-based IVIM model fitting network to map b-value images to the IVIM parameter maps and a domain discriminator to improve the accuracy of the IVIM parameter maps in the real data. The proposed method was compared with previously reported fitting methods, including the nonlinear least squares (NLSs), IVIM-NEToptim , and self-supervised U-network methods. The IVIM parameter-fitting performance was assessed by measuring the DWI reconstruction performance and testing the robustness of each method against noise using noise-corrupted data. RESULTS The DWI reconstruction performance demonstrates that the proposed method has better reconstruction accuracy for DWI with a low signal-to-noise ratio, which implies that the proposed method improves the fitting accuracy of the IVIM parameters. Noise-corrupt experiments show that the proposed method is more robust against noise-corrupted signals. With the proposed method, no outliers were found in Dt , and outliers were reduced for Fp in the abnormal regions (proposed method: 1.85%; NLS: 5.90%; IVIM-NEToptim : 6.61%; and self-U-net: 25.36%). Moreover, experiments show that the proposed method has a more stable parameter estimation performance than the existing methods in the absence of real data. CONCLUSIONS IVIM parameters can be estimated using a synthetic-to-real domain-adaptation framework with deep learning, and the proposed method outperforms previously reported methods.
Collapse
Affiliation(s)
- Haoyuan Huang
- School of Medical Information Engineering, Guangzhou University of Chinese Medicine, Guangzhou, China
| | - Baoer Liu
- Department of Medical Imaging Center, Nanfang Hospital, Southern Medical University, Guangzhou, China
| | - Yikai Xu
- Department of Medical Imaging Center, Nanfang Hospital, Southern Medical University, Guangzhou, China
| | - Wu Zhou
- School of Medical Information Engineering, Guangzhou University of Chinese Medicine, Guangzhou, China
| |
Collapse
|
12
|
Abdolalizadeh A, Ohadi MAD, Ershadi ASB, Aarabi MH. Graph theoretical approach to brain remodeling in multiple sclerosis. Netw Neurosci 2023; 7:148-159. [PMID: 37334009 PMCID: PMC10270718 DOI: 10.1162/netn_a_00276] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2021] [Accepted: 09/05/2022] [Indexed: 03/21/2024] Open
Abstract
Multiple sclerosis (MS) is a neuroinflammatory disorder damaging structural connectivity. Natural remodeling processes of the nervous system can, to some extent, restore the damage caused. However, there is a lack of biomarkers to evaluate remodeling in MS. Our objective is to evaluate graph theory metrics (especially modularity) as a biomarker of remodeling and cognition in MS. We recruited 60 relapsing-remitting MS and 26 healthy controls. Structural and diffusion MRI, plus cognitive and disability evaluations, were done. We calculated modularity and global efficiency from the tractography-derived connectivity matrices. Association of graph metrics with T2 lesion load, cognition, and disability was evaluated using general linear models adjusting for age, gender, and disease duration wherever applicable. We showed that MS subjects had higher modularity and lower global efficiency compared with controls. In the MS group, modularity was inversely associated with cognitive performance but positively associated with T2 lesion load. Our results indicate that modularity increase is due to the disruption of intermodular connections in MS because of the lesions, with no improvement or preserving of cognitive functions.
Collapse
Affiliation(s)
- AmirHussein Abdolalizadeh
- Students’ Scientific Research Program, Tehran University of Medical Sciences, Tehran, Iran
- Interdisciplinary Neuroscience Research Program, Tehran University of Medical Sciences, Tehran, Iran
| | - Mohammad Amin Dabbagh Ohadi
- Students’ Scientific Research Program, Tehran University of Medical Sciences, Tehran, Iran
- Interdisciplinary Neuroscience Research Program, Tehran University of Medical Sciences, Tehran, Iran
| | - Amir Sasan Bayani Ershadi
- Students’ Scientific Research Program, Tehran University of Medical Sciences, Tehran, Iran
- Interdisciplinary Neuroscience Research Program, Tehran University of Medical Sciences, Tehran, Iran
| | - Mohammad Hadi Aarabi
- Department of Neuroscience, Padova Neuroscience Center, University of Padova, Padova, Italy
| |
Collapse
|
13
|
Howard AF, Cottaar M, Drakesmith M, Fan Q, Huang SY, Jones DK, Lange FJ, Mollink J, Rudrapatna SU, Tian Q, Miller KL, Jbabdi S. Estimating axial diffusivity in the NODDI model. Neuroimage 2022; 262:119535. [PMID: 35931306 PMCID: PMC9802007 DOI: 10.1016/j.neuroimage.2022.119535] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2022] [Revised: 07/20/2022] [Accepted: 08/01/2022] [Indexed: 01/03/2023] Open
Abstract
To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d∥=1.7μm2/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5μm2/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.
Collapse
Affiliation(s)
- Amy Fd Howard
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom.
| | - Michiel Cottaar
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Mark Drakesmith
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom; Neuroscience and Mental Health Research Institute, Cardiff University, Cardiff, United Kingdom
| | - Qiuyun Fan
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States; Harvard Medical School, Boston, Massachusetts, United States; Department of Biomedical Engineering, College of Precision Instruments and Optoelectronics Engineering, Tianjin University, Tianjin, China; Academy of Medical Engineering and Translational Medicine, Medical College, Tianjin University, Tianjin, China
| | - Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States; Harvard Medical School, Boston, Massachusetts, United States; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom; Neuroscience and Mental Health Research Institute, Cardiff University, Cardiff, United Kingdom
| | - Frederik J Lange
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Jeroen Mollink
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Suryanarayana Umesh Rudrapatna
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, United Kingdom; Philips Innovation Campus, Bangalore, India
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States; Harvard Medical School, Boston, Massachusetts, United States
| | - Karla L Miller
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| | - Saad Jbabdi
- FMRIB Centre, Wellcome Centre for Integrative Neuroimaging, Nuffield Department of Clinical Neurosciences, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
14
|
Huang HM. An unsupervised convolutional neural network method for estimation of intravoxel incoherent motion parameters. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac9a1f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2022] [Accepted: 10/13/2022] [Indexed: 11/07/2022]
Abstract
Abstract
Objective. Intravoxel incoherent motion (IVIM) imaging obtained by fitting a biexponential model to multiple b-value diffusion-weighted magnetic resonance imaging (DW-MRI) has been shown to be a promising tool for different clinical applications. Recently, several deep neural network (DNN) methods were proposed to generate IVIM imaging. Approach. In this study, we proposed an unsupervised convolutional neural network (CNN) method for estimation of IVIM parameters. We used both simulated and real abdominal DW-MRI data to evaluate the performance of the proposed CNN-based method, and compared the results with those obtained from a non-linear least-squares fit (TRR, trust-region reflective algorithm) and a feed-forward backward-propagation DNN-based method. Main results. The simulation results showed that both the DNN- and CNN-based methods had lower coefficients of variation than the TRR method, but the CNN-based method provided more accurate parameter estimates. The results obtained from real DW-MRI data showed that the TRR method produced many biased IVIM parameter estimates that hit the upper and lower parameter bounds. In contrast, both the DNN- and CNN-based methods yielded less biased IVIM parameter estimates. Overall, the perfusion fraction and diffusion coefficient obtained from the DNN- and CNN-based methods were close to literature values. However, compared with the CNN-based method, both the TRR and DNN-based methods tended to yield increased pseudodiffusion coefficients (55%–180%). Significance. Our preliminary results suggest that it is feasible to estimate IVIM parameters using CNN.
Collapse
|
15
|
Powell E, Schneider T, Battiston M, Grussu F, Toosy A, Clayden JD, Wheeler‐Kingshott CAMG. SENSE EPI reconstruction with 2D phase error correction and channel-wise noise removal. Magn Reson Med 2022; 88:2157-2166. [PMID: 35877787 PMCID: PMC9545987 DOI: 10.1002/mrm.29349] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Revised: 05/16/2022] [Accepted: 05/16/2022] [Indexed: 11/23/2022]
Abstract
PURPOSE To develop a robust reconstruction pipeline for EPI data that enables 2D Nyquist phase error correction using sensitivity encoding without incurring major noise artifacts in low SNR data. METHODS SENSE with 2D phase error correction (PEC-SENSE) was combined with channel-wise noise removal using Marcenko-Pastur principal component analysis (MPPCA) to simultaneously eliminate Nyquist ghost artifacts in EPI data and mitigate the noise amplification associated with phase correction using parallel imaging. The proposed pipeline (coined SPECTRE) was validated in phantom DW-EPI data using the accuracy and precision of diffusion metrics; ground truth values were obtained from data acquired with a spin echo readout. Results from the SPECTRE pipeline were compared against PEC-SENSE reconstructions with three alternate denoising strategies: (i) no denoising; (ii) denoising of magnitude data after image formation; (iii) denoising of complex data after image formation. SPECTRE was then tested using highb $$ b $$ -value (i.e., low SNR) diffusion data (up tob = 3000 $$ b=3000 $$ s/mm2 $$ {}^2 $$ ) in four healthy subjects. RESULTS Noise amplification associated with phase error correction incurred a 23% bias in phantom mean diffusivity (MD) measurements. Phantom MD estimates using the SPECTRE pipeline were within 8% of the ground truth value. In healthy volunteers, the SPECTRE pipeline visibly corrected Nyquist ghost artifacts and reduced associated noise amplification in highb $$ b $$ -value data. CONCLUSION The proposed reconstruction pipeline is effective in correcting low SNR data, and improves the accuracy and precision of derived diffusion metrics.
Collapse
Affiliation(s)
- Elizabeth Powell
- Queen Square MS Centre, UCL Institute of NeurologyUniversity College LondonLondonUK
- Centre for Medical Image Computing, Department of Medical Physics and Biomedical EngineeringUniversity College LondonLondonUK
| | | | - Marco Battiston
- Queen Square MS Centre, UCL Institute of NeurologyUniversity College LondonLondonUK
| | - Francesco Grussu
- Queen Square MS Centre, UCL Institute of NeurologyUniversity College LondonLondonUK
- Radiomics GroupVall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital CampusBarcelonaSpain
| | - Ahmed Toosy
- Queen Square MS Centre, UCL Institute of NeurologyUniversity College LondonLondonUK
| | - Jonathan D. Clayden
- Developmental Imaging and Biophysics Section, Great Ormond Street Institute of Child HealthUniversity College LondonLondonUK
| | - Claudia A. M. Gandini Wheeler‐Kingshott
- Queen Square MS Centre, UCL Institute of NeurologyUniversity College LondonLondonUK
- Department of Brain and Behavioural SciencesUniversity of PaviaPaviaItaly
- Brain MRI 3T CenterIRCCS Mondino FoundationPaviaItaly
| |
Collapse
|
16
|
Wood TC, Cash D, MacNicol E, Simmons C, Kim E, Lythgoe DJ, Zelaya F, Turkheimer F. Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents. Wellcome Open Res 2022; 6:109. [PMID: 36081865 PMCID: PMC9428501 DOI: 10.12688/wellcomeopenres.16734.4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/16/2022] [Indexed: 11/20/2022] Open
Abstract
Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.
Collapse
Affiliation(s)
- Tobias C Wood
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - Diana Cash
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - Eilidh MacNicol
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - Camilla Simmons
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - Eugene Kim
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - David J Lythgoe
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - Fernando Zelaya
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| | - Federico Turkheimer
- Department of Neuroimaging, Institute of Psychiatry, Psychology & Neuroscience, King's College London, SE5 8AF, UK
| |
Collapse
|
17
|
Wood TC, Cash D, MacNicol E, Simmons C, Kim E, Lythgoe DJ, Zelaya F, Turkheimer F. Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents. Wellcome Open Res 2022; 6:109. [DOI: 10.12688/wellcomeopenres.16734.3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/08/2022] [Indexed: 11/20/2022] Open
Abstract
Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.
Collapse
|
18
|
Fan Q, Eichner C, Afzali M, Mueller L, Tax CMW, Davids M, Mahmutovic M, Keil B, Bilgic B, Setsompop K, Lee HH, Tian Q, Maffei C, Ramos-Llordén G, Nummenmaa A, Witzel T, Yendiki A, Song YQ, Huang CC, Lin CP, Weiskopf N, Anwander A, Jones DK, Rosen BR, Wald LL, Huang SY. Mapping the human connectome using diffusion MRI at 300 mT/m gradient strength: Methodological advances and scientific impact. Neuroimage 2022; 254:118958. [PMID: 35217204 PMCID: PMC9121330 DOI: 10.1016/j.neuroimage.2022.118958] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 01/27/2022] [Accepted: 01/31/2022] [Indexed: 12/20/2022] Open
Abstract
Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength diffusion MRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for diffusion MRI and where the field is headed in the coming years.
Collapse
Affiliation(s)
- Qiuyun Fan
- Department of Biomedical Engineering, College of Precision Instruments and Optoelectronics Engineering, Tianjin University, Tianjin, China; Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Cornelius Eichner
- Max Planck Institute for Human Cognitive and Brain Sciences, Department of Neuropsychology, Leipzig, Germany
| | - Maryam Afzali
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, Wales, UK; Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds LS2 9JT, UK
| | - Lars Mueller
- Leeds Institute of Cardiovascular and Metabolic Medicine, University of Leeds, Leeds LS2 9JT, UK
| | - Chantal M W Tax
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, Wales, UK; Image Sciences Institute, University Medical Center (UMC) Utrecht, Utrecht, the Netherlands
| | - Mathias Davids
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Mirsad Mahmutovic
- Institute of Medical Physics and Radiation Protection (IMPS), TH-Mittelhessen University of Applied Sciences (THM), Giessen, Germany
| | - Boris Keil
- Institute of Medical Physics and Radiation Protection (IMPS), TH-Mittelhessen University of Applied Sciences (THM), Giessen, Germany
| | - Berkin Bilgic
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Kawin Setsompop
- Department of Radiology, Stanford University, Stanford, CA, USA; Department of Electrical Engineering, Stanford University, Stanford, CA, USA
| | - Hong-Hsi Lee
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Qiyuan Tian
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Chiara Maffei
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Gabriel Ramos-Llordén
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Aapo Nummenmaa
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | | | - Anastasia Yendiki
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA
| | - Yi-Qiao Song
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA USA
| | - Chu-Chung Huang
- Key Laboratory of Brain Functional Genomics (MOE & STCSM), Affiliated Mental Health Center (ECNU), School of Psychology and Cognitive Science, East China Normal University, Shanghai, China; Shanghai Changning Mental Health Center, Shanghai, China
| | - Ching-Po Lin
- Institute of Neuroscience, National Yang Ming Chiao Tung University, Taipei, Taiwan; Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai, China
| | - Nikolaus Weiskopf
- Department of Neurophysics, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany; Felix Bloch Institute for Solid State Physics, Faculty of Physics and Earth Sciences, Leipzig University, Leipzig, Germany
| | - Alfred Anwander
- Max Planck Institute for Human Cognitive and Brain Sciences, Department of Neuropsychology, Leipzig, Germany
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, Wales, UK
| | - Bruce R Rosen
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Lawrence L Wald
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Susie Y Huang
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, USA; Harvard Medical School, Boston, MA, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
19
|
Olesen JL, Østergaard L, Shemesh N, Jespersen SN. Diffusion time dependence, power-law scaling, and exchange in gray matter. Neuroimage 2022; 251:118976. [PMID: 35168088 PMCID: PMC8961002 DOI: 10.1016/j.neuroimage.2022.118976] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2021] [Revised: 12/24/2021] [Accepted: 02/04/2022] [Indexed: 12/27/2022] Open
Abstract
Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins' interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b-1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.
Collapse
Affiliation(s)
- Jonas L Olesen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Leif Østergaard
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Sune N Jespersen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark.
| |
Collapse
|
20
|
Grussu F, Bernatowicz K, Casanova-Salas I, Castro N, Nuciforo P, Mateo J, Barba I, Perez-Lopez R. Diffusion MRI signal cumulants and hepatocyte microstructure at fixed diffusion time: Insights from simulations, 9.4T imaging, and histology. Magn Reson Med 2022; 88:365-379. [PMID: 35181943 PMCID: PMC9303340 DOI: 10.1002/mrm.29174] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2021] [Revised: 12/21/2021] [Accepted: 01/07/2022] [Indexed: 11/09/2022]
Abstract
Purpose Relationships between diffusion‐weighted MRI signals and hepatocyte microstructure were investigated to inform liver diffusion MRI modeling, focusing on the following question: Can cell size and diffusivity be estimated at fixed diffusion time, realistic SNR, and negligible contribution from extracellular/extravascular water and exchange? Methods Monte Carlo simulations were performed within synthetic hepatocytes for varying cell size/diffusivity L/D0, and clinical protocols (single diffusion encoding; maximum b‐value: {1000, 1500, 2000} s/mm2; 5 unique gradient duration/separation pairs; SNR = {∞, 100, 80, 40, 20}), accounting for heterogeneity in (D0,L) and perfusion contamination. Diffusion (D) and kurtosis (K) coefficients were calculated, and relationships between (D0,L) and (D,K) were visualized. Functions mapping (D,K) to (D0,L) were computed to predict unseen (D0,L) values, tested for their ability to classify discrete cell‐size contrasts, and deployed on 9.4T ex vivo MRI‐histology data of fixed mouse livers Results Relationships between (D,K) and (D0,L) are complex and depend on the diffusion encoding. Functions mapping D,K to (D0,L) captures salient characteristics of D0(D,K) and L(D,K) dependencies. Mappings are not always accurate, but they enable just under 70% accuracy in a three‐class cell‐size classification task (for SNR = 20, bmax = 1500 s/mm2, δ = 20 ms, and Δ = 75 ms). MRI detects cell‐size contrasts in the mouse livers that are confirmed by histology, but overestimates the largest cell sizes. Conclusion Salient information about liver cell size and diffusivity may be retrieved from minimal diffusion encodings at fixed diffusion time, in experimental conditions and pathological scenarios for which extracellular, extravascular water and exchange are negligible.
Collapse
Affiliation(s)
- Francesco Grussu
- Radiomics Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Kinga Bernatowicz
- Radiomics Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Irene Casanova-Salas
- Prostate Cancer Translational Research Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Natalia Castro
- Prostate Cancer Translational Research Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Paolo Nuciforo
- Molecular Oncology Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Joaquin Mateo
- Prostate Cancer Translational Research Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Ignasi Barba
- NMR Lab, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain
| | - Raquel Perez-Lopez
- Radiomics Group, Vall d'Hebron Institute of Oncology, Vall d'Hebron Barcelona Hospital Campus, Barcelona, Spain.,Department of Radiology, Hospital Universitari Vall d'Hebron, Barcelona, Spain
| |
Collapse
|
21
|
Zhou X, Wang X, Liu E, Zhang L, Zhang H, Zhang X, Zhu Y, Kuai Z. An Unsupervised Deep Learning Approach for
Dynamic‐Exponential
Intravoxel Incoherent Motion
MRI
Modeling and Parameter Estimation in the Liver. J Magn Reson Imaging 2022; 56:848-859. [PMID: 35064945 DOI: 10.1002/jmri.28074] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Revised: 01/09/2022] [Accepted: 01/10/2022] [Indexed: 12/18/2022] Open
Affiliation(s)
- Xin‐Xiang Zhou
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| | - Xin‐Yu Wang
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| | - En‐Hui Liu
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| | - Lan Zhang
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| | - Hong‐Xia Zhang
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| | - Xiu‐Shi Zhang
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| | - Yue‐Min Zhu
- CREATIS CNRS UMR 5220‐INSERM U1206‐University Lyon 1‐INSA Lyon‐University Jean Monnet Saint‐Etienne Lyon France
| | - Zi‐Xiang Kuai
- Imaging Center Harbin Medical University Cancer Hospital Harbin China
| |
Collapse
|
22
|
Shukla V, Khandekar P, Khaparde A. Noise estimation in 2D MRI using DWT coefficients and optimized neural network. Biomed Signal Process Control 2022. [DOI: 10.1016/j.bspc.2021.103225] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
23
|
Nonparametric D-R 1-R 2 distribution MRI of the living human brain. Neuroimage 2021; 245:118753. [PMID: 34852278 DOI: 10.1016/j.neuroimage.2021.118753] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Revised: 11/17/2021] [Accepted: 11/22/2021] [Indexed: 11/23/2022] Open
Abstract
Diffusion-relaxation correlation NMR can simultaneously characterize both the microstructure and the local chemical composition of complex samples that contain multiple populations of water. Recent developments on tensor-valued diffusion encoding and Monte Carlo inversion algorithms have made it possible to transfer diffusion-relaxation correlation NMR from small-bore scanners to clinical MRI systems. Initial studies on clinical MRI systems employed 5D D-R1 and D-R2 correlation to characterize healthy brain in vivo. However, these methods are subject to an inherent bias that originates from not including R2 or R1 in the analysis, respectively. This drawback can be remedied by extending the concept to 6D D-R1-R2 correlation. In this work, we present a sparse acquisition protocol that records all data necessary for in vivo 6D D-R1-R2 correlation MRI across 633 individual measurements within 25 min-a time frame comparable to previous lower-dimensional acquisition protocols. The data were processed with a Monte Carlo inversion algorithm to obtain nonparametric 6D D-R1-R2 distributions. We validated the reproducibility of the method in repeated measurements of healthy volunteers. For a post-therapy glioblastoma case featuring cysts, edema, and partially necrotic remains of tumor, we present representative single-voxel 6D distributions, parameter maps, and artificial contrasts over a wide range of diffusion-, R1-, and R2-weightings based on the rich information contained in the D-R1-R2 distributions.
Collapse
|
24
|
Kurokawa R, Kamiya K, Koike S, Nakaya M, Uematsu A, Tanaka SC, Kamagata K, Okada N, Morita K, Kasai K, Abe O. Cross-scanner reproducibility and harmonization of a diffusion MRI structural brain network: A traveling subject study of multi-b acquisition. Neuroimage 2021; 245:118675. [PMID: 34710585 DOI: 10.1016/j.neuroimage.2021.118675] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Revised: 09/26/2021] [Accepted: 10/21/2021] [Indexed: 01/18/2023] Open
Abstract
Characterization of brain networks by diffusion MRI (dMRI) has rapidly evolved, and there are ongoing movements toward data sharing and multi-center studies. To extract meaningful information from multi-center data, methods to correct for the bias caused by scanner differences, that is, harmonization, are urgently needed. In this work, we report the cross-scanner differences in structural network analyses using data from nine traveling subjects (four males and five females, 21-49 years-old) who underwent scanning using four 3T scanners (public database available from the Brain/MINDS Beyond Human Brain MRI project (http://mriportal.umin.jp/)). The reliability and reproducibility were compared to those of data from another set of four subjects (all males, 29-42 years-old) who underwent scan-rescan (interval, 105-147 days) with the same scanner as well as scan-rescan data from the Human Connectome Project database. The results demonstrated that the reliability of the edge weights and graph theory metrics was lower for data including different scanners, compared to the scan-rescan with the same scanner. Besides, systematic differences between scanners were observed, indicating the risk of bias in comparing networks obtained from different scanners directly. We further demonstrate that it is feasible to reduce inter-scanner variabilities while preserving the inter-subject differences among healthy individuals by modeling the scanner effects at the level of network matrices, when traveling-subject data are available for calibration between scanners. The present data and results are expected to serve as a basis for developing and evaluating novel harmonization methods.
Collapse
Affiliation(s)
- Ryo Kurokawa
- Department of Radiology, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan.
| | - Kouhei Kamiya
- Department of Radiology, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan; Department of Radiology, Toho University, Tokyo, Japan; Department of Radiology, Juntendo University, Tokyo, Japan.
| | - Shinsuke Koike
- Center for Evolutionary Cognitive Sciences (ECS), Graduate School of Art and Sciences, The University of Tokyo, Tokyo, Japan; University of Tokyo Institute for Diversity & Adaptation of Human Mind (UTIDAHM), Tokyo, Japan; University of Tokyo Center for Integrative Science of Human Behavior (CiSHuB), Tokyo, Japan; The International Research Center for Neurointelligence (WPI-IRCN), Institutes for Advanced Study (UTIAS), The University of Tokyo, Tokyo, Japan.
| | - Moto Nakaya
- Department of Radiology, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan.
| | - Akiko Uematsu
- Center for Evolutionary Cognitive Sciences (ECS), Graduate School of Art and Sciences, The University of Tokyo, Tokyo, Japan.
| | - Saori C Tanaka
- Brain Information Communication Research Laboratory Group, Advanced Telecommunications Research Institutes International (ATR), Kyoto, Japan.
| | - Koji Kamagata
- Department of Radiology, Juntendo University, Tokyo, Japan.
| | - Naohiro Okada
- University of Tokyo Institute for Diversity & Adaptation of Human Mind (UTIDAHM), Tokyo, Japan; The International Research Center for Neurointelligence (WPI-IRCN), Institutes for Advanced Study (UTIAS), The University of Tokyo, Tokyo, Japan; Department of Neuropsychiatry, The University of Tokyo, Tokyo, Japan.
| | - Kentaro Morita
- Department of Neuropsychiatry, The University of Tokyo, Tokyo, Japan.
| | - Kiyoto Kasai
- University of Tokyo Institute for Diversity & Adaptation of Human Mind (UTIDAHM), Tokyo, Japan; University of Tokyo Center for Integrative Science of Human Behavior (CiSHuB), Tokyo, Japan; The International Research Center for Neurointelligence (WPI-IRCN), Institutes for Advanced Study (UTIAS), The University of Tokyo, Tokyo, Japan; Department of Neuropsychiatry, The University of Tokyo, Tokyo, Japan.
| | - Osamu Abe
- Department of Radiology, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan.
| |
Collapse
|
25
|
Kaandorp MPT, Barbieri S, Klaassen R, van Laarhoven HWM, Crezee H, While PT, Nederveen AJ, Gurney‐Champion OJ. Improved unsupervised physics-informed deep learning for intravoxel incoherent motion modeling and evaluation in pancreatic cancer patients. Magn Reson Med 2021; 86:2250-2265. [PMID: 34105184 PMCID: PMC8362093 DOI: 10.1002/mrm.28852] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Revised: 04/30/2021] [Accepted: 05/03/2021] [Indexed: 12/12/2022]
Abstract
PURPOSE Earlier work showed that IVIM-NETorig , an unsupervised physics-informed deep neural network, was faster and more accurate than other state-of-the-art intravoxel-incoherent motion (IVIM) fitting approaches to diffusion-weighted imaging (DWI). This study presents a substantially improved version, IVIM-NEToptim , and characterizes its superior performance in pancreatic cancer patients. METHOD In simulations (signal-to-noise ratio [SNR] = 20), the accuracy, independence, and consistency of IVIM-NET were evaluated for combinations of hyperparameters (fit S0, constraints, network architecture, number of hidden layers, dropout, batch normalization, learning rate), by calculating the normalized root-mean-square error (NRMSE), Spearman's ρ, and the coefficient of variation (CVNET ), respectively. The best performing network, IVIM-NEToptim was compared to least squares (LS) and a Bayesian approach at different SNRs. IVIM-NEToptim 's performance was evaluated in an independent dataset of 23 patients with pancreatic ductal adenocarcinoma. Fourteen of the patients received no treatment between two repeated scan sessions and nine received chemoradiotherapy between the repeated sessions. Intersession within-subject standard deviations (wSD) and treatment-induced changes were assessed. RESULTS In simulations (SNR = 20), IVIM-NEToptim outperformed IVIM-NETorig in accuracy (NRMSE(D) = 0.177 vs 0.196; NMRSE(f) = 0.220 vs 0.267; NMRSE(D*) = 0.386 vs 0.393), independence (ρ(D*, f) = 0.22 vs 0.74), and consistency (CVNET (D) = 0.013 vs 0.104; CVNET (f) = 0.020 vs 0.054; CVNET (D*) = 0.036 vs 0.110). IVIM-NEToptim showed superior performance to the LS and Bayesian approaches at SNRs < 50. In vivo, IVIM-NEToptim showed significantly less noisy parameter maps with lower wSD for D and f than the alternatives. In the treated cohort, IVIM-NEToptim detected the most individual patients with significant parameter changes compared to day-to-day variations. CONCLUSION IVIM-NEToptim is recommended for accurate, informative, and consistent IVIM fitting to DWI data.
Collapse
Affiliation(s)
- Misha P. T. Kaandorp
- Department of Radiology and Nuclear MedicineCancer Center Amsterdam, Amsterdam UMC, University of AmsterdamAmsterdamthe Netherlands
- Department of Radiology and Nuclear MedicineSt. Olav’s University HospitalTrondheimNorway
- Department of Circulation and Medical ImagingNTNU – Norwegian University of Science and TechnologyTrondheimNorway
| | | | - Remy Klaassen
- Department of Medical OncologyCancer Center Amsterdam, Amsterdam UMC, University of AmsterdamAmsterdamthe Netherlands
| | - Hanneke W. M. van Laarhoven
- Department of Medical OncologyCancer Center Amsterdam, Amsterdam UMC, University of AmsterdamAmsterdamthe Netherlands
| | - Hans Crezee
- Department of Radiology and Nuclear MedicineCancer Center Amsterdam, Amsterdam UMC, University of AmsterdamAmsterdamthe Netherlands
| | - Peter T. While
- Department of Radiology and Nuclear MedicineSt. Olav’s University HospitalTrondheimNorway
- Department of Circulation and Medical ImagingNTNU – Norwegian University of Science and TechnologyTrondheimNorway
| | - Aart J. Nederveen
- Department of Radiology and Nuclear MedicineCancer Center Amsterdam, Amsterdam UMC, University of AmsterdamAmsterdamthe Netherlands
| | - Oliver J. Gurney‐Champion
- Department of Radiology and Nuclear MedicineCancer Center Amsterdam, Amsterdam UMC, University of AmsterdamAmsterdamthe Netherlands
| |
Collapse
|
26
|
Gyori NG, Palombo M, Clark CA, Zhang H, Alexander DC. Training data distribution significantly impacts the estimation of tissue microstructure with machine learning. Magn Reson Med 2021; 87:932-947. [PMID: 34545955 DOI: 10.1002/mrm.29014] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Revised: 08/30/2021] [Accepted: 08/30/2021] [Indexed: 12/18/2022]
Abstract
PURPOSE Supervised machine learning (ML) provides a compelling alternative to traditional model fitting for parameter mapping in quantitative MRI. The aim of this work is to demonstrate and quantify the effect of different training data distributions on the accuracy and precision of parameter estimates when supervised ML is used for fitting. METHODS We fit a two- and three-compartment biophysical model to diffusion measurements from in-vivo human brain, as well as simulated diffusion data, using both traditional model fitting and supervised ML. For supervised ML, we train several artificial neural networks, as well as random forest regressors, on different distributions of ground truth parameters. We compare the accuracy and precision of parameter estimates obtained from the different estimation approaches using synthetic test data. RESULTS When the distribution of parameter combinations in the training set matches those observed in healthy human data sets, we observe high precision, but inaccurate estimates for atypical parameter combinations. In contrast, when training data is sampled uniformly from the entire plausible parameter space, estimates tend to be more accurate for atypical parameter combinations but may have lower precision for typical parameter combinations. CONCLUSION This work highlights that estimation of model parameters using supervised ML depends strongly on the training-set distribution. We show that high precision obtained using ML may mask strong bias, and visual assessment of the parameter maps is not sufficient for evaluating the quality of the estimates.
Collapse
Affiliation(s)
- Noemi G Gyori
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK.,Great Ormond Street Institute of Child Health, University College London, London, UK
| | - Marco Palombo
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK
| | - Christopher A Clark
- Great Ormond Street Institute of Child Health, University College London, London, UK
| | - Hui Zhang
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK
| | - Daniel C Alexander
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK
| |
Collapse
|
27
|
Kuczera S, Alipoor M, Langkilde F, Maier SE. Optimized bias and signal inference in diffusion-weighted image analysis (OBSIDIAN). Magn Reson Med 2021; 86:2716-2732. [PMID: 34278590 DOI: 10.1002/mrm.28773] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Revised: 01/29/2021] [Accepted: 02/24/2021] [Indexed: 12/25/2022]
Abstract
PURPOSE Correction of Rician signal bias in magnitude MR images. METHODS A model-based, iterative fitting procedure is used to simultaneously estimate true signal and underlying Gaussian noise with standard deviation σ g on a pixel-by-pixel basis in magnitude MR images. A precomputed function that relates absolute residuals between measured signals and model fit to σ g is used to iteratively estimate σ g . The feasibility of the method is evaluated and compared to maximum likelihood estimation (MLE) for diffusion signal decay simulations and diffusion-weighted images of the prostate considering 21 linearly spaced b-values from 0 to 3000 s/mm2 . A multidirectional analysis was performed with publically available brain data. RESULTS Model simulations show that the Rician bias correction algorithm is fast, with an accuracy and precision that is on par to model-based MLE and direct fitting in the case of pure Gaussian noise. Increased accuracy in parameter prediction in a low signal-to-noise ratio (SNR) scenario is ideally achieved by using a composite of multiple signal decays from neighboring voxels as input for the algorithm. For patient data, good agreement with high SNR reference data of diffusion in prostate is achieved. CONCLUSIONS OBSIDIAN is a novel, alternative, simple to implement approach for rapid Rician bias correction applicable in any case where differences between true signal decay and underlying model function can be considered negligible in comparison to noise. The proposed composite fitting approach permits accurate parameter estimation even in typical clinical scenarios with low SNR, which significantly simplifies comparison of complex diffusion parameters among studies.
Collapse
Affiliation(s)
- Stefan Kuczera
- Institute of Clinical Sciences, Sahlgrenska Academy, Gothenburg University, Gothenburg, Sweden.,MedTech West, Sahlgrenska University Hospital, Gothenburg, Sweden
| | - Mohammad Alipoor
- Institute of Clinical Sciences, Sahlgrenska Academy, Gothenburg University, Gothenburg, Sweden
| | - Fredrik Langkilde
- Institute of Clinical Sciences, Sahlgrenska Academy, Gothenburg University, Gothenburg, Sweden
| | - Stephan E Maier
- Institute of Clinical Sciences, Sahlgrenska Academy, Gothenburg University, Gothenburg, Sweden.,Department of Radiology, Brigham Women's Hospital, Harvard Medical School, Boston, MA, USA
| |
Collapse
|
28
|
Lee HH, Novikov DS, Fieremans E. Removal of partial Fourier-induced Gibbs (RPG) ringing artifacts in MRI. Magn Reson Med 2021; 86:2733-2750. [PMID: 34227142 DOI: 10.1002/mrm.28830] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 03/29/2021] [Accepted: 04/16/2021] [Indexed: 01/03/2023]
Abstract
PURPOSE To investigate and remove Gibbs-ringing artifacts caused by partial Fourier (PF) acquisition and zero filling interpolation in MRI data. THEORY AND METHODS Gibbs ringing of fully sampled data, leading to oscillations around tissue boundaries, is caused by the symmetric truncation of k-space. Such ringing can be removed by conventional methods, with the local subvoxel shifts method being the state-of-the-art. However, the asymmetric truncation of k-space in routinely used PF acquisitions leads to additional ringings of wider intervals in the PF sampling dimension that cannot be corrected solely based on magnitude images reconstructed via zero filling. Here, we develop a pipeline for the Removal of PF-induced Gibbs ringing (RPG) to remove ringing patterns of different periods by applying the conventional method twice. The proposed pipeline is validated on numerical phantoms, demonstrated on in vivo diffusion MRI measurements, and compared with the conventional method and neural network-based approach. RESULTS For PF = 7/8 and 6/8, Gibbs-ringings and subsequent bias in diffusion metrics induced by PF acquisition and zero filling are robustly removed by using the proposed RPG pipeline. For PF = 5/8, however, ringing removal via RPG leads to excessive image blurring due to the interplay of image phase and convolution kernel. CONCLUSIONS RPG corrects Gibbs-ringing artifacts in magnitude images of PF acquired data and reduces the bias in quantitative MR metrics. Considering the benefit of PF acquisition and the feasibility of ringing removal, we suggest applying PF = 6/8 when PF acquisition is necessary.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, New York University School of Medicine, New York, NY, USA
| | - Dmitry S Novikov
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, New York University School of Medicine, New York, NY, USA
| | - Els Fieremans
- Center for Advanced Imaging Innovation and Research (CAI2R), Department of Radiology, New York University School of Medicine, New York, NY, USA
| |
Collapse
|
29
|
Cai LY, Yang Q, Hansen CB, Nath V, Ramadass K, Johnson GW, Conrad BN, Boyd BD, Begnoche JP, Beason-Held LL, Shafer AT, Resnick SM, Taylor WD, Price GR, Morgan VL, Rogers BP, Schilling KG, Landman BA. PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images. Magn Reson Med 2021; 86:456-470. [PMID: 33533094 PMCID: PMC8387107 DOI: 10.1002/mrm.28678] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Revised: 12/19/2020] [Accepted: 12/22/2020] [Indexed: 12/31/2022]
Abstract
PURPOSE Diffusion weighted MRI imaging (DWI) is often subject to low signal-to-noise ratios (SNRs) and artifacts. Recent work has produced software tools that can correct individual problems, but these tools have not been combined with each other and with quality assurance (QA). A single integrated pipeline is proposed to perform DWI preprocessing with a spectrum of tools and produce an intuitive QA document. METHODS The proposed pipeline, built around the FSL, MRTrix3, and ANTs software packages, performs DWI denoising; inter-scan intensity normalization; susceptibility-, eddy current-, and motion-induced artifact correction; and slice-wise signal drop-out imputation. To perform QA on the raw and preprocessed data and each preprocessing operation, the pipeline documents qualitative visualizations, quantitative plots, gradient verifications, and tensor goodness-of-fit and fractional anisotropy analyses. RESULTS Raw DWI data were preprocessed and quality checked with the proposed pipeline and demonstrated improved SNRs; physiologic intensity ratios; corrected susceptibility-, eddy current-, and motion-induced artifacts; imputed signal-lost slices; and improved tensor fits. The pipeline identified incorrect gradient configurations and file-type conversion errors and was shown to be effective on externally available datasets. CONCLUSIONS The proposed pipeline is a single integrated pipeline that combines established diffusion preprocessing tools from major MRI-focused software packages with intuitive QA.
Collapse
Affiliation(s)
- Leon Y. Cai
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Qi Yang
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN, USA
| | - Colin B. Hansen
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN, USA
| | - Vishwesh Nath
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN, USA
| | - Karthik Ramadass
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN, USA
| | - Graham W. Johnson
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Benjamin N. Conrad
- Neuroscience Graduate Program, Vanderbilt Brain Institute, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Psychology and Human Development, Peabody College, Vanderbilt University, Nashville, TN, USA
| | - Brian D. Boyd
- Department of Psychiatry and Behavioral Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Center for Cognitive Medicine, Vanderbilt University Medical Center, Nashville, TN, USA
| | - John P. Begnoche
- Department of Psychiatry and Behavioral Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Center for Cognitive Medicine, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Lori L. Beason-Held
- Laboratory of Behavioral Neuroscience, National Institute on Aging, National Institutes of Health, Baltimore, MD, USA
| | - Andrea T. Shafer
- Laboratory of Behavioral Neuroscience, National Institute on Aging, National Institutes of Health, Baltimore, MD, USA
| | - Susan M. Resnick
- Laboratory of Behavioral Neuroscience, National Institute on Aging, National Institutes of Health, Baltimore, MD, USA
| | - Warren D. Taylor
- Department of Psychiatry and Behavioral Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Center for Cognitive Medicine, Vanderbilt University Medical Center, Nashville, TN, USA
| | - Gavin R. Price
- Department of Psychology and Human Development, Peabody College, Vanderbilt University, Nashville, TN, USA
| | - Victoria L. Morgan
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Neurological Surgery, Vanderbilt University Medical Center, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, TN, USA
| | - Baxter P. Rogers
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Department of Psychiatry and Behavioral Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, TN, USA
| | - Kurt G. Schilling
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, TN, USA
| | - Bennett A. Landman
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
- Department of Electrical Engineering and Computer Science, Vanderbilt University, Nashville, TN, USA
- Department of Psychiatry and Behavioral Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Department of Radiology and Radiological Sciences, Vanderbilt University Medical Center, Nashville, TN, USA
- Vanderbilt University Institute of Imaging Science, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
30
|
Vukovic N, Hansen B, Lund TE, Jespersen S, Shtyrov Y. Rapid microstructural plasticity in the cortical semantic network following a short language learning session. PLoS Biol 2021; 19:e3001290. [PMID: 34125828 PMCID: PMC8202930 DOI: 10.1371/journal.pbio.3001290] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 05/17/2021] [Indexed: 01/22/2023] Open
Abstract
Despite the clear importance of language in our life, our vital ability to quickly and effectively learn new words and meanings is neurobiologically poorly understood. Conventional knowledge maintains that language learning—especially in adulthood—is slow and laborious. Furthermore, its structural basis remains unclear. Even though behavioural manifestations of learning are evident near instantly, previous neuroimaging work across a range of semantic categories has largely studied neural changes associated with months or years of practice. Here, we address rapid neuroanatomical plasticity accompanying new lexicon acquisition, specifically focussing on the learning of action-related language, which has been linked to the brain’s motor systems. Our results show that it is possible to measure and to externally modulate (using transcranial magnetic stimulation (TMS) of motor cortex) cortical microanatomic reorganisation after mere minutes of new word learning. Learning-induced microstructural changes, as measured by diffusion kurtosis imaging (DKI) and machine learning-based analysis, were evident in prefrontal, temporal, and parietal neocortical sites, likely reflecting integrative lexico-semantic processing and formation of new memory circuits immediately during the learning tasks. These results suggest a structural basis for the rapid neocortical word encoding mechanism and reveal the causally interactive relationship of modal and associative brain regions in supporting learning and word acquisition. This combined neuroimaging and brain stimulation study reveals rapid and distributed microstructural plasticity after a single immersive language learning session, demonstrating the causal relevance of the motor cortex in encoding the meaning of novel action words.
Collapse
Affiliation(s)
- Nikola Vukovic
- Department of Psychiatry, University of California San Francisco, San Francisco, United States of America
- * E-mail:
| | - Brian Hansen
- Center of Functionally Integrative Neuroscience, Aarhus University, Aarhus, Denmark
| | | | - Sune Jespersen
- Center of Functionally Integrative Neuroscience, Aarhus University, Aarhus, Denmark
- Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Yury Shtyrov
- Center of Functionally Integrative Neuroscience, Aarhus University, Aarhus, Denmark
- Centre for Cognition and Decision making, HSE University, Moscow, Russia
| |
Collapse
|
31
|
Olesen JL, Østergaard L, Shemesh N, Jespersen SN. Beyond the diffusion standard model in fixed rat spinal cord with combined linear and planar encoding. Neuroimage 2021; 231:117849. [PMID: 33582270 DOI: 10.1016/j.neuroimage.2021.117849] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 01/20/2021] [Accepted: 02/04/2021] [Indexed: 10/22/2022] Open
Abstract
Information about tissue on the microscopic and mesoscopic scales can be accessed by modelling diffusion MRI signals, with the aim of extracting microstructure-specific biomarkers. The standard model (SM) of diffusion, currently the most broadly adopted microstructural model, describes diffusion in white matter (WM) tissues by two Gaussian components, one of which has zero radial diffusivity, to represent diffusion in intra- and extra-axonal water, respectively. Here, we reappraise these SM assumptions by collecting comprehensive double diffusion encoded (DDE) MRI data with both linear and planar encodings, which was recently shown to substantially enhance the ability to estimate SM parameters. We find however, that the SM is unable to account for data recorded in fixed rat spinal cord at an ultrahigh field of 16.4 T, suggesting that its underlying assumptions are violated in our experimental data. We offer three model extensions to mitigate this problem: first, we generalize the SM to accommodate finite radii (axons) by releasing the constraint of zero radial diffusivity in the intra-axonal compartment. Second, we include intracompartmental kurtosis to account for non-Gaussian behaviour. Third, we introduce an additional (third) compartment. The ability of these models to account for our experimental data are compared based on parameter feasibility and Bayesian information criterion. Our analysis identifies the three-compartment description as the optimal model. The third compartment exhibits slow diffusion with a minor but non-negligible signal fraction (∼12%). We demonstrate how failure to take the presence of such a compartment into account severely misguides inferences about WM microstructure. Our findings bear significance for microstructural modelling at large and can impact the interpretation of biomarkers extracted from the standard model of diffusion.
Collapse
Affiliation(s)
- Jonas L Olesen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
| | - Leif Østergaard
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
| | - Noam Shemesh
- Champalimaud Research, Champalimaud Centre for the Unknown, Lisbon, Portugal
| | - Sune N Jespersen
- Center of Functionally Integrative Neuroscience (CFIN) and MINDLab, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark.
| |
Collapse
|
32
|
Thaler C, Kyselyova AA, Faizy TD, Nawka MT, Jespersen S, Hansen B, Stellmann JP, Heesen C, Stürner KH, Stark M, Fiehler J, Bester M, Gellißen S. Heterogeneity of multiple sclerosis lesions in fast diffusional kurtosis imaging. PLoS One 2021; 16:e0245844. [PMID: 33539364 PMCID: PMC7861404 DOI: 10.1371/journal.pone.0245844] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Accepted: 01/09/2021] [Indexed: 12/14/2022] Open
Abstract
Background Mean kurtosis (MK), one of the parameters derived from diffusion kurtosis imaging (DKI), has shown increased sensitivity to tissue microstructure damage in several neurological disorders. Methods Thirty-seven patients with relapsing-remitting MS and eleven healthy controls (HC) received brain imaging on a 3T MR scanner, including a fast DKI sequence. MK and mean diffusivity (MD) were measured in the white matter of HC, normal-appearing white matter (NAWM) of MS patients, contrast-enhancing lesions (CE-L), FLAIR lesions (FLAIR-L) and black holes (BH). Results Overall 1529 lesions were analyzed, including 30 CE-L, 832 FLAIR-L and 667 BH. Highest MK values were obtained in the white matter of HC (0.814 ± 0.129), followed by NAWM (0.724 ± 0.137), CE-L (0.619 ± 0.096), FLAIR-L (0.565 ± 0.123) and BH (0.549 ± 0.12). Lowest MD values were obtained in the white matter of HC (0.747 ± 0.068 10−3mm2/sec), followed by NAWM (0.808 ± 0.163 10−3mm2/sec), CE-L (0.853 ± 0.211 10−3mm2/sec), BH (0.957 ± 0.304 10−3mm2/sec) and FLAIR-L (0.976 ± 0.35 10−3mm2/sec). While MK differed significantly between CE-L and non-enhancing lesions, MD did not. Conclusion MK adds predictive value to differentiate between MS lesions and might provide further information about diffuse white matter injury and lesion microstructure.
Collapse
Affiliation(s)
- Christian Thaler
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- * E-mail:
| | - Anna A. Kyselyova
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Tobias D. Faizy
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Marie T. Nawka
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Sune Jespersen
- Department of Clinical Medicine - Center of Functionally Integrative Neuroscience, Aarhus University, Aarhus, Denmark
| | - Brian Hansen
- Department of Clinical Medicine - Center of Functionally Integrative Neuroscience, Aarhus University, Aarhus, Denmark
| | - Jan-Patrick Stellmann
- Department of Neurology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Institute for Neuroimmunology and Clinical MS Research, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- APHM, Hospital de la Timone, CEMEREM, Marseille, France
- Aix Marseille University, CNRS, CRMBM, UMR 7339, Marseille, France
| | - Christoph Heesen
- Department of Neurology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Institute for Neuroimmunology and Clinical MS Research, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Klarissa H. Stürner
- Department of Neurology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Institute for Neuroimmunology and Clinical MS Research, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Department of Neurology, University Hospital Schleswig-Holstein, Kiel, Germany
| | - Maria Stark
- Institute of Medical Biometry and Epidemiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Jens Fiehler
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Maxim Bester
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Susanne Gellißen
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| |
Collapse
|
33
|
Data Preparation Protocol for Low Signal-to-Noise Ratio Fluorine-19 MRI. Methods Mol Biol 2021. [PMID: 33476033 DOI: 10.1007/978-1-0716-0978-1_43] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
Fluorine-19 MRI shows great promise for a wide range of applications including renal imaging, yet the typically low signal-to-noise ratios and sparse signal distribution necessitate a thorough data preparation.This chapter describes a general data preparation workflow for fluorine MRI experiments. The main processing steps are: (1) estimation of noise level, (2) correction of noise-induced bias and (3) background subtraction. The protocol is supplemented by an example script and toolbox available online.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers. This analysis protocol chapter is complemented by two separate chapters describing the basic concept and experimental procedure.
Collapse
|
34
|
Afzali M, Pieciak T, Newman S, Garyfallidis E, Özarslan E, Cheng H, Jones DK. The sensitivity of diffusion MRI to microstructural properties and experimental factors. J Neurosci Methods 2021; 347:108951. [PMID: 33017644 PMCID: PMC7762827 DOI: 10.1016/j.jneumeth.2020.108951] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2020] [Revised: 08/27/2020] [Accepted: 09/15/2020] [Indexed: 12/13/2022]
Abstract
Diffusion MRI is a non-invasive technique to study brain microstructure. Differences in the microstructural properties of tissue, including size and anisotropy, can be represented in the signal if the appropriate method of acquisition is used. However, to depict the underlying properties, special care must be taken when designing the acquisition protocol as any changes in the procedure might impact on quantitative measurements. This work reviews state-of-the-art methods for studying brain microstructure using diffusion MRI and their sensitivity to microstructural differences and various experimental factors. Microstructural properties of the tissue at a micrometer scale can be linked to the diffusion signal at a millimeter-scale using modeling. In this paper, we first give an introduction to diffusion MRI and different encoding schemes. Then, signal representation-based methods and multi-compartment models are explained briefly. The sensitivity of the diffusion MRI signal to the microstructural components and the effects of curvedness of axonal trajectories on the diffusion signal are reviewed. Factors that impact on the quality (accuracy and precision) of derived metrics are then reviewed, including the impact of random noise, and variations in the acquisition parameters (i.e., number of sampled signals, b-value and number of acquisition shells). Finally, yet importantly, typical approaches to deal with experimental factors are depicted, including unbiased measures and harmonization. We conclude the review with some future directions and recommendations on this topic.
Collapse
Affiliation(s)
- Maryam Afzali
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom.
| | - Tomasz Pieciak
- AGH University of Science and Technology, Kraków, Poland; LPI, ETSI Telecomunicación, Universidad de Valladolid, Valladolid, Spain.
| | - Sharlene Newman
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN 47405, USA; Program of Neuroscience, Indiana University, Bloomington, IN 47405, USA.
| | - Eleftherios Garyfallidis
- Program of Neuroscience, Indiana University, Bloomington, IN 47405, USA; Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN 47408, USA.
| | - Evren Özarslan
- Department of Biomedical Engineering, Linköping University, Linköping, Sweden; Center for Medical Image Science and Visualization, Linköping University, Linköping, Sweden.
| | - Hu Cheng
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN 47405, USA; Program of Neuroscience, Indiana University, Bloomington, IN 47405, USA.
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Psychology, Cardiff University, Cardiff, United Kingdom.
| |
Collapse
|
35
|
Lee HH, Papaioannou A, Novikov DS, Fieremans E. In vivo observation and biophysical interpretation of time-dependent diffusion in human cortical gray matter. Neuroimage 2020; 222:117054. [PMID: 32585341 PMCID: PMC7736473 DOI: 10.1016/j.neuroimage.2020.117054] [Citation(s) in RCA: 43] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2019] [Revised: 06/08/2020] [Accepted: 06/10/2020] [Indexed: 12/25/2022] Open
Abstract
The dependence of the diffusion MRI signal on the diffusion time t is a hallmark of tissue microstructure at the scale of the diffusion length. Here we measure the time-dependence of the mean diffusivity D(t) and mean kurtosis K(t) in cortical gray matter and in 25 gray matter sub-regions, in 10 healthy subjects. Significant diffusivity and kurtosis time-dependence is observed for t=21.2-100 ms, and is characterized by a power-law tail ∼t-ϑ with dynamical exponent ϑ. To interpret our measurements, we systematize the relevant scenarios and mechanisms for diffusion time-dependence in the brain. Using the effective medium theory formalism, we derive an exact relation between the power-law tails in D(t) and K(t). The estimated dynamical exponent ϑ≃1/2 in both D(t) and K(t) is consistent with one-dimensional diffusion in the presence of randomly positioned restrictions along neurites. We analyze the short-range disordered statistics of synapses on axon collaterals in the cortex, and perform one-dimensional Monte Carlo simulations of diffusion restricted by permeable barriers with a similar randomness in their placement, to confirm the ϑ=1/2 exponent. In contrast, the Kärger model of exchange is less consistent with the data since it does not capture the diffusivity time-dependence, and the estimated exchange time from K(t) falls below our measured t-range. Although we cannot exclude exchange as a contributing factor, we argue that structural disorder along neurites is mainly responsible for the observed time-dependence of diffusivity and kurtosis. Our observation and theoretical interpretation of the t-1/2 tail in D(t) and K(t) altogether establish the sensitivity of a macroscopic MRI signal to micrometer-scale structural heterogeneities along neurites in human gray matter in vivo.
Collapse
Affiliation(s)
- Hong-Hsi Lee
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA.
| | - Antonios Papaioannou
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, NY, USA; Center for Advanced Imaging Innovation and Research (CAI2R), New York University School of Medicine, New York, NY, USA
| |
Collapse
|
36
|
Pestryaev EM. Chain Heterogeneity in Simulated Polymer Melts: NMR Free Induction Decay and Absorption Line. POLYMER SCIENCE SERIES A 2020. [DOI: 10.1134/s0965545x20060097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
37
|
Pestryaev EM. Chain Heterogeneity in Simulated Polymer Melts: Segment Orientational Autocorrelation Function. POLYMER SCIENCE SERIES A 2020. [DOI: 10.1134/s0965545x20060085] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
38
|
Sudeep P, Palanisamy P, Kesavadas C, Rajan J. An improved nonlocal maximum likelihood estimation method for denoising magnetic resonance images with spatially varying noise levels. Pattern Recognit Lett 2020. [DOI: 10.1016/j.patrec.2018.02.007] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
39
|
Grussu F, Battiston M, Veraart J, Schneider T, Cohen-Adad J, Shepherd TM, Alexander DC, Fieremans E, Novikov DS, Gandini Wheeler-Kingshott CAM. Multi-parametric quantitative in vivo spinal cord MRI with unified signal readout and image denoising. Neuroimage 2020; 217:116884. [PMID: 32360689 PMCID: PMC7378937 DOI: 10.1016/j.neuroimage.2020.116884] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 03/18/2020] [Accepted: 04/23/2020] [Indexed: 12/11/2022] Open
Abstract
Multi-parametric quantitative MRI (qMRI) of the spinal cord is a promising non-invasive tool to probe early microstructural damage in neurological disorders. It is usually performed in vivo by combining acquisitions with multiple signal readouts, which exhibit different thermal noise levels, geometrical distortions and susceptibility to physiological noise. This ultimately hinders joint multi-contrast modelling and makes the geometric correspondence of parametric maps challenging. We propose an approach to overcome these limitations, by implementing state-of-the-art microstructural MRI of the spinal cord with a unified signal readout in vivo (i.e. with matched spatial encoding parameters across a range of imaging contrasts). We base our acquisition on single-shot echo planar imaging with reduced field-of-view, and obtain data from two different vendors (vendor 1: Philips Achieva; vendor 2: Siemens Prisma). Importantly, the unified acquisition allows us to compare signal and noise across contrasts, thus enabling overall quality enhancement via multi-contrast image denoising methods. As a proof-of-concept, here we provide a demonstration with one such method, known as Marchenko-Pastur (MP) Principal Component Analysis (PCA) denoising. MP-PCA is a singular value (SV) decomposition truncation approach that relies on redundant acquisitions, i.e. such that the number of measurements is large compared to the number of components that are maintained in the truncated SV decomposition. Here we used in vivo and synthetic data to test whether a unified readout enables more efficient MP-PCA denoising of less redundant acquisitions, since these can be denoised jointly with more redundant ones. We demonstrate that a unified readout provides robust multi-parametric maps, including diffusion and kurtosis tensors from diffusion MRI, myelin metrics from two-pool magnetisation transfer, and T1 and T2 from relaxometry. Moreover, we show that MP-PCA improves the quality of our multi-contrast acquisitions, since it reduces the coefficient of variation (i.e. variability) by up to 17% for mean kurtosis, 8% for bound pool fraction (myelin-sensitive), and 13% for T1, while enabling more efficient denoising of modalities limited in redundancy (e.g. relaxometry). In conclusion, multi-parametric spinal cord qMRI with unified readout is feasible and provides robust microstructural metrics with matched resolution and distortions, whose quality benefits from multi-contrast denoising methods such as MP-PCA.
Collapse
Affiliation(s)
- Francesco Grussu
- Queen Square MS Centre, UCL Queen Square Institute of Neurology, Faculty of Brain Sciences, University College London, London, UK; Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK.
| | - Marco Battiston
- Queen Square MS Centre, UCL Queen Square Institute of Neurology, Faculty of Brain Sciences, University College London, London, UK
| | - Jelle Veraart
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | | | - Julien Cohen-Adad
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montreal, Montreal, Canada; Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montreal, Canada
| | - Timothy M Shepherd
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | - Daniel C Alexander
- Centre for Medical Image Computing, Department of Computer Science, University College London, London, UK
| | - Els Fieremans
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of Radiology, New York University School of Medicine, New York, USA
| | - Claudia A M Gandini Wheeler-Kingshott
- Queen Square MS Centre, UCL Queen Square Institute of Neurology, Faculty of Brain Sciences, University College London, London, UK; Brain MRI 3T Research Centre, IRCCS Mondino Foundation, Pavia, Italy; Department of Brain and Behavioural Sciences, University of Pavia, Pavia, Italy
| |
Collapse
|
40
|
St-Jean S, De Luca A, Tax CMW, Viergever MA, Leemans A. Automated characterization of noise distributions in diffusion MRI data. Med Image Anal 2020; 65:101758. [PMID: 32599491 DOI: 10.1016/j.media.2020.101758] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Revised: 06/11/2020] [Accepted: 06/16/2020] [Indexed: 02/07/2023]
Abstract
Knowledge of the noise distribution in magnitude diffusion MRI images is the centerpiece to quantify uncertainties arising from the acquisition process. The use of parallel imaging methods, the number of receiver coils and imaging filters applied by the scanner, amongst other factors, dictate the resulting signal distribution. Accurate estimation beyond textbook Rician or noncentral chi distributions often requires information about the acquisition process (e.g., coils sensitivity maps or reconstruction coefficients), which is usually not available. We introduce two new automated methods using the moments and maximum likelihood equations of the Gamma distribution to estimate noise distributions as they explicitly depend on the number of coils, making it possible to estimate all unknown parameters using only the magnitude data. A rejection step is used to make the framework automatic and robust to artifacts. Simulations using stationary and spatially varying noncentral chi noise distributions were created for two diffusion weightings with SENSE or GRAPPA reconstruction and 8, 12 or 32 receiver coils. Furthermore, MRI data of a water phantom with different combinations of parallel imaging were acquired on a 3T Philips scanner along with noise-only measurements. Finally, experiments on freely available datasets from a single subject acquired on a 3T GE scanner are used to assess reproducibility when limited information about the acquisition protocol is available. Additionally, we demonstrated the applicability of the proposed methods for a bias correction and denoising task on an in vivo dataset acquired on a 3T Siemens scanner. A generalized version of the bias correction framework for non integer degrees of freedom is also introduced. The proposed framework is compared with three other algorithms with datasets from three vendors, employing different reconstruction methods. Simulations showed that assuming a Rician distribution can lead to misestimation of the noise distribution in parallel imaging. Results on the acquired datasets showed that signal leakage in multiband can also lead to a misestimation of the noise distribution. Repeated acquisitions of in vivo datasets show that the estimated parameters are stable and have lower variability than compared methods. Results for the bias correction and denoising task show that the proposed methods reduce the appearance of noise at high b-value. The proposed algorithms herein can estimate both parameters of the noise distribution automatically, are robust to signal leakage artifacts and perform best when used on acquired noise maps.
Collapse
Affiliation(s)
- Samuel St-Jean
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, the Netherlands.
| | - Alberto De Luca
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, the Netherlands.
| | - Chantal M W Tax
- Cardiff University Brain Research Imaging Centre (CUBRIC), School of Physics and Astronomy, Cardiff University, Cardiff, United Kingdom.
| | - Max A Viergever
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, the Netherlands.
| | - Alexander Leemans
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, the Netherlands.
| |
Collapse
|
41
|
Goodburn RJ, Barrett T, Patterson I, Gallagher FA, Lawrence EM, Gnanapragasam VJ, Kastner C, Priest AN. Removing rician bias in diffusional kurtosis of the prostate using real-data reconstruction. Magn Reson Med 2020; 83:2243-2252. [PMID: 31737935 PMCID: PMC7065237 DOI: 10.1002/mrm.28080] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2018] [Revised: 10/22/2019] [Accepted: 10/23/2019] [Indexed: 12/24/2022]
Abstract
PURPOSE To compare prostate diffusional kurtosis imaging (DKI) metrics generated using phase-corrected real data with those generated using magnitude data with and without noise compensation (NC). METHODS Diffusion-weighted images were acquired at 3T in 16 prostate cancer patients, measuring 6 b-values (0-1500 s/mm2 ), each acquired with 6 signal averages along 3 diffusion directions, with noise-only images acquired to allow NC. In addition to conventional magnitude averaging, phase-corrected real data were averaged in an attempt to reduce rician noise-bias, with a range of phase-correction low-pass filter (LPF) sizes (8-128 pixels) tested. Each method was also tested using simulations. Pixelwise maps of apparent diffusion (D) and apparent kurtosis (K) were calculated for magnitude data with and without NC and phase-corrected real data. Average values were compared in tumor, normal transition zone (NTZ), and normal peripheral zone (NPZ). RESULTS Simulations indicated LPF size can strongly affect K metrics, where 64-pixel LPFs produced accurate metrics. Relative to metrics estimated from magnitude data without NC, median NC K were lower (P < 0.0001) by 6/11/8% in tumor/NPZ/NTZ, 64-LPF real-data K were lower (P < 0.0001) by 4/10/7%, respectively. CONCLUSION Compared with magnitude data with NC, phase-corrected real data can produce similar K, although the choice of phase-correction LPF should be chosen carefully.
Collapse
Affiliation(s)
- Rosie J. Goodburn
- Department of Medical PhysicsCambridge University Hospitals NHS Foundation TrustCambridgeUnited Kingdom
- Division of Radiotherapy and ImagingThe Institute of Cancer ResearchLondon
| | - Tristan Barrett
- Department of RadiologySchool of Clinical MedicineUniversity of CambridgeCambridgeUnited Kingdom
| | - Ilse Patterson
- Department of RadiologyCambridge University Hospitals NHS Foundation TrustCambridgeUnited Kingdom
| | - Ferdia A. Gallagher
- Department of RadiologySchool of Clinical MedicineUniversity of CambridgeCambridgeUnited Kingdom
| | - Edward M. Lawrence
- Department of RadiologySchool of Clinical MedicineUniversity of CambridgeCambridgeUnited Kingdom
| | | | - Christof Kastner
- Department of UrologyCambridge University Hospitals NHS Foundation TrustCambridgeUnited Kingdom
| | - Andrew N. Priest
- Department of RadiologySchool of Clinical MedicineUniversity of CambridgeCambridgeUnited Kingdom
- Department of RadiologyCambridge University Hospitals NHS Foundation TrustCambridgeUnited Kingdom
| |
Collapse
|
42
|
Das P, Pal C, Chakrabarti A, Acharyya A, Basu S. Adaptive denoising of 3D volumetric MR images using local variance based estimator. Biomed Signal Process Control 2020. [DOI: 10.1016/j.bspc.2020.101901] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
|
43
|
Denoise magnitude diffusion magnetic resonance images via variance-stabilizing transformation and optimal singular-value manipulation. Neuroimage 2020; 215:116852. [PMID: 32305566 PMCID: PMC7292796 DOI: 10.1016/j.neuroimage.2020.116852] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2019] [Revised: 04/07/2020] [Accepted: 04/10/2020] [Indexed: 12/12/2022] Open
Abstract
Although shown to have a great utility for a wide range of neuroscientific and clinical applications, diffusion-weighted magnetic resonance imaging (dMRI) faces a major challenge of low signal-to-noise ratio (SNR), especially when pushing the spatial resolution for improved delineation of brain's fine structure or increasing the diffusion weighting for increased angular contrast or both. Here, we introduce a comprehensive denoising framework for denoising magnitude dMRI. The framework synergistically combines the variance stabilizing transform (VST) with optimal singular value manipulation. The purpose of VST is to transform the Rician data to Gaussian-like data so that an asymptotically optimal singular value manipulation strategy tailored for Gaussian data can be used. The output of the framework is the estimated underlying diffusion signal for each voxel in the image domain. The usefulness of the proposed framework for denoising magnitude dMRI is demonstrated using both simulation and real-data experiments. Our results show that the proposed denoising framework can significantly improve SNR across the entire brain, leading to substantially enhanced performances for estimating diffusion tensor related indices and for resolving crossing fibers when compared to another competing method. More encouragingly, the proposed method when used to denoise a single average of 7 Tesla Human Connectome Project-style diffusion acquisition provided comparable performances relative to those achievable with ten averages for resolving multiple fiber populations across the brain. As such, the proposed denoising method is expected to have a great utility for high-quality, high-resolution whole-brain dMRI, desirable for many neuroscientific and clinical applications.
Collapse
|
44
|
Veraart J, Nunes D, Rudrapatna U, Fieremans E, Jones DK, Novikov DS, Shemesh N. Nonivasive quantification of axon radii using diffusion MRI. eLife 2020; 9:e49855. [PMID: 32048987 PMCID: PMC7015669 DOI: 10.7554/elife.49855] [Citation(s) in RCA: 95] [Impact Index Per Article: 23.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Accepted: 01/07/2020] [Indexed: 12/13/2022] Open
Abstract
Axon caliber plays a crucial role in determining conduction velocity and, consequently, in the timing and synchronization of neural activation. Noninvasive measurement of axon radii could have significant impact on the understanding of healthy and diseased neural processes. Until now, accurate axon radius mapping has eluded in vivo neuroimaging, mainly due to a lack of sensitivity of the MRI signal to micron-sized axons. Here, we show how - when confounding factors such as extra-axonal water and axonal orientation dispersion are eliminated - heavily diffusion-weighted MRI signals become sensitive to axon radii. However, diffusion MRI is only capable of estimating a single metric, the effective radius, representing the entire axon radius distribution within a voxel that emphasizes the larger axons. Our findings, both in rodents and humans, enable noninvasive mapping of critical information on axon radii, as well as resolve the long-standing debate on whether axon radii can be quantified.
Collapse
Affiliation(s)
- Jelle Veraart
- Champalimaud ResearchChampalimaud Centre for the UnknownLisbonPortugal
- Center for Biomedical Imaging, Department of RadiologyNew York University School of MedicineNew YorkUnited States
- imec-Vision Lab, Department of PhysicsUniversity of AntwerpAntwerpBelgium
| | - Daniel Nunes
- Champalimaud ResearchChampalimaud Centre for the UnknownLisbonPortugal
| | - Umesh Rudrapatna
- CUBRIC, School of PsychologyCardiff UniversityCardiffUnited Kingdom
| | - Els Fieremans
- Center for Biomedical Imaging, Department of RadiologyNew York University School of MedicineNew YorkUnited States
| | - Derek K Jones
- CUBRIC, School of PsychologyCardiff UniversityCardiffUnited Kingdom
- Mary MacKillop Institute for Health ResearchAustralian Catholic UniversityMelbourneAustralia
| | - Dmitry S Novikov
- Center for Biomedical Imaging, Department of RadiologyNew York University School of MedicineNew YorkUnited States
| | - Noam Shemesh
- Champalimaud ResearchChampalimaud Centre for the UnknownLisbonPortugal
| |
Collapse
|
45
|
Chen K, Lin X, Hu X, Wang J, Zhong H, Jiang L. An enhanced adaptive non-local means algorithm for Rician noise reduction in magnetic resonance brain images. BMC Med Imaging 2020; 20:2. [PMID: 31906873 PMCID: PMC6945655 DOI: 10.1186/s12880-019-0407-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Accepted: 12/27/2019] [Indexed: 11/10/2022] Open
Abstract
Background The Rician noise formed in magnetic resonance (MR) imaging greatly reduced the accuracy and reliability of subsequent analysis, and most of the existing denoising methods are suitable for Gaussian noise rather than Rician noise. Aiming to solve this problem, we proposed fuzzy c-means and adaptive non-local means (FANLM), which combined the adaptive non-local means (NLM) with fuzzy c-means (FCM), as a novel method to reduce noise in the study. Method The algorithm chose the optimal size of search window automatically based on the noise variance which was estimated by the improved estimator of the median absolute deviation (MAD) for Rician noise. Meanwhile, it solved the problem that the traditional NLM algorithm had to use a fixed size of search window. Considering the distribution characteristics for each pixel, we designed three types of search window sizes as large, medium and small instead of using a fixed size. In addition, the combination with the FCM algorithm helped to achieve better denoising effect since the improved the FCM algorithm divided the membership degrees of images and introduced the morphological reconstruction to preserve the image details. Results The experimental results showed that the proposed algorithm (FANLM) can effectively remove the noise. Moreover, it had the highest peak signal-noise ratio (PSNR) and structural similarity (SSIM), compared with other three methods: non-local means (NLM), linear minimum mean square error (LMMSE) and undecimated wavelet transform (UWT). Using the FANLM method, the image details can be well preserved with the noise being mostly removed. Conclusion Compared with the traditional denoising methods, the experimental results showed that the proposed approach effectively suppressed the noise and the edge details were well retained. However, the FANLM method took an average of 13 s throughout the experiment, and its computational cost was not the shortest. Addressing these can be part of our future research.
Collapse
Affiliation(s)
- Kaixin Chen
- Shanghai Key Lab of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, 200093, P, Shanghai, R, China
| | - Xiao Lin
- Shanghai Key Lab of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, 200093, P, Shanghai, R, China
| | - Xing Hu
- Shanghai Key Lab of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, 200093, P, Shanghai, R, China
| | - Jiayao Wang
- Shanghai Key Lab of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, 200093, P, Shanghai, R, China
| | - Han Zhong
- Shanghai Key Lab of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, 200093, P, Shanghai, R, China
| | - Linhua Jiang
- Shanghai Key Lab of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, 200093, P, Shanghai, R, China. .,School of Medicine, Stanford University, 269 Campus Drive, Stanford, CA, 94305, USA.
| |
Collapse
|
46
|
Chen G, Dong B, Zhang Y, Lin W, Yap PT. Denoising of Diffusion MRI Data via Graph Framelet Matching in x-q Space. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:2838-2848. [PMID: 31071025 PMCID: PMC8325050 DOI: 10.1109/tmi.2019.2915629] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Diffusion magnetic resonance imaging (DMRI) suffers from lower signal-to-noise-ratio (SNR) due to MR signal attenuation associated with the motion of water molecules. To improve SNR, the non-local means (NLM) algorithm has demonstrated state-of-the-art performance in noise reduction. However, existing NLM algorithms do not take into account explicitly the fact that DMRI signal can vary significantly with local fiber orientations. Applying NLM naïvely can hence blur subtle structures and aggravate partial volume effects. To overcome this limitation, we improve NLM by performing neighborhood matching in non-flat domains and removing noise with information from both x -space (spatial domain) and q -space (wavevector domain). Specifically, we first encode the q -space sampling domain using a graph. We then perform graph framelet transforms to extract robust rotation-invariant features for each sampling point in x-q space. The resulting features are employed for robust neighborhood matching to locate recurrent information. Finally, we remove noise via an NLM framework. To adapt to the various types of noise in multi-coil MR imaging, we transform the signal before denoising so that it is Gaussian-distributed, allowing noise removal to be carried out in an unbiased manner. Our method is able to more effectively locate recurrent information in white matter structures with different orientations, avoiding the blurring effects caused by naïvely applying NLM. Experiments on synthetic, repetitively-acquired, and infant DMRI data demonstrate that our method is able to preserve subtle structures while effectively removing noise.
Collapse
Affiliation(s)
- Geng Chen
- Department of Radiology and Biomedical Research Imaging Center (BRIC), University of North Carolina at Chapel Hill, NC, U.S.A. D. Shen is also with the Department of Brain and Cognitive Engineering, Korea University, Seoul, Korea
| | - Bin Dong
- Beijing International Center for Mathematical Research, Peking University, Beijing, China
| | - Yong Zhang
- Vancouver Research Center, Huawei, Burnaby, Canada
| | | | | |
Collapse
|
47
|
Cordero-Grande L, Christiaens D, Hutter J, Price AN, Hajnal JV. Complex diffusion-weighted image estimation via matrix recovery under general noise models. Neuroimage 2019; 200:391-404. [PMID: 31226495 PMCID: PMC6711461 DOI: 10.1016/j.neuroimage.2019.06.039] [Citation(s) in RCA: 162] [Impact Index Per Article: 32.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Revised: 03/31/2019] [Accepted: 06/17/2019] [Indexed: 11/28/2022] Open
Abstract
We propose a patch-based singular value shrinkage method for diffusion magnetic resonance image estimation targeted at low signal to noise ratio and accelerated acquisitions. It operates on the complex data resulting from a sensitivity encoding reconstruction, where asymptotically optimal signal recovery guarantees can be attained by modeling the noise propagation in the reconstruction and subsequently simulating or calculating the limit singular value spectrum. Simple strategies are presented to deal with phase inconsistencies and optimize patch construction. The pertinence of our contributions is quantitatively validated on synthetic data, an in vivo adult example, and challenging neonatal and fetal cohorts. Our methodology is compared with related approaches, which generally operate on magnitude-only data and use data-based noise level estimation and singular value truncation. Visual examples are provided to illustrate effectiveness in generating denoised and debiased diffusion estimates with well preserved spatial and diffusion detail.
Collapse
Affiliation(s)
- Lucilio Cordero-Grande
- Centre for the Developing Brain and Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, King's Health Partners, St Thomas' Hospital, London, SE1 7EH, UK.
| | - Daan Christiaens
- Centre for the Developing Brain and Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, King's Health Partners, St Thomas' Hospital, London, SE1 7EH, UK
| | - Jana Hutter
- Centre for the Developing Brain and Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, King's Health Partners, St Thomas' Hospital, London, SE1 7EH, UK
| | - Anthony N Price
- Centre for the Developing Brain and Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, King's Health Partners, St Thomas' Hospital, London, SE1 7EH, UK
| | - Jo V Hajnal
- Centre for the Developing Brain and Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, King's Health Partners, St Thomas' Hospital, London, SE1 7EH, UK
| |
Collapse
|
48
|
Maximov II, Alnæs D, Westlye LT. Towards an optimised processing pipeline for diffusion magnetic resonance imaging data: Effects of artefact corrections on diffusion metrics and their age associations in UK Biobank. Hum Brain Mapp 2019; 40:4146-4162. [PMID: 31173439 PMCID: PMC6865652 DOI: 10.1002/hbm.24691] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2019] [Revised: 05/14/2019] [Accepted: 05/27/2019] [Indexed: 12/30/2022] Open
Abstract
Increasing interest in the structural and functional organisation of the human brain encourages the acquisition of big data sets comprising multiple neuroimaging modalities, often accompanied by additional information obtained from health records, cognitive tests, biomarkers and genotypes. Diffusion weighted magnetic resonance imaging data enables a range of promising imaging phenotypes probing structural connections as well as macroanatomical and microstructural properties of the brain. The reliability and biological sensitivity and specificity of diffusion data depend on processing pipeline. A state-of-the-art framework for data processing facilitates cross-study harmonisation and reduces pipeline-related variability. Using diffusion magnetic resonance imaging (MRI) data from 218 individuals in the UK Biobank, we evaluate the effects of different processing steps that have been suggested to reduce imaging artefacts and improve reliability of diffusion metrics. In lack of a ground truth, we compared diffusion metric sensitivity to age between pipelines. By comparing distributions and age sensitivity of the resulting diffusion metrics based on different approaches (diffusion tensor imaging, diffusion kurtosis imaging and white matter tract integrity), we evaluate a general pipeline comprising seven postprocessing blocks: noise correction; Gibbs ringing correction; evaluation of field distortions; susceptibility, eddy-current and motion-induced distortion corrections; bias field correction; spatial smoothing and final diffusion metric estimations. Based on this evaluation, we suggest an optimised processing pipeline for diffusion weighted MRI data.
Collapse
Affiliation(s)
- Ivan I. Maximov
- Department of PsychologyUniversity of OsloOsloNorway
- Department of Mental Health and AddictionNorwegian Centre for Mental Disorders Research spiepr132 (NORMENT), Oslo University Hospital & Institute of Clinical Medicine, University of OsloOsloNorway
| | - Dag Alnæs
- Department of Mental Health and AddictionNorwegian Centre for Mental Disorders Research spiepr132 (NORMENT), Oslo University Hospital & Institute of Clinical Medicine, University of OsloOsloNorway
| | - Lars T. Westlye
- Department of PsychologyUniversity of OsloOsloNorway
- Department of Mental Health and AddictionNorwegian Centre for Mental Disorders Research spiepr132 (NORMENT), Oslo University Hospital & Institute of Clinical Medicine, University of OsloOsloNorway
| |
Collapse
|
49
|
Barbieri S, Gurney‐Champion OJ, Klaassen R, Thoeny HC. Deep learning how to fit an intravoxel incoherent motion model to diffusion‐weighted MRI. Magn Reson Med 2019; 83:312-321. [DOI: 10.1002/mrm.27910] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2019] [Revised: 06/22/2019] [Accepted: 06/26/2019] [Indexed: 12/29/2022]
Affiliation(s)
| | - Oliver J. Gurney‐Champion
- Joint Department of Physics The Institute of Cancer Research London United Kingdom
- The Royal Marsden NHS Foundation Trust London United Kingdom
| | - Remy Klaassen
- Cancer Center Amsterdam, Department of Medical Oncology and LEXOR (Laboratory for Experimental Oncology and Radiobiology) Academic Medical Center Amsterdam The Netherlands
| | - Harriet C. Thoeny
- Department of Radiology HFR Fribourg‐Hôpital Cantonal Fribourg Switzerland
| |
Collapse
|
50
|
Zhong X, Dale BM, Nickel MD, Kannengiesser SAR, Kiefer B, Bashir M. Improved accuracy of apparent diffusion coefficient quantification using a fully automatic noise bias compensation method: Preliminary evaluation in prostate diffusion weighted imaging. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2019; 305:22-30. [PMID: 31158792 DOI: 10.1016/j.jmr.2019.05.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/15/2019] [Revised: 05/11/2019] [Accepted: 05/20/2019] [Indexed: 06/09/2023]
Abstract
Noise in diffusion magnetic resonance imaging can introduce bias in apparent diffusion coefficient (ADC) quantification. Previous studies proposed methods that are site-specific techniques as research tools with limited availability and typically require manual intervention, not completely ready to use in the clinical environment. The purpose of this study was to develop a fully automatic computational method to correct noise bias in ADC quantification and perform a preliminary evaluation in the clinical prostate diffusion weighted imaging (DWI). Using a pseudo replica approach for the noise map calculation as well as a direct mapping and a stepwise Chebychev polynomial modelling approach for the ADC fitting, a fully automatic noise-bias-compensated ADC calculation method was proposed and implemented both on the scanner and offline. The proposed method was validated in a computer simulation and a standardized diffusion phantom with ground-truth values. Two in vivo studies were performed to evaluate the proposed method in the clinical environment. The first in vivo study performed acquisitions using a clinically routine prostate DWI protocol on 29 subjects to evaluate the consistency between simulated and empirical results. In the second in vivo study, prostate ADC values of 14 subjects were compared between data acquired with external coils only and reconstructed with the proposed method vs. acquired with external combined with endorectal coils and reconstructed with the conventional method. In statistical analyses, p < 0.05 was regarded as significantly different. In the computer simulation, the proposed method showed smaller error percentage than the other methods and was significantly different (p < 2.2 × 10-16). With low signal-to-noise ratio (SNR), the conventional method underestimated ADC values compared to the ground truth values of the diffusion phantom, while the results of the proposed method were more consistent with the ground truth values. Statistical analyses showed no significant differences between measured and simulated results in the first in vivo study (p = 0.5618). Data from the second in vivo study showed that agreement between ADC measured with external coils only and combined coils was improved for the proposed method (mean bias: 0.04 × 10-3 mm2/s, 95% confidence interval (CI) = [-0.01, 0.09] × 10-3 mm2/s, p = 0.187), compared to the conventional method (mean bias: -0.12 × 10-3 mm2/s, 95% CI = [-0.17, -0.06] × 10-3 mm2/s, p < 0.0001). The proposed method compensates noise bias in low-SNR diffusion-weighted acquisitions and results show improved ADC quantification accuracy in the prostate. This method may be suitable for both clinical imaging and research utilizing ADC quantification.
Collapse
Affiliation(s)
- Xiaodong Zhong
- MR R&D Collaborations, Siemens Healthcare, Los Angeles, CA, United States.
| | - Brian M Dale
- MR R&D Collaborations, Siemens Healthcare, Cary, NC, United States
| | - Marcel D Nickel
- MR Application Development, Siemens Healthcare GmbH, Erlangen, Germany
| | | | - Berthold Kiefer
- MR Application Development, Siemens Healthcare GmbH, Erlangen, Germany
| | - Mustafa Bashir
- Department of Radiology, Duke University Medical Center, Durham, NC, United States; Center for Advanced Magnetic Resonance Development, Duke University Medical Center, Durham, NC, United States
| |
Collapse
|