Search

Objectives The aim of this study was to assess the effect of a deep learning (DL)–based reconstruction algorithm on late gadolinium enhancement (LGE) image quality and to evaluate its influence on scar quantification. Methods Sixty patients (46 ± 17 years, 50% male) with suspected or known cardiomyopathy underwent CMR. Short-axis LGE images were reconstructed using the conventional reconstruction and a DL network (DLRecon) with tunable noise reduction (NR) levels from 0 to 100%. Image quality of standard LGE images and DLRecon images with 75% NR was scored using a 5-point scale (poor to excellent). In 30 patients with LGE, scar size was quantified using thresholding techniques with different standard deviations (SD) above remote myocardium, and using full width at half maximum (FWHM) technique in images with varying NR levels. Results DLRecon images were of higher quality than standard LGE images (subjective quality score 3.3 ± 0.5 vs. 3.6 ± 0.7, p p values p = 0.06). Conclusions LGE image quality improved significantly using a DL-based reconstruction algorithm. However, this algorithm has an important impact on scar quantification depending on which quantification technique is used. The FWHM method is preferred because of its independency of NR. Clinicians should be aware of this impact on scar quantification, as DL-based reconstruction algorithms are being used. Key Points • The image quality based on (subjective) visual assessment and image sharpness of late gadolinium enhancement images improved significantly using a deep learning–based reconstruction algorithm that aims to reconstruct high signal-to-noise images using a denoising technique. • Special care should be taken when scar size is quantified using thresholding techniques with different standard deviations above remote myocardium because of the large impact of these advanced image enhancement algorithms. • The full width at half maximum method is recommended to quantify scar size when deep learning algorithms based on noise reduction are used, as this method is the least sensitive to the level of noise and showed the best agreement with visual late gadolinium enhancement assessment.

PURPOSE To develop and evaluate an efficient precontrast T1 mapping technique suitable for quantitative high-resolution whole-brain dynamic contrast-enhanced-magnetic resonance imaging (DCE-MRI). METHODS Variable flip angle (VFA) T1 mapping was considered that provides 1 × 1 × 2 mm3 resolution to match a recent high-resolution whole-brain DCE-MRI protocol. Seven FAs were logarithmically spaced from 1.5° to 15°. T1 and M0 maps were estimated using model-based reconstruction. This approach was evaluated using an anatomically realistic brain tumor digital reference object (DRO) with noise-mimicking 3T neuroimaging and fully sampled data acquired from one healthy volunteer. Methods were also applied on fourfold prospectively undersampled VFA data from 13 patients with high-grade gliomas. RESULTS T1 -mapping precision decreased with undersampling factor R, althoughwhereas bias remained small before a critical R. In the noiseless DRO, T1 bias was

The U-net is a deep-learning network model that has been used to solve a number of inverse problems. In this work, the concatenation of two-element U-nets, termed the W-net, operating in k-space (K) and image (I) domains, were evaluated for multi-channel magnetic resonance (MR) image reconstruction. The two-element network combinations were evaluated for the four possible image-k-space domain configurations: a) W-net II, b) W-net KK, c) W-net IK, and d) W-net KI. Selected four element (WW-nets) and six element (WWW-nets) networks were also examined. Two configurations of each network were compared: 1) each coil channel was processed independently, and 2) all channels were processed simultaneously. One hundred and eleven volumetric, T1-weighted, 12-channel coil k-space datasets were used in the experiments. Normalized root mean squared error, peak signal-to-noise ratio and visual information fidelity were used to assess the reconstructed images against the fully sampled reference images. Our results indicated that networks that operate solely in the image domain were better when independently processing individual channels of multi-channel data. Dual-domain methods were better when simultaneously reconstructing all channels of multi-channel data. In addition, the best cascade of U-nets performed better (p 0.01) than the previously published, state-of-the-art Deep Cascade and Hybrid Cascade models in three out of four experiments.

Purpose To apply tracer kinetic models as temporal constraints during reconstruction of under‐sampled brain tumor dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI). Methods A library of concentration vs time profiles is simulated for a range of physiological kinetic parameters. The library is reduced to a dictionary of temporal bases, where each profile is approximated by a sparse linear combination of the bases. Image reconstruction is formulated as estimation of concentration profiles and sparse model coefficients with a fixed sparsity level. Simulations are performed to evaluate modeling error, and error statistics in kinetic parameter estimation in presence of noise. Retrospective under‐sampling experiments are performed on a brain tumor DCE digital reference object (DRO), and 12 brain tumor in‐vivo 3T datasets. The performances of the proposed under‐sampled reconstruction scheme and an existing compressed sensing‐based temporal finite‐difference (tFD) under‐sampled reconstruction were compared against the fully sampled inverse Fourier Transform‐based reconstruction. Results Simulations demonstrate that sparsity levels of 2 and 3 model the library profiles from the Patlak and extended Tofts‐Kety (ETK) models, respectively. Noise sensitivity analysis showed equivalent kinetic parameter estimation error statistics from noisy concentration profiles, and model approximated profiles. DRO‐based experiments showed good fidelity in recovery of kinetic maps from 20‐fold under‐sampled data. In‐vivo experiments demonstrated reduced bias and uncertainty in kinetic mapping with the proposed approach compared to tFD at under‐sampled reduction factors >= 20. Conclusions Tracer kinetic models can be applied as temporal constraints during brain tumor DCE‐MRI reconstruction. The proposed under‐sampled scheme resulted in model parameter estimates less biased with respect to conventional fully sampled DCE MRI reconstructions and parameter estimation. The approach is flexible, can use nonlinear kinetic models, and does not require tuning of regularization parameters.

