Anisotropic phantom to calibrate high-q diffusion MRI methods

A silicon oil-filled glass capillary array is proposed as an anisotropic diffusion MRI phantom. Together with a computational/theoretical pipeline these provide a gold standard for calibrating and validating high-q diffusion MRI experiments. The phantom was used to test high angular resolution diffusion imaging (HARDI) and double pulsed-field gradient (d-PFG) MRI acquisition schemes. MRI-based predictions of microcapillary diameter using both acquisition schemes were compared with results from optical microscopy. This phantom design can be used for quality control and quality assurance purposes and for testing and validating proposed microstructure imaging experiments and the processing pipelines used to analyze them.

Analysis of the effects of noise, DWI sampling, and value of assumed parameters in diffusion MRI models

Purpose

This study was a systematic evaluation across different and prominent diffusion MRI models to better understand the ways in which scalar metrics are influenced by experimental factors, including experimental design (diffusion‐weighted imaging [DWI] sampling) and noise.

Methods

Four diffusion MRI models—diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), mean apparent propagator MRI (MAP‐MRI), and neurite orientation dispersion and density imaging (NODDI)—were evaluated by comparing maps and histogram values of the scalar metrics generated using DWI datasets obtained in fixed mouse brain with different noise levels and DWI sampling complexity. Additionally, models were fit with different input parameters or constraints to examine the consequences of model fitting procedures.

Results

Experimental factors affected all models and metrics to varying degrees. Model complexity influenced sensitivity to DWI sampling and noise, especially for metrics reporting non‐Gaussian information. DKI metrics were highly susceptible to noise and experimental design. The influence of fixed parameter selection for the NODDI model was found to be considerable, as was the impact of initial tensor fitting in the MAP‐MRI model.

Conclusion

Across DTI, DKI, MAP‐MRI, and NODDI, a wide range of dependence on experimental factors was observed that elucidate principles and practical implications for advanced diffusion MRI. Magn Reson Med 78:1767–1780, 2017. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Prefrontal GABA Levels Measured With Magnetic Resonance Spectroscopy in Patients With Psychosis and Unaffected Siblings

OBJECTIVE

The authors sought to compare GABA levels in treated and untreated patients with psychosis with levels in their unaffected siblings and healthy control subjects, and to assess the effects of antipsychotic medications on GABA levels.

METHOD

GABA+ levels (i.e., including signal from unrelated proteins or macromolecules) referenced to creatine or water were studied with J-edited proton spectroscopy in the dorsal anterior cingulate cortex of 289 individuals: 184 healthy control subjects, 83 treated patients with psychosis, 25 untreated patients, 31 unaffected siblings, and 17 patients studied both while off all medications and while on a single antipsychotic.

RESULTS

GABA+ levels in the dorsal anterior cingulate did not differ between untreated patients and healthy controls. For treated patients, levels were modestly lower for GABA+/creatine but did not differ for GABA+/water compared with healthy controls. For both GABA+ measures, unaffected siblings had significantly lower levels compared with controls. GABA+/creatine showed a modest degree of familiality (intraclass correlation=0.36). Antipsychotic dosage was negatively correlated with GABA+ levels, but the on-off medication studies indicated no difference in GABA+ levels on antipsychotics compared with off antipsychotics.

CONCLUSIONS

GABA+/creatine in the dorsal anterior cingulate may constitute an intermediate phenotype with low effect size for psychosis, but GABA+/water measures do not fully support this conclusion. Low GABA+ levels in unaffected siblings could suggest a genetic association, but the failure to find consistent evidence of this phenotype in the patients themselves weakens genetic inference on risk for psychosis. Replication in independent samples of siblings is warranted to confirm the potential genetic risk association.

Investigation of anatomical thalamo-cortical connectivity and FMRI activation in schizophrenia

