T1 Estimation Research Articles

Motion correction for native myocardial T1 mapping using self-supervised deep learning registration with contrast separation.

In myocardial T1 mapping, undesirable motion poses significant challenges because uncorrected motion can affect T1 estimation accuracy and cause incorrect diagnosis. In this study, we propose and evaluate a motion correction method for myocardial T1 mapping using self-supervised deep learning based registration with contrast separation (SDRAP). A sparse coding based method was first proposed to separate the contrast component from T1 -weighted (T1w) images. Then, a self-supervised deep neural network with cross-correlation (SDRAP-CC) or mutual information as the registration similarity measurement was developed to register contrast separated images, after which signal fitting was performed on the motion corrected T1w images to generate motion corrected T1 maps. The registration network was trained and tested in 80 healthy volunteers with images acquired using the modified Look-Locker inversion recovery (MOLLI) sequence. The proposed SDRAP was compared with the free form deformation (FFD) registration method regarding (1) Dice similarity coefficient (DSC) and mean boundary error (MBE) of myocardium contours, (2) T1 value and standard deviation (SD) of T1 fitting, (3) subjective evaluation score for overall image quality and motion artifact level, and (4) computation time. Results showed that SDRAP-CC achieved the highest DSC of 85.0 ± 3.9% and the lowest MBE of 0.92 ± 0.25 mm among the methods compared. Additionally, SDRAP-CC performed the best by resulting in lower SD value (28.1 ± 17.6 ms) and higher subjective image quality scores (3.30 ± 0.79 for overall quality and 3.53 ± 0.68 for motion artifact) evaluated by a cardiologist. The proposed SDRAP took only 0.52 s to register one slice of MOLLI images, achieving about sevenfold acceleration over FFD (3.7 s/slice).

Application of a novel T1 retrospective quantification using internal references (T1-REQUIRE) algorithm to derive quantitative T1 relaxation maps of the brain.

Most MRI sequences used clinically are qualitative or weighted. While such images provide useful information for clinicians to diagnose and monitor disease progression, they lack the ability to quantify tissue damage for more objective assessment. In this study, an algorithm referred to as the T1-REQUIRE is presented as a proof-of-concept which uses nonlinear transformations to retrospectively estimate T1 relaxation times in the brain using T1-weighted MRIs, the appropriate signal equation, and internal, healthy tissues as references. T1-REQUIRE was applied to two T1-weighted MR sequences, a spin-echo and a MPRAGE, and validated with a reference standard T1 mapping algorithm in vivo. In addition, a multiscanner study was run using MPRAGE images to determine the effectiveness of T1-REQUIRE in conforming the data from different scanners into a more uniform way of analyzing T1-relaxation maps. The T1-REQUIRE algorithm shows good agreement with the reference standard (Lin's concordance correlation coefficients of 0.884 for the spin-echo and 0.838 for the MPRAGE) and with each other (Lin's concordance correlation coefficient of 0.887). The interscanner studies showed improved alignment of cumulative distribution functions after T1-REQUIRE was performed. T1-REQUIRE was validated with a reference standard and shown to be an effective estimate of T1 over a clinically relevant range of T1 values. In addition, T1-REQUIRE showed excellent data conformity across different scanners, providing evidence that T1-REQUIRE could be a useful addition to big data pipelines.

Bowel wall MRI T1 relaxation estimates for assessment of intestinal inflammation in pediatric Crohn's disease.