Accurate and reliable quantification of brain metabolites measured in vivo using 1H magnetic resonance spectroscopy (MRS) is a topic of continued interest. Aside from differences in the basic approach to quantification, the quantification of metabolite data acquired at different sites and on different platforms poses an additional methodological challenge. In this study, spectrally edited γ-aminobutyric acid (GABA) MRS data were analyzed and GABA levels were quantified relative to an internal tissue water reference. Data from 284 volunteers scanned across 25 research sites were collected using GABA+ (GABA + co-edited macromolecules (MM)) and MM-suppressed GABA editing. The unsuppressed water signal from the volume of interest was acquired for concentration referencing. Whole-brain T1-weighted structural images were acquired and segmented to determine gray matter, white matter and cerebrospinal fluid voxel tissue fractions. Water-referenced GABA measurements were fully corrected for tissue-dependent signal relaxation and water visibility effects. The cohort-wide coefficient of variation was 17% for the GABA + data and 29% for the MM-suppressed GABA data. The mean within-site coefficient of variation was 10% for the GABA + data and 19% for the MM-suppressed GABA data. Vendor differences contributed 53% to the total variance in the GABA + data, while the remaining variance was attributed to site- (11%) and participant-level (36%) effects. For the MM-suppressed data, 54% of the variance was attributed to site differences, while the remaining 46% was attributed to participant differences. Results from an exploratory analysis suggested that the vendor differences were related to the unsuppressed water signal acquisition. Discounting the observed vendor-specific effects, water-referenced GABA measurements exhibit similar levels of variance to creatine-referenced GABA measurements. It is concluded that quantification using internal tissue water referencing is a viable and reliable method for the quantification of in vivo GABA levels. ispartof: NEUROIMAGE vol:191 pages:537-548 ispartof: location:United States status: published

Background Advances in sub-Nyquist-sampled dynamic contrast-enhanced (DCE) MRI enable monitoring of brain tumors with millimeter resolution and whole-brain coverage. Such undersampled quantitative methods need careful characterization regarding achievable test-retest reproducibility. Purpose To demonstrate a fully automated high-resolution whole-brain DCE MRI pipeline with 30-fold sparse undersampling and estimate its reproducibility on the basis of reference regions of stable tissue types during multiple posttreatment time points by using longitudinal clinical images of high-grade glioma. Materials and Methods Two methods for sub-Nyquist-sampled DCE MRI were extended with automatic estimation of vascular input functions. Continuously acquired three-dimensional k-space data with ramped-up flip angles were partitioned to yield high-resolution, whole-brain tracer kinetic parameter maps with matched precontrast-agent T1 and M

PURPOSE To demonstrate the feasibility of integrating the magnetization transfer (MT) preparations required for inhomogeneous MT (ihMT) within an MPRAGE-style acquisition. Such a sequence allows for reduced power deposition and easy inclusion of other modules. METHODS An ihMT MPRAGE-style sequence (ihMTRAGE) was initially simulated to investigate acquisition of the 3D ihMT data sequentially, or in an interleaved manner. The ihMTRAGE sequence was implemented on a 3T clinical scanner to acquire ihMT data from the brain and spine. RESULTS Both simulations and in vivo data provided an ihMT signal that was significantly greater using a sequential ihMTRAGE acquisition, compared with an interleaved implementation. Comparison with a steady-state ihMT acquisition (defined as having one MT RF pulse between successive acquisition modules) demonstrated how ihMTRAGE allows for a reduction in average power deposition, or greater ihMT signal at equal average power deposition. Inclusion of a prospective motion-correction module did not significantly affect the ihMT signal obtained from regions of interest in the brain. The ihMTRAGE acquisition allowed combination with a spatial saturation module to reduce phase wrap artifacts in a cervical spinal cord acquisition. CONCLUSIONS Use of preparations necessary for ihMT experiments within an MPRAGE-style sequence provides a useful alternative for acquiring 3D ihMT data. Compared with our steady-state implementation, ihMTRAGE provided reduced power deposition, while allowing use of the maximum intensity from off-resonance RF pulses. The 3D ihMTRAGE acquisition allowed combination of other modules with the preparation necessary for ihMT experiments, specifically motion compensation and spatial saturation modules.

Background Arterial Spin Labeling (ASL) MRI can provide quantitative images that are sensitive to both time averaged blood flow and its temporal fluctuations. 3D image acquisitions for ASL are desirable because they are more readily compatible with background suppression to reduce noise, can reduce signal loss and distortion, and provide uniform flow sensitivity across the brain. However, single-shot 3D acquisition for maximal temporal resolution typically involves degradation of image quality through blurring or noise amplification by parallel imaging. Here, we report a new approach to accelerate a common stack of spirals 3D image acquisition by pseudo golden-angle rotation and compressed sensing reconstruction without any degradation of time averaged blood flow images. Methods 28 healthy volunteers were imaged at 3T with background-suppressed unbalanced pseudo-continuous ASL combined with a pseudo golden-angle Stack-of-Spirals 3D RARE readout. A fully-sampled perfusion-weighted volume was reconstructed by 3D non-uniform Fast Fourier Transform (nuFFT) followed by sum-of-squares combination of the 32 individual channels. Coil sensitivities were estimated followed by reconstruction of the 39 single-shot volumes using an L1-wavelet Compressed-Sensing reconstruction. Finally, brain connectivity analyses were performed in regions where BOLD signal suffers from low signal-to-noise ratio and susceptibility artifacts. Results Image quality, assessed with a non-reference 3D blurring metric, of full time averaged blood flow was comparable to a conventional interleaved acquisition. The temporal resolution provided by the acceleration enabled identification and quantification of resting-state networks even in inferior regions such as the amygdala and inferior frontal lobes, where susceptibility artifacts can degrade conventional resting-state fMRI acquisitions. Conclusion This approach can provide measures of blood flow modulations and resting-state networks for free within any research or clinical protocol employing ASL for resting blood flow.

