Routine clinical brain MRI sequences for use at 3.0 Tesla

Black-Blood Magnetization Prepared 2 Rapid Acquisition Gradient Echoes: A Fast and Three-Dimensional MR Black-Blood T1 Mapping Technique for Quantitative Assessment of Atherosclerosis and Venous Thrombosis

BACKGROUND

Blood flow signals may be a confounder in quantifying T1 values of plaque or thrombus and how to realize black-blood T1 mapping remains a challenge task.

PURPOSE

To develop a fast and three-dimensional black-blood T1 mapping technique for quantitative assessment of atherosclerosis and venous thrombosis.

STUDY TYPE

Sequence development and optimization via phantoms and volunteers as well as pilot prospective.

PHANTOM AND SUBJECTS

Numerical simulations, a standard phantom, 8 healthy volunteers (mean age, 22 ± 1 years; 5 males), and 19 patients (mean age, 57 ± 14 years; 13 males) with atherosclerosis or venous thrombosis.

FIELD STRENGTH/SEQUENCE

3T/inversion recovery spin-echo sequence (IR-SE), magnetization prepared 2 rapid acquisition gradient echoes (MP2RAGE), and black-blood prepared MP2RAGE (BB-MP2RAGE).

ASSESSMENT

The black-blood preparation (i.e., delay alternating with nutation for tailored excitation, DANTE) was incorporated into MP2RAGE for black-blood T1 mapping. The BB-MP2RAGE was optimized numerically based on the Bloch equation, and then the phantom study was performed to verify the accuracy of T1 mapping by BB-MP2RAGE against IR-SE and MP2RAGE. Preliminary clinical validation was prospectively performed to assess the flow suppression effect and its potential application in plaque and thrombosis identification.

STATISTICAL TESTS

Pearson correlation test, Bland-Altman analysis, paired t-test, and intraclass correlation coefficient. A P value <0.05 indicates a statistically significant difference.

RESULTS

Phantom experiments showed comparable accuracy of T1 maps by BB-MP2RAGE with IR-SE and MP2RAGE (all r2 > 0.99); Compared to MP2RAGE, BB-MP2RAGE effectively nulled the blood flow signals, and had a significant improvement in contrast-to-noise ratio between static tissue and blood (250.5 ± 66.6 vs. 91.9 ± 35.9). BB-MP2RAGE can quantify plaque or thrombus T1 relaxation time with blood flow signal suppression.

DATA CONCLUSION

Accurate T1 mapping with sufficient blood flow suppression was achieved by BB-MP2RAGE. BB-MP2RAGE has the potential to quantitatively characterize atherosclerosis and venous thrombosis.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.

Retrospective analysis of Braak stage- and APOE4 allele-dependent associations between MR spectroscopy and markers of tau and neurodegeneration in cognitively unimpaired elderly

PURPOSE

The pathological hallmarks of Alzheimer's disease (AD), amyloid, tau, and associated neurodegeneration, are present in the cortical gray matter (GM) years before symptom onset, and at significantly greater levels in carriers of the apolipoprotein E4 (APOE4) allele. Their respective biomarkers, A/T/N, have been found to correlate with aspects of brain biochemistry, measured with magnetic resonance spectroscopy (MRS), indicating a potential for MRS to augment the A/T/N framework for staging and prediction of AD. Unfortunately, the relationships between MRS and A/T/N biomarkers are unclear, largely due to a lack of studies examining them in the context of the spatial and temporal model of T/N progression. Advanced MRS acquisition and post-processing approaches have enabled us to address this knowledge gap and test the hypotheses, that glutamate-plus-glutamine (Glx) and N-acetyl-aspartate (NAA), metabolites reflecting synaptic and neuronal health, respectively, measured from regions on the Braak stage continuum, correlate with: (i) cerebrospinal fluid (CSF) p-tau181 level (T), and (ii) hippocampal volume or cortical thickness of parietal lobe GM (N). We hypothesized that these correlations will be moderated by Braak stage and APOE4 genotype.

METHODS

We conducted a retrospective imaging study of 34 cognitively unimpaired elderly individuals who received APOE4 genotyping and lumbar puncture from pre-existing prospective studies at the NYU Grossman School of Medicine between October 2014 and January 2019. Subjects returned for their imaging exam between April 2018 and February 2020. Metabolites were measured from the left hippocampus (Braak II) using a single-voxel semi-adiabatic localization by adiabatic selective refocusing sequence; and from the bilateral posterior cingulate cortex (PCC; Braak IV), bilateral precuneus (Braak V), and bilateral precentral gyrus (Braak VI) using a multi-voxel echo-planar spectroscopic imaging sequence. Pearson and Spearman correlations were used to examine the relationships between absolute levels of choline, creatine, myo-inositol, Glx, and NAA and CSF p-tau181, and between these metabolites and hippocampal volume or parietal cortical thicknesses. Covariates included age, sex, years of education, Fazekas score, and months between CSF collection and MRI exam.

RESULTS

There was a direct correlation between hippocampal Glx and CSF p-tau181 in APOE4 carriers (Pearson's r = 0.76, p = 0.02), but not after adjusting for covariates. In the entire cohort, there was a direct correlation between hippocampal NAA and hippocampal volume (Spearman's r = 0.55, p = 0.001), even after adjusting for age and Fazekas score (Spearman's r = 0.48, p = 0.006). This relationship was observed only in APOE4 carriers (Pearson's r = 0.66, p = 0.017), and was also retained after adjustment (Pearson's r = 0.76, p = 0.008; metabolite-by-carrier interaction p = 0.03). There were no findings in the PCC, nor in the negative control (late Braak stage) regions of the precuneus and precentral gyrus.

CONCLUSIONS

Our findings are in line with the spatially- and temporally-resolved Braak staging model of pathological severity in which the hippocampus is affected earlier than the PCC. The correlations, between MRS markers of synaptic and neuronal health and, respectively, T and N pathology, were found exclusively within APOE4 carriers, suggesting a connection with AD pathological change, rather than with normal aging. We therefore conclude that MRS has the potential to augment early A/T/N staging, with the hippocampus serving as a more sensitive MRS target compared to the PCC.

Compartmental anisotropy of filtered exchange imaging (FEXI) in human white matter: What is happening in FEXI?

PURPOSE

To investigate the effects of compartmental anisotropy on filtered exchange imaging (FEXI) in white matter (WM).

THEORY AND METHODS

FEXI signals were measured using multiple combinations of diffusion filter and detection directions in five healthy volunteers. Additional filters, including a trace-weighted diffusion filter with trapezoidal gradients, a spherical b-tensor encoded diffusion filter, and a T2 filter, were tested with trace-weighted diffusion detection.

RESULTS

A large range of apparent exchange rates (AXR) and both positive and negative filter efficiencies (σ) were found depending on the mutual orientation of the filter and detection gradients relative to WM fiber orientation. The data demonstrated that the fast-diffusion compartment suppressed by diffusional filtering is not exclusively extra-cellular, but also intra-cellular. While not comprehensive, a simple two-compartment diffusion tensor model with water exchange was able to account qualitatively for the trends in positive and negative filtering efficiencies, while standard model imaging (SMI) without exchange could not. This two-compartment diffusion tensor model also demonstrated smaller AXR variances across subjects. When employing trace-weighted diffusion detection, AXR values were on the order of the R1 (=1/T1) of water at 3T for crossing fibers, while being less than R1 for parallel fibers.

CONCLUSION

Orientation-dependent AXR and σ values were observed when using multi-orientation filter and detection gradients in FEXI, indicating that WM FEXI models need to account for compartmental anisotropy. When using trace-weighted detection, AXR values were on the order of or less than R1, complicating the interpretation of FEXI results in WM in terms of biological exchange properties. These findings may contribute toward better understanding of FEXI results in WM.

VICTR: Venous transit time imaging by changes in T1 relaxation

PURPOSE

Abnormalities in cerebral veins are a common finding in many neurological diseases, yet there is a scarcity of MRI techniques to assess venous hemodynamic function. The present study aims to develop a noncontrast technique to measure a novel blood flow circulatory measure, venous transit time (VTT), which denotes the time it takes for water to travel from capillary to major veins.

METHODS

The proposed sequence, venous transit time imaging by changes in T1 relaxation (VICTR), is based on the notion that as water molecules transition from the tissue into the veins, they undergo a change in T1 relaxation time. The validity of the measured VTT was tested by studying the VTT along the anatomically known flow trajectory of venous vessels as well as using a physiological vasoconstrictive challenge of caffeine ingestion. Finally, we compared the VTT measured with VICTR MRI to a bolus-tracking method using gadolinium-based contrast agent.

RESULTS

VTT was measured to be 3116.3 ± 326.0 ms in the posterior superior sagittal sinus (SSS), which was significantly longer than 2865.0 ± 390.8 ms at the anterior superior sagittal sinus (p = 0.004). The test-retest assessment showed an interclass correlation coefficient of 0.964. VTT was significantly increased by 513.8 ± 239.3 ms after caffeine ingestion (p < 0.001). VTT measured with VICTR MRI revealed a strong correlation (R = 0.84, p = 0.002) with that measured with the contrast-based approach. VTT was found inversely correlated to cerebral blood flow and venous oxygenation across individuals.

CONCLUSION

A noncontrast MRI technique, VICTR MRI, was developed to measure the VTT of the brain.

Association between hearing ability and cortical morphology in the elderly: multiparametric mapping, cognitive relevance, and neurobiological underpinnings

BACKGROUND

Hearing impairment is a common condition in the elderly. However, a comprehensive understanding of its neural correlates is still lacking.

METHODS

We recruited 284 elderly adults who underwent structural MRI, magnetic resonance spectroscopy, audiometry, and cognitive assessments. Individual hearing abilities indexed by pure tone average (PTA) were correlated with multiple structural MRI-derived cortical morphological indices. For regions showing significant correlations, mediation analyses were performed to examine their role in the relationship between hearing ability and cognitive function. Finally, the correlation maps between hearing ability and cortical morphology were linked with publicly available connectomic gradient, transcriptomic, and neurotransmitter maps.

FINDINGS

Poorer hearing was related to cortical thickness (CT) reductions in widespread regions and gyrification index (GI) reductions in the right Area 52 and Insular Granular Complex. The GI in the right Area 52 mediated the relationship between hearing ability and executive function. This mediating effect was further modulated by glutamate and N-acetylaspartate levels in the right auditory region. The PTA-CT correlation map followed microstructural connectomic hierarchy, were related to genes involved in certain biological processes (e.g., glutamate metabolic process), cell types (e.g., excitatory neurons and astrocytes), and developmental stages (i.e., childhood to young adulthood), and covaried with dopamine receptor 1, dopamine transporter, and fluorodopa. The PTA-GI correlation map was related to 5-hydroxytryptamine receptor 2a.

INTERPRETATION

Poorer hearing is associated with cortical thinning and folding reductions, which may be engaged in the relationship between hearing impairment and cognitive decline in the elderly and have different neurobiological substrates.

FUNDING

See the Acknowledgements section.

Brain extended and closed forms glutathione levels decrease with age and extended glutathione is associated with visuospatial memory

During aging, the brain is subject to greater oxidative stress (OS), which is thought to play a critical role in cognitive impairment. Glutathione (GSH), as a major antioxidant in the brain, can be used to combat OS. However, how brain GSH levels vary with age and their associations with cognitive function is unclear. In this study, we combined point-resolved spectroscopy and edited spectroscopy sequences to investigate extended and closed forms GSH levels in the anterior cingulate cortex (ACC), posterior cingulate cortex (PCC), and occipital cortex (OC) of 276 healthy participants (extended form, 166 females, age range 20-70 years) and 15 healthy participants (closed form, 7 females, age range 26-56 years), and examined their relationships with age and cognitive function. The results revealed decreased extended form GSH levels with age in the PCC among 276 participants. Notably, the timecourse of extended form GSH level changes in the PCC and ACC differed between males and females. Additionally, positive correlations were observed between extended form GSH levels in the PCC and OC and visuospatial memory. Additionally, a decreased trend of closed form GSH levels with age was also observed in the PCC among 15 participants. Taken together, these findings enhance our understanding of the brain both closed and extended form GSH time course during normal aging and associations with sex and memory, which is an essential first step for understanding the neurochemical underpinnings of healthy aging.

Reduced cerebrospinal fluid motion in patients with Parkinson's disease revealed by magnetic resonance imaging with low b-value diffusion weighted imaging

BACKGROUND

Parkinson's disease is characterized by dopamine-responsive symptoms as well as aggregation of α-synuclein protofibrils. New diagnostic methods assess α-synuclein aggregation characteristics from cerebrospinal fluid (CSF) and recent pathophysiologic mechanisms suggest that CSF circulation disruptions may precipitate α-synuclein retention. Here, diffusion-weighted MRI with low-to-intermediate diffusion-weightings was applied to test the hypothesis that CSF motion is reduced in Parkinson's disease relative to healthy participants.

METHODS

Multi-shell diffusion weighted MRI (spatial resolution = 1.8 × 1.8 × 4.0 mm) with low-to-intermediate diffusion weightings (b-values = 0, 50, 100, 200, 300, 700, and 1000 s/mm2) was applied over the approximate kinetic range of suprasellar cistern fluid motion at 3 Tesla in Parkinson's disease (n = 27; age = 66 ± 6.7 years) and non-Parkinson's control (n = 32; age = 68 ± 8.9 years) participants. Wilcoxon rank-sum tests were applied to test the primary hypothesis that the noise floor-corrected decay rate of CSF signal as a function of b-value, which reflects increasing fluid motion, is reduced within the suprasellar cistern of persons with versus without Parkinson's disease and inversely relates to choroid plexus activity assessed from perfusion-weighted MRI (significance-criteria: p < 0.05).

RESULTS

Consistent with the primary hypothesis, CSF decay rates were higher in healthy (D = 0.00673 ± 0.00213 mm2/s) relative to Parkinson's disease (D = 0.00517 ± 0.00110 mm2/s) participants. This finding was preserved after controlling for age and sex and was observed in the posterior region of the suprasellar cistern (p < 0.001). An inverse correlation between choroid plexus perfusion and decay rate in the voxels within the suprasellar cistern (Spearman's-r=-0.312; p = 0.019) was observed.

