Faculty Bibliography

A novel mega-analytical approach that reduced methodological variance was evaluated using a multisite diffusion tensor imaging (DTI) fractional anisotropy (FA) data by comparing white matter integrity in people with schizophrenia to controls. Methodological variance was reduced through regression of variance captured from quality assurance (QA) and by using Marchenko-Pastur Principal Component Analysis (MP-PCA) denoising. N = 192 (119 patients/73 controls) data sets were collected at three sites equipped with 3T MRI systems: GE MR750, GE HDx, and Siemens Trio. DTI protocol included five b = 0 and 60 diffusion-sensitized gradient directions (b = 1,000 s/mm(2) ). In-house DTI QA protocol data was acquired weekly using a uniform phantom; factor analysis was used to distil into two orthogonal QA factors related to: SNR and FA. They were used as site-specific covariates to perform mega-analytic data aggregation. The effect size of patient-control differences was compared to these reported by the enhancing neuro imaging genetics meta-analysis (ENIGMA) consortium before and after regressing QA variance. Impact of MP-PCA filtering was evaluated likewise. QA-factors explained approximately 3-4% variance in the whole-brain average FA values per site. Regression of QA factors improved the effect size of schizophrenia on whole brain average FA values-from Cohen's d = .53 to .57-and improved the agreement between the regional pattern of FA differences observed in this study versus ENIGMA from r = .54 to .70. Application of MP-PCA-denoising further improved the agreement to r = .81. Regression of methodological variances captured by routine QA and advanced denoising that led to a better agreement with a large mega-analytic study.

Biophysical modeling of macroscopic diffusion-weighted MRI (dMRI) signal in terms of microscopic cellular parameters holds the promise of quantifying the integrity of white matter. Unfortunately, even fairly simple multi-compartment models of proton diffusion in the white matter do not provide a unique, biophysically plausible solution. Here we report a nontrivial diffusion MRI signal dependence on echo time (TE) in human white matter in vivo. We demonstrate that such TE dependence originates from compartment-specific T2 values and that it is a promising "orthogonal measure" able to break the degeneracy in parameter estimation, and to yield important relaxation metrics robustly. We thereby enable the precise estimation of the intra- and extra-axonal water T2 relaxation values, which is precluded by a limited signal-to-noise ratio when using multi-echo relaxometry alone.

BACKGROUND:Segmentation of gliomas in multi-parametric (MP-)MR images is challenging due to their heterogeneous nature in terms of size, appearance and location. Manual tumor segmentation is a time-consuming task and clinical practice would benefit from (semi-) automated segmentation of the different tumor compartments. METHODS:We present a semi-automated framework for brain tumor segmentation based on non-negative matrix factorization (NMF) that does not require prior training of the method. L1-regularization is incorporated into the NMF objective function to promote spatial consistency and sparseness of the tissue abundance maps. The pathological sources are initialized through user-defined voxel selection. Knowledge about the spatial location of the selected voxels is combined with tissue adjacency constraints in a post-processing step to enhance segmentation quality. The method is applied to an MP-MRI dataset of 21 high-grade glioma patients, including conventional, perfusion-weighted and diffusion-weighted MRI. To assess the effect of using MP-MRI data and the L1-regularization term, analyses are also run using only conventional MRI and without L1-regularization. Robustness against user input variability is verified by considering the statistical distribution of the segmentation results when repeatedly analyzing each patient's dataset with a different set of random seeding points. RESULTS:Using L1-regularized semi-automated NMF segmentation, mean Dice-scores of 65%, 74 and 80% are found for active tumor, the tumor core and the whole tumor region. Mean Hausdorff distances of 6.1 mm, 7.4 mm and 8.2 mm are found for active tumor, the tumor core and the whole tumor region. Lower Dice-scores and higher Hausdorff distances are found without L1-regularization and when only considering conventional MRI data. CONCLUSIONS:Based on the mean Dice-scores and Hausdorff distances, segmentation results are competitive with state-of-the-art in literature. Robust results were found for most patients, although careful voxel selection is mandatory to avoid sub-optimal segmentation.

Non-negative matrix factorization (NMF) has become a widely used tool for additive parts-based analysis in a wide range of applications. As NMF is a non-convex problem, the quality of the solution will depend on the initialization of the factor matrices. In this study, the successive projection algorithm (SPA) is proposed as an initialization method for NMF. SPA builds on convex geometry and allocates endmembers based on successive orthogonal subspace projections of the input data. SPA is a fast and reproducible method, and it aligns well with the assumptions made in near-separable NMF analyses. SPA was applied to multi-parametric magnetic resonance imaging (MRI) datasets for brain tumor segmentation using different NMF algorithms. Comparison with common initialization methods shows that SPA achieves similar segmentation quality and it is competitive in terms of convergence rate. Whereas SPA was previously applied as a direct endmember extraction tool, we have shown improved segmentation results when using SPA as an initialization method, as it allows further enhancement of the sources during the NMF iterative procedure.

