Functional Quantitative Susceptibility Mapping (fQSM)...olunteers were sc ed, which alterna tapping...

1
1 H Notting Introd Functio suscept more re depend gradien signific functio induced Materi Two v present finger- EPI da followi the slic space d modulu was the in the t of itera ppm. (v iterativ betwee derivat Result Tempo subseq variatio for suc relative perform maps a shows that ac |Δφ|-m negativ similar a result highest phase c Conclu Quanti extende activati change possibi a basi suscept pulse-s algorith if phas there i modific Refere MRM Neuroi (2011) Panchu 103:25 47:344 in pres 46:662 JCBFM (2005) (2009) David Z Ba High-Field MR Cen gham, Nottingham duction onal activation ma tibility) and imagi eliable at the same dent. Its measurem nt effects. Quantit cant development onal QSM (fQSM) d susceptibility-ch ials and Methods volunteers were sc ted, which alterna tapping with the l ata (TE=25ms, TR= ing order for optim ce direction were c data as a nuisance us data were motio en performed with time-series [8], 8 th ations set to 100 [ viii) These maps s ve deconvolution a en 1 and 30 TRs) tive, convolved wi ts oral stability was i quent spatial filteri ons and their appli ccessful quantificat e functional suscep med and the dipole and maps of the p the estimated |Δφ ctivation-induced B maps additionally s ve susceptibility S rity, in spite of diff t of the BOLD-ch t BOLD changes changes in the pare usion tative susceptibi ed for use in fun ion-induced BO es were mapped a ility of BOLD-fM ic physiological tibility), limiting sequence parame hm can be applie se information has s no need for har cations. ences: 1. de Roc 63:194; 2.Wharton image 53:515; 3 Neuroimage 54 uelo et al., (20 544; 5. Pfeuffer e 4; 6. Hagberg et al. s; 7. Hussein et al. 23; 8.Tomasi an M 27:33; 9. Ma Concepts MR 25 NeuroImage 47:5 alla 1,2 , Rosa M San ntre, Max Planck In m, United Kingdom aps based on BOL ing (e.g. TE) param e field strength, fa ment requires the a tative susceptibili in recent years, ) as a quantitative hanges. s canned on a 7T P ated 12s blocks of left hand, somatos =3s) with 1mm iso mal temporal stabil corrected by rewin e regressor [6]. (iv on-corrected using h one of the follow order polynomial 1,2]. (vii) The resu served as inputs dir algorithm [2] with in the time doma ith a double-gamm improved significa ing step. After init ication did not intr tion of baseline su ptibility maps cou e-fit reconstruction phase change estim |-map overlaid on BOLD susceptibil served for optimiz SPMs overlaid on ferences between T hange in blood sus in modulus, phase enchyma (1 o predi ility mapping nctional imaging OLD susceptibi at 7T. This offers MRI related directly quantity (magn the dependence eters. The propo d on any time-ser s been collected, rdware or acquisi chefort et al., (20 n and Bowtell, (20 3. Schweser et 4:2789; 4. Sanch 10) J Neurophy et al., (2002) M ., (2011) Neuroima ., (2007) Appl. Op d Caparelli, (20 arques and Bow B:65; 10. Feng et 540 Functional nchez-Panchuelo 2 , nstitute for Biolog m, 3 Physiology of C Neuroimaging a LD-fMRI provide meters. As an atte st T2*-mapping ca acquisition of two ty mapping (QSM partly due to the e approach for BO Philips scanner du task/stimulation w sensory stimulus: 6 otropic resolution lity and statistical nding the phase in v) Spatial unwrap g 3dAllineate (AFN wing methods: hom background field ulting phase-imag rectly for the calcu a maximum of 15 ain before uncorre ma HRF and filtere antly by correcting tial global filtering roduce temporal m usceptibility maps. uld be produced w n method was use mated using forwa the filtered phase ity changes are no zation of the recon n the negative QS T2* (Fig.1) and ne ceptibility, but thi e and susceptibilit icted by [10]) and was and ility the y to netic on osed ries, and tion 010) 010) al., hez- ysiol RM age, ptics 007) wtell t al., l Quantitative Sam Wharton 2 , Gi gical Cybernetics, Cognitive Processe and MR-Physics, U qualitative topolog mpt to relate BOL an serve as an alte o echoes per repeti M) uses phase inf emergence of ad OLD-fMRI, and p uring a multi-task with 18s blocks of 60 Hz vibrotactile [4] spanning eithe power. (i) Phase k-space [5]. (iii) G pping was perform NI), and the resulti modyne Gaussian w removal, SHARP ges, except for the ulation of phase-SP 5 iterations. Filtere cted activation ma d in time with the g for global intens g in the time doma modulation. Remo The relative phas with this method. F ed. Filtered phase i ard field calculatio image. The measu ot exclusively enc nstruction process SM. Comparison egative susceptibil s negative change ty were observed to quantify suscep e Susceptibility isela E Hagberg 3 , K Tuebingen, Germa es, Max Planck Ins University of Tuebi gies of task-relate LD changes to an ernative to T2*-we ition, is sensitive formation from M dvanced backgrou proposes a multi-s k fMRI experimen f rest; visual stimu e stimulation appli er the visual or the maps were maske Global temporal ph med using a fast 3 ing correction mat with kernel size 12 P with a threshold o relative phase filt PMs or for the QS ed modulus and ph aps (p<10 -10 ) were same parameters a sity variations. Th ain, all of the app ving slow phase m se filtering only ret Figures 1-5 present images revealed fi on [9], by replacin ured |Δφ|-map is p coded in the T2*- s. Fig.4 presents a of activation patt lity contrast (Fig.4 e coexists with pos around veins. Filt ptibility changes w y Mapping (fQ Klaus Scheffler 1,4 , any, 2 Sir Peter Ma stitute for Biologic ingen, Tuebingen, ed modulus intensi absolute quantity eighted imaging. H to unwanted dyna MR images to cal und field removal step filtering algor nt. Paradigm: A v ulus: a flickering ied to the fingertip e sensorimotor cor ed using binarised hase variations we 3D algorithm [7], trices were then ap 2mm, relative pha of 0.1 and a kerne tered data-sets, we SM calculations. Su hase images and su e calculated with as the functional d his was essential f lied spatial filters modulations of ext etained information ted below result fr fine structural brai ng susceptibility w presented in Fig.3 -weighted image i an SPM derived f terns in modulus 4). As expected, su sitive modulation tering out this vas with better spatial c QSM) , Sue T Francis 2 , a ansfield Magnetic R cal Cybernetics, Tu Germany ity changes which and to render the However, T2* is an amic effects and t lculate local tissu methods [1-3]. T rithm for the calc visual, motor or s checkerboard ring ps of the left hand rtex. Analysis: Pro modulus images. ere corrected by u and was followed pplied to the phase ase filtering by com el size of 2mm [3] ere then divided b usceptibility maps usceptibility maps FSL. The GLM i data. Pre-whitening for efficient reduc were successful i ternal origin from n about dynamic c from an fMRI-data in contrast (Fig.2). with the modulus . Significant differ intensity. The com from the magnitud (Fig.1) and QSM usceptibility decrea in neighbouring v scular contribution correlation to neura and Richard W Bow Resonance Centre uebingen, German h depend on physio comparison of BO n MR-specific par to static, but non-l ue susceptibility. I This study introdu culation of maps s somatosensory blo g on a grey backgr d. Protocol: Zoom ocessing steps wer (ii) Global intensi using the phase at t d by temporal unw e data. (vi) High-pa mplex division wi or dipole-fit with by γB 0 TE×10 -6 to y s (QSM) were calc s were smoothed (b included the block g of the data was n ction of the backgr in removing low-f the region of inte changes in phase a aset where a finge . Activation induc %-change, were i rences between Fi mparison of estima de data along wit M data (Fig.4-5) r ases (-Δχ) up to 7% voxels (Fig.5). Als n will be necessar al activity. wtell 2 e, University of ny, 4 Department of ological (magnetic OLD-fMRI studies rameter and is field local susceptibility It has experienced uces time-resolved showing activation ock-paradigm wa round, motor task ed multi-slice GE re performed in the ty variations along the origin of the k wrapping. (v) The ass spatial filtering ith the first volume maximum numbe yield field maps in culated by using an band-pass filtering k-paradigm and it not performed. round fields in the frequency intensity erest was importan and, therefore, only er tapping task wa ced Δφ (%-change investigated. Fig.2 igs.2 and 3 sugges ated and measured th the positive and reveals remarkable % on activation a o, as expected, the ry for detection o f c s d y d d n s : - e g - e g e r n n g s e y t y s ) 2 t d d e s e f 325 Proc. Intl. Soc. Mag. Reson. Med. 20 (2012)

Transcript of Functional Quantitative Susceptibility Mapping (fQSM)...olunteers were sc ed, which alterna tapping...

Page 1: Functional Quantitative Susceptibility Mapping (fQSM)...olunteers were sc ed, which alterna tapping with the l ta (TE=25ms, TR= ng order for optim e direction were c ata as a nuisance

1HNotting

IntrodFunctiosusceptmore redependgradiensignificfunctioinduced

MateriTwo vpresentfinger-EPI dafollowithe slicspace dmoduluwas thein the tof iterappm. (viterativbetweederivat

ResultTemposubseqvariatiofor sucrelativeperformmaps ashows that ac|Δφ|-mnegativsimilara resulthighestphase c

ConcluQuantiextendeactivatichangepossibia basisusceptpulse-salgorithif phasthere imodific

RefereMRM Neuroi(2011) Panchu103:2547:344in pres46:662JCBFM(2005) (2009)

David Z BaHigh-Field MR Cen

gham, Nottingham

duction onal activation matibility) and imagieliable at the same

dent. Its measuremnt effects. Quantitcant development onal QSM (fQSM)d susceptibility-chials and Methods

volunteers were scted, which alternatapping with the l

ata (TE=25ms, TR=ing order for optimce direction were cdata as a nuisanceus data were motioen performed withtime-series [8], 8th ations set to 100 [viii) These maps s

ve deconvolution aen 1 and 30 TRs) tive, convolved wits oral stability was iquent spatial filterions and their appliccessful quantificate functional suscepmed and the dipoleand maps of the pthe estimated |Δφ

ctivation-induced Bmaps additionally sve susceptibility Srity, in spite of difft of the BOLD-cht BOLD changes changes in the pareusion tative susceptibied for use in funion-induced BO

es were mapped aility of BOLD-fMic physiological tibility), limiting sequence paramehm can be appliese information hass no need for harcations.

ences: 1. de Roc63:194; 2.Whartonimage 53:515; 3

Neuroimage 54uelo et al., (20544; 5. Pfeuffer e4; 6. Hagberg et al.s; 7. Hussein et al.

23; 8.Tomasi anM 27:33; 9. Ma

Concepts MR 25NeuroImage 47:5

alla1,2, Rosa M Sanntre, Max Planck Inm, United Kingdom

aps based on BOLing (e.g. TE) parame field strength, fa

ment requires the atative susceptibiliin recent years,

) as a quantitativehanges. s canned on a 7T Pated 12s blocks of left hand, somatos=3s) with 1mm iso

mal temporal stabilcorrected by rewine regressor [6]. (ivon-corrected usingh one of the follow

order polynomial1,2]. (vii) The resu

served as inputs diralgorithm [2] with in the time doma

ith a double-gamm

improved significaing step. After initication did not intrtion of baseline suptibility maps coue-fit reconstruction

phase change estim|-map overlaid on BOLD susceptibilserved for optimizSPMs overlaid onferences between T

hange in blood susin modulus, phaseenchyma (1o predi

ility mapping nctional imaging OLD susceptibiat 7T. This offers

MRI related directlyquantity (magnthe dependence

eters. The propod on any time-sers been collected, rdware or acquisi

chefort et al., (20n and Bowtell, (203. Schweser et 4:2789; 4. Sanch10) J Neurophyet al., (2002) M., (2011) Neuroima., (2007) Appl. Opd Caparelli, (20arques and BowB:65; 10. Feng et

540

Functionalnchez-Panchuelo2, Institute for Biologm, 3Physiology of C

Neuroimaging a

LD-fMRI provide meters. As an attest T2*-mapping caacquisition of twoty mapping (QSMpartly due to the

e approach for BO

Philips scanner duf task/stimulation wsensory stimulus: 6otropic resolution lity and statistical

nding the phase in v) Spatial unwrap

g 3dAllineate (AFNwing methods: hom background field ulting phase-imagrectly for the calcua maximum of 15

ain before uncorrema HRF and filtere

antly by correctingtial global filteringroduce temporal m

usceptibility maps.uld be produced wn method was usemated using forwathe filtered phase ity changes are no

zation of the reconn the negative QST2* (Fig.1) and neceptibility, but thie and susceptibiliticted by [10]) and

was and ility the

y to netic

on osed ries, and tion

010) 010)

al., hez-ysiol RM age,

ptics 007) wtell t al.,

l QuantitativeSam Wharton2, Gi

gical Cybernetics, Cognitive Processeand MR-Physics, U

qualitative topologmpt to relate BOLan serve as an alte

o echoes per repetiM) uses phase inf

emergence of adOLD-fMRI, and p

uring a multi-taskwith 18s blocks of60 Hz vibrotactile[4] spanning eithepower. (i) Phase k-space [5]. (iii) G

pping was performNI), and the resultimodyne Gaussian w