CONCLUSIONS

Multi-shell diffusion MRI was applied to identify reduced CSF motion at the level of the suprasellar cistern in adults with versus without Parkinson's disease; the strengths and limitations of this methodology are discussed in the context of the growing literature on CSF flow.

Uncertainty propagation in absolute metabolite quantification for in vivo MRS of the human brain

PURPOSE

Absolute spectral quantification is the standard method for deriving estimates of the concentration from metabolite signals measured using in vivo proton MRS (1 H-MRS). This method is often reported with minimum variance estimators, specifically the Cramér-Rao lower bound (CRLB) of the metabolite signal amplitude's scaling factor from linear combination modeling. This value serves as a proxy for SD and is commonly reported in MRS experiments. Characterizing the uncertainty of absolute quantification, however, depends on more than simply the CRLB. The uncertainties of metabolite-specific (T1m , T2m ), reference-specific (T1ref , T2ref ), and sequence-specific (TR , TE ) parameters are generally ignored, potentially leading to an overestimation of precision. In this study, the propagation of uncertainty is used to derive a comprehensive estimate of the overall precision of concentrations from an internal reference.

METHODS

The propagated uncertainty is calculated using analytical derivations and Monte Carlo simulations and subsequently analyzed across a set of commonly measured metabolites and macromolecules. The effect of measurement error from experimentally obtained quantification parameters is estimated using published uncertainties and CRLBs from in vivo 1 H-MRS literature.

RESULTS

The additive effect of propagated measurement uncertainty from applied quantification correction factors can result in up to a fourfold increase in the concentration estimate's coefficient of variation compared to the CRLB alone. A case study analysis reveals similar multifold increases across both metabolites and macromolecules.

CONCLUSION

The precision of absolute metabolite concentrations derived from 1 H-MRS experiments is systematically overestimated if the uncertainties of commonly applied corrections are neglected as sources of error.

Investigations on correlations between changes of optical tissue properties and NMR relaxation times

BACKGROUND

Accurate light dosimetry is a complex remaining challenge in interstitial photodynamic therapy (iPDT) for malignant gliomas. The light dosimetry should ideally be based on the tissue morphology and the individual optical tissue properties of each tissue type in the target region. First investigations are reported on using NMR information to estimate changes of individual optical tissue properties.

METHODS

Porcine brain tissue and optical tissue phantoms were investigated. To the porcine brain, supplements were added to simulate an edema or high blood content. The tissue phantoms were based on agar, Lipoveneous, ink, blood and gadobutrol (Gd-based MRI contrast agent). The concentrations of phantom ingredients and tissue additives are varied to compare concentration-dependent effects on optical and NMR properties. A 3-tesla whole-body MRI system was used to determine T1 and T2 relaxation times. Optical tissue properties, i.e., the spectrally resolved absorption and reduced scattering coefficient, were obtained using a single integrating sphere setup. The observed changes of NMR and optical properties were compared to each other.

RESULTS

By adjusting the NMR relaxation times and optical tissue properties of the tissue phantoms to literature values, recipes for human brain tumor, white matter and grey matter tissue phantoms were obtained that mimic these brain tissues simultaneously in both properties. For porcine brain tissue, it was observed that with increasing water concentration in the tissue, both NMR-relaxation times increased, while µa decreased and µs' increased at 635 nm. The addition of blood to porcine brain samples showed a constant T1, while T2 shortened and the absorption coefficient at 635 nm increased.

CONCLUSIONS

In this investigation, by changing sample contents, notable changes of both NMR relaxation times and optical tissue properties have been observed and their relations examined. The developed dual NMR/optical tissue phantoms can be used in iPDT research, clinical training and demonstrations.

Technical Note: Novel imaging method to obtain gray matter-attenuated inversion recovery image using low-field magnetic resonance imaging systems

INTRODUCTION

The double inversion recovery (DIR) technique suppresses two types of tissue signals with different T1 values by applying two inversion recovery (IR) pulses with different inversion times (TI). In contrast, the double tissue suppression with multi-echo acquisition and single TI combining HIRE (DOMUST-HIRE) method, is a technique enabling the white-matter-attenuated inversion recovery (WAIR) images by setting one inversion time (TI) in a sequence based on the multi-echo method and subtracting the second echo image from the first echo image. Here, we propose a new sequence that can provide the gray-matter-attenuated inversion recovery image based on the DOMUST-HIRE method.

METHODS

In this small clinical study, we performed determination of optimal TI and physical evaluation by imaging a subject's head with T1WI and our proposed method for GAIR images.

RESULTS

Our proposed method could increase the contrast ratio and the contrast-to-noise ratio between white matter (WM) and gray matter (GM), whereas the signal-to-noise ratio WM and GM decreased than with T1WI method.

CONCLUSIONS

Our proposed method can be used to suppress GM and CSF signals.

IMPLICATIONS FOR PRACTICE

The use of our proposed method in low-field MRI systems could provide GAIR image.

Optimal control in a magnetization-prepared rapid acquisition gradient-echo sequence

This article proposes a numerical framework to determine the optimal magnetization preparation in a three-dimensional magnetization-prepared rapid gradient-echo (MP-RAGE) sequence to obtain the best achievable contrast between target tissues based on differences in their relaxation times. The benefit lies in the adaptation of the algorithm of optimal control, GRAdient Ascent Pulse Engineering (GRAPE), to the optimization of magnetization preparation in a cyclic sequence without full recovery between each cycle. This numerical approach optimizes magnetization preparation of an arbitrary number of radio frequency pulses to enhance contrast, taking into account the establishment of a steady state in the longitudinal component of the magnetization. The optimal control preparation offers an optimized mixed T 1 / T 2 contrast in this traditional T 1 -weighted sequence. To show the versatility of the proposed method, numerical and in vitro results are described. Examples of contrasts acquired on brain regions of a healthy volunteer are presented for potential applications at 3 T.

Quantification of T1 and T2 of subarachnoid CSF: Implications for water exchange between CSF and brain tissues

PURPOSE

To quantify the T1 and T2 values of CSF in the subarachnoid space (SAS) at 3 T and interpret them in the context of water exchange between CSF and brain tissues.

METHODS

CSF T1 was measured using inversion recovery, and CSF T2 was assessed using T2 -preparation. T1 and T2 values in the SAS were compared with those in the frontal horns of lateral ventricles, which have less brain-CSF exchange. Phantom experiments were performed to examine whether there were spatial variations in T1 and T2 that were unrelated to brain-CSF exchange. Simulations were conducted to investigate the relationship between the brain-CSF exchange rate and the apparent T1 and T2 values of SAS CSF.

RESULTS

The CSF T1 and T2 values were 4308.7 ± 146.9 ms and 1885.5 ± 67.9 ms, respectively, in the SAS and were 4454.0 ± 187.9 ms and 2372.9 ± 72.0 ms in the frontal horns. The SAS CSF had shorter T1 (p = 0.006) and T2 (p < 0.0001) than CSF in the frontal horns. Phantom experiments showed negligible (< 6 ms for T1 ; < 1 ms for T2 ) spatial variations in T1 and T2 , suggesting that the T1 and T2 differences between SAS and frontal horns were largely attributed to physiological reasons. Simulations revealed that faster brain-CSF exchange rates lead to shorter apparent T1 and T2 of SAS CSF. However, the experimentally observed T2 difference between SAS and frontal horns was greater than that attributable to typical exchange effect, suggesting that the T2 shortening in SAS may reflect a combined effect of exchange and deoxyhemoglobin susceptibility.

CONCLUSION

Quantification of SAS CSF relaxation times may be useful to assess the brain-CSF exchange.