To compare bowel wall T1 relaxation estimates in young patients with newly diagnosed ileal CD to healthy control participants, characterize their change over time in response to biologic medical therapy, and evaluate their associations with clinical markers of intestinal inflammation. Patients with newly diagnosed ileal CD and healthy control participants were prospectively recruited between December 2018 and October 2021. Patients underwent research MRI examinations of the bowel at baseline and at 6-weeks and 6-months into biologic medical treatment; control participants underwent single MRI examinations. MRI examinations included native T1 relaxometry of the terminal ileum using a modified Look-Locker inversion recovery (MOLLI) sequence. T1 estimates were measured on scanner-generated parametric maps. Clinical markers of intestinal inflammation were recorded at each visit. Group differences were assessed using the Mann-Whitney U test; the Friedman test was used to assess longitudinal changes in T1 estimates. Spearman correlation was used to evaluate associations between T1 estimates and inflammatory markers. Nineteen participants with CD (12 males; median age 14years) and 15 control participants (7 males; median age 17years) were included in the study. Bowel wall T1 estimates in CD patients (median 1302ms) were significantly longer compared to control participants (median 1159ms) (p < 0.001). In CD patients, T1 estimates changed over time after treatment (p = 0.001), with largest reductions between baseline and 6-weeks (p < 0.001). T1 estimates correlated with inflammatory markers, including erythrocyte sedimentation rate (ρ = 0.35; p = 0.01), c-reactive protein level (ρ = 0.34; p = 0.02), and weighted Pediatric Crohn's Disease Activity Index (ρ = 0.39; p = 0.005). T1 estimates did not correlate with serum albumin (ρ = - 0.28; p = 0.051) and fecal calprotectin (ρ = 0.07; p = 0.63). Bowel wall T1 estimates are abnormally increased in newly diagnosed ileal CD patients and decrease in response to medical therapy.

Cancellation of streak artifacts in radial abdominal imaging using interference null space projection.

In radial abdominal imaging, it has been commonly observed that signal from the arms cause streaks due to system imperfections. We previously introduced a streak removal technique (B-STAR), which is inherently spatially variant and limited to work in image space. In this work, we propose a spatially invariant streak cancellation technique (CACTUS), which can be applied in either image space or k-space and is compatible with iterative reconstructions. Streak sources are typically spatially localized and can be represented using a low-dimensional subspace. CACTUS identifies the streak subspace by leveraging the spatial redundancy of receiver coils and projects the data onto the streak null space to eliminate the streaks. When applied in k-space, CACTUS can be combined with iterative reconstructions. CACTUS was tested in phantoms and in vivo abdominal imaging using a radial turbo spin-echo pulse sequence. In phantoms, CACTUS improved T2 estimation in comparison to previous de-streaking methods. In vivo experiments showed that CACTUS reduced streaks and yielded T2 estimation, in regions affected by streaks, closer to a streak-free reference. Evaluation using a clinical abdominal dataset (n=20) showed that CACTUS is comparable to B-STAR and yields significantly better signal preservation and streak cancellation than coil removal and suppression methods. CACTUS provides superior signal preservation and streak reduction performance compared to coil removal and suppression methods. As a clear advantage over B-STAR, CACTUS can be integrated with iterative reconstruction methods. In abdominal T2 mapping, CACTUS improves the accuracy of parameter estimation in areas affected by streaks.

SAturation-recovery and Variable-flip-Angle-based three-dimensional free-breathing cardiovascular magnetic resonance T1 mapping at 3T.

The purpose of the current study was to develop and validate a three-dimensional (3D) free-breathing cardiac T1 -mapping sequence using SAturation-recovery and Variable-flip-Angle (SAVA). SAVA sequentially acquires multiple electrocardiogram-triggered volumes using a multishot spoiled gradient-echo sequence. The first volume samples the equilibrium signal of the longitudinal magnetization, where a flip angle of 2° is used to reduce the time for the magnetization to return to equilibrium. The succeeding three volumes are saturation prepared with variable delays, and are acquired using a 15° flip angle to maintain the signal-to-noise ratio. A diaphragmatic navigator is used to compensate the respiratory motion. T1 is calculated using a saturation-recovery model that accounts for the flip angle. We validated SAVA by simulations, phantom, and human subject experiments at 3T. SAVA was compared with modified Look-Locker inversion recovery (MOLLI) and saturation-recovery single-shot acquisition (SASHA) in vivo. In phantoms, T1 by SAVA had good agreement with the reference (R2 = 0.99). In vivo 3D T1 mapping by SAVA could achieve an imaging resolution of 1.25 × 1.25 × 8mm3 . Both global and septal T1 values by SAVA (1347 ± 37 and 1332 ± 42 ms) were in between those by SASHA (1612 ± 63 and 1618 ± 51 ms) and MOLLI (1143 ± 59 and 1188 ± 65 ms). According to the standard deviation (SD) and coefficient of variation (CV), T1 precision measured by SAVA (SD: 99 ± 14 and 60 ± 8ms; CV: 7.4% ± 0.9% and 4.5% ± 0.6%) was comparable with MOLLI (SD: 99 ± 25 and 46 ± 12 ms; CV: 8.8% ± 2.5% and 3.9% ± 1.1%) and superior to SASHA (SD: 222 ± 89 and 132 ± 33 ms; CV: 13.8% ± 5.5% and 8.1% ± 2.0%). It was concluded that the proposed free-breathing SAVA sequence enables more efficient 3D whole-heart T1 estimation with good accuracy and precision.