removal, SHARPges, except for the ulation of phase-SP5 iterations. Filterected activation mad in time with the

g for global intensg in the time domamodulation. Remo The relative phas

with this method. Fed. Filtered phase iard field calculatio

image. The measuot exclusively encnstruction process

SM. Comparison egative susceptibils negative changety were observed to quantify suscep

e Susceptibilityisela E Hagberg3, KTuebingen, Germaes, Max Planck InsUniversity of Tuebi

gies of task-relateLD changes to an ernative to T2*-weition, is sensitive formation from Mdvanced backgrouproposes a multi-s

k fMRI experimenf rest; visual stimu

e stimulation applier the visual or themaps were maske

Global temporal phmed using a fast 3ing correction matwith kernel size 12

P with a threshold orelative phase filt

PMs or for the QSed modulus and phaps (p<10-10) weresame parameters a

sity variations. Thain, all of the appving slow phase m

se filtering only retFigures 1-5 presentimages revealed fion [9], by replacinured |Δφ|-map is pcoded in the T2*-s. Fig.4 presents aof activation pattlity contrast (Fig.4e coexists with pos

around veins. Filtptibility changes w

y Mapping (fQKlaus Scheffler1,4,any, 2Sir Peter Mastitute for Biologicingen, Tuebingen,

ed modulus intensiabsolute quantity

eighted imaging. Hto unwanted dyna

MR images to calund field removal step filtering algor

nt. Paradigm: A vulus: a flickering ied to the fingertipe sensorimotor cored using binarised hase variations we

3D algorithm [7], trices were then ap2mm, relative phaof 0.1 and a kernetered data-sets, we

SM calculations. Suhase images and sue calculated with as the functional d

his was essential flied spatial filters

modulations of extetained informationted below result fr

fine structural braing susceptibility wpresented in Fig.3-weighted image ian SPM derived fterns in modulus 4). As expected, susitive modulation tering out this vas

with better spatial c

QSM) , Sue T Francis2, a

