Documentation‎ > ‎

Group connectometry analysis

Introduction


What is connectometry?
Diffusion MRI connectometry uses permutation test to find out the association of white matter pathways with any study factor [1]. It can be combined with a simple linear regression that includes all relevant variables (e.g. sex, site difference) in the model. The results will show tracks that exhibit positive or negative correlation with a variable of interest.

Why is connectometry a powerful test?
Most studies used either track-based or region-based analysis to compare diffusion data along the tracks or within a given region. These approaches require users to assign pre-defined tracks or regions to sample diffusion indices for further analysis. The assigned regions or tracks, however, will inevitably include regions that have low SNR (near gray matter) or show no correlation at all (irrelevant branches), consequently bringing high variance and making results insignificant.

Connectometry adopts an entirely different mindset to avoid including irrelevant regions. The paradigm behind it can be summarized as "tracking the difference". This is fundamentally different from the conventional paradigm of "finding the difference in tracks or regions". The analysis first identifies voxels that have strong correlations and tracks along axonal fiber directions to identify the consecutive fiber segment that shows continuous positive and negative correlation. The tracking is only limit to regions with high correlation, thus minimizing the noise contributed from other non-correlated regions.

What is the biophysics behind connnectometry?

Connectometry compares spin distribution function (SDF) between subjects, a density-based measurement of diffusion at different orientations. SDF is different from diffusivity measurements such as FA, ADC, RD. Diffusivity measures how fast water diffuses, whereas SDF measures the density of diffusing water. A recent study [4] showed that SDF provides a unique structural characterization that can reliably identify individuals (termed local connectome fingerprint). Its reproducibility and uniqueness 

is higher than diffusivity-based measurements. It is noteworthy that since SDF reveals high individuality, the inter-subject variance of SDF can be very high (this does not mean that SDF is unreliable). As a result, connectometry is most suitable for longitudinal study, though cross-sectional study can also be benefited by connectometry analysis.

What can connectometry do in a brain study? 
Two types of connectometry analysis are available in DSI Studio:

One is group connectometry [pdf], which identifies tracks associated with group difference or correlated with a study variable. It can also be used in a longitudinal study to study the difference between paired data (e.g. pre- and post- treatment). Group connectometry use a linear regression model and can include other variables to consider in the regression.

Another is individual connectometry (see Individual Connectometry for instruction), which identifies damaged/enhanced pathways of an individual by comparing one subject with a group of normal subjects. 

You may refer to the publication at the bottom of this page to see how connectometry can be used in other studies.

The following documentation shows how to do group connectometry in DSI Studio:

STEP 1: Data Reconstruction


Create a connectometry database containing ALL subjects. At the end of this step, you will have a connectometry database, which is a file name with extension *.db.fib.gz

*You can add more subjects or remove subjects from the database later using the [Edit a Connectometry Database] function provided in the main window tab.

*Quality control: please check the R2 value in the connectometry database using the "source data" tab shown in the following figure. A substantially lower R2 indicates reconstruction error. You may contact me at frank.yeh@gmail.com about fixing the error.

STEP 2: Handle Longitudinal Scans (if applicable)

Skipping this if your study is a cross-sectional study without repeat scans of the same subjects.

If your study is a longitudinal study (e.g. comparing pre- and post-treatment of the same subject), then the scans of the same individuals needs to be matched to calculate their difference. The following are the steps:

a) Click "Open a Connectometry Database" and select the database created in STEP1:

b) In the top menu, select [Tools][Longitudinal scans...] to bring up the subject match window:

DSI Studio will try to match subjects automatically using "connectome fingerprint". 
The baseline (b) is the baseline scans, whereas the study scans (a) are the follow-up scans. 

Alternatively, you can load your own matching table using a text file. The text file should be the consecutive numbers of the matching pairs. The first number in a pair will be used as the baseline, and the second number will be the study scans. For example, a text file with "0 1 2 3 4 5" will match scan 0 (base) and scan 1(study), scan 2 and scan 3, scan 4 and scan 5. Please note that the first scan in the database is labeled by 0, not 1. 

The button "Match consecutive scans" will match scan 0 (base) and scan 1 (study), 2 and 3, 4 and 5, through out the entire database. 

