Normative volumes and relaxation times at 3T during brain development

While research has unveiled and quantified brain markers of abnormal neurodevelopment, clinicians still work with qualitative metrics for MRI brain investigation. The purpose of the current article is to bridge the knowledge gap between case-control cohort studies and individual patient care. Here, we provide a unique dataset of seventy-three 3-to-17 years-old healthy subjects acquired with a 6-minute MRI protocol encompassing T1 and T2 relaxation quantitative sequence that can be readily implemented in the clinical setting; MP2RAGE for T1 mapping and the prototype sequence GRAPPATINI for T2 mapping. White matter and grey matter volumes were automatically quantified. We further provide normative developmental curves based on these two imaging sequences; T1, T2 and volume normative ranges with respect to age were computed, for each ROI of a pediatric brain atlas. This open-source dataset provides normative values allowing to position individual patients acquired with the same protocol on the brain maturation curve and as such provides potentially useful quantitative biomarkers facilitating precise and personalized care.

Tract-wise microstructural analysis informs on current and future disability in early multiple sclerosis

OBJECTIVES

Microstructural characterization of patients with multiple sclerosis (MS) has been shown to correlate better with disability compared to conventional radiological biomarkers. Quantitative MRI provides effective means to characterize microstructural brain tissue changes both in lesions and normal-appearing brain tissue. However, the impact of the location of microstructural alterations in terms of neuronal pathways has not been thoroughly explored so far. Here, we study the extent and the location of tissue changes probed using quantitative MRI along white matter (WM) tracts extracted from a connectivity atlas.

METHODS

We quantified voxel-wise T1 tissue alterations compared to normative values in a cohort of 99 MS patients. For each WM tract, we extracted metrics reflecting tissue alterations both in lesions and normal-appearing WM and correlated these with cross-sectional disability and disability evolution after 2 years.

RESULTS

In early MS patients, T1 alterations in normal-appearing WM correlated better with disability evolution compared to cross-sectional disability. Further, the presence of lesions in supratentorial tracts was more strongly associated with cross-sectional disability, while microstructural alterations in infratentorial pathways yielded higher correlations with disability evolution. In progressive patients, all major WM pathways contributed similarly to explaining disability, and correlations with disability evolution were generally poor.

CONCLUSIONS

We showed that microstructural changes evaluated in specific WM pathways contribute to explaining future disability in early MS, hence highlighting the potential of tract-wise analyses in monitoring disease progression. Further, the proposed technique allows to estimate WM tract-specific microstructural characteristics in clinically compatible acquisition times, without the need for advanced diffusion imaging.

Precision Brain Morphometry Using Cluster Scanning

Measurement error limits the statistical power to detect group differences and longitudinal change in structural MRI morphometric measures (e.g., hippocampal volume, prefrontal thickness). Recent advances in scan acceleration enable extremely fast T1-weighted scans (~1 minute) to achieve morphometric errors that are close to the errors in longer traditional scans. As acceleration allows multiple scans to be acquired in rapid succession, it becomes possible to pool estimates to increase measurement precision, a strategy known as "cluster scanning." Here we explored brain morphometry using cluster scanning in a test-retest study of 40 individuals (12 younger adults, 18 cognitively unimpaired older adults, and 10 adults diagnosed with mild cognitive impairment or Alzheimer's Dementia). Morphometric errors from a single compressed sensing (CS) 1.0mm scan with 6x acceleration (CSx6) were, on average, 12% larger than a traditional scan using the Alzheimer's Disease Neuroimaging Initiative (ADNI) protocol. Pooled estimates from four clustered CSx6 acquisitions led to errors that were 34% smaller than ADNI despite having a shorter total acquisition time. Given a fixed amount of time, a gain in measurement precision can thus be achieved by acquiring multiple rapid scans instead of a single traditional scan. Errors were further reduced when estimates were pooled from eight CSx6 scans (51% smaller than ADNI). Neither pooling across a break nor pooling across multiple scan resolutions boosted this benefit. We discuss the potential of cluster scanning to improve morphometric precision, boost statistical power, and produce more sensitive disease progression biomarkers.

Implementation of a piritramide based patient-controlled analgesia (PCA) as a standard of care for pain control in late abortion induction: A prospective cohort study from a patient perspective

Objective

To assess whether the implementation of patient-controlled analgesia (PCA) with piritramide using an automatic pump system under routine conditions is effective to reduce pain in late abortion inductions.

Study design

Prospective observational cohort study.

Setting

Patients requiring medically indicated abortion induction from 14 weeks of pregnancy onwards between July 2019 and July 2020 at the department of Obstetrics and Prenatal Medicine of the Bonn University Hospital in Germany.

Methods

Evaluation of pain management after implementation of a PCA system compared with previous nurse-controlled tramadol-based standard under routine conditions. Patients answered a validated pain questionnaire and requirement of rescue analgesics was assessed. Pain intensity and satisfaction were measured on a ten-point numeric rating scale. Main Outcome Measure Maximal pain intensity.

Results

Forty patients were included. Patients using Piritramide-PCA complained of higher pain sores than those in the standard group (6.90 (± 2.34) vs. 4.83 (± 2.87), (p < 0.05)). In both groups the level of satisfaction with the analgesia received was comparable (8.00 (± 2.45) vs 7.67 (± 2.62), (p = 0.7)). Patients in the PCA group suffered more nausea (63.2 % vs 30 % respectively, OR 4.0, 95 % CI 1.05-15.20, p < 0.05) and expressed more the desire for more analgesic support compared to the control group (OR 5.7 (1-33.25), p = 0.05).

Conclusion

Women with abortion induction after 14 weeks of gestation suffer from relevant severe pain, which requires adequate therapy. However, addition of PCA does not seem to bring any advantage in patients undergoing this procedure.

SPARCQ: A new approach for fat fraction mapping using asymmetries in the phase-cycled balanced SSFP signal profile

PURPOSE

To develop SPARCQ (Signal Profile Asymmetries for Rapid Compartment Quantification), a novel approach to quantify fat fraction (FF) using asymmetries in the phase-cycled balanced SSFP (bSSFP) profile.

METHODS

SPARCQ uses phase-cycling to obtain bSSFP frequency profiles, which display asymmetries in the presence of fat and water at certain TRs. For each voxel, the measured signal profile is decomposed into a weighted sum of simulated profiles via multi-compartment dictionary matching. Each dictionary entry represents a single-compartment bSSFP profile with a specific off-resonance frequency and relaxation time ratio. Using the results of dictionary matching, the fractions of the different off-resonance components are extracted for each voxel, generating quantitative maps of water and FF and banding-artifact-free images for the entire image volume. SPARCQ was validated using simulations, experiments in a water-fat phantom and in knees of healthy volunteers. Experimental results were compared with reference proton density FFs obtained with 1 H-MRS (phantoms) and with multiecho gradient-echo MRI (phantoms and volunteers). SPARCQ repeatability was evaluated in six scan-rescan experiments.

RESULTS

Simulations showed that FF quantification is accurate and robust for SNRs greater than 20. Phantom experiments demonstrated good agreement between SPARCQ and gold standard FFs. In volunteers, banding-artifact-free quantitative maps and water-fat-separated images obtained with SPARCQ and ME-GRE demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ FF was 0.0512.

CONCLUSION

SPARCQ demonstrates potential for fat quantification using asymmetries in bSSFP profiles and may be a promising alternative to conventional FF quantification techniques.

Morphological assessment of cartilage and osteoarthritis in clinical practice and research: Intermediate-weighted fat-suppressed sequences and beyond

Magnetic resonance imaging (MRI) is widely regarded as the primary modality for the morphological assessment of cartilage and all other joint tissues involved in osteoarthritis. 2D fast spin echo fat-suppressed intermediate-weighted (FSE FS IW) sequences with a TE between 30 and 40ms have stood the test of time and are considered the cornerstone of MRI protocols for clinical practice and trials. These sequences offer a good balance between sensitivity and specificity and provide appropriate contrast and signal within the cartilage as well as between cartilage, articular fluid, and subchondral bone. Additionally, FS IW sequences enable the evaluation of menisci, ligaments, synovitis/effusion, and bone marrow edema-like signal changes. This review article provides a rationale for the use of FSE FS IW sequences in the morphological assessment of cartilage and osteoarthritis, along with a brief overview of other clinically available sequences for this indication. Additionally, the article highlights ongoing research efforts aimed at improving FSE FS IW sequences through 3D acquisitions with enhanced resolution, shortened examination times, and exploring the potential benefits of different magnetic field strengths. While most of the literature on cartilage imaging focuses on the knee, the concepts presented here are applicable to all joints. KEY POINTS: 1. MRI is currently considered the modality of reference for a "whole-joint" morphological assessment of osteoarthritis. 2. Fat-suppressed intermediate-weighted sequences remain the keystone of MRI protocols for the assessment of cartilage morphology, as well as other structures involved in osteoarthritis. 3. Trends for further development in the field of cartilage and joint imaging include 3D FSE imaging, faster acquisition including AI-based acceleration, and synthetic imaging providing multi-contrast sequences.

Brain morphometry in older adults with and without dementia using extremely rapid structural scans

T1-weighted structural MRI is widely used to measure brain morphometry (e.g., cortical thickness and subcortical volumes). Accelerated scans as fast as one minute or less are now available but it is unclear if they are adequate for quantitative morphometry. Here we compared the measurement properties of a widely adopted 1.0 mm resolution scan from the Alzheimer's Disease Neuroimaging Initiative (ADNI = 5'12'') with two variants of highly accelerated 1.0 mm scans (compressed-sensing, CSx6 = 1'12''; and wave-controlled aliasing in parallel imaging, WAVEx9 = 1'09'') in a test-retest study of 37 older adults aged 54 to 86 (including 19 individuals diagnosed with a neurodegenerative dementia). Rapid scans produced highly reliable morphometric measures that largely matched the quality of morphometrics derived from the ADNI scan. Regions of lower reliability and relative divergence between ADNI and rapid scan alternatives tended to occur in midline regions and regions with susceptibility-induced artifacts. Critically, the rapid scans yielded morphometric measures similar to the ADNI scan in regions of high atrophy. The results converge to suggest that, for many current uses, extremely rapid scans can replace longer scans. As a final test, we explored the possibility of a 0'49'' 1.2 mm CSx6 structural scan, which also showed promise. Rapid structural scans may benefit MRI studies by shortening the scan session and reducing cost, minimizing opportunity for movement, creating room for additional scan sequences, and allowing for the repetition of structural scans to increase precision of the estimates.