Data-driven electrical conductivity brain imaging using 3 T MRI

Magnetic resonance electrical properties tomography (MR-EPT) is a non-invasive measurement technique that derives the electrical properties (EPs, e.g., conductivity or permittivity) of tissues in the radiofrequency range (64 MHz for 1.5 T and 128 MHz for 3 T MR systems). Clinical studies have shown the potential of tissue conductivity as a biomarker. To date, model-based conductivity reconstructions rely on numerical assumptions and approximations, leading to inaccuracies in the reconstructed maps. To address such limitations, we propose an artificial neural network (ANN)-based non-linear conductivity estimator trained on simulated data for conductivity brain imaging. Network training was performed on 201 synthesized T2-weighted spin-echo (SE) data obtained from the finite-difference time-domain (FDTD) electromagnetic (EM) simulation. The dataset was composed of an approximated T2-w SE magnitude and transceive phase information. The proposed method was tested three in-silico and in-vivo on two volunteers and three patients' data. For comparison purposes, various conventional phase-based EPT reconstruction methods were used that ignoreB 1 + magnitude information, such as Savitzky-Golay kernel combined with Gaussian filter (S-G Kernel), phase-based convection-reaction EPT (cr-EPT), magnitude-weighted polynomial-fitting phase-based EPT (Poly-Fit), and integral-based phase-based EPT (Integral-based). From the in-silico experiments, quantitative analysis showed that the proposed method provides more accurate and improved quality (e.g., high structural preservation) conductivity maps compared to conventional reconstruction methods. Representatively, in the healthy brain in-silico phantom experiment, the proposed method yielded mean conductivity values of 1.97 ± 0.20 S/m for CSF, 0.33 ± 0.04 S/m for WM, and 0.52 ± 0.08 S/m for GM, which were closer to the ground-truth conductivity (2.00, 0.30, 0.50 S/m) than the integral-based method (2.56 ± 2.31, 0.39 ± 0.12, 0.68 ± 0.33 S/m). In-vivo ANN-based conductivity reconstructions were also of improved quality compared to conventional reconstructions and demonstrated network generalizability and robustness to in-vivo data and pathologies. The reported in-vivo brain conductivity values were in agreement with literatures. In addition, the proposed method was observed for various SNR levels (SNR levels = 10, 20, 40, and 58) and repeatability conditions (the eight acquisitions with the number of signal averages = 1). The preliminary investigations on brain tumor patient datasets suggest that the network trained on simulated dataset can generalize to unforeseen in-vivo pathologies, thus demonstrating its potential for clinical applications.

Advanced software for MRgFUS treatment planning

BACKGROUND AND OBJECTIVES

Herein, a user-friendly software platform for 3-dimensional Focused Ultrasound treatment planning based on Magnetic Resonance Imaging (MRI) images is presented.

METHODS

The software directly retrieves and loads MRI images. Various design tools can be used on the MRI images to define the treatment area and the sonication parameters. Based on the treatment plan, the software controls the robotic motion and motion pattern of Magnetic Resonance guided Focused Ultrasound (MRgFUS) robotic systems to execute the treatment procedure. Real-time treatment monitoring is achieved through MRI images and thermometry. The software's functionality and performance were evaluated in both laboratory and MRI environments. Different treatment plans were designed on MRI images and sonications were executed on agar-based phantoms and polymer films.

RESULTS

Magnetic Resonance (MR) thermometry maps were acquired in the agar-based phantoms. An exceptional agreement was observed between the software-planned treatment area and the lesions produced on the polymer films.

CONCLUSIONS

The developed software was successfully integrated with the MRI and robotic system controls for performing accurate treatment planning and real-time monitoring during sonications. The software provides an extremely user-friendly interface, while in the future it could be enhanced by providing dynamic modulation of the ultrasonic parameters during the treatment process.

Reduced suprasellar cistern cerebrospinal fluid motion in patients with Parkinson's disease revealed by magnetic resonance imaging with dynamic cycling of diffusion weightings

BACKGROUND

Parkinson's disease is characterized by dopamine-responsive symptoms as well as aggregation and accumulation of a-synuclein protofibrils. New diagnostic methods assess a-synuclein aggregation characteristics from cerebrospinal fluid and recent pathophysiologic mechanisms suggest that cerebrospinal fluid circulation disruptions may precipitate a-synuclein retention. Here, we test the hypothesis that cerebrospinal fluid motion at the level of the suprasellar cistern is reduced in Parkinson's disease relative to healthy participants and this reduction relates to choroid plexus perfusion.

METHODS

Diffusion weighted imaging (spatial resolution=1.8×1.8×4 mm) magnetic resonance imaging with cycling of diffusion weightings (b-values=0, 50, 100, 200, 300, 700, and 1000 s/mm2) over the approximate kinetic range of suprasellar cistern neurofluid motion was applied at 3-Tesla in Parkinson's disease (n=27; age=66±6.7 years) and healthy (n=32; age=68±8.9 years) participants. Wilcoxon rank-sum tests were applied to test the primary hypothesis that the decay rate of cerebrospinal fluid signal as a function of b-value, which reflects increasing fluid motion, is reduced in persons with versus without Parkinson's disease and inversely relates to choroid plexus activity assessed from perfusion-weighted magnetic resonance imaging (Spearman rank-order correlation; significance-criteria: p<0.05).

RESULTS

Consistent with the primary hypothesis, decay rates were higher in healthy (D=0.00328±0.00123mm2/s) relative to Parkinson's disease (D=0.00256±0.0094mm2/s) participants (p=0.016). This finding was preserved after controlling for age and sex. An inverse correlation between choroid plexus perfusion and decay rate (p=0.011) was observed in Parkinson's disease participants.

CONCLUSIONS

Cerebrospinal fluid motion at the level of the suprasellar cistern is often reduced in adults with versus without Parkinson's disease and this reduction correlates on average with choroid plexus perfusion.

Brain imaging with portable low-field MRI

The advent of portable, low-field MRI (LF-MRI) heralds new opportunities in neuroimaging. Low power requirements and transportability have enabled scanning outside the controlled environment of a conventional MRI suite, enhancing access to neuroimaging for indications that are not well suited to existing technologies. Maximizing the information extracted from the reduced signal-to-noise ratio of LF-MRI is crucial to developing clinically useful diagnostic images. Progress in electromagnetic noise cancellation and machine learning reconstruction algorithms from sparse k-space data as well as new approaches to image enhancement have now enabled these advancements. Coupling technological innovation with bedside imaging creates new prospects in visualizing the healthy brain and detecting acute and chronic pathological changes. Ongoing development of hardware, improvements in pulse sequences and image reconstruction, and validation of clinical utility will continue to accelerate this field. As further innovation occurs, portable LF-MRI will facilitate the democratization of MRI and create new applications not previously feasible with conventional systems.

Risk factors and clinical features associated with basal ganglia manganese deposition in patients with hereditary hemorrhagic telangiectasia

BACKGROUND

T1-hyperintensity of the basal ganglia (BG) due to manganese deposition is a known radiologic finding in patients with hereditary hemorrhagic telangiectasia (HHT), but risk factors and associated clinical manifestations are unclear. This study conducted a quantitative analysis of the association of T1-hyperintensity in HHT patients with specific risk factors, signs, and symptoms.

METHODS

Patients seen at our center between 2005 and 2020 with a definitive diagnosis of HHT who had an available non-contrast T1-weighted brain MRI were included. Hyperintensity was evaluated using oval regions of interest measurements. The BG: thalamus intensity ratio was used to quantitatively evaluate T1-hyperintensity. Patient laboratory values and clinical findings were collected from electronic medical records. Hyperintensity was analyzed for its association with laboratory values, and clinical findings. Variables were analyzed through regression analysis.

RESULTS

A total of 239 patients were included in this study. On 1.5 T scanners, values that were significant on multivariable regression analysis were age (p < .001), hepatic AVMs (p < .001), iron deficiency anemia (p = .0021), and cirrhosis (p = .016). On 3 T scanners, values that were significant on multivariable analysis were hepatic AVMs (p = .0024) and cirrhosis (p = .0056). On 3 T scanners, hyperintensity was significantly associated with tremor (OR = 1.17, p = .033), restless leg syndrome (OR = 1.22, p = .0086), and memory problems (OR = 1.17, p = .046).

CONCLUSIONS

BG hyperintensity due to manganese deposition is significantly associated with hepatic risk factors on 1.5 T and 3 T scanners and iron deficiency anemia on 1.5 T scanners. On 3 T scanners, T1-hyperintensity is associated with neuropsychiatric signs and symptoms, such as tremor, restless leg syndrome, and memory problems.

Non-invasive flow mapping of parasagittal meningeal lymphatics using 2D interslice flow saturation MRI

The clearance pathways of brain waste products in humans are still under debate in part due to the lack of noninvasive imaging techniques for meningeal lymphatic vessels (mLVs). In this study, we propose a new noninvasive mLVs imaging technique based on an inter-slice blood perfusion MRI called alternate ascending/descending directional navigation (ALADDIN). ALADDIN with inversion recovery (IR) at single inversion time of 2300 ms (single-TI IR-ALADDIN) clearly demonstrated parasagittal mLVs around the human superior sagittal sinus (SSS) with better detectability and specificity than the previously suggested noninvasive imaging techniques. While in many studies it has been difficult to detect mLVs and confirm their signal source noninvasively, the detection of mLVs in this study was confirmed by their posterior to anterior flow direction and their velocities and morphological features, which were consistent with those from the literature. In addition, IR-ALADDIN was compared with contrast-enhanced black blood imaging to confirm the detection of mLVs and its similarity. For the quantification of flow velocity of mLVs, IR-ALADDIN was performed at three inversion times of 2000, 2300, and 2600 ms (three-TI IR-ALADDIN) for both a flow phantom and humans. For this preliminary result, the flow velocity of the dorsal mLVs in humans ranged between 2.2 and 2.7 mm/s. Overall, (i) the single-TI IR-ALADDIN can be used as a novel non-invasive method to visualize mLVs in the whole brain with scan time of ~ 17 min and (ii) the multi-TI IR-ALADDIN can be used as a way to quantify the flow velocity of mLVs with a scan time of ~ 10 min (or shorter) in a limited coverage. Accordingly, the suggested approach can be applied to noninvasively studying meningeal lymphatic flows in general and also understanding the clearance pathways of waste production through mLVs in humans, which warrants further investigation.

Red cell exchange transfusions increase cerebral capillary transit times and may alter oxygen extraction in sickle cell disease

