Diffusion Tensor Imaging I. Jennifer Campbell

Transcription

1 Diffusion Tensor Imaging I Jennifer Campbell

2 Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications

3 Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications

4 Diffusion MRI measures Brownian motion of water molecules Path of diffusing water molecule Water displacement distribution

5 Diffusion MRI measures Brownian motion of water molecules P(r r 0,! d) = 1 $ (r # r 0) T D #1 (r # r 0) ' exp 3 D(4"! d) % & 4! d ( ) ;< r 2 >= 6D! d

6 Diffusion MRI measures Brownian motion of water molecules Hindered diffusion: Apparent Diffusion Coefficient (ADC)

7 Diffusion MRI measures Brownian motion of water molecules Hindered diffusion: Apparent Diffusion Coefficient (ADC)

8 Diffusion MRI measures Brownian motion of water molecules P(r r 0,! d) = 1 % (r $ r 0) T (ADC) $1 (r $ r 0) ( exp 3 ADC"(4#! d) & ' 4! d ) *

9 Tissue structures determine which directions of motion are most probable White matter fibre bundle: oriented structure Water molecules prefer to travel parallel to fibre direction

10 Diffusion MRI Tractography

11 Tissue structures determine which directions of motion are most probable P(r r 0,! d) = 1 $ (r # r 0) T D #1 (r # r 0) ' exp 3 D (4"! d) % & 4! d ( )

12 Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications

13 The diffusion tensor! # # # " Dxx Dxy Dxz Dyx Dyy Dyz Dzx Dzy Dzz $ & & & % P(r r 0,! d) = 1 $ (r # r 0) T D #1 (r # r 0) ' exp 3 D (4"! d) % & 4! d ( )

14 The diffusion tensor "! % $ $ $ 0! ! 3 ' ' ' λ 1 λ 3 λ 2 e 1 # & P(r r 0,! d) = 1 $ (r # r 0) T D #1 (r # r 0) ' exp 3 D (4"! d) % & 4! d ( )

15 The diffusion tensor "! % $ $ 0! 2 0 ' ' λ 1 λ 3 λ 2 e 1 $ 0 0! 3 ' # &

16 The diffusion tensor "! % $ $ 0! 2 0 ' ' λ 1 λ 3 λ 2 e 1 $ 0 0! 3 ' # &

17 The diffusion tensor "! % $ $ 0! 2 0 ' ' λ 1 λ 3 λ 2 e 1 $ 0 0! 3 ' # &

18 Measuring the diffusion tensor: diffusion weighted images (DWIs)

19 Computing the diffusion tensor b value describes the extent of diffusion weighting in the MRI sequence S(b) = Soe!bD b matrix used in 3D tensor calculation S(b) S(b) = Soexp(! " " b ij D ij ) 3 i=1 3 j=1 b

20 Computing the diffusion tensor b value describes the extent of diffusion weighting in the MRI sequence S(b) = Soe!bD b matrix used in 3D tensor calculation log(s(b)) S(b) = Soexp(! " " b ij D ij ) 3 i=1 3 j=1 b

21 Computing the diffusion tensor λ 1 λ 3 λ 2 e 1

22 Maps obtainable from diffusion tensor FA = 3 2 FA: anisotropy index (!1 "!) 2 + (! 2 "!) 2 + (! 3 "!) 2 (!1 2 +! 3 2 +! 3 2 ) trace of diffusion tensor MD =! = trace / 3 =!1 +! 2 +! 3 3 RGB plot: principal eigenvector (e 1 ) direction, scaled by FA

23 Diffusion tensor imaging (DTI) λ 1 λ 2 λ 3 cylindrical symmetry

24 Diffusion tensor imaging (DTI) λ, axial λ, radial T cylindrical symmetry

25 Diffusion tensor imaging (DTI) normal white matter

26 Diffusion tensor imaging (DTI) myelin degeneration: radial diffusivity increases FA decreases

27 Diffusion tensor imaging (DTI) normal white matter

28 Diffusion tensor imaging (DTI) axonal fragmentation: axial diffusivity decreases FA decreases

29 Diffusion tensor imaging (DTI) normal white matter

30 Diffusion tensor imaging (DTI) decreased axon density: radial diffusivity increases FA decreases

31 Diffusion tensor imaging (DTI) crossing fibres: FA lower than for single direction

32 Diffusion tensor imaging (DTI) degeneration of one fibre population in crossing case: FA increases

33 Diffusion tensor imaging (DTI) curvature: FA lower than for single direction

34 Diffusion tensor imaging (DTI) splay: FA lower than for single direction

35 Diffusion tensor imaging (DTI) normal tissue

36 Diffusion tensor imaging (DTI) cellular swelling: MD decreases

37 Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications

38 Diffusion weighted MR sequence TE S(b) = Soe!bD echo b = t! k(t) 2 dt ; k(t) = " G(t ')dt ' 0 t! 0 G G b =! 2 G 2 " 2 (# $ " / 3) δ Δ Δ Stejskal-Tanner diffusion weighted sequence: pulsed-gradient spin-echo (PGSE)

39 Gradients, phase, and motion 180 echo G G B 0 + G B 0, B180 o 0 B 0 + G B 0

40 Gradients, phase, and motion 180 echo G G B 0 + G B 0, 180 o B 0 + G B 0 B 0

41 Acquisition parameters for diffusion weighted images (DWIs) series of DWIs: encoding directions full brain coverage: minutes per 100 directions b value: s/mm 2 ; b~1000 s/mm 2 best for tensor voxel size: roughly 2mm x 2mm x 2mm isotropic voxels typically used adjust for population being studied, e.g. infants, patients

42 Sequence modifications: isotropic diffusion weighting RF Gx G G Gy G G Gz G G

43 Sequence design and artifacts diffusion MRI acquisitions generally use fast imaging techniques, e.g., single-shot EPI: reduce sensitivity to motion reduce acquisition time these sequences are prone to artifacts such as magnetic susceptibility induced distortion these artifacts confound registration with standard structural scans and make investigation of diffusion parameters difficult

44 Sequence design and artifacts to improve spatial resolution and reduce distortion, multi-shot imaging possible, but must use navigator echo to avoid phase errors can use other types of sequences, e.g., STEAM, SSFP,

45 Sequence design and artifacts Heidemann et al single-shot EPI readout segmented EPI readout

46 Eddy-current induced distortion eddy currents are induced by the rapidly changing magnetic field gradients required for diffusion sensitization these eddy currents induce magnetic fields which distort the image. artifacts: shifting, sheering, stretching, compression, signal loss

47 Eddy-current induced distortion without correction/ compensation, can have image shifts up to, e.g., 5mm affine registration can correct for some, but not all, of this artifact problem largely alleviated with twice refocused balanced echo (TRBE) sequence

48 Eddy-current induced distortion TE TE echo echo G G G G G G δ Δ δ Δ Δ Twice refocused balanced echo (TRBE)

49 Signal to noise ratio (SNR); spatial resolution (b=1000 s/mm 2 ) (b=3000 s/mm 2 )

50 Signal to noise ratio (SNR); spatial resolution near noise floor, FA and trace are decreased. however, note that in regions of isotropic diffusion, noise increases FA tradeoff between SNR, spatial resolution, achievable b value Jones and Basser 2004.

51 Signal to noise ratio (SNR); spatial resolution Higher field strength (7T) Heidemann et al

52 Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications

53 Diffusion weighted signal profile Diffusion weighted signal intensity: high value in directions perpendicular to fibres (lower diffusivity)

54 Probability of water displacement orientation distribution function (ODF) Diffusion tensor model

55 High Angular Resolution Diffusion Imaging (HARDI) q-ball imaging

56 High Angular Resolution Diffusion Imaging (HARDI) QBI DTI

57 Why use high angular resolution diffusion MRI? Reduce false positive and false negative tractography results due to crossings In voxels where fibres cross, single fibre approaches can yield: ambiguous fibre direction incorrect fibre direction only one fibre direction (that of fascicle with largest volume fraction) Examples of such regions: basis pontis, subcortical white matter, superior longitudinal fasciculus, acoustic radiations, projection from subgenual white matter to amygdala, optic chiasm, caudate nucleus, corpus callosum, cortical spinal tract, cingulate bundle.

58 Crossing fibre reconstruction approaches multi-tensor approaches (Alexander et al., Parker et al., others) Mixture models (Gaussian (Tuch et al.); Wishart (Jian et al.)) ball and multi-stick (Behrens et al.) Composite hindered and restricted model of diffusion (CHARMED) (Assaf et al.) diffusion spectrum imaging (DSI) (Wedeen et al.) q-ball imaging (QBI) (Tuch et al.) spherical deconvolution (Tournier et al., Anderson, others) other variants

59 Multi-tensor model Tuch et al Hosey et al. 2008

60 Behrens ball and stick model Behrens et al. 2007

61 q-space analysis of diffusion MRI data the diffusion displacement distribution is the Fourier transform of the diffusion signal as a function of q, where q =! G"

62 Diffusion Spectrum Imaging (DSI) 3D diffusion pdf 2D diffusion Orientation Distribution Function (ODF) Wedeen et al. ; Canales-Rodrigue et al.

63 q-ball imaging direct calculation of the 2D diffusion ODF calculation uses Funk- Radon transform Tuch et al. MRM 2004

64 Crossing fibre detection: QBI vs. Deconvolution Diffusion ODF Deconvolved ODF

65 Composite hindered and restricted model of diffusion (CHARMED) hindered restricted

66 Composite hindered and restricted model of diffusion (CHARMED) restricted pool fraction F fractional anisotropy FA

67 Inferring higher angular resolution from limited datasets DTI QBI regularized DTI Curve Inference using neighbourhood information

68 Beyond crossing: other complex subvoxel geometries Curve Inference to distinguish fanning from bending fibres Savadjiev et al. NeuroImage subvoxel fanning of fibres Campbell et al. ISMRM 2011.

69 Effect of complex subvoxel geometry on simple scalar measures FA lower when there are multiple subvoxel orientations trace lower when there are multiple subvoxel orientations (Vos et al. ISMRM 2011) The way multiple populations share one voxel affects the signal profile, and hence, mean diffusivity and kurtosis Peled et al. ISMRM 2011.

70 Fibre directions and uncertainty

71 Fibre directions and uncertainty bootstrap analysis using tensors calculated from subsets of data cones of uncertainty subtending 95% confidence intervals Jones et al.

72 Fibre directions and uncertainty

73 Residual bootstrap deconvolution 1. fit diffusion weighted signal profile to spherical harmonic basis 2. perform N iterations of residual bootstrap 2 generation of synthetic signal profiles 3.calculate fibre orientation distribution function using spherical deconvolution for all iterations. 1 Savadjiev et al., Jones et al., 2003; Berman et al., 2008; Jeurisson et al., 2010; Haroon et al., 2007

74 Assignment of confidence that a fibre exists in a given direction 4. calculate mean orientation and uncertainty for the fibres single fibre double crossing triple crossing Campbell et al. ISMRM 2011

75 Fibre directions and uncertainty Behrens et al Bayesian inference on marginal distribution of diffusion parameters performed in MCMC framework Hosey et al. 2008

76 Diffusion Imaging Molecular diffusion The diffusion tensor Diffusion weighting in MRI Alternatives to the tensor Overview of applications

77 Example: Stroke Moseley (1990): ADC decreases in ischemia. Causes: cell membrane permeability decreases lead to cellular swelling Early diagnosis: diffusion MRI can identify regions of ischemia within 30 minutes of arterial occlusion, while T 2 images show changes only after 2 hours Chronic stroke can also cause Wallerian degeneration, leading to reduced FA

78 Example: Multiple Sclerosis T 2 Image FA Image Reduced FA values suggest a possible breakdown in the structural integrity of axons passing through there

79 Example: Multiple Sclerosis FA decrease in NAWM: significant correlation with NA/Cr value (a measure of axonal integrity): r 2 =0.53, p=0.02 insignificant correlation with lateralization index in fmri finger tapping task (a measure of cerebro-functional compensation): r 2 =0.35, p=0.07 insignificant correlation with EDSS ( a measure of clinical disability): r 2 =0.38, p=0.06 Caramanos et al, ISMRM 2002

80 Tract-based spatial statistics (TBSS) project highest FA value in region ( centre of tract ) onto average FA skeleton reduces problems associated with, e.g., partial volume averaging, misalignment Interpretation may be difficult in the presence of pathology Smith et al., NIMG 2006

81 See Diffusion Tensor Imaging II for more applications