Feasibility of accelerated T2 mapping for the preoperative assessment of endometrial carcinoma

Objective

The application value of T2 mapping in evaluating endometrial carcinoma (EMC) features remains unclear. The aim of the study was to determine the quantitative T2 values in EMC using a novel accelerated T2 mapping, and evaluate them for detection, classification,and grading of EMC.

Materials and methods

Fifty-six patients with pathologically confirmed EMC and 17 healthy volunteers were prospectively enrolled in this study. All participants underwent pelvic magnetic resonance imaging, including DWI and accelerated T2 mapping, before treatment. The T2 and apparent diffusion coefficient (ADC) values of different pathologic EMC features were extracted and compared. Receiver operating characteristic (ROC) curve analysis was performed to analyze the diagnostic efficacy of the T2 and ADC values in distinguishing different pathological features of EMC.

Results

The T2 values and ADC values were significantly lower in EMC than in normal endometrium (bothl p < 0.05). The T2 and ADC values were significantly different between endometrioid adenocarcinoma (EA) and non-EA (both p < 0.05) and EMC tumor grades (all p < 0.05) but not for EMC clinical types (both p > 0.05) and depth of myometrial invasion (both p > 0.05). The area under the ROC curve (AUC) was higher for T2 values than for ADC values in predicting grade 3 EA (0.939 vs. 0.764, p = 0.048). When combined T2 and ADC values, the AUC for predicting grade 3 EA showed a significant increase to 0.947 (p = 0.03) compared with those of ADC values. The T2 and ADC values were negatively correlated with the tumor grades (r = -0.706 and r = -0.537, respectively).

Conclusion

Quantitative T2 values demonstrate potential suitability in discriminating between EMC and normal endometrium, EA and non-EA, grade 3 EA and grade 1/2 EA. Combining T2 and ADC values performs better in predicting the histological grades of EA in comparison with ADC values alone.

Hi-Fi fMRI: High-resolution, fast-sampled and sub-second whole-brain functional MRI at 3T in humans

Functional magnetic resonance imaging (fMRI) is a methodological cornerstone of neuroscience. Most studies measure blood-oxygen-level-dependent (BOLD) signal using echo-planar imaging (EPI), Cartesian sampling, and image reconstruction with a one-to-one correspondence between the number of acquired volumes and reconstructed images. However, EPI schemes are subject to trade-offs between spatial and temporal resolutions. We overcome these limitations by measuring BOLD with a gradient recalled echo (GRE) with 3D radial-spiral phyllotaxis trajectory at a high sampling rate (28.24ms) on standard 3T field-strength. The framework enables the reconstruction of 3D signal time courses with whole-brain coverage at simultaneously higher spatial (1mm ³ ) and temporal (up to 250ms) resolutions, as compared to optimized EPI schemes. Additionally, artifacts are corrected before image reconstruction; the desired temporal resolution is chosen after scanning and without assumptions on the shape of the hemodynamic response. By showing activation in the calcarine sulcus of 20 participants performing an ON-OFF visual paradigm, we demonstrate the reliability of our method for cognitive neuroscience research.

Clinical feasibility and validation of the accelerated T2 mapping sequence GRAPPATINI in brain imaging

Rationale and objectives

To prospectively evaluate feasibility and robustness of an accelerated T2 mapping sequence (GRAPPATINI) in brain imaging and to assess its synthetic T2-weighted images (sT2w) in comparison with a standard T2-weighted sequence (T2 TSE).

Material and methods

Volunteers were included to evaluate the robustness and consecutive patients for morphological evaluation. They were scanned on a 3 T MR-scanner. Healthy volunteers underwent GRAPPATINI of the brain three times (day 1: scan/rescan; day 2: follow-up). Patients between the ages of 18 and 85 years who were able to provide written informed consent and who had no MRI contraindications were included. For morphological comparison two radiologists with 5 and 7 years of experience in brain MRI evaluated image quality using a Likert scale (1 being poor, 4 being excellent) in a blinded and randomized fashion.

Results

Images were successfully acquired in ten volunteers with a mean age of 25 years (ranging from 22 to 31 years) and 52 patients (23 men/29 women) with a mean age of 55 years (range of 22-83 years). Most brain regions showed repeatable and reproducible T2 values (rescan: CoV 0.75%-2.06%, ICC 69%-92.3%; follow-up: CoV 0.41%-1.59%, ICC 79.4%-95.8%), except for the caudate nucleus (rescan: CoV 7.25%, ICC 66.3%; follow-up: CoV 4.78%, ICC 80.9%). Image quality of sT2w was rated inferior to T2 TSE (median for T2 TSE: 3; sT2w: 1-2), but measurements revealed good interrater reliability of sT2w (lesion counting: ICC 0.85; diameter measure: ICC 0.68 and 0.67).

Conclusion

GRAPPATINI is a feasible and robust T2 mapping sequence of the brain on intra- and intersubject level. The resulting sT2w depict brain lesions comparable to T2 TSE despite its inferior image quality.

Submillimeter T1 atlas for subject-specific abnormality detection at 7T

PURPOSE

Studies at 3T have shown that T₁ relaxometry enables characterization of brain tissues at the single-subject level by comparing individual physical properties to a normative atlas. In this work, an atlas of normative T₁ values at 7T is introduced with 0.6 mm isotropic resolution and its clinical potential is explored in comparison to 3T.

METHODS

T₁ maps were acquired in two separate healthy cohorts scanned at 3T and 7T. Using transfer learning, a template-based brain segmentation algorithm was adapted to ultra-high field imaging data. After segmenting brain tissues, volumes were normalized into a common space, and an atlas of normative T₁ values was established by modeling the T₁ inter-subject variability. A method for single-subject comparisons restricted to white matter and subcortical structures was developed by computing Z-scores. The comparison was applied to eight patients scanned at both field strengths for proof of concept.

RESULTS

The proposed method for morphometry delivered segmentation masks without statistically significant differences from those derived with the original pipeline at 3T and achieved accurate segmentation at 7T. The established normative atlas allowed characterizing tissue alterations in single-subject comparisons at 7T, and showed greater anatomical details compared with 3T results.

CONCLUSION

A high-resolution quantitative atlas with an adapted pipeline was introduced and validated. Several case studies on different clinical conditions showed the feasibility, potential and limitations of high-resolution single-subject comparisons based on quantitative MRI atlases. This method in conjunction with 7T higher resolution broadens the range of potential applications of quantitative MRI in clinical practice.

Assessment of idiopathic inflammatory myopathy using a deep learning method for muscle T2 mapping segmentation

OBJECTIVE

To investigate the utility of an automatic deep learning (DL) method for segmentation of T2 maps in patients with idiopathic inflammatory myopathy (IIM) against healthy controls, and also the association of quantitative T2 values in patients with laboratory and pulmonary findings.

METHODS

Structural MRI and T2 mapping of bilateral thigh muscles from patients with IIM and healthy volunteers were segmented using dedicated software based on a pre-trained convolutional neural network. Incremental and federated learning were implemented for continuous adaptation and improvement. Muscle T2 values derived from DL segmentation were compared between patients and healthy controls, and T2 values of patients were further analyzed with serum muscle enzymes, and interstitial lung disease (ILD) which was diagnosed and graded based on chest HRCT.

RESULTS

