Publications

2019

Jimenez-Carretero D, aez DB-P, Nardelli P, Fraga P, Fraile E, epar R ul SJ e E, Ledesma-Carbayo MJ. A graph-cut approach for pulmonary artery-vein segmentation in noncontrast CT images. Med Image Anal. 2019;52:144–159. doi:10.1016/j.media.2018.11.011
Lung vessel segmentation has been widely explored by the biomedical image processing community; however, the differentiation of arterial from venous irrigation is still a challenge. Pulmonary artery-vein (AV) segmentation using computed tomography (CT) is growing in importance owing to its undeniable utility in multiple cardiopulmonary pathological states, especially those implying vascular remodelling, allowing the study of both flow systems separately. We present a new framework to approach the separation of tree-like structures using local information and a specifically designed graph-cut methodology that ensures connectivity as well as the spatial and directional consistency of the derived subtrees. This framework has been applied to the pulmonary AV classification using a random forest (RF) pre-classifier to exploit the local anatomical differences of arteries and veins. The evaluation of the system was performed using 192 bronchopulmonary segment phantoms, 48 anthropomorphic pulmonary CT phantoms, and 26 lungs from noncontrast CT images with precise voxel-based reference standards obtained by manually labelling the vessel trees. The experiments reveal a relevant improvement in the accuracy ( \~ 20%) of the vessel particle classification with the proposed framework with respect to using only the pre-classification based on local information applied to the whole area of the lung under study. The results demonstrated the accurate differentiation between arteries and veins in both clinical and synthetic cases, specifically when the image quality can guarantee a good airway segmentation, which opens a huge range of possibilities in the clinical study of cardiopulmonary diseases.
Seitz J, Kubicki M, Jacobs EG, Cherkerzian S, Weiss BK, Papadimitriou G, Mouradian P, Buka S, Goldstein JM, Makris N. Impact of sex and reproductive status on memory circuitry structure and function in early midlife using structural covariance analysis. Hum Brain Mapp. 2019;40(4):1221–1233. doi:10.1002/hbm.24441
Research on age-related memory alterations traditionally targets individuals aged >=65 years. However, recent studies emphasize the importance of early aging processes. We therefore aimed to characterize variation in brain gray matter structure in early midlife as a function of sex and menopausal status. Subjects included 94 women (33 premenopausal, 29 perimenopausal, and 32 postmenopausal) and 99 demographically comparable men from the New England Family Study. Subjects were scanned with a high-resolution T1 sequence on a 3 T whole body scanner. Sex and reproductive-dependent structural differences were evaluated using Box’s M test and analysis of covariances (ANCOVAs) for gray matter volumes. Brain regions of interest included dorsolateral prefrontal cortex (DLPFC), inferior parietal lobule (iPAR), anterior cingulate cortex (ACC), hippocampus (HIPP), and parahippocampus. While we observed expected significant sex differences in volume of hippocampus with women of all groups having higher volumes than men relative to cerebrum size, we also found significant differences in the covariance matrices of perimenopausal women compared with postmenopausal women. Associations between ACC and HIPP/iPAR/DLPFC were higher in postmenopausal women and correlated with better memory performance. Findings in this study underscore the importance of sex and reproductive status in early midlife for understanding memory function with aging.
Topgaard D. Diffusion tensor distribution imaging. NMR Biomed. 2019;32(5):e4066. doi:10.1002/nbm.4066
Conventional diffusion MRI yields voxel-averaged parameters that suffer from ambiguities for heterogeneous anisotropic materials such as brain tissue. Using principles from solid-state NMR spectroscopy, we have previously introduced the shape of the diffusion encoding tensor as a separate acquisition dimension that disentangles isotropic and anisotropic contributions to the observed diffusivities, thereby allowing for unconstrained data inversion into diffusion tensor distributions with "size," "shape," and orientation dimensions. Here we combine our recent non-parametric data inversion algorithm and data acquisition protocol with an imaging pulse sequence to demonstrate spatial mapping of diffusion tensor distributions using a previously developed composite phantom with multiple isotropic and anisotropic components. We propose a compact format for visualizing two-dimensional arrays of the distributions, new scalar parameters quantifying intra-voxel heterogeneity, and a binning procedure giving maps of all relevant parameters for each of the components resolved in the multidimensional distribution space.
Tobisch A, Schultz T, Stirnberg R, Varela-Mattatall G, Knutsson H, azaval PI, Stöcker T. Comparison of basis functions and q-space sampling schemes for robust compressed sensing reconstruction accelerating diffusion spectrum imaging. NMR Biomed. 2019;32(3):e4055. doi:10.1002/nbm.4055
Time constraints placed on magnetic resonance imaging often restrict the application of advanced diffusion MRI (dMRI) protocols in clinical practice and in high throughput research studies. Therefore, acquisition strategies for accelerated dMRI have been investigated to allow for the collection of versatile and high quality imaging data, even if stringent scan time limits are imposed. Diffusion spectrum imaging (DSI), an advanced acquisition strategy that allows for a high resolution of intra-voxel microstructure, can be sufficiently accelerated by means of compressed sensing (CS) theory. CS theory describes a framework for the efficient collection of fewer samples of a data set than conventionally required followed by robust reconstruction to recover the full data set from sparse measurements. For an accurate recovery of DSI data, a suitable acquisition scheme for sparse q-space sampling and the sensing and sparsifying bases for CS reconstruction need to be selected. In this work we explore three different types of q-space undersampling schemes and two frameworks for CS reconstruction based on either Fourier or SHORE basis functions. After CS recovery, diffusion and microstructural parameters and orientational information are estimated from the reconstructed data by means of state-of-the-art processing techniques for dMRI analysis. By means of simulation, diffusion phantom and in vivo DSI data, an isotropic distribution of q-space samples was found to be optimal for sparse DSI. The CS reconstruction results indicate superior performance of Fourier-based CS-DSI compared to the SHORE-based approach. Based on these findings we outline an experimental design for accelerated DSI and robust CS reconstruction of the sparse measurements that is suitable for the application within time-limited studies.
en MN, Schremm A, Nilsson M, Horne M, Roll M. Cortical thickness of Broca’s area and right homologue is related to grammar learning aptitude and pitch discrimination proficiency. Brain Lang. 2019;188:42–47. doi:10.1016/j.bandl.2018.12.002
Aptitude for and proficiency in acquiring new languages varies in the human population but their neural bases are largely unknown. We investigated the influence of cortical thickness on language learning predictors measured by the LLAMA tests and a pitch-change discrimination test. The LLAMA tests are first language-independent assessments of language learning aptitude for vocabulary, phonetic working memory, sound-symbol correspondence (not used in this study), and grammatical inferencing. Pitch perception proficiency is known to predict aptitude for learning new phonology. Results show a correlation between scores in a grammatical meaning-inferencing aptitude test and cortical thickness of Broca’s area (r(30) = 0.65, p = 0.0202) and other frontal areas (r(30) = 0.66, p = 0.0137). Further, a correlation was found between proficiency in discriminating pitch-change direction and cortical thickness of the right Broca homologue (r(30) = 0.57, p = 0.0006). However, no correlations were found for aptitude for vocabulary learning or phonetic working memory. Results contribute to locating cortical regions important for language-learning aptitude.
Del Re EC, Bouix S, Fitzsimmons J, Blokland G ella AM, Mesholam-Gately R, Wojcik J, Kikinis Z, Kubicki M, Petryshen T, Pasternak O, et al. Diffusion abnormalities in the corpus callosum in first episode schizophrenia: Associated with enlarged lateral ventricles and symptomatology. Psychiatry Res. 2019;277:45–51. doi:10.1016/j.psychres.2019.02.038
INTRODUCTION: Abnormalities in the corpus callosum (CC) and the lateral ventricles (LV) are hallmark features of schizophrenia. These abnormalities have been reported in chronic and in first episode schizophrenia (FESZ). Here we explore further associations between CC and LV in FESZ using diffusion tensor imaging (DTI). METHODS: . Sixteen FESZ patients and 16 healthy controls (HC), matched on age, gender, and handedness participated in the study. Diffusion and structural imaging scans were acquired on a 3T GE Signa magnet. Volumetric measures for LV and DTI measures for five CC subdivisions were completed in both groups. In addition, two-tensor tractography, the latter corrected for free-water (FA), was completed for CC. Correlations between LV and DTI measures of the CC were examined in both groups, while correlations between DTI and clinical measures were examined in only FESZ. RESULTS: Results from two-tensor tractography demonstrated decreased FA and increased trace and radial diffusivity (RD) in the five CC subdivisions in FESZ compared to HC. Central CC diffusion measures in FESZ were significantly correlated with volume of the LV, i.e., decreased FA values were associated with larger LV volume, while increased RD and trace values were associated with larger LV volume. In controls, correlations were also significant, but they were in the opposite direction from FESZ. In addition, decreased FA in FESZ was associated with more positive symptoms. DISCUSSION: Partial volume corrected FA, RD, and trace abnormalities in the CC in FESZ suggest possible de- or dys-myelination, or changes in axonal diameters, all compatible with neurodevelopmental theories of schizophrenia. Correlational findings between the volume of LV and diffusion measures in FESZ reinforce the concept of a link between abnormalities in the LV and CC in early stages of schizophrenia and are also compatible with neurodevelopmental abnormalities in this population.
Hahn A, Strandberg TO, Stomrud E, Nilsson M, van Westen D, Palmqvist S, Ossenkoppele R, Hansson O. Association Between Earliest Amyloid Uptake and Functional Connectivity in Cognitively Unimpaired Elderly. Cereb Cortex. 2019;29(5):2173–2182. doi:10.1093/cercor/bhz020
Alterations in cognitive performance have been noted in nondemented subjects with elevated accumulation of amyloid-β (Aβ) fibrils. However, it is not yet understood whether brain function is already influenced by Aβ deposition during the very earliest stages of the disease. We therefore investigated associations between [18F]Flutemetamol PET, resting-state functional connectivity, gray and white matter structure and cognitive performance in 133 cognitively normal elderly that exhibited normal global Aβ PET levels. [18F]Flutemetamol uptake in regions known to accumulate Aβ fibrils early in preclinical AD (i.e., mainly certain parts of the default-mode network) was positively associated with dynamic but not static functional connectivity (r = 0.77). Dynamic functional connectivity was further related to better cognitive performance (r = 0.21-0.72). No significant associations were found for Aβ uptake with gray matter volume or white matter diffusivity. The findings demonstrate that the earliest accumulation of Aβ fibrils is associated with increased functional connectivity, which occurs before any structural alterations. The enhanced functional connectivity may reflect a compensatory mechanism to maintain high cognitive performance in the presence of increasing amyloid accumulation during the earliest phases of AD.
Synn AJ, Zhang C, Washko GR, epar R ul SJ e E, O\textquoterightConnor GT, Li W, Mittleman MA, Rice MB. Cigarette Smoke Exposure and Radiographic Pulmonary Vascular Morphology in the Framingham Heart Study. Ann Am Thorac Soc. 2019;16(6):698–706. doi:10.1513/AnnalsATS.201811-795OC
Cigarette smoke exposure is a risk factor for many lung diseases, and histologic studies suggest that tobacco-related vasoconstriction and vessel loss plays a role in the development of emphysema. However, it remains unclear how tobacco affects the pulmonary vasculature in general populations with a typical range of tobacco exposure, and whether these changes are detectable by radiographic methods. To determine whether tobacco exposure in a generally healthy population manifests as lower pulmonary blood vessel volumes and vascular pruning on imaging. A total of 2,410 Framingham Heart Study participants with demographic data and smoking history underwent volumetric whole-lung computed tomography from 2008 to 2011. Automated algorithms calculated the total blood volume of all intrapulmonary vessels (TBV), smaller peripheral vessels (defined as cross-sectional area
Herberthson M, Yolcu C, Knutsson H, Westin C-F, Özarslan E. Orientationally-averaged diffusion-attenuated magnetic resonance signal for locally-anisotropic diffusion. Sci Rep. 2019;9(1):4899. doi:10.1038/s41598-019-41317-8
Diffusion-attenuated MR signal for heterogeneous media has been represented as a sum of signals from anisotropic Gaussian sub-domains to the extent that this approximation is permissible. Any effect of macroscopic (global or ensemble) anisotropy in the signal can be removed by averaging the signal values obtained by differently oriented experimental schemes. The resulting average signal is identical to what one would get if the micro-domains are isotropically (e.g., randomly) distributed with respect to orientation, which is the case for "powdered" specimens. We provide exact expressions for the orientationally-averaged signal obtained via general gradient waveforms when the microdomains are characterized by a general diffusion tensor possibly featuring three distinct eigenvalues. This extends earlier results which covered only axisymmetric diffusion as well as measurement tensors. Our results are expected to be useful in not only multidimensional diffusion MR but also solid-state NMR spectroscopy due to the mathematical similarities in the two fields.