The neurophysiological effects of pediatric concussion several years after injury remain inadequately characterized. The objective of this study was to determine if a history of concussion was associated with BOLD response differences during an n-back working memory task in youth.Observational, cross-sectional.Participants include 52 children and adolescents (M = 15.1 years, 95%CI = 14.4-15.8, range = 9-19) with past concussion (n = 33) or orthopedic injury (OI; n = 19). Mean time since injury was 2.5 years (95%CI = 2.0-3.0). Measures included postconcussion symptom ratings, neuropsychological testing, and blood-oxygen-dependent-level (BOLD) functional magnetic resonance imaging (fMRI) during an n-back working memory task.Groups did not differ on accuracy or speed during the three n-back conditions. They also did not differ in BOLD signal change for the 1- vs. 0-back or 2- vs. 0-back contrasts (controlling for task performance).This study does not support group differences in BOLD response during an n-back working memory task in youth who are on average 2.5 years post-concussion. The findings are encouraging from the perspective of understanding recovery after pediatric concussion.

The role of white matter in reading has been established by diffusion tensor imaging (DTI), but DTI cannot identify specific microstructural features driving these relationships. Neurite orientation dispersion and density imaging (NODDI), inhomogeneous magnetization transfer (ihMT) and multicomponent driven equilibrium single-pulse observation of T1/T2 (mcDESPOT) can be used to link more specific aspects of white matter microstructure and reading due to their sensitivity to axonal packing and fiber coherence (NODDI) and myelin (ihMT and mcDESPOT). We applied principal component analysis (PCA) to combine DTI, NODDI, ihMT and mcDESPOT measures (10 in total), identify major features of white matter structure, and link these features to both reading and age. Analysis was performed for nine reading-related tracts in 46 neurotypical 6-16 year olds. We identified three principal components (PCs) which explained 79.5% of variance in our dataset. PC1 probed tissue complexity, PC2 described myelin and axonal packing, while PC3 was related to axonal diameter. Mixed effects regression models did not identify any significant relationships between principal components and reading skill. Bayes factor analysis revealed that the absence of relationships was not due to low power. Increasing PC1 in the left arcuate fasciculus with age suggest increases in tissue complexity, while increases of PC2 in the bilateral arcuate, inferior longitudinal, inferior fronto-occipital fasciculi, and splenium suggest increases in myelin and axonal packing with age. Multimodal white matter imaging and PCA provide microstructurally informative, powerful principal components which can be used by future studies of development and cognition. Our findings suggest major features of white matter undergo development during childhood and adolescence, but changes are not linked to reading during this period in our typically-developing sample.

Background The hardware and software differences between MR vendors and individual sites influence the quantification of MR spectroscopy data. An analysis of a large data set may help to better understand sources of the total variance in quantified metabolite levels. Purpose To compare multisite quantitative brain MR spectroscopy data acquired in healthy participants at 26 sites by using the vendor-supplied single-voxel point-resolved spectroscopy (PRESS) sequence. Materials and Methods An MR spectroscopy protocol to acquire short-echo-time PRESS data from the midparietal region of the brain was disseminated to 26 research sites operating 3.0-T MR scanners from three different vendors. In this prospective study, healthy participants were scanned between July 2016 and December 2017. Data were analyzed by using software with simulated basis sets customized for each vendor implementation. The proportion of total variance attributed to vendor-, site-, and participant-related effects was estimated by using a linear mixed-effects model. P values were derived through parametric bootstrapping of the linear mixed-effects models (denoted Pboot). Results In total, 296 participants (mean age, 26 years ± 4.6; 155 women and 141 men) were scanned. Good-quality data were recorded from all sites, as evidenced by a consistent linewidth of N-acetylaspartate (range, 4.4-5.0 Hz), signal-to-noise ratio (range, 174-289), and low Cramér-Rao lower bounds (≤5%) for all of the major metabolites. Among the major metabolites, no vendor effects were found for levels of myo-inositol (Pboot > .90), N-acetylaspartate and N-acetylaspartylglutamate (Pboot = .13), or glutamate and glutamine (Pboot = .11). Among the smaller resonances, no vendor effects were found for ascorbate (Pboot = .08), aspartate (Pboot > .90), glutathione (Pboot > .90), or lactate (Pboot = .28). Conclusion Multisite multivendor single-voxel MR spectroscopy studies performed at 3.0 T can yield results that are coherent across vendors, provided that vendor differences in pulse sequence implementation are accounted for in data analysis. However, the site-related effects on variability were more profound and suggest the need for further standardization of spectroscopic protocols. © RSNA, 2020 Online supplemental material is available for this article. ispartof: RADIOLOGY vol:295 issue:1 pages:171-180 ispartof: location:United States status: published

The long-term effects of concussion in youth remain poorly understood. The objective of this study was to determine the association between history of concussion and cerebral blood flow (CBF) in youth.A total of 53 children and adolescents with a history of concussion (n = 37) or orthopaedic injury (OI; n = 16) were considered. Measures included pseudo-continuous arterial spin labelling magnetic resonance imaging to quantify CBF, post-concussion symptoms, psychological symptoms, and cognitive testing.Participants (mean age: 14.4 years, 95% CI = 13.8-15.4, range = 8-19) were on average 2.7 years (95% CI = 2.2-3.1) post-injury. Youth with a history of concussion had higher parent-reported physical, cognitive, anxiety, and depression symptoms than children with OI, but the groups did not differ on self-reported symptoms (post-concussive or psychological) or cognitive testing. Global CBF did not differ between groups. Regional CBF analyses suggested that youth with a history of concussion had hypoperfusion in posterior and inferior regions and hyperperfusion in anterior/frontal/temporal regions as compared to those with OI. However, neither global nor regional CBF were significantly associated with demographics, pre-injury functioning, number of concussions, time since injury, post-concussive symptoms, psychological symptoms, or cognitive abilities.Youth with a history of concussion demonstrate differences in regional CBF (not global CBF), but without clear clinical expression.

