Publications by Year: 2021

Diffusion tensor imaging (DTI) consistently detects increased mean diffusivity and decreased fractional anisotropy with advancing age in regions of primarily single white matter (WM) fiber populations, but findings have been inconsistent in regions of more complex fiber architecture. Given that DTI remains more common for characterizing aging WM than advanced diffusion MRI models due to DTI’s simplicity, robustness, and efficiency, it is critical to strive to maximize the information extracted from DTI across the entire WM. The present study uses an orthogonal diffusion tensor decomposition based on the 3 eigenvalue moments (mean diffusivity, norm of anisotropy, and mode of anisotropy), yielding clear voxelwise degeneration patterns across the WM, including regions of complex fiber architecture. This indicates that the previous challenges of DTI in these regions were due to the choice of tensor decomposition rather than the DTI model itself. This study therefore presents a revised view of DTI of aging WM and indicates how age-related degeneration in complex fiber architecture can manifest in forms other than decreased fractional anisotropy.

Segmentation of brain tissue types from diffusion MRI (dMRI) is an important task, required for quantification of brain microstructure and for improving tractography. Current dMRI segmentation is mostly based on anatomical MRI (e.g., T1- and T2-weighted) segmentation that is registered to the dMRI space. However, such inter-modality registration is challenging due to more image distortions and lower image resolution in dMRI as compared with anatomical MRI. In this study, we present a deep learning method for diffusion MRI segmentation, which we refer to as DDSeg. Our proposed method learns tissue segmentation from high-quality imaging data from the Human Connectome Project (HCP), where registration of anatomical MRI to dMRI is more precise. The method is then able to predict a tissue segmentation directly from new dMRI data, including data collected with different acquisition protocols, without requiring anatomical data and inter-modality registration. We train a convolutional neural network (CNN) to learn a tissue segmentation model using a novel augmented target loss function designed to improve accuracy in regions of tissue boundary. To further improve accuracy, our method adds diffusion kurtosis imaging (DKI) parameters that characterize non-Gaussian water molecule diffusion to the conventional diffusion tensor imaging parameters. The DKI parameters are calculated from the recently proposed mean-kurtosis-curve method that corrects implausible DKI parameter values and provides additional features that discriminate between tissue types. We demonstrate high tissue segmentation accuracy on HCP data, and also when applying the HCP-trained model on dMRI data from other acquisitions with lower resolution and fewer gradient directions.

Younger age at first exposure (AFE) to repetitive head impacts while playing American football increases the risk for later-life neuropsychological symptoms and brain alterations. However, it is not known whether AFE is associated with cortical thickness in American football players. Sixty-three former professional National Football League players (55.5 ± 7.7 years) with cognitive, behavioral, and mood symptoms underwent neuroimaging and neuropsychological testing. First, the association between cortical thickness and AFE was tested. Second, the relationship between clusters of decreased cortical thickness and verbal and visual memory, and composite measures of mood/behavior and attention/psychomotor speed was assessed. AFE was positively correlated with cortical thickness in the right superior frontal cortex (cluster-wise P value [CWP] = 0.0006), the left parietal cortex (CWP = 0.0003), and the occipital cortices (right: CWP = 0.0023; left: CWP = 0.0008). A positive correlation was found between cortical thickness of the right superior frontal cortex and verbal memory (R = 0.333, P = 0.019), and the right occipital cortex and visual memory (R = 0.360, P = 0.012). In conclusion, our results suggest an association between younger AFE and decreased cortical thickness, which in turn is associated with worse neuropsychological performance. Furthermore, an association between younger AFE and signs of neurodegeneration later in life in symptomatic former American football players seems likely.

INTRODUCTION: Recent network-based analyses suggest that schizophrenia symptoms are intricately connected and interdependent, such that central symptoms can activate adjacent symptoms and increase global symptom burden. Here, we sought to identify key clinical and neurobiological factors that relate to symptom organization in established schizophrenia. METHODS: A symptom comorbidity network was mapped for a broad constellation of symptoms measured in 642 individuals with a schizophrenia-spectrum disorder. Centrality analyses were used to identify hub symptoms. The extent to which each patient’s symptoms formed clusters in the comorbidity network was quantified with cluster analysis and used to predict (1) clinical features, including illness duration and psychosis (positive symptom) severity and (2) brain white matter microstructure, indexed by the fractional anisotropy (FA), in a subset (n = 296) of individuals with diffusion-weighted imaging (DWI) data. RESULTS: Global functioning, substance use, and blunted affect were the most central symptoms within the symptom comorbidity network. Symptom profiles for some patients formed highly interconnected clusters, whereas other patients displayed unrelated and disconnected symptoms. Stronger clustering among an individual’s symptoms was significantly associated with shorter illness duration (t = 2.7; P = .0074), greater psychosis severity (ie, positive symptoms expression) (t = -5.5; P < 0.0001) and lower fractional anisotropy in fibers traversing the cortico-cerebellar-thalamic-cortical circuit (r = .59, P < 0.05). CONCLUSION: Symptom network structure varies over the course of schizophrenia: symptom interactions weaken with increasing illness duration and strengthen during periods of high positive symptom expression. Reduced white matter coherence relates to stronger symptom clustering, and thus, may underlie symptom cascades and global symptomatic burden in individuals with schizophrenia.

