Publications by Year: 2020

BACKGROUND: Post-traumatic stress disorder (PTSD) is a psychiatric disorder that afflicts many individuals, yet the neuropathological mechanisms that contribute to this disorder remain to be fully determined. Moreover, it is unclear how exposure to mild traumatic brain injury (mTBI), a condition that is often comorbid with PTSD, particularly among military personnel, affects the clinical and neurological presentation of PTSD. To address these issues, the present study explores relationships between PTSD symptom severity and the microstructure of limbic and paralimbic gray matter brain regions, as well as the impact of mTBI comorbidity on these relationships. METHODS: Structural and diffusion MRI data were acquired from 102 male veterans who were diagnosed with current PTSD. Diffusion data were analyzed with free-water imaging to quantify average CSF-corrected fractional anisotropy (FA) and mean diffusivity (MD) in 18 limbic and paralimbic gray matter regions. Associations between PTSD symptom severity and regional average dMRI measures were examined with repeated measures linear mixed models. Associations were studied separately in veterans with PTSD only, and in veterans with PTSD and a history of military mTBI. RESULTS: Analyses revealed that in the PTSD only cohort, more severe symptoms were associated with higher FA in the right amygdala-hippocampus complex, lower FA in the right cingulate cortex, and lower MD in the left medial orbitofrontal cortex. In the PTSD and mTBI cohort, more severe PTSD symptoms were associated with higher FA bilaterally in the amygdala-hippocampus complex, with higher FA bilaterally in the nucleus accumbens, with lower FA bilaterally in the cingulate cortex, and with higher MD in the right amygdala-hippocampus complex. CONCLUSIONS: These findings suggest that the microstructure of limbic and paralimbic brain regions may influence PTSD symptomatology. Further, given the additional associations observed between microstructure and symptom severity in veterans with head trauma, we speculate that mTBI may exacerbate the impact of brain microstructure on PTSD symptoms, especially within regions of the brain known to be vulnerable to chronic stress. A heightened sensitivity to the microstructural environment of the brain could partially explain why individuals with PTSD and mTBI comorbidity experience more severe symptoms and poorer illness prognoses than those without a history of brain injury. The relevance of these microstructural findings to the conceptualization of PTSD as being a disorder of stress-induced neuronal connectivity loss is discussed.

We investigated brain wiring in chronic schizophrenia and healthy controls in frontostriatal circuits using diffusion magnetic resonance imaging tractography in a novel way. We extracted diffusion streamlines in 27 chronic schizophrenia and 26 healthy controls connecting 4 frontal subregions to the striatum. We labeled the projection zone striatal surface voxels into 2 subtypes: dominant-input from a single cortical subregion, and, functionally integrative, with mixed-input from diverse cortical subregions. We showed: 1) a group difference for total striatal surface voxel number (P = .045) driven by fewer mixed-input voxels in the left (P = .007), but not right, hemisphere; 2) a group by hemisphere interaction for the ratio quotient between voxel subtypes (P = .04) with a left (P = .006), but not right, hemisphere increase in schizophrenia, also reflecting fewer mixed-input voxels; and 3) fewer mixed-input voxel counts in schizophrenia (P = .045) driven by differences in left hemisphere limbic (P = .007) and associative (P = .01), but not sensorimotor, striatum. These results demonstrate a less integrative pattern of frontostriatal structural connectivity in chronic schizophrenia. A diminished integrative pattern yields a less complex input pattern to the striatum from the cortex with less circuit integration at the level of the striatum. Further, as brain wiring occurs during early development, aberrant brain wiring could serve as a developmental biomarker for schizophrenia.

Joint relaxation-diffusion measurements can provide new insight about the tissue microstructural properties. Most recent methods have focused on inverting the Laplace transform to recover the joint distribution of relaxation-diffusion. However, as is well-known, this problem is notoriously ill-posed and numerically unstable. In this work, we address this issue by directly computing the joint moments of transverse relaxation rate and diffusivity, which can be robustly estimated. To zoom into different parts of the joint distribution, we further enhance our method by applying multiplicative filters to the joint probability density function of relaxation and diffusion and compute the corresponding moments. We propose an approach to use these moments to compute several novel scalar indices to characterize specific properties of the underlying tissue microstructure. Furthermore, for the first time, we propose an algorithm to estimate diffusion signals that are independent of echo time based on the moments of the marginal probability density function of diffusion. We demonstrate its utility in extracting tissue information not contaminated with multiple intra-voxel relaxation rates. We compare the performance of four types of filters that zoom into tissue components with different relaxation and diffusion properties and demonstrate it on an in-vivo human dataset. Experimental results show that these filters are able to characterize heterogeneous tissue microstructure. Moreover, the filtered diffusion signals are also able to distinguish fiber bundles with similar orientations but different relaxation rates. The proposed method thus allows to characterize the neural microstructure information in a robust and unique manner not possible using existing techniques.