The purpose of this study was to examine measures of anatomical connectivity between the thalamus and lateral prefrontal cortex (LPFC) in schizophrenia and to assess their functional implications. We measured thalamocortical connectivity with diffusion tensor imaging (DTI) and probabilistic tractography in 15 patients with schizophrenia and 22 age- and sex-matched controls. The relationship between thalamocortical connectivity and prefrontal cortical blood-oxygenation-level-dependent (BOLD) functional activity as well as behavioral performance during working memory was examined in a subsample of 9 patients and 18 controls. Compared with controls, schizophrenia patients showed reduced total connectivity of the thalamus to only one of six cortical regions, the LPFC. The size of the thalamic region with at least 25% of model fibers reaching the LPFC was also reduced in patients compared with controls. The total thalamocortical connectivity to the LPFC predicted working memory task performance and also correlated with LPFC BOLD activation. Notably, the correlation with BOLD activation was accentuated in patients as compared with controls in the ventral LPFC. These results suggest that thalamocortical connectivity to the LPFC is altered in schizophrenia with functional consequences on working memory processing in LPFC.

Reproducibility of prefrontal γ-aminobutyric acid measurements with J-edited spectroscopy

γ-Aminobutyric acid (GABA) is the chief inhibitory neurotransmitter of the human brain, and GABA-ergic dysfunction has been implicated in a variety of neuropsychiatric disorders. Recent MRS techniques have allowed the quantification of GABA concentrations in vivo, and could therefore provide biologically relevant information. Few reports have formally characterized the reproducibility of these techniques, and differences in field strength, acquisition and processing parameters may result in large differences in measured GABA values. Here, we used a J-edited, single-voxel spectroscopy method of measurement of GABA + macromolecules (GABA + ) in the anterior cingulate cortex (ACC) and right frontal white matter (rFWM) at 3 T. We measured the coefficient of variation within subjects (CVw) and intra-class correlation coefficients on two repeated scans obtained from 10 healthy volunteers with processing procedures developed in-house for the quantification of GABA + and other major metabolites. In addition, by segmenting the spectroscopic voxel into cerebrospinal fluid, gray matter and white matter, and employing a linear regression technique to extrapolate metabolite values to pure gray and white matter, we determined metabolite differences between gray and white matter in ACC and rFWM. CVw values for GABA + /creatine, GABA + /H(2) O, GABA + , creatine, partially co-edited glutamate + glutamine (Glx)/creatine, partially co-edited Glx and N-acetylaspartic acid (NAA)/creatine were all below 12% in both ACC and rFWM. After extrapolation to pure gray and pure white matter, CVw values for all metabolites were below 16%. We found metabolite ratios between gray and white matter for GABA + /creatine, GABA + , creatine, partially co-edited Glx and NAA/creatine to be 0.88 ± 0.21 (standard deviation), 1.52 ± 0.32, 1.77 ± 0.4, 2.69 ± 0.74 and 0.70 ± 0.05, respectively. This study validates a reproducible method for the quantification of brain metabolites, and provides information on gray/white matter differences that may be important in the interpretation of results in clinical populations.

Genetic modulation of GABA levels in the anterior cingulate cortex by GAD1 and COMT

Gamma-aminobutyric acid (GABA)-ergic transmission is critical for normal cortical function and is likely abnormal in a variety of neuropsychiatric disorders. We tested the in vivo effects of variations in two genes implicated in GABA function on GABA concentrations in prefrontal cortex of living subjects: glutamic acid decarboxylase 1 (GAD1), which encodes GAD67, and catechol-o-methyltransferase (COMT), which regulates synaptic dopamine in the cortex. We studied six single nucleotide polymorphisms (SNPs) in GAD1 previously associated with risk for schizophrenia or cognitive dysfunction and the val158met polymorphism in COMT in 116 healthy volunteers using proton magnetic resonance spectroscopy. Two of the GAD1 SNPs (rs1978340 (p=0.005) and rs769390 (p=0.004)) showed effects on GABA levels as did COMT val158met (p=0.04). We then tested three SNPs in GAD1 (rs1978340, rs11542313, and rs769390) for interaction with COMT val158met based on previous clinical results. In this model, rs11542313 and COMT val158met showed significant main effects (p=0.001 and 0.003, respectively) and a trend toward a significant interaction (p=0.05). Interestingly, GAD1 risk alleles for schizophrenia were associated with higher GABA/Cre, and Val-Val homozygotes had high GABA/Cre levels when on a GAD1 risk genotype background (N=6). These results support the importance of genetic variation in GAD1 and COMT in regulating prefrontal cortical GABA function. The directionality of the effects, however, is inconsistent with earlier evidence of decreased GABA activity in schizophrenia.