Persons with sickle cell disease (SCD) suffer from chronic hemolytic anemia, reduced blood oxygen content, and lifelong risk of silent and overt stroke. Major conventional stroke risk factors are absent in most individuals with SCD, yet nearly 50% have evidence of brain infarcts by the age of 30 years, indicating alternative etiologies for ischemia. We investigated whether radiological evidence of accelerated blood water transit through capillaries, visible on arterial spin labeling (ASL) magnetic resonance imaging, reduces following transfusion-induced increases in hemoglobin and relates to oxygen extraction fraction (OEF). Neurological evaluation along with anatomical and hemodynamic imaging with cerebral blood flow (CBF)-weighted pseudocontinuous ASL and OEF imaging with T2 -relaxation-under-spin-tagging were applied in sequence before and after blood transfusion therapy (n = 32) and in a comparator cohort of nontransfused SCD participants on hydroxyurea therapy scanned at two time points to assess stability without interim intervention (n = 13). OEF was calculated separately using models derived from human hemoglobin-F, hemoglobin-A, and hemoglobin-S. Gray matter CBF and dural sinus signal, indicative of rapid blood transit, were evaluated at each time point and compared with OEF using paired statistical tests (significance: two-sided p < 0.05). No significant change in sinus signal was observed in nontransfused participants (p = 0.650), but a reduction was observed in transfused participants (p = 0.034), consistent with slower red cell transit following transfusion. The dural sinus signal intensity was inversely associated with OEF pretransfusion (p = 0.011), but not posttransfusion. Study findings suggest that transfusion-induced increases in total hemoglobin may lengthen blood transit times through cerebral capillaries and alter cerebral OEF in SCD.

Silent infarction in sickle cell disease is associated with brain volume loss in excess of infarct volume

Introduction

Sickle cell disease (SCD) increases cerebral infarct risk, but reported effects on brain volume have varied. More detailed information using larger cohorts and contemporary methods could motivate the use of longitudinal brain volume assessment in SCD as an automated marker of disease stability or future progression. The purpose of this study was to rigorously evaluate whether children and young adults with SCD have reduced gray matter volume (GMV) and white matter volume (WMV) compared to healthy controls using high-resolution MRI. We tested the hypotheses that (i) elevated CBF, a marker of cerebral hemodynamic compensation in SCD, is associated with global and regional brain atrophy, and (ii) silent cerebral infarct burden is associated with brain atrophy in excess of infarct volume.

Methods

Healthy controls (n = 49) and SCD participants without overt stroke (n = 88) aged 7-32 years completed 3 T brain MRI; pseudocontinuous arterial spin labeling measured CBF. Multivariable linear regressions assessed associations of independent variables with GMV, WMV, and volumes of cortical/subcortical regions.

Results

Reduced hemoglobin was associated with reductions in both GMV (p = 0.032) and WMV (p = 0.005); reduced arterial oxygen content (CaO2) was also associated with reductions in GMV (p = 0.035) and WMV (p = 0.006). Elevated gray matter CBF was associated with reduced WMV (p = 0.018). Infarct burden was associated with reductions in WMV 30-fold greater than the infarct volume itself (p = 0.005). Increased GM CBF correlated with volumetric reductions of the insula and left and right caudate nuclei (p = 0.017, 0.017, 0.036, respectively). Infarct burden was associated with reduced left and right nucleus accumbens, right thalamus, and anterior corpus callosum volumes (p = 0.002, 0.002, 0.009, 0.002, respectively).

Discussion

We demonstrate that anemia and decreased CaO2 are associated with reductions in GMV and WMV in SCD. Increased CBF and infarct burden were also associated with reduced volume in subcortical structures. Global WMV deficits associated with infarct burden far exceed infarct volume itself. Hemodynamic compensation via increased cerebral blood flow in SCD seems inadequate to prevent brain volume loss. Our work highlights that silent cerebral infarcts are just a portion of the brain injury that occurs in SCD; brain volume is another potential biomarker of brain injury in SCD.

Neurochemical and functional reorganization of the cognitive-ear link underlies cognitive impairment in presbycusis

Recent studies suggest that the interaction between presbycusis and cognitive impairment may be partially explained by the cognitive-ear link. However, the underlying neurophysiological mechanisms remain largely unknown. In this study, we combined magnetic resonance spectroscopy (MRS) and resting-state functional magnetic resonance imaging (fMRI) to investigate auditory gamma-aminobutyric acid (GABA) and glutamate (Glu) levels, intra- and inter-network functional connectivity, and their relationships with auditory and cognitive function in 51 presbycusis patients and 51 well-matched healthy controls. Our results confirmed reorganization of the cognitive-ear link in presbycusis, including decreased auditory GABA and Glu levels and aberrant functional connectivity involving auditory networks (AN) and cognitive-related networks, which were associated with reduced speech perception or cognitive impairment. Moreover, mediation analyses revealed that decreased auditory GABA levels and dysconnectivity between the AN and default mode network (DMN) mediated the association between hearing loss and impaired information processing speed in presbycusis. These findings highlight the importance of AN-DMN dysconnectivity in cognitive-ear link reorganization leading to cognitive impairment, and hearing loss may drive reorganization via decreased auditory GABA levels. Modulation of GABA neurotransmission may lead to new treatment strategies for cognitive impairment in presbycusis patients.

Silicone-based materials with tailored MR relaxation characteristics for use in reduced coil visibility and in tissue-mimicking phantom design

BACKGROUND

The development of materials with tailored signal intensity in MR imaging is critically important both for the reduction of signal from non-tissue hardware, as well as for the construction of tissue-mimicking phantoms. Silicone-based phantoms are becoming more popular due to their structural stability, stretchability, longer shelf life, and ease of handling, as well as for their application in dynamic imaging of physiology in motion. Moreover, silicone can be also used for the design of stretchable receive radio-frequency (RF) coils.

PURPOSE

Fabrication of materials with tailored signal intensity for MRI requires knowledge of precise T₁ and T₂ relaxation times of the materials used. In order to increase the range of possible relaxation times, silicone materials can be doped with gadolinium (Gd). In this work, we aim to systematically evaluate relaxation properties of Gd-doped silicone material at a broad range of Gd concentrations and at three clinically relevant magnetic field strengths (1.5 T, 3 T, and 7 T). We apply the findings for rendering silicone substrates of stretchable receive RF coils less visible in MRI. Moreover, we demonstrate early stage proof-of-concept applicability in tissue-mimicking phantom development.

MATERIALS AND METHODS

Ten samples of pure and Gd-doped Ecoflex silicone polymer samples were prepared with various Gd volume ratios ranging from 1:5000 to 1:10, and studied using 1.5 T and 3 T clinical and 7 T preclinical scanners. T₁ and T₂ relaxation times of each sample were derived by fitting the data to Bloch signal intensity equations. A receive coil made from Gd-doped Ecoflex silicone polymer was fabricated and evaluated in vitro at 3 T.

RESULTS

With the addition of a Gd-based contrast agent, it is possible to significantly change T₂ relaxation times of Ecoflex silicone polymer (from 213 ms to 20 ms at 1.5 T; from 135 ms to 17 ms at 3 T; and from 111.4 ms to 17.2 ms at 7 T). T₁ relaxation time is less affected by the introduction of the contrast agent (changes from 608 ms to 579 ms; from 802.5 ms to 713 ms at 3 T; from 1276 ms to 979 ms at 7 T). First results also indicate that liver, pancreas, and white matter tissues can potentially be closely mimicked using this phantom preparation technique. Gd-doping reduces the appearance of the silicone-based coil substrate during the MR scan by up to 81%.

CONCLUSIONS

Gd-based contrast agents can be effectively used to create Ecoflex silicone polymer-based phantoms with tailored T₂ relaxation properties. The relative low cost, ease of preparation, stretchability, mechanical stability, and long shelf life of Ecoflex silicone polymer all make it a good candidate for "MR invisible" coil development and bears promise for tissue-mimicking phantom development applicability.

Replicability of proton MR spectroscopic imaging findings in mild traumatic brain injury: Implications for clinical applications

PURPOSE

Proton magnetic resonance spectroscopy (1H MRS) offers biomarkers of metabolic damage after mild traumatic brain injury (mTBI), but a lack of replicability studies hampers clinical translation. In a conceptual replication study design, the results reported in four previous publications were used as the hypotheses (H1-H7), specifically: abnormalities in patients are diffuse (H1), confined to white matter (WM) (H2), comprise low N-acetyl-aspartate (NAA) levels and normal choline (Cho), creatine (Cr) and myo-inositol (mI) (H3), and correlate with clinical outcome (H4); additionally, a lack of findings in regional subcortical WM (H5) and deep gray matter (GM) structures (H6), except for higher mI in patients' putamen (H7).

METHODS

