Publications

Diffusion MRI may enable non-invasive mapping of axonal microstructure. Most approaches infer axon diameters from effects of time-dependent diffusion on the diffusion-weighted MR signal by modeling axons as straight cylinders. Axons do not, however, propagate in straight trajectories, and so far the impact of the axonal trajectory on diameter estimation has been insufficiently investigated. Here, we employ a toy model of axons, which we refer to as the undulating thin fiber model, to analyze the impact of undulating trajectories on the time dependence of diffusion. We study time-dependent diffusion in the frequency domain and characterize the diffusion spectrum by its height, width, and low-frequency behavior (power law exponent). Results show that microscopic orientation dispersion of the thin fibers is the main parameter that determines the characteristics of the diffusion spectra. At lower frequencies (longer diffusion times), straight cylinders and undulating thin fibers can have virtually identical spectra. If the straight-cylinder assumption is used to interpret data from undulating thin axons, the diameter is overestimated by an amount proportional to the undulation amplitude and microscopic orientation dispersion of the fibers. At higher frequencies (shorter diffusion times), spectra from cylinders and undulating thin fibers differ. The low-frequency behavior of the spectra from the undulating thin fibers may also differ from that of cylinders, because the power law exponent of undulating fibers can reach values below 2 for experimentally relevant frequency ranges. In conclusion, we argue that the non-straight nature of axonal trajectories should not be overlooked when analyzing and interpreting diffusion MRI data.

Motion is a major confound in diffusion-weighted imaging (DWI) in the body, and it is a common cause of image artefacts. The effects are particularly severe in cardiac applications, due to the nonrigid cyclical deformation of the myocardium. Spin echo-based DWI commonly employs gradient moment-nulling techniques to desensitise the acquisition to velocity and acceleration, ie, nulling gradient moments up to the 2nd order (M2-nulled). However, current M2-nulled DWI scans are limited to encode diffusion along a single direction at a time. We propose a method for designing b-tensors of arbitrary shapes, including planar, spherical, prolate and oblate tensors, while nulling gradient moments up to the 2nd order and beyond. The design strategy comprises initialising the diffusion encoding gradients in two encoding blocks about the refocusing pulse, followed by appropriate scaling and rotation, which further enables nulling undesired effects of concomitant gradients. Proof-of-concept assessment of in vivo mean diffusivity (MD) was performed using linear and spherical tensor encoding (LTE and STE, respectively) in the hearts of five healthy volunteers. The results of the M2-nulled STE showed that (a) the sequence was robust to cardiac motion, and (b) MD was higher than that acquired using standard M2-nulled LTE, where diffusion-weighting was applied in three orthogonal directions, which may be attributed to the presence of restricted diffusion and microscopic diffusion anisotropy. Provided adequate signal-to-noise ratio, STE could significantly shorten estimation of MD compared with the conventional LTE approach. Importantly, our theoretical analysis and the proposed gradient waveform design may be useful in microstructure imaging beyond diffusion tensor imaging where the effects of motion must be suppressed.

To date, few studies have evaluated the contribution of early life experiences to neurocognitive abnormalities observed in posttraumatic stress disorder (PTSD). Childhood maltreatment is common among individuals with PTSD and is thought to catalyze stress-related biobehavioral changes that might impact both brain structure and function in adulthood. The current study examined differences in brain morphology (brain volume, cortical thickness) and neuropsychological performance in individuals with PTSD characterized by low or high self-reported childhood maltreatment, compared with healthy comparison participants. Data were drawn from the INjury and TRaUmatic STress (INTRuST) Clinical Consortium imaging repository, which contains MRI and self-report data for individuals classified as PTSD positive (with and without a history of mild traumatic brain injury [mTBI]), individuals with mTBI only, and healthy comparison participants. The final sample included 36 individuals with PTSD without childhood maltreatment exposure (PTSD, n = 30 with mTBI), 31 individuals with PTSD and childhood maltreatment exposure (PTSD + M, n = 26 with mTBI), and 114 healthy comparison participants without history of childhood maltreatment exposure (HC). The PTSD + M and PTSD groups demonstrated cortical thinning in prefrontal and occipital regions, and poorer verbal memory and processing speed compared to the HC group. PTSD + M participants demonstrated cortical thinning in frontal and cingulate regions, and poorer executive functioning relative to the PTSD and HC groups. Thus, neurocognitive features varied between individuals with PTSD who did versus did not have exposure to childhood maltreatment, highlighting the need to assess developmental history of maltreatment when examining biomarkers in PTSD.

We present an open-source framework for pulmonary fissure completeness assessment. Fissure incompleteness has been shown to associate with emphysema treatment outcomes, motivating the development of tools that facilitate completeness estimation. Generally, the task of fissure completeness assessment requires accurate detection of fissures and definition of the boundary surfaces separating the lung lobes. The framework we describe acknowledges a) the modular nature of fissure detection and lung lobe segmentation (lobe boundary detection), and b) that methods to address these challenges are varied and continually developing. It is designed to be readily deployable on existing lung lobe segmentation and fissure detection data sets. The framework consists of multiple components: a flexible quality control module that enables rapid assessment of lung lobe segmentations, an interactive lobe segmentation tool exposed through 3D Slicer for handling challenging cases, a flexible fissure representation using particles-based sampling that can handle fissure feature-strength or binary fissure detection images, and a module that performs fissure completeness estimation using voxel counting and a novel surface area estimation approach. We demonstrate the usage of the proposed framework by deploying on 100 cases exhibiting various levels of fissure completeness. We compare the two completeness level approaches and also compare to visual reads. The code is available to the community via github as part of the Chest Imaging Platform and a 3D Slicer extension module.

The central nervous system is known to have limited regenerative capacity. Not only does this halt the human body’s reparative processes after central nervous system lesions, but it also impedes the establishment of effective and safe therapeutic options for such patients. Despite the high prevalence of stroke and spinal cord injury in the general population, these conditions remain incurable and place a heavy burden on patients’ families and on society more broadly. Neuroregeneration and neural engineering are diverse biomedical fields that attempt reparative treatments, utilizing stem cells-based strategies, biologically active molecules, nanotechnology, exosomes and highly tunable biodegradable systems (e.g., certain hydrogels). Although there are studies demonstrating promising preclinical results, safe clinical translation has not yet been accomplished. A key gap in clinical translation is the absence of an ideal animal or ex vivo model that can perfectly simulate the human microenvironment, and also correspond to all the complex pathophysiological and neuroanatomical factors that affect functional outcomes in humans after central nervous system injury. Such an ideal model does not currently exist, but it seems that the nonhuman primate model is uniquely qualified for this role, given its close resemblance to humans. This review considers some regenerative therapies for central nervous system repair that hold promise for future clinical translation. In addition, it attempts to uncover some of the main reasons why clinical translation might fail without the implementation of nonhuman primate models in the research pipeline.

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.