COPD is a heterogeneous syndrome. Many COPD subtypes have been proposed, but there is not yet consensus on how many COPD subtypes there are and how they should be defined. The COPD Genetic Epidemiology Study (COPDGene), which has generated 10-year longitudinal chest imaging, spirometry, and molecular data, is a rich resource for relating COPD phenotypes to underlying genetic and molecular mechanisms. In this article, we place COPDGene clustering studies in context with other highly cited COPD clustering studies, and summarize the main COPD subtype findings from COPDGene. First, most manifestations of COPD occur along a continuum, which explains why continuous aspects of COPD or disease axes may be more accurate and reproducible than subtypes identified through clustering methods. Second, continuous COPD-related measures can be used to create subgroups through the use of predictive models to define cut-points, and we review COPDGene research on blood eosinophil count thresholds as a specific example. Third, COPD phenotypes identified or prioritized through machine learning methods have led to novel biological discoveries, including novel emphysema genetic risk variants and systemic inflammatory subtypes of COPD. Fourth, trajectory-based COPD subtyping captures differences in the longitudinal evolution of COPD, addressing a major limitation of clustering analyses that are confounded by disease severity. Ongoing longitudinal characterization of subjects in COPDGene will provide useful insights about the relationship between lung imaging parameters, molecular markers, and COPD progression that will enable the identification of subtypes based on underlying disease processes and distinct patterns of disease progression, with the potential to improve the clinical relevance and reproducibility of COPD subtypes.

The middle longitudinal fascicle (MdLF) is a recently delineated association cortico-cortical fiber pathway in humans, connecting superior temporal gyrus and temporal pole principally with the angular gyrus, and is likely to be involved in language processing. However, the MdLF has not been studied in language disorders as primary progressive aphasia (PPA). We hypothesized that the MdLF will exhibit evidence of neurodegeneration in PPA patients. In this study, 20 PPA patients and 25 healthy controls were recruited in the Primary Progressive Aphasia program in the Massachusetts General Hospital Frontotemporal Disorders Unit. We used diffusion tensor imaging (DTI) tractography to reconstruct the MdLF and extract tract-specific DTI metrics (fractional anisotropy (FA), radial diffusivity (RD), mean diffusivity (MD) and axial diffusivity (AD)) to assess white matter changes in PPA and their relationship with language impairments. We found severe WM damage in the MdLF in PPA patients, which was principally pronounced in the left hemisphere. Moreover, the WM alterations in the MdLF in the dominant hemisphere were significantly correlated with impairments in word comprehension and naming, but not with articulation and fluency. In addition, asymmetry analysis revealed that the DTI metrics of controls were similar for each hemisphere, whereas PPA patients had clear laterality differences in MD, AD and RD. These findings add new insight into the localization and severity of white matter fiber bundle neurodegeneration in PPA, and provide evidence that degeneration of the MdLF contribute to impairment in semantic processing and lexical retrieval in PPA.

Cerebrospinal fluid partial volume effect is a known bias in the estimation of Diffusion Tensor Imaging (DTI) parameters from diffusion MRI data. The Free-Water Imaging model for diffusion MRI data adds a second compartment to the DTI model, which explicitly accounts for the signal contribution of extracellular free-water, such as cerebrospinal fluid. As a result the DTI parameters obtained through the free-water model are corrected for partial volume effects, and thus better represent tissue microstructure. In addition, the model estimates the fractional volume of free-water, and can be used to monitor changes in the extracellular space. Under certain assumptions, the model can be estimated from single-shell diffusion MRI data. However, by using data from multi-shell diffusion acquisitions, these assumptions can be relaxed, and the fit becomes more robust. Nevertheless, fitting the model to multi-shell data requires high computational cost, with a non-linear iterative minimization, which has to be initialized close enough to the global minimum to avoid local minima and to robustly estimate the model parameters. Here we investigate the properties of the main initialization approaches that are currently being used, and suggest new fast approaches to improve the initial estimates of the model parameters. We show that our proposed approaches provide a fast and accurate initial approximation of the model parameters, which is very close to the final solution. We demonstrate that the proposed initializations improve the final outcome of non-linear model fitting.

