Deep brain stimulation (DBS) has the potential to improve the quality of life of people with a variety of neurological diseases. A key challenge in DBS is in the placement of a stimulation electrode in the anatomical location that maximizes efficacy and minimizes side effects. Pre-operative localization of the optimal stimulation zone can reduce surgical times and morbidity. Current methods of producing efficacy probability maps follow an anatomical guidance on magnetic resonance imaging (MRI) to identify the areas with the highest efficacy in a population. In this work, we propose to revisit this problem as a classification problem, where each voxel in the MRI is a sample informed by the surrounding anatomy. We use a patch-based convolutional neural network to classify a stimulation coordinate as having a positive reduction in symptoms during surgery. We use a cohort of 187 patients with a total of 2,869 stimulation coordinates, upon which 3D patches were extracted and associated with an efficacy score. We compare our results with a registration-based method of surgical planning. We show an improvement in the classification of intraoperative stimulation coordinates as a positive response in reduction of symptoms with AUC of 0.670 compared to a baseline registration-based approach, which achieves an AUC of 0.627 (p < 0.01). Although additional validation is needed, the proposed classification framework and deep learning method appear well-suited for improving pre-surgical planning and personalize treatment strategies.
Diffusion weighted MRI (DW-MRI) depends on accurate quantification signal intensities that reflect directional apparent diffusion coefficients (ADC). Signal drift and fluctuations during imaging can cause systematic non-linearities that manifest as ADC changes if not corrected. Here, we present a case study on a large longitudinal dataset of typical diffusion tensor imaging. We investigate observed variation in the cerebral spinal fluid (CSF) regions of the brain, which should represent compartments with isotropic diffusivity. The study contains 3949 DW-MRI acquisitions of the human brain with 918 subjects and 542 with repeated scan sessions. We provide an analysis of the inter-scan, inter-session, and intra-session variation and an analysis of the associations with the applied diffusion gradient directions. We investigate a hypothesis that CSF models could be used in lieu of an interspersed minimally diffusion-weighted image (b0) correction. Variation in CSF signal is not largely attributable to within-scan dynamic anatomical changes (3.6%), but rather has substantial variation across scan sessions (10.6%) and increased variation across individuals (26.6%). Unfortunately, CSF intensity is not solely explained by a main drift model or a gradient model, but rather has statistically significant associations with both possible explanations. Further exploration is necessary for CSF drift to be used as an effective harmonization technique.
KEYWORDS: Diffusion, Gold, Signal to noise ratio, Data modeling, Magnetic resonance imaging, Biological research, Brain, Spherical lenses, In vivo imaging, Data acquisition
The diffusion tensor model is nonspecific in regions where micrometer structural patterns are inconsistent at the millimeter scale (i.e., brain regions with pathways that cross, bend, branch, fan, etc.). Numerous models have been proposed to represent crossing fibers and complex intravoxel structure from in vivo diffusion weighted magnetic resonance imaging (e.g., high angular resolution diffusion imaging—HARDI). Here, we present an empirical comparison of two HARDI approaches—persistent angular structure MRI (PAS-MRI) and Q-ball—using a newly acquired reproducibility dataset. Briefly, a single subject was scanned 11 times with 96 diffusion weighted directions and 10 reference volumes for each of two b values (1000 and 3000 s / mm2 for a total of 2144 volumes). Empirical reproducibility of intravoxel fiber fractions (number/strength of peaks), angular orientation, and fractional anisotropy was compared with metrics from a traditional tensor analysis approach, focusing on b values of 1000 and 3000 s / mm2. PAS-MRI is shown to be more reproducible than Q-ball and offers advantages at low b values. However, there are substantial and biologically meaningful differences between the intravoxel structures estimated both in terms of analysis method as well as by b value. The two methods suggest a fundamentally different microarchitecture of the human brain; therefore, it is premature to perform meta-analysis or combine results across HARDI studies using a different analysis model or acquisition sequences.
The choice of surface template plays an important role in cross-sectional subject analyses involving cortical brain surfaces because there is a tendency toward registration bias given variations in inter-individual and inter-group sulcal and gyral patterns. In order to account for the bias and spatial smoothing, we propose a feature-based unbiased average template surface. In contrast to prior approaches, we factor in the sample population covariance and assign weights based on feature information to minimize the influence of covariance in the sampled population. The mean surface is computed by applying the weights obtained from an inverse covariance matrix, which guarantees that multiple representations from similar groups (e.g., involving imaging, demographic, diagnosis information) are down-weighted to yield an unbiased mean in feature space. Results are validated by applying this approach in two different applications. For evaluation, the proposed unbiased weighted surface mean is compared with un-weighted means both qualitatively and quantitatively (mean squared error and absolute relative distance of both the means with baseline). In first application, we validated the stability of the proposed optimal mean on a scan-rescan reproducibility dataset by incrementally adding duplicate subjects. In the second application, we used clinical research data to evaluate the difference between the weighted and unweighted mean when different number of subjects were included in control versus schizophrenia groups. In both cases, the proposed method achieved greater stability that indicated reduced impacts of sampling bias. The weighted mean is built based on covariance information in feature space as opposed to spatial location, thus making this a generic approach to be applicable to any feature of interest.
An understanding of the bias and variance of diffusion weighted magnetic resonance imaging (DW-MRI) acquisitions across scanners, study sites, or over time is essential for the incorporation of multiple data sources into a single clinical study. Studies that combine samples from various sites may be introducing confounding factors due to site-specific artifacts and patterns. Differences in bias and variance across sites may render the scans incomparable, and, without correction, inferences obtained from these data may be misleading. We present an analysis of the bias and variance of scans of the same subjects across different sites and evaluate their impact on statistical analyses. In previous work, we presented a simulation extrapolation (SIMEX) technique for bias estimation as well as a wild bootstrap technique for variance estimation in metrics obtained from a Q-ball imaging (QBI) reconstruction of empirical high angular resolution diffusion imaging (HARDI) data. We now apply those techniques to data acquired from 5 healthy volunteers on 3 independent scanners under closely matched acquisition protocols. The bias and variance of GFA measurements were estimated on a voxel-wise basis for each scan and compared across study sites to identify site-specific differences. Further, we provide model recommendations that can be used to determine the extent of the impact of bias and variance as well as aspects of the analysis to account for these differences. We include a decision tree to help researchers determine if model adjustments are necessary based on the bias and variance results.
High Angular Resolution Diffusion Imaging (HARDI) models are used to capture complex intra-voxel microarchitectures. The magnetic resonance imaging sequences that are sensitized to diffusion are often highly accelerated and prone to motion, physiologic, and imaging artifacts. In diffusion tensor imaging, robust statistical approaches have been shown to greatly reduce these adverse factors without human intervention. Similar approaches would be possible with HARDI methods, but robust versions of each distinct HARDI approach would be necessary. To avoid the computational and pragmatic burdens of creating individual robust HARDI analysis variants, we propose a robust outlier imputation model to mitigate outliers prior to traditional HARDI analysis. This model uses a weighted spherical harmonic fit of diffusion weighted magnetic resonance imaging scans to estimate the values which had been corrupted during acquisition to restore them. Briefly, spherical harmonics of 6th order were used to generate basis function which were weighted by diffusion signal for detection of outliers. For validation, a single healthy volunteer was scanned for a single session comprising of two scans one without head movement and the other with deliberate head movement at a b-value of 3000 s/mm2 with 64 diffusion weighted directions with a single b0 (5 averages) per scan. The deliberate motion from the volunteer created natural artifacts in the acquisition of one of the scans. The imputation model shows reduction in root mean squared error of the raw signal intensities and improvement for the HARDI method Q-ball in terms of the Angular Correlation Coefficient. The results reveal that there is quantitative and qualitative improvement. The proposed model can be used as general pre-processing model before implementing any HARDI model in general to restore the artifacts which are created because of the outlier diffusion signal in certain gradient volumes.
Pathologic complete response following neoadjuvant therapy (NAT) is used as a short-term surrogate marker of eventual outcome in patients with breast cancer. Analyzing voxel-level heterogeneity in MRI-derived parametric maps, obtained before and after the first cycle of NAT (n=33), in conjunction with receptor status, may improve the predictive accuracy of tumor response to NAT. Toward that end, we incorporated two MRI-derived parameters, the apparent diffusion coefficient and efflux rate constant, with receptor status in a logistic ridge-regression model. The area under the curve (AUC) and Brier score of the model computed via 10-fold cross validation were 0.94 (95% CI: 0.85, 0.99) and 0.11 (95% CI: 0.06, 0.16), respectively. These two statistics strongly support the hypothesis that our proposed model outperforms the other models that we investigated (namely, models without either receptor information or voxel-level information). The contribution of the receptor information was manifested by an 8% to 15% increase in AUC and a 14% to 21% decrease in Brier score. These data indicate that combining multiparametric MRI with hormone receptor status has a high likelihood of improved prediction of pathologic response to NAT in breast cancer.
KEYWORDS: Magnetic resonance imaging, Breast cancer, Tumors, Diagnostics, Temporal resolution, Data modeling, Statistical analysis, Receptors, Error analysis, Cancer
This meta-analysis assesses the prognostic value of quantitative dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and diffusion-weighted MRI (DW-MRI) performed during neoadjuvant therapy (NAT) of locally advanced breast cancer. A systematic literature search was conducted to identify studies of quantitative DCE-MRI and DW-MRI performed during breast cancer NAT that report the sensitivity and specificity for predicting pathological complete response (pCR). Details of the study population and imaging parameters were extracted from each study for subsequent meta-analysis. Metaregression analysis, subgroup analysis, study heterogeneity, and publication bias were assessed. Across 10 studies that met the stringent inclusion criteria for this meta-analysis (out of 325 initially identified studies), we find that MRI had a pooled sensitivity of 0.91 [95% confidence interval (CI), 0.80 to 0.96] and specificity of 0.81(95% CI, 0.68 to 0.89) when adjusted for covariates. Quantitative DCE-MRI exhibits greater specificity for predicting pCR than semiquantitative DCE-MRI (p<0.001). Quantitative DCE-MRI and DW-MRI are able to predict, early in the course of NAT, the eventual response of breast tumors, with a high level of specificity and sensitivity. However, there is a high degree of heterogeneity in published studies highlighting the lack of standardization in the field.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.