{"id":4253,"date":"2024-07-01T00:00:10","date_gmt":"2024-07-01T05:00:10","guid":{"rendered":"https:\/\/my.vanderbilt.edu\/masi\/?p=4253"},"modified":"2025-12-08T16:10:38","modified_gmt":"2025-12-08T21:10:38","slug":"characterizing-patterns-of-diffusion-tensor-imaging-variance-in-aging-brains","status":"publish","type":"post","link":"https:\/\/my.vanderbilt.edu\/masi\/2024\/07\/characterizing-patterns-of-diffusion-tensor-imaging-variance-in-aging-brains\/","title":{"rendered":"Characterizing patterns of diffusion tensor imaging variance in aging brains"},"content":{"rendered":"<p>Chenyu Gao, Qi Yang, Michael E Kim, Nazirah Mohd Khairi, Leon Y Cai, Nancy R Newlin, Praitayini Kanakaraj, Lucas W Remedios, Aravind R Krishnan, Xin Yu, Tianyuan Yao, Panpan Zhang, Kurt G Schilling, Daniel Moyer, Derek B Archer, Susan M Resnick, Bennett A Landman, for the Alzheimer\u2019s Disease Neuroimaging Initiative, the BIOCARD Study team. &#8220;Characterizing patterns of diffusion tensor imaging variance in aging brains&#8221;. Journal of Medical Imaging 11 (4), 044007-044007. https:\/\/doi.org\/10.1117\/1.JMI.11.4.044007<\/p>\n<h2>Abstract<\/h2>\n<h3>Purpose<\/h3>\n<p>As large analyses merge data across sites, a deeper understanding of variance in statistical assessment across the sources of data becomes critical for valid analyses. Diffusion tensor imaging (DTI) exhibits spatially varying and correlated noise, so care must be taken with distributional assumptions. Here, we characterize the role of physiology, subject compliance, and the interaction of the subject with the scanner in the understanding of DTI variability, as modeled in the spatial variance of derived metrics in homogeneous regions.<\/p>\n<h3>Approach<\/h3>\n<p>We analyze DTI data from 1035 subjects in the Baltimore Longitudinal Study of Aging, with ages ranging from 22.4 to 103 years old. For each subject, up to 12 longitudinal sessions were conducted. We assess the variance of DTI scalars within regions of interest (ROIs) defined by four segmentation methods and investigate the relationships between the variance and covariates, including baseline age, time from the baseline (referred to as \u201cinterval\u201d), motion, sex, and whether it is the first scan or the second scan in the session.<\/p>\n<h3>Results<\/h3>\n<p>Covariate effects are heterogeneous and bilaterally symmetric across ROIs. Inter-session interval is positively related ( p \u226a 0.001 ) to FA variance in the cuneus and occipital gyrus, but negatively ( p \u226a 0.001 ) in the caudate nucleus. Males show significantly ( p \u226a 0.001 ) higher FA variance in the right putamen, thalamus, body of the corpus callosum, and cingulate gyrus. In 62 out of 176 ROIs defined by the Eve type-1 atlas, an increase in motion is associated ( p &lt; 0.05 ) with a decrease in FA variance. Head motion increases during the rescan of DTI ( \u0394 \u03bc = 0.045 mm per volume).<\/p>\n<h3>Conclusions<\/h3>\n<p>The effects of each covariate on DTI variance and their relationships across ROIs are complex. Ultimately, we encourage researchers to include estimates of variance when sharing data and consider models of heteroscedasticity in analysis. This work provides a foundation for study planning to account for regional variations in metric variance.<\/p>\n<figure id=\"attachment_4254\" aria-describedby=\"caption-attachment-4254\" style=\"width: 650px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4254 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f001-650x437.png\" alt=\"\" width=\"650\" height=\"437\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f001-650x437.png 650w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f001-300x201.png 300w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f001-768x516.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f001-1536x1032.png 1536w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f001.png 1544w\" sizes=\"auto, (max-width: 650px) 100vw, 650px\" \/><figcaption id=\"caption-attachment-4254\" class=\"wp-caption-text\">Fig. 1 Simulation shows that applying the whitening matrix to the standard linear regression equation reduces the number of false positives (FP) and false negatives (FN) under heteroskedasticity. In the top row, the population truth has zero slope. In data sampled from the synthetic population data, ordinary least square (OLS) regression using the standard equation generates FP, while the solution with whitening, W, does not falsely reject the null hypothesis (the horizontal line). After 10,000 experiments, the FP frequency is lower with whitening, centering at 5 per 100. In the second row, the population truth has a positive slope. In data sampled from the synthetic population data, OLS regression using the standard equation generates FN, while the solution with whitening, W, does not. After 10,000 experiments, the FN frequency with whitening is half that of the one without whitening.<\/figcaption><\/figure>\n<figure id=\"attachment_4255\" aria-describedby=\"caption-attachment-4255\" style=\"width: 650px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4255 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f002-650x612.png\" alt=\"\" width=\"650\" height=\"612\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f002-650x612.png 650w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f002-300x282.png 300w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f002-768x723.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f002.png 1139w\" sizes=\"auto, (max-width: 650px) 100vw, 650px\" \/><figcaption id=\"caption-attachment-4255\" class=\"wp-caption-text\">Fig. 2 We observe that the noise (approximated by the difference between the scan and rescan acquired within the same imaging session) in DTI scalar images, such as FA images, generally increases with age. (Subjects\u2019 ages are grouped into 5-year bins to respect privacy.) But motion is also considered to increase with age. We would like to know the following: Which factor is associated with DTI variance? Where and how does this association manifest?<\/figcaption><\/figure>\n<figure id=\"attachment_4256\" aria-describedby=\"caption-attachment-4256\" style=\"width: 650px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4256 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f003-650x363.png\" alt=\"\" width=\"650\" height=\"363\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f003-650x363.png 650w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f003-300x167.png 300w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f003-768x429.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f003.png 1446w\" sizes=\"auto, (max-width: 650px) 100vw, 650px\" \/><figcaption id=\"caption-attachment-4256\" class=\"wp-caption-text\">Fig. 3 The BLSA dataset we use has a slight imbalance between the number of females and males, but it is well-matched and appropriate for our research objectives in other aspects: (i) the age ranges of females and males align closely; (ii) rescan DTI data were acquired in most sessions, enabling inter-scan comparisons; and (iii) the distributions of sessions of females and males align closely.<\/figcaption><\/figure>\n<figure id=\"attachment_4257\" aria-describedby=\"caption-attachment-4257\" style=\"width: 650px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4257 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f004-650x355.png\" alt=\"\" width=\"650\" height=\"355\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f004-650x355.png 650w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f004-300x164.png 300w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f004-768x420.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f004-1536x840.png 1536w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f004.png 1556w\" sizes=\"auto, (max-width: 650px) 100vw, 650px\" \/><figcaption id=\"caption-attachment-4257\" class=\"wp-caption-text\">Fig. 4 Brain segmentation labels are obtained using the SLANT segmentation of the target subject\u2019s T1w image and using three types of manual parcellations provided by the Eve atlas. To generate transformation matrices for transferring these labels to DTI scalar images, intra- and inter-modality registrations are performed. Standard deviations of DTI scalars within each ROI are computed.<\/figcaption><\/figure>\n<figure id=\"attachment_4258\" aria-describedby=\"caption-attachment-4258\" style=\"width: 526px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4258 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f005-526x650.png\" alt=\"\" width=\"526\" height=\"650\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f005-526x650.png 526w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f005-243x300.png 243w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f005-768x949.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f005-1243x1536.png 1243w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f005-1657x2048.png 1657w\" sizes=\"auto, (max-width: 526px) 100vw, 526px\" \/><figcaption id=\"caption-attachment-4258\" class=\"wp-caption-text\">Fig. 5 Covariate effects on FA standard deviation (standardized) are region-specific. Motion and interval exhibit opposite effect directions in many ROIs. Gender differences exist in multiple ROIs. Counterintuitively, motion is negatively related to FA standard deviation in many ROIs. The lookup table for the abbreviation of the ROI name is in the Supplementary Material 2.<\/figcaption><\/figure>\n<figure id=\"attachment_4259\" aria-describedby=\"caption-attachment-4259\" style=\"width: 526px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4259 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f006-526x650.png\" alt=\"\" width=\"526\" height=\"650\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f006-526x650.png 526w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f006-243x300.png 243w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f006-768x949.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f006-1243x1536.png 1243w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f006-1657x2048.png 1657w\" sizes=\"auto, (max-width: 526px) 100vw, 526px\" \/><figcaption id=\"caption-attachment-4259\" class=\"wp-caption-text\">Fig. 6 The region-specific and bidirectional patterns of covariate effects are similarly observed in the results derived from SLANT segmentation, despite its differing definitions and delineations of ROIs compared to Eve type-1 segmentation (Fig. 5).<\/figcaption><\/figure>\n<figure id=\"attachment_4260\" aria-describedby=\"caption-attachment-4260\" style=\"width: 631px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4260 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f007-631x650.png\" alt=\"\" width=\"631\" height=\"650\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f007-631x650.png 631w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f007-291x300.png 291w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f007-768x791.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f007-1491x1536.png 1491w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f007.png 1500w\" sizes=\"auto, (max-width: 631px) 100vw, 631px\" \/><figcaption id=\"caption-attachment-4260\" class=\"wp-caption-text\">Fig. 7 Despite the different definitions and delineations of ROIs between Eve type-1 and SLANT segmentations, results based on the two segmentation methods are largely similar (comparable regions are colored similarly) and both show that the effects of motion and interval on FA variance vary across ROIs.<\/figcaption><\/figure>\n<figure id=\"attachment_4261\" aria-describedby=\"caption-attachment-4261\" style=\"width: 650px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4261 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f008-650x326.png\" alt=\"\" width=\"650\" height=\"326\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f008-650x326.png 650w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f008-300x151.png 300w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f008-768x385.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f008-1536x771.png 1536w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f008-2048x1028.png 2048w\" sizes=\"auto, (max-width: 650px) 100vw, 650px\" \/><figcaption id=\"caption-attachment-4261\" class=\"wp-caption-text\">Fig. 8 SNR of the FA images across eight Eve type-1 atlas-defined ROIs, including white matter (WM) regions (body of the corpus callosum, precuneus WM, cuneus WM), gray matter regions (insular, putamen, middle occipital gyrus, thalamus), and mixed regions (medulla), in four different scanners of BLSA (where scanner A is the 1.5 Tesla Philips Intera scanner and scanners B\/C\/D are the 3 Tesla Philips Achieva scanners). \u201cL.\u201d stands for the left hemisphere of the brain.<\/figcaption><\/figure>\n<figure id=\"attachment_4262\" aria-describedby=\"caption-attachment-4262\" style=\"width: 526px\" class=\"wp-caption alignnone\"><img loading=\"lazy\" decoding=\"async\" class=\"wp-image-4262 size-large\" src=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f009-526x650.png\" alt=\"\" width=\"526\" height=\"650\" srcset=\"https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f009-526x650.png 526w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f009-243x300.png 243w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f009-768x950.png 768w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f009-1242x1536.png 1242w, https:\/\/cdn.vanderbilt.edu\/t2-my\/my-prd\/wp-content\/uploads\/sites\/2304\/2025\/12\/JMI_11_4_044007_f009-1656x2048.png 1656w\" sizes=\"auto, (max-width: 526px) 100vw, 526px\" \/><figcaption id=\"caption-attachment-4262\" class=\"wp-caption-text\">Fig. 9 To assess the generalizability of our findings, we include two additional datasets, ADNI and BIOCARD, and fit the linear mixed-effects models. The coefficients and p-values show similar patterns to those from BLSA alone, despite that the effect sizes and the hierarchical clustering are different, partly due to the omission of the rescan covariate.<\/figcaption><\/figure>\n","protected":false},"excerpt":{"rendered":"<p>Chenyu Gao, Qi Yang, Michael E Kim, Nazirah Mohd Khairi, Leon Y Cai, Nancy R Newlin, Praitayini Kanakaraj, Lucas W Remedios, Aravind R Krishnan, Xin Yu, Tianyuan Yao, Panpan Zhang, Kurt G Schilling, Daniel Moyer, Derek B Archer, Susan M Resnick, Bennett A Landman, for the Alzheimer\u2019s Disease Neuroimaging Initiative, the BIOCARD Study team. &#8220;Characterizing&#8230;<\/p>\n","protected":false},"author":9582,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[27,201,33,9,116,144,46,44],"tags":[],"class_list":["post-4253","post","type-post","status-publish","format-standard","hentry","category-big-data","category-characterization","category-diffusion-tensor-imaging","category-diffusion-weighted-mri","category-image-processing","category-longitudinal","category-multi-atlas-segmentation","category-noise-estimation"],"_links":{"self":[{"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/posts\/4253","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/users\/9582"}],"replies":[{"embeddable":true,"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/comments?post=4253"}],"version-history":[{"count":1,"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/posts\/4253\/revisions"}],"predecessor-version":[{"id":4263,"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/posts\/4253\/revisions\/4263"}],"wp:attachment":[{"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/media?parent=4253"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/categories?post=4253"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/my.vanderbilt.edu\/masi\/wp-json\/wp\/v2\/tags?post=4253"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}