26 mTBI patients (20 female, age 36.5 ± 12.5 [mean ± standard deviation] years), within two months from injury and 21 age-, sex-, and education-matched healthy controls were scanned at 3 Tesla with 3D echo-planar spectroscopic imaging. To test H1-H3, global analysis using linear regression was used to obtain metabolite levels of GM and WM in each brain lobe. For H4, patients were stratified into non-recovered and recovered subgroups using the Glasgow Outcome Scale Extended. To test H5-H7, regional analysis using spectral averaging estimated metabolite levels in four GM and six WM structures segmented from T1-weighted MRI. The Mann-Whitney U test and weighted least squares analysis of covariance were used to examine mean group differences in metabolite levels between all patients and all controls (H1-H3, H5-H7), and between recovered and non-recovered patients and their respectively matched controls (H4). Replicability was defined as the support or failure to support the null hypotheses in accordance with the content of H1-H7, and was further evaluated using percent differences, coefficients of variation, and effect size (Cohen's d).

RESULTS

Patients' occipital lobe WM Cho and Cr levels were 6.0% and 4.6% higher than controls', respectively (Cho, d = 0.37, p = 0.04; Cr, d = 0.63, p = 0.03). The same findings, i.e., higher patients' occipital lobe WM Cho and Cr (both p = 0.01), but with larger percent differences (Cho, 8.6%; Cr, 6.3%) and effect sizes (Cho, d = 0.52; Cr, d = 0.88) were found in the comparison of non-recovered patients to their matched controls. For the lobar WM Cho and Cr comparisons without statistical significance (frontal, parietal, temporal), unidirectional effect sizes were observed (Cho, d = 0.07 - 0.37; Cr, d = 0.27 - 0.63). No differences were found in any metabolite in any lobe in the comparison between recovered patients and their matched controls. In the regional analyses, no differences in metabolite levels were found in any GM or WM region, but all WM regions (posterior, frontal, corona radiata, and the genu, body, and splenium of the corpus callosum) exhibited unidirectional effect sizes for Cho and Cr (Cho, d = 0.03 - 0.34; Cr, d = 0.16 - 0.51).

CONCLUSIONS

We replicated findings of diffuse WM injury, which correlated with clinical outcome (supporting H1-H2, H4). These findings, however, were among the glial markers Cho and Cr, not the neuronal marker NAA (not supporting H3). No differences were found in regional GM and WM metabolite levels (supporting H5-H6), nor in putaminal mI (not supporting H7). Unidirectional effect sizes of higher patients' Cho and Cr within all WM analyses suggest widespread injury, and are in line with the conclusion from the previous publications, i.e., that detection of WM injury may be more dependent upon sensitivity of the 1H MRS technique than on the selection of specific regions. The findings lend further support to the corollary that clinic-ready 1H MRS biomarkers for mTBI may best be achieved by using high signal-to-noise-ratio single-voxels placed anywhere within WM. The biochemical signature of the injury, however, may differ and therefore absolute levels, rather than ratios may be preferred. Future replication efforts should further test the generalizability of these findings.

White matter microstructure and longitudinal relaxation time anisotropy in human brain at 3 and 7 T

A high degree of structural order by white matter (WM) fibre tracts creates a physicochemical environment where water relaxations are rendered anisotropic. Recently, angularly dependent longitudinal relaxation has been reported in human WM. We have characterised interrelationships between T1 relaxation and diffusion MRI microstructural indices at 3 and 7 T. Eleven volunteers consented to participate in the study. Multishell diffusion MR images were acquired with b-values of 0/1500/3000 and 0/1000/2000 s/mm² at 1.5 and 1.05 mm³ isotropic resolutions at 3 and 7 T, respectively. DTIFIT was used to compute DTI indices; the fibre-to-field angle (θ_FB ) maps were obtained using the principal eigenvector images. The orientations and volume fractions of multiple fibre populations were estimated using BedpostX in FSL, and the orientation dispersion index (ODI) was estimated using the NODDI protocol. MP2RAGE was used to acquire images for T1 maps at 1.0 and 0.9 mm³ isotropic resolutions at 3 and 7 T, respectively. At 3 T, T1 as a function of θ_FB in WM with high fractional anisotropy and one-fibre orientation volume fraction or low ODI shows a broad peak centred at 50^o , but a flat baseline at 0^o and 90^o . The broad peak amounted up to 7% of the mean T1. At 7 T, the broad peak appeared at 40^o and T1 in fibres running parallel to B0 was longer by up to 75 ms (8.3% of the mean T1) than in those perpendicular to the field. The peak at 40^o was approximately 5% of mean T1 (i.e., proportionally smaller than that at 54^o at 3 T). The data demonstrate T1 anisotropy in WM with high microstructural order at both fields. The angular patterns are indicative of the B0-dependency of T1 anisotropy. Thus myelinated WM fibres influence T1 contrast both by acting as a T1 contrast agent and rendering T1 dependent on fibre orientation with B0.

CEST MR fingerprinting (CEST-MRF) for brain tumor quantification using EPI readout and deep learning reconstruction

PURPOSE

To develop a clinical CEST MR fingerprinting (CEST-MRF) method for brain tumor quantification using EPI acquisition and deep learning reconstruction.

METHODS

A CEST-MRF pulse sequence originally designed for animal imaging was modified to conform to hardware limits on clinical scanners while keeping scan time under 2 min. Quantitative MRF reconstruction was performed using a deep reconstruction network (DRONE) to yield the water relaxation and chemical exchange parameters. The feasibility of the six parameter DRONE reconstruction was tested in simulations using a digital brain phantom. A healthy subject was scanned with the CEST-MRF sequence, conventional MRF and CEST sequences for comparison. Reproducibility was assessed via test-retest experiments and the concordance correlation coefficient calculated for white matter and gray matter. The clinical utility of CEST-MRF was demonstrated on four patients with brain metastases in comparison to standard clinical imaging sequences. Tumors were segmented into edema, solid core, and necrotic core regions and the CEST-MRF values compared to the contra-lateral side.

RESULTS

DRONE reconstruction of the digital phantom yielded a normalized RMS error of ≤7% for all parameters. The CEST-MRF parameters were in good agreement with those from conventional MRF and CEST sequences and previous studies. The mean concordance correlation coefficient for all six parameters was 0.98 ± 0.01 in white matter and 0.98 ± 0.02 in gray matter. The CEST-MRF values in nearly all tumor regions were significantly different (P = 0.05) from each other and the contra-lateral side.

CONCLUSION

Combination of EPI readout and deep learning reconstruction enabled fast, accurate and reproducible CEST-MRF in brain tumors.

A New Combination of Radio-Frequency Coil Configurations Using High-Permittivity Materials and Inductively Coupled Structures for Ultrahigh-Field Magnetic Resonance Imaging

In ultrahigh-field (UHF) magnetic resonance imaging (MRI) system, the RF power required to excite the nuclei of the target object increases. As the strength of the main magnetic field (B0 field) increases, the improvement of the RF transmit field (B1+ field) efficiency and receive field (B1- field) sensitivity of radio-frequency (RF) coils is essential to reduce their specific absorption rate and power deposition in UHF MRI. To address these problems, we previously proposed a method to simultaneously improve the B1+ field efficiency and B1- field sensitivity of 16-leg bandpass birdcage RF coils (BP-BC RF coils) by combining a multichannel wireless RF element (MCWE) and segmented cylindrical high-permittivity material (scHPM) comprising 16 elements in 7.0 T MRI. In this work, we further improved the performance of transmit/receive RF coils. A new combination of RF coil with wireless element and HPM was proposed by comparing the BP-BC RF coil with the MCWE and the scHPM proposed in the previous study and the multichannel RF coils with a birdcage RF coil-type wireless element (BCWE) and the scHPM proposed in this study. The proposed 16-ch RF coils with the BCWE and scHPM provided excellent B1+ field efficiency and B1- field sensitivity improvement.

Effect of Normal Breathing on the Movement of CSF in the Spinal Subarachnoid Space

BACKGROUND AND PURPOSE

Forced respirations reportedly have an effect on CSF movement in the spinal canal. We studied respiratory-related CSF motion during normal respiration.

MATERIALS AND METHODS

Six healthy subjects breathed at their normal rate with a visual guide to ensure an unchanging rhythm. Respiratory-gated phase-contrast MR flow images were acquired at 5 selected axial planes along the spine. At each spinal level, we computed the flow rate voxelwise in the spinal canal, together with the associated stroke volume. From these data, we computed the periodic volume changes of spinal segments. A phantom was used to quantify the effect of respiration-related magnetic susceptibility changes on the velocity data measured.

RESULTS

At each level, CSF moved cephalad during inhalation and caudad during expiration. While the general pattern of fluid movement was the same in the 6 subjects, the flow rates, stroke volumes, and spine segment volume changes varied among subjects. Peak flow rates ranged from 0.60 to 1.59 mL/s in the cervical region, 0.46 to 3.17 mL/s in the thoracic region, and 0.75 to 3.64 mL/s in the lumbar region. The differences in flow rates along the canal yielded cyclic volume variations of spine segments that were largest in the lumbar spine, ranging from 0.76 to 3.07 mL among subjects. In the phantom study, flow velocities oscillated periodically during the respiratory cycle by up to 0.02 cm/s or 0.5%.

CONCLUSIONS

Respiratory-gated measurements of the CSF motion in the spinal canal showed cyclic oscillatory movements of spinal fluid correlated to the breathing pattern.

Synthetic MRI improves radiomics-based glioblastoma survival prediction

Glioblastoma is an aggressive and fast-growing brain tumor with poor prognosis. Predicting the expected survival of patients with glioblastoma is a key task for efficient treatment and surgery planning. Survival predictions could be enhanced by means of a radiomic system. However, these systems demand high numbers of multicontrast images, the acquisitions of which are time consuming, giving rise to patient discomfort and low healthcare system efficiency. Synthetic MRI could favor deployment of radiomic systems in the clinic by allowing practitioners not only to reduce acquisition time, but also to retrospectively complete databases or to replace artifacted images. In this work we analyze the replacement of an actually acquired MR weighted image by a synthesized version to predict survival of glioblastoma patients with a radiomic system. Each synthesized version was realistically generated from two acquired images with a deep learning synthetic MRI approach based on a convolutional neural network. Specifically, two weighted images were considered for the replacement one at a time, a T2w and a FLAIR, which were synthesized from the pairs T1w and FLAIR, and T1w and T2w, respectively. Furthermore, a radiomic system for survival prediction, which can classify patients into two groups (survival >480 days and ≤ 480 days), was built. Results show that the radiomic system fed with the synthesized image achieves similar performance compared with using the acquired one, and better performance than a model that does not include this image. Hence, our results confirm that synthetic MRI does add to glioblastoma survival prediction within a radiomics-based approach.

A relaxometry method that emphasizes practicality and availability

PURPOSE

Although many methods have been proposed to quantitatively map the main MRI parameters (e.g., T₁ , T₂ , C × M₀ ), these methods often involve special sequences not readily available on clinical scanners and/or may require long scan times. In contrast, the proposed method can readily run on most scanners, offer flexible tradeoffs between scan time and image quality, and map MRI parameters jointly to ensure spatial alignment.

METHODS

The approach is based on the multi-shot spin-echo (SE) EPI sequence. The corresponding signal equation was derived and strategies for solving it were developed. As usual with multi-shot EPI, scan time can readily be traded-off against image quality by adjusting the echo train length. Validation was performed against reference relaxometry methods, in gel phantoms with varying concentrations of gadobutrol and gadoterate meglumine contrast agents. In vivo examples are further presented, from 3 neuroradiology patients.

RESULTS

Bland-Altman analysis was performed: for T₂ , as compared to 2D SE, bias was 0.29 ms and the 95% limits of agreement ranged from -1.15 to +1.73 ms. For T₁ , compared to inversion-recovery SE (and MOLLI), bias was -20.2 ms (and -14.5 ms) and the limits of agreement ranged from -62.4 to +22.0 ms (and -53.8 to +24.9 ms). The mean relative T₁ error between the proposed method and each of the 2 reference methods was similar to that of the reference methods among themselves.

CONCLUSION

In the constellation of existing relaxometry methods, the proposed method is meant to stand out in terms of its practicality and availability.

Interleaved fluid-attenuated inversion recovery (FLAIR) MRI and deuterium metabolic imaging (DMI) on human brain in vivo

PURPOSE

To integrate deuterium metabolic imaging (DMI) with clinical MRI through an interleaved MRI and DMI acquisition workflow. Interleaved MRI-DMI was enabled with hardware and pulse sequence modifications, and the performance was demonstrated using fluid-attenuated inversion recovery (FLAIR) MRI as an example.

METHODS

Interleaved FLAIR-DMI was developed by interleaving the ² H excitation and acquisition time windows into the intrinsic delay periods presented in the FLAIR method. All ² H MR signals were up-converted to the ¹ H Larmor frequency using a custom-built hardware unit, which also achieved frequency and phase locking of the output signal in real-time. The interleaved measurements were compared with direct measurements both in phantom and in the human brain in vivo.

RESULTS

The interleaved MRI-DMI acquisition strategy allowed simultaneous detection of FLAIR MRI and DMI in the same scan time as a FLAIR-only MRI acquisition. Both phantom and in vivo data showed that the MR image quality, DMI sensitivity as well as information content were preserved using interleaved MRI-DMI.

CONCLUSION

The interleaved MRI-DMI technology can be used to extend clinical MRI protocols with DMI, thereby offering a metabolic component to the MR imaging contrasts without a penalty on patient comfort or scan time.

In vivo measurements of lamination patterns in the human cortex

The laminar composition of the cerebral cortex is tightly connected to the development and connectivity of the brain, as well as to function and pathology. Although most of the research on the cortical layers is done with the aid of ex vivo histology, there have been recent attempts to use magnetic resonance imaging (MRI) with potential in vivo applications. However, the high-resolution MRI technology and protocols required for such studies are neither common nor practical. In this article, we present a clinically feasible method for assessing the laminar properties of the human cortex using standard pulse sequence available on any common MRI scanner. Using a series of low-resolution inversion recovery (IR) MRI scans allows us to calculate multiple T1 relaxation time constants for each voxel. Based on the whole-brain T1 -distribution, we identify six different gray matter T1 populations and their variation across the cortex. Based on this, we show age-related differences in these population and demonstrate that this method is able to capture the difference in laminar composition across varying brain areas. We also provide comparison to ex vivo high-resolution MRI scans. We show that this method is feasible for the estimation of layer variability across large population cohorts, which can lead to research into the links between the cortical layers and function, behavior and pathologies that was heretofore unexplorable.

[Material Investigation for the Development of Non-rigid Phantoms for CT-MRI Image Registration]

PURPOSE

In radiotherapy, deformable image registration (DIR) has been frequently used in different imaging examinations in recent years. However, no phantom has been established for quality assurance for DIR. In order to develop a non-rigid phantom for accuracy control between CT and MRI images, we investigated the suitability of 3D printing materials and gel materials in this study.

METHODS

We measured CT values, T₁ values, T₂ values, and the proton densities of 31 3D printer materials-purchased from three manufacturers-and one gel material. The dice coefficient after DIR was calculated for the CT-MRI images using a prototype phantom made of a gel material compatible with CT-MRI.

RESULTS

The CT number of the 3D printing materials ranged from -6.8 to 146.4 HU. On MRI, T₁ values were not measurable in most cases, whereas T₂ values were not measurable in all cases; proton density (PD) ranged from 2.51% to 4.9%. The gel material had a CT number of 111.16 HU, T₁ value of 813.65 ms, and T₂ value of 27.19 ms. The prototype phantom was flexible, and the usefulness of DIR with CT and MRI images was demonstrated using this phantom.

CONCLUSION

The CT number and T₁ and T₂ values of the gel material are close to those of the human body and may therefore be developed as a DIR verification phantom between CT and MRI. These findings may contribute to the development of non-rigid phantoms for DIR in the future.

Measurement of CSF pulsation from EPI-based human fMRI

It is recently discovered that the glymphatic system and meningeal lymphatic system are the primary routes for the clearance of brain waste products. The CSF flow is part of these systems, facilitating the clearance procedure. Nonetheless, the relationship between CSF flow and brain functional activity has been underexplored. To investigate CSF dynamics and functional brain activity simultaneously, recent studies have proposed a CSF inflow index measured on edge slices (CSFedge) of echo-planar imaging (EPI) based functional magnetic resonance imaging (fMRI), however, it lacks the quantitative aspect of the CSF pulsation. We proposed a new method for quantifying CSF pulsation (CSFpulse) based on an interslice CSF pulsation model in the 4^th ventricle of EPI-based fMRI. The proposed CSFpulse successfully detected the higher CSF flow during the resting state than the typical task states (visual and motor) (p<.05), which is consistent with previous studies based on phase contrast (PC) MRI and CSF volume MRI, while it was not detected in CSFedge based indices or baseline CSF signals in various regions of interest (ROIs). Moreover, CSFpulse demonstrated dynamic functional changes in CSF pulsation: it decreased during the activation-on blocks while it increased during the activation-off blocks. CSFpulse significantly correlated with stroke volume measured using PC MRI, a standard method for CSF pulsation quantification, under the same functional state, while CSFedge based indices or CSF ROIs showed no correlation with the PC MRI stroke volume. Lastly, the correlation of CSFpulse with global BOLD was weaker than that of CSFedge, suggesting that CSFpulse may reflect distinct CSF physiological information that is less affected by global BOLD effects. Based on these results, the proposed CSFpulse provides CSF pulsatility information more accurately in a quantitative manner than CSFedge based indices from the recent CSF studies or the conventional ROI-based analysis. In addition to the high correlation with PC MRI, CSFpulse is much faster than PC MRI and provides information of functional brain activations simultaneously, advantageous over PC MRI or CSF volume MRI. Accordingly, the suggested CSFpulse can be used for investigating intra-subject functional changes in BOLD and CSF pulsation simultaneously and inter-subject CSF pulsation variations based on conventional EPI-based fMRI, which warrants further investigation.

Impact of tissue properties on time-dependent alterations in apparent diffusion coefficient: a phantom study using oscillating-gradient spin-echo and pulsed-gradient spin-echo sequences

Purpose

The purpose of this study was to investigate whether the changes in apparent diffusion coefficients (ADCs) due to differences in diffusion time reflect tissue properties in actual measurements of phantoms.

Materials and methods

Various n-alkane phantoms and sucrose/collagen phantoms with various collagen densities were set up with and without polyvinyl alcohol (PVA) foam with an average pore diameter of 300 μm. Thus, n-alkanes or sucrose/collagen represented substrate viscosity and the presence of PVA foam represented tissue structure with septum. Diffusion-weighted images with various diffusion times (7.71–60 ms) were acquired using pulsed-gradient spin-echo (PGSE) and oscillating-gradient spin-echo (OGSE) sequences. The ADCs of the phantoms with and without PVA foam were calculated.

Results

The ADCs of some of the phantoms without PVA decreased with diffusion times decreased. In the n-alkane phantoms, only C₈H₁₈ showed significantly different ADCs depending on the use of PVA foam in the OGSE sequence. On the other hand, sucrose/collagen phantoms showed significant differences according to diffusion time. The ADCs of the phantoms decreased as the molecular size of the n-alkanes or collagen density of the sucrose/collagen phantom increased. Compared to phantoms without PVA foam, the ADC of the phantoms with PVA foam decreased as the diffusion time increased.

Conclusion

Changes in ADCs due to differences in diffusion time reflect tissue properties in actual measurements of phantoms. These changes in ADCs can be used for tissue characterization in vivo.

Three-dimensional simultaneous brain mapping of T1, T2, T2∗ and magnetic susceptibility with MR Multitasking

PURPOSE

To develop a new technique that enables simultaneous quantification of whole-brain T1 , T2 , T 2 ∗ , as well as susceptibility and synthesis of six contrast-weighted images in a single 9.1-minute scan.

METHODS

The technique uses hybrid T2 -prepared inversion-recovery pulse modules and multi-echo gradient-echo readouts to collect k-space data with various T1, T2, and T 2 ∗ weightings. The underlying image is represented as a six-dimensional low-rank tensor consisting of three spatial dimensions and three temporal dimensions corresponding to T1 recovery, T2 decay, and multi-echo behaviors, respectively. Multiparametric maps were fitted from reconstructed image series. The proposed method was validated on phantoms and healthy volunteers, by comparing quantitative measurements against corresponding reference methods. The feasibility of generating six contrast-weighted images was also examined.

RESULTS

High quality, co-registered T1 , T2 , and T 2 ∗ susceptibility maps were generated that closely resembled the reference maps. Phantom measurements showed substantial consistency (R2 > 0.98) with the reference measurements. Despite the significant differences of T1 (p < .001), T2 (p = .002), and T 2 ∗ (p = 0.008) between our method and the references for in vivo studies, excellent agreement was achieved with all intraclass correlation coefficients greater than 0.75. No significant difference was found for susceptibility (p = .900). The framework is also capable of synthesizing six contrast-weighted images.

CONCLUSION

The MR Multitasking-based 3D brain mapping of T1 , T2 , T 2 ∗ , and susceptibility agrees well with the reference and is a promising technique for multicontrast and quantitative imaging.

Avoiding data loss: Synthetic MRIs generated from diffusion imaging can replace corrupted structural acquisitions for freesurfer-seeded tractography

Magnetic Resonance Imaging (MRI) motion artefacts frequently complicate structural and diffusion MRI analyses. While diffusion imaging is easily ‘scrubbed’ of motion affected volumes, the same is not true for T1w or T2w ‘structural’ images. Structural images are critical to most diffusion-imaging pipelines thus their corruption can lead to disproportionate data loss. To enable diffusion-image processing when structural images are missing or have been corrupted, we propose a means by which synthetic structural images can be generated from diffusion MRI. This technique combines multi-tissue constrained spherical deconvolution, which is central to many existing diffusion analyses, with the Bloch equations that allow simulation of MRI intensities for given scanner parameters and magnetic resonance (MR) tissue properties. We applied this technique to 32 scans, including those acquired on different scanners, with different protocols and with pathology present. The resulting synthetic T1w and T2w images were visually convincing and exhibited similar tissue contrast to acquired structural images. These were also of sufficient quality to drive a Freesurfer-based tractographic analysis. In this analysis, probabilistic tractography connecting the thalamus to the primary sensorimotor cortex was delineated with Freesurfer, using either real or synthetic structural images. Tractography for real and synthetic conditions was largely identical in terms of both voxels encountered (Dice 0.88–0.95) and mean fractional anisotropy (intrasubject absolute difference 0.00–0.02). We provide executables for the proposed technique in the hope that these may aid the community in analysing datasets where structural image corruption is common, such as studies of children or cognitively impaired persons.

Comparative Analysis of Blood T2 Values Measured by T2 -TRIR and TRUST

BACKGROUND

Venous blood oxygenation (Yv), which can be derived from venous blood T₂ (T₂ b), combined with oxygen-extraction fraction (OEF) and cerebral metabolic rate of oxygen, is considered indicative for tissue viability and brain functioning and frequently assessed in patients with sickle cell disease. Recently, T₂ -Prepared-Blood-Relaxation-Imaging-with-Inversion-Recovery (T₂ -TRIR) was introduced allowing for simultaneous measurements of blood T₂ and T₁ (T₁ b), potentially improving Yv estimation by overcoming the need to estimate hematocrit.

PURPOSE

To optimize and compare T₂ -TRIR with T₂ -relaxation-under-spin-tagging (TRUST) sequence.

STUDY TYPE

Prospective.

POPULATION

A total of 12 healthy volunteers (six female, 27 ± 3 years old) and 7 patients with sickle cell disease (five female, 32 ± 12 years old).

FIELD STRENGTH/SEQUENCE

3 T; turbo field echo planar imaging (TFEPI), echo planar imaging (EPI), and fast field echo (FFE).

ASSESSMENT

T₂ b, Yv, and OEF from TRUST and T₂ -TRIR were compared and T₂ -TRIR-derived T₁ b was assessed. Within- and between-session repeatability was quantified in the controls, whereas sensitivity to hemodynamic changes after acetazolamide (ACZ) administration was assessed in the patients.

STATISTICAL TESTS

Shapiro-Wilk, one-sample and paired-sample t-test, repeated measures ANOVA, mixed linear model, Bland-Altman analysis and correlation analysis. Sidak multiple-comparison correction was performed. Significance level was 0.05.

RESULTS

In controls, T₂ b from T₂ -TRIR (70 ± 11 msec) was higher compared to TRUST (60 ± 8 msec). In patients, T₂ b values were lower pre- compared to post-ACZ administration (TRUST: 80 ± 15 msec and 106 ± 23 msec and T₂ -TRIR: 95 ± 21 msec and 125 ± 36 msec). Consequently, Yv and OEF were lower and higher pre- compared to post-ACZ administration (TRUST Yv: 68% ± 7% and 77% ± 8%, T₂ -TRIR Yv: 74% ± 8% and 80% ± 6%, TRUST OEF: 30% ± 7% and 21% ± 8%, and T₂ -TRIR OEF: 25% ± 8% and 18% ± 6%).

DATA CONCLUSION

TRUST and T₂ -TRIR are reproducible, but T₂ -TRIR-derived T₂ b values are significantly higher compared to TRUST, resulting in higher Yv and lower OEF estimates. This bias might be considered when evaluating cerebral oxygen homeostasis.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Measuring Arterial Pulsatility With Dynamic Inflow Magnitude Contrast

The pulsatility of blood flow through cerebral arteries is clinically important, as it is intrinsically associated with cerebrovascular health. In this study we outline a new MRI approach to measuring the real-time pulsatile flow in cerebral arteries, which is based on the inflow phenomenon associated with fast gradient-recalled-echo acquisitions. Unlike traditional phase-contrast techniques, this new method, which we dub dynamic inflow magnitude contrast (DIMAC), does not require velocity-encoding gradients as sensitivity to flow velocity is derived purely from the inflow effect. We achieved this using a highly accelerated single slice EPI acquisition with a very short TR (15 ms) and a 90° flip angle, thus maximizing inflow contrast. We simulate the spoiled GRE signal in the presence of large arteries and perform a sensitivity analysis. The sensitivity analysis demonstrates that in the regime of high inflow contrast, DIMAC shows much greater sensitivity to flow velocity over blood volume changes. We support this theoretical prediction with in-vivo data collected in two separate experiments designed to demonstrate the utility of the DIMAC signal contrast. We perform a hypercapnia challenge experiment in order to experimentally modulate arterial tone within subjects, and thus modulate the arterial pulsatile flow waveform. We also perform a thigh-cuff release challenge, designed to induce a transient drop in blood pressure, and demonstrate that the continuous DIMAC signal captures the complex transient change in the pulsatile and non-pulsatile components of flow. In summary, this study proposes a new role for a well-established source of MR image contrast and demonstrates its potential for measuring both steady-state and dynamic changes in arterial tone.

Compressed SENSE in Pediatric Brain Tumor MR Imaging : Assessment of Image Quality, Examination Time and Energy Release

Purpose

To compare the image quality, examination time, and total energy release of a standardized pediatric brain tumor magnetic resonance imaging (MRI) protocol performed with and without compressed sensitivity encoding (C-SENSE). Recently introduced as an acceleration technique in MRI, we hypothesized that C‑SENSE would improve image quality, reduce the examination time and radiofrequency-induced energy release compared with conventional examination in a pediatric brain tumor protocol.

Methods

This retrospective study included 22 patients aged 2.33–18.83 years with different brain tumor types who had previously undergone conventional MRI examination and underwent follow-up C‑SENSE examination. Both examinations were conducted with a 3.0-Tesla device and included pre-contrast and post-contrast T1-weighted turbo-field-echo, T2-weighted turbo-spin-echo, and fluid-attenuated inversion recovery sequences. Image quality was assessed in four anatomical regions of interest (tumor area, cerebral cortex, basal ganglia, and posterior fossa) using a 5-point scale. Reader preference between the standard and C‑SENSE images was evaluated. The total examination duration and energy deposit were compared based on scanner log file analysis.

Results

Relative to standard examinations, C‑SENSE examinations were characterized by shorter total examination times (26.1 ± 3.93 vs. 22.18 ± 2.31 min; P = 0.001), reduced total energy deposit (206.0 ± 19.7 vs. 92.3 ± 18.2 J/kg; P < 0.001), and higher image quality (overall P < 0.001).

Conclusion

C‑SENSE contributes to the improvement of image quality, reduction of scan times and radiofrequency-induced energy release relative to the standard protocol in pediatric brain tumor MRI.

Supplementary Information

The online version of this article (10.1007/s00062-021-01112-3) contains supplementary material, which is available to authorized users.

Relaxation time of brain tissue in the elderly assessed by synthetic MRI

BACKGROUND

Synthetic MRI (SyMRI) is a quantitative technique that allows measurements of T1 and T2 relaxation times (RTs). Brain RT evolution across lifespan is well described for the younger population. The aim was to study RTs of brain parenchyma in a healthy geriatric population in order to define the normal value of structures in this group population. Normal values for geriatric population could help find biomarker for age-related brain disease.

MATERIALS AND METHODS

Fifty-four normal-functioning individuals (22 females, 32 males) with mean age of 83 years (range 56-98) underwent SyMRI. RT values in manually defined ROIs (centrum semiovale, middle cerebellar peduncles, thalamus, and insular cortex) and in segmented whole-brain components (brain parenchyma, gray matter, white matter, myelin, CSF, and stromal structures) were extracted from the SyMRI segmentation software. Patients' results were combined into the group age. Main ROI-based and whole-brain results were compared for the all dataset and for age group results as well.

RESULTS

For white matter, RTs between ROI-based analyses and whole-brain results for T2 and for T1 were statistically different and a trend of increasing T1 in centrum semiovale and cerebellar peduncle was observed. For gray matter, thalamic T1 was statistically different from insular T1. A difference was also found between left and right insula (p < .0001). T1 RTs of ROI-based and whole-brain-based analyses were statistically different (p < .0001). No significant difference in T1 and T2 was found between age groups on ROI-based analysis, but T1 in centrum semiovale and thalamus increased with age. No statistical difference between age groups was found for the various segmented volumes except for myelin between 65-74 years of age and the 95-105 years of age groups (p = .038).

CONCLUSIONS

SyMRI is a new tool that allows faster imaging and permits to obtain quantitative T1 and T2. By defining RT values of different brain components of normal-functioning elderly individuals, this technique may be used as a biomarker for clinical disorders like dementia.

Predictive accuracy of T2-FLAIR mismatch sign for the IDH-mutant, 1p/19q noncodeleted low-grade glioma: An updated systematic review and meta-analysis

BACKGROUND

The T2-fluid-attenuated inversion recovery (FLAIR) mismatch sign, has been considered a highly specific imaging biomarker of IDH-mutant, 1p/19q noncodeleted low-grade glioma. This systematic review and meta-analysis aimed to evaluate the diagnostic performance of T2-FLAIR mismatch sign for prediction of a patient with IDH-mutant, 1p/19q noncodeleted low-grade glioma, and identify the causes responsible for the heterogeneity across the included studies.

METHODS

A systematic literature search in the Ovid-MEDLINE and EMBASE databases was performed for studies reporting the relevant topic before November 17, 2020. The pooled sensitivity and specificity values with their 95% confidence intervals were calculated using bivariate random-effects modeling. Meta-regression analyses were also performed to determine factors influencing heterogeneity.

RESULTS

For all the 10 included cohorts from 8 studies, the pooled sensitivity was 40% (95% confidence interval [CI] 28-53%), and the pooled specificity was 100% (95% CI 95-100%). In the hierarchic summary receiver operating characteristic curve, the difference between the 95% confidence and prediction regions was relatively large, indicating heterogeneity among the studies. Higgins I2 statistics demonstrated considerable heterogeneity in sensitivity (I2 = 83.5%) and considerable heterogeneity in specificity (I2 = 95.83%). Among the potential covariates, it seemed that none of factors was significantly associated with study heterogeneity in the joint model. However, the specificity was increased in studies with all the factors based on the differences in the composition of the detailed tumors.

CONCLUSIONS

The T2-FLAIR mismatch sign is near-perfect specific marker of IDH mutation and 1p/19q noncodeletion.

A low-cost and shielding-free ultra-low-field brain MRI scanner

Magnetic resonance imaging is a key diagnostic tool in modern healthcare, yet it can be cost-prohibitive given the high installation, maintenance and operation costs of the machinery. There are approximately seven scanners per million inhabitants and over 90% are concentrated in high-income countries. We describe an ultra-low-field brain MRI scanner that operates using a standard AC power outlet and is low cost to build. Using a permanent 0.055 Tesla Samarium-cobalt magnet and deep learning for cancellation of electromagnetic interference, it requires neither magnetic nor radiofrequency shielding cages. The scanner is compact, mobile, and acoustically quiet during scanning. We implement four standard clinical neuroimaging protocols (T1- and T2-weighted, fluid-attenuated inversion recovery like, and diffusion-weighted imaging) on this system, and demonstrate preliminary feasibility in diagnosing brain tumor and stroke. Such technology has the potential to meet clinical needs at point of care or in low and middle income countries.

Quantification of cerebral blood volume changes caused by visual stimulation at 3 T using DANTE-prepared dual-echo EPI

PURPOSE

We investigate the influence of moving blood-attenuation effects when using "delay alternating with nutation for tailored excitation" (DANTE) pulses in conjunction with blood oxygen level dependent (BOLD) of functional MRI (fMRI) at 3 T. Based on the effects of including DANTE pulses, we propose quantification of cerebral blood volume (CBV) changes following functional stimulation.

METHODS

Eighteen volunteers in total underwent fMRI scans at 3 T. Seven volunteers were scanned to investigate the effects of DANTE pulses on the fMRI signal. CBV changes in response to visual stimulation were quantified in 11 volunteers using a DANTE-prepared dual-echo EPI sequence.

RESULTS

The inflow effects from flowing blood in arteries and draining vein effects from flowing blood in large veins can be suppressed by use of a DANTE preparation module. Using DANTE-prepared dual-echo EPI, we quantitatively measured intravascular-weighted microvascular CBV changes of 25.4%, 29.8%, and 32.6% evoked by 1, 5, and 10 Hz visual stimulation, respectively. The extravascular fraction (∆S/S)_extra at TE = 30 ms in total BOLD signal was determined to be 64.8 ± 3.4%, which is in line with previous extravascular component estimation at 3 T. Results show that the microvascular CBV changes are linearly dependent on total BOLD changes at TE = 30 ms with a slope of 0.113, and this relation is independent of stimulation frequency and subject.

CONCLUSION

The DANTE preparation pulses can be incorporated into a standard EPI fMRI sequence for the purpose of minimizing inflow effects and reducing draining veins effects in large vessels. Additionally, the DANTE-prepared dual-echo EPI sequence is a promising fast imaging tool for quantification of intravascular-weighted CBV change in the microvascular space at 3 T.

Neurophysiological modulations in the (pre)motor-motor network underlying age-related increases in reaction time and the role of GABA levels - a bimodal TMS-MRS study

It has been argued that age-related changes in the neurochemical and neurophysiological properties of the GABAergic system may underlie increases in reaction time (RT) in older adults. However, the role of GABA levels within the sensorimotor cortices (SMC) in mediating interhemispheric interactions (IHi) during the processing stage of a fast motor response, as well as how both properties explain interindividual differences in RT, are not yet fully understood. In this study, edited magnetic resonance spectroscopy (MRS) was combined with dual-site transcranial magnetic stimulation (dsTMS) for probing GABA+ levels in bilateral SMC and task-related neurophysiological modulations in corticospinal excitability (CSE), and primary motor cortex (M1)-M1 and dorsal premotor cortex (PMd)-M1 IHi, respectively. Both CSE and IHi were assessed during the preparatory and premotor period of a delayed choice RT task. Data were collected from 25 young (aged 18-33 years) and 28 older (aged 60-74 years) healthy adults. Our results demonstrated that older as compared to younger adults exhibited a reduced bilateral CSE suppression, as well as a reduced magnitude of long latency M1-M1 and PMd-M1 disinhibition during the preparatory period, irrespective of the direction of the IHi. Importantly, in older adults, the GABA+ levels in bilateral SMC partially accounted for task-related neurophysiological modulations as well as individual differences in RT. In contrast, in young adults, neither task-related neurophysiological modulations, nor individual differences in RT were associated with SMC GABA+ levels. In conclusion, this study contributes to a comprehensive initial understanding of how age-related differences in neurochemical properties and neurophysiological processes are related to increases in RT.

Post mortem brain temperature and its influence on quantitative MRI of the brain

Objective

MRI temperature sensitivity presents a major issue in in situ post mortem MRI (PMMRI), as the tissue temperatures differ from living persons due to passive cooling of the deceased. This study aims at computing brain temperature effects on the MRI parameters to correct for temperature in PMMRI, laying the foundation for future projects on post mortem validation of in vivo MRI techniques.

Materials and methods

Brain MRI parameters were assessed in vivo and in situ post mortem using a 3 T MRI scanner. Post mortem brain temperature was measured in situ transethmoidally. The temperature effect was computed by fitting a linear model to the MRI parameters and the corresponding brain temperature.

Results

Linear positive temperature correlations were observed for T₁, T₂* and mean diffusivity in all tissue types. A significant negative correlation was observed for T₂ in white matter. Fractional anisotropy revealed significant correlations in all gray matter regions except for the thalamus.

Discussion

The linear models will allow to correct for temperature in post mortem MRI. Comparing in vivo to post mortem conditions, the mean diffusivity, in contrast to T₁ and T₂, revealed additional effects besides temperature, such as cessation of perfusion and active diffusion.

Separating spin compartments in arterial spin labeling using delays alternating with nutation for tailored excitation (DANTE) pulse: A validation study using T2 -relaxometry and application to arterial cerebral blood volume imaging

PURPOSE

To clarify the type of spin compartment in arterial spin labeling (ASL) that is eliminated by delays alternating with nutation for tailored excitation (DANTE) pulse using T₂ -relaxometry, and to demonstrate the feasibility of arterial cerebral blood volume (CBV_a ) imaging using DANTE-ASL in combination with a simplified two-compartment model.

METHOD

The DANTE and T₂ -preparation modules were combined into a single ASL sequence. T₂ values under the application of DANTE were determined to evaluate changes in T₂ , along with the post-labeling delay (PLD) and the relationship between transit time without DANTE (TT_noVS ) and T₂ . The reference tissue T₂ (T_{2_ref} ) was also obtained. Subsequently, the DANTE module was embedded into the Hadamard-encoded ASL. Cerebral blood flow (CBF) and CBV_a were computed using two Hadamard-encoding datasets (with and without DANTE) in a rest and breath-holding (BH) task.

RESULTS

While T₂ without DANTE (T_{2_noVS} ) decreased as the PLD increased, T₂ with DANTE (T_{2_DANTE} ) was equivalent to T_{2_ref} and did not change with the PLD. Although there was a significant positive correlation between TT_noVS and T_{2_noVS} with short PLD, T_{2_DANTE} was not correlated with TT_noVS nor PLD. Baseline CBV_a values obtained at rest were 0.64 ± 0.12, 0.64 ± 0.11, and 0.58 ± 0.15 mL/100 g for anterior, middle, and posterior cerebral arteries, respectively. Significant CBF and CBV_a elevations were observed in the BH task.

CONCLUSION

Microvascular compartment signals were eliminated from the total ASL signals by DANTE. CBV_a can be measured using Hadamard-encoded DANTE-ASL in combination with a simplified two-compartment model.

A deep learning approach for synthetic MRI based on two routine sequences and training with synthetic data

BACKGROUND AND OBJECTIVE

Synthetic magnetic resonance imaging (MRI) is a low cost procedure that serves as a bridge between qualitative and quantitative MRI. However, the proposed methods require very specific sequences or private protocols which have scarcely found integration in clinical scanners. We propose a learning-based approach to compute T1, T2, and PD parametric maps from only a pair of T1- and T2-weighted images customarily acquired in the clinical routine.

METHODS

Our approach is based on a convolutional neural network (CNN) trained with synthetic data; specifically, a synthetic dataset with 120 volumes was constructed from the anatomical brain model of the BrainWeb tool and served as the training set. The CNN learns an end-to-end mapping function to transform the input T1- and T2-weighted images to their underlying T1, T2, and PD parametric maps. Then, conventional weighted images unseen by the network are analytically synthesized from the parametric maps. The network can be fine tuned with a small database of actual weighted images and maps for better performance.

RESULTS

This approach is able to accurately compute parametric maps from synthetic brain data achieving normalized squared error values predominantly below 1%. It also yields realistic parametric maps from actual MR brain acquisitions with T1, T2, and PD values in the range of the literature and with correlation values above 0.95 compared to the T1 and T2 maps obtained from relaxometry sequences. Further, the synthesized weighted images are visually realistic; the mean square error values are always below 9% and the structural similarity index is usually above 0.90. Network fine tuning with actual maps improves performance, while training exclusively with a small database of actual maps shows a performance degradation.

CONCLUSIONS

These results show that our approach is able to provide realistic parametric maps and weighted images out of a CNN that (a) is trained with a synthetic dataset and (b) needs only two inputs, which are in turn obtained from a common full-brain acquisition that takes less than 8 min of scan time. Although a fine tuning with actual maps improves performance, synthetic data is crucial to reach acceptable performance levels. Hence, we show the utility of our approach for both quantitative MRI in clinical viable times and for the synthesis of additional weighted images to those actually acquired.

DeepVolume: Brain Structure and Spatial Connection-Aware Network for Brain MRI Super-Resolution

Thin-section magnetic resonance imaging (MRI) can provide higher resolution anatomical structures and more precise clinical information than thick-section images. However, thin-section MRI is not always available due to the imaging cost issue. In multicenter retrospective studies, a large number of data are often in thick-section manner with different section thickness. The lack of thin-section data and the difference in section thickness bring considerable difficulties in the study based on the image big data. In this article, we introduce DeepVolume, a two-step deep learning architecture to address the challenge of accurate thin-section MR image reconstruction. The first stage is the brain structure-aware network, in which the thick-section MR images in axial and sagittal planes are fused by a multitask 3-D U-net with prior knowledge of brain volume segmentation, which encourages the reconstruction result to have correct brain structure. The second stage is the spatial connection-aware network, in which the preliminary reconstruction results are adjusted slice-by-slice by a recurrent convolutional network embedding convolutional long short-term memory (LSTM) block, which enhances the precision of the reconstruction by utilizing the previously unassessed sagittal information. We used 305 paired brain MRI samples with thickness of 1.0 mm and 6.5 mm in this article. Extensive experiments illustrate that DeepVolume can produce the state-of-the-art reconstruction results by embedding more anatomical knowledge. Furthermore, considering DeepVolume as an intermediate step, the practical and clinical value of our method is validated by applying the brain volume estimation and voxel-based morphometry. The results show that DeepVolume can provide much more reliable brain volume estimation in the normalized space based on the thick-section MR images compared with the traditional solutions.

Characterization of the 1H-MRS Metabolite Spectra in Transgender Men with Gender Dysphoria and Cisgender People

Much research has been conducted on sexual differences of the human brain to determine whether and to what extent a brain gender exists. Consequently, a variety of studies using different neuroimaging techniques attempted to identify the existence of a brain phenotype in people with gender dysphoria (GD). However, to date, brain sexual differences at the metabolite level using magnetic resonance spectroscopy (1H-MRS) have not been explored in transgender people. In this study, 28 cisgender men (CM) and 34 cisgender women (CW) and 29 transgender men with GD (TMGD) underwent 1H-MRS at 3 Tesla MRI to characterize common brain metabolites. Specifically, levels of N-acetyl aspartate (NAA), choline (Cho), creatine (Cr), glutamate and glutamine (Glx), and myo-inositol + glycine (mI + Gly) were assessed in two brain regions, the amygdala-anterior hippocampus and the lateral parietal cortex. The results indicated a sex-assigned at birth pattern for Cho/Cr in the amygdala of TMGD. In the parietal cortex, a sex-assigned at birth and an intermediate pattern were found. Though assessed post-hoc, exploration of the age of onset of GD in TMGD demonstrated within-group differences in absolute NAA and relative Cho/Cr levels, suggestive for a possible developmental trend. While brain metabolite levels in TMGD resembled those of CW, some interesting findings, such as modulation of metabolite concentrations by age of onset of GD, warrant future inquiry.

Three-dimensional whole-brain mapping of cerebral blood volume and venous cerebral blood volume using Fourier transform-based velocity-selective pulse trains

PURPOSE

To develop 3D MRI methods for cerebral blood volume (CBV) and venous cerebral blood volume (vCBV) estimation with whole-brain coverage using Fourier transform-based velocity-selective (FT-VS) pulse trains.

METHODS

For CBV measurement, FT-VS saturation pulse trains were used to suppress static tissue, whereas CSF contamination was corrected voxel-by-voxel using a multi-readout acquisition and a fast CSF T₂ scan. The vCBV mapping was achieved by inserting an arterial-nulling module that included a FT-VS inversion pulse train. Using these methods, CBV and vCBV maps were obtained on 6 healthy volunteers at 3 T.

RESULTS

The mean CBV and vCBV values in gray matter and white matter in different areas of the brain showed high correlation (r = 0.95 and P < .0001). The averaged CBV and vCBV values of the whole brain were 5.4 ± 0.6 mL/100 g and 2.5 ± 0.3 mL/100 g in gray matter, and 2.6 ± 0.5 mL/100 g and 1.5 ± 0.2 mL/100 g in white matter, respectively, comparable to the literature.

CONCLUSION

The feasibility of FT-VS-based CBV and vCBV estimation was demonstrated for 3D acquisition with large spatial coverage.