Body size interacts with the structure of the central nervous system: A multi-center in vivo neuroimaging study.

Clinical research emphasizes the implementation of rigorous and reproducible study designs that rely on between-group matching or controlling for sources of biological variation such as subject's sex and age. However, corrections for body size (i.e. height and weight) are mostly lacking in clinical neuroimaging designs. This study investigates the importance of body size parameters in their relationship with spinal cord (SC) and brain magnetic resonance imaging (MRI) metrics. Data were derived from a cosmopolitan population of 267 healthy human adults (age 30.1±6.6 years old, 125 females). We show that body height correlated strongly or moderately with brain gray matter (GM) volume, cortical GM volume, total cerebellar volume, brainstem volume, and cross-sectional area (CSA) of cervical SC white matter (CSA-WM; 0.44≤r≤0.62). In comparison, age correlated weakly with cortical GM volume, precentral GM volume, and cortical thickness (-0.21≥r≥-0.27). Body weight correlated weakly with magnetization transfer ratio in the SC WM, dorsal columns, and lateral corticospinal tracts (-0.20≥r≥-0.23). Body weight further correlated weakly with the mean diffusivity derived from diffusion tensor imaging (DTI) in SC WM (r=-0.20) and dorsal columns (-0.21), but only in males. CSA-WM correlated strongly or moderately with brain volumes (0.39≤r≤0.64), and weakly with precentral gyrus thickness and DTI-based fractional anisotropy in SC dorsal columns and SC lateral corticospinal tracts (-0.22≥r≥-0.25). Linear mixture of sex and age explained 26±10% of data variance in brain volumetry and SC CSA. The amount of explained variance increased at 33±11% when body height was added into the mixture model. Age itself explained only 2±2% of such variance. In conclusion, body size is a significant biological variable. Along with sex and age, body size should therefore be included as a mandatory variable in the design of clinical neuroimaging studies examining SC and brain structure.

Nonaqueous magnetization following adiabatic and selective pulses in brain: T1 and cross-relaxation dynamics.

Inversion pulses are commonly employed in MRI for T 1 -weighted contrast and relaxation measurements. In the brain, it is often assumed that adiabatic pulses saturate the nonaqueous magnetization. We investigated this assumption using solid-state NMR to monitor the nonaqueous signal directly following adiabatic inversion and compared this with signals following hard and soft inversion pulses. The effects of the different preparations on relaxation dynamics were explored. Inversion recovery experiments were performed on ex vivo bovine and porcine brains using 360-MHz (8.4 T) and 200-MHz (4.7 T) NMR spectrometers, respectively, using broadband rectangular, adiabatic, and sinc inversion pulses as well as a long rectangular saturation pulse. Analogous human brain MRI experiments were performed at 3 T using single-slice echo-planar imaging. Relaxation data were fitted by mono- and biexponential decay models. Further fitting analysis was performed using only two inversion delay times. Adiabatic and sinc inversion left much of the nonaqueous magnetization along B 0 and resulted in biexponential relaxation. Saturation of both aqueous and nonaqueous magnetization components led to effectively monoexponential T 1 relaxation. Typical adiabatic inversion pulses do not, as has been widely assumed, saturate the nonaqueous proton magnetization in white matter. Unequal magnetization states in aqueous and nonaqueous 1 H reservoirs prepared by soft and adiabatic pulses result in biexponential T 1 relaxation. Both pools must be prepared in the same magnetization state (e.g., saturated or inverted) in order to observe consistent monoexponential relaxation.

Predictive Value of Acute Neuroplastic Response to rTMS in Treatment Outcome in Depression: A Concurrent TMS-fMRI Trial.