Overall, 64 patients (27 patients with dermatomyositis, 29 with polymyositis, and 8 with antisynthetase syndrome (ASS)) and 10 healthy controls were included. By using DL-based muscle segmentation, T2 values generated from T2 maps accurately differentiated patients from those of controls (p < 0.001) with a cutoff value of 36.4 ms (sensitivity 96.9%, and specificity 100%). In patients with IIM, muscle T2 values positively correlated with all the serum muscle enzymes (all p < 0.05). ILD score of patients with ASS was markedly higher than that of those without ASS (p = 0.011), while dissociation between the severity of muscular involvement and ILD was observed (p = 0.080).

CONCLUSION

Automatic DL could be used to segment thigh muscles and help quantitatively assess muscular inflammation of IIM through T2 mapping.

KEY POINTS

• Muscle T2 mapping automatically segmented by deep learning can differentiate IIM from healthy controls. • T2 value, an indicator of active muscle inflammation, positively correlates with serum muscle enzymes. • T2 mapping can detect muscle disease in patients with normal muscle enzyme levels.

Variability and reproducibility of multi-echo T2 relaxometry: Insights from multi-site, multi-session and multi-subject MRI acquisitions

Quantitative magnetic resonance imaging (qMRI) can increase the specificity and sensitivity of conventional weighted MRI to underlying pathology by comparing meaningful physical or chemical parameters, measured in physical units, with normative values acquired in a healthy population. This study focuses on multi-echo T₂ relaxometry, a qMRI technique that probes the complex tissue microstructure by differentiating compartment-specific T₂ relaxation times. However, estimation methods are still limited by their sensitivity to the underlying noise. Moreover, estimating the model's parameters is challenging because the resulting inverse problem is ill-posed, requiring advanced numerical regularization techniques. As a result, the estimates from distinct regularization strategies are different. In this work, we aimed to investigate the variability and reproducibility of different techniques for estimating the transverse relaxation time of the intra- and extra-cellular space (T2IE) in gray (GM) and white matter (WM) tissue in a clinical setting, using a multi-site, multi-session, and multi-run T₂ relaxometry dataset. To this end, we evaluated three different techniques for estimating the T₂ spectra (two regularized non-negative least squares methods and a machine learning approach). Two independent analyses were performed to study the effect of using raw and denoised data. For both the GM and WM regions, and the raw and denoised data, our results suggest that the principal source of variance is the inter-subject variability, showing a higher coefficient of variation (CoV) than those estimated for the inter-site, inter-session, and inter-run, respectively. For all reconstruction methods studied, the CoV ranged between 0.32 and 1.64%. Interestingly, the inter-session variability was close to the inter-scanner variability with no statistical differences, suggesting that T2IE is a robust parameter that could be employed in multi-site neuroimaging studies. Furthermore, the three tested methods showed consistent results and similar intra-class correlation (ICC), with values superior to 0.7 for most regions. Results from raw data were slightly more reproducible than those from denoised data. The regularized non-negative least squares method based on the L-curve technique produced the best results, with ICC values ranging from 0.72 to 0.92.

Data-driven myelin water imaging based on T1 and T2 relaxometry

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T₁ and T₂ relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T₁ /T₂ data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

Evaluating reproducibility and subject-specificity of microstructure-informed connectivity

Tractography enables identifying and evaluating the healthy and diseased brain's white matter pathways from diffusion-weighted magnetic resonance imaging data. As previous evaluation studies have reported significant false-positive estimation biases, recent microstructure-informed tractography algorithms have been introduced to improve the trade-off between specificity and sensitivity. However, a major limitation for characterizing the performance of these techniques is the lack of ground truth brain data. In this study, we compared the performance of two relevant microstructure-informed tractography methods, SIFT2 and COMMIT, by assessing the subject specificity and reproducibility of their derived white matter pathways. Specifically, twenty healthy young subjects were scanned at eight different time points at two different sites. Subject specificity and reproducibility were evaluated using the whole-brain connectomes and a subset of 29 white matter bundles. Our results indicate that although the raw tractograms are more vulnerable to the presence of false-positive connections, they are highly reproducible, suggesting that the estimation bias is subject-specific. This high reproducibility was preserved when microstructure-informed tractography algorithms were used to filter the raw tractograms. Moreover, the resulting track-density images depicted a more uniform coverage of streamlines throughout the white matter, suggesting that these techniques could increase the biological meaning of the estimated fascicles. Notably, we observed an increased subject specificity by employing connectivity pre-processing techniques to reduce the underlaying noise and the data dimensionality (using principal component analysis), highlighting the importance of these tools for future studies. Finally, no strong bias from the scanner site or time between measurements was found. The largest intraindividual variance originated from the sole repetition of data measurements (inter-run).

Periventricular gradient of T1 tissue alterations in multiple sclerosis

•

T₁ relaxation times alterations were assessed based on a study-specific atlas.

•

T₁ alterations depend on distance from lateral ventricles (“periventricular gradient”).

•

Gradient parameters correlate better with disability compared to conventional MRI.

Objective

Pathology in multiple sclerosis is not homogenously distributed. Recently, it has been shown that structures adjacent to CSF are more severely affected. A gradient of brain tissue involvement was shown with more severe pathology in periventricular areas and in proximity to brain surfaces such as the subarachnoid spaces and ependyma, and hence termed the “surface–in” gradient. Here, we study whether (i) the surface-in gradient of periventricular tissue alteration measured by T₁ relaxometry is already present in early multiple sclerosis patients, (ii) how it differs between early and progressive multiple sclerosis patients, and (iii) whether the gradient-derived metrics in normal-appearing white matter and lesions correlate better with physical disability than conventional MRI-based metrics.

Methods

Forty-seven patients with early multiple sclerosis, 52 with progressive multiple sclerosis, and 92 healthy controls were included in the study. Isotropic 3D T₁ relaxometry maps were obtained using the Magnetization-Prepared 2 Rapid Acquisition Gradient Echoes sequence at 3 T. After spatially normalizing the T₁ maps into a study-specific common space, T₁ inter-subject variability within the healthy cohort was modelled voxel-wise, yielding a normative T₁ atlas. Individual comparisons of each multiple sclerosis patient against the atlas were performed by computing z-scores. Equidistant bands of voxels were defined around the ventricles in the supratentorial white matter; the z-scores in these bands were analysed and compared between the early and progressive multiple sclerosis cohorts. Correlations between both conventional and z-score-gradient-derived MRI metrics and the Expanded Disability Status Scale were assessed.

Results

Patients with early and progressive multiple sclerosis demonstrated a periventricular gradient of T₁ relaxation time z-scores. In progressive multiple sclerosis, z-score-derived metrics reflecting the gradient of tissue abnormality in normal-appearing white matter were more strongly correlated with disability (maximal rho = 0.374) than the conventional lesion volume and count (maximal rho = 0.189 and 0.21 respectively). In early multiple sclerosis, the gradient of normal-appearing white matter volume with z-scores > 2 at baseline correlated with clinical disability assessed at two years follow-up.

Conclusion

Our results suggest that the surface-in white matter gradient of tissue alteration is detectable with T1 relaxometry and is already present at clinical disease onset. The periventricular gradients correlate with clinical disability. The periventricular gradient in normal-appearing white matter may thus qualify as a promising biomarker for monitoring of disease activity from an early stage in all phenotypes of multiple sclerosis.

Simultaneous 3D acquisition of 1 H MRF and 23 Na MRI

PURPOSE

To develop a 3D MR technique to simultaneously acquire proton multiparametric maps (T₁ , T₂ , and proton density) and sodium density weighted images over the whole brain.

METHODS

We implemented a 3D stack-of-stars MR pulse sequence which consists of interleaved proton (¹ H) and sodium (²³ Na) excitations, tailored slice encoding gradients that can encode the same slice for both nuclei, and simultaneous readout with different radial trajectories (¹ H, full-radial; ²³ Na, center-out radial). The receive chain of our 7T scanner was modified to enable simultaneous acquisition of ¹ H and ²³ Na signal. A heuristically optimized flip angle train was implemented for proton MR fingerprinting (MRF). The SNR and the accuracy of proton T₁ and T₂ were evaluated in phantoms. Finally, in vivo application of the method was demonstrated in five healthy subjects.

RESULTS

The SNR for the simultaneous measurement was almost identical to that for the single-nucleus measurements (<2% change). The proton T₁ and T₂ maps remained similar to the results from a reference 2D MRF technique (normalized RMS error in T₁ ≈ 4.2% and T₂ ≈ 11.3%). Measurements in healthy subjects corroborated these results and demonstrated the feasibility of our method for in vivo application. The in vivo T₁ values measured using our method were lower than the results measured by other conventional techniques.

CONCLUSIONS

With the 3D simultaneous implementation, we were able to acquire sodium and proton density weighted images in addition to proton T₁ , T₂ , and B1+ from ¹ H MRF that covers the whole brain volume within 21 min.

Value of T2 Mapping MRI for Prostate Cancer Detection and Classification

BACKGROUND

Currently, multi-parametric prostate MRI (mpMRI) consists of a qualitative T₂ , diffusion weighted, and dynamic contrast enhanced imaging. Quantification of T₂ imaging might further standardize PCa detection and support artificial intelligence solutions.

PURPOSE