This study aimed to evaluate the usefulness of transit time corrected cerebral blood flow (CBF) maps based on multi-phase arterial spin labeling MR perfusion imaging (ASL-MRP). The Institutional Review Board of our hospital approved this retrospective study. Written informed consent was waived. Conventional and multi-phase ASL-MRPs and dynamic susceptibility contrast MR perfusion imaging (DSC-MRP) were acquired for 108 consecutive patients. Vascular territory-based volumes of interest were applied to CBF and time to peak (TTP) maps obtained from DSC-MRP and CBF maps obtained from conventional and multi-phase ASL-MRPs. The concordances between normalized CBF (nCBF) from DSC-MRP and nCBF from conventional and transition time corrected CBF maps from multi-phase ASL-MRP were evaluated using Bland-Altman analysis. In addition, the dependence of difference between nCBF (ΔnCBF) values obtained from DSC-MRP and conventional ASL-MRP (or multi-phase ASL-MRP) on TTP obtained from DSC-MRP was also analyzed using regression analysis. The values of nCBFs from conventional and multi-phase ASL-MRPs had lower values than nCBF based on DSC-MRP (mean differences, 0.08 and 0.07, respectively). The values of ΔnCBF were dependent on TTP values from conventional ASL-MRP technique (F = 5.5679, P = 0.0384). No dependency of ΔnCBF on TTP values from multi-phase ASL-MRP technique was revealed (F = 0.1433, P > 0.05). The use of transit time corrected CBF maps based on multi-phase ASL-MRP technique can overcome the effect of delayed transit time on perfusion maps based on conventional ASL-MRP.

Purpose To develop and evaluate a model-based reconstruction framework for joint arterial input function (AIF) and kinetic parameter estimation from undersampled brain tumor dynamic contrast-enhanced MRI (DCE-MRI) data. Methods The proposed method poses the tracer-kinetic (TK) model as a model consistency constraint, enabling the flexible inclusion of different TK models and TK solvers, and the joint estimation of the AIF. The proposed method is evaluated using an anatomic realistic digital reference object (DRO), and nine retrospectively down-sampled brain tumor DCE-MRI datasets. We also demonstrate application to 30-fold prospectively undersampled brain tumor DCE-MRI. Results In DRO studies with up to 60-fold undersampling, the proposed method provided TK maps with low error that were comparable to fully sampled data and were demonstrated to be compatible with a third-party TK solver. In retrospective undersampling studies, this method provided patient-specific AIF with normalized root mean-squared-error (normalized by the 90th percentile value) less than 8% at up to 100-fold undersampling. In the 30-fold undersampled prospective study, the proposed method provided high-resolution whole-brain TK maps and patient-specific AIF. Conclusion The proposed model-based DCE-MRI reconstruction enables the use of different TK solvers with a model consistency constraint and enables joint estimation of patient-specific AIF. TK maps and patient-specific AIF with high fidelity can be reconstructed at up to 100-fold undersampling in k,t-space. Magn Reson Med, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Intraspinal microstimulation (ISMS) is currently under investigation for its ability to restore function following spinal cord injury and aid in addressing basic investigations of the spinal cord in feline and murine (rat) models. In this report we describe the procedures for fabricating and implanting intraspinal microwires, with special emphasis on the rat model. We also report our results on targeting success and long-term stability and functionality of the implants. Early targeting with implants fabricated based on general “average” dimensions of the spinal cord was approximately 50% successful in reaching the proper targets within the ventral grey matter in cats. Improvements in insertion technique and the use of multiple contact electrodes have raised the targeting success to 100%. Furthermore, the manufacturing of ISMS arrays has been improved by the use of magnetic resonance imaging to create subject-specific implants for cats and track the location of the arrays post-implant. In the rat, our procedures have produced desirable targeting of all recovered microwires. We speculate this is due to the different targeting parameters and the shorter depth of insertion in the rat spinal cord. Although there is a heightened mechanical mismatch between the $30~\mu \mathrm {m}$ -diameter microwires and the small rat spinal cord, chronic implantation and stimulation produce limited histological damage and do not compromise function. Furthermore, despite the increased difficulties of implanting into the smaller rat spinal cord, ISMS is effective in activating spinal cord networks in the lumbosacral enlargement in a manner that is safe, stable and reproducible.

In vivo quantification of glutamate (Glu) and γ-aminobutyric acid (GABA) using MRS is often achieved using two separate sequences: a short-echo point resolved spectroscopy (PRESS) acquisition for Glu and a Mescher-Garwood PRESS (MEGA-PRESS) acquisition for GABA. The purpose of this study was to examine the agreement of Glu and Glx (the combined signal of glutamate + glutamine) quantified from two different GABA-edited MEGA-PRESS acquisitions (GABA plus macromolecules, GABA+, TE = 68 ms, and macromolecule suppressed, MMSup, TE = 80 ms) with Glu and Glx quantified from a short-echo PRESS (PRESS-35, TE = 35 ms) acquisition. Fifteen healthy male volunteers underwent a single scan session, in which data were acquired using the three acquisitions (GABA+, MMSup and PRESS-35) in both the sensorimotor and anterior cingulate cortices using a voxel size of 3 × 3 × 3 cm3 . Glx and Glu were quantified from the MEGA-PRESS data using both the OFF sub-spectra and the difference (DIFF) spectra. Agreement was assessed using correlation analyses, Bland-Altman plots and intraclass correlation coefficients. Glx quantified from the OFF sub-spectra from both the GABA+ and MMSup acquisitions showed poor agreement with PRESS-35 in both brain regions. In the sensorimotor cortex, Glu quantified from the OFF sub-spectra of GABA+ showed moderate agreement with PRESS-35 data, but this finding was not replicated in the anterior cingulate cortex. Glx and Glu quantified using the DIFF spectra of either MEGA-PRESS sequence were in poor agreement with the PRESS-35 data in both brain regions. In conclusion, Glx and Glu measured from MEGA-PRESS data generally showed poor agreement with Glx and Glu measured using PRESS-35.