Saturated multi‐delay renal arterial spin labeling technique for simultaneous perfusion and T1 quantification in kidneys

To propose a free-breathing simultaneous multi-delay arterial spin labeling (ASL) and T1 mapping technique with a stepwise kinetic model for renal assessment in a single 4-min scan at 3T. The proposed saturated multi-delay renal arterial spin labeling (SAMURAI) sequence used flow-sensitive alternating inversion recovery (FAIR) preparation, followed by acquisition of 9 images with Look-Locker spoiled gradient recalled echo (SPGR). Pre-saturation at the imaging slice was used to achieve saturation-based T1 mapping. A 4-step 2-compartment kinetic model was proposed to characterize water transition through artery- and tissue-compartment. The impact of the Look-Locker sampling scheme on the ASL signal was corrected in this model. T1 estimation with dictionary searching method and perfusion quantification based on the proposed kinetic model fitting were conducted after groupwise registration of the acquired images. The feasibility and repeatability of SAMURAI were validated in healthy subjects (n=11) and patients with different renal diseases (n=4). The proposed SAMURAI technique can provide accurate T1 map with strong correlation (R2 =0.98) with inversion recovery spin echo (IR-SE) on phantom. SAMURAI provided equally reliable whole kidney and cortical ASL and T1 quantification results compared with multi-TI FAIR (intraclass correlation coefficient [ICC], 0.880-0.958) and IR-SPGR (ICC, 0.875-0.912), respectively. Low renal blood flow and increased T1 were detected by SAMURAI in the affected kidneys of the patients. SAMURAI had excellent scan-rescan repeatability (ICC, 0.905-0.992) and significantly reduced scan time (4 min 6 s vs. 45 min for 9 TIs) compared to multi-TI FAIR. The proposed SAMURAI technique is feasible and repeatable for simultaneously quantifying T1 and perfusion of kidneys with high time-efficiency.

Using Variable Flip Angle (VFA) and Modified Look-Locker Inversion Recovery (MOLLI) T1 mapping in clinical OE-MRI

Background and purposeThe imaging technique known as Oxygen-Enhanced MRI is under development as a noninvasive technique for imaging hypoxia in tumours and pulmonary diseases. While promising results have been shown in preclinical experiments, clinical studies have mentioned experiencing difficulties with patient motion, image registration, and the limitations of single-slice images compared to 3D volumes. As clinical studies begin to assess feasibility of using OE-MRI in patients, it is important for researchers to communicate about the practical challenges experienced when using OE-MRI on patients to help the technique advance. Materials and methodsWe report on our experience with using two types of T1 mapping (MOLLI and VFA) for a recently completed OE-MRI clinical study on oropharyngeal squamous cell carcinoma. ResultsWe report: (1) the artefacts and practical difficulties encountered in this study; (2) the difference in estimated T1 from each method used - the VFA T1 estimation was higher than the MOLLI estimation by 27% on average; (3) the standard deviation within the tumour ROIs - there was no significant difference in the standard deviation seen within the tumour ROIs from the VFA versus MOLLI; and (4) the OE-MRI response collected from either method. Lastly, we collated the MRI acquisition details from over 45 relevant manuscripts as a convenient reference for researchers planning future studies. ConclusionWe have reported our practical experience from an OE-MRI clinical study, with the aim that sharing this is helpful to researchers planning future studies. In this study, VFA was a more useful technique for using OE-MRI in tumours than MOLLI T1 mapping.

Effects of T2 * on accuracy and precision of dynamic T1 measurements using the single reference variable flip angle method: a simulation study.