OBJECTIVE: To assess the association between peripheral lipid/fat profiles and cerebral gray matter (GM) and white matter (WM) in healthy Old Order Amish (OOA). METHODS: Blood lipids, abdominal adiposity, liver lipid contents, and cerebral microstructure were assessed in OOA (N = 64, 31 males/33 females, ages 18-77). Orthogonal factors were extracted from lipid and imaging adiposity measures. GM assessment used the Human Connectome Project protocol to measure whole-brain average cortical thickness. Diffusion-weighted imaging was used to derive WM fractional anisotropy and kurtosis anisotropy measurements. RESULTS: Lipid/fat measures were captured by three orthogonal factors explaining 80% of the variance. Factor one loaded on cholesterol and/or low-density lipoprotein cholesterol measurements; factor two loaded on triglyceride/liver measurements; and factor three loaded on abdominal fat measurements. A two-stage regression including age/sex (first stage) and the three factors (second stage) examined the peripheral lipid/fat effects. Factors two and three significantly contributed to WM measures after Bonferroni corrections (P < 0.007). No factor significantly contributed to GM. Blood pressure (BP) inclusion did not meaningfully alter the lipid/fat-WM relationship. CONCLUSIONS: Peripheral lipid/fat indicators were significantly and negatively associated with cerebral WM rather than with GM, independent of age and BP level. Dissecting the fat/lipid components contributing to different brain imaging parameters may open a new understanding of the body-brain connection through lipid metabolism.

OBJECTIVE: Prior studies in prostate diffusion-weighted magnetic resonance imaging (MRI) have largely explored the impact of b-value and diffusion directions on estimated diffusion coefficient D. Here we suggest varying diffusion time, t, to study time-dependent D(t) in prostate cancer, thereby adding an extra dimension in the development of prostate cancer biomarkers. METHODS: Thirty-eight patients with peripheral zone prostate cancer underwent 3-T MRI using an external-array coil and a diffusion-weighted image sequence acquired for b = 0, as well as along 12 noncollinear gradient directions for b = 500 s/mm using stimulated echo acquisition mode (STEAM) diffusion tensor imaging (DTI). For this sequence, 6 diffusion times ranging from 20.8 to 350 milliseconds were acquired. Tumors were classified as low-grade (Gleason score [GS] 3 + 3; n = 11), intermediate-grade (GS 3 + 4; n = 16), and high-grade (GS >/=4 + 3; n = 11). Benign peripheral zone and transition zone were also studied. RESULTS: Apparent diffusion coefficient (ADC) D(t) decreased with increasing t in all zones of the prostate, though the rate of decay in D(t) was different between sampled zones. Analysis of variance and area under the curve analyses suggested better differentiation of tumor grades at shorter t. Fractional anisotropy (FA) increased with t for all regions of interest. On average, highest FA was observed within GS 3 + 3 tumors. CONCLUSIONS: There is a measurable time dependence of ADC in prostate cancer, which is dependent on the underlying tissue and Gleason score. Therefore, there may be an optimal selection of t for prediction of tumor grade using ADC. Controlling t should allow ADC to achieve greater reproducibility between different sites and vendors. Intentionally varying t enables targeted exploration of D(t), a previously overlooked biophysical phenomenon in the prostate. Its further microstructural understanding and modeling may lead to novel diffusion-derived biomarkers.

The time dependence of the diffusion coefficient is a hallmark of tissue complexity at the micrometer level. Here we demonstrate how biophysical modeling, combined with a specifically tailored diffusion MRI acquisition performing diffusion tensor imaging (DTI) for varying diffusion times, can be used to determine fiber size and membrane permeability of muscle fibers in vivo. We describe the random permeable barrier model (RPBM) and its assumptions, as well as the details of stimulated echo DTI acquisition, signal processing steps, and potential pitfalls. We illustrate the RPBM method on a few pilot examples involving human subjects (previously published as well as new), such as revealing myofiber size derived from RPBM increase after training in a calf muscle, and size decrease with atrophy in shoulder rotator cuff muscle. Finally, we comment on the potential clinical relevance of our results

Tumor segmentation is a particularly challenging task in high-grade gliomas (HGGs), as they are among the most heterogeneous tumors in oncology. An accurate delineation of the lesion and its main subcomponents contributes to optimal treatment planning, prognosis and follow-up. Conventional MRI (cMRI) is the imaging modality of choice for manual segmentation, and is also considered in the vast majority of automated segmentation studies. Advanced MRI modalities such as perfusion-weighted imaging (PWI), diffusion-weighted imaging (DWI) and magnetic resonance spectroscopic imaging (MRSI) have already shown their added value in tumor tissue characterization, hence there have been recent suggestions of combining different MRI modalities into a multi-parametric MRI (MP-MRI) approach for brain tumor segmentation. In this paper, we compare the performance of several unsupervised classification methods for HGG segmentation based on MP-MRI data including cMRI, DWI, MRSI and PWI. Two independent MP-MRI datasets with a different acquisition protocol were available from different hospitals. We demonstrate that a hierarchical non-negative matrix factorization variant which was previously introduced for MP-MRI tumor segmentation gives the best performance in terms of mean Dice-scores for the pathologic tissue classes on both datasets.