To evaluate the value of T₂ mapping to detect prostate cancer (PCa) and to differentiate PCa aggressiveness.

STUDY TYPE

Retrospective single center cohort study.

POPULATION

Forty-four consecutive patients (mean age 67 years; median PSA 7.9 ng/mL) with mpMRI and verified PCa by subsequent targeted plus systematic MR/ultrasound (US)-fusion biopsy from February 2019 to December 2019.

FIELD STRENGTH/SEQUENCE

Standardized mpMRI at 3 T with an additionally acquired T₂ mapping sequence.

ASSESSMENT

Primary endpoint was the analysis of quantitative T₂ values and contrast differences/ratios (CD/CR) between PCa and benign tissue. Secondary objectives were the correlation between T₂ values, ISUP grade, apparent diffusion coefficient (ADC) value, and PI-RADS, and the evaluation of thresholds for differentiating PCa and clinically significant PCa (csPCa).

STATISTICAL TESTS

Mann-Whitney test, Spearman's rank (r_s ) correlation, receiver operating curves, Youden's index (J), and AUC were performed. Statistical significance was defined as P < 0.05.

RESULTS

Median quantitative T₂ values were significantly lower for PCa in PZ (85 msec) and PCa in TZ (75 msec) compared to benign PZ (141 msec) or TZ (97 msec) (P < 0.001). CD/CR between PCa and benign PZ (51.2/1.77), respectively TZ (19.8/1.29), differed significantly (P < 0.001). The best T₂ -mapping threshold for PCa/csPCa detection was for TZ 81/86 msec (J = 0.929/1.0), and for PZ 110 msec (J = 0.834/0.905). Quantitative T₂ values of PCa did not correlate significantly with the ISUP grade (r_s  = 0.186; P = 0.226), ADC value (r_s  = 0.138; P = 0.372), or PI-RADS (r_s  = 0.132; P = 0.392).

DATA CONCLUSION

Quantitative T₂ values could differentiate PCa in TZ and PZ and might support standardization of mpMRI of the prostate. Different thresholds seem to apply for PZ and TZ lesions. However, in the present study quantitative T₂ values were not able to indicate PCa aggressiveness.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Optimization of magnetization transfer contrast for EPI FLAIR brain imaging

PURPOSE

To evaluate the impact of magnetization transfer (MT) on brain tissue contrast in turbo-spin-echo (TSE) and EPI fluid-attenuated inversion recovery (FLAIR) images, and to optimize an MT-prepared EPI FLAIR pulse sequence to match the tissue contrast of a clinical reference TSE FLAIR protocol.

METHODS

Five healthy volunteers underwent 3T brain MRI, including single slice TSE FLAIR, multi-slice TSE FLAIR, EPI FLAIR without MT-preparation, and MT-prepared EPI FLAIR with variations of the MT-preparation parameters, including number of preparation pulses, pulse amplitude, and resonance offset. Automated co-registration and gray matter (GM) versus white matter (WM) segmentation was performed using a T1-MPRAGE acquisition, and the GM versus WM signal intensity ratio (contrast ratio) was calculated for each FLAIR acquisition.

RESULTS

Without MT preparation, EPI FLAIR showed poor tissue contrast (contrast ratio = 0.98), as did single slice TSE FLAIR. Multi-slice TSE FLAIR provided high tissue contrast (contrast ratio = 1.14). MT-prepared EPI FLAIR closely approximated the contrast of the multi-slice TSE FLAIR images for two combinations of the MT-preparation parameters (contrast ratio = 1.14). Optimized MT-prepared EPI FLAIR provided a 50% reduction in scan time compared to the reference TSE FLAIR acquisition.

CONCLUSION

Optimized MT-prepared EPI FLAIR provides comparable brain tissue contrast to the multi-slice TSE FLAIR images used in clinical practice.

Revisiting the T2 spectrum imaging inverse problem: Bayesian regularized non-negative least squares

Multi-echo T₂ magnetic resonance images contain information about the distribution of T₂ relaxation times of compartmentalized water, from which we can estimate relevant brain tissue properties such as the myelin water fraction (MWF). Regularized non-negative least squares (NNLS) is the tool of choice for estimating non-parametric T₂ spectra. However, the estimation is ill-conditioned, sensitive to noise, and highly affected by the employed regularization weight. The purpose of this study is threefold: first, we want to underline that the apparently innocuous use of two alternative parameterizations for solving the inverse problem, which we called the standard and alternative regularization forms, leads to different solutions; second, to assess the performance of both parameterizations; and third, to propose a new Bayesian regularized NNLS method (BayesReg). The performance of BayesReg was compared with that of two conventional approaches (L-curve and Chi-square (X²) fitting) using both regularization forms. We generated a large dataset of synthetic data, acquired in vivo human brain data in healthy participants for conducting a scan-rescan analysis, and correlated the myelin content derived from histology with the MWF estimated from ex vivo data. Results from synthetic data indicate that BayesReg provides accurate MWF estimates, comparable to those from L-curve and X², and with better overall stability across a wider signal-to-noise range. Notably, we obtained superior results by using the alternative regularization form. The correlations reported in this study are higher than those reported in previous studies employing the same ex vivo and histological data. In human brain data, the estimated maps from L-curve and BayesReg were more reproducible. However, the T₂ spectra produced by BayesReg were less affected by over-smoothing than those from L-curve. These findings suggest that BayesReg is a good alternative for estimating T₂ distributions and MWF maps.

Q-Dixon and GRAPPATINI T2 Mapping Parameters: A Whole Spinal Assessment of the Relationship Between Osteoporosis and Intervertebral Disc Degeneration

BACKGROUND

The relationship between osteoporosis and intervertebral disc (IVD) degeneration remains controversial. Novel quantitative Dixon (Q-Dixon) and GRAPPATINI T2 mapping techniques have shown potential for evaluating the biochemical components of the spine.

PURPOSE

To investigate the correlation of osteoporosis with IVD degeneration in postmenopausal women.

STUDY TYPE

Prospective.

SUBJECTS

A total of 105 postmenopausal females (mean age, 65 years; mean body mass index, 26 kg/m² ).

FIELD STRENGTH/SEQUENCE

3 T; sagittal; 6-echo Q-Dixon, multiecho spin-echo GRAPPATINI T2 mapping, turbo spin echo (TSE) T1-weighted and TSE T2-weighted sequences.

ASSESSMENT

The subjects were divided into normal (N = 47), osteopenia (N = 28), and osteoporosis (N = 30) groups according to quantitative computed tomography examination. The Pfirrmann grade of each IVD was obtained. Region of interest analysis was performed separately by two radiologists (X.L., with 10 years of experience, and S.C., with 20 years of experience) on a fat fraction map and T2 map to calculate the bone marrow fat fraction (BMFF) from the L1 to L5 vertebrae and the T2 values of each adjacent IVD separately.

STATISTICAL TESTS

One-way analysis of variance, post-hoc comparisons, and Kruskal-Wallis H tests were performed to evaluate the differences in the magnetic resonance imaging parameters between the groups. The relationships between BMFF and the IVD features were analyzed using the Spearman correlation analysis and linear regression models.

RESULTS

There were significant differences in BMFF among the three groups. The osteoporosis group had higher BMFF values (64.5 ± 5.9%). No significant correlation was found between BMFF and Pfirrmann grade (r = 0.251, P = 0.06). BMFF was significantly negatively correlated with the T2 of the adjacent IVD from L1 to L3 (r = -0.731; r = -0.637; r = -0.547), while significant weak correlations were found at the L4 to L5 levels (r = -0.337; r = -0.278).

DATA CONCLUSION

This study demonstrated that osteoporosis is associated with IVD degeneration.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 4.

Synthetic T2-weighted images of the lumbar spine derived from an accelerated T2 mapping sequence: Comparison to conventional T2w turbo spin echo

OBJECTIVES

To evaluate the diagnostic usefulness of synthetic T₂-weighted images of the lumbar spine derived from ten-fold undersampled k-space data using GRAPPATINI, a combination of a model-based approach for rapid T₂ and M₀ quantification (MARTINI) extended by generalized autocalibrating partial parallel acquistion (GRAPPA).

MATERIALS AND METHODS

Overall, 58 individuals (26 female, mean age 23.3 ± 8.1 years) were examined at 3 Tesla with sagittal and axial T₂w turbo spin echo (TSE) sequences compared to synthetic T_2-weighted contrasts derived at identical effective echo times and spatial resolutions. Two blinded readers graded disk degeneration and evaluated the lumbar intervertebral disks for present herniation or annular tear. One reader reassessed all studies after four weeks. Weighted kappa statistics were calculated to assess inter-rater and intra-rater agreement. Also, all studies were segmented manually by one reader to compute contrast ratios (CR) and contrast-to-noise ratios (CNR) of the nucleus pulposus and the annulus fibrosus.

RESULTS