To study in simulation and in theory the accuracy and precision of dynamic T1 measurements obtained using the previously published single-reference variable flip angle (SR-VFA) technique, with a focus on the effects of dynamic changes in T2 * on the calculation. Monte Carlo simulations were performed over 1000 noisy iterations for the VFA method, the SR-VFA method, and a proposed method, SR-VFA with a T2 * correction (SR-VFA-T2 *). Dynamic T1 estimates were calculated analytically for each method, with signals modeled by the steady-state spoiled gradient echo equation. The mean and standard deviation of these estimates were calculated and compared to truth, while varying repetition time (TR), baseline and dynamic T1 , echo time (TE), baseline and dynamic T2 *, flip angles, and the number of averages on baseline scans. Additionally, the variance of T1 in the SR-VFA and SR-VFA-T2 * methods was derived analytically based on the theory of propagation of errors. This equation was used to produce an inverse-variance weighted linear combination to improve T1 mapping precision in the SR-VFA-T2 * method. Flip angle sensitivity of dynamic T1 precision in the SR-VFA and SR-VFA-T2 * methods was also performed. Substantial bias can be produced by the SR-VFA method when the ratio of the T2 * decay of the dynamic signal versus that of the baseline signals deviates from 1, with a 0.01 deviation leading to approximately a 1% bias in cases of high SNR and TR ≫ T1 . This bias can be corrected by estimating the baseline and dynamic T2 * values in this ratio via multiecho measurements. The bias and precision of the SR-VFA-T2 * method, when normalized to scan time, is found to rival and sometimes improve upon the two flip angle VFA method when an inverse variance weighted linear combination is applied across its multiecho T1 maps. The analytic variance equation presented is found to be accurate within 1% relative to the Monte Carlo simulations over a broad parameter space. Flip angle ranges that maximize SR-VFA and SR-VFA-T2 *T1 precision over a broad parameter space are given, and each is defined relative to TR and T1 . Multiecho SR-VFA-T2 * T1 mapping is found in simulation and theory to be a promising alternative to the VFA method that maintains speed of the SR-VFA method with accuracy and precision similar to the VFA method.

Tailored magnetic resonance fingerprinting for simultaneous non-synthetic and quantitative imaging: A repeatability study.

The goals of this study include: (a) generating tailored magnetic resonance fingerprinting (TMRF) based non-synthetic imaging; (b) assessing the repeatability of TMRF and deep learning-based mapping of in vitro ISMRM/NIST phantom and in vivo brain data of healthy human subjects. We have acquired qualitative images obtained from the vendor-supplied gold standard (GS), MRF (synthetic), and TMRF (non-synthetic) on one representative healthy human brain. We also acquired 30 datasets on the ISMRM/NIST phantom for the in vitro repeatability study on a GE Discovery 3T MR750w scanner using the TMRF sequence. We compared T1 and T2 maps generated from 30 ISMRM/NIST phantom datasets to the spin-echo (SE) based GS method as part of the in vitro repeatability study. R-squared coefficient of determination in a simple linear regression and Bland-Altman analysis were computed for 30 datasets of ISMRM/NIST phantom to assess the accuracy of in vitro quantitative TMRF data. The repeatability of T1 and T2 estimates by TMRF was evaluated by calculating the standard deviation (SD) divided by the average of 30 datasets for each sphere, respectively. We acquired 10 volunteers for the in vivo repeatability study on the same scanner using the same TMRF sequence. These volunteers were imaged five times with two runs per repetition, resulting in 100 in vivo datasets. Five contrasts, T1 and T2 maps of 10 human volunteers acquired over five repetitions, were evaluated in the in vivo repeatability study. We computed the intraclass correlation coefficient (ICC) of the signal-to-noise ratio (SNR), signal intensities, T1 and T2 relaxation times in white matter (WM), and gray matter (GM). The synthetic images generated from MRF show partial volume and flow artifacts compared to non-synthetic images obtained from TMRF images and the GS. In vitro studies show that TMRF estimates have less than 5% variations except sphere 14 in the T2 array (6.36%). TMRF and SE relaxometry measurements were strongly correlated; R2 values were 0.9958 and 0.9789 for T1 and T2 estimates, respectively. Based on the ICC values, SNR, mean intensity values, and relaxation times of WM and GM for the in vivo studies were consistent. T1 and T2 values of WM and GM were similar to previously published values. The mean ± SD of T1 and T2 for WM for ten subjects and five repeats are 992 ± 41ms and 99 ± 6ms, while the corresponding values for T1 and T2 for GM are 1598 ± 73ms and 152 ± 14ms. TMRF and deep learning-based reconstruction produce repeatable, non-synthetic multi-contrast images, and parametric maps simultaneously.