Detection and characterization of neurotoxicity in cancer patients using proton MR spectroscopy

OBJECTIVE

The study objective was to detect abnormalities and identify relationships between brain metabolic ratios determined by proton magnetic resonance spectroscopic imaging ((1)H-MRSI) and neuropsychological (NP) function in cancer patients at risk for neurotoxicity.

METHODS

Thirty-two patients received (1)H-MRSI using a multi-slice, multi-voxel technique on a 1.5T magnet. Cho/NAA, NAA/Cr, and Cho/Cr ratios were identified in seven pre-determined sites without tumor involvement. A battery of age-appropriate NP tests was administered within 7 days of imaging. Relationships were examined between test scores and metabolite ratios.

CONCLUSIONS

This study identifies relationships between brain metabolite ratios and cognitive functioning in cancer patients. (1)H-MRSI may be useful in early detection of neurotoxic effects, but prospective longitudinal studies in a homogeneous population are recommended to determine the prognostic value.

Three-dimensional mapping of lingual myoarchitecture by diffusion tensor MRI

This study was performed to assess the feasibility of investigating the complex lingual myoarchitecture through segmentation of muscles from diffusion tensor imaging (DTI) data. The primary eigenvectors were found to be adequate for delineating the superior and inferior longitudinalis, genioglossus, and hyoglossus. The tertiary eigenvector orientations effectively revealed the homogeneous and systematic change of muscle orientation in the tongue core. In the longitudinalis near the tongue tip, the secondary eigenvectors were oriented in the radial direction. Lingual muscles were segmented using two methods: modified directional correlation (DC) and tensor coherence (TC) methods. The DC method, based on one eigenvector, was found to be inadequate for lingual muscle segmentation, whereas the TC method, based on the tensor shape and orientation, was used successfully to segment most lingual muscles. The segmentation result was used to report the diffusion tensor properties of individual lingual muscles. Also found was a continuous change in skewness of the intrinsic tongue core from negative in the anterior region to positive in the posterior region. DTI and the proposed segmentation method provide an adequate means of imaging and visualizing the complex, compartmentalized musculature of the tongue. The potential for in vivo research and clinical applications is demonstrated.

Regional distribution of measurement error in diffusion tensor imaging

The characterization of measurement error is critical in assessing the significance of diffusion tensor imaging (DTI) findings in longitudinal and cohort studies of psychiatric disorders. We studied 20 healthy volunteers, each one scanned twice (average interval between scans of 51 +/- 46.8 days) with a single shot echo planar DTI technique. Intersession variability for fractional anisotropy (FA) and Trace (D) was represented as absolute variation (standard deviation within subjects: SDw), percent coefficient of variation (CV) and intra-class correlation coefficient (ICC). The values from the two sessions were compared for statistical significance with repeated measures analysis of variance or a non-parametric equivalent of a paired t-test. The results showed good reproducibility for both FA and Trace (CVs below 10% and ICCs at or above 0.70 in most regions of interest) and evidence of systematic global changes in Trace between scans. The regional distribution of reproducibility described here has implications for the interpretation of regional findings and for rigorous pre-processing. The regional distribution of reproducibility measures was different for SDw, CV and ICC. Each one of these measures reveals complementary information that needs to be taken into consideration when performing statistical operations on groups of DT images.

Dependence on diffusion time of apparent diffusion tensor of ex vivo calf tongue and heart

The time dependence of the apparent diffusion tensor of ex vivo calf heart and tongue was measured for diffusion times (tau(d)) between 32 and 810 ms. The results showed evidence of restricted diffusion in the muscle tissues of both organs. In regions where the myofibers are parallel, the largest eigenvalue (lambda(1)) of the diffusion tensor remained the same for all diffusion times measured, while the other eigenvalues (lambda(2), lambda(3)) decreased by 29-36% between tau(d) = 32 ms and tau(d) = 400 ms. In regions where the fibers cross, the lambda(1) also changed, decreasing by 17% between tau(d) = 32 ms and tau(d) = 400 ms. The restricting compartment size and volume fraction were effectively estimated by fitting the time courses of the eigenvalues to a model consisting of a nonrestricted compartment and a cylindrically restricted compartment. To our knowledge, this study is the first demonstrating diffusion time dependence of measured water diffusion tensor in muscular tissue. With improvement in scanning technology, future studies may permit noninvasive, in vivo detection of changes in muscle myoarchitecture due to disease, treatment, and exercise.

