# Publications

2022
Chakwizira A, Westin C-F, Brabec J, Lasič S, Knutsson L, Szczepankiewicz F, Nilsson M. Diffusion MRI With Pulsed and Free Gradient Waveforms: Effects of Restricted Diffusion and Exchange. NMR Biomed. 2022 :e4827.Abstract
Monitoring time dependence with diffusion MRI yields observables sensitive to compartment sizes (restricted diffusion) and membrane permeability (water exchange). However, restricted diffusion and exchange have opposite effects on the diffusion-weighted signal, which can lead to errors in parameter estimates. In this work, we propose a signal representation that incorporates the effects of both restricted diffusion and exchange up to second order in b-value and is compatible with gradient waveforms of arbitrary shape. The representation features mappings from a gradient waveform to two scalars that separately control the sensitivity to restriction and exchange. We demonstrate that these scalars span a two-dimensional space that can be used to choose waveforms that selectively probe restricted diffusion or exchange, eliminating the correlation between the two phenomena. We found that waveforms with specific but unconventional shapes provide an advantage over conventional pulsed and oscillating gradient acquisitions. We also show that parametrization of waveforms into a two-dimensional space can be used to understand protocols from other approaches that probe restricted diffusion and exchange. For example, we found that the variation of mixing time in filter-exchange imaging corresponds to variation of our exchange-weighting scalar at a fixed value of the restriction-weighting scalar. The proposed signal representation was evaluated using Monte Carlo simulations in identical parallel cylinders with hexagonal and random packing as well as parallel cylinders with gamma-distributed radii. Results showed that the approach is sensitive to sizes in the interval 4-12 μm $$\upmu \mathrm{m}$$ and exchange rates in the simulated range of 0 to 20 s - 1 $${\mathrm{s}}^{-1}$$ , but also that there is a sensitivity to the extracellular geometry. The presented theory constitutes a simple and intuitive description of how restricted diffusion and exchange influence the signal as well as a guide to protocol design capable of separating the two effects.
Mayer C, Nägele FL, Petersen M, Frey BM, Hanning U, Pasternak O, Petersen E, Gerloff C, Thomalla G, Cheng B. Free-Water Diffusion MRI Detects Structural Alterations Surrounding White Matter Hyperintensities in the Early Stage of Cerebral Small Vessel Disease. J Cereb Blood Flow Metab. 2022;42 (9) :1707-18.Abstract
In cerebral small vessel disease (CSVD), both white matter hyperintensities (WMH) of presumed vascular origin and the normal-appearing white matter (NAWM) contain microstructural brain alterations on diffusion-weighted MRI (DWI). Contamination of DWI-derived metrics by extracellular free-water can be corrected with free-water (FW) imaging. We investigated the alterations in FW and FW-corrected fractional anisotropy (FA-t) in WMH and surrounding tissue and their association with cerebrovascular risk factors. We analysed 1,000 MRI datasets from the Hamburg City Health Study. DWI was used to generate FW and FA-t maps. WMH masks were segmented on FLAIR and T1-weighted MRI and dilated repeatedly to create 8 NAWM masks representing increasing distance from WMH. Linear models were applied to compare FW and FA-t across WMH and NAWM masks and in association with cerebrovascular risk. Median age was 64 ± 14 years. FW and FA-t were altered 8 mm and 12 mm beyond WMH, respectively. Smoking was significantly associated with FW in NAWM (p = 0.008) and FA-t in WMH (p = 0.008) and in NAWM (p = 0.003) while diabetes and hypertension were not. Further research is necessary to examine whether FW and FA-t alterations in NAWM are predictors for developing WMH.
Laurent D, Riek J, Sinclair CDJ, Houston P, Roubenoff R, Papanicolaou DA, Nagy A, Pieper S, Yousry TA, Hanna MG, et al. Longitudinal Changes in MRI Muscle Morphometry and Composition in People With Inclusion Body Myositis. Neurology. 2022;99 (9) :e865-e876.Abstract
BACKGROUND AND OBJECTIVES: Limited data suggest that quantitative MRI (qMRI) measures have potential to be used as trial outcome measures in sporadic inclusion body myositis (sIBM) and as a noninvasive assessment tool to study sIBM muscle pathologic processes. Our aim was to evaluate changes in muscle structure and composition using a comprehensive multiparameter set of qMRI measures and to assess construct validity and responsiveness of qMRI measures in people with sIBM. METHODS: This was a prospective observational cohort study with assessments at baseline (n = 30) and 1 year (n = 26). qMRI assessments include thigh muscle volume (TMV), inter/intramuscular adipose tissue (IMAT), muscle fat fraction (FF), muscle inflammation (T2 relaxation time), IMAT from T2* relaxation (T2*-IMAT), intermuscular connective tissue from T2* relaxation (T2*-IMCT), and muscle macromolecular structure from the magnetization transfer ratio (MTR). Physical performance assessments include sIBM Physical Functioning Assessment (sIFA), 6-minute walk distance, and quantitative muscle testing of the quadriceps. Correlations were assessed using the Spearman correlation coefficient. Responsiveness was assessed using the standardized response mean (SRM). RESULTS: After 1 year, we observed a reduction in TMV (6.8%, p < 0.001) and muscle T2 (6.7%, p = 0.035), an increase in IMAT (9.7%, p < 0.001), FF (11.2%, p = 0.030), connective tissue (22%, p = 0.995), and T2*-IMAT (24%, p < 0.001), and alteration in muscle macromolecular structure (ΔMTR = -26%, p = 0.002). A decrease in muscle T2 correlated with an increase in T2*-IMAT (r = -0.47, p = 0.008). Deposition of connective tissue and IMAT correlated with deterioration in sIFA (r = 0.38, p = 0.032; r = 0.34, p = 0.048; respectively), whereas a decrease in TMV correlated with a decrease in quantitative muscle testing (r = 0.36, p = 0.035). The most responsive qMRI measures were T2*-IMAT (SRM = 1.50), TMV (SRM = -1.23), IMAT (SRM = 1.20), MTR (SRM = -0.83), and T2 relaxation time (SRM = -0.65). DISCUSSION: Progressive deterioration in muscle quality measured by qMRI is associated with a decline in physical performance. Inflammation may play a role in triggering fat infiltration into muscle. qMRI provides valid and responsive measures that might prove valuable in sIBM experimental trials and assessment of muscle pathologic processes. CLASSIFICATION OF EVIDENCE: This study provides Class I evidence that qMRI outcome measures are associated with physical performance measures in patients with sIBM.
Bridge CP, Gorman C, Pieper S, Doyle SW, Lennerz JK, Kalpathy-Cramer J, Clunie DA, Fedorov AY, Herrmann MD. Highdicom: a Python Library for Standardized Encoding of Image Annotations and Machine Learning Model Outputs in Pathology and Radiology. J Digit Imaging. 2022.Abstract
Machine learning (ML) is revolutionizing image-based diagnostics in pathology and radiology. ML models have shown promising results in research settings, but the lack of interoperability between ML systems and enterprise medical imaging systems has been a major barrier for clinical integration and evaluation. The DICOM® standard specifies information object definitions (IODs) and services for the representation and communication of digital images and related information, including image-derived annotations and analysis results. However, the complexity of the standard represents an obstacle for its adoption in the ML community and creates a need for software libraries and tools that simplify working with datasets in DICOM format. Here we present the highdicom library, which provides a high-level application programming interface (API) for the Python programming language that abstracts low-level details of the standard and enables encoding and decoding of image-derived information in DICOM format in a few lines of Python code. The highdicom library leverages NumPy arrays for efficient data representation and ties into the extensive Python ecosystem for image processing and machine learning. Simultaneously, by simplifying creation and parsing of DICOM-compliant files, highdicom achieves interoperability with the medical imaging systems that hold the data used to train and run ML models, and ultimately communicate and store model outputs for clinical use. We demonstrate through experiments with slide microscopy and computed tomography imaging, that, by bridging these two ecosystems, highdicom enables developers and researchers to train and evaluate state-of-the-art ML models in pathology and radiology while remaining compliant with the DICOM standard and interoperable with clinical systems at all stages. To promote standardization of ML research and streamline the ML model development and deployment process, we made the library available free and open-source at https://github.com/herrmannlab/highdicom .
Pace DF, Dalca AV, Brosch T, Geva T, Powell AJ, Weese J, Moghari MH, Golland P. Learned Iterative Segmentation of Highly Variable Anatomy From Limited Data: Applications to Whole Heart Segmentation for Congenital Heart Disease. Med Image Anal. 2022;80 :102469.Abstract
Training deep learning models that segment an image in one step typically requires a large collection of manually annotated images that captures the anatomical variability in a cohort. This poses challenges when anatomical variability is extreme but training data is limited, as when segmenting cardiac structures in patients with congenital heart disease (CHD). In this paper, we propose an iterative segmentation model and show that it can be accurately learned from a small dataset. Implemented as a recurrent neural network, the model evolves a segmentation over multiple steps, from a single user click until reaching an automatically determined stopping point. We develop a novel loss function that evaluates the entire sequence of output segmentations, and use it to learn model parameters. Segmentations evolve predictably according to growth dynamics encapsulated by training data, which consists of images, partially completed segmentations, and the recommended next step. The user can easily refine the final segmentation by examining those that are earlier or later in the output sequence. Using a dataset of 3D cardiac MR scans from patients with a wide range of CHD types, we show that our iterative model offers better generalization to patients with the most severe heart malformations.
Zhang F, Wells WM, O'Donnell LJ. Deep Diffusion MRI Registration (DDMReg): A Deep Learning Method for Diffusion MRI Registration. IEEE Trans Med Imaging. 2022;41 (6) :1454-67.Abstract
In this paper, we present a deep learning method, DDMReg, for accurate registration between diffusion MRI (dMRI) datasets. In dMRI registration, the goal is to spatially align brain anatomical structures while ensuring that local fiber orientations remain consistent with the underlying white matter fiber tract anatomy. DDMReg is a novel method that uses joint whole-brain and tract-specific information for dMRI registration. Based on the successful VoxelMorph framework for image registration, we propose a novel registration architecture that leverages not only whole brain information but also tract-specific fiber orientation information. DDMReg is an unsupervised method for deformable registration between pairs of dMRI datasets: it does not require nonlinearly pre-registered training data or the corresponding deformation fields as ground truth. We perform comparisons with four state-of-the-art registration methods on multiple independently acquired datasets from different populations (including teenagers, young and elderly adults) and different imaging protocols and scanners. We evaluate the registration performance by assessing the ability to align anatomically corresponding brain structures and ensure fiber spatial agreement between different subjects after registration. Experimental results show that DDMReg obtains significantly improved registration performance compared to the state-of-the-art methods. Importantly, we demonstrate successful generalization of DDMReg to dMRI data from different populations with varying ages and acquired using different acquisition protocols and different scanners.
Berger M, Pirpamer L, Hofer E, Ropele S, Duering M, Gesierich B, Pasternak O, Enzinger C, Schmidt R, Koini M. Free Water Diffusion MRI and Executive Function With a Speed Component in Healthy Aging. Neuroimage. 2022;257 :119303.Abstract
Extracellular free water (FW) increases are suggested to better provide pathophysiological information in brain aging than conventional biomarkers such as fractional anisotropy. The aim of the present study was to determine the relationship between conventional biomarkers, FW in white matter hyperintensities (WMH), FW in normal appearing white matter (NAWM) and in white matter tracts and executive functions (EF) with a speed component in elderly persons. We examined 226 healthy elderly participants (median age 69.83 years, IQR: 56.99-74.42) who underwent brain MRI and neuropsychological examination. FW in WMH and in NAWM as well as FW corrected diffusion metrics and measures derived from conventional MRI (white matter hyperintensities, brain volume, lacunes) were used in partial correlation (adjusted for age) to assess their correlation with EF with a speed component. Random forest analysis was used to assess the relative importance of these variables as determinants. Lastly, linear regression analyses of FW in white matter tracts corrected for risk factors of cognitive and white matter deterioration, were used to examine the role of specific tracts on EF with a speed component, which were then ranked with random forest regression. Partial correlation analyses revealed that almost all imaging metrics showed a significant association with EF with a speed component (r = -0.213 - 0.266). Random forest regression highlighted FW in WMH and in NAWM as most important among all diffusion and structural MRI metrics. The fornix (R2=0.421, p = 0.018) and the corpus callosum (genu (R2 = 0.418, p = 0.021), prefrontal (R2 = 0.416, p = 0.026), premotor (R2 = 0.418, p = 0.021)) were associated with EF with a speed component in tract based regression analyses and had highest variables importance. In a normal aging population FW in WMH and NAWM is more closely related to EF with a speed component than standard DTI and brain structural measures. Higher amounts of FW in the fornix and the frontal part of the corpus callosum leads to deteriorating EF with a speed component.
Wang S, Zhang F, Huang P, Hong H, Jiaerken Y, Yu X, Zhang R, Zeng Q, Zhang Y, Kikinis R, et al. Superficial White Matter Microstructure Affects Processing Speed in Cerebral Small Vessel Disease. Hum Brain Mapp. 2022.Abstract
White matter hyperintensities (WMH) are a typical feature of cerebral small vessel disease (CSVD), which contributes to about 50% of dementias worldwide. Microstructural alterations in deep white matter (DWM) have been widely examined in CSVD. However, little is known about abnormalities in superficial white matter (SWM) and their relevance for processing speed, the main cognitive deficit in CSVD. In 141 CSVD patients, processing speed was assessed using Trail Making Test Part A. White matter abnormalities were assessed by WMH burden (volume on T2-FLAIR) and diffusion MRI measures. SWM imaging measures had a large contribution to processing speed, despite a relatively low SWM WMH burden. Across all imaging measures, SWM free water (FW) had the strongest association with processing speed, followed by SWM mean diffusivity (MD). SWM FW was the only marker to significantly increase between two subgroups with the lowest WMH burdens. When comparing two subgroups with the highest WMH burdens, the involvement of WMH in the SWM was accompanied by significant differences in processing speed and white matter microstructure. Mediation analysis revealed that SWM FW fully mediated the association between WMH volume and processing speed, while no mediation effect of MD or DWM FW was observed. Overall, results suggest that the SWM has an important contribution to processing speed, while SWM FW is a sensitive imaging marker associated with cognition in CSVD. This study extends the current understanding of CSVD-related dysfunction and suggests that the SWM, as an understudied region, can be a potential target for monitoring pathophysiological processes.
Pujol S, Cabeen RP, Yelnik J, François C, Fernandez Vidal S, Karachi C, Bardinet E, Cosgrove RG, Kikinis R. Somatotopic Organization of Hyperdirect Pathway Projections From the Primary Motor Cortex in the Human Brain. Front Neurol. 2022;13 :791092.Abstract
Background: The subthalamic nucleus (STN) is an effective neurosurgical target to improve motor symptoms in Parkinson's Disease (PD) patients. MR-guided Focused Ultrasound (MRgFUS) subthalamotomy is being explored as a therapeutic alternative to Deep Brain Stimulation (DBS) of the STN. The hyperdirect pathway provides a direct connection between the cortex and the STN and is likely to play a key role in the therapeutic effects of MRgFUS intervention in PD patients. Objective: This study aims to investigate the topography and somatotopy of hyperdirect pathway projections from the primary motor cortex (M1). Methods: We used advanced multi-fiber tractography and high-resolution diffusion MRI data acquired on five subjects of the Human Connectome Project (HCP) to reconstruct hyperdirect pathway projections from M1. Two neuroanatomy experts reviewed the anatomical accuracy of the tracts. We extracted the fascicles arising from the trunk, arm, hand, face and tongue area from the reconstructed pathways. We assessed the variability among subjects based on the fractional anisotropy (FA) and mean diffusivity (MD) of the fibers. We evaluated the spatial arrangement of the different fascicles using the Dice Similarity Coefficient (DSC) of spatial overlap and the centroids of the bundles. Results: We successfully reconstructed hyperdirect pathway projections from M1 in all five subjects. The tracts were in agreement with the expected anatomy. We identified hyperdirect pathway fascicles projecting from the trunk, arm, hand, face and tongue area in all subjects. Tract-derived measurements showed low variability among subjects, and similar distributions of FA and MD values among the fascicles projecting from different M1 areas. We found an anterolateral somatotopic arrangement of the fascicles in the corona radiata, and an average overlap of 0.63 in the internal capsule and 0.65 in the zona incerta. Conclusion: Multi-fiber tractography combined with high-resolution diffusion MRI data enables the identification of the somatotopic organization of the hyperdirect pathway. Our preliminary results suggest that the subdivisions of the hyperdirect pathway projecting from the trunk, arm, hand, face, and tongue motor area are intermixed at the level of the zona incerta and posterior limb of the internal capsule, with a predominantly overlapping topographical organization in both regions. Subject-specific knowledge of the hyperdirect pathway somatotopy could help optimize target definition in MRgFUS intervention.
Chauvin L, Kumar K, Desrosiers C, Wells W, Toews M. Efficient Pairwise Neuroimage Analysis Using the Soft Jaccard Index and 3D Keypoint Sets. IEEE Trans Med Imaging. 2022;41 (4) :836-45.Abstract
We propose a novel pairwise distance measure between image keypoint sets, for the purpose of large-scale medical image indexing. Our measure generalizes the Jaccard index to account for soft set equivalence (SSE) between keypoint elements, via an adaptive kernel framework modeling uncertainty in keypoint appearance and geometry. A new kernel is proposed to quantify the variability of keypoint geometry in location and scale. Our distance measure may be estimated between O (N 2) image pairs in [Formula: see text] operations via keypoint indexing. Experiments report the first results for the task of predicting family relationships from medical images, using 1010 T1-weighted MRI brain volumes of 434 families including monozygotic and dizygotic twins, siblings and half-siblings sharing 100%-25% of their polymorphic genes. Soft set equivalence and the keypoint geometry kernel improve upon standard hard set equivalence (HSE) and appearance kernels alone in predicting family relationships. Monozygotic twin identification is near 100%, and three subjects with uncertain genotyping are automatically paired with their self-reported families, the first reported practical application of image-based family identification. Our distance measure can also be used to predict group categories, sex is predicted with an AUC = 0.97. Software is provided for efficient fine-grained curation of large, generic image datasets.
Abulnaga MS, Abaci Turk E, Bessmeltsev M, Grant EP, Solomon J, Golland P. Volumetric Parameterization of the Placenta to a Flattened Template. IEEE Trans Med Imaging. 2022;41 (4) :925-36.Abstract
We present a volumetric mesh-based algorithm for parameterizing the placenta to a flattened template to enable effective visualization of local anatomy and function. MRI shows potential as a research tool as it provides signals directly related to placental function. However, due to the curved and highly variable in vivo shape of the placenta, interpreting and visualizing these images is difficult. We address interpretation challenges by mapping the placenta so that it resembles the familiar ex vivo shape. We formulate the parameterization as an optimization problem for mapping the placental shape represented by a volumetric mesh to a flattened template. We employ the symmetric Dirichlet energy to control local distortion throughout the volume. Local injectivity in the mapping is enforced by a constrained line search during the gradient descent optimization. We validate our method using a research study of 111 placental shapes extracted from BOLD MRI images. Our mapping achieves sub-voxel accuracy in matching the template while maintaining low distortion throughout the volume. We demonstrate how the resulting flattening of the placenta improves visualization of anatomy and function. Our code is freely available at https://github.com/mabulnaga/placenta-flattening.
Zheng J, Yang Q, Makris N, Huang K, Liang J, Ye C, Yu X, Tian M, Ma T, Mou T, et al. Three-Dimensional Digital Reconstruction of the Cerebellar Cortex: Lobule Thickness, Surface Area Measurements, and Layer Architecture. Cerebellum. 2022.Abstract
The cerebellum is ontogenetically one of the first structures to develop in the central nervous system; nevertheless, it has been only recently reconsidered for its significant neurobiological, functional, and clinical relevance in humans. Thus, it has been a relatively under-studied compared to the cerebrum. Currently, non-invasive imaging modalities can barely reach the necessary resolution to unfold its entire, convoluted surface, while only histological analyses can reveal local information at the micrometer scale. Herein, we used the BigBrain dataset to generate area and point-wise thickness measurements for all layers of the cerebellar cortex and for each lobule in particular. We found that the overall surface area of the cerebellar granular layer (including Purkinje cells) was 1,732 cm2 and the molecular layer was 1,945 cm2. The average thickness of the granular layer is 0.88 mm (± 0.83) and that of the molecular layer is 0.32 mm (± 0.08). The cerebellum (both granular and molecular layers) is thicker at the depth of the sulci and thinner at the crowns of the gyri. Globally, the granular layer is thicker in the lateral-posterior-inferior region than the medial-superior regions. The characterization of individual layers in the cerebellum achieved herein represents a stepping-stone for investigations interrelating structural and functional connectivity with cerebellar architectonics using neuroimaging, which is a matter of considerable relevance in basic and clinical neuroscience. Furthermore, these data provide templates for the construction of cerebellar topographic maps and the precise localization of structural and functional alterations in diseases affecting the cerebellum.
Seitz-Holland J, Seethaler M, Makris N, Rushmore J, Cho K-IK, Rizzoni E, Vangel M, Sahin OS, Heller C, Pasternak O, et al. The Association of Matrix Metalloproteinase 9 (MMP9) With Hippocampal Volume in Schizophrenia: A Preliminary MRI Study. Neuropsychopharmacology. 2022;47 (2) :524-30.Abstract
Matrix metalloproteinases 9 (MMP9) are enzymes involved in regulating neuroplasticity in the hippocampus. This, combined with evidence for disrupted hippocampal structure and function in schizophrenia, has prompted our current investigation into the relationship between MMP9 and hippocampal volumes in schizophrenia. 34 healthy individuals (mean age = 32.50, male = 21, female = 13) and 30 subjects with schizophrenia (mean age = 33.07, male = 19, female = 11) underwent a blood draw and T1-weighted magnetic resonance imaging. The hippocampus was automatically segmented utilizing FreeSurfer. MMP9 plasma levels were measured with ELISA. ANCOVAs were conducted to compare MMP9 plasma levels (corrected for age and sex) and hippocampal volumes between groups (corrected for age, sex, total intracranial volume). Spearman correlations were utilized to investigate the relationship between symptoms, medication, duration of illness, number of episodes, and MMP9 plasma levels in patients. Last, we explored the correlation between MMP9 levels and hippocampal volumes in patients and healthy individuals separately. Patients displayed higher MMP9 plasma levels than healthy individuals (F(1, 60) = 21.19, p < 0.0001). MMP9 levels correlated with negative symptoms in patients (R = 0.39, p = 0.035), but not with medication, duration of illness, or the number of episodes. Further, patients had smaller left (F(1,59) = 9.12, p = 0.0040) and right (F(1,59) = 6.49, p = 0.013) hippocampal volumes. Finally, left (R = -0.39, p = 0.034) and right (R = -0.37, p = 0.046) hippocampal volumes correlated negatively with MMP9 plasma levels in patients. We observe higher MMP9 plasma levels in SCZ, associated with lower hippocampal volumes, suggesting involvement of MMP9 in the pathology of SCZ. Future studies are needed to investigate how MMP9 influences the pathology of SCZ over the lifespan, whether the observed associations are specific for schizophrenia, and if a therapeutic modulation of MMP9 promotes neuroprotective effects in SCZ.
Bayat A, Pace DF, Sekuboyina A, Payer C, Stern D, Urschler M, Kirschke JS, Menze BH. Anatomy-Aware Inference of the 3D Standing Spine Posture from 2D Radiographs. Tomography. 2022;8 (1) :479-96.Abstract
An important factor for the development of spinal degeneration, pain and the outcome of spinal surgery is known to be the balance of the spine. It must be analyzed in an upright, standing position to ensure physiological loading conditions and visualize load-dependent deformations. Despite the complex 3D shape of the spine, this analysis is currently performed using 2D radiographs, as all frequently used 3D imaging techniques require the patient to be scanned in a prone position. To overcome this limitation, we propose a deep neural network to reconstruct the 3D spinal pose in an upright standing position, loaded naturally. Specifically, we propose a novel neural network architecture, which takes orthogonal 2D radiographs and infers the spine's 3D posture using vertebral shape priors. In this work, we define vertebral shape priors using an atlas and a spine shape prior, incorporating both into our proposed network architecture. We validate our architecture on digitally reconstructed radiographs, achieving a 3D reconstruction Dice of 0.95, indicating an almost perfect 2D-to-3D domain translation. Validating the reconstruction accuracy of a 3D standing spine on real data is infeasible due to the lack of a valid ground truth. Hence, we design a novel experiment for this purpose, using an orientation invariant distance metric, to evaluate our model's ability to synthesize full-3D, upright, and patient-specific spine models. We compare the synthesized spine shapes from clinical upright standing radiographs to the same patient's 3D spinal posture in the prone position from CT.
O'Donnell LJ. Editorial for "Early-Onset Micromorphological Changes of Neuronal Fiber Bundles During Radiotherapy". J Magn Reson Imaging. 2022;56 (1) :219-220.
Zhang F, Daducci A, He Y, Schiavi S, Seguin C, Smith RE, Yeh C-H, Zhao T, O'Donnell LJ. Quantitative Mapping of the Brain's Structural Connectivity Using Diffusion MRI Tractography: A Review. Neuroimage. 2022;249 :118870.Abstract
Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.
Zaks N, Velikonja T, Parvaz MA, Zinberg J, Done M, Mathalon DH, Addington J, Cadenhead K, Cannon T, Cornblatt B, et al. Sleep Disturbance in Individuals at Clinical High Risk for Psychosis. Schizophr Bull. 2022;48 (1) :111-21.Abstract
INTRODUCTION: Disturbed sleep is a common feature of psychotic disorders that is also present in the clinical high risk (CHR) state. Evidence suggests a potential role of sleep disturbance in symptom progression, yet the interrelationship between sleep and CHR symptoms remains to be determined. To address this knowledge gap, we examined the association between disturbed sleep and CHR symptoms over time. METHODS: Data were obtained from the North American Prodrome Longitudinal Study (NAPLS)-3 consortium, including 688 CHR individuals and 94 controls (mean age 18.25, 46% female) for whom sleep was tracked prospectively for 8 months. We used Cox regression analyses to investigate whether sleep disturbances predicted conversion to psychosis up to >2 years later. With regressions and cross-lagged panel models, we analyzed longitudinal and bidirectional associations between sleep (the Pittsburgh Sleep Quality Index in conjunction with additional sleep items) and CHR symptoms. We also investigated the independent contribution of individual sleep characteristics on CHR symptom domains separately and explored whether cognitive impairments, stress, depression, and psychotropic medication affected the associations. RESULTS: Disturbed sleep at baseline did not predict conversion to psychosis. However, sleep disturbance was strongly correlated with heightened CHR symptoms over time. Depression accounted for half of the association between sleep and symptoms. Importantly, sleep was a significant predictor of CHR symptoms but not vice versa, although bidirectional effect sizes were similar. DISCUSSION: The critical role of sleep disturbance in CHR symptom changes suggests that sleep may be a promising intervention target to moderate outcome in the CHR state.
Ji Y, Hoge WS, Gagoski B, Westin C-F, Rathi Y, Ning L. Accelerating Joint Relaxation-Diffusion MRI by Integrating Time Division Multiplexing and Simultaneous Multi-Slice (TDM-SMS) Strategies. Magn Reson Med. 2022;87 (6) :2697-709.Abstract
PURPOSE: To accelerate the acquisition of relaxation-diffusion imaging by integrating time-division multiplexing (TDM) with simultaneous multi-slice (SMS) for EPI and evaluate imaging quality and diffusion measures. METHODS: The time-division multiplexing (TDM) technique and SMS method were integrated to achieve a high slice-acceleration (e.g., 6×) factor for acquiring relaxation-diffusion MRI. Two variants of the sequence, referred to as TDM3e-SMS and TDM2s-SMS, were developed to simultaneously acquire slice groups with three distinct TEs and two slice groups with the same TE, respectively. Both sequences were evaluated on a 3T scanner with in vivo human brains and compared with standard single-band (SB) -EPI and SMS-EPI using diffusion measures and tractography results. RESULTS: Experimental results showed that the TDM3e-SMS sequence with total slice acceleration of 6 (multiplexing factor (MP) = 3 × multi-band factor (MB) = 2) provided similar image intensity and microstructure measures compared to standard SMS-EPI with MB = 2, and yielded less bias in intensity compared to standard SMS-EPI with MB = 4. The three sequences showed a similar positive correlation between TE and mean kurtosis (MK) and a negative correlation between TE and mean diffusivity (MD) in white matter. Multi-fiber tractography also shows consistency of results in TE-dependent measures between different sequences. The TDM2s-SMS sequence (MP = 2, MB = 2) also provided imaging measures similar to standard SMS-EPI sequences (MB = 2) for single-TE diffusion imaging. CONCLUSIONS: The TDM-SMS sequence can provide additional 2× to 3× acceleration to SMS without degrading imaging quality. With the significant reduction in scan time, TDM-SMS makes joint relaxation-diffusion MRI a feasible technique in neuroimaging research to investigate new markers of brain disorders.
McGarry SD, Brehler M, Bukowy JD, Lowman AK, Bobholz SA, Duenweg SR, Banerjee A, Hurrell SL, Malyarenko D, Chenevert TL, et al. Multi-Site Concordance of Diffusion-Weighted Imaging Quantification for Assessing Prostate Cancer Aggressiveness. J Magn Reson Imaging. 2022;55 (6) :1745-58.Abstract
BACKGROUND: Diffusion-weighted imaging (DWI) is commonly used to detect prostate cancer, and a major clinical challenge is differentiating aggressive from indolent disease. PURPOSE: To compare 14 site-specific parametric fitting implementations applied to the same dataset of whole-mount pathologically validated DWI to test the hypothesis that cancer differentiation varies with different fitting algorithms. STUDY TYPE: Prospective. POPULATION: Thirty-three patients prospectively imaged prior to prostatectomy. FIELD STRENGTH/SEQUENCE: 3 T, field-of-view optimized and constrained undistorted single-shot DWI sequence. ASSESSMENT: Datasets, including a noise-free digital reference object (DRO), were distributed to the 14 teams, where locally implemented DWI parameter maps were calculated, including mono-exponential apparent diffusion coefficient (MEADC), kurtosis (K), diffusion kurtosis (DK), bi-exponential diffusion (BID), pseudo-diffusion (BID*), and perfusion fraction (F). The resulting parametric maps were centrally analyzed, where differentiation of benign from cancerous tissue was compared between DWI parameters and the fitting algorithms with a receiver operating characteristic area under the curve (ROC AUC). STATISTICAL TEST: Levene's test, P < 0.05 corrected for multiple comparisons was considered statistically significant. RESULTS: The DRO results indicated minimal discordance between sites. Comparison across sites indicated that K, DK, and MEADC had significantly higher prostate cancer detection capability (AUC range = 0.72-0.76, 0.76-0.81, and 0.76-0.80 respectively) as compared to bi-exponential parameters (BID, BID*, F) which had lower AUC and greater between site variation (AUC range = 0.53-0.80, 0.51-0.81, and 0.52-0.80 respectively). Post-processing parameters also affected the resulting AUC, moving from, for example, 0.75 to 0.87 for MEADC varying cluster size. DATA CONCLUSION: We found that conventional diffusion models had consistent performance at differentiating prostate cancer from benign tissue. Our results also indicated that post-processing decisions on DWI data can affect sensitivity and specificity when applied to radiological-pathological studies in prostate cancer. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 3.
Brabec J, Szczepankiewicz F, Lennartsson F, Englund E, Pebdani H, Bengzon J, Knutsson L, Westin C-F, Sundgren PC, Nilsson M. Histogram Analysis of Tensor-Valued Diffusion MRI in Meningiomas: Relation to Consistency, Histological Grade and Type. Neuroimage Clin. 2022;33 :102912.Abstract
BACKGROUND: Preoperative radiological assessment of meningioma characteristics is of value for pre- and post-operative patient management, counselling, and surgical approach. PURPOSE: To investigate whether tensor-valued diffusion MRI can add to the preoperative prediction of meningioma consistency, grade and type. MATERIALS AND METHODS: 30 patients with intracranial meningiomas (22 WHO grade I, 8 WHO grade II) underwent MRI prior to surgery. Diffusion MRI was performed with linear and spherical b-tensors with b-values up to 2000 s/mm2. The data were used to estimate mean diffusivity (MD), fractional anisotropy (FA), mean kurtosis (MK) and its components-the anisotropic and isotropic kurtoses (MKA and MKI). Meningioma consistency was estimated for 16 patients during resection based on ultrasonic aspiration intensity, ease of resection with instrumentation or suction. Grade and type were determined by histopathological analysis. The relation between consistency, grade and type and dMRI parameters was analyzed inside the tumor ("whole-tumor") and within brain tissue in the immediate periphery outside the tumor ("rim") by histogram analysis. RESULTS: Lower 10th percentiles of MK and MKA in the whole-tumor were associated with firm consistency compared with pooled soft and variable consistency (n = 7 vs 9; U test, p = 0.02 for MKA 10 and p = 0.04 for MK10) and lower 10th percentile of MD with variable against soft and firm (n = 5 vs 11; U test, p = 0.02). Higher standard deviation of MKI in the rim was associated with lower grade (n = 22 vs 8; U test, p = 0.04) and in the MKI maps we observed elevated rim-like structure that could be associated with grade. Higher median MKA and lower median MKI distinguished psammomatous type from other pooled meningioma types (n = 5 vs 25; U test; p = 0.03 for MKA 50 and p = 0.03 and p = 0.04 for MKI 50). CONCLUSION: Parameters from tensor-valued dMRI can facilitate prediction of consistency, grade and type.