Efficacy of D-cycloserine augmented brief intensive cognitive-behavioural therapy for paediatric obsessive-compulsive disorder: A randomised clinical trial.

ObjectiveTo examine the efficacy of weight‐adjusted D‐cycloserine (DCS) (35 or 70 mg) relative to placebo augmentation of intensive exposure therapy for youth with obsessive‐compulsive disorder (OCD) in a double‐blind, randomised controlled trial, and examine whether antidepressant medication or patient age moderated outcomes.MethodsYouth (n = 100, 7–17 years) with OCD were randomised in a 1:1 ratio to either DCS + exposure (n = 49) or placebo + exposure (n = 51). Assessments occurred posttreatment, 1 month later, and at 3 and 6 months. Pills were ingested immediately before sessions.ResultsSignificant improvements on all outcomes were observed at posttreatment, and to 6‐month follow‐up. Treatment arms did not differ across time, with no significant time‐by‐medication interactions on symptom severity (T1 to T2 estimate: 9.3, 95% confidence interval [CI]: −11.2 to −7.4, and estimate −10.7, 95% CI: −12.6 to −8.7), diagnostic severity (T1 to T2 estimate: −2.0, 95% CI: −2.4 to −1.5 and estimate −2.5, 95% CI: −3.0 to −2.0) or global functioning (T1 to T2 estimate: 13.8, 95% CI: 10.6 to 17.0, and estimate 16.6, 95% CI: 13.2 to 19.9). Neither antidepressants at baseline nor age moderated primary outcomes. There were significantly fewer responders/remitters at 1‐ and 6‐month follow‐up among youth in the DCS condition stabilised on SSRIs, relative to youth not taking SSRIs.ConclusionsDCS augmented intensive exposure therapy did not result in overall additional benefits relative to placebo. Intensive exposure proved effective in reducing symptoms for the overall sample.

MP3RAGE: Simultaneous mapping of T1 and in human brain at 7T.

To map T1 and the local flip angle ( ) in human brain using a single MP3RAGE sequence with 3 rapid acquisitions of gradient echoes (RAGEs). A third RAGE with a relatively high flip angle was appended to an MP2RAGE sequence. Through curve fitting and a rational approximation for small flip angles and short TR, closed form solutions for T1 and were derived. The influence of different k-space encoding schemes on precision and whether edge enhancement artifacts could be reduced with a saturation pulse applied prior to the third RAGE were explored. Validation of T1 estimates was performed using single-slice inversion recovery (IR) and a subsequent region-of-interest-based comparison, whereas validation of was performed using a whole brain pixelwise comparison to a DREAM flip angle mapping protocol. Lastly, MP3RAGE was compared to T1 -mapping by MP2RAGE with separate correction. Whole brain maps of T1 and at 1 mm isotropic resolution were obtained with MP3RAGE in 06:37 min. A linear-reverse centric-reverse centric phase-encoding order of the 3 RAGEs improved precision, and artifacts were successfully reduced with the saturation pulse. Estimations of T1 and deviated +2.5 ± 3.1% and -1.7 ± 8.6% from their respective references. T1 and can be mapped simultaneously using MP3RAGE. The approach can be thought of as combining MP2RAGE with a dual flip angle T1 -mapping protocol. Both maps can be solved for analytically and will be inherently co-registered at the high resolution associated with MPRAGE.

Comparison of T2 Quantification Strategies in the Abdominal-Pelvic Region for Clinical Use.