As the COVID-19 pandemic impacts global populations, computed tomography (CT) lung imaging is being used in many countries to help manage patient care as well as to rapidly identify potentially useful quantitative COVID-19 CT imaging biomarkers. Quantitative COVID-19 CT imaging applications, typically based on computer vision modeling and artificial intelligence algorithms, include the potential for better methods to assess COVID-19 extent and severity, assist with differential diagnosis of COVID-19 versus other respiratory conditions, and predict disease trajectory. To help accelerate the development of robust quantitative imaging algorithms and tools, it is critical that CT imaging is obtained following best practices of the quantitative lung CT imaging community. Toward this end, the Radiological Society of North America’s (RSNA) Quantitative Imaging Biomarkers Alliance (QIBA) CT Lung Density Profile Committee and CT Small Lung Nodule Profile Committee developed a set of best practices to guide clinical sites using quantitative imaging solutions and to accelerate the international development of quantitative CT algorithms for COVID-19. This guidance document provides quantitative CT lung imaging recommendations for COVID-19 CT imaging, including recommended CT image acquisition settings for contemporary CT scanners. Additional best practice guidance is provided on scientific publication reporting of quantitative CT imaging methods and the importance of contributing COVID-19 CT imaging datasets to open science research databases.

The brain hemispheres can be divided into an upper dorsal and a lower ventral system. Each system consists of distinct cortical regions connected via long association tracts. The tracts cross the central sulcus or the limen insulae to connect the frontal lobe with the posterior brain. The dorsal stream is associated with sensorimotor mapping. The ventral stream serves structural analysis and semantics in different domains, as visual, acoustic or space processing. How does the prefrontal cortex, regarded as the platform for the highest level of integration, incorporate information from these different domains? In the current view, the ventral pathway consists of several separate tracts, related to different modalities. Originally the assumption was that the ventral path is a continuum, covering all modalities. The latter would imply a very different anatomical basis for cognitive and clinical models of processing. To further define the ventral connections, we used cutting-edge in vivo global tractography on high-resolution diffusion tensor imaging (DTI) data from 100 normal subjects from the human connectome project and ex vivo preparation of fiber bundles in the extreme capsule of 8 humans using the Klingler technique. Our data showed that ventral stream tracts, traversing through the extreme capsule, form a continuous band of fibers that fan out anteriorly to the prefrontal cortex, and posteriorly to temporal, occipital and parietal cortical regions. Introduction of additional volumes of interest in temporal and occipital lobes differentiated between the inferior fronto-occipital fascicle (IFOF) and uncinate fascicle (UF). Unequivocally, in both experiments, in all subjects a connection between the inferior frontal and middle-to-posterior temporal cortical region, otherwise known as the temporo-frontal extreme capsule fascicle (ECF) from nonhuman primate brain-tracing experiments was identified. In the human brain, this tract connects the language domains of Broca’s area and Wernicke’s area. The differentiation in the three tracts, IFOF, UF and ECF seems arbitrary, all three pass through the extreme capsule. Our data show that the ventral pathway represents a continuum. The three tracts merge seamlessly and streamlines showed considerable overlap in their anterior and posterior course. Terminal maps identified prefrontal cortex in the frontal lobe and association cortices in temporal, occipital and parietal lobes as streamline endings. This anatomical substrate potentially facilitates the prefrontal cortex to integrate information across different domains and modalities.

Purpose: The aim of this study is to assess the role of traction bronchiectasis/bronchiolectasis and its progression as a predictor for early fibrosis in interstitial lung abnormalities (ILA). Methods: Three hundred twenty-seven ILA participants out of 5764 in the Age, Gene/Environment Susceptibility (AGES)-Reykjavik Study who had undergone chest CT twice with an interval of approximately five-years were enrolled in this study. Traction bronchiectasis/bronchiolectasis index (TBI) was classified on a four-point scale: 0, ILA without traction bronchiectasis/bronchiolectasis; 1, ILA with bronchiolectasis but without bronchiectasis or architectural distortion; 2, ILA with mild to moderate traction bronchiectasis; 3, ILA and severe traction bronchiectasis and/or honeycombing. Traction bronchiectasis (TB) progression was classified on a five-point scale: 1, Improved; 2, Probably improved; 3, No change; 4, Probably progressed; 5, Progressed. Overall survival (OS) among participants with different TB Progression Score and between the TB progression group and No TB progression group was also investigated. Hazard radio (HR) was estimated with Cox proportional hazards model.