The "metric" group box allows you to define how the difference is calculated. "a-b" calculates the absolute difference by "the study scan (a) - the baseline (b)". "(a-b)/(a+b)" calculates the difference as a ratio. The absolute difference is recommended.

Click "Ok" to calculate the difference between scans and go back to the database window. Use [File][Save DB as..]  to save the differences as a new connectometry db.db.connectometry db.

The new connectometry db will contain only the difference between the paired scans.

STEP 3: Group Connectometry


Open the db.fib.gz file using the "Group Connectometry Analysis" button in the main window (Under the connectometry tab) to bring up the following window. 
If the dialog popups a warning about poor image quality (low R2), you may need to check the data at the "Source Data" tab to find out which subject has a registration problem or motion artifact. The possible solutions are mentioned in the documentation for creating a connectometry database.

There are four analysis models available: multiple regression, group difference, paired group difference, individual difference. 

STEP 3.1: Load demographics

Prepare a text file that records any scalar values you would like to include in the regression model. The format is shown in the following figure. The first row is the name of the scalar values. The second row is the value for the first subject, and the rest follows. If there is any missing data, assign 9999 and check the "missing data labeled by" checkbox (as shown in the interface). DSI Studio will ignore subjects with missing data.


Click on the button labeled "open subjects demographics". Load the txt file that records all the scalar values that will be included in the regression model. In the demographic table, you need to check whether the values matches well with the subject IDs. The "X" button is for removing unwanted cases (e.g. case with missing data). There is no undo for this function, and you may need to reload the connectometry database if you mistakenly delete a case.

STEP 3.2 Regression model


Please choose variables to be included in the regression model. This can include any confounding factors like sex and handedness. There are several tips in choosing the variables in the model.

TIPS:
1.If you have multiples study variables, do you include all of them in the model unless you have enough sample size and really want to regress out the dependent effects between them. Otherwise, include basic variables like the gender and only one study variable at a time.
2. If your study variables are highly correlated to each other, consider using a PCA to get the principle components/

STEP 3.3 Setup parameters

Choose the study variable to run the connectometry analysis. DSI Studio will find the tracks correlated with this study variable. Assign the T-score threshold. Higher values give more specific results, whereas lower values give more sensitive results. You can run separate analyses by different threshold (e.g., t=1, 2, and 3) to map different levels of correlation. Each threshold is viewed as a different hypothesis and will receive an FDR. The FDR value gives us an idea of how reliable the findings are, and higher thresholds may not necessarily reduce FDR. My experience is that an FDR of 0.05 is highly confirmative while 0.05~0.2 suggests a high possibility of positive findings. Any findings with FDR lower than 0.5 may be worth reporting since the results have more "true positive" findings than "false positive" findings.


Advanced options (OPTIONAL)
Seeding density determines the total number of seeds used in the fiber tracking. A larger number takes more computation time but may give better track coverage.
Track trimming helps remove noisy findings. A Higher value may need higher seeding density. Track trimming helps to improve FDR.
Permutation count determines the total number of permutation applied to the FDR analysis. Higher values give smoother and more accurate estimation of FDR curves, but it requires more computation time. 
Length threshold: This threshold determines the minimum length of the affected tracks.
Normalize SDF: Recommend unchecked unless the directionality of the results is wrong. This SDF normalization will reduce the systemic error and difference, but it assumes that the maximum anisotropy is the same for all subjects to stabilize the quality. 
Tracking region: The default study region for connectometry is "whole brain", which means that connectometry will look at whole brain region. You can choose to use an ROI in the MNI space by loading a NIFTI file of the mask. You can also choose a region from the atlas. The region type can be ROI (select tracts passing the region), ROA (discard any track passing the region), End (select track ending in the region), Seed (Only start the seeds from the region). "Whole brain" is most suitable for the exploratory purpose. If you have specific region of interest, limiting the computation to an ROI can give more specific result (e.g. assign cerebral peduncle as the ROI can investigate whether CST affected). Using only one ROI is good for most cases unless you have specific need in the tracking routine.


STEP 3.4: Run connectometry