The aim of the study was to compare different magnetic resonance imaging (MRI) acquisition strategies appropriate for T2 quantification in the abdominal-pelvic area. The different techniques targeted in the study were chosen according to 2 main considerations: performing T2 measurement in an acceptable time for clinical use and preventing/correcting respiratory motion. Acquisitions were performed at 3 T. To select sequences for in vivo measurements, a phantom experiment was conducted, for which the T2 values obtained with the different techniques of interest were compared with the criterion standard (single-echo SE sequence, multiple acquisitions with varying echo time). Repeatability and temporal reproducibility studies for the different techniques were also conducted on the phantom. Finally, an in vivo study was conducted on 12 volunteers to compare the techniques that offer acceptable acquisition time for clinical use and either address or correct respiratory motion. For the phantom study, the DESS and T2-preparation techniques presented the lowest precision (ρ2 = 0.9504 and ρ2 = 0.9849 respectively), and showed a poor repeatability/reproducibility compared with the other techniques. The strategy relying on SE-EPI showed the best precision and accuracy (ρ2 = 0.9994 and Cb = 0.9995). GRAPPATINI exhibited a very good precision (ρ2 = 0.9984). For the technique relying on radial TSE, the precision was not as good as GRAPPATINI (ρ2 = 0.9872). The in vivo study demonstrated good respiratory motion management for all of the selected techniques. It also showed that T2 estimate ranges were different from one method to another. For GRAPPATINI and radial TSE techniques, there were significant differences between all the different types of organs of interest. To perform T2 measurement in the abdominal-pelvic region, one should favor a technique with acceptable acquisition time for clinical use, with proper respiratory motion management, with good repeatability, reproducibility, and precision. In this study, the techniques relying respectively on SE-EPI, radial TSE, and GRAPPATINI appeared as good candidates.

Accelerated cardiac T1 mapping in four heartbeats with inline MyoMapNet: a deep learning-based T1 estimation approach

PurposeTo develop and evaluate MyoMapNet, a rapid myocardial T1 mapping approach that uses fully connected neural networks (FCNN) to estimate T1 values from four T1-weighted images collected after a single inversion pulse in four heartbeats (Look-Locker, LL4).MethodWe implemented an FCNN for MyoMapNet to estimate T1 values from a reduced number of T1-weighted images and corresponding inversion-recovery times. We studied MyoMapNet performance when trained using native, post-contrast T1, or a combination of both. We also explored the effects of number of T1-weighted images (four and five) for native T1. After rigorous training using in-vivo modified Look-Locker inversion recovery (MOLLI) T1 mapping data of 607 patients, MyoMapNet performance was evaluated using MOLLI T1 data from 61 patients by discarding the additional T1-weighted images. Subsequently, we implemented a prototype MyoMapNet and LL4 on a 3 T scanner. LL4 was used to collect T1 mapping data in 27 subjects with inline T1 map reconstruction by MyoMapNet. The resulting T1 values were compared to MOLLI.ResultsMyoMapNet trained using a combination of native and post-contrast T1-weighted images had excellent native and post-contrast T1 accuracy compared to MOLLI. The FCNN model using four T1-weighted images yields similar performance compared to five T1-weighted images, suggesting that four T1 weighted images may be sufficient. The inline implementation of LL4 and MyoMapNet enables successful acquisition and reconstruction of T1 maps on the scanner. Native and post-contrast myocardium T1 by MOLLI and MyoMapNet was 1170 ± 55 ms vs. 1183 ± 57 ms (P = 0.03), and 645 ± 26 ms vs. 630 ± 30 ms (P = 0.60), and native and post-contrast blood T1 was 1820 ± 29 ms vs. 1854 ± 34 ms (P = 0.14), and 508 ± 9 ms vs. 514 ± 15 ms (P = 0.02), respectively.ConclusionA FCNN, trained using MOLLI data, can estimate T1 values from only four T1-weighted images. MyoMapNet enables myocardial T1 mapping in four heartbeats with similar accuracy as MOLLI with inline map reconstruction.

Fast B1 Mapping Based on Double-Angle Method with T1 Correction Using Standard Pulse Sequence.

Radiofrequency (RF) field (B1) mapping by combining the double-angle method (DAM) and T1 correction was investigated. The signal intensities S1 and S2 acquired by flip angle (FA) α and double FA 2α at short repetition time (TR) were converted to a signal intensity at TR=∞ by T1 correction. Then, these were used for DAM calculation. The T1 values are measured from two different images acquired with different TRs based on the saturation recovery (SR) method preliminarily. The effects of imaging parameters for T1 estimation and measured FA were investigated using CuSO4-doped water phantoms. A two-dimensional gradient echo type echo planar imaging pulse sequence was used. T1 values obtained by the 2-SR method were underestimated compared to the multipoint inversion recovery method. FA error was less than 5% when the appropriate imaging parameters were used. The acquisition time could be shortened to under 25 s by the use of T1-corrected DAM.