OBJECTIVE: The study objective was to investigate the predictive value of functional connectivity changes induced by acute repetitive transcranial magnetic stimulation (rTMS) for clinical response in treatment-resistant depression. METHODS: Cross-sectional changes in functional connectivity induced by a single concurrent rTMS-fMRI session were assessed in 38 outpatients with treatment-resistant depression (26 of them female; mean age, 41.87 years) who subsequently underwent a 4-week course of rTMS. rTMS was delivered at 1 Hz over the right dorsolateral prefrontal cortex. Acute rTMS-induced functional connectivity changes were computed and subjected to connectome-based predictive modeling to test their association with changes in score on the Montgomery-Åsberg Depression Rating Scale (MADRS) after rTMS treatment. RESULTS: TMS-fMRI induced widespread, acute, and transient alterations in functional connectivity. The rTMS-induced connectivity changes predicted about 30% of the variance of improvement in the MADRS score. The most robust predictive associations involved connections between prefrontal regions and motor, parietal, and insular cortices and between bilateral regions of the thalamus. CONCLUSIONS: Acute rTMS-induced connectivity changes in patients with treatment-resistant depression may index macro-level neuroplasticity, relevant to interindividual variability in rTMS treatment response. Large-scale network phenomena occurring during rTMS might be used to inform prospective clinical trials.

Open-access quantitative MRI data of the spinal cord and reproducibility across participants, sites and manufacturers.

In a companion paper by Cohen-Adad et al. we introduce the spine generic quantitative MRI protocol that provides valuable metrics for assessing spinal cord macrostructural and microstructural integrity. This protocol was used to acquire a single subject dataset across 19 centers and a multi-subject dataset across 42 centers (for a total of 260 participants), spanning the three main MRI manufacturers: GE, Philips and Siemens. Both datasets are publicly available via git-annex. Data were analysed using the Spinal Cord Toolbox to produce normative values as well as inter/intra-site and inter/intra-manufacturer statistics. Reproducibility for the spine generic protocol was high across sites and manufacturers, with an average inter-site coefficient of variation of less than 5% for all the metrics. Full documentation and results can be found at https://spine-generic.rtfd.io/ . The datasets and analysis pipeline will help pave the way towards accessible and reproducible quantitative MRI in the spinal cord.

Generic acquisition protocol for quantitative MRI of the spinal cord.

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols . The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition.

Frequency drift in MR spectroscopy at 3T.

PURPOSE: Heating of gradient coils and passive shim components is a common cause of instability in the B0 field, especially when gradient intensive sequences are used. The aim of the study was to set a benchmark for typical drift encountered during MR spectroscopy (MRS) to assess the need for real-time field-frequency locking on MRI scanners by comparing field drift data from a large number of sites. METHOD: A standardized protocol was developed for 80 participating sites using 99 3T MR scanners from 3 major vendors. Phantom water signals were acquired before and after an EPI sequence. The protocol consisted of: minimal preparatory imaging; a short pre-fMRI PRESS; a ten-minute fMRI acquisition; and a long post-fMRI PRESS acquisition. Both pre- and post-fMRI PRESS were non-water suppressed. Real-time frequency stabilization/adjustment was switched off when appropriate. Sixty scanners repeated the protocol for a second dataset. In addition, a three-hour post-fMRI MRS acquisition was performed at one site to observe change of gradient temperature and drift rate. Spectral analysis was performed using MATLAB. Frequency drift in pre-fMRI PRESS data were compared with the first 5:20 minutes and the full 30:00 minutes of data after fMRI. Median (interquartile range) drifts were measured and showed in violin plot. Paired t-tests were performed to compare frequency drift pre- and post-fMRI. A simulated in vivo spectrum was generated using FID-A to visualize the effect of the observed frequency drifts. The simulated spectrum was convolved with the frequency trace for the most extreme cases. Impacts of frequency drifts on NAA and GABA were also simulated as a function of linear drift. Data from the repeated protocol were compared with the corresponding first dataset using Pearson's and intraclass correlation coefficients (ICC). RESULTS: Of the data collected from 99 scanners, 4 were excluded due to various reasons. Thus, data from 95 scanners were ultimately analyzed. For the first 5:20 min (64 transients), median (interquartile range) drift was 0.44 (1.29) Hz before fMRI and 0.83 (1.29) Hz after. This increased to 3.15 (4.02) Hz for the full 30 min (360 transients) run. Average drift rates were 0.29 Hz/min before fMRI and 0.43 Hz/min after. Paired t-tests indicated that drift increased after fMRI, as expected (p < 0.05). Simulated spectra convolved with the frequency drift showed that the intensity of the NAA singlet was reduced by up to 26%, 44 % and 18% for GE, Philips and Siemens scanners after fMRI, respectively. ICCs indicated good agreement between datasets acquired on separate days. The single site long acquisition showed drift rate was reduced to 0.03 Hz/min approximately three hours after fMRI. DISCUSSION: This study analyzed frequency drift data from 95 3T MRI scanners. Median levels of drift were relatively low (5-min average under 1 Hz), but the most extreme cases suffered from higher levels of drift. The extent of drift varied across scanners which both linear and nonlinear drifts were observed.

