1
|
Fennessy FM, Maier SE. Quantitative diffusion MRI in prostate cancer: Image quality, what we can measure and how it improves clinical assessment. Eur J Radiol 2023; 167:111066. [PMID: 37651828 PMCID: PMC10623580 DOI: 10.1016/j.ejrad.2023.111066] [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/05/2023] [Revised: 08/19/2023] [Accepted: 08/24/2023] [Indexed: 09/02/2023]
Abstract
Diffusion-weighted imaging is a dependable method for detection of clinically significant prostate cancer. In prostate tissue, there are several compartments that can be distinguished from each other, based on different water diffusion decay signals observed. Alterations in cell architecture, such as a relative increase in tumor infiltration and decrease in stroma, will influence the observed diffusion signal in a voxel due to impeded random motion of water molecules. The amount of restricted diffusion can be assessed quantitatively by measuring the apparent diffusion coefficient (ADC) value. This is traditionally calculated using a monoexponential decay formula represented by the slope of a line produced between the logarithm of signal intensity decay plotted against selected b-values. However, the choice and number of b-values and their distribution, has a significant effect on the measured ADC values. There have been many models that attempt to use higher-order functions to better describe the observed diffusion signal decay, requiring an increased number and range of b-values. While ADC can probe heterogeneity on a macroscopic level, there is a need to optimize advanced diffusion techniques to better interrogate prostate tissue microstructure. This could be of benefit in clinical challenges such as identifying sparse tumors in normal prostate tissue or better defining tumor margins. This paper reviews the principles of diffusion MRI and novel higher order diffusion signal analysis techniques to improve the detection of prostate cancer.
Collapse
Affiliation(s)
- Fiona M Fennessy
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, United States.
| | - Stephan E Maier
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, United States; Department of Radiology, Institute of Clinical Sciences, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden
| |
Collapse
|
2
|
Springer CS, Baker EM, Li X, Moloney B, Wilson GJ, Pike MM, Barbara TM, Rooney WD, Maki JH. Metabolic activity diffusion imaging (MADI): I. Metabolic, cytometric modeling and simulations. NMR IN BIOMEDICINE 2023; 36:e4781. [PMID: 35654608 DOI: 10.1002/nbm.4781] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 05/25/2022] [Accepted: 05/26/2022] [Indexed: 06/15/2023]
Abstract
Evidence mounts that the steady-state cellular water efflux (unidirectional) first-order rate constant (kio [s-1 ]) magnitude reflects the ongoing, cellular metabolic rate of the cytolemmal Na+ , K+ -ATPase (NKA), c MRNKA (pmol [ATP consumed by NKA]/s/cell), perhaps biology's most vital enzyme. Optimal 1 H2 O MR kio determinations require paramagnetic contrast agents (CAs) in model systems. However, results suggest that the homeostatic metabolic kio biomarker magnitude in vivo is often too large to be reached with allowable or possible CA living tissue distributions. Thus, we seek a noninvasive (CA-free) method to determine kio in vivo. Because membrane water permeability has long been considered important in tissue water diffusion, we turn to the well-known diffusion-weighted MRI (DWI) modality. To analyze the diffusion tensor magnitude, we use a parsimoniously primitive model featuring Monte Carlo simulations of water diffusion in virtual ensembles comprising water-filled and -immersed randomly sized/shaped contracted Voronoi cells. We find this requires two additional, cytometric properties: the mean cell volume (V [pL]) and the cell number density (ρ [cells/μL]), important biomarkers in their own right. We call this approach metabolic activity diffusion imaging (MADI). We simulate water molecule displacements and transverse MR signal decays covering the entirety of b-space from pure water (ρ = V = 0; kio undefined; diffusion coefficient, D0 ) to zero diffusion. The MADI model confirms that, in compartmented spaces with semipermeable boundaries, diffusion cannot be described as Gaussian: the nanoscopic D (Dn ) is diffusion time-dependent, a manifestation of the "diffusion dispersion". When the "well-mixed" (steady-state) condition is reached, diffusion becomes limited, mainly by the probabilities of (1) encountering (ρ, V), and (2) permeating (kio ) cytoplasmic membranes, and less so by Dn magnitudes. Importantly, for spaces with large area/volume (A/V; claustrophobia) ratios, this can happen in less than a millisecond. The model matches literature experimental data well, with implications for DWI interpretations.
Collapse
Affiliation(s)
- Charles S Springer
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
- Knight Cancer Institute, Oregon Health & Science University, Portland, Oregon, USA
- Brenden-Colson Center for Pancreatic Care, Oregon Health & Science University, Portland, Oregon, USA
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
- Department of Chemical Physiology and Biochemistry, Oregon Health & Science University, Portland, Oregon, USA
| | - Eric M Baker
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
| | - Xin Li
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
- Brenden-Colson Center for Pancreatic Care, Oregon Health & Science University, Portland, Oregon, USA
| | - Brendan Moloney
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
| | - Gregory J Wilson
- Department of Radiology, University of Washington, Seattle, Washington, USA
- Bayer Healthcare, Radiology, New Jersey, USA
| | - Martin M Pike
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
| | - Thomas M Barbara
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
| | - William D Rooney
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon, USA
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, Oregon, USA
- Department of Behavioral Neuroscience, Oregon Health & Science University, Portland, Oregon, USA
| | - Jeffrey H Maki
- Anschutz Medical Center Department of Radiology, University of Colorado, Aurora, Colorado, USA
| |
Collapse
|
3
|
Li X, Mangia S, Lee JH, Bai R, Springer CS. NMR shutter-speed elucidates apparent population inversion of 1 H 2 O signals due to active transmembrane water cycling. Magn Reson Med 2019; 82:411-424. [PMID: 30903632 PMCID: PMC6593680 DOI: 10.1002/mrm.27725] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Revised: 02/01/2019] [Accepted: 02/11/2019] [Indexed: 12/13/2022]
Abstract
Purpose The desire to quantitatively discriminate the extra‐ and intracellular tissue 1H2O MR signals has gone hand‐in‐hand with the continual, historic increase in MRI instrument magnetic field strength [B0]. However, recent studies have indicated extremely valuable, novel metabolic information can be readily accessible at ultra–low B0. The two signals can be distinguished, and the homeostatic activity of the cell membrane sodium/potassium pump (Na+,K+,ATPase) detected. The mechanism allowing 1H2O MRI to do this is the newly discovered active transmembrane water cycling (AWC) phenomenon, which we found using paramagnetic extracellular contrast agents at clinical B0 values. AWC is important because Na+,K+,ATPase can be considered biology’s most vital enzyme, and its in vivo steady‐state activity has not before been measurable, let alone amenable to mapping with high spatial resolution. Recent reports indicate AWC correlates with neuronal firing rate, with malignant tumor metastatic potential, and inversely with cellular reducing equivalent fraction. We wish to systematize the ways AWC can be precisely measured. Methods We present a theoretical longitudinal relaxation analysis of considerable scope: it spans the low‐ and high–field situations. Results We show the NMR shutter‐speed organizing principle is pivotal in understanding how trans–membrane steady–state water exchange kinetics are manifest throughout the range. Our findings illuminate an aspect, apparent population inversion, which is crucial in understanding ultra‐low field results. Conclusions Without an appreciation of apparent population inversion, significant misinterpretations of future data are likely. These could have unfortunate diagnostic consequences.
Collapse
Affiliation(s)
- Xin Li
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon
| | - Silvia Mangia
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota
| | - Jing-Huei Lee
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, Ohio
| | - Ruiliang Bai
- Interdisciplinary Institute of Neuroscience and Technology, Qiushi Academy for Advanced Studies, Key Laboratory of Biomedical Engineering of Ministry of Education, College of Biomedical Engineering and Instrument Science, Zhejiang University, Hangzhou, China
| | - Charles S Springer
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, Oregon
| |
Collapse
|
4
|
Reci A, Ainte MI, Sederman AJ, Mantle MD, Gladden LF. Optimising sampling patterns for bi-exponentially decaying signals. Magn Reson Imaging 2018; 56:14-18. [PMID: 30413334 DOI: 10.1016/j.mri.2018.09.029] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2018] [Revised: 09/25/2018] [Accepted: 09/26/2018] [Indexed: 11/29/2022]
Abstract
A recently reported method, based on the Cramér-Rao Lower Bound theory, for optimising sampling patterns for a wide range of nuclear magnetic resonance (NMR) experiments is applied to the problem of optimising sampling patterns for bi-exponentially decaying signals. Sampling patterns are optimised by minimizing the percentage error in estimating the most difficult to estimate parameter of the bi-exponential model, termed the objective function. The predictions of the method are demonstrated in application to pulsed field gradient NMR data recorded for the two-component diffusion of a binary mixture of methane/ethane in a zeolite. It is shown that the proposed method identifies an optimal sampling pattern with the predicted objective function being within 10% of that calculated from the experiment dataset. The method is used to advise on the number of sampled points and the noise level needed to resolve two-component systems characterised by a range of ratios of populations and diffusion coefficients. It is subsequently illustrated how the method can be used to reduce the experiment acquisition time while still being able to resolve a given two-component system.
Collapse
Affiliation(s)
- A Reci
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - M I Ainte
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - A J Sederman
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom.
| | - M D Mantle
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - L F Gladden
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| |
Collapse
|
5
|
Nilsson M, Englund E, Szczepankiewicz F, van Westen D, Sundgren PC. Imaging brain tumour microstructure. Neuroimage 2018; 182:232-250. [DOI: 10.1016/j.neuroimage.2018.04.075] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Revised: 04/27/2018] [Accepted: 04/30/2018] [Indexed: 01/18/2023] Open
|
6
|
Reci A, Sederman AJ, Gladden LF. Optimising magnetic resonance sampling patterns for parametric characterisation. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2018; 294:35-43. [PMID: 30005192 DOI: 10.1016/j.jmr.2018.06.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2018] [Revised: 05/29/2018] [Accepted: 06/18/2018] [Indexed: 06/08/2023]
Abstract
Sampling strategies are often central to experimental design. Choosing efficiently which data to acquire can improve the estimation of parameters and reduce the acquisition time. This work is focused on designing optimal sampling patterns for Nuclear Magnetic Resonance (NMR) applications, illustrated with respect to the best estimate of the parameters characterising a lognormal distribution. Lognormal distributions are commonly used as fitting models for distributions of spin-lattice relaxation time constants, spin-spin relaxation time constants and diffusion coefficients. A method for optimising the choice of points to be sampled is presented which is based on the Cramér-Rao Lower Bound (CRLB) theory. The method's capabilities are demonstrated experimentally by applying it to the problem of estimating the emulsion droplet size distribution from a pulsed field gradient (PFG) NMR diffusion experiment. A difference of <5% is observed between the predictions of CRLB theory and the PFG NMR experimental results. It is shown that CLRB theory is stable down to signal-to-noise ratios of ∼10. A sensitivity analysis for the CRLB theory is also performed. The method of optimizing sampling patterns is easily adapted to distributions other than lognormal and to other aspects of experimental design; case studies of optimising the sampling scheme for a fixed acquisition time and determining the potential for reduction in acquisition time for a fixed parameter estimation accuracy are presented. The experimental acquisition time is typically reduced by a factor of 3 using the proposed method compared to a constant gradient increment approach that would usually be used.
Collapse
Affiliation(s)
- A Reci
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| | - A J Sederman
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom.
| | - L F Gladden
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| |
Collapse
|
7
|
Springer CS. Using 1H 2O MR to measure and map sodium pump activity in vivo. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2018; 291:110-126. [PMID: 29705043 DOI: 10.1016/j.jmr.2018.02.018] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2017] [Revised: 02/16/2018] [Accepted: 02/26/2018] [Indexed: 05/26/2023]
Abstract
The cell plasma membrane Na+,K+-ATPase [NKA] is one of biology's most [if not the most] significant enzymes. By actively transporting Na+ out [and K+ in], it maintains the vital trans-membrane ion concentration gradients and the membrane potential. The forward NKA reaction is shown in the Graphical Abstract [which is elaborated in the text]. Crucially, NKA does not operate in isolation. There are other transporters that conduct K+ back out of [II, Graphical Abstract] and Na+ back into [III, Graphical Abstract] the cell. Thus, NKA must function continually. Principal routes for ATP replenishment include mitochondrial oxidative phosphorylation, glycolysis, and creatine kinase [CrK] activity. However, it has never been possible to measure, let alone map, this integrated, cellular homeostatic NKA activity in vivo. Active trans-membrane water cycling [AWC] promises a way to do this with 1H2O MR. Inthe Graphical Abstract, the AWC system is characterized by active contributions totheunidirectional rate constants for steady-state water efflux and influx, respectively, kio(a) and koi(a). The discovery, validation, and initial exploration of active water cycling are reviewed here. Promising applications in cancer, cardiological, and neurological MRI are covered. This initial work employed paramagnetic Gd(III)chelate contrast agents [CAs]. However, the significant problems associated with in vivo CA use are also reviewed. A new analysis of water diffusion-weighted MRI [DWI] is presented. Preliminary results suggest a non-invasive way to measure the cell number density [ρ (cells/μL)], the mean cell volume [V (pL)], and the cellular NKA metabolic rate [cMRNKA(fmol(ATP)/s/cell)] with high spatial resolution. These crucial cell biology properties have not before been accessible invivo. Furthermore, initial findings indicate their absolute values can be determined.
Collapse
Affiliation(s)
- Charles S Springer
- Advanced Imaging Research Center, Oregon Health & Science University, Portland, OR 97239, United States.
| |
Collapse
|
8
|
Kim D, Doyle EK, Wisnowski JL, Kim JH, Haldar JP. Diffusion-relaxation correlation spectroscopic imaging: A multidimensional approach for probing microstructure. Magn Reson Med 2017; 78:2236-2249. [PMID: 28317261 PMCID: PMC5605406 DOI: 10.1002/mrm.26629] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2016] [Revised: 12/19/2016] [Accepted: 01/10/2017] [Indexed: 12/16/2022]
Abstract
PURPOSE To propose and evaluate a novel multidimensional approach for imaging subvoxel tissue compartments called Diffusion-Relaxation Correlation Spectroscopic Imaging. THEORY AND METHODS Multiexponential modeling of MR diffusion or relaxation data is commonly used to infer the many different microscopic tissue compartments that contribute signal to macroscopic MR imaging voxels. However, multiexponential estimation is known to be difficult and ill-posed. Observing that this ill-posedness is theoretically reduced in higher dimensions, diffusion-relaxation correlation spectroscopic imaging uses a novel multidimensional imaging experiment that jointly encodes diffusion and relaxation information, and then uses a novel constrained reconstruction technique to generate a multidimensional diffusion-relaxation correlation spectrum for every voxel. The peaks of the multidimensional spectrum are expected to correspond to the distinct tissue microenvironments that are present within each macroscopic imaging voxel. RESULTS Using numerical simulations, experiment data from a custom-built phantom, and experiment data from a mouse model of traumatic spinal cord injury, diffusion-relaxation correlation spectroscopic imaging is demonstrated to provide substantially better multicompartment resolving power compared to conventional diffusion- and relaxation-based methods. CONCLUSION The diffusion-relaxation correlation spectroscopic imaging approach provides powerful new capabilities for resolving the different components of multicompartment tissue models, and can be leveraged to significantly expand the insights provided by MRI in studies of tissue microstructure. Magn Reson Med 78:2236-2249, 2017. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Daeun Kim
- Electrical Engineering, University of Southern California, Los Angeles, CA, USA
| | - Eamon K. Doyle
- Biomedical Engineering, University of Southern California, Los Angeles, CA, USA
- Cardiology, Children’s Hospital Los Angeles, Los Angeles, CA, USA
| | | | - Joong Hee Kim
- Neurology and Radiology, Washington University, St. Louis, MO, USA
| | - Justin P. Haldar
- Electrical Engineering, University of Southern California, Los Angeles, CA, USA
- Biomedical Engineering, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
9
|
Langkilde F, Kobus T, Fedorov A, Dunne R, Tempany C, Mulkern RV, Maier SE. Evaluation of fitting models for prostate tissue characterization using extended-range b-factor diffusion-weighted imaging. Magn Reson Med 2017; 79:2346-2358. [PMID: 28718517 DOI: 10.1002/mrm.26831] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2017] [Revised: 06/16/2017] [Accepted: 06/19/2017] [Indexed: 11/08/2022]
Abstract
PURPOSE To compare the fitting and tissue discrimination performance of biexponential, kurtosis, stretched exponential, and gamma distribution models for high b-factor diffusion-weighted images in prostate cancer. METHODS Diffusion-weighted images with 15 b-factors ranging from b = 0 to 3500 s/mm2 were obtained in 62 prostate cancer patients. Pixel-wise signal decay fits for each model were evaluated with the Akaike Information Criterion (AIC). Parameter values for each model were determined within normal prostate and the index lesion. Their potential to differentiate normal from cancerous tissue was investigated through receiver operating characteristic analysis and comparison with Gleason score. RESULTS The biexponential slow diffusion fraction fslow , the apparent kurtosis diffusion coefficient ADCK , and the excess kurtosis factor K differ significantly among normal peripheral zone (PZ), normal transition zone (TZ), tumor PZ, and tumor TZ. Biexponential and gamma distribution models result in the lowest AIC, indicating a superior fit. Maximum areas under the curve (AUCs) of all models ranged from 0.93 to 0.96 for the PZ and from 0.95 to 0.97 for the TZ. Similar AUCs also result from the apparent diffusion coefficient (ADC) of a monoexponential fit to a b-factor sub-range up to 1250 s/mm2 . For kurtosis and stretched exponential models, single parameters yield the highest AUCs, whereas for the biexponential and gamma distribution models, linear combinations of parameters produce the highest AUCs. Parameters with high AUC show a trend in differentiating low from high Gleason score, whereas parameters with low AUC show no such ability. CONCLUSION All models, including a monoexponential fit to a lower-b sub-range, achieve similar AUCs for discrimination of normal and cancer tissue. The biexponential model, which is favored statistically, also appears to provide insight into disease-related microstructural changes. Magn Reson Med 79:2346-2358, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Fredrik Langkilde
- Institute of Clinical Sciences, Sahlgrenska Academy, Gothenburg University, Gothenburg, Sweden
| | - Thiele Kobus
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA.,Department of Radiology and Nuclear Medicine, Radboud university medical center, Nijmegen, The Netherlands
| | - Andriy Fedorov
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Ruth Dunne
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Clare Tempany
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Robert V Mulkern
- Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA.,Department of Radiology, Boston Children's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Stephan E Maier
- Institute of Clinical Sciences, Sahlgrenska Academy, Gothenburg University, Gothenburg, Sweden.,Department of Radiology, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, USA
| |
Collapse
|