Chapter #1: Functional Connectivity Demonstration¶
Before we begin using the CONN toolbox, a short demonstration of functional connectivity using a simpler, easy-to-understand method may help prepare you for using a more sophisticated package. I prefer to use the fMRI analysis package FSL for this kind of demonstration, although both AFNI and SPM are able to do functional connectivity analyses as well. Before going on, I recommend that you work through the FSL tutorials to familiarize yourself with how FSL works, what a time-series is, and how to extract data from regions of interest (ROIs). If you already know how to do a functional connectivity analysis in FSL and simply want to get started with the CONN toolbox, click the
Next button down below.
This tutorial will use a sample dataset that you can download here. The file size is about 20MB.
Functional Connectivity Analysis with FSL¶
As you saw in the FSL tutorial, each voxel of an fMRI dataset contains a time-series, or sequence of BOLD activation measurements that are concatenated together. The sampling rate of the time-series is the Time to Repetition (TR), which for most fMRI studies is around 2-3 seconds. For each voxel, therefore, the BOLD activation is measured at each TR, and this continues for a period of time specified by the user. This period of time when the scanner is measuring the BOLD signal is called a run (or session, in SPM), and a typical experiment is composed of several runs.
To examine the time-series of the sample dataset, navigate to where you downloaded the
TimeSeries_ExampleData.nii.gz file, and type:
This will open the file in the fsleyes viewer. From the menu at the top of your window, click on
View -> Time series. This will open a time-series panel at the bottom of the fsleyes viewer. To make the numbers easier to understand, and to enable us to compare time-series from distant voxels, click on the “Plotting mode” dropdown menu and select “Normalised”. This will scale the time-series relative to a mean of 0.
Now, click on a random voxel, and then from the menu at the top of the screen select
Tools -> Seed correlation (Pearson). The voxel you have currently selected will be a seed voxel; the time-series of this seed voxel will be used as a reference. This will generate a connectivity map between your seed voxel and every other voxel in the brain, with hotter colors (i.e., more intense shades of red and yellow) representing greater correlation. Cooler colors (i.e., shades of blue), by contrast, represent where there is negative correlation.
You will also notice that this has created a new file in the Overlay list,
TimeSeries_ExampleData/correlation. If you highlight this image by clicking on it, you can then set the “Min.” and “Max.” thresholds to only show those voxels that correlate with the seed region at or higher than a certain level. Set the “Min” field to
0.4, and take a look at which voxels remain after this threshold.
If you’ve read the Cluster Correction appendix, you will remember that one of the justifications for cluster correction is that a given voxel is not completely independent from its neighbors. How do you see that principle being reflected in this correlation map?
Now let’s see whether we can replicate one of the most well-known functional connectivity networks, the so-called Default Mode Network. This network is a pattern of correlated regions, primarily the ventromedial prefrontal cortex and the posterior cingulate cortex. Use the “Voxel location” fields to enter values of X=32, Y=53, and Z=16, and use this voxel as a seed for another functional connectivity analysis. What do you notice about the resulting correlation map?
Functional Connectivity Analysis with Neurosynth¶
Another useful tool for understanding functional connectivity is Neurosynth. Neurosynth is primarily a meta-analysis tool, but it can also be used to quickly generate functional connectivity maps for seed regions in a template brain. If you click on the “Locations” tab, a new window opens that says “Functional connectivity and coactivation maps”. Click on any coordinate in the brain, and then click on “What’s here?” After a few moments, a correlation map will be generated.
For example, let’s say that we wanted to generate a functional connectivity map for the default mode network (described above). If you click on the coordinates [0, 44, -14] and click the “What’s here?” button, you should see significant clusters in both the region that you clicked, and also in the posterior cingulate cortex. These maps are generated based on a sample of a thousand subjects collected by Yeo et al. (2011), and contain the typical networks you would see in any sample.
As an exercise, try setting the coordinates to one of the nodes of the salience network, for example - [0 26 32], which is roughly within the anterior cingulate cortex - and see whether it correlates with the bilateral insula. Try setting the seed in other regions, such as the motor cortex or auditory cortex. What other regions would you expect would correlate with those seed regions? Why?
Click here for a brief overview of the history of functional connectivity, and how to do a simple functional connectivity analysis in FSL.
Now that you have an idea of what functional connectivity is and how we generate correlation maps, you are ready to begin learning about the CONN toolbox. Click
Next to learn about how to first download the data we will be analyzing with CONN.