Click on the "STEP 4:Run" button and wait until the computation is finished.

TIPS:
1. If the direction of the correlation is wrong, try checking normalize SDF to fix the problem.
2. If you have very different FDR results at different permutation counts, increase the permutation count until the FDR converges. 
3. You can run analysis multiple times with different parameters. The result is still reliable since connectometry already considers multiple comparisons (as long as the permutation count is large enough)With a large permutation count (usually greater than 2000), it is almost impossible for users to manipulate the result by tuning the parameters. 


You may click on "View Results in 3D" button to view the tracks with a substantial difference. 

DSI Studio will output the following files:

demo.txt is the demography text file prepared in STEP0. 
demo.txt.length40.s1.p2000.bmi.t20.dist_value.txt records the distribution of track length.
demo.txt.length40.s1.p2000.bmi.t20.fdr_value.txt records the FDR with respect to length.
demo.txt.length40.s1.p2000.bmi.t20.greater.dist.bmp shows the histogram of track length. "greater" here means positive correlation in multiple regression, or group 0 greater than group 1 if group analysis was selected.
demo.txt.length40.s1.p2000.bmi.t20.greater.fdr.bmp shows the FDR with respect to track length.
demo.txt.length40.s1.p2000.bmi.t20.greater.fib.gz is a FIB file which can be opened in "STEP3 fiber tracking" to observe the how local connectome correlates with the study variable.
demo.txt.length40.s1.p2000.bmi.t20.greater.trk.gz is a track file of the connectometry finding.
demo.txt.length40.s1.p2000.bmi.t20.report.txt is a method report.

You may also check out the report to find out the FDR of the finding.

How to interpret and report the result?

First, you may visualize the tracks that are correlated with the study variable at different T threshold (see the next section about visualization). The FDR value will be different for different track length. You can either report the FDR values at 40 mm (default approach in the automatic report), or report which track length has FDR values lower than 0.05 (check out the raw data results). 

An FDR < 0.05 indicates a highly confirmative finding. An FDR > 0.5 suggest that there are false findings than true findings. The result is thus inconclusive and not reliable. An FDR between 0.05 and 0.5 cover a wide spectrum of reliability. 

FDR not good enough?

There are several ways to improve the FDR, with the expense of losing the sensitivity. You may increase the length threshold (e.g. 50 mm), increase "t" threshold (e.g. 3), increase track trimming (e.g. 2 rounds, in the advanced options). The permutation test has corrected for multiple comparison, and it is okay to run the analysis multiple time to find the best parameters that separates true findings from false findings. 


STEP 4) Visualize results and report finding


Tracks with surface rendering


1) Click on the "Show Results" button to view the tracks picked up by connectometry analysis. Two set of tracks will be presented. One is for positive correlation, and another is for negative correlation in multiple comparison.   

Another way to show the results is first openning the connectometry database (the *.db.fib.gz file) in "STEP3: fiber tracking" and loading the trk files under the [Tracts] menu. You may find the trk files placed under same directory of the demographic files.

2) Switch to "wm" as shown in the figure and click "+isosurface". You may choose full, or right half, right half...etc. to visualize full white matter or only part of it. The cut off is defined by the slice position. You can switch slices image to "T1w" to see the T1-weighted images.



3) Change the opacity of the surface in the Options window (right upper corner) under the "Surface Rendering" node.




Tracks in mosaic T1w


1) Switch slice image to "T1w".

2) In the main menu, add tracks to ROI by [Tracks][Tracks to ROI]

3) In the Options window (right upper corner), under the "Region Window" node. Change the slice layout to Mosaic 2

4) Change the contrast of the slices using the slider on the top of the region window to enhance the tracks regions.


Get SDF values along the track

Once you get tracks that show significant difference or correlation. You can open the connectometry db at [STEP 3 Fiber tracking] and load the resulting tracks. Click on track statistics and you will get SDF values along the tracks for "each subject". This allows you to plot the difference between the groups.