Overall, the CR_T2w was 4.45 ± 1.80 and CR_T2synth was 4.71 ± 2.14. Both correlated (rsp = 0.768;p < 0.001) and differed (0.26 ± 1.38;p = 0.002) significantly. The CNR_T2w was 1.73 ± 0.52 and CNR_T2synth was 1.63 ± 0.50. Both correlated (rsp = 0.875;p < 0.001) and differed (-0.10 ± 0.25;p < 0.001) significantly. The inter-rater agreement was substantial to almost perfect (κ = 0.808-0.925) with the intra-rater agreement also substantial to almost perfect (κ = 0.862-0.963). The area under the curve of the receiver operating characteristics assessing disk herniation or annular tear ranged from 0.787 to 0.892.

CONCLUSIONS

This study concludes that synthetic images derived by GRAPPATINI can be used for clinical routine assessment with inter-rater and intra-rater agreements comparable to conventional T₂w TSE.

B subgroup detection in a small hospital transfusion service

The ABO blood group system includes phenotypes, or subgroups, that differ in the amount of A and B antigens present on the red blood cells (RBCs). These subgroups also differ in the A, B, or H substances present in secretions (for individuals who have the secretor phenotype). B subgroups are very rare and are less frequently reported than A subgroups. Usually, B subgroups are discovered during serologic testing when there is a discrepancy between RBC and serum grouping results. Subgroups of B are usually identified by a reference laboratory using molecular and adsorption-elution methods. This report details a case of a young, healthy, pregnant woman with a B subgroup detected by a small transfusion service using adsorption-elution methods. Serology and genotyping of the ABO gene was performed at a reference laboratory where the serology was consistent with a B subgroup, but no changes were identified in ABO gene sequencing. It is important to correctly identify B subgroups in donors and recipients to help resolve ABO discrepancies and potentially prevent ABO incompatibility in blood transfusion, thus minimizing transfusion reactions.

The ABO blood group system includes phenotypes, or subgroups, that differ in the amount of A and B antigens present on the red blood cells (RBCs). These subgroups also differ in the A, B, or H substances present in secretions (for individuals who have the secretor phenotype). B subgroups are very rare and are less frequently reported than A subgroups. Usually, B subgroups are discovered during serologic testing when there is a discrepancy between RBC and serum grouping results. Subgroups of B are usually identified by a reference laboratory using molecular and adsorption-elution methods. This report details a case of a young, healthy, pregnant woman with a B subgroup detected by a small transfusion service using adsorption-elution methods. Serology and genotyping of the ABO gene was performed at a reference laboratory where the serology was consistent with a B subgroup, but no changes were identified in ABO gene sequencing. It is important to correctly identify B subgroups in donors and recipients to help resolve ABO discrepancies and potentially prevent ABO incompatibility in blood transfusion, thus minimizing transfusion reactions.

Clinical implementation of accelerated T2 mapping: Quantitative magnetic resonance imaging as a biomarker for annular tear and lumbar disc herniation

OBJECTIVES

This study evaluates GRAPPATINI, an accelerated T2 mapping sequence combining undersampling and model-based reconstruction to facilitate the clinical implementation of T2 mapping of the lumbar intervertebral disc.

METHODS

Fifty-eight individuals (26 females, 32 males, age 23.3 ± 8.0 years) were prospectively examined at 3 T. This cohort study consisted of 19 patients, 20 rowers, and 19 volunteers. GRAPPATINI was conducted with the same parameters as a conventional 2D multi-echo spin-echo (MESE) sequence in 02:27 min instead of 13:18 min. Additional T2 maps were calculated after discarding the first echo (T2-WO1ST) and only using even echoes (T2-EVEN). Segmentation was done on the four most central slices. The resulting T2 values were compared for all four measurements.

RESULTS

T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of the nucleus pulposus of normal discs differed significantly from those of bulging discs or herniated discs (all p < 0.001). For the posterior annular region, only T2-GRAPPATINI showed a significant difference (p = 0.011) between normal and herniated discs. There was a significant difference between T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of discs with and without an annular tear for the nucleus pulposus (all p < 0.001). The nucleus pulposus' T2 at different degeneration states showed significant differences between all group comparisons of Pfirrmann grades for T2-GRAPPATINI (p = 0.000-0.018), T2-MESE (p = 0.000-0.015), T2-EVEN (p = 0.000-0.019), and T2-WO1ST (p = 0.000-0.015).

CONCLUSIONS

GRAPPATINI facilitates the use of T2 values as quantitative imaging biomarkers to detect disc pathologies such as degeneration, lumbar disc herniation, and annular tears while simultaneously shortening the acquisition time from 13:18 to 2:27 min.

KEY POINTS

• T2-GRAPPATINI, T2-MESE, T2-EVEN, and T2-WO1ST of the nucleus pulposus of normal discs differed significantly from those of discs with bulging or herniation (all p < 0.001). • The investigated T2 mapping techniques differed significantly in discs with and without annular tearing (all p < 0.001). • The nucleus pulposus' T2 showed significant differences between different stages of degeneration in all group comparisons for T2-GRAPPATINI (p = 0.000-0.018), T2-MESE (p = 0.000-0.015), T2-EVEN (p = 0.000-0.019), and T2-WO1ST (p = 0.000-0.015).

Differentiation between benign and malignant vertebral compression fractures using qualitative and quantitative analysis of a single fast spin echo T2-weighted Dixon sequence

Objectives

To determine and compare the qualitative and quantitative diagnostic performance of a single sagittal fast spin echo (FSE) T2-weighted Dixon sequence in differentiating benign and malignant vertebral compression fractures (VCF), using multiple readers and different quantitative methods.

Methods

From July 2014 to June 2020, 95 consecutive patients with spine MRI performed prior to cementoplasty for acute VCFs were retrospectively included. VCFs were categorized as benign (n = 63, mean age = 76 ± 12 years) or malignant (n = 32, mean age = 63 ± 12 years) with a best valuable comparator as a reference. Qualitative analysis was independently performed by four radiologists by categorizing each VCF as either benign or malignant using only the image sets provided by FSE T2-weighted Dixon sequences. Quantitative analysis was performed using two different regions of interest (ROI1-2) and three methods (signal drop, fat fraction (FF) from ROIs, FF maps). Diagnostic performance was compared using ROC curves analyses. Interobserver agreement was assessed using kappa statistics and intraclass correlation coefficients (ICC).

Results

The qualitative diagnostic performance ranged from area under the curve (AUC) = 0.97 (95% CI: 0.91–1.00) to AUC = 0.99 (95% CI: 0.95–1.0). The quantitative diagnostic performance ranged from AUC = 0.82 (95% CI: 0.73–0.89) to AUC = 0.97 (95% CI: 0.91–0.99). Pairwise comparisons showed no statistical difference in diagnostic performance (all p > 0.0013, Bonferroni-corrected p < 0.0011). All five cases with disagreement among the readers were correctly diagnosed at quantitative analysis using ROI2. Interobserver agreement was excellent for both qualitative and quantitative analyses.

Conclusions

A single FSE T2-weighted Dixon sequence can be used to differentiate benign and malignant VCF with high diagnostic performance using both qualitative and quantitative analyses, which can provide complementary information.

Key Points

• Qualitative analysis of a single FSE T2-weighted Dixon sequence yields high diagnostic performance and excellent observer agreement for differentiating benign and malignant compression fractures.

• The same FSE T2-weighted Dixon sequence allows quantitative assessment with high diagnostic performance.

• Quantitative data can readily be extracted from the FSE T2-weighted Dixon sequence and may provide complementary information to the qualitative analysis, which may be useful in doubtful cases.

Comparison of non-parametric T2 relaxometry methods for myelin water quantification

Multi-component T₂ relaxometry allows probing tissue microstructure by assessing compartment-specific T₂ relaxation times and water fractions, including the myelin water fraction. Non-negative least squares (NNLS) with zero-order Tikhonov regularization is the conventional method for estimating smooth T₂ distributions. Despite the improved estimation provided by this method compared to non-regularized NNLS, the solution is still sensitive to the underlying noise and the regularization weight. This is especially relevant for clinically achievable signal-to-noise ratios. In the literature of inverse problems, various well-established approaches to promote smooth solutions, including first-order and second-order Tikhonov regularization, and different criteria for estimating the regularization weight have been proposed, such as L-curve, Generalized Cross-Validation, and Chi-square residual fitting. However, quantitative comparisons between the available reconstruction methods for computing the T₂ distribution, and between different approaches for selecting the optimal regularization weight, are lacking. In this study, we implemented and evaluated ten reconstruction algorithms, resulting from the individual combinations of three penalty terms with three criteria to estimate the regularization weight, plus non-regularized NNLS. Their performance was evaluated both in simulated data and real brain MRI data acquired from healthy volunteers through a scan-rescan repeatability analysis. Our findings demonstrate the need for regularization. As a result of this work, we provide a list of recommendations for selecting the optimal reconstruction algorithms based on the acquired data. Moreover, the implemented methods were packaged in a freely distributed toolbox to promote reproducible research, and to facilitate further research and the use of this promising quantitative technique in clinical practice.

Compressed sensing with signal averaging for improved sensitivity and motion artifact reduction in fluorine-19 MRI