Cross-validation study between the HRRT and the PET component of the SIGNA PET/MRI system with focus on neuroimaging.

BACKGROUND: The Siemens high-resolution research tomograph (HRRT - a dedicated brain PET scanner) is to this day one of the highest resolution PET scanners; thus, it can serve as useful benchmark when evaluating performance of newer scanners. Here, we report results from a cross-validation study between the HRRT and the whole-body GE SIGNA PET/MR focusing on brain imaging. Phantom data were acquired to determine recovery coefficients (RCs), % background variability (%BG), and image voxel noise (%). Cross-validation studies were performed with six healthy volunteers using [11C]DTBZ, [11C]raclopride, and [18F]FDG. Line profiles, regional time-activity curves, regional non-displaceable binding potentials (BPND) for [11C]DTBZ and [11C]raclopride scans, and radioactivity ratios for [18F]FDG scans were calculated and compared between the HRRT and the SIGNA PET/MR. RESULTS: Phantom data showed that the PET/MR images reconstructed with an ordered subset expectation maximization (OSEM) algorithm with time-of-flight (TOF) and TOF + point spread function (PSF) + filter revealed similar RCs for the hot spheres compared to those obtained on the HRRT reconstructed with an ordinary Poisson-OSEM algorithm with PSF and PSF + filter. The PET/MR TOF + PSF reconstruction revealed the highest RCs for all hot spheres. Image voxel noise of the PET/MR system was significantly lower. Line profiles revealed excellent spatial agreement between the two systems. BPND values revealed variability of less than 10% for the [11C]DTBZ scans and 19% for [11C]raclopride (based on one subject only). Mean [18F]FDG ratios to pons showed less than 12% differences. CONCLUSIONS: These results demonstrated comparable performances of the two systems in terms of RCs with lower voxel-level noise (%) present in the PET/MR system. Comparison of in vivo human data confirmed the comparability of the two systems. The whole-body GE SIGNA PET/MR system is well suited for high-resolution brain imaging as no significant performance degradation was found compared to that of the reference standard HRRT.

An International Standardized Magnetic Resonance Imaging Protocol for Diagnosis and Follow-up of Patients with Multiple Sclerosis: Advocacy, Dissemination, and Implementation Strategies.

Standardized magnetic resonance imaging (MRI) protocols are important for the diagnosis and monitoring of patients with multiple sclerosis (MS). The Consortium of Multiple Sclerosis Centers (CMSC) convened an international panel of MRI experts to review and update the current guidelines. The objective was to update the standardized MRI protocol and clinical guidelines for diagnosis and follow-up of MS and develop strategies for advocacy, dissemination, and implementation. Conference attendees included neurologists, radiologists, technologists, and imaging scientists with expertise in MS. Representatives from the CMSC, Magnetic Resonance Imaging in MS (MAGNIMS), North American Imaging in Multiple Sclerosis Cooperative, US Department of Veteran Affairs, National Multiple Sclerosis Society, Multiple Sclerosis Association of America, MRI manufacturers, and commercial image analysis companies were present. Before the meeting, CMSC members were surveyed about standardized MRI protocols, gadolinium use, need for diffusion-weighted imaging, and the central vein sign. The panel worked to make the CMSC and MAGNIMS MRI protocols similar so that the updated guidelines could ultimately be accepted by international consensus. Advocacy efforts will promote the importance of standardized MS MRI protocols. Dissemination will include publications, meeting abstracts, educational programming, webinars, "meet the expert" teleconferences, and examination cards. Implementation will require comprehensive and coordinated efforts to make the protocol easy to access and use. The ultimate vision, and goal, is for the guidelines to be universally useful, usable, and used as the standard of care for patients with MS.