ansfield Magnetic Rcal Cybernetics, TuGermany

ity changes whichand to render the

However, T2* is anamic effects and tlculate local tissu

methods [1-3]. Trithm for the calc

visual, motor or scheckerboard ring

ps of the left handrtex. Analysis: Promodulus images. ere corrected by uand was followed

pplied to the phasease filtering by comel size of 2mm [3] ere then divided busceptibility mapsusceptibility mapsFSL. The GLM i

data. Pre-whitening

for efficient reducwere successful i

ternal origin from n about dynamic cfrom an fMRI-datain contrast (Fig.2).with the modulus . Significant differintensity. The comfrom the magnitud(Fig.1) and QSM

usceptibility decreain neighbouring vscular contributioncorrelation to neura

and Richard W BowResonance Centre

Tuebingen, German

h depend on physiocomparison of BO

n MR-specific parto static, but non-lue susceptibility. IThis study introduculation of maps s

somatosensory blog on a grey backgrd. Protocol: Zoomocessing steps wer(ii) Global intensi

using the phase at td by temporal unw

e data. (vi) High-pamplex division wior dipole-fit with

by γB0TE×10-6 to ys (QSM) were calcs were smoothed (bincluded the blockg of the data was n

ction of the backgrin removing low-fthe region of inte

changes in phase aaset where a finge. Activation induc%-change, were i

rences between Fimparison of estimade data along wit

M data (Fig.4-5) rases (-Δχ) up to 7%

voxels (Fig.5). Alsn will be necessaral activity.

wtell2 e, University of ny, 4Department of

ological (magneticOLD-fMRI studiesrameter and is fieldlocal susceptibilityIt has experienceduces time-resolvedshowing activation

ock-paradigm waround, motor tasked multi-slice GE

re performed in thety variations alongthe origin of the kwrapping. (v) Theass spatial filteringith the first volumemaximum numbe

yield field maps inculated by using anband-pass filteringk-paradigm and itnot performed.

round fields in thefrequency intensityerest was importanand, therefore, onlyer tapping task waced Δφ (%-changeinvestigated. Fig.2igs.2 and 3 suggesated and measuredth the positive andreveals remarkable% on activation ao, as expected, thery for detection o

f

c s d y d d n

s k: -e g -e g e r n n g s

e y t y s ) 2 t d d e s e f

325Proc. Intl. Soc. Mag. Reson. Med. 20 (2012)