PURPOSE: To develop an accelerated, robust, and accurate diffusion MRI acquisition and reconstruction technique for submillimeter whole human brain in vivo scan on a clinical scanner. METHODS: We extend the ultra-high resolution diffusion MRI acquisition technique, gSlider, by allowing undersampling in q-space and radiofrequency (RF)-encoding space, thereby dramatically reducing the total acquisition time of conventional gSlider. The novel method, termed gSlider-SR, compensates for the lack of acquired information by exploiting redundancy in the dMRI data using a basis of spherical ridgelets (SR), while simultaneously enhancing the signal-to-noise ratio. Using Monte Carlo simulation with realistic noise levels and several acquisitions of in vivo human brain dMRI data (acquired on a Siemens Prisma 3T scanner), we demonstrate the efficacy of our method using several quantitative metrics. RESULTS: For high-resolution dMRI data with realistic noise levels (synthetically added), we show that gSlider-SR can reconstruct high-quality dMRI data at different acceleration factors preserving both signal and angular information. With in vivo data, we demonstrate that gSlider-SR can accurately reconstruct 860 μm diffusion MRI data (64 diffusion directions at ), at comparable quality as that obtained with conventional gSlider with four averages, thereby providing an eight-fold reduction in scan time (from 1 hour 20 to 10 minutes). CONCLUSIONS: gSlider-SR enables whole-brain high angular resolution dMRI at a submillimeter spatial resolution with a dramatically reduced acquisition time, making it feasible to use the proposed scheme on existing clinical scanners.

Age-related memory impairments have been linked to differences in structural brain parameters, including cerebral white matter (WM) microstructure and hippocampal (HC) volume, but their combined influences are rarely investigated. In a population-based sample of 337 older participants aged 61-82 years (Mage = 69.66, SDage = 3.92 years), we modeled the independent and joint effects of limbic WM microstructure and HC subfield volumes on verbal learning. Participants completed a verbal learning task of recall over five repeated trials and underwent magnetic resonance imaging (MRI), including structural and diffusion scans. We segmented three HC subregions on high-resolution MRI data and sampled mean fractional anisotropy (FA) from bilateral limbic WM tracts identified via deterministic fiber tractography. Using structural equation modeling, we evaluated the associations between learning rate and latent factors representing FA sampled from limbic WM tracts, and HC subfield volumes, and their latent interaction. Results showed limbic WM and the interaction of HC and WM-but not HC volume alone-predicted verbal learning rates. Model decomposition revealed HC volume is only positively associated with learning rate in individuals with higher WM anisotropy. We conclude that the structural characteristics of limbic WM regions and HC volume jointly contribute to verbal learning in older adults.

BACKGROUND: In previous work, we applied novel in vivo imaging methods to reveal that white matter pathology in patients with first-episode psychosis (FEP) is mainly characterized by excessive extracellular free-water, and to a lesser extent by cellular processes, such as demyelination. Here, we apply a back-translational approach to evaluate whether or not a rodent model of maternal immune activation (MIA) induces patterns of white matter pathology that we observed in patients with FEP. To this end, we examined free-water and tissue-specific white matter alterations in rats born to mothers exposed to the viral mimic polyriboinosinic-polyribocytidylic acid (Poly-I:C) in pregnancy, which is widely used to produce alterations relevant to schizophrenia and is characterized by a robust neuroinflammatory response. METHOD: Pregnant dams were injected on gestational day 15 with the viral mimic Poly-I:C (4 mg/kg) or saline. Diffusion-weighted magnetic resonance images were acquired from 17 male offspring (9 Poly-I:C and 8 saline) on postnatal day 90, after the emergence of brain structural and behavioral abnormalities. The free-water fraction (FW) and tissue-specific fractional anisotropy (FA), as well as conventional fractional anisotropy (FA) were computed across voxels traversing a white matter skeleton. Voxel-wise and whole-brain averaged white matter were tested for significant microstructural alterations in immune-challenged, relative to saline-exposed offspring.