Repeatability and robustness of MP-GRASP T1 mapping.

To demonstrate the repeatability of fast 3D T1 mapping using Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI and its robustness to variation of imaging parameters including flip angle and spatial resolution in phantoms and the brain. Multiple imaging experiments were performed to (1) assess the robustness of MP-GRASP T1 mapping to B1 inhomogeneity using a single tube phantom filled with uniform MnCl2 liquid; (2) compare the repeatability of T1 mapping between MP-GRASP and inversion recovery-based spin-echo (IR-SE; over 12 scans), using a commercial T1MES phantom; (3) evaluate the longitudinal variation of T1 estimation using MP-GRASP with varying imaging parameters, including spatial resolution, flip angle, TR/TE, and acceleration rate, using the T1MES phantom (106 scans performed over a period of 12 months); and (4) evaluate the variation of T1 estimation using MP-GRASP with varying imaging parameters in the brain (24 scans in a single visit). In addition, the accuracy of MP-GRASP T1 mapping was also validated against IR-SE by performing linear correlation and calculating the Lin's concordance correlation coefficient (CCC). MP-GRASP demonstrates good robustness to B1 inhomogeneity, with intra-slice variability below 1% in the single tube phantom experiment. The longitudinal variability is good both in the phantom (below 2.5%) and in the brain (below 2%) with varying imaging parameters. The T1 values estimated from MP-GRASP are accurate compared to that from the IR-SE imaging (R2 = 0.997, Lin's CCC = 0.996). MP-GRASP shows excellent repeatability of T1 estimation over time, and it is also robust to variation of different imaging parameters evaluated in this study.

Feasibility of Magnetic Resonance Fingerprinting on Aging MRI Hardware.

The purpose of this work is to evaluate the feasibility of performing magnetic resonance fingerprinting (MRF) on older and lower-performance MRI hardware as a means to bring advanced imaging to the aging MRI install base. Phantom and in vivo experiments were performed on a 1.5T Siemens Aera (installed 2015) and 1.5T Siemens Symphony (installed 2002). A 2D spiral MRF sequence for simultaneous T1/T2/M0 mapping was implemented on both scanners with different gradient trajectories to accommodate system specifications. In phantom, for T1/T2 values in a physiologically relevant range (T1: 195–1539 ms; T2: 20–267 ms), scanners had strong correlation (R2 > 0.999) with average absolute percent difference of 8.1% and 10.1%, respectively. Comparison of the two trajectories on the newer scanner showed differences of 2.6% (T1) and 10.9% (T2), suggesting a partial explanation of the observed inter-scanner bias. Inter-scanner agreement was better when the same trajectory was used, with differences of 6.0% (T1) and 4.0% (T2). Intra-scanner coefficient of variation (CV) of T1 and T2 estimates in phantom were <2.0% and in vivo were ≤3.5%. In vivo inter-scanner white matter CV was 4.8% (T1) and 5.1% (T2). White matter measurements on the aging scanner after two months were consistent, with differences of 1.9% (T1) and 3.9% (T2). In conclusion, MRF is feasible on an aging MRI scanner and required only changes to the gradient trajectory.

Multi-exponential MRI T2 maps: A tool to classify and characterize fruit tissues

The estimation of multi-exponential relaxation time T2 and their associated amplitudes A0 at the voxel level has been made possible by recent developments in the field of image processing. These data are of great interest for the characterization of biological tissues, such as fruit tissues. However, they represent a high number of information, not easily interpretable. Moreover, the non-uniformity of the MRI images, which mainly directly impacts A0, could induce interpretation errors. In this paper, we propose a post-processing scheme that clusters similar voxels according to the multi-exponential relaxation parameters in order to reduce the complexity of the information while avoiding the problems associated with intensity non-uniformity. We also suggest a data representation suitable for the visualization of the multi-T2 distribution within each tissue. We illustrate this work with results for different fruits, demonstrating the great potential of multi-T2 information to shed new light on fruit characterization.

Axonal T2 estimation using the spherical variance of the strongly diffusion-weighted MRI signal

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T2 that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T2, as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T2 based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T2 estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system.