Fluorine-19 (19 F) MRI of injected perfluorocarbon emulsions (PFCs) allows for the non-invasive quantification of inflammation and cell tracking, but suffers from a low signal-to-noise ratio and extended scan time. To address this limitation, we tested the hypotheses that a 19 F MRI pulse sequence that combines a specific undersampling regime with signal averaging has both increased sensitivity and robustness against motion artifacts compared with a non-averaged fully sampled pulse sequence, when both datasets are reconstructed with compressed sensing. As a proof of principle, numerical simulations and phantom experiments were performed on selected variable ranges to characterize the point spread function of undersampling patterns, as well as the vulnerability to noise of undersampling and reconstruction parameters with paired numbers of x signal averages and acceleration factor x (NAx-AFx). The numerical simulations demonstrated that a probability density function that uses 25% of the samples to fully sample the k-space central area allowed for an optimal balance between limited blurring and artifact incoherence. At all investigated noise levels, the Dice similarity coefficient (DSC) strongly depended on the regularization parameters and acceleration factor. In phantoms, the motion robustness of an NA8-AF8 undersampling pattern versus NA1-AF1 was evaluated with simulated and real motion patterns. Differences were assessed with the DSC, which was consistently higher for the NA8-AF8 compared with the NA1-AF1 strategy, for both simulated and real cyclic motion patterns (P < 0.001). Both strategies were validated in vivo in mice (n = 2) injected with perfluoropolyether. Here, the images displayed a sharper delineation of the liver with the NA8-AF8 strategy than with the NA1-AF1 strategy. In conclusion, we validated the hypotheses that in 19 F MRI the combination of undersampling and averaging improves both the sensitivity and the robustness against motion artifacts.

Model-informed machine learning for multi-component T2 relaxometry

Recovering the T₂ distribution from multi-echo T₂ magnetic resonance (MR) signals is challenging but has high potential as it provides biomarkers characterizing the tissue micro-structure, such as the myelin water fraction (MWF). In this work, we propose to combine machine learning and aspects of parametric (fitting from the MRI signal using biophysical models) and non-parametric (model-free fitting of the T₂ distribution from the signal) approaches to T₂ relaxometry in brain tissue by using a multi-layer perceptron (MLP) for the distribution reconstruction. For training our network, we construct an extensive synthetic dataset derived from biophysical models in order to constrain the outputs with a priori knowledge of in vivo distributions. The proposed approach, called Model-Informed Machine Learning (MIML), takes as input the MR signal and directly outputs the associated T₂ distribution. We evaluate MIML in comparison to a Gaussian Mixture Fitting (parametric) and Regularized Non-Negative Least Squares algorithms (non-parametric) on synthetic data, an ex vivo scan, and high-resolution scans of healthy subjects and a subject with Multiple Sclerosis. In synthetic data, MIML provides more accurate and noise-robust distributions. In real data, MWF maps derived from MIML exhibit the greatest conformity to anatomical scans, have the highest correlation to a histological map of myelin volume, and the best unambiguous lesion visualization and localization, with superior contrast between lesions and normal appearing tissue. In whole-brain analysis, MIML is 22 to 4980 times faster than the non-parametric and parametric methods, respectively.

Comparison of 2D simultaneous multi-slice and 3D GRASE readout schemes for pseudo-continuous arterial spin labeling of cerebral perfusion at 3 T

OBJECTIVE

In this perfusion magnetic resonance imaging study, the performances of different pseudo-continuous arterial spin labeling (PCASL) sequences were compared: two-dimensional (2D) single-shot readout with simultaneous multislice (SMS), 2D single-shot echo-planar imaging (EPI) and multishot three-dimensional (3D) gradient and spin echo (GRASE) sequences combined with a background-suppression (BS) module.

MATERIALS AND METHODS

Whole-brain PCASL images were acquired from seven healthy volunteers. The performance of each protocol was evaluated by extracting regional cerebral blood flow (rCBF) measures using an inline morphometric segmentation prototype. Image data postprocessing and subsequent statistical analyses enabled comparisons at the regional and sub-regional levels.

RESULTS

The main findings were as follows: (i) Mean global CBF obtained across methods was were highly correlated, and these correlations were significantly higher among the same readout sequences. (ii) Temporal signal-to-noise ratio and gray-matter-to-white-matter CBF ratio were found to be equivalent for all 2D variants but lower than those of 3D-GRASE.

DISCUSSION

Our study demonstrates that the accelerated SMS readout can provide increased acquisition efficiency and/or a higher temporal resolution than conventional 2D and 3D readout sequences. Among all of the methods, 3D-GRASE showed the lowest variability in CBF measurements and thus highest robustness against noise.

Probing myelin content of the human brain with MRI: A review

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Normal volumetric and T1 relaxation time values at 1.5 T in segmented pediatric brain MRI using a MP2RAGE acquisition

OBJECTIVES

This study introduced a tailored MP2RAGE-based brain acquisition for a comprehensive assessment of the normal maturing brain.

METHODS

Seventy normal patients (35 girls and 35 boys) from 1 to 16 years of age were recruited within a prospective monocentric study conducted from a single University Hospital. Brain MRI examinations were performed at 1.5 T using a 20-channel head coil and an optimized 3D MP2RAGE sequence with a total acquisition time of 6:36 min. Automated 38 region segmentation was performed using the MorphoBox (template registration, bias field correction, brain extraction, and tissue classification) which underwent a major adaptation of three age-group T1-weighted templates. Volumetry and T1 relaxometry reference ranges were established using a logarithmic model and a modified Gompertz growth respectively.

RESULTS

Detailed automated brain segmentation and T1 mapping were successful in all patients. Using these data, an age-dependent model of normal brain maturation with respect to changes in volume and T1 relaxometry was established. After an initial rapid increase until 24 months of life, the total intracranial volume was found to converge towards 1400 mL during adolescence. The expected volumes of white matter (WM) and cortical gray matter (GM) showed a similar trend with age. After an initial major decrease, T1 relaxation times were observed to decrease progressively in all brain structures. The T1 drop in the first year of life was more pronounced in WM (from 1000-1100 to 650-700 ms) than in GM structures.

CONCLUSION

The 3D MP2RAGE sequence allowed to establish brain volume and T1 relaxation time normative ranges in pediatrics.

KEY POINTS

• The 3D MP2RAGE sequence provided a reliable quantitative assessment of brain volumes and T1 relaxation times during childhood. • An age-dependent model of normal brain maturation was established. • The normative ranges enable an objective comparison to a normal cohort, which can be useful to further understand, describe, and identify neurodevelopmental disorders in children.

Novel T2 Mapping for Evaluating Cervical Cancer Features by Providing Quantitative T2 Maps and Synthetic Morphologic Images: A Preliminary Study

BACKGROUND

The application value of T₂ mapping in evaluating cervical cancer (CC) features remains unclear.

PURPOSE

To investigate the role of T₂ values in evaluating CC classification, grade, and lymphovascular space invasion (LVSI) in comparison to apparent diffusion coefficient (ADC), and to compare synthetic T₂ -weighted (T₂ W) images calculated from T₂ values to conventional T₂ W images for CC staging.

STUDY TYPE

Retrospective.

POPULATION

Sixty-three patients with histopathologically confirmed CC.

FIELD STRENGTH/SEQUENCE

3T, conventional T₂ W turbo spin-echo, diffusion-weighted echo-planar, and accelerated T₂ mapping sequence.

ASSESSMENT

T₂ and ADC values between different pathological features of CC were compared. The diagnostic accuracies of conventional and synthetic T₂ W images in staging were also compared.

STATISTICAL TESTS

Parameters were compared using an independent t-test, Wilcoxon signed-rank test, and the chi-square test. Receiver operating characteristic analysis was performed.

RESULTS

The T₂ values varied significantly between well/moderately differentiated and poorly differentiated tumors ([92.8 ± 9.5 msec] vs. [83.8 ± 9.5 msec], P < 0.05) and between LVSI-positive and LVSI-negative CC ([82.2 ± 8.2 msec] vs. [93.9 ± 9.1 msec], P < 0.05). The ADC values showed a significant difference for grade ([0.76 ± 0.10 × 10^-3 mm² /s] vs. [0.65 ± 0.11 × 10^-3 mm² /s], P < 0.05) and no difference for LVSI status ([0.71 ± 0.11× 10^-3 mm² /s] vs. [0.73 ± 0.12× 10^-3 mm² /s], P = 0.472). There was no significant difference in T₂ and ADC values between squamous cell carcinoma and adenocarcinoma (P = 0.378 and P = 0.661, respectively). In MRI staging, the conventional and synthetic T₂ W images resulted in a similar accuracy (71% vs. 68%, P = 0.698).

DATA CONCLUSION

The accelerated T₂ mapping sequence may facilitate grading and staging of CC by providing quantitative T₂ maps and synthetic T₂ W images in one acquisition. T₂ values may be superior to ADC in predicting LVSI.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2 J. MAGN. RESON. IMAGING 2020;52:1859-1869.

Fast and high-resolution myelin water imaging: Accelerating multi-echo GRASE with CAIPIRINHA