Publications using group connectometry analysis

  • Monnot, Cyril, Xiaoqun Zhang, Sahar Nikkhou-Aski, Peter Damberg, and Per Svenningsson. "Asymmetric dopaminergic degeneration and levodopa alter functional corticostriatal connectivity bilaterally in experimental parkinsonism." Experimental Neurology (2017).
  • Rahmani, Farzaneh, and Mohammad Hadi Aarabi. "Does apolipoprotein A1 predict microstructural changes in subgenual cingulum in early Parkinson?." Journal of Neurology: 1-10.
  • Delaparte, Lauren, Fang-Cheng Yeh, Phil Adams, Ashley Malchow, Madhukar H. Trivedi, Maria A. Oquendo, Thilo Deckersbach et al. "A comparison of structural connectivity in anxious depression versus non-anxious depression." Journal of Psychiatric Research (2017). 
  • Rahmani, Farzaneh, Mina Ansari, Atefeh Pooyan, Mehdi M. Mirbagheri, and Mohammad Hadi Aarabi. "Differences in white matter microstructure between Parkinson's disease patients with and without REM sleep behavior disorder." In Engineering in Medicine and Biology Society (EMBC), 2016 IEEE 38th Annual International Conference of the, pp. 1124-1126. IEEE, 2016.
  • Ansari, Mina, Farzaneh Rahmani, Mahsa Dolatshahi, Atefe Pooyan, and Mohammad Hadi Aarabi. "Brain pathway differences between Parkinson’s disease patients with and without REM sleep behavior disorder." Sleep and Breathing (2016): 1-7.
  • Lichenstein, Sarah D., Timothy Verstynen, and Erika E. Forbes. "Adolescent brain development and depression: A case for the importance of connectivity of the anterior cingulate cortex." Neuroscience & Biobehavioral Reviews 70 (2016): 271-287.
  • Wen, Ming-Ching, Hannah SE Heng, Samuel YE Ng, Louis CS Tan, Ling Ling Chan, and Eng King Tan. "White matter microstructural characteristics in newly diagnosed Parkinson’s disease: An unbiased whole-brain study." Scientific Reports 6 (2016): 35601.
  • Olvet, Doreen M., et al. "A comprehensive examination of white matter tracts and connectometry in major depressive disorder." Depression and anxiety(2015).
  • FC Yeh, D Badre, T Verstynen, "Connectometry: A statistical approach harnessing the analytical potential of the local connectome", - NeuroImage, 2015
  • Yeh, Fang-Cheng, Pei-Fang Tang, and Wen-Yih Isaac Tseng. "Diffusion MRI connectometry automatically reveals affected fiber pathways in individuals with chronic stroke." NeuroImage: Clinical 2 (2013): 912-921.
  • Faraji, Amir H., et al. "Longitudinal evaluation of corticospinal tract in patients with resected brainstem cavernous malformations using high-definition fiber tractography and diffusion connectometry analysis: preliminary experience."Journal of neurosurgery (2015): 1-12.
  • Abhinav, Kumar, et al. "Use of diffusion spectrum imaging in preliminary longitudinal evaluation of amyotrophic lateral sclerosis: development of an imaging biomarker." Frontiers in human neuroscience 8 (2014).

Reference:

[1] F.-C. Yeh, D. Badre, and T. Verstynen, “Connectometry: a statistical approach harnessing the analytical potential of the local connectome”, Neuroimage, accepted , 2015. (pdf)(link)
[2] F-C. Yeh, P-F. Tang, and W-Y. I. Tseng, “Diffusion MRI Connectometry Automatically Reveals Affected Fiber Pathways in Individuals with Chronic Stroke”, Neuroimage: Clinical 2 (2013) 912-921, (pdf)
[3] F. C. Yeh and W. Y. Tseng, "NTU-90: a high angular resolution brain atlas constructed by q-space diffeomorphic reconstruction," Neuroimage 58, 91-99 (2011). (pdf)
[4] Yeh F-C, Vettel JM, Singh A, Poczos B, Grafton ST, Erickson KI, et al. (2016) Quantifying Differences and Similarities in Whole-Brain White Matter Architecture Using Local Connectome Fingerprints. PLoS Comput Biol 12(11): e1005203. doi:10.1371/journal.pcbi.1005203

Ċ
Fang-Cheng Yeh,
Oct 22, 2015, 6:38 PM
Comments