Arterial spin labeling (ASL) MRI can provide seizure onset zone (SOZ) localizing information in up to 80% of patients. Clinical implementation of this technique is limited by the need to obtain two scans per patient: a postictal scan that is subtracted from an interictal scan. We aimed to determine whether it is possible to limit the number of ASL scans to one per patient by comparing patient postictal ASL scans to baseline scans of 100 healthy controls. Eighteen patients aged 20-55 years underwent ASL MRI

Objectives We attempted to demonstrate the clinical applicability and utility of a three-dimensional multidelay arterial spin labeling magnetic resonance imaging technique in pediatric neuroimaging through a series of case studies. Methods Whole-brain three-dimensional multidelay arterial spin labeling data were acquired in five pediatric patients with different neurological conditions using 3 mm to 4 mm slices and a scan time of six to seven minutes. Results Three-dimensional multidelay arterial spin labeling provided complementary diagnostic information via quantitative cerebral blood flow and arterial transit time maps. Conclusions Three-dimensional multidelay arterial spin labeling sequence provides simultaneous quantification of cerebral blood flow and arterial transit time and is feasible for pediatric patients.

Adequate cerebral blood flow (CBF) is essential to proper brain development and function. Detailed characterization of CBF developmental trajectories will lead to better understanding of the development of cognitive, motor, and sensory functions, as well as behaviour in children. Previous studies have shown CBF increases during infancy and decreases during adolescence; however, the trajectories during childhood, and in particular the timing of peak CBF, remain unclear. Here, we used arterial spin labeling to map age-related changes of CBF across a large longitudinal sample that included 279 scans on 96 participants (46 girls and 50 boys) aged 2–7 years. CBF maps were analyzed using hierarchical linear regression for every voxel inside the grey matter mask, controlling for multiple comparisons. The results revealed a significant positive linear association between CBF and age in distributed brain regions including prefrontal, temporal, parietal, and occipital cortex, and in the cerebellum. There were no differences in developmental trajectories between males and females. Our findings show that CBF continues to increase until the age of 7 years, likely supporting ongoing improvements in behaviour, cognition, motor, and sensory functions in early childhood.

Adequate cerebral blood flow (CBF) is essential to proper brain development and function. Detailed characterization of CBF developmental trajectories will lead to better understanding of the development of cognitive, motor, and sensory functions, as well as behaviour in children. Previous studies have shown CBF increases during infancy and decreases during adolescence; however, the trajectories during childhood, and in particular the timing of peak CBF, remain unclear. Here, we used arterial spin labeling to map age-related changes of CBF across a large longitudinal sample that included 279 scans on 96 participants (46 girls and 50 boys) aged 2-7 years. CBF maps were analyzed using hierarchical linear regression for every voxel inside the grey matter mask, controlling for multiple comparisons. The results revealed a significant positive linear association between CBF and age in distributed brain regions including prefrontal, temporal, parietal, and occipital cortex, and in the cerebellum. There were no differences in developmental trajectories between males and females. Our findings suggest that CBF continues to increase at least until the age of 7 years, likely supporting the ongoing improvements in behaviour, cognition, motor, and sensory functions of the children in early childhood.Highlights*We mapped cerebral blood flow development longitudinally in 96 typically developing children aged 2-7 years*Cerebral blood flow increased over the full age range*No significant differences between males and females were foundConflicting of InterestR.M.L. is an employee of GE.

Background Smaller hippocampal volumes, as assessed by magnetic resonance imaging (MRI), and proton magnetic resonance spectroscopy ( 1 H-MRS) indexed alterations in brain metabolites have been identified in adults with major depressive disorder (MDD). Our group has found similar effects in MDD youth. However, this has not been studied in youth with treatment resistant MDD (TRD), nor has the interaction between regional N-acetyl-aspartate and volume deficits. N-acetyl-aspartate is an amino acid in the synthesis pathway of glutamate, and serves a marker of neuronal viability/number. Methods Fifteen typically developing youth (16–22 years of age; 7 males, 8 females) and eighteen youth with TRD (14–22 years of age; 8 males, 10 females) underwent 1 H-MRS and MRI on a 3 T scanner. A short echo PRESS protocol was used with voxels in the right and left hippocampi (6 mL each). Hippocampal volume was evaluated using FreeSurfer. Results Compared with the typically developing group, youth with TRD had lower concentrations of N-acetyl-aspartate in the left hippocampus (p=0.004), and a trend for smaller left hippocampal volume (p=0.067). In TRD subjects, hippocampal N-acetyl-aspartate was inversely correlated with left (r=−0.68, p=0.003) but not right hippocampal volume. Right hippocampal glutamate+glutamine was greater in TRD youth compared to typically developing controls (p=0.007). Conclusions These results suggest a neurochemical and structural deficit in the hippocampi of youth with TRD. These findings fit with the role of N-acetyl-aspartate in glutamate neurotransmission and the effect of glutamate on brain morphology.

Purpose The purpose of this work was to develop and evaluate a T1-weighted dynamic contrast enhanced (DCE) MRI methodology where tracer-kinetic (TK) parameter maps are directly estimated from undersampled (k,t)-space data. Theory and Methods The proposed reconstruction involves solving a nonlinear least squares optimization problem that includes explicit use of a full forward model to convert parameter maps to (k,t)-space, utilizing the Patlak TK model. The proposed scheme is compared against an indirect method that creates intermediate images by parallel imaging and compressed sensing before to TK modeling. Thirteen fully sampled brain tumor DCE-MRI scans with 5-second temporal resolution are retrospectively undersampled at rates R = 20, 40, 60, 80, and 100 for each dynamic frame. TK maps are quantitatively compared based on root mean-squared-error (rMSE) and Bland-Altman analysis. The approach is also applied to four prospectively R = 30 undersampled whole-brain DCE-MRI data sets. Results In the retrospective study, the proposed method performed statistically better than indirect method at R ≥ 80 for all 13 cases. This approach provided restoration of TK parameter values with less errors in tumor regions of interest, an improvement compared to a state-of-the-art indirect method. Applied prospectively, the proposed method provided whole-brain, high-resolution TK maps with good image quality. Conclusion Model-based direct estimation of TK maps from k,t-space DCE-MRI data is feasible and is compatible up to 100-fold undersampling. Magn Reson Med 78:1566–1578, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

