resting-state fmri analysis

Hello there:
i am tring to do the resting-state fmri analysis. i did the preprocessing step using afni_proc.py then i did the first level analysis to detect the DMN network but when i underlay the resulted anatomic dataset and overlay the resulted the correlation maps i get a correlatin outside the brain borders(The areas outlined in blue).
Here is the afni proc py used in the preprocessing:

#!/bin/tcsh

afni_proc.py
-subj_id s01
-blocks despike tshift align tlrc volreg blur mask
scale regress
-copy_anat anat_shft+orig
-dsets func+orig.HEAD
-tcat_remove_first_trs 2
-align_opts_aea -cost lpc+ZZ -giant_move
-tlrc_base TT_N27+tlrc
-tlrc_NL_warp
-volreg_align_to MIN_OUTLIER
-volreg_align_e2a
-volreg_tlrc_warp
-volreg_warp_dxyz 2.5
-blur_size 4
-mask_segment_anat yes
-mask_segment_erode yes
-mask_import Tvent template_ventricle_2.5mm+tlrc
-mask_intersect Svent CSFe Tvent
-mask_epi_anat yes
-regress_motion_per_run
-regress_ROI_PC Svent 3
-regress_ROI_PC_per_run Svent
-regress_make_corr_vols WMe Svent
-regress_anaticor_fast
-regress_censor_motion 0.2
-regress_censor_outliers 0.05
-regress_apply_mot_types demean deriv
-regress_est_blur_epits
-regress_est_blur_errts
-regress_run_clustsim yes

Hi, Osman-

That is a good thing to check. There is even an automatically generated image of seedbased correlation for the DMN in the afni_proc.py generated QC HTML. From your output/results directory, copy+paste:


afni_open -aw QC_*/index.html

… and jump to the “vstat” section. How does that look?

The BOLD effect can be observed outside the brain (meaning here, the final anatomical volume or template brain used for alignment) for various reasons:

  • B0 inhomogeneity geometrically stretches/compresses/distorts the EPI volume, so it doesn’t perfectly match the subject anatomical
    → check this in the ‘vorig’ section of the QC HTML
  • alignment can be imperfect between EPI-anatomical, and/or anatomical-template
    → check this in the ‘ve2a’ and ‘va2t’ sections of the QC HTML, mentioned above
  • the BOLD information measured with FMRI does not measure neuronal activity, but instead blood oxygenation affecting the magnetic field. Blood vessels around the brain can affect the magnetic field greatly (because there some big ones there), so one can see some effects occurring outside the brain at times
  • there can be noise—well, there is noise in the data
  • if the EPI data is blurred during processing (as in your case), some signal from inside the brain can get smeared outside a bit

That all being noted, the amount of correlation leaking/overlapping the brain there doesn’t seem too terrible to me. But it would be good for you to verify features of your data and processing, noted above.

–pt

And blood vessels can carry BOLD signal from inside the brain to outside of it. So any blood vessel departing area of activity will likely show activation along the vessel, including as it goes outside of the brain.

  • rick

thank you for you valuable information. i attached the “vstat” section. how can we comment this incident?

Hi, Osman-

To me, this looks like a reasonable single-subject seedbased-correlation pattern for a PCC-placed seed (default mode network), with standard resting state data.

The QC HTML there will show piece-by-piece useful aspects of processing to verify–e.g., the alignment images just preceding the vstat–so check everything make sure that all is well (or at least that there aren’t glaring issues…).

–pt