Kinetic modelling of dissolution dynamic nuclear polarisation 13 C magnetic resonance spectroscopy data for analysis of pyruvate delivery and fate in tumours.

Dissolution dynamic nuclear polarisation (dDNP) of 13 C-labelled pyruvate in magnetic resonance spectroscopy/imaging (MRS/MRSI) has the potential for monitoring tumour progression and treatment response. Pyruvate delivery, its metabolism to lactate and efflux were investigated in rat P22 sarcomas following simultaneous intravenous administration of hyperpolarised 13 C-labelled pyruvate (13 C1 -pyruvate) and urea (13 C-urea), a nonmetabolised marker. A general mathematical model of pyruvate-lactate exchange, incorporating an arterial input function (AIF), enabled the losses of pyruvate and lactate from tumour to be estimated, in addition to the clearance rate of pyruvate signal from blood into tumour, Kip , and the forward and reverse fractional rate constants for pyruvate-lactate signal exchange, kpl and klp . An analogous model was developed for urea, enabling estimation of urea tumour losses and the blood clearance parameter, Kiu . A spectral fitting procedure to blood time-course data proved superior to assuming a gamma-variate form for the AIFs. Mean arterial blood pressure marginally correlated with clearance rates. Kiu equalled Kip , indicating equivalent permeability of the tumour vasculature to urea and pyruvate. Fractional loss rate constants due to effluxes of pyruvate, lactate and urea from tumour tissue into blood (kpo , klo and kuo , respectively) indicated that T1 s and the average flip angle, θ, obtained from arterial blood were poor surrogates for these parameters in tumour tissue. A precursor-product model, using the tumour pyruvate signal time-course as the input for the corresponding lactate signal time-course, was modified to account for the observed delay between them. The corresponding fractional rate constant, kavail , most likely reflected heterogeneous tumour microcirculation. Loss parameters, estimated from this model with different TRs, provided a lower limit on the estimates of tumour T1 for lactate and urea. The results do not support use of hyperpolarised urea for providing information on the tumour microcirculation over and above what can be obtained from pyruvate alone. The results also highlight the need for rigorous processes controlling signal quantitation, if absolute estimations of biological parameters are required.

Detection of doxorubicin-induced cardiotoxicity using myocardial T1 and T2 relaxation times in childhood acute lymphoblastic leukemia survivors.

Doxorubicin leads to dose-dependent cardiotoxicity in childhood acute lymphoblastic leukemia (ALL) survivors. The first aim was to propose a contour-based estimation of T1 and T2 relaxation times based on the myocardial area, while our second aim was to evaluate native T1, post-gadolinium T1 and T2 relaxation time sensitivity to detect myocardial changes. A total of 84 childhood ALL survivors were stratified in regard to their prognostic risk groups: standard risk (SR), n = 20), high-risk with and without dexrazoxane (HR + DEX, n = 39 and HR, n = 25). Survivors' mean age was of 22.0 ± 6.9years, with a mean age at cancer diagnosis of 8.0 ± 5.2years. CMR acquisitions were performed on a 3T MRI system and included an ECG-gated 3(3)3(3)5 MOLLI sequence for T1 mapping and an ECG-gated T2-prepared TrueFISP sequence for T2 mapping. Myocardial contours were semi-automatically segmented using an interactive implementation of cubic Bezier curves. We found excellent repeatability between operators for native T1 (ICC = 0.91), and good repeatability between operators for post-gadolinium T1 (ICC = 0.84) and T2 (ICC = 0.79). Bland and Altman tests demonstrated a strong agreement between our contour-based method and images analyzed using the CVI42 software on the measure of native T1, post-gadolinium T1, and T2. No significant differences between survivors' prognostic risk groups in native T1 were reported, while we observed significant differences between survivors' prognostic risk groups in post-gadolinium T1 and T2. Significant differences were observed between male and female survivors. Differences between groups were also observed in partition coefficients, but no significant differences were observed between male and female survivors. The use of CMR parameters with native T1, post-gadolinium T1, and T2 allowed to show that survivors at a high-risk prognostic were more exposed to doxorubicin-related cardiotoxicity than those who were at a standard risk prognostic or who received dexrazoxane treatments.