To develop and evaluate a novel 3D Cartesian sampling scheme which is well suited for time-resolved 3D MRI using parallel imaging and compressed sensing.The proposed sampling scheme, termed GOlden-angle CArtesian Randomized Time-resolved (GOCART) 3D MRI, is based on golden angle (GA) Cartesian sampling, with random sampling of the ky-kz phase encode locations along each Cartesian radial spoke. This method was evaluated in conjunction with constrained reconstruction of retrospectively and prospectively undersampled in-vivo dynamic contrast enhanced (DCE) MRI data and simulated phantom data.In in-vivo retrospective studies and phantom simulations, images reconstructed from phase encodes defined by GOCART were equal to or superior to those with Poisson disc or GA sampling schemes. Typical GOCART sampling tables were generated in100ms. GOCART has also been successfully utilized prospectively to produce clinically valuable whole-brain DCE-MRI images.GOCART is a practical and efficient sampling scheme for time-resolved 3D MRI. It shows great potential for highly accelerated DCE-MRI and is well suited to modern reconstruction methods such as parallel imaging and compressed sensing.

Background: Deep gray matter (DGM) is affected in relapsing–remitting multiple sclerosis (RRMS) and may be studied using short-term longitudinal MRI. Objective: To investigate two-year changes in spin-echo transverse relaxation rate (R2) and atrophy in DGM, and its relationship with disease severity in RRMS patients. Methods: Twenty six RRMS patients and 26 matched controls were imaged at 4.7 T. Multiecho spin-echo R2 maps and atrophy measurements were obtained in DGM at baseline and two-year follow-up. Differences between MRI measures and correlations to disease severity were examined. Results: After two years, mean R2 values in the globus pallidus and pulvinar increased by ~4% ( p2 showed significant correlation to disease severity in the globus pallidus, pulvinar, substantia nigra, and thalamus. Multiple regression of the two-year R2 difference using these four DGM structures as variables, yielded high correlation with disease severity ( r=0.83, p2 showed significant correlation only for the globus pallidus in multiple sclerosis (MS) ( pConclusions: Two-year difference R2 measurements in DGM correlate to disease severity in MS. R2 mapping and atrophy measurements over two years can be used to identify changes in DGM in MS.

Aim Structural, functional, and metabolic changes in the dorsolateral prefrontal cortex (DLPFC) are implicated in the pathogenesis of major depressive disorder (MDD). We used proton magnetic resonance spectroscopy ((1) H-MRS) to examine the metabolite choline (glycerophosphocholine plus phosphocholine), which is used as an index of membrane integrity in the left DLPFC, in adolescents and young adults with MDD who were treatment-resistant and had a positive family history compared to healthy controls. Differences in the choline resonance indicate an imbalance between synthesis and degradation activity of neuronal and glia membrane phospholipids. Methods Seventeen adolescents with MDD and 11 healthy controls underwent (1) H-MRS. A short echo point-resolved spectroscopy (echo time = 30 ms, repetition time = 2000 ms) protocol was used with a voxel (4.5cm(3) , 128 averages) placed within the left DLPFC. Results There were significantly increased choline (P = 0.04) and creatine concentrations (P = 0.005) in the left DLPFC of the MDD group compared to controls. In MDD participants, choline concentration correlated with scores on the Beck Depression Inventory (r = 0.41, P = 0.03). Conclusion Increased left DLPFC choline and creatine levels in depressed adolescents may be biomarkers for the disorder. The increased choline levels may indicate abnormalities in neuronal membrane integrity, and the increased creatine could be reflective of altered energy demands and metabolism.

Purpose: To clinically evaluate a highly accelerated T1-weighted dynamic contrast-enhanced (DCE) MRI technique that provides high spatial resolution and whole-brain coverage via undersampling and constrained reconstruction with multiple sparsity constraints. Methods: Conventional (rate-2 SENSE) and experimental DCE-MRI (rate-30) scans were performed 20 minutes apart in 15 braintumor patients. The conventional clinical DCE-MRI had voxel dimensions 0.9 × 1.3 × 7.0 mm3, FOV 22 × 22 × 4.2 cm3, and the experimental DCE-MRI had voxel dimensions 0.9 × 0.9 × 1.9 mm3, and broader coverage 22 × 22 × 19 cm3. Temporal resolution was 5 s for both protocols. Time-resolvedimages and blood–brain barrier permeability maps were qualitatively evaluated by two radiologists. Results: The experimental DCE-MRI scans showed no loss of qualitative information in any of the cases, while achieving substantially higher spatial resolution and whole-brain spatial coverage. Average qualitative scores (from 0 to 3) were 2.1 for the experimental scans and 1.1 for the conventional clinical scans. Conclusions: The proposed DCE-MRI approach provides clinically superior image quality with higher spatial resolution and coverage than currently available approaches. These advantages may allow comprehensive permeability mapping in the brain, which is especially valuable in the setting of large lesions or multiple lesions spread throughout the brain.

Objective To construct a clinical diagnostic biomarker using state-of-the-art microstructural MRI in the motor cortex of people with amyotrophic lateral sclerosis (ALS). Methods Clinical and MRI data were obtained from 21 ALS patients (aged 54 ± 14 years, 33% female) and 63 age- and gender-matched controls (aged 48 ± 18 years, 43% female). MRI was acquired at 3T and included T1-weighted scan (for volumetrics), arterial spin labelling (for cerebral blood flow), susceptibility-weighted angiography (for iron deposition) and multiband diffusion kurtosis imaging (for tissue microstructure). Group differences in imaging measures in the motor cortex were tested by general linear model and relationships to clinical variables by linear regression. Results The ALS group had mild-to-moderate impairment (disease duration: 1.8 ± 0.8 years; ALS functional rating scale 40.2 ± 6.0; forced vital capacity 83% ± 22%). No age or gender differences were present between groups. We found significant group differences in diffusion kurtosis metrics (apparent, mean, radial and axial kurtosis: p, Highlights • A biomarker based on diffusion kurtosis imaging achieved an accuracy of 83%. • Kurtosis-based measures were more abnormal in ALS than tensor-based measures. • Motor cortex in the symptomatic hemisphere was smaller and had greater iron concentration. • There was a 1 mL volume loss per year in ALS motor cortex.

BACKGROUND AND PURPOSE:Multidelay arterial spin-labeling is a promising emerging method in clinical practice. The effect of imaging parameters in multidelay arterial spin-labeling on estimated cerebral blood flow measurements remains unknown. We directly compared 3-delay versus 7-delay sequences, assessing the difference in the estimated transit time and blood flow. MATERIALS AND METHODS:This study included 87 cognitively healthy controls (78.7 3.8 years of age; 49 women). We assessed delay and transit time-uncorrected and transit time-corrected CBF maps. Data analysis included voxelwise permutation-based between-sequence comparisons of 3-delay versus 7-delay, within-sequence comparison of transit time-uncorrected versus transit time-corrected maps, and average CBF calculations in regions that have been shown to differ. RESULTS:The 7-delay sequence estimated a higher CBF value than the 3-delay for the transit time-uncorrected and transit time-corrected maps in regions corresponding to the watershed areas (transit time-uncorrected = 27.62 +/- 12.23 versus 24.58 +/- 11.70 mL/min/100 g, Cohen's d = 0.25; transit time-corrected = 33.48 +/- 14.92 versus 30.16 +/- 14.32 mL/min/100 g, Cohen's d = 0.23). In the peripheral regions of the brain, the estimated delay was found to be longer for the 3-delay sequence (1.52408 +/- 0.25236 seconds versus 1.47755 +/- 0.24242 seconds, Cohen's d = 0.19), while the inverse was found in the center of the brain (1.39388 +/- 0.22056 seconds versus 1.42565 +/- 0.21872 seconds, Cohen's d = 0.14). Moreover, 7-delay had lower hemispheric asymmetry. CONCLUSIONS:The results of this study support the necessity of standardizing acquisition parameters in multidelay arterial spin-labeling and identifying basic parameters as a confounding factor in CBF quantification studies. Our findings conclude that multidelay arterial spin-labeling sequences with a high number of delays estimate higher CBF values than those with a lower number of delays.

Purpose We propose the use of a retrospectively gated cine fast spin echo (FSE) sequence for characterization of carotid artery dynamics. The aim of this study was to compare cine FSE measures of carotid dynamics with measures obtained on prospectively gated FSE images. Methods The common carotid arteries in 10 volunteers were imaged using two temporally resolved sequences: (i) cine FSE and (ii) prospectively gated FSE. Three raters manually traced a common carotid artery area for all cardiac phases on both sequences. Measured areas and systolic–diastolic area changes were calculated and compared. Inter- and intra-rater reliability were assessed for both sequences. Results No significant difference between cine FSE and prospectively gated FSE areas were observed (P = 0.36). Both sequences produced repeatable cross-sectional area measurements: inter-rater intraclass correlation coefficient (ICC) = 0.88 on cine FSE images and 0.87 on prospectively gated FSE images. Minimum detectable difference (MDD) in systolic–diastolic area was 4.9 mm2 with cine FSE and 6.4 mm2 with prospectively gated FSE. Conclusion This cine FSE method produced repeatable dynamic carotid artery measurements with less artifact and greater temporal efficiency compared with prospectively gated FSE. Magn Reson Med 74:1103–1109, 2015. © 2014 Wiley Periodicals, Inc.

A novel sequence has been introduced that combines multiband imaging with a multi-echo acquisition for simultaneous high spatial resolution pseudo-continuous arterial spin labeling (ASL) and blood-oxygenation-level dependent (BOLD) echo-planar imaging (MBME ASL/BOLD). Resting-state connectivity in healthy adult subjects was assessed using this sequence. Four echoes were acquired with a multiband acceleration of four, in order to increase spatial resolution, shorten repetition time, and reduce slice-timing effects on the ASL signal. In addition, by acquiring four echoes, advanced multi-echo independent component analysis (ME-ICA) denoising could be employed to increase the signal-to-noise ratio (SNR) and BOLD sensitivity. Seed-based and dual-regression approaches were utilized to analyze functional connectivity. Cerebral blood flow (CBF) and BOLD coupling was also evaluated by correlating the perfusion-weighted timeseries with the BOLD timeseries. These metrics were compared between single echo (E2), multi-echo combined (MEC), multi-echo combined and denoised (MECDN), and perfusion-weighted (PW) timeseries. Temporal SNR increased for the MECDN data compared to the MEC and E2 data. Connectivity also increased, in terms of correlation strength and network size, for the MECDN compared to the MEC and E2 datasets. CBF and BOLD coupling was increased in major resting-state networks, and that correlation was strongest for the MECDN datasets. These results indicate our novel MBME ASL/BOLD sequence, which collects simultaneous high-resolution ASL/BOLD data, could be a powerful tool for detecting functional connectivity and dynamic neurovascular coupling during the resting state. The collection of more than two echoes facilitates the use of ME-ICA denoising to greatly improve the quality of resting state functional connectivity MRI.

Magnetic resonance spectroscopy (MRS) is the only biomedical imaging method that can noninvasively detect endogenous signals from the neurotransmitter γ-aminobutyric acid (GABA) in the human brain. Its increasing popularity has been aided by improvements in scanner hardware and acquisition methodology, as well as by broader access to pulse sequences that can selectively detect GABA, in particular J-difference spectral editing sequences. Nevertheless, implementations of GABA-edited MRS remain diverse across research sites, making comparisons between studies challenging. This large-scale multi-vendor, multi-site study seeks to better understand the factors that impact measurement outcomes of GABA-edited MRS. An international consortium of 24 research sites was formed. Data from 272 healthy adults were acquired on scanners from the three major MRI vendors and analyzed using the Gannet processing pipeline. MRS data were acquired in the medial parietal lobe with standard GABA+ and macromolecule- (MM-) suppressed GABA editing. The coefficient of variation across the entire cohort was 12% for GABA+ measurements and 28% for MM-suppressed GABA measurements. A multilevel analysis revealed that most of the variance (72%) in the GABA+ data was accounted for by differences between participants within-site, while site-level differences accounted for comparatively more variance (20%) than vendor-level differences (8%). For MM-suppressed GABA data, the variance was distributed equally between site- (50%) and participant-level (50%) differences. The findings show that GABA+ measurements exhibit strong agreement when implemented with a standard protocol. There is, however, increased variability for MM-suppressed GABA measurements that is attributed in part to differences in site-to-site data acquisition. This study's protocol establishes a framework for future methodological standardization of GABA-edited MRS, while the results provide valuable benchmarks for the MRS community. This work was supported by NIH grants R01 EB016089, R01 EB023963 and P41 EB015909. Data collection was supported by the Shandong Provincial Key Research and Development Plan of China (2016ZDJS07A16) and the National Natural Science Foundation of China for Young Scholars (no. 81601479). IDW thanks Mrs. J. Bigley of the University of Sheffield MRI Unit for her assistance with data acquisition. JJP was supported by NIAAA grant K23 AA020842. MPS was supported by NIH grant F32 EY025121. NAJP receives salary support from NIH grant K99 MH107719. The authors acknowledge implementation contributions from a number of employees of Siemens Medical Solutions, including Dr. Keith Heberlein and Dr. Sinyeob Ahn, to the Siemens WIP sequences, which are shared with several research sites under sequencespecific agreements.

Transverse relaxation (T2) mapping has many applications, including imaging of iron accumulation in grey matter. Using the typical multiecho spin-echo sequence with long echo trains, stimulated echo compensation can enable T2 fitting under conditions of variable radio frequency homogeneity arising from slice profile and in-plane radio frequency variation. Substantial reduction in the number of refocusing pulses could enable use at high magnetic fields where specific absorption rate is a major limitation, and enable multislice use with reduced incidental magnetization transfer at all field strengths. We examine the effect of reduced echo train lengths and multislice imaging on T2 fitting using stimulated echo compensation applied to iron-rich subcortical grey matter in human brain at 4.7 T. Our findings indicate that reducing from 20 echoes to as few as four echoes can maintain consistent T2 values when using stimulated echo compensation in grey and white matter, but not for cerebrospinal fluid. All territories produce marginal results when using standard exponential fitting. Savings from reduced echoes can be used to substantially increase slice coverage. In multislice mode, the resulting incidental magnetization transfer decreased brain signal but had minimal effect on measured T2 values. Magn Reson Med 70:1340–1346, 2013. © 2013 Wiley Periodicals, Inc.

Purpose To investigate the brain iron dependence of transverse relaxation rate difference methods derived from spin echo and gradient echo measurements from two field strengths. Methods Transverse relaxation rates R 2 and R 2 ⁎ were measured in human brain in 17 healthy subjects at 1.5 T and 4.7 T using multi-slice, multiecho spin echo and gradient echo sequences. R 2 quantification used stimulated echo compensation and R 2 ⁎ quantification used linear background gradient correction at 4.7 T only. Subtraction of R 2 from R 2 ⁎ within each field strength yielded R 2 ’ , and R 2 subtraction across fields yielded Field Dependent R 2 Increase (FDRI). All transverse relaxation measures were then correlated with published post-mortem iron concentrations using linear regression analysis. Regional differences were tested using paired t-tests. Phantom measurements of FDRI were also performed. Results In deep grey matter, all transverse relaxation rates (R 2 , R 2 ⁎ , R 2 ’ ) at both 1.5 T and 4.7 T, and FDRI had moderate to strong correlations (r > 0.71, p 0.9) than 1.5 T measures. R 2 , R 2 ⁎ , R 2 ’ at 4.7 T and FDRI had slopes 0.49, 1.96, 1.48 and 0.33 [s − 1 /mg Fe/100 g wt. tissue] and intercepts 14.40, 16.87, 2.47 and 3.21 [s − 1 ] respectively. Even though FDRI yielded a zero intercept in phantom, in vivo FDRI was found to be ineffective at fully removing non-iron contributions and yielded a large intercept. The slope for R 2 ’ was 3.4 times greater at 4.7 T than 1.5 T. For white matter fiber tracts oriented predominantly perpendicular versus parallel to B 0 , R 2 ’ increased by ~ 50% at 4.7 T and ~ 30% at 1.5 T, while R 2 and FDRI in white matter was insensitive to its orientation with respect to B 0 . Conclusion The transverse relaxation difference methods FDRI and R 2 ’ at 4.7 T had high correlations to predicted iron content similar to R 2 and R 2 ⁎ at 4.7 T. Although R 2 ’ had smaller y-intercept with estimated iron concentration than FDRI, in white matter R 2 ’ demonstrated strongest dependence on fiber orientation with respect to B 0. These results suggest that for brain iron correlation, there is minimal value of transverse difference methods over a single R 2 ⁎ measurement at highest available field, which was 4.7 T.