OBJECTIVE: To investigate the relative predictive value of circulating immune cell markers for cardiovascular mortality in ambulatory adults without cardiovascular disease. METHODS: We analyzed data of participants enrolled in the National Health and Nutrition Examination Survey from January 1, 1999, to December 31, 2010, with the total leukocyte count within a normal range (4000-11,000 cells/μL [to convert to cells ×10/L, multiply by 0.001]) and without cardiovascular disease. The relative predictive value of circulating immune cell markers measured at enrollment-including total leukocyte count, absolute neutrophil count, absolute lymphocyte count, absolute monocyte count, monocyte-lymphocyte ratio (MLR), neutrophil-lymphocyte ratio, and C-reactive protein-for cardiovascular mortality was evaluated. The marker with the best predictive value was added to the 10-year atherosclerotic cardiovascular disease (ASCVD) risk score to estimate net risk reclassification indices for 10-year cardiovascular mortality. RESULTS: Among 21,599 participants eligible for this analysis, the median age was 47 years (interquartile range, 34-63 years); 10,651 (49.2%) participants were women, and 10,713 (49.5%) were self-reported non-Hispanic white. During a median follow-up of 9.6 years (interquartile range, 6.8-13.1 years), there were 627 cardiovascular deaths. MLR had the best predictive value for cardiovascular mortality. The addition of elevated MLR (>=0.3) to the 10-year ASCVD risk score improved the classification by 2.7%±1.4% (P=.04). Elevated MLR had better predictive value than C-reactive protein and several components of the 10-year ASCVD risk score. CONCLUSION: Among ambulatory US adults without preexisting cardiovascular disease, we found that MLR had the best predictive value for cardiovascular mortality among circulating immune markers. The addition of MLR to the 10-year risk score significantly improved the risk classification of participants.

PURPOSE: To introduce, develop, and evaluate a novel denoising technique for diffusion MRI that leverages nonlinear redundancy in the data to boost the SNR while preserving signal information. METHODS: We exploit nonlinear redundancy of the dMRI data by means of kernel principal component analysis (KPCA), a nonlinear generalization of PCA to reproducing kernel Hilbert spaces. By mapping the signal to a high-dimensional space, a higher level of redundant information is exploited, thereby enabling better denoising than linear PCA. We implement KPCA with a Gaussian kernel, with parameters automatically selected from knowledge of the noise statistics, and validate it on realistic Monte Carlo simulations as well as with in vivo human brain submillimeter and low-resolution dMRI data. We also demonstrate KPCA denoising on multi-coil dMRI data. RESULTS: SNR improvements up to 2.7 were obtained in real in vivo datasets denoised with KPCA, in comparison to SNR gains of up to 1.8 using a linear PCA denoising technique called Marchenko-Pastur PCA (MPPCA). Compared to gold-standard dataset references created from averaged data, we showed that lower normalized root mean squared error was achieved with KPCA compared to MPPCA. Statistical analysis of residuals shows that anatomical information is preserved and only noise is removed. Improvements in the estimation of diffusion model parameters such as fractional anisotropy, mean diffusivity, and fiber orientation distribution functions were also demonstrated. CONCLUSION: Nonlinear redundancy of the dMRI signal can be exploited with KPCA, which allows superior noise reduction/SNR improvements than the MPPCA method, without loss of signal information.

PURPOSE: Reperfusion therapy enables effective treatment of ischemic stroke presenting within 4-6 hours. However, tissue progression from ischemia to infarction is variable, and some patients benefit from treatment up until 24 hours. Improved imaging techniques are needed to identify these patients. Here, it was hypothesized that time dependence in diffusion MRI may predict tissue outcome in ischemic stroke. METHODS: Diffusion MRI data were acquired with multiple diffusion times in five non-reperfused patients at 2, 9, and 100 days after stroke onset. Maps of "rate of kurtosis change" (k), mean kurtosis, ADC, and fractional anisotropy were derived. The ADC maps defined lesions, normal-appearing tissue, and the lesion tissue that would either be infarcted or remain viable by day 100. Diffusion parameters were compared (1) between lesions and normal-appearing tissue, and (2) between lesion tissue that would be infarcted or remain viable. RESULTS: Positive values of k were observed within stroke lesions on day 2 (P = .001) and on day 9 (P = .023), indicating diffusional exchange. On day 100, high ADC values indicated infarction of 50 ± 20% of the lesion volumes. Tissue infarction was predicted by high k values both on day 2 (P = .026) and on day 9 (P = .046), by low mean kurtosis values on day 2 (P = .043), and by low fractional anisotropy values on day 9 (P = .029), but not by low ADC values. CONCLUSIONS: Diffusion time dependence predicted tissue outcome in ischemic stroke more accurately than the ADC, and may be useful for predicting reperfusion benefit.