Multi-spin echo T2 relaxation imaging with compressed sensing (METRICS) for rapid myelin water imaging.

PURPOSE: Myelin water imaging (MWI) provides a valuable biomarker for myelin, but clinical application has been restricted by long acquisition times. Accelerating the standard multi-echo T2 acquisition with gradient echoes (GRASE) or by 2D multi-slice data collection results in image blurring, contrast changes, and other issues. Compressed sensing (CS) can vastly accelerate conventional MRI. In this work, we assessed the use of CS for in vivo human MWI, using a 3D multi spin-echo sequence. METHODS: We implemented multi-echo T2 relaxation imaging with compressed sensing (METRICS) and METRICS with partial Fourier acceleration (METRICS-PF). Scan-rescan data were acquired from 12 healthy controls for assessment of repeatability. MWI data were acquired for METRICS in 9 m:58 s and for METRICS-PF in 7 m:25 s, both with 1.5 × 2 × 3 mm3 voxels, 56 echoes, 7 ms ΔTE, and 240 × 240 × 170 mm3 FOV. METRICS was compared with a novel multi-echo spin-echo gold-standard (MSE-GS) MWI acquisition, acquired for a single additional subject in 2 h:2 m:40 s. RESULTS: METRICS/METRICS-PF myelin water fraction had mean: repeatability coefficient 1.5/1.1, coefficient of variation 6.2/4.5%, and intra-class correlation coefficient 0.79/0.84. Repeatability metrics comparing METRICS with METRICS-PF were similar, and both sequences agreed with reference values from literature. METRICS images and quantitative maps showed excellent qualitative agreement with those of MSE-GS. CONCLUSION: METRICS and METRICS-PF provided highly repeatable MWI data without the inherent disadvantages of GRASE or 2D multi-slice acquisition. CS acceleration allows MWI data to be acquired rapidly with larger FOV, higher estimated SNR, more isotropic voxels and more echoes than with previous techniques. The approach introduced here generalizes to any multi-component T2 mapping application.

Myelin Water Fraction and Intra/Extracellular Water Geometric Mean T2 Normative Atlases for the Cervical Spinal Cord from 3T MRI.