PURPOSE

Although several MRI methods have been explored to achieve in vivo myelin quantification, imaging the whole brain in clinically acceptable times and sufficiently high resolution remains challenging. To address this problem, this work investigates the acceleration of multi-echo T₂ acquisitions based on the multi-echo gradient and spin echo (GRASE) sequence using CAIPIRINHA undersampling and adapted k-space reordering patterns.

METHODS

A prototype multi-echo GRASE sequence supporting CAIPIRINHA parallel imaging was implemented. Multi-echo T₂ data were acquired from 12 volunteers using the implemented sequence (1.6 × 1.6 × 1.6 mm³ , 84 slices, acquisition time [TA] = 10:30 min) and a multi-echo spin echo (MESE) sequence as reference (1.6 × 1.6 × 3.2 mm³ , single-slice, TA = 5:41 min). Myelin water fraction (MWF) maps derived from both acquisitions were compared via correlation and Bland-Altman analyses. In addition, scan-rescan datasets were acquired to evaluate the repeatability of the derived maps.

RESULTS

Resulting maps from the MESE and multi-echo GRASE sequences were found to be correlated (r = 0.83). The Bland-Altman analysis revealed a mean bias of -0.2% (P = .24) with the limits of agreement ranging from -3.7% to 3.3%. The Pearson's correlation coefficient among MWF values obtained from the scan-rescan datasets was found to be 0.95 and the mean bias equal to 0.11% (P = .32), indicating good repeatability of the retrieved maps.

CONCLUSION

By combining a 3D multi-echo GRASE sequence with CAIPIRINHA sampling, whole-brain MWF maps were obtained in 10:30 min with 1.6 mm isotropic resolution. The good correlation with conventional MESE-based maps demonstrates that the implemented sequence may be a promising alternative to time-consuming MESE acquisitions.

Magnetization transfer in magnetic resonance fingerprinting

PURPOSE

To study the effects of magnetization transfer (MT, in which a semi-solid spin pool interacts with the free pool), in the context of magnetic resonance fingerprinting (MRF).

METHODS

Simulations and phantom experiments were performed to study the impact of MT on the MRF signal and its potential influence on T1 and T2 estimation. Subsequently, an MRF sequence implementing off-resonance MT pulses and a dictionary with an MT dimension, generated by incorporating a two-pool model, were used to estimate the fractional pool size in addition to the B 1 + , T1 , and T2 values. The proposed method was evaluated in the human brain.

RESULTS

Simulations and phantom experiments showed that an MRF signal obtained from a cross-linked bovine serum sample is influenced by MT. Using a dictionary based on an MT model, a better match between simulations and acquired MR signals can be obtained (NRMSE 1.3% vs. 4.7%). Adding off-resonance MT pulses can improve the differentiation of MT from T1 and T2 . In vivo results showed that MT affects the MRF signals from white matter (fractional pool-size ~16%) and gray matter (fractional pool-size ~10%). Furthermore, longer T1 (~1060 ms vs. ~860 ms) and T2 values (~47 ms vs. ~35 ms) can be observed in white matter if MT is accounted for.

CONCLUSION

Our experiments demonstrated a potential influence of MT on the quantification of T1 and T2 with MRF. A model that encompasses MT effects can improve the accuracy of estimated relaxation parameters and allows quantification of the fractional pool size.

MRI T2 Mapping of the Knee Providing Synthetic Morphologic Images: Comparison to Conventional Turbo Spin-Echo MRI

Background Use of a T2 mapping sequence in addition to the conventional knee MRI protocol increases sensitivity to early cartilage lesions but is time consuming. Purpose To test the in vitro validity of quantitative data from an accelerated parallel T2 mapping sequence (combined generalized autocalibrating partially parallel acquisition and model-based accelerated relaxometry by iterative nonlinear inversion [GRAPPATINI]) of the knee and to compare in vivo synthetic images generated with this sequence with those generated with conventional morphologic sequences. Materials and Methods T2 estimations with GRAPPATINI were validated in vitro in comparison with T2 estimations with routine multisection multiecho and reference standard single-section single-echo spin-echo T2 mapping sequences by using a Bland-Altman plot. Synthetic morphologic images (intermediate-weighted sequence, 34-msec echo time; T2-weighted sequence, 80-msec echo time) were compared in vivo with corresponding conventional morphologic turbo spin-echo 3-T sequences by three readers in consecutive patients recruited retrospectively from February to May 2018. Synthetic and conventional morphologic images were compared by using rates of interreader agreement, κ statistics, and rates of findings. Results T2 values with GRAPPATINI were accurate compared with those obtained with the reference single-section single-echo sequence, with slight T2 overestimation (2.7 msec). Sixty-one patients (mean age, 43 years ± 16 [standard deviation]; 32 men) were included. The rate of agreement when one reader used synthetic morphologic images and the other used conventional sequences was not inferior to the rate of agreement when all readers used conventional sequences (upper bounds of 95% confidence intervals of differences between rates of agreement ≤ 4.8%). Interreader agreement was similar for the conventional set alone, the synthetic set alone, and when readers used different sets (two-by-two differences between κ values for all items ≤ 0.15). The rates of findings were not different between synthetic and conventional image sets (all P ≥ .07) except for two items (femoral trochlear cartilage [3.0% vs 0.3%, P = .006] and joint effusion [0.3% vs 2.7%, P = .005]). Conclusion This T2 mapping sequence yields, in one acquisition, accurate T2 values and synthetic morphologic images that are comparable with those obtained with conventional turbo spin-echo sequences. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Fritz in this issue.

Model-based super-resolution reconstruction of T2 maps

PURPOSE

High-resolution isotropic T₂ mapping of the human brain with multi-echo spin-echo (MESE) acquisitions is challenging. When using a 2D sequence, the resolution is limited by the slice thickness. If used as a 3D acquisition, specific absorption rate limits are easily exceeded due to the high power deposition of nonselective refocusing pulses. A method to reconstruct 1-mm³ isotropic T₂ maps is proposed based on multiple 2D MESE acquisitions. Data were undersampled (10-fold) to compensate for the prolonged scan time stemming from the super-resolution acquisition.

THEORY AND METHODS

The proposed method integrates a classical super-resolution with an iterative model-based approach to reconstruct quantitative maps from a set of undersampled low-resolution data. The method was tested on numerical and multipurpose phantoms, and in vivo data. T₂ values were assessed with a region-of-interest analysis using a single-slice spin-echo and a fully sampled MESE acquisition in a phantom, and a MESE acquisition in healthy volunteers.

RESULTS

Numerical simulations showed that the best trade-off between acceleration and number of low-resolution datasets is 10-fold acceleration with 4 acquisitions (acquisition time = 18 min). The proposed approach showed improved resolution over low-resolution images for both phantom and brain. Region-of-interest analysis of the phantom compartments revealed that at shorter T₂ , the proposed method was comparable with the fully sampled MESE. For the volunteer data, the T₂ values found in the brain structures were consistent across subjects (8.5-13.1 ms standard deviation).

CONCLUSION

The proposed method addresses the inherent limitations associated with high-resolution T₂ mapping and enables the reconstruction of 1 mm³ isotropic relaxation maps with a 10 times faster acquisition.

Quantitative brain relaxation atlases for personalized detection and characterization of brain pathology

PURPOSE

To exploit the improved comparability and hardware independency of quantitative MRI, databases of MR physical parameters in healthy tissue are required, to which tissue properties of patients can be compared. In this work, normative values for longitudinal and transverse relaxation times in the brain were established and tested in single-subject comparisons for detection of abnormal relaxation times.

METHODS

Relaxometry maps of the brain were acquired from 52 healthy volunteers. After spatially normalizing the volumes into a common space, T₁ and T₂ inter-subject variability within the healthy cohort was modeled voxel-wise. A method for a single-subject comparison against the atlases was developed by computing z-scores with respect to the established healthy norms. The comparison was applied to two multiple sclerosis and one clinically isolated syndrome cases for a proof of concept.

RESULTS

The established atlases exhibit a low variation in white matter structures (median RMSE of models equal to 32 ms for T₁ and 4 ms for T₂ ), indicating that relaxation times are in a narrow range for normal tissues. The proposed method for single-subject comparison detected relaxation time deviations from healthy norms in the example patient data sets. Relaxation times were found to be increased in brain lesions (mean z-scores >5). Moreover, subtle and confluent differences (z-scores ~2-4) were observed in clinically plausible regions (between lesions, corpus callosum).

CONCLUSIONS

Brain T₁ and T₂ quantitative norms were derived voxel-wise with low variability in healthy tissue. Example patient deviation maps demonstrated good sensitivity of the atlases for detecting relaxation time alterations.

Quantitative T2 mapping accelerated by GRAPPATINI for evaluation of muscles in patients with myositis

OBJECTIVE

Dermatomyositis (DM) and polymyositis (PM) make up the largest group of potentially treatable myopathies and require early diagnosis. This study investigates whether the edema of thigh muscles in DM/PM can be quantitatively assessed by a novel accelerated T2 mapping technique-GRAPPATINI.

