1
|
Canales-Rodríguez EJ, Pizzolato M, Zhou FL, Barakovic M, Thiran JP, Jones DK, Parker GJM, Dyrby TB. Pore size estimation in axon-mimicking microfibers with diffusion-relaxation MRI. Magn Reson Med 2024; 91:2579-2596. [PMID: 38192108 DOI: 10.1002/mrm.29991] [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: 10/17/2023] [Revised: 12/04/2023] [Accepted: 12/12/2023] [Indexed: 01/10/2024]
Abstract
PURPOSE This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T2-based pore size estimation technique. THEORY AND METHODS A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images. RESULTS Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T2-based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes. CONCLUSION The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T2-based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research.
Collapse
Affiliation(s)
- Erick J Canales-Rodríguez
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| | - Marco Pizzolato
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
- Department of Applied Mathematics and Computer Science, Technical University of Denmark (DTU), Kongens Lyngby, Denmark
| | - Feng-Lei Zhou
- Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London (UCL), London, UK
- MicroPhantoms Limited, Cambridge, UK
| | - Muhamed Barakovic
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
- Centre d'Imagerie Biomédicale (CIBM), EPFL, Lausanne, Switzerland
| | - Derek K Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, UK
| | - Geoffrey J M Parker
- Centre for Medical Image Computing, Department of Medical Physics and Biomedical Engineering, University College London (UCL), London, UK
- Department of Neuroinflammation, Queen Square Institute of Neurology, University College London (UCL), London, UK
- Bioxydyn Limited, Manchester, UK
| | - Tim B Dyrby
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
- Department of Applied Mathematics and Computer Science, Technical University of Denmark (DTU), Kongens Lyngby, Denmark
| |
Collapse
|
2
|
Mallick H, Porwal A, Saha S, Basak P, Svetnik V, Paul E. An integrated Bayesian framework for multi-omics prediction and classification. Stat Med 2024; 43:983-1002. [PMID: 38146838 DOI: 10.1002/sim.9953] [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: 04/21/2022] [Revised: 10/06/2023] [Accepted: 10/24/2023] [Indexed: 12/27/2023]
Abstract
With the growing commonality of multi-omics datasets, there is now increasing evidence that integrated omics profiles lead to more efficient discovery of clinically actionable biomarkers that enable better disease outcome prediction and patient stratification. Several methods exist to perform host phenotype prediction from cross-sectional, single-omics data modalities but decentralized frameworks that jointly analyze multiple time-dependent omics data to highlight the integrative and dynamic impact of repeatedly measured biomarkers are currently limited. In this article, we propose a novel Bayesian ensemble method to consolidate prediction by combining information across several longitudinal and cross-sectional omics data layers. Unlike existing frequentist paradigms, our approach enables uncertainty quantification in prediction as well as interval estimation for a variety of quantities of interest based on posterior summaries. We apply our method to four published multi-omics datasets and demonstrate that it recapitulates known biology in addition to providing novel insights while also outperforming existing methods in estimation, prediction, and uncertainty quantification. Our open-source software is publicly available at https://github.com/himelmallick/IntegratedLearner.
Collapse
Affiliation(s)
- Himel Mallick
- Division of Biostatistics, Department of Population Health Sciences, Weill Cornell Medicine, Cornell University, New York, 10065, New York, USA
- Department of Statistics and Data Science, Cornell University, Ithaca, New York, USA
| | - Anupreet Porwal
- Department of Statistics, University of Washington, Seattle, Washington, USA
| | - Satabdi Saha
- Department of Biostatistics, University of Texas MD Anderson Cancer Center, Houston, Texas, USA
| | - Piyali Basak
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| | - Vladimir Svetnik
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| | - Erina Paul
- Biostatistics and Research Decision Sciences, Merck & Co., Inc., Rahway, New Jersey, USA
| |
Collapse
|
3
|
Barakovic M, Pizzolato M, Tax CMW, Rudrapatna U, Magon S, Dyrby TB, Granziera C, Thiran JP, Jones DK, Canales-Rodríguez EJ. Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology. Front Neurosci 2023; 17:1209521. [PMID: 37638307 PMCID: PMC10457121 DOI: 10.3389/fnins.2023.1209521] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Accepted: 07/24/2023] [Indexed: 08/29/2023] Open
Abstract
Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement.
Collapse
Affiliation(s)
- Muhamed Barakovic
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
- Department of Neurology, University Hospital Basel, Basel, Switzerland
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Roche Pharma Research and Early Development, Neuroscience and Rare Diseases, Roche Innovation Center, Basel, Switzerland
| | - Marco Pizzolato
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Chantal M. W. Tax
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
- Image Sciences Institute, University Medical Center Utrecht, Utrecht, Netherlands
| | - Umesh Rudrapatna
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
| | - Stefano Magon
- Roche Pharma Research and Early Development, Neuroscience and Rare Diseases, Roche Innovation Center, Basel, Switzerland
| | - Tim B. Dyrby
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
- Danish Research Centre for Magnetic Resonance (DRCMR), Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| | - Cristina Granziera
- Translational Imaging in Neurology (ThINk) Basel, Department of Biomedical Engineering, University Hospital Basel and University of Basel, Basel, Switzerland
- Department of Neurology, University Hospital Basel, Basel, Switzerland
- Research Center for Clinical Neuroimmunology and Neuroscience Basel (RC2NB), University Hospital Basel and University of Basel, Basel, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Radiology Department, Centre Hospitalier Universitaire Vaudois and University of Lausanne, Lausanne, Switzerland
- Centre d’Imagerie Biomédicale (CIBM), EPFL, Lausanne, Switzerland
| | - Derek K. Jones
- Cardiff University Brain Research Imaging Centre, Cardiff University, Cardiff, Wales, United Kingdom
| | - Erick J. Canales-Rodríguez
- Signal Processing Laboratory 5 (LTS5), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
4
|
Endt S, Engel M, Naldi E, Assereto R, Molendowska M, Mueller L, Mayrink Verdun C, Pirkl CM, Palombo M, Jones DK, Menzel MI. In Vivo Myelin Water Quantification Using Diffusion-Relaxation Correlation MRI: A Comparison of 1D and 2D Methods. APPLIED MAGNETIC RESONANCE 2023; 54:1571-1588. [PMID: 38037641 PMCID: PMC10682074 DOI: 10.1007/s00723-023-01584-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Revised: 07/18/2023] [Accepted: 07/25/2023] [Indexed: 12/02/2023]
Abstract
Multidimensional Magnetic Resonance Imaging (MRI) is a versatile tool for microstructure mapping. We use a diffusion weighted inversion recovery spin echo (DW-IR-SE) sequence with spiral readouts at ultra-strong gradients to acquire a rich diffusion-relaxation data set with sensitivity to myelin water. We reconstruct 1D and 2D spectra with a two-step convex optimization approach and investigate a variety of multidimensional MRI methods, including 1D multi-component relaxometry, 1D multi-component diffusometry, 2D relaxation correlation imaging, and 2D diffusion-relaxation correlation spectroscopic imaging (DR-CSI), in terms of their potential to quantify tissue microstructure, including the myelin water fraction (MWF). We observe a distinct spectral peak that we attribute to myelin water in multi-component T1 relaxometry, T1-T2 correlation, T1-D correlation, and T2-D correlation imaging. Due to lower achievable echo times compared to diffusometry, MWF maps from relaxometry have higher quality. Whilst 1D multi-component T1 data allows much faster myelin mapping, 2D approaches could offer unique insights into tissue microstructure and especially myelin diffusion.
Collapse
Affiliation(s)
- Sebastian Endt
- Technical University of Munich, Munich, Germany
- AImotion Bavaria, Technische Hochschule Ingolstadt, Ingolstadt, Germany
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
| | - Maria Engel
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
| | - Emanuele Naldi
- Technische Universität Braunschweig, Braunschweig, Germany
| | | | - Malwina Molendowska
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
- Medical Radiation Physics, Lund University, Lund, Sweden
| | - Lars Mueller
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
- Biomedical Imaging Science Department, Leeds Institute of Cardiovascular and Metabolic Medicine (LICAMM), University of Leeds, Leeds, United Kingdom
| | - Claudio Mayrink Verdun
- Technical University of Munich, Munich, Germany
- Munich Center for Machine Learning, Munich, Germany
| | | | - Marco Palombo
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
| | - Derek K. Jones
- Cardiff University Brain Research Imaging Centre (CUBRIC), Cardiff University, Cardiff, United Kingdom
| | - Marion I. Menzel
- Technical University of Munich, Munich, Germany
- AImotion Bavaria, Technische Hochschule Ingolstadt, Ingolstadt, Germany
- GE HealthCare, Munich, Germany
| |
Collapse
|
5
|
Ben-Atya H, Freiman M. P 2T 2: A physically-primed deep-neural-network approach for robust T 2 distribution estimation from quantitative T 2-weighted MRI. Comput Med Imaging Graph 2023; 107:102240. [PMID: 37224742 DOI: 10.1016/j.compmedimag.2023.102240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2023] [Revised: 04/27/2023] [Accepted: 04/27/2023] [Indexed: 05/26/2023]
Abstract
Estimating T2 relaxation time distributions from multi-echo T2-weighted MRI (T2W) data can provide valuable biomarkers for assessing inflammation, demyelination, edema, and cartilage composition in various pathologies, including neurodegenerative disorders, osteoarthritis, and tumors. Deep neural network (DNN) based methods have been proposed to address the complex inverse problem of estimating T2 distributions from MRI data, but they are not yet robust enough for clinical data with low Signal-to-Noise ratio (SNR) and are highly sensitive to distribution shifts such as variations in echo-times (TE) used during acquisition. Consequently, their application is hindered in clinical practice and large-scale multi-institutional trials with heterogeneous acquisition protocols. We propose a physically-primed DNN approach, called P2T2, that incorporates the signal decay forward model in addition to the MRI signal into the DNN architecture to improve the accuracy and robustness of T2 distribution estimation. We evaluated our P2T2 model in comparison to both DNN-based methods and classical methods for T2 distribution estimation using 1D and 2D numerical simulations along with clinical data. Our model improved the baseline model's accuracy for low SNR levels (SNR<80) which are common in the clinical setting. Further, our model achieved a ∼35% improvement in robustness against distribution shifts in the acquisition process compared to previously proposed DNN models. Finally, Our P2T2 model produces the most detailed Myelin-Water fraction maps compared to baseline approaches when applied to real human MRI data. Our P2T2 model offers a reliable and precise means of estimating T2 distributions from MRI data and shows promise for use in large-scale multi-institutional trials with heterogeneous acquisition protocols. Our source code is available at: https://github.com/Hben-atya/P2T2-Robust-T2-estimation.git.
Collapse
Affiliation(s)
- Hadas Ben-Atya
- Faculty of Biomedical Engineering, Technion - Israel Institute of Technology, Haifa, Israel.
| | - Moti Freiman
- Faculty of Biomedical Engineering, Technion - Israel Institute of Technology, Haifa, Israel
| |
Collapse
|
6
|
Pizzolato M, Canales-Rodríguez EJ, Andersson M, Dyrby TB. Axial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI. Med Image Anal 2023; 86:102767. [PMID: 36867913 DOI: 10.1016/j.media.2023.102767] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 12/13/2022] [Accepted: 02/08/2023] [Indexed: 02/18/2023]
Abstract
We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities.
Collapse
Affiliation(s)
- Marco Pizzolato
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kgs. Lyngby, Denmark; Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark.
| | | | - Mariam Andersson
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| | - Tim B Dyrby
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kgs. Lyngby, Denmark; Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Copenhagen, Denmark
| |
Collapse
|
7
|
Drenthen GS, Jansen JFA, van der MM, Voorter PHM, Backes WH. An optimized b-value sampling for the quantification of interstitial fluid using diffusion-weighted MRI, a genetic algorithm approach. Magn Reson Med 2023; 90:194-201. [PMID: 36744716 DOI: 10.1002/mrm.29612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 01/17/2023] [Accepted: 01/21/2023] [Indexed: 02/07/2023]
Abstract
PURPOSE Multi-b-value diffusion-weighted MRI techniques can simultaneously measure the parenchymal diffusivity, microvascular perfusion, and a third, intermediate diffusion component. This component is related to the interstitial fluid in the brain parenchyma. However, simultaneously estimating three diffusion components from multi-b-value data is difficult and has strong dependence on SNR and chosen b-values. As the number of acquired b-values is limited due to scanning time, it is important to know which b-values are most effective to be included. Therefore, this study evaluates an optimized b-value sampling for interstitial fluid estimation. METHOD The optimized b-value sampling scheme is determined using a genetic algorithm. Subsequently, the performance of this optimized sampling is assessed by comparing it with a linear, logarithmic, and previously proposed sampling scheme, in terms of the RMS error (RMSE) for the intermediate component estimation. The in vivo performance of the optimized sampling is assessed using 7T data with 101 equally spaced b-values ranging from 0 to 1000 s/mm2 . In this case, the RMSE was determined by comparing the fit that includes all b-values. RESULTS The optimized b-value sampling for estimating the intermediate component was reported to be [0, 30, 90, 210, 280, 350, 580, 620, 660, 680, 720, 760, 980, 990, 1000] s/mm2 . For computer simulations, the optimized sampling had a lower RMSE, compared with the other samplings for varying levels of SNR. For the in vivo data, the voxel-wise RMSE of the optimized sampling was lower compared with other sampling schemes. CONCLUSION The genetic algorithm-optimized b-value scheme improves the quantification of the diffusion component related to interstitial fluid in terms of a lower RMSE.
Collapse
Affiliation(s)
- Gerhard S Drenthen
- Department of Radiology and Nuclear Medicine, Maastricht University Medical Center, Maastricht, The Netherlands.,School for Mental Health and Neuroscience, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Jacobus F A Jansen
- Department of Radiology and Nuclear Medicine, Maastricht University Medical Center, Maastricht, The Netherlands.,School for Mental Health and Neuroscience, Maastricht University Medical Center, Maastricht, The Netherlands.,Department of Electrical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Merel M van der
- Department of Radiology and Nuclear Medicine, Maastricht University Medical Center, Maastricht, The Netherlands.,School for Mental Health and Neuroscience, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Paulien H M Voorter
- Department of Radiology and Nuclear Medicine, Maastricht University Medical Center, Maastricht, The Netherlands.,School for Mental Health and Neuroscience, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Walter H Backes
- Department of Radiology and Nuclear Medicine, Maastricht University Medical Center, Maastricht, The Netherlands.,School for Mental Health and Neuroscience, Maastricht University Medical Center, Maastricht, The Netherlands
| |
Collapse
|
8
|
Fischi-Gomez E, Girard G, Koch PJ, Yu T, Pizzolato M, Brügger J, Piredda GF, Hilbert T, Cadic-Melchior AG, Beanato E, Park CH, Morishita T, Wessel MJ, Schiavi S, Daducci A, Kober T, Canales-Rodríguez EJ, Hummel FC, Thiran JP. Variability and reproducibility of multi-echo T2 relaxometry: Insights from multi-site, multi-session and multi-subject MRI acquisitions. FRONTIERS IN RADIOLOGY 2022; 2:930666. [PMID: 37492668 PMCID: PMC10365099 DOI: 10.3389/fradi.2022.930666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/28/2022] [Accepted: 06/30/2022] [Indexed: 07/27/2023]
Abstract
Quantitative magnetic resonance imaging (qMRI) can increase the specificity and sensitivity of conventional weighted MRI to underlying pathology by comparing meaningful physical or chemical parameters, measured in physical units, with normative values acquired in a healthy population. This study focuses on multi-echo T2 relaxometry, a qMRI technique that probes the complex tissue microstructure by differentiating compartment-specific T2 relaxation times. However, estimation methods are still limited by their sensitivity to the underlying noise. Moreover, estimating the model's parameters is challenging because the resulting inverse problem is ill-posed, requiring advanced numerical regularization techniques. As a result, the estimates from distinct regularization strategies are different. In this work, we aimed to investigate the variability and reproducibility of different techniques for estimating the transverse relaxation time of the intra- and extra-cellular space (T2IE) in gray (GM) and white matter (WM) tissue in a clinical setting, using a multi-site, multi-session, and multi-run T2 relaxometry dataset. To this end, we evaluated three different techniques for estimating the T2 spectra (two regularized non-negative least squares methods and a machine learning approach). Two independent analyses were performed to study the effect of using raw and denoised data. For both the GM and WM regions, and the raw and denoised data, our results suggest that the principal source of variance is the inter-subject variability, showing a higher coefficient of variation (CoV) than those estimated for the inter-site, inter-session, and inter-run, respectively. For all reconstruction methods studied, the CoV ranged between 0.32 and 1.64%. Interestingly, the inter-session variability was close to the inter-scanner variability with no statistical differences, suggesting that T2IE is a robust parameter that could be employed in multi-site neuroimaging studies. Furthermore, the three tested methods showed consistent results and similar intra-class correlation (ICC), with values superior to 0.7 for most regions. Results from raw data were slightly more reproducible than those from denoised data. The regularized non-negative least squares method based on the L-curve technique produced the best results, with ICC values ranging from 0.72 to 0.92.
Collapse
Affiliation(s)
- Elda Fischi-Gomez
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Translational Machine Learning Lab, Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| | - Gabriel Girard
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| | - Philipp J. Koch
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
- Department of Neurology, University of Lübeck, Lübeck, Germany
- Center of Brain, Behavior and Metabolism (CBBM), University of Lübeck, Lübeck, Germany
| | - Thomas Yu
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland
| | - Marco Pizzolato
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Julia Brügger
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
| | - Gian Franco Piredda
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland
| | - Tom Hilbert
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland
| | - Andéol G. Cadic-Melchior
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
| | - Elena Beanato
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
| | - Chang-Hyun Park
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
| | - Takuya Morishita
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
| | - Maximilian J. Wessel
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
- Department of Neurology, University Hospital and Julius-Maximilians-University, Wuerzburg, Germany
| | - Simona Schiavi
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Health (DINOGMI), University of Genoa, Genoa, Italy
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy
| | - Alessandro Daducci
- Diffusion Imaging and Connectivity Estimation (DICE) Lab, Department of Computer Science, University of Verona, Verona, Italy
| | - Tobias Kober
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
- Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland
| | - Erick J. Canales-Rodríguez
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Friedhelm C. Hummel
- Defitech Chair for Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Defitech Chair of Clinical Neuroengineering, Neuro-X Institute (NIX) and Brain Mind Institute (BMI), École Polytechnique Fédérale de Lausanne (EPFL Valais), Clinique Romande de Réadaptation, Sion, Switzerland
- Clinical Neuroscience, University Hospital of Geneva (HUG), Geneva, Switzerland
| | - Jean-Philippe Thiran
- Signal Processing Laboratory 5 (LTS5), École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- CIBM Center for Biomedical Imaging, Lausanne, Switzerland
- Department of Radiology, Centre Hospitalier Universitaire Vaudois, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|