BACKGROUND AND PURPOSE: Acquiring and interpreting quantitative myelin-specific MRI data at an individual level is challenging because of technical difficulties and natural myelin variation in the population. To overcome these challenges, we used multiecho T2 myelin water imaging (MWI) to create T2 metric healthy population atlases that depict the mean and variation of myelin water fraction (MWF), and intra- and extracellular water mobility as described by geometric mean T2 (IEGMT2 ). METHODS: Cervical cord MWI was performed at 3T on 20 healthy individuals (10M/10F, mean age: 36 years) and 3 relapsing remitting multiple sclerosis (RRMS) participants (1M/2F, age: 39/42/37 years). Anatomical data were collected for the purpose of image segmentation and registration. Atlases were created by coregistering and averaging T2 metrics from all controls. Voxel-wise z-score maps from 3 RRMS participants were produced to demonstrate the preliminary utility of the MWF and IEGMT2 atlases. RESULTS: The average MWF atlas provides a representation of myelin in the spinal cord consistent with well-known spinal cord anatomical characteristics. The IEGMT2 atlas also depicted structural variations in the spinal cord. Z-score analysis illustrated distinct abnormalities in MWF and IEGMT2 in the 3 RRMS cases. CONCLUSIONS: Our findings highlight the potential for using a quantitative T2 relaxation metric atlas to visualize and detect pathology in spinal cord. Our MWF and IEGMT2 atlases (URL: https://sourceforge.net/projects/mwi-spinal-cord-atlases/) can serve as normative references in the cervical spinal cord for other studies.

Rapid myelin water imaging for the assessment of cervical spinal cord myelin damage.

BACKGROUND: Rapid myelin water imaging (MWI) using a combined gradient and spin echo (GRASE) sequence can produce myelin specific metrics for the human brain. Spinal cord MWI could be similarly useful, but technical challenges have hindered routine application. GRASE rapid MWI was recently successfully implemented for imaging of healthy cervical spinal cord and may complement other advanced imaging methods, such as diffusion tensor imaging (DTI) and quantitative T1 (qT1). OBJECTIVE: To demonstrate the feasibility of cervical cord GRASE rapid MWI in multiple sclerosis (MS), primary lateral sclerosis (PLS) and neuromyelitis optica spectrum disorder (NMO), with comparison to DTI and qT1 metrics. METHODS: GRASE MWI, DTI and qT1 data were acquired in 2 PLS, 1 relapsing-remitting MS (RRMS), 1 primary-progressive MS (PPMS) and 2 NMO subjects, as well as 6 age (±3 yrs) and sex matched healthy controls (HC). Internal cord structure guided template registrations, used for region of interest (ROI) analysis. Z score maps were calculated for the difference between disease subject and mean HC metric values. RESULTS: PLS subjects had low myelin water fraction (MWF) in the lateral funiculi compared to HC. RRMS subject MWF was heterogeneous within the cord. The PPMS subject showed no trends in ROI results but had a region of low MWF Z score corresponding to a focal lesion. The NMO subject with a longitudinally extensive transverse myelitis lesion had low values for whole cord mean MWF of 12.8% compared to 24.3% (standard deviation 2.2%) for HC. The NMO subject without lesions also had low MWF compared to HC. DTI and qT1 metrics showed similar trends, corroborating the MWF results and providing complementary information. CONCLUSION: GRASE is sufficiently sensitive to detect decreased myelin within MS spinal cord plaques, NMO lesions, and PLS diffuse spinal cord injury. Decreased MWF in PLS is consistent with demyelination secondary to motor neuron degeneration. GRASE MWI is a feasible method for rapid assessment of myelin content in the cervical spinal cord and provides complementary information to that of DTI and qT1 measures.

The Killer Fly Hunger Games: Target Size and Speed Predict Decision to Pursuit.

Predatory animals have evolved to optimally detect their prey using exquisite sensory systems such as vision, olfaction and hearing. It may not be so surprising that vertebrates, with large central nervous systems, excel at predatory behaviors. More striking is the fact that many tiny insects, with their miniscule brains and scaled down nerve cords, are also ferocious, highly successful predators. For predation, it is important to determine whether a prey is suitable before initiating pursuit. This is paramount since pursuing a prey that is too large to capture, subdue or dispatch will generate a substantial metabolic cost (in the form of muscle output) without any chance of metabolic gain (in the form of food). In addition, during all pursuits, the predator breaks its potential camouflage and thus runs the risk of becoming prey itself. Many insects use their eyes to initially detect and subsequently pursue prey. Dragonflies, which are extremely efficient predators, therefore have huge eyes with relatively high spatial resolution that allow efficient prey size estimation before initiating pursuit. However, much smaller insects, such as killer flies, also visualize and successfully pursue prey. This is an impressive behavior since the small size of the killer fly naturally limits the neural capacity and also the spatial resolution provided by the compound eye. Despite this, we here show that killer flies efficiently pursue natural (Drosophila melanogaster) and artificial (beads) prey. The natural pursuits are initiated at a distance of 7.9 ± 2.9 cm, which we show is too far away to allow for distance estimation using binocular disparities. Moreover, we show that rather than estimating absolute prey size prior to launching the attack, as dragonflies do, killer flies attack with high probability when the ratio of the prey's subtended retinal velocity and retinal size is 0.37. We also show that killer flies will respond to a stimulus of an angular size that is smaller than that of the photoreceptor acceptance angle, and that the predatory response is strongly modulated by the metabolic state. Our data thus provide an exciting example of a loosely designed matched filter to Drosophila, but one which will still generate successful pursuits of other suitable prey.