Synopsis
Quantitative imaging metrics represent an important area of development in our field, but also present a number of challenges in bringing them to fruition. Several advanced diffusion MRI methods will be discussed and prospective or candidate quantitative imaging biomarkers will be described.
What
is a qualitative imaging?
The
term “quantitative imaging” may seem like an oxymoron to some or at least an
unattainable goal. Neuroradiological practice often entails making a collective
assessment of a set of qualitative “weighted” MR images and providing a concise
narrative explaining their clinical significance. However, T1- T2-,
D-, … weighed images are susceptible to a myriad of experimental confounds and
artifacts, and often cannot be quantitatively compared in longitudinal and
particularly, in multi-site studies.
What
is the goal of quantitative imaging?
A
goal of quantitative MRI is to transform the radiologist’s MRI scanner into a scientific
instrument that produces robust, accurate, precise, and stable measurements of
quantities or parameters having a physical meaning, and hopefully a concomitant
biological significance that can improve diagnostic sensitivity and/or
specificity, enabling a more informed clinical assessment.
What
characteristics and properties are shared by quantitative imaging metrics?
Quantitative
imaging parameters should share several common characteristics or features.
They should either be intrinsic physical quantities or constants, which are
defined in terms of a well-accepted physical or chemical model; or parameters,
which represent a meaningful ratio of physical quantities or constants that
measure an intrinsic feature or quality of a physical process. Image-derived physical
quantities should be measured in physical units, robustly, accurately and
precisely. One should be able to relate these quantities to a measured MR
signal using well-established physical-mathematical models. It is prudent for some
scalar parameters to be unitless or dimensionless. For directional processes,
like anisotropic diffusion, scalar parameters should be rotationally invariant (1,
2), rendering them independent of
the choice of the laboratory frame of reference or the position and orientation
of the subject or subject’s tissue being scanned.
What
are examples of quantitative metrics?
Common
quantitative metrics in diffusion MRI applications include the Trace of the
apparent diffusion tensor (1,
3) or one-third of it -- the mean
apparent diffusion coefficient (mADC), which characterize the isotropic part of
the tissue water diffusion process. Dimensionless scalar parameters that
measure the degree of diffusion anisotropy, like the Fractional Anisotropy (FA) (4,
5) (or the
Relative Anisotropy) represent ratios of different features of the anisotropic diffusion
process. It is critical that such scalar parameters, like the FA, are also designed to be rotationally
invariant so that they reveal intrinsic
features of the underlying physical process independent of the experimental
design, details and laboratory or subject coordinate system.
What
is the value of quantitative metrics?
From
a bench, preclinical, and clinical perspective, our ability to more sensitively
and selectively extract useful microstructural features, independent of
experimental design, hardware used, … should enable us to better follow changes
in normal and abnormal development, disease progression, tissue degeneration,
sequelae of brain trauma and effects of normal and pathological aging
processes. In an era of precision medicine, they should improve our ability to
make genotype-phenotype correlations, as in the ENIGMA DTI project (6).
Challenges and requirements of measuring and
mapping quantitative imaging metrics: Generally,
“It takes a pipeline” to measure and map a quantitative imaging metric.
Quantitative imaging is significantly more challenging than weighted imaging,
and much more costly and labor intensive. One must have a robust experimental
design to enable the measurement robustly, accurately and precisely. There must
be a means to calibrate the scanner for a measurement to assess artifacts,
image quality, etc. There needs to be an accounting of all known or foreseeable
experimental, biological, and imaging artifacts and a means to correct or at
least mitigate their effects. There should be a quantitative assessment of the
effect of the image reconstruction process on the generation of the MR signal
of interest. There has to be a physico-mathematical framework in place to
estimate quantities or parameters from the imaging data. There has to be framework
for mapping and analyzing the data using a rigorous statistical approach, such
as hypothesis testing. There should be methods in place for segmenting and
clustering imaging, and a means to do meaningful group (population) analysis.
What
are candidate higher-order or “advanced” quantitative imaging parameters or
metrics?
In
the field of diffusion MRI, the advent of higher-order tensor (HOT) models of diffusion,
such as Generalized Diffusion Tensor Imaging (7,
8) Diffusion
Kurtosis Imaging (7,
9) and Mean
Apparent Propagator (MAP)-MRI (10) have enlarged
the family of potential quantitative imaging quantities and parameters to
characterize features of diffusion processes occurring in tissue that do not necessarily
obey the classical Einstein equation. These models can be viewed as
“non-parametric” because they do not explicitly contain material or microstructural
parameters to estimate or measure. MAP-MRI is an implementation of Callaghan et
al.’s “k-space and q-space imaging” (11) that requires a vastly reduced
amount of DWI data to construct the “average propagator” or MAP. Several new candidate quantitative metrics can
be derived from extracting features of the size, shape, and orientation of these
MAPs in each voxel. One such feature is the Non-Gaussianity (NG). This represents the fraction of the
MAP that cannot be fit by a DTI or Gaussian diffusion model. Another is the
Propagator Anisotropy (PA). This
quantity reports the fraction of the propagator that cannot be described by an
isotropic displacement distribution. In some ways, this is a generalization of
the FA in DTI. Another quantitative parameter of interest is
the return-to-origin probability (RTOP).
In a restricted diffusion process, RTOP
can be shown (in the limit of long diffusion times) to be inversely related to
the volume of a confining pore. RTOP’s
dependence on diffusion time may make it less robust as an imaging biomarker
than other MAP-MRI derived parameters
What
candidate higher-order or “advanced” imaging parameters or metrics arise from
parametric modeling approaches?
There have been several parametric models of diffusion
within a voxel presented in the past decade. Parametric models, such as NODDI (12), CHARMED (13), ActiveAx (14), and AxCaliber (15), all ascribe the total MR signal
to a sum of contributions from different distinct, non-exchanging compartments,
with different assumed morphologies. The features that these models yield, such
as the intra-axonal fraction, neurite density, extracellular fraction, … are
potentially useful in characterizing tissue microstructure and their possible
alterations in pathology.
Overhanging issues with proposed quantitative parameters
from “advanced” imaging methods?
For non-parametric models like DKI, one must ask about the
possible robustness of kurtosis tensor-derived parameters, including the mean
Kurtosis. These parameters are measured from DWIs using a model that is
quadratic as opposed to a linear with respect to “b”. Diffusion Kurtosis
Imaging (DKI) parameters are therefore sensitive to the range of b-values used
and will invariably will change if higher order terms are included in the diffusion
model. DKI parameters also inherently sensitive to the timing parameters of the
diffusion gradients (i.e., the pulse duration and diffusion time), used to
acquire the DWIs.
For parametric modeling approaches, one must inquire about
whether the proposed assumed microstructure is consistent with the ground biological
truth. More testing and vetting are required to assess these parametric models before
parameters or features derived from them are promulgated from bench to bedside
as quantitative metrics or biomarkers.
Other candidate quantitative parameters:
Another class of potential quantitative parameters could be
based on the hierarchical character of many biological tissues. Structural
tissues, such as tendon, ligaments, and skeletal muscle are characterized by a “bundles-within-bundles”
fractal-like architecture. This architectural paradigm is also present in white
matter pathways in the CNS and nerve fascicles in the PNS. Extracellular matrix
and brain parenchyma also have a hierarchical structure roughly consisting of tubes,
sticks, beads, and plates organized at different hierarchical length scales (16).
Despite the fact that we cannot directly probe sub-micron displacements
in tissues with current MR hardware, we can infer some of these tissue’s
topological features at finer length scales. By measuring the average
propagator over a range of diffusion times, and using models of fractal diffusion
appropriate for each tissue type, we can extrapolate into a high-resolution
regime, inferring architectural features of tissue too small to be probed
directly by conventional PFG MR (17). Although more data intensive
than conventional PFG MRI, these approaches may one day provide sensitive, early
biomarkers of changes in tissue microstructure associated with disease and
development. There are several newly proposed parameters that characterize the
fractal structure of complex, hierarchically organized media, which possess
many of the features of quantitative imaging parameters.
Acknowledgements
PJB is supported by funding from the NIH Intramural Research Program and the Eunice Kennedy Shriver National Institute of Child Health and Human Development.References
1. Basser
PJ & Le Bihan D (1992) Fiber orientation mapping in an anisotropic medium
with NMR diffusion spectroscopy. 11th
Annual Meeting of the SMRM, p 1221.
2. Basser
PJ, Mattiello J, & LeBihan D (1993) MR imaging of fiber-tract direction and
diffusion in anisotropic tissues. 12th Annual
Meeting of the Society of Magnetic Resonance in Medicine, p 288.
3. Basser
PJ, Mattiello J, & LeBihan D (1994) MR diffusion tensor spectroscopy and
imaging. Biophys J 66(1):259-267.
4. Basser
PJ (1995) Inferring microstructural features and the physiological state of
tissues from diffusion-weighted images. NMR
Biomed 8(7-8):333-344.
5. Pierpaoli
C & Basser PJ (1996) Toward a quantitative assessment of diffusion
anisotropy. Magn Reson Med
36(6):893-906.
6. Jahanshad
N, et al. (2013) Multi-site genetic
analysis of diffusion images and voxelwise heritability analysis: a pilot
project of the ENIGMA-DTI working group. Neuroimage
81:455-469.
7. Liu CL,
Bammer R, & Moseley ME (2003) Generalized diffusion tensor imaging (GDTI):
A method for characterizing and imaging diffusion anisotropy caused by
non-Gaussian diffusion. ISRAEL J CHEM
43(1-2):145-154.
8. Ozarslan
E & Mareci TH (2003) Generalized diffusion tensor imaging and analytical
relationships between diffusion tensor imaging and high angular resolution diffusion
imaging. Magn Reson Med
50(5):955-965.
9. Jensen
JH, Helpern JA, Ramani A, Lu H, & Kaczynski K (2005) Diffusional kurtosis
imaging: the quantification of non-gaussian water diffusion by means of
magnetic resonance imaging. Magn Reson
Med 53(6):1432-1440.
10. Ozarslan E, et al. (2013) Mean apparent propagator
(MAP) MRI: a novel diffusion imaging method for mapping tissue microstructure. Neuroimage 78:16-32.
11. Callaghan
PT, Eccles CD, & Xia Y (1988) NMR microscopy of dynamic displacements:
k-space and q-space imaging. J. Phys. E:
Sci. Instrum. 21:820-822.
12. Zhang H,
Schneider T, Wheeler-Kingshott CA, & Alexander DC (2012) NODDI: practical
in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage 61(4):1000-1016.
13. Assaf Y
& Basser PJ (2005) Composite hindered and restricted model of diffusion
(CHARMED) MR imaging of the human brain. Neuroimage
27(1):48-58.
14. Dyrby TB,
Sogaard LV, Hall MG, Ptito M, & Alexander DC (2013) Contrast and stability
of the axon diameter index from microstructure imaging with diffusion MRI. Magn Reson Med 70(3):711-721.
15. Assaf Y,
Blumenfeld-Katzir T, Yovel Y, & Basser PJ (2008) AxCaliber: a method for
measuring axon diameter distribution from diffusion MRI. Magn Reson Med 59(6):1347-1354.
16. Kinney JP, et al. (2013) Extracellular sheets and
tunnels modulate glutamate diffusion in hippocampal neuropil. J Comp Neurol 521(2):448-464.
17. Ozarslan
E, Shepherd TM, Koay CG, Blackband SJ, & Basser PJ (2012) Temporal scaling
characteristics of diffusion as a new MRI contrast: findings in rat
hippocampus. Neuroimage
60(2):1380-1393.