Estimating intensity variance due to noise in registered images: Applications to diffusion tensor MRI

Image registration techniques which require image interpolation are widely used in neuroimaging research. We show that signal variance in interpolated images differs significantly from the signal variance of the original images in native space. We describe a simple approach to compute the signal variance in registered images based on the signal variance and covariance of the original images, the spatial transformations computed by the registration procedure, and the interpolation or approximation kernel chosen. The method is general and could handle various sources of signal variability, such as thermal noise and physiological noise, provided that their effects can be assessed in the original images. Our approach is applied to diffusion tensor (DT) MRI data, assuming only thermal noise as the source of variability in the data. We show that incorrect noise variance estimates in registered diffusion-weighted images can affect DT parameters, as well as indices of goodness of fit such as chi-square maps. In addition to DT-MRI, we believe that this methodology would be useful any time parameter extraction methods are applied to registered or interpolated data, such as in relaxometry and functional MRI studies.

Comprehensive approach for correction of motion and distortion in diffusion-weighted MRI

Patient motion and image distortion induced by eddy currents cause artifacts in maps of diffusion parameters computed from diffusion-weighted (DW) images. A novel and comprehensive approach to correct for spatial misalignment of DW imaging (DWI) volumes acquired with different strengths and orientations of the diffusion sensitizing gradients is presented. This approach uses a mutual information-based registration technique and a spatial transformation model containing parameters that correct for eddy current-induced image distortion and rigid body motion in three dimensions. All parameters are optimized simultaneously for an accurate and fast solution to the registration problem. The images can also be registered to a normalized template with a single interpolation step without additional computational cost. Following registration, the signal amplitude of each DWI volume is corrected to account for size variations of the object produced by the distortion correction, and the b-matrices are properly recalculated to account for any rotation applied during registration. Both qualitative and quantitative results show that this approach produces a significant improvement of diffusion tensor imaging (DTI) data acquired in the human brain.

Functional anatomy of cognitive development: fMRI of verbal fluency in children and adults

OBJECTIVE

To identify age-dependent activation patterns of verbal fluency with functional MRI (fMRI).

BACKGROUND

Few fMRI language studies have been performed in children, and none provide comparison data to adult studies. Normative data are important for interpretation of similar studies in patients with epilepsy.

METHODS

A total of 10 normal children (5 boys, 5 girls; mean age, 10.7 years; range, 8.1 to 13.1 years) and 10 normal adults (5 men, 5 women; mean age, 28.7 years; range, 19.3 to 48 years) were studied on a 1.5-T Signa MRI scanner using BOLD echo planar imaging of the frontal lobes with a verbal fluency paradigm, covert word generation to letters. Studies were analyzed with a cross-correlation algorithm (r = 0.7). A region-of-interest analysis was used to determine the extent, magnitude, and laterality of brain activation.

RESULTS