METHODS

Three conventional MR sequences and GRAPPATINI accelerated T2 mapping of bilateral thighs from 20 patients (7 DM and 13 PM) and 10 healthy volunteers were prospectively carried out on a 3 T MR scanner. Afterwards, T2 values of 477 thigh muscles from the patients and the healthy controls were manually measured. In addition, the correlations between T2 values and serum muscle enzymes in patients were also analyzed.

RESULTS

The new GRAPPATINI technique made quantitative T2 mapping of bilateral thighs feasible with a scanning time of only 2 min 18 s. Moreover, GRAPPATINI-generated T2 values of muscles from patients were markedly higher than those from healthy subjects (p < 0.001). GRAPPATINI accelerated T2 mapping appeared a more sensitive technique in that some DM/PM muscles appearing normal per conventional MRI had increased T2 relaxation time. Furthermore, GRAPPATINI-generated T2 values of DM/PM thigh muscles positively correlated with serum enzyme levels (p < 0.001), which reflected the severity of myopathy.

CONCLUSION

GRAPPATINI can significantly shorten acquisition time of T2 mapping and may potentially be applied clinically in DM and PM.

ADVANCES IN KNOWLEDGE

GRAPPATINI acceleration makes T2 mapping feasible in clinical practice in providing quantitative information regarding thigh muscle inflammation in DM and PM.

Patient respiratory-triggered quantitative T2 mapping in the pancreas

Background

Long acquisition times and motion sensitivity limit T₂ mapping in the abdomen. Accelerated mapping at 3 T may allow for quantitative assessment of diffuse pancreatic disease in patients during free‐breathing.

Purpose

To test the feasibility of respiratory‐triggered quantitative T₂ analysis in the pancreas and correlate T₂‐values with age, body mass index, pancreatic location, main pancreatic duct dilatation, and underlying pathology.

Study Type

Retrospective single‐center pilot study.

Population

Eighty‐eight adults.

Field Strength/Sequence

Ten‐fold accelerated multiecho‐spin‐echo 3 T MRI sequence to quantify T₂ at 3 T.

Assessment

Two radiologists independently delineated three regions of interest inside the pancreatic head, body, and tail for each acquisition. Means and standard deviations for T₂ values in these regions were determined. T₂‐value variation with demographic data, intraparenchymal location, pancreatic duct dilation, and underlying pancreatic disease was assessed.

Statistical Tests

Interreader reliability was determined by calculating the interclass coefficient (ICCs). T₂ values were compared for different pancreatic locations by analysis of variance (ANOVA). Interpatient associations between T₂ values and demographical, clinical, and radiological data were calculated (ANOVA).

Results

The accelerated T₂ mapping sequence was successfully performed in all participants (mean acquisition time, 2:48 ± 0:43 min). Low T₂ value variability was observed across all patients (intersubject) (head: 60.2 ± 8.3 msec, body: 63.9 ± 11.5 msec, tail: 66.8 ± 16.4 msec). Interreader agreement was good (ICC, 0.82, 95% confidence interval: 0.77–0.86). T₂‐values differed significantly depending on age (P < 0.001), location (P < 0.001), main pancreatic duct dilatation (P < 0.001), and diffuse pancreatic disease (P < 0.03).

Data Conclusion

The feasibility of accelerated T₂ mapping at 3 T in moving abdominal organs was demonstrated in the pancreas, since T₂ values were stable and reproducible. In the pancreatic parenchyma, T₂‐values were significantly dependent on demographic and clinical parameters.

Level of Evidence: 4

Technical Efficacy: Stage 1

J. Magn. Reson. Imaging 2019;50:410–416.

Personalized pathology maps to quantify diffuse and focal brain damage

Background and objectives

Quantitative MRI (qMRI) permits the quantification of brain changes compatible with inflammation, degeneration and repair in multiple sclerosis (MS) patients. In this study, we propose a new method to provide personalized maps of tissue alterations and longitudinal brain changes based on different qMRI metrics, which provide complementary information about brain pathology.

Methods

We performed baseline and two-years follow-up on (i) 13 relapsing-remitting MS patients and (ii) four healthy controls. A group consisting of up to 65 healthy controls was used to compute the reference distribution of qMRI metrics in healthy tissue. All subjects underwent 3T MRI examinations including T1, T2, T2* relaxation and Magnetization Transfer Ratio (MTR) imaging. We used a recent partial volume estimation algorithm to estimate the concentration of different brain tissue types on T1 maps; then, we computed a deviation map (z-score map) for each contrast at both time-points. Finally, we subtracted those deviation maps only for voxels showing a significant difference with healthy tissue in one of the time points, to obtain a difference map for each subject.

Results and conclusion

Control subjects did not show any significant z-score deviations or longitudinal z-score changes. On the other hand, MS patients showed brain regions with cross-sectional and longitudinal concomitant increase in T1, T2, T2* z-scores and decrease of MTR z-scores, suggesting brain tissue degeneration/loss. In the lesion periphery, we observed areas with cross-sectional and longitudinal decreased T1/T2 and slight decrease in T2* most likely related to iron accumulation. Moreover, we measured longitudinal decrease in T1, T2 - and to a lesser extent in T2* - as well as a concomitant increase in MTR, suggesting remyelination/repair.

In summary, we have developed a method that provides whole-brain personalized maps of cross-sectional and longitudinal changes in MS patients, which are computed in patient space. These maps may open new perspectives to complement and support radiological evaluation of brain damage for a given patient.

Accelerated T2 mapping combining parallel MRI and model-based reconstruction: GRAPPATINI

BACKGROUND

Quantitative T₂ measurements are sensitive to intra- and extracellular water accumulation and myelin loss. Therefore, quantitative T₂ promises to be a good biomarker of disease. However, T₂ measurements require long acquisition times.

PURPOSE

To accelerate T₂ quantification and subsequent generation of synthetic T₂ -weighted (T₂ -w) image contrast for clinical research and routine. To that end, a recently developed model-based approach for rapid T₂ and M₀ quantification (MARTINI) based on undersampling k-space, was extended by parallel imaging (GRAPPA) to enable high-resolution T₂ mapping with access to T₂ -w images in less than 2 minutes acquisition time for the entire brain.

STUDY TYPE

Prospective cross-sectional study.

SUBJECTS/PHANTOM

Fourteen healthy subjects and a multipurpose phantom.

FIELD STRENGTH/SEQUENCE

Carr-Purcell-Meiboom-Gill sequence at a 3T scanner.

ASSESSMENT

The accuracy and reproducibility of the accelerated T₂ quantification was assessed. Validations comprised MRI studies on a phantom as well as the brain, knee, prostate, and liver from healthy volunteers. Synthetic T₂ -w images were generated from computed T₂ and M₀ maps and compared to conventional fast spin-echo (SE) images.

STATISTICAL TESTS

Root mean square distance (RMSD) to the reference method and region of interest analysis.

RESULTS

The combination of MARTINI and GRAPPA (GRAPPATINI) lead to a 10-fold accelerated T₂ mapping protocol with 1:44 minutes acquisition time and full brain coverage. The RMSD of GRAPPATINI increases less (4.3%) than a 10-fold MARTINI reconstruction (37.6%) in comparison to the reference. Reproducibility tests showed low standard deviation (SD) of T₂ values in regions of interest between scan and rescan (<0.4 msec) and across subjects (<4 msec).

DATA CONCLUSION

GRAPPATINI provides highly reproducible and fast whole-brain T₂ maps and arbitrary synthetic T₂ -w images in clinically compatible acquisition times of less than 2 minutes. These abilities are expected to support more widespread clinical applications of quantitative T₂ mapping.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2018;48:359-368.

An in vivo study of the orientation-dependent and independent components of transverse relaxation rates in white matter

Diffusion-weighted imaging (DWI) provides information that allows the estimation of white-matter (WM) fibre orientation and distribution, but it does not provide information about myelin density, fibre concentration or fibre size within each voxel. On the other hand, quantitative relaxation contrasts (like the apparent transverse relaxation, R2∗) offer iron and myelin-related contrast, but their dependence on the orientation of microstructure with respect to the applied magnetic field, B₀ , is often neglected. The aim of this work was to combine the fibre orientation information retrieved from the DWI acquisition and the sensitivity to microstructural information from quantitative relaxation parameters. The in vivo measured quantitative transverse relaxation maps (R₂ and R2∗) were decomposed into their orientation-dependent and independent components, using the DWI fibre orientation information as prior knowledge. The analysis focused on major WM fibre bundles such as the forceps major (FMj), forceps minor (FMn), cingulum (CG) and corticospinal tracts (CST). The orientation-dependent R₂ parameters, despite their small size (0-1.5 Hz), showed higher variability across different fibre populations, while those derived from R2∗, although larger (3.1-4.5 Hz), were mostly bundle-independent. With this article, we have, for the first time, attempted the in vivo characterization of the orientation-(in)dependent components of the transverse relaxation rates and demonstrated that the orientation of WM fibres influences both R₂ and R2∗ contrasts.