Children and adults activated similar regions, predominantly in left inferior frontal cortex (Broca's area) and left middle frontal gyrus (dorsolateral prefrontal cortex). Children had, on average, 60% greater extent of activation than adults, with a trend for greater magnitude of activation. Children also had significantly more right hemisphere and inferior frontal gyrus activation than adults.

CONCLUSIONS

In a test of verbal fluency, children tended to activate cortex more widely than adults, but activation patterns for fluency appear to be established by middle childhood. Thus, functional MRI using verbal fluency paradigms may be applied to pediatric patient populations for determining language dominance in anterior brain regions. The greater activation found in children, including the right inferior frontal gyrus, may reflect developmental plasticity for the ongoing organization of neural networks, which underlie language capacity.

Comparative MR imaging study of brain maturation in kittens with T1, T2, and the trace of the diffusion tensor

PURPOSE

To assess the time-course of the relaxation times and the orientationally averaged water diffusion coefficient Doav in postnatal brain development.

MATERIALS AND METHODS

Multisection maps of T1, T2, and the trace of the diffusion tensor (Trace[D] = 3 x Doav) were obtained in four kittens at eight time points.

RESULTS

In the adult, Doav was about 700 micron 2/sec in both white and gray matter. In the newborn, Doav was 1,100-1,350 micron 2/sec in white matter and 1,000 micron 2/sec in gray matter. For all anatomic regions and time points, the correlation between Doav and 1/T2 was high (R2 = 0.87, P << .001). T1 showed a lower correlation with Doav and a higher sensitivity to myelinization than did T2.

CONCLUSION

Although Doav shows dramatic changes in the maturing brain, the high correlation between Doav and T2 indicates that little additional information can be obtained by measuring this diffusion parameter during normal brain development. This contrasts with previous findings in brain ischemia, where Doav and T2 appear to be uncorrelated. After including the authors' data and published iontophoretic measurements in a simple model of diffusion in tissues, the authors suggest that the underlying mechanisms of Doav reduction in brain maturation and ischemia are different. Doav changes during development are mainly affected by events occurring in the cellular compartment, while changes in extracellular volume fraction and tortuosity, which are thought to determine the reduction in Doav during ischemia, are probably of secondary importance.

Characterization of and correction for eddy current artifacts in echo planar diffusion imaging

Magnetic resonance diffusion imaging is potentially an important tool for the noninvasive characterization of normal and pathological tissue. The technique, however, is prone to a number of artifacts that can severely affect its ability to provide clinically useful information. In this study, the problem of eddy current-induced geometric distortions that occur in diffusion images acquired with echo planar sequences was addressed. These geometric distortions produce artifacts in computed maps of diffusion parameters and are caused by misalignments in the individual diffusion-weighted images that comprise the diffusion data set. A new approach is presented to characterize and calibrate the eddy current effects, enabling the eddy current distortions to be corrected in sets of interleaved (or snapshot) echo planar diffusion images. Correction is achieved by acquiring one-dimensional field maps in the read and phase encode direction for each slice and each diffusion step. The method is then demonstrated through the correction of distortions in diffusion images of the human brain. It is shown that by using the eddy current correction scheme outlined, the eddy current-induced artifacts in the diffusion-weighted images are almost completely eliminated. In addition, there is a significant improvement in the quality of the resulting diffusion tensor maps.

Regionally specific pattern of neurochemical pathology in schizophrenia as assessed by multislice proton magnetic resonance spectroscopic imaging

OBJECTIVE

Several single-voxel proton magnetic resonance spectroscopy (1H-MRS) studies of patients with schizophrenia have found evidence of reductions of N-acetyl-aspartate (NAA) concentrations in the temporal lobes. Multislice proton magnetic resonance spectroscopy imaging (1H-MRSI) permits simultaneous acquisition and mapping of NAA, choline-containing compounds (CHO), and creatine/phosphocreatine (CRE) signal intensities from multiple whole brain slices consisting of 1.4-ml single-volume elements. We have used 1H-MRSI to assess the regional specificity of previously reported changes of metabolite signal intensities in schizophrenia. Hippocampal volume was also measured to test the relationship between 1H-MRSI findings and tissue volume in this region.

METHOD

Ratios of areas under the metabolite peaks of the proton spectra were determined (i.e., NAA/CRE, NAA/CHO, CHO/CRE) for multiple cortical and subcortical regions in 10 inpatients with schizophrenia.

RESULTS

Patients showed significant reductions of NAA/CRE and NAA/CHO bilaterally in the hippocampal region and in the dorsolateral prefrontal cortex. There were no significant changes in CHO/CRE or in NAA ratios in any other area sampled. No significant correlation was found between metabolite ratios in the hippocampal region and its volume.

CONCLUSIONS

NAA-relative signal intensity reductions in schizophrenia appear to be remarkably localized, involving primarily the hippocampal region and the dorsolateral prefrontal cortex, two regions implicated prominently in the pathophysiology of this disorder.

Proton magnetic resonance spectroscopic imaging in patients with cerebellar degeneration

Using proton magnetic resonance spectroscopic imaging, we studied the cerebellum of 9 patients with cerebellar degeneration and of 9 age-matched normal control subjects. This technique permits the simultaneous measurement of N-acetylaspartate, choline-containing compounds, creatine/phosphocreatine, and lactate signal intensities from four 15-mm slices divided into 0.84-ml single-volume elements. Because patients with cerebellar degeneration often show substantial atrophy on magnetic resonance imaging (MRI), we specifically chose to analyze the spectroscopic signals only from tissue that did not have an atrophic appearance on the MRI. The spectroscopic findings showed a significant reduction of N-acetylaspartate in all parts of the cerebellum, a significant correlation with MRI scores of cerebellar atrophy, and a significant correlation with clinical rating scores of cerebellar disturbance. Our method of analysis suggests the presence of a neurodegenerative process in cerebellar areas that do not appear to be atrophic on the MRI. Some limitations of proton magnetic resonance spectroscopic imaging in the present study were related to the partial field inhomogeneity characteristics of the posterior fossa, the anatomical location of the cerebellum, and the particularly severe cerebellar atrophy in some of the patients.

Reproducibility of proton MR spectroscopic imaging findings

PURPOSE

To study the intraindividual, interindividual, and intraregional reproducibility of multisection proton MR spectroscopic imaging in healthy adults.

METHODS

Six subjects were studied three times with proton MR spectroscopic imaging. Multisection long-echo-time proton MR spectroscopic imaging permits simultaneous acquisition of N-acetylaspartate (NAA), choline-containing compounds (Cho), and creatine plus phosphocreatine (Cr) signal intensities from four 15-mm-thick sections divided into 0.84-mL single-volume elements. Regions of interest were the frontal cortex, the occipital cortex, the parietal cortex, the insular cortex, the cingulate gyrus, the centrum semiovale, the thalamus, and the caudate. Statistical evaluation was performed by analyses of variance and components of variance method.

RESULTS

The ratio NAA/Cr showed the lowest overall coefficient of variation (CV, %) in most of the regions of interest (range, 8.9 to 26.1). Interregional differences in the overall CV were present. Interindividual CVs ranged from 4.2 to 8.7 for NAA/Cr, from 6.8 to 17.4 for NAA/Cho, and from 5.0 to 13.6 for Cho/Cr. Intraindividual CVs ranged from 8.2 to 22.2 for NAA/Cr, from 12.8 to 25.8 for NAA/Cho, and from 4.5 to 21.0 for Cho/Cr. Intraregional CVs ranged from 12.3 to 21.2 for NAA/Cr, from 13.0 to 20.4 for NAA/Cho, and from 12.2 to 18.9 for Cho/Cr.

CONCLUSIONS

Proton MR spectroscopic imaging showed good overall reproducibility. The finding of interregional variations of CV indicates that care is needed when using this imaging technique for follow-up studies.

Blue blood or black blood: R1 effects in gradient-echo echo-planar functional neuroimaging

Changes in the longitudinal relaxation rate (R1) may play a role in the MRI signal intensity increases that have been associated with physiological brain activation. We used gradient-echo echo-planar MRI (GRE-EPI) to test whether physiological activations associated with hypercapnia in dogs were dependent on the delay (TR) between successive images in a time-series. Our results show that, in addition to activation-induced changes in the R2 (transverse relaxation including inhomogeneity effects), activation-induced changes in R1 are significant under certain pulsing conditions. In our paradigm, the R1 contribution became significant at TR values of 1 s or less.

Cerebral white matter in the centrum semiovale exhibits a larger N-acetyl signal than does gray matter in long echo time 1H-magnetic resonance spectroscopic imaging

Long echo time (272 ms) 1H magnetic resonance spectroscopic imaging was used to measure the relative magnitudes of the N-acetylaspartate (NAA) signal in a variety of anatomically defined brain structures (centrum semiovale, thalamus, medial frontal cortex, and genu of the corpus callosum) composed primarily of gray matter or white matter. Six normal young adult humans aged 30-40 were studied. With a 95% level of statistical confidence, the white matter in the centrum semiovale (CSO) produced a more intense NAA signal than did the gray matter in the thalamus and the frontal cortex. Differences between the white matter regions were also noted. The CSO white matter's NAA signal yielded a larger NAA signal than did the white matter of the genu of the corpus callosum. Possible reasons for the anatomical variation in the